Multiconfiguration perturbation theory: size consistency at second order

A modified version of a previously elaborated multiconfiguration perturbation theory (MCPT) [Rolik et al. J. Chem. Phys. 119, 1922 (2003)] is presented. In the modified formulation size consistency is ensured at second order in energy, by omitting projectors from the zero order Hamiltonian operator. This MCPT formulation is abbreviated as SC2-MCPT (size consistent at second order). To ensure proper separability, we also require that energy denominators are constructed as differences of some one-particle energies. A similar choice for energy denominators also renders the well-known multireference Moller-Plesset (MRMP) energy size consistent at second order. The same thing applies to the related multireference perturbation theory by Witek, Nakano, and Hirao.     

两种哈密顿量划分的单参考态微扰展开式的大小一致性研究

利用一个含N个氢分子的超分子模型, 从解析的角度仔细地分析了相应于两种哈密顿量划分方式的单参考态微扰展开式前四级能量的大小一致性. 对于Epstein-Nesbet划分, 微扰展开式是大小一致的. 而对于另一种划分(Chen, F; Davidson, E. R.; Iwata, S. Int. J. Quantum Chem., 2002, 86: 256, 见公式(24)的讨论), 其展开式大小不一致. 但将它的分母作泰勒展开, 经过重新整理, 将那些大小不一致的项消去后, 得到的新微扰展开式是逐项大小一致的.     

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.     

A method of incorporating quadruple correction in the scheme of multi-reference singly and doubly excited configuration interaction — a CSF based coupled pair approximation

Summary We develop an approximate size consistent method within a framework of the multi-reference configuration interaction scheme. The Rayleigh-Schrödinger perturbation theory is employed with a specific selection of the unperturbed part of the electronic Hamiltonian. The second order energy is obtained by a set of equations similar to the quasidegenerate variational perturbation theory of Cave and Davidson. The approximate fourth order energy is obtained by solving a set of equations similar to the coupled electron pair approximation (CEPA). The method has been tested for two simple systems, BeH₂ and N₂, and the results are quite encouraging.     

New perspectives in multireference perturbation theory: the n-electron valence state approach

The n-electron valence state perturbation theory (NEVPT) is a form of multireference perturbation theory which is based on a zero order reference wavefunction of CAS-CI type (complete active space configuration interaction) and which is characterized by the utilization of correction functions (zero order wavefunctions external to the CAS) of multireference nature, obtained through the diagonalization of a suitable two-electron model Hamiltonian (Dyall’s Hamiltonian) in some well defined determinant spaces. A review of the NEVPT approach is presented, starting from the original second order state-specific formulation, going through the quasidegenerate multi-state extension and arriving at the recent implementations of the third order in the energy and of the internally contracted configuration interaction. The chief properties of NEVPT—size consistence and absence of intruder states—are analyzed. Finally, an application concerning the calculation of the vertical spectrum of the biologically important free base porphin molecule, is presented.     

Localized electron pair theory for the calculation of ground state energies of large molecules

The quantum mechanical energy is examined in which groups of one, two, three, and four localized electron pairs found within a molecule are separately computed. From these results, the interaction energies of the electron pairs taken one, two, three, and four at a time form the terms of a convergent molecular mechanics like expansion of the molecular ground state energy. This procedure can be used with any size consistent quantum mechanical method. The computational time for large molecules depends chiefly upon the order needed in the energy expansion to obtain sufficient convergence and not on the particular quantum mechanical method used. Preliminary results within the framework of a semiempirical CNDO/2 model Hamiltonian show at the Hartree–Fock and Møller–Plesset perturbation levels that relative energies converge to within a few tenths of a kcal/mol of the exact values at the four body level for molecules that have little delocalization. In strained ring and aromatic systems, convergence is however not nearly as rapid. Results can be improved somewhat by using larger interacting fragments containing two or more electron pairs over three or more atomic centers. © 1992 by John Wiley & Sons, Inc.     

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.     

Perturbational approach to the electron correlation cusp applied to helium‐like atoms

A recently proposed perturbational approach to the electron correlation cusp problem 1 is tested in the context of three spherically symmetrical two‐electron systems: helium atom, hydride anion, and a solvable model system. The interelectronic interaction is partitioned into long‐ and short‐range components. The long‐range interaction, lacking the singularities responsible for the electron correlation cusp, is included in the reference Hamiltonian. Accelerated convergence of orbital‐based methods for this smooth reference Hamiltonian is shown by a detailed partial wave analysis. Contracted orbital basis sets constructed from atomic natural orbitals are shown to be significantly better for the new Hamiltonian than standard basis sets of the same size. The short‐range component becomes the perturbation. The low‐order perturbation equations are solved variationally using basis sets of correlated Gaussian geminals. Variational energies and low‐order perturbation wave functions for the model system are shown to be in excellent agreement with highly accurate numerical solutions for that system. Approximations of the reference wave functions, described by fewer basis functions, are tested for use in the perturbation equations and shown to provide significant computational advantages with tolerable loss of accuracy. Lower bounds for the radius of convergence of the resulting perturbation expansions are estimated. The proposed method is capable of achieving sub‐μHartree accuracy for all systems considered here. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

A real-space perturbation theory for electronic correlation

We introduce a method for treating electronic correlation in which a correlation factor is optimized using Hylleraas variational perturbation theory. The factor is defined by a number of parameters which grows only linearly with system size. We test the theory on two-electron atoms using the shielded-nucleus Hamiltonian at zeroth order, obtaining −2.9035 a.u. for helium. The convergence of the method is investigated, and energies and intracule densities are compared with accurate variational results. The application of the theory to an N-electron problem with a Hartree–Fock Hamiltonian at zeroth order is discussed.     

On the use of UHF orbitals in ionization energy calculations

The UHF Hamiltonian and simple Löwdin-like annihilators are formulated in the second quantization formalism. The so formulated Hamiltonian was employed in many-body Rayleigh-Schrödinger perturbation theory to evaluate the corrections to the UHF orbital energies.     

Combination of perturbative and variational methods for calculating molecular spectra: calculation of the upsilon=3-5 CH stretch overtone spectrum of CHF3

Highly excited states of the CHF3 molecule belonging to the third, fourth, and fifth Fermi polyad are calculated using a combination of the Van Vleck perturbation theory and a variational treatment. The perturbation theory preconditions the Hamiltonian matrix by transforming away all couplings except those between nearly degenerate states. This transformation is implemented so that eigenvalues can be found with significantly smaller matrices than that which would be needed in the original normal mode representation. Even with preconditioning, at the energies as high as 3-5 quanta in the CH stretch, it is not possible to directly diagonalize the Hamiltonian matrix due to the large basis sets required. Iterative methods, particularly the block-Davidson method, are explored for finding the eigenvalues. The methods are compared and the advantages discussed.     

Second order explicitly correlated R12 theory revisited: a second quantization framework for treatment of the operators' partitionings

Second order R12 theory is presented and derived alternatively using the second quantized hole-particle formalism. We have shown that in order to ensure the strong orthogonality between the R12 and the conventional part of the wave function, the explicit use of projection operators can be easily avoided by an appropriate partitioning of the involved operators to parts which are fully describable within the computational orbital basis and complementary parts that involve imaginary orbitals from the complete orbital basis. Various Hamiltonian splittings are discussed and computationally investigated for a set of nine molecules and their atomization energies. If no generalized Brillouin condition is assumed, with all relevant partitionings the one-particle contribution arising in the explicitly correlated part of the first order wave function has to be considered and has a significant role when smaller atomic orbital basis sets are used. The most appropriate Hamiltonian splitting results if one follows the conventional perturbation theory for a general non-Hartree-Fock reference. Then, no couplings between the R12 part and the conventional part arise within the first order wave function. The computationally most favorable splitting when the whole complementary part of the Hamiltonian is treated as a perturbation fails badly. These conclusions also apply to MP2-F12 approaches with different correlation factors.     

On the characterization of three state conical intersections: a quasianalytic theory using a group homomorphism approach

In this work, degenerate perturbation theory through second order is used to characterize the vicinity of a three state conical intersection. This report extends our recent demonstration that it is possible to describe the branching space (in which the degeneracy is lifted linearly) and seam space (in which the degeneracy is preserved) in the vicinity of a two state conical intersection using second order perturbation theory. The general analysis developed here is based on a group homomorphism approach. Second order perturbation theory, in conjunction with high quality ab initio electronic structure data, produces an approximately diabatic Hamiltonian whose eigenenergies and eigenstates can accurately describe the three adiabatic potential energy surfaces, the interstate derivative couplings, and the branching and seam spaces in their full dimensionality. The application of this approach to the minimum energy three state conical intersection of the pyrazolyl radical demonstrates the potential of this method. A Hamiltonian comprised of the ten characteristic (linear) parameters and over 300 second order parameters is constructed to describe the branching space associated with a point of conical intersection. The second order parameters are determined using data at only 30 points. In the vicinity of the conical intersection the energy and derivative couplings are well reproduced and the singularity in the derivative coupling is analyzed.     

New relativistic ANO basis sets for transition metal atoms

New basis sets of the atomic natural orbital (ANO) type have been developed for the first, second, and third row transition metal atoms. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive and negative ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies, electron affinities, and excitation energies for all atoms and polarizabilities for spherically symmetric atoms. These calculations include spin-orbit coupling using a variation-perturbation approach. Computed ionization energies have an accuracy better than 0.2 eV in most cases. The accuracy of computed electron affinities is the same except in cases where the experimental values are smaller than 0.5 eV. Accurate results are obtained for the polarizabilities of atoms with spherical symmetry. Multiplet levels are presented for some of the third row transition metals.     

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.     

A second-quantization framework for the unified treatment of relativistic and nonrelativistic molecular perturbations by response theory

A formalism is presented for the calculation of relativistic corrections to molecular electronic energies and properties. After a discussion of the Dirac and Breit equations and their first-order Foldy-Wouthuysen [Phys. Rev. 78, 29 (1950)] transformation, we construct a second-quantization electronic Hamiltonian, valid for all values of the fine-structure constant alpha. The resulting alpha-dependent Hamiltonian is then used to set up a perturbation theory in orders of alpha(2), using the general framework of time-independent response theory, in the same manner as for geometrical and magnetic perturbations. Explicit expressions are given to second order in alpha(2) for the Hartree-Fock model. However, since all relativistic considerations are contained in the alpha-dependent Hamiltonian operator rather than in the wave function, the same approach may be used for other wave-function models, following the general procedure of response theory. In particular, by constructing a variational Lagrangian using the alpha-dependent electronic Hamiltonian, relativistic corrections can be calculated for nonvariational methods as well.     

Correlated one-particle method: numerical results

In a previous paper a correlated one-particle method was formulated, where the effective Hamiltonian was composed of the Fock operator and a correlation potential. The objective was to define a correlated one-particle theory that would give all properties that can be obtained from a one-particle theory. The Fock-space coupled-cluster method was used to construct the infinite-order correlation potential, which yields correct ionization potentials (IP's) and electron affinities (EA's) as the negative of the eigenvalues. The model, however, was largely independent of orbital choice. To exploit the degree of freedom of improving the orbitals, the Brillouin-Brueckner condition is imposed, which leads to an effective Brueckner Hamiltonian. To assess its numerical properties, the effective Brueckner Hamiltonian is approximated through second order in perturbation. Its eigenvalues are the negative of IP's and EA's correct through second order, and its eigenfunctions are second-order Brueckner orbitals. We also give expressions for its energy and density matrix. Different partitioning schemes of the Hamiltonian are used and the intruder state problem is discussed. The results for ionization potentials, electron affinities, dipole moments, energies, and potential curves are given for some sample molecules.     

Localized bond model for molecular energy to fourth order in perturbation theory

The localized bond model of Malrieu, Diner, and Claverie is extended to fourth order in perturbation theory. Single, double, triple, and quadruple replacements from the doubly occupied bonding reference function are included utilizing a symmetric form of diagrammatic perturbation theory. The fourth order theory derived executes on a computer as quickly as does the third order theory. Results are examined utilizing the Pariser–Parr–Pople and CNDO/2 model Hamiltonians, and are compared with third order results and with either exact results where they are known, or with a configuration interaction of all singles and doubles. The influence of the initial hybridization, localization, and bond polarization is discussed. In general, the fourth order corrections are of comparable size to third order. Improvement in results appears to be marginal in the Nesbet–Epstein scheme in passing to fourth order because of the oscillating nature of the series; for Moller–Plesset theory errors are approximately halved. The relative energies as a function of modest geometry change about minima is about the same at third order as it is at fourth for most cases examined.