Multi‐state multi‐reference Møller–Plesset second‐order perturbation theory for molecular calculations

This work presents multi‐state multi‐reference Møller–Plesset second‐order perturbation theory as a variant of multi‐reference perturbation theory to treat electron correlation in molecules. An effective Hamiltonian is constructed from the first‐order wave operator to treat several strongly interacting electronic states simultaneously. The wave operator is obtained by solving the generalized Bloch equation within the first‐order interaction space using a multi‐partitioning of the Hamiltonian based on multi‐reference Møller–Plesset second‐order perturbation theory. The corresponding zeroth‐order Hamiltonians are nondiagonal. To reduce the computational effort that arises from the nondiagonal generalized Fock operator, a selection procedure is used that divides the configurations of the first‐order interaction space into two sets based on the strength of the interaction with the reference space. In the weaker interacting set, only the projected diagonal part of the zeroth‐order Hamiltonian is taken into account. The justification of the approach is demonstrated in two examples: the mixing of valence Rydberg states in ethylene, and the avoided crossing of neutral and ionic potential curves in LiF. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Relativistic and electron‐correlation effects on magnetizabilities investigated by the Douglas‐Kroll‐Hess method and the second‐order Møller‐Plesset perturbation theory

Isotropic and anisotropic magnetizabilities for noble gas atoms and a series of singlet and triplet molecules were calculated using the second‐order Douglas‐Kroll‐Hess (DKH2) Hamiltonian containing the vector potential A and in part using second‐order generalized unrestricted Møller‐Plesset (GUMP2) theory. The DKH2 Hamiltonian was resolved into three parts (spin‐free terms, spin‐dependent terms, and magnetic perturbation terms), and the magnetizabilities were decomposed into diamagnetic and paramagnetic terms to investigate the relativistic and electron‐correlation effects in detail. For Ne, Kr, and Xe, the calculated magnetizabilities approached the experimental values, once relativistic and electron‐correlation effects were included. For the IF molecule, the magnetizability was strongly affected by the spin‐orbit interaction, and the total relativistic contribution amounted to 22%. For group 17, 16, 15, and 14 hydrides, the calculated relativistic effects were small (less than 3%), and trends were observed in relativistic and electron‐correlation effects across groups and periods. The magnetizability anisotropies of triplet molecules were generally larger than those of similar singlet molecules. The so‐called relativistic‐correlation interference for the magnetizabilities computed using the relativistic GUMP2 method can be neglected for the molecules evaluated, with exception of triplet SbH. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Construction of Giant‐Spin Hamiltonians from Many‐Spin Hamiltonians by Third‐Order Perturbation Theory and Application to an Fe3Cr Single‐Molecule Magnet

A general giant‐spin Hamiltonian (GSH) describing an effective spin multiplet of an exchange‐coupled metal cluster with dominant Heisenberg interactions was derived from a many‐spin Hamiltonian (MSH) by treating anisotropic interactions at the third order of perturbation theory. Going beyond the existing second‐order perturbation treatment allows irreducible tensor operators of rank six (or corresponding Stevens operator equivalents) in the GSH to be obtained. Such terms were found to be of crucial importance for the fitting of high‐field EPR spectra of a number of single‐molecule magnets (SMMs). Also, recent magnetization measurements on trigonal and tetragonal SMMs have found the inclusion of such high‐rank axial and transverse terms to be necessary to account for experimental data in terms of giant‐spin models. While mixing of spin multiplets by local zero‐field splitting interactions was identified as the major origin of these contributions to the GSH, a direct and efficient microscopic explanation had been lacking. The third‐order approach developed in this work is used to illustrate the mapping of an MSH onto a GSH for an trigonal Fe₃Cr complex that was recently investigated by high‐field EPR spectroscopy. Comparisons between MSH and GSH consider the simulation of EPR data with both Hamiltonians, as well as locations of diabolical points (conical intersections) in magnetic‐field space. The results question the ability of present high‐field EPR techniques to determine high‐rank zero‐field splitting terms uniquely, and lead to a revision of the experimental GSH parameters of the Fe₃Cr SMM. Indeed, a bidirectional mapping between MSH and GSH effectively constrains the number of free parameters in the GSH. This notion may in the future facilitate spectral fitting for highly symmetric SMMs.     

The quasiparticle concept in vibrational–electronic problems in molecules. I. Partitioning of the vibrational–electronic Hamiltonian

The partitioning of the vibrational–electronic Hamiltonian is presented. This partitioning is based on a new quasiparticle transformation that is constructed in such a way that the adiabatic approximation is included into the unperturbed Hamiltonian; nonadiabacity, anharmonicity, and electron correlation are treated as perturbations. We also present the second quantization treatment for bosons. The many body perturbation theory expansion for the vibrational–electronic Hamiltonian is suggested. A comparison of this approach is made with gradient techniques.     

Perturbative calculation of the Hartree–Fock interaction energy using orthogonalized orbitals

A many‐body perturbation theory based on the partitioning of the dimer Hamiltonian, formulated in an orthogonalized basis set, is used for the calculation of interaction energies at the Hartree–Fock (HF) level. Numerical results for the (HF)₂ and (H₂O)₂ systems in selected geometries are presented. The interaction‐energy components are compared with the results obtained from the standard supermolecular approach and the intermolecular perturbation theory based on the biorthogonal basis set. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 81–88, 1999     

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     

Green's function calculations using non-Hartree–Fock orbitals

The single-particle Green's function is used to generate a new zero-order Hamiltonian. The idea to generate a new zero order from the previous zero order by incorporating perturbative corrections up to certain order is attractive since it allows an iterative procedure to repeatedly improve the results by decreasing the perturbation. In particular, in those cases where the Hartree–Fock Hamiltonian is not a good approximation to the full Hamiltonian and where perturbation theory usually does not produce sufficiently accurate results, one might hope that such a repetitive procedure ultimately yields an improved zero order and accurate perturbative corrections from this newly generated zero order. Two such approaches are investigated: first, one in which the ω-independent part of the self-energy is fully incorporated in the zero order and, second, one in which the correlation energy is incorporated in a one-electron potential in an average way. Numerical calculations are reported.     

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     

Exact size consistency of multireference Møller–Plesset perturbation theory

Single‐reference closed‐shell Møller–Plesset perturbation theory is well known for its size consistency, a quality that is essential for consistent comparisons of calculations on molecules of different size. However, it is far from obvious whether this quality can be retained in the multireference case. In this work it is shown that an exactly size consistently generalization to multireference perturbation theory can be constructed. The central result is that the zeroth‐order Hamiltonian should be constructed using separate projection operators for each excitation level, i.e., it should contain no couplings between different excitation levels. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 549–558, 1999     

Geminal approach to vibronic models of superconductivity in crystals. I. The Jahn–Teller effect

Electronic geminals constructed as linear combinations of binary products of site functions are used to formulate a vibronic model of superconductivity in crystals that is based upon the approximation of independent correlated electron pairs obtained variationally from an electron‐pair Hamiltonian and the Jahn–Teller effect. The cyclic symmetry of the system is taken into account and the geminals are sorted into doubly degenerate pairs. The Herzberg–Teller expansion of the pair Hamiltonian in terms of vibrational modes leads directly to the Jahn–Teller effect. A contact transformation of the vibronic Hamiltonian containing only linear terms lowers the energy of the system by a second‐order term associated with the Jahn–Teller stabilization energy. A possible model for superconductivity in solids is proposed on the basis of the Jahn–Teller effect. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

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.     

Superconducting first‐order pairing: Finite‐temperature simulations

A recently developed first‐order mechanism for superconducting pairing has been extended from T = 0 K to finite temperatures. On the basis of quantum statistical considerations, we have suggested a direct pairing interaction that does not necessarily involve second‐order elements, such as the electron–phonon coupling or specific magnetic interactions submitted by spin fluctuations. The driving force for the (energy‐driven) first‐order pairing is an attenuation of the destabilizing influence of the Pauli antisymmetry principle (PAP). Only the moves of unpaired fermions are controlled by the PAP, while the moves of superconducting Cooper pairs are not. The quantum statistics of Cooper pairs is of a mixed type, as it combines fermionic on‐site and bosonic intersite properties. The strong correlation between the strength of PAP constraints and system topology in combination with the electron number has been discussed for some larger clusters. Detailed finite‐temperature simulations on first‐order pairing have been performed for four‐center–four‐electron clusters with different topologies. A canonical ensemble statistics has been employed to derive the electronic energy, the electronic configuration entropy, and the free energy of paired and unpaired states in thermal equilibrium. The simulations show that pairing can be caused by either the electronic energy or the electronic configuration entropy. The coexistence of two different sets of quantum particles in paired states (i.e., the Cooper pairs and the unpaired electrons) can lead to an enhanced configuration entropy. In this context, we discuss the possibility of an entropy‐driven high‐temperature superconductor emerging from a low‐temperature unpaired state. The charge and spin degrees of freedom of the four‐center–four‐electron systems have been studied with the help of the charge and spin fluctuations. The spin fluctuations are helpful in judging the validity of pairing theories based on magnetic interactions. The charge fluctuations are a measure for the carrier delocalization in unpaired and paired states. The well‐known proximity between Jahn–Teller activity and superconductivity is analyzed in the zero‐temperature limit. It is demonstrated that both processes compete in their ability to reduce PAP constraints. All theoretical results have been derived within the framework of the simple Hubbard Hamiltonian. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Spin–orbit coupling of spin‐frustrated systems

We have investigated the effects of spin–orbit (SO) interactions on noncollinear molecular magnetism by combining the classical Dzyaloshinsky–Moriya (DM) model and ab initio generalized spin orbital (GSO) method. We have derived an estimation scheme of the magnetic anisotropy energy (MAE) and the Dzyaloshinsky vector based on the SO first‐order perturbation theory (SOPT1) for GSO Hartree–Fock (GHF) solutions. We found that the fundamental results of GHF‐SOPT1 method can be reproduced by diagonalizing the core Hamiltonian plus SO terms, and that the spin topologies of odd‐ring systems can be determined by the topological indices of the singly occupied molecular orbitals. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Spin‐orbit ZORA and four‐component Dirac–Coulomb estimation of relativistic corrections to isotropic nuclear shieldings and chemical shifts of noble gas dimers

Hartree–Fock and density functional theory with the hybrid B3LYP and general gradient KT2 exchange‐correlation functionals were used for nonrelativistic and relativistic nuclear magnetic shielding calculations of helium, neon, argon, krypton, and xenon dimers and free atoms. Relativistic corrections were calculated with the scalar and spin‐orbit zeroth‐order regular approximation Hamiltonian in combination with the large Slater‐type basis set QZ4P as well as with the four‐component Dirac–Coulomb Hamiltonian using Dyall's acv4z basis sets. The relativistic corrections to the nuclear magnetic shieldings and chemical shifts are combined with nonrelativistic coupled cluster singles and doubles with noniterative triple excitations [CCSD(T)] calculations using the very large polarization‐consistent basis sets aug‐pcSseg‐4 for He, Ne and Ar, aug‐pcSseg‐3 for Kr, and the AQZP basis set for Xe. For the dimers also, zero‐point vibrational (ZPV) corrections are obtained at the CCSD(T) level with the same basis sets were added. Best estimates of the dimer chemical shifts are generated from these nuclear magnetic shieldings and the relative importance of electron correlation, ZPV, and relativistic corrections for the shieldings and chemical shifts is analyzed. © 2015 Wiley Periodicals, Inc.     

A two‐scale approach to electron correlation in multiconfigurational perturbation theory

We present a new approach for the calculation of dynamic electron correlation effects in large molecular systems using multiconfigurational second‐order perturbation theory (CASPT2). The method is restricted to cases where partitioning of the molecular system into an active site and an environment is meaningful. Only dynamic correlation effects derived from orbitals extending over the active site are included at the CASPT2 level of theory, whereas the correlation effects of the environment are retrieved at lower computational costs. For sufficiently large systems, the small errors introduced by this approximation are contrasted by the substantial savings in both storage and computational demands compared to the full CASPT2 calculation. Provided that static correlation effects are correctly taken into account for the whole system, the proposed scheme represent a hierarchical approach to the electron correlation problem, where two molecular scales are treated each by means of the most suitable level of theory. © 2014 Wiley Periodicals, Inc.     

Ab initio correlation corrections to the Hartree-Fock quasi band-structure of periodic systems employing Wannier-type orbitals

A size-consistent ab initio formalism to calculate correlation corrections to ionization potentials as well as electron affinities of periodic systems is presented. Our approach is based on a Hartree-Fock scheme which directly yields local orbitals without any a posteriori localization step. The use of local orbitals implies non-zero off-diagonal matrix elements of the Fock operator, which are treated as an additional perturbation and give rise to localization diagrams. Based on the obtained local orbitals, an effective Bloch Hamiltonian is constructed to second order of perturbation theory with all third-order localization diagrams included. In addition, the summation of certain classes of diagrams up to infinite order in the off-diagonal Fock elements as well as the Epstein-Nesbet partitioning of the full Hamiltonian are discussed. The problem of intruder states, frequently encountered in many-body perturbation theory, is dealt with by employing the theory of intermediate Hamiltonians. As model systems we have chosen cyclic periodic structures up to an oligoethylene ring in double-zeta basis; however, the theory presented here straightforwardly carries over to infinite periodic systems. Received: 30 April 1998 / Accepted: 27 July 1998 /  Published online: 7 October 1998     

Heisenberg behavior of some carbon‐beryllium compounds: How well truncated‐CI approaches work

This works tries to establish the performance of truncated CI calculations on the evaluation of magnetic coupling parameters with respect to available FCI estimates on a set of carbon‐beryllium clusters. First‐, second‐ and third‐neighbor magnetic coupling constants have been evaluated and many body effective parameters as the cyclic terms. They result from the fitting of the low‐lying states to the eigenvalues of an extended Heisenberg Hamiltonian, involving not only two‐body isotropic terms but also cyclic terms. SDCI and DDCI calculations have been carried out and their performance compared with FCI ones. The impact of the basis set choice and size‐consistency errors have been explored. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Studies in perturbation theory XIII. Treatment of constants of motion in resolvent method,partitioning technique,and perturbation theory

After a brief survey of some basic concepts in the theory of linear spaces, the eigenvalue problem is formulated in the resolvent technique based on the introduction of a reference function φ and a complex variable ?. This leads to a series of fundamental concepts including the trial wave function, the inhomogeneous equation, and finally the transition and expectation values of the Hamiltonian, of which the former renders a “bracketing function” for the energy. In order to avoid the explicit limiting procedures in this approach, the eigenvalue problem is then reformulated in terms of the partitioning technique which, in turn, leads to a closed form of infinite-order perturbation theory. The eigenvalue problem is greatly simplified if the Hamiltonian H has a constant of motion Λ or has symmetry properties characterized by the group G = {g}, and the question is now how these simplifications can be incorporated into the partitioning technique and into perturbation theory. In both cases, there exists a set of projection operators {Q_k} which lead to a splitting of the Hilbert space into subspaces which have virtually nothing to do with each other. It is shown that, in the partitioning technique, it is sufficient to consider one of these subspaces at a time, and the results are then generalized to perturbation theory. It turns out that the finite-order expansions are no longer unique, and the commutation rules connecting the various forms are derived. The infinite-order results are finally presented in such a form that they are later suitable for the evaluation of upper and lower bounds to the energy eigenvalues.     

Two‐component calculations of spin–orbit effects for a van der Waals molecule Rn2

Electronic structures of the weakly bound Rn₂ were calculated by the two‐component Møller–Plesset second‐order perturbation and coupled‐cluster methods with relativistic effective core potentials including spin–orbit operators. The calculated spin–orbit effects are small, but depend strongly on the size of basis sets and the amount of electron correlations. Magnitudes of spin–orbit effects on D_e (0.7–3.0 meV) and R_e (−0.4∼−2.2 Å) of Rn₂ are comparable to previously reported values based on configuration interaction calculations. A two‐component approach seems to be a promising tool to investigate spin–orbit effects for the weak‐bonded systems containing heavy elements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 139–143, 1999