A constant denominator perturbation theory for molecular energy

A constant denominator perturbation theory is developed based on a zeroth order Hamiltonian characterized by degenerate subsets of orbitals. Such a formulation allows for a decoupling of the numerators of the perturbation sequence, allowing for much more rapid evaluation of the resultant sums. For example, the full fourth order theory can be evaluated as an N ⁶ step rather than N ⁷, where N is proportional to the basis set.Although the theory is general, a constant denominator is chosen for this study as the difference between the average occupied and average virtual orbital energies scaled so that the first order wavefunction yields the lowest possible variational bound. The third order correction then appears naturally as a scaled Langhoff-Davidson correction. The full fourth order with this partitioning is developed. Results are presented within the localized bond model utilizing both the Pariser-Parr-Pople and CNDO/2 model Hamiltonians. The second order theory presents a useful bound, usually containing a good deal of the basis set correlation. In all cases examined the fourth order theory shows remarkable stability, even in those cases in which the Nesbet-Epstein partitioning seems poorly convergent, and the Moller-Plesset theory uncertain.     

On the use of spin graphs for spin adapting many-body perturbation theory

In connection with spin adaptation in many-body perturbation theory, this paper reexamines the use of spin graphs in view of a Hugenholtz-like representation where both the orbital and the spin parts coexist. Together with the idea of essentially distinct diagrams, this representation leads to an economic handling of spin adaptation. As a side issue, the appropriate generalization of the Epstein–Nesbet partitioning for spin-adapted perturbation theory is obtained. Compact formulas up to fourth order of the ground-state energy, and up to third order for excitation energies and ionization potentials are given.     

Bond length alternation in cyclic polyenes. V. Local finite-order perturbation theory approach

The finite-order many-body perturbation theory using the localized Wannier orbital basis is applied to the problem of bond length alternation in the Pariser–Parr–Pople model of cyclic polyenes C_N H_N, N = 4v + 2, which may be regarded as a simplified model of polyacetylene. Both the Møller–Plesset and the Epstein–Nesbet-type partitionings of the model Hamiltonian are employed. The localized orbital basis enables an efficient truncation of the perturbation theory summations over the intermediate states as well as an elimination of energetically unimportant diagrams, thus enabling one to obtain the fourth-order Møller–Plesset-perturbation energies with a relatively small computational effort even for large polyenes. The results obtained with the second-, third-, and fourth-order Møller–Plesset and with the third-order Epstein–Nesbet perturbation theories yield very similar bond length distortions (about 0.05 Å) and stabilization energies per site (about 0.04 eV) as obtained earlier with the RHF , one-parameter AMO , and delocalized orbital perturbation theories. The effects of truncation and diagram elimination in the fourth-order Møller–Plesset perturbation theory and the abnormal behavior of the second-order Epstein–Nesbet perturbation theory results in the localized Wannier basis near the instability threshold of the RHF solutions are discussed.     

A nondiagonal quasidegenerate fourth-order perturbation theory

A canonical quasidegenerate Rayleigh-Schrödinger perturbation theory, correct through fourth order in the energy, is explored for a block-diagonal unperturbed Hamiltonian. The theory is developed completely within a Lie Algebra in Hilbert space. Explicit equations forn-particle transition elements in terms of solutions of simultaneous linear equations are presented. A two-dimensional anisotropic anharmonic oscillator is used to provide numerical results. The perturbation theory is shown to be stable under small separation of model and complement spaces. An iterative variant of the fourth order perturbation theory is introduced; the iterative variant is related to the non-iterative one in much the same way as nondegenerate coupled cluster theories are related to nondegenerate perturbation theory. The quasidegenerate coupled cluster theory appears to be stable in the presence of multiple intruder states.     

Second order of perturbation theory correction to the radial distribution function of a liquid with the interaction potential of hard spheres plus a square-well

A liquid with the interaction potential of hard spheres plus a square-well is analyzed using the Monte-Carlo technique. Numerical results for the perturbation theory series over a square-well potential are obtained in the form of the Barker and Henderson discrete representation. Approximating expressions for the correction to a liquid radial distribution function in the second order of perturbation theory are presented. The obtained results allow us to define this correction with a root-mean-square deviation of about 0.007. It is shown that the given approach provides a complete calculation in the second order of perturbation theory, and also the determination of the third order correction to the free energy for a liquid interacting with the potential of the Lennard-Jones type.     

Perturbation theory to determine the stable isomer of triatomic linear molecules

The electronegativity perturbation theory up to the fourth order has been used to study the preferred isomer of linear triatomic molecules. Systems with 10–16 valence–electron systems are discussed and our results are in good agreement with the experiment. The only exception is the isomers of PNO. The usefulness of the perturbation theory is briefly discussed. © 1993 John Wiley & Sons, Inc.     

The calculation of vibrational energy levels of polyatomic molecules including anharmonic effect using contact transformation perturbation method

Using contact transformation perturbation method based on the Taylor expansion of the potential energy function in terms of dimensionless normal coordinates up to sixth‐order, the vibrational energy levels in terms of force constants are derived. The contact transformation theory has been applied to simplify the calculation of perturbation effects. To calculate the second‐order vibrational energy correction, the third and fourth‐order terms of potential function have been placed in the first‐order perturbation Hamiltonian and the second‐order Hamiltonian contains hexatic ones. We present expressions which give relations between the fourth‐ and sixth‐order terms in dimensionless normal coordinates of the potential and the anharmonicity coefficients. For illustration, a set of vibrational energies levels of SO₂, and H₂O molecules including anharmonic effects has been calculated. © 2013 Wiley Periodicals, Inc.     

Application of the many-body perturbation theory by using localized orbitals

Diagrammatic formulation of the many-body perturbation theory is investigated when both the occupied orbitals and the virtual ones are localized, i.e., they are unitary transforms of the canonical Hartree–Fock orbitals. All diagrams representing ground state correlation energy can be generated through fifth order. For cyclic polyenes C₆H₆ and C₁₀H₁₀ as model systems, the energy corrections are calculated in the Pariser–Parr–Pople approximation for a wide range of the coupling constant β^?1, through fourth order including some fifth order terms. The results are compared to those obtained by other methods: perturbation theory by using canonical orbitals and full CI. The effect of neglecting contributions from orbitals localized into neighboring sites is also studied.     

Structure, force fields, and vibrational spectra of cerium tetrahalides

The geometrical structure, force fields, and vibrational spectra of CeX₄ (X = F, Cl, Br, I) were investigated by second, third, and fourth order Möller-Plesset perturbation theory, CISD+Q configuration interaction method, and the CCSD(T) coupled cluster method. Calculations on CeF₄ were also performed by multiconfiguration second order perturbation theory MCQDPT2/CASSCF. The wave function of the ground state of CeX₄ molecules was found to be appreciably non-one-configurational; this property increases from cerium fluorides to iodides and leads to the divergence of the series of Möller-Plesset perturbation theory. The calculated data point to a tetrahedral equilibrium nuclear configuration in CeX₄ molecules. The energy barriers to the inversion of the tetrahedral CeX₄ molecules via the square configurations are high enough, 74–89 kJ/mol. The calculated vibration frequencies, effective internuclear distances, and mean amplitudes of nuclear vibrations in CeF₄ agree with IR and Raman spectroscopic and high-temperature gas-phase electron diffraction data.     

Perturbation theory for multiconfiguration reference states

Vector method procedures are adapted to evaluate Rayleigh-Schrödinger perturbation corrections to a multiconfiguration zeroth order function. If this function is sufficiently flexible, this perturbation theory can be applied to low lying excited states. The effectiveness of our theory is demonstrated on the ground state of F₂ and the low lying excited states of Mg₂. Energies calculated through fourth order are compared with appropriate CI results.     

Energy analysis of weak electron-donor-acceptor complexes and water clusters with the perturbation theory based on the locally projected molecular orbitals: charge-transfer and dispersion terms

The third order single excitation perturbation theory corrected with the dispersion energy based on the locally projected molecular orbital was applied to study the weak electron-donor-acceptor (charge-transfer) complexes and the hydrogen bonds in the water clusters. In the weak electron-donor-acceptor complexes, the dispersion energy is larger than the charge-transfer energy in absolute value. The dispersion energy is as large as the charge-transfer energy in the hydrogen bond. The cage form of (H(2)O)(6) is the most stable among eight isomers examined, because the dispersion energy is the largest among them.     

A new size extensive multireference perturbation theory

A new multireference perturbation series is derived based on the Rayleigh–Schrödinger perturbation theory. It is orbitally invariant. Its computational cost is comparable to the single reference Møller–Plesset perturbation theory. It is demonstrated numerically that the present multireference second‐ and third‐order energies are size extensive by two types of supermolecules composed of H₂ and BH monomers. Spectroscopic constants of as well as the ground state energies of H₂O, NH₂, and CH₂ at three bond lengths have been calculated with the second multireference perturbation theory. The dissociation behaviors of CH₄ and HF have also been investigated. Comparisons with other approximate theoretical models as well as the experimental data have been carried out to show their relative performances. © 2013 Wiley Periodicals, Inc.     

Non-fused polyaromatic hydrocarbons: interactions of aromatic and antiaromatic rings through a CC bond

The interaction of the moieties of benzene, cyclobutadiene, cyclopentadinyl anion, and the cyclopentadianide cation upon each other and upon a CC bond connecting pairs of these rings is investigated computationally. The resulting non-fused bicycles include biphenyl, phenylcyclobutadiene, phenylcyclopentadienylium, phenylcyclopentadienide, pentafulvalene, cyclobutadienyl–cyclopentadienylium, cyclobutadienyl–cyclopentadienide, and bicyclobutadiene. The relative stability and aromaticity are assessed from hydrogenation energies, aromatic stabilization energies, ring separation energies, nucleus-independent chemical-shift, harmonic oscillator model of aromaticity, and natural bond orbital analysis. Calculations are performed with density functional theory (B3LYP) and Møller–Plesset perturbation theory of second order (MP2). Enthalpy quantities are also determined by G3. When both rings are aromatic in character, the bridging bond is mostly σ in character. When one or both of the rings is antiaromatic, the bridging bond has significant π character. Systems with contrasting aromaticities have CC bridging bonds of lengths between CC single bond lengths and CC double bond lengths and where the systems were charged, the charge is evenly distributed between the rings.     

CCSD(T) expectation value calculations of first-order properties

An expectation value approach to calculations of first-order properties using the non-iterative, triple-excitation amplitudes in the coupled cluster wave function is exploited. Three methods are suggested and analysed using the many body perturbation theory (MBPT) expansion arguments. The first method, in which non-iterative triple-excitation amplitudes are used in the expression for the expectation values, makes the wave function accurate through the second order of MBPT. In the second method, which is an extension of the first, effects of triple-excitation amplitudes are coupled with single- and double-excitation amplitudes. The correlated density matrix equivalent through the fourth order to that obtained when CCSDT-la amplitudes are used is employed in the third method. The suggested methods are tested on dipole moment and polarizability calculations for several diatomic closed-shell molecules and are compared to other related approaches. Received: 15 May 1997 / Accepted: 5 June 1997     

Coordination bonding in dicopper and dichromium tetrakis(μ-acetato)-diaqua complexes: Nature,strength, length,and topology

Geometry optimization, energetics, electronic structure, and topology of electron density of dicopper ( I ) and dichromium ( II ) tetrakis(μ-acetato)-diaqua complexes are studied focusing on the metal–metal interactions. The performance of broken symmetry (BS) single-determinant ab initio (Hartree–Fock, Møller–Plesset perturbation theory to the second and third orders, coupled clusters singles and doubles) and density functional theory (BLYP, B3LYP, B3LYP-D3, B2PLYP, MPW2PLYP) methods is compared to multideterminant ab initio (CASSCF, NEVPT2) methods as well as to the multipole model of charge density from a single-crystal X-ray diffraction experiment (Herich et al., Acta Cryst. 2018, B74, 681–692). In vacuo DFT geometry optimizations (improper axial water ligand orientation) are compared against the periodic ones. The singlet state is found to be energetically preferred. J coupling of ( I ) becomes underestimated for all ab initio methods used, when compared to experiment. It is concluded that the strength of the direct M─M interactions correlates closely with the J coupling magnitude at a given level of theory. The double potential well character of (II) and of the dehydrated form of (II) are considered with respect to the Cr─Cr distance. The physical effective bond order of the metal–metal interaction is small (below 0.1 e) in ( I ) and moderate (0.4 e) in ( II ). The CASSCF results overestimate the electron density of the metal–metal bond critical point by 20% and 50% in ( I ) and ( II ), respectively, when compared to the multipole model. © 2019 Wiley Periodicals, Inc.     

Improvement on macroscopic compressibility approximation and beyond

A numerical procedure is proposed to extend the thermodynamic perturbation expansion (TPE) to a higher order. It is shown that the present second order term is superior to that due to a macroscopic compressibility approximation (MCA), a local compressibility approximation, and a superposition approximation by Barker and Henderson [Rev. Mod. Phys. 48, 587 (1976)]. Extensive model calculation and comparison with simulation data available in literature and supplied in the present report indicate that the present third order TPE is superior to a previous second order TPE based on the MCA, two previous perturbation theories, which are respectively based on an analytical mean spherical approximation for an Ornstein-Zernike equation, and an assumed explicit functional form for the Laplace transform of radial distribution function multiplied by radial distance, and a recent generalized van der Waals theory. The present critical temperature for a hard core attractive Yukawa fluid of varying range is in very good agreement with that due to a hierarchical reference theory. The present third order TPE is computationally far more modest than the self-consistent integral equation theory, and therefore is a viable alternative to use of the latter.     

Further test of third order + second-order perturbation DFT approach: hard core repulsive Yukawa fluid subjected to diverse external fields

Grand canonical Monte Carlo simulation is used to investigate density profiles of hard-core repulsive Yukawa (HCRY) model fluid under the influence of various external fields and radial distribution function (RDF) of the bulk HCRY system. The aim of these extensive simulations is to provide exact data for purely repulsive interaction potential against which the validity of a third order + second-order perturbation DFT approach can be tested. It is found that a semiempirical parametrized bridge function due to Malijevsky and Labik performs very well for the RDF of the bulk HCRY fluid. Incorporation of a bulk second-order direct correlation function (DCF) of the HCRY fluid based on the Malijevsky-Labik bridge function into the third order + second-order perturbation DFT approach yields the resulting theoretical predictions for the density profiles of inhomogeneous HCRY fluid that are in a very good agreement with the simulation data, an exception being somewhat larger deviations appearing for the structure of the fluid around the center of a hard spherical cavity. Both theory and simulation predict layering transition and gas-liquid coexistence phenomena occurring with the HCRY model fluid under confined conditions. For the case of an inverse sixth-power repulsive potential under the influence of a flat stationary wall defined by an inverse twelfth-power repulsive potential, the present third order + second-order perturbation DFT approach is found to be superior to several existing weighted density approximations (WDA) and partitioned WDA.     

Frequency dependent nonlinear optical properties of molecules: Formulation and implementation in the HONDO program

This article summarizes the detailed equations for the time-dependent Hartree–Fock treatment of nonlinear properties for perturbations made up of a static electric field and an oscillating field. Explicit expressions for all nonlinear processes up to third order are obtained in terms of the density matrices at the same order. For processes at second and third order in perturbation, expressions in terms of lower order quantities are also obtained by applying the (2n + 1) theorem of perturbation theory. The corresponding computer implementation in the HONDO program is described.     

Tests of Monte Carlo perturbation theory for the free energy of liquid copper

Monte Carlo perturbation theory, in which terms in the thermodynamic perturbation series are evaluated by Monte Carlo averaging, has potentially large advantages in efficiency for calculating free energies of liquids from ab initio potential surfaces. In order to test the accuracy of perturbation theory for liquid metals, a series of calculations has been done on liquid copper, modeled by an embedded atom potential. A simple 1/r(12) pair potential is used as the reference system. The free energy is calculated to third order in perturbation theory, and the results are compared to an exact formula. It is found that for optimal reference potential parameters, second order perturbation theory is essentially exact. Second and third order theories give accurate results for significantly nonoptimal reference parameters. The relation between perturbation theory and reweighting is discussed, and an approximate formula is derived that shows an exponential dependence of the efficiency of reweighting on the second order free energy correction. Finally, techniques for application to ab initio potentials are discussed. It is shown that with samples of 100 configurations, it is possible to obtain accuracy and precision at the level of approximately 1 meV/atom.