Calculation of vertical ionization potentials of H2O and Ne by many-body Rayleigh-Schrödinger perturbation theory

The vertical ionization potentials of H₂O and Ne are calculated by many-body Rayleigh-Schrödinger perturbation theory up to third order. The comparison of the present method with the other approaches is done.     

Assessing the Accuracy of SAPT(DFT) Interaction Energies by Comparison with Experimentally Derived Noble Gas Potentials and Molecular Crystal Lattice Energies

The density functional version of symmetry‐adapted perturbation theory, SAPT(DFT), is a computationally efficient method for calculating intermolecular interaction energies. We evaluate its accuracy by comparison with experimentally determined noble gas interaction potentials and sublimation enthalpies, most of which have not been previously calculated using this method. In order to compare the results with wavefunction methods, we also calculate these quantities using MP2 and, for noble gas dimers, using CCSD(T). For the crystal lattice energy calculations, we include corrections to the dispersion, electrostatic, and induction energies that account for the finite interaction distance cutoff and higher‐order induction contributions. Overall, the energy values extrapolated to the complete basis set limit show that SAPT(DFT) achieves significantly better agreement with experiment than MP2.     

The hydrogen perturbation in molecular connectivity computations

A new algorithm for the delta(v) number, the basic parameter of molecular connectivity indices, is proposed. The new algorithm, which is centered on graph concepts like complete graphs and general graphs, encodes the information of the bonded hydrogen on different atoms through a perturbation parameter that makes use of no new graph concepts. The model quality of the new algorithm is tested with 13 properties of seven different classes of compounds, as well as with composite classes of compounds with the same property and with composite properties of the same class of compounds. Chosen properties and classes of compounds display different percentage of bonded hydrogen atoms, which allow a checking of the importance of this parameter. A comparison is drawn with previous results with zero contribution for the hydrogen perturbation as well as among results obtained by changing the number of compounds of a property but keeping constant the percentage of hydrogen atoms. Results underline the importance of the property as well as the importance of the number of compounds in determining the level of the hydrogen perturbation. Molecular connectivity terms are in some cases more critical than the combination of indices in detecting the perturbation introduced by the hydrogen atoms.     

Multi-reference-state perturbation theory for computation of potential-energy surfaces

The implementation of multi-reference-state Rayleigh-Schrödinger perturbation theory for the evaluation of potential-energy surfaces is reviewed. The organization of the computation and the basic algorithms are outlined. A major difficulty is the quadratic increase in computing sources as the reference space is expanded. Truncating the space of nonreference determinants is suggested as a solution for the problem. The method applies a numerical cutoff criterion for the weights of the first-order wave function. The truncation reduces the computing time significantly with negligible sacrifice in the quality of the results.This paper was presented at the international conference on The Impact of Supercomputers on Chemistry, held at the University of London, London, UK, 13–16 April 1987     

Anharmonic Corrections to Structural Experiment Data

Anharmonic corrections to internuclear distances observed in diffraction and microwave experiments are calculated at the level of first-order perturbation theory. The results obtained for a simple system with the use of a quantum-mechanical potential function are compared with the values calculated according to currently used approximations. It is shown that the diatomic approximation enormously exaggerates anharmonic corrections to distances between nonhydrogen atoms.     

The potential of mean force on halide ions near the Cu(100) surface

This paper reports an investigation of the phenomenon of specific adsorption of halide ions on a Cu(100) surface using Monte Carlo simulations. The system was modeled by considering each ion in a water lamina placed between two copper walls. The potentials used in simulations were constructed by fitting to results of quantum calculations. The solvent contribution to the potential of mean force (pmf) was calculated for each of the four halide ions using the free energy perturbation method. Given the difficulty of finding a reliable ion–metal potential, several alternatives, representing extremal models, were tested in combination with the solvent mean force on the ions, F⁻, Cl⁻, Br⁻ or I⁻. The results for the pmf on an ion near the metal surface are discussed in the light of the experimental data available. The sensitivity of the results to the type of ion–metal potential used in the simulations is stressed.     

氯化钠水溶液的monte Carlo分子模拟研究

采用MonteCarlo计算机分子模拟方法,研究了氯化钠水溶液在常温和高温情况下粒子的径向分布函数和粒子间各种作用位能的热力学性质。模拟过程采用了NTV正则系综,粒子间的作用能包括离子库仑静电作用,偶极子作用以及色散作用,这些作用构成了电解质溶液的基本框架,模拟结果与微扰理论和平均球近似积分方程理论的预测值进行了比较。     

Dimers of Nineteen‐Electron Sandwich Compounds: Crystal and Electronic Structures,and Comparison of Reducing Strengths

The dimers of some Group 8 metal cyclopentadienyl/arene complexes and Group 9 metallocenes can be handled in air, yet are strongly reducing, making them useful n‐dopants in organic electronics. In this work, the X‐ray molecular structures are shown to resemble those of Group 8 metal cyclopentadienyl/pentadienyl or Group 9 metal cyclopentadienyl/diene model compounds. Compared to those of the model compounds, the DFT HOMOs of the dimers are significantly destabilized by interactions between the metal and the central C?C σ‐bonding orbital, accounting for the facile oxidation of the dimers. The lengths of these C?C bonds (X‐ray or DFT) do not correlate with DFT dissociation energies, the latter depending strongly on the monomer stabilities. Ru and Ir monomers are more reducing than their Fe and Rh analogues, but the corresponding dimers also exhibit much higher dissociation energies, so the estimated monomer cation/neutral dimer potentials are, with the exception of that of [RhCp₂]₂, rather similar (?1.97 to ?2.15 V vs. FeCp₂^+/0 in THF). The consequences of the variations in bond strength and redox potentials for the reactivity of the dimers are discussed.     

Direct calculation of ionization potentials of closed-shell atoms and molecules

Vertical ionization potentials, electron affinities and information about quasi-particles can be obtained by using the technique of the single-particle propagator. The expansion of the self-energy part up to third order perturbation theory can be evaluated numerically, but does not lead, in most cases, to satisfying results. A theoretical and numerical analysis of the diagrammatic expansion of the self-energy part requires the introduction of a renormalized interaction and renormalized hole and particle lines.     

Electrochemical and electrochromic properties of rare-earth metal diphthalocyanine complexes

The electrochemical and spectroelectrochemical properties of several new rare-earth metal diphthalocyanine complexes with different substituents containing both electron-withdrawing and -donating substituents in the phthalocyanine rings were studied. The influence of structural modification of the phthalocyanine complexes (viz., the nature of the central metal atom and substituents in the phthalocyanine rings and the number of phthalocyanine rings in the complexes) on the total number of redox transitions, their potentials, and spectral characteristics of anionic and cationic forms of the complexes was examined. The potentials of the first anodic and cathodic redox transitions of the diphthalocyanine complexes are in a good linear correlation with the ionic radii of lanthanides. The potentials of the observed redox transitions of the complexes under study correlate well with the sums of Hammett constants for the substituents and the minimum molecular electrostatic potential of the benzene rings in the phthalocyanine moiety, which serves as a measure of electron perturbations introduced by the substituents. The revealed regularities allow the prediction of the redox properties and structure of rare-earth element diphthalocyanine complexes, which are redox-active in a specified potential range.     

Ab initio quasirelativistic calculations on electronic transitions in ICl by the multireference many‐body perturbation theory

A quasirelativistic perturbative method of ab initio calculations on ground and excited molecular electronic states and transition properties within the relativistic effective core potential approximation is presented and discussed. The method is based on the construction of a state‐selective many‐electron effective Hamiltonian in the model space spanned by an appropriate set of Slater determinants by means of the second‐order many‐body multireference perturbation theory. The neglect of effective spin–orbit interactions outside of the model space allows the exploitation of relatively high nonrelativistic symmetry during the evaluation of perturbative corrections and therefore dramatic reduction of the cost of computations without any contraction of the model‐space functions. One‐electron transition properties are evaluated via the perturbative construction of spin‐free transition density matrices. Illustrative calculations on the X0⁺ ? A1, B0⁺, and (ii)1 transitions in the ICl molecule are reported. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Geometry optimizations with explicit inclusion of electron correlation

Automatic techniques for geometry optimization are applied in conjunction with configuration interaction and perturbation treatments of electron correlation. The computational effort and numerical accuracy of the optimizations are discussed, as well as problems with approximate correlation methods concerning the continuity of the potential surface. The optimized geometries of fourteen molecules obtained with different correlation treatments (MNDO SCF MOs) are compared. The configuration interaction results are reproduced satisfactorily by simple perturbation approaches. The largest change of the optimized SCF geometry is found for hydrogen peroxide.     

Cooperative and Diminutive Effects of Pnicogen Bonds and Cation–π Interactions

The interplay between pnicogen bonds and cation–π interactions has been investigated at the MP2/aug‐cc‐pVDZ level. Interesting cooperative and diminutive effects are observed when pnicogen bonds and cation–π interactions coexist in the same complex. These effects have been analyzed in terms of the structural, energetic, and charge‐transfer properties of the complexes. The variations in electron density at critical points of the intermolecular bond have been used to analyze bond strengthening or weakening. The nature of the interactions and the mechanisms of cooperative and diminutive effects have been studied by means of symmetry‐adapted perturbation theory and molecular electrostatic potentials.     

A Green's function approach to the photoelectron spectrum of bis(π-allyl)nickel[1]

The vertical ionization potentials of bis(-allyl)nickel (see (1) in Fig. 1) are calculated by means of the Green's function approach within a semiempirical INDO extension to the first transition metal series. The computed ionization potentials are in good agreement with an experimentally deduced assignment. In contrast to earlier theoretical and experimental studies, the 7a _u () level is predicted on top of the levels corresponding to the Ni 3d orbitals. Our approach leads to a complete assignment of the PE spectrum of (1) in the outer valence region.     

Uranyl-Bound Tetra-Dentate Non-Innocent Ligands: Prediction of Structure and Redox Behaviour Using Density Functional Theory

Computational methods have been applied to understand the reduction potentials of [UO₂-salmnt-L] complexes (L=pyridine, DMSO, DMF and TPPO), and their redox behavior is compared with previous experiments in dichloromethane solution. Since the experimental results were inconclusive regarding the influence of the uranyl-bound tetra-dentate ‘salmnt’ ligand, here we will show that salmnt acts as a redox-active ligand and exhibits non-innocent behavior to interfere with the otherwise expected one-electron metal (U) reduction. We have employed two approaches to determine the uranyl (VI/V) reduction potentials, using a direct study of one-electron reduction processes and an estimation of the overall reduction using isodesmic reactions. Hybrid density functional theory (DFT) methods were combined with the Conductor-like Polarizable Continuum Model (CPCM) to account for solvation effects. The computationally predicted one-electron reduction potentials for the range of [UO₂-salmnt-L] complexes are in excellent agreement with shoulder peaks (∼1.4 eV) observed in the cyclic voltammetry experiments and clearly correlate with ligand reduction. Highly conjugated pi-bonds stabilize the ligand based delocalized orbital relative to the localized U f-orbitals, and as a consequence, the ligand traps the incoming electron. A second reduction step results in metal U(VI) to U(V) reduction, in good agreement with the experimentally assigned uranyl (VI/V) reduction potentials.     

Thermodynamic perturbation theory of simple liquids

It was proven that after averaging over the canonical Gibbs ensemble, the mean perturbation energy was singled out of the classical partition function before the expansion in a series of perturbation theory. Therefore, the term that formally coincides with first order perturbation theory in a decomposition of the Helmholtz free energy bears no relationship to perturbation theory. Then the proper series of the thermodynamic perturbation theory always starts with a second order infinitesimal. Therefore, the wellknown condition of applicability of the thermodynamic perturbation theory, “...the requirement that the perturbation energy per particle be small compared with T...” (L. D. Landau and E. M. Livshits, Statistical Physics, Vol. V, Pt. I), can be substantially weakened. The most important factor for applicability of thermodynamic perturbation theory is the value of many-particle correlations in an unperturbed system, but not the smallness of the perturbation potential.     

Hydrogen Bonding between Metal‐Ion Complexes and Noncoordinated Water: Electrostatic Potentials and Interaction Energies

The hydrogen bonding of noncoordinated water molecules to each other and to water molecules that are coordinated to metal‐ion complexes has been investigated by means of a search of the Cambridge Structural Database (CSD) and through quantum chemical calculations. Tetrahedral and octahedral complexes that were both charged and neutral were studied. A general conclusion is that hydrogen bonds between noncoordinated water and coordinated water are much stronger than those between noncoordinated waters, whereas hydrogen bonds of water molecule in tetrahedral complexes are stronger than in octahedral complexes. We examined the possibility of correlating the computed interaction energies with the most positive electrostatic potentials on the interacting hydrogen atoms prior to interaction and obtained very good correlation. This study illustrates the fact that electrostatic potentials computed for ground‐state molecules, prior to interaction, can provide considerable insight into the interactions.     

Theoretical study on the size consistency of the second and third order energies of the multireference perturbation theory

The size consistency of the second and third order energies of the multireference perturbation theory(Chen F, Davidson E, Iwata S. Int J Quant Chem, 2002, 86: 256) is investigated theoretically with a su-per-molecular model composed of N-hydrogen molecules separated by a large distance. It is found that the two perturbation series corresponding to two Hamiltonian partitions are not size consistent at the second and third order. However, two size consistent forms are suggested for two Hamiltonian parti-tions at the second order, if some approximations to the denominators of the original second order energies are assumed.     

Atom-type assignment in molecules and clusters by perturbation theory-A complement to X-ray structure analysis

An approach to distinguish elements with similar atomic numbers in molecules and clusters is presented and applied to experimentally synthesized and structurally characterized mixed-metallic compounds. By first treating a homogenized reference system constructed from the original compound and applying first-order perturbation theory it is possible to find the most stable distribution of the atom types to the atomic sites in a very efficient way. This work is focused on the appropriate choice of homogenized reference systems and on applications treating experimentally synthesized compounds. With these examples is shown that the method is a helpful complement to X-ray crystal structure analysis.     

Free‐energy differences between states with different conformational ensembles

Multiple conformations separated by high‐energy barriers represent a challenging problem in free‐energy calculations due to the difficulties in achieving adequate sampling. We present an application of thermodynamic integration (TI) in conjunction with the local elevation umbrella sampling (LE/US) method to improve convergence in alchemical free‐energy calculations. TI‐LE/US was applied to the guanosine triphosphate (GTP) to 8‐Br‐GTP perturbation, molecules that present high‐energy barriers between the anti and syn states and that have inverted preferences for those states. The convergence and reliability of TI‐LE/US was assessed by comparing with previous results using the enhanced‐sampling one‐step perturbation (OSP) method. A linear interpolation of the end‐state biasing potentials was sufficient to dramatically improve sampling along the chosen reaction coordinate. Conformational free‐energy differences were also computed for the syn and anti states and compared to experimental and theoretical results. Additionally, a coupled OSP with LE/US was carried out, allowing the calculation of conformational and alchemical free energies of GTP and 8‐substituted GTP analogs. © 2013 Wiley Periodicals, Inc.