Applications of the MBPT in the localized representation

Diagrammatic formulation of the MBPT is applied when the occupied and the virtual canonical orbitals are separately localized by unitary transformations. In this localized representation, due to the off-diagonal Fock matrix elements, the perturbation operator contains extra terms generating the so-called localization corrections. These corrections enter the perturbation energy in third and higher orders. Their magnitude depends on the type of localization, but they represent only a small fraction of the canonical corrections. The calculation of the localization corrections, however, does not need a significant amount of extra computer time. It is shown that by introducing an “order of neighborhood” local and nonlocal effects of the electron correlation can be separated and the contribution of the nonlocal effects can be neglected to a good approximation. Ab initio calculations have been carried out for the normal saturated hydrocarbons: C_2n+1H_4n+4 and for the all-trans conjugated polyenes C_2n+2H_2n+4. As to the ratio of the local and nonlocal corrections, it is shown that there is only a quantitative difference for these two kinds of systems (strongly or weakly localizable). Neglecting nonlocal effects, considerable amount of computer time can be saved.     

Localized orbitals for the description of molecular interaction

The intermolecular interaction between the molecules CH₂O and NH₃ was investigated by the supermolecule method. The interaction energies were first calculated at the ab initio SCF level, and the electron correlation was included via second-order Møller-Plesset perturbation theory (MP 2). The basis set superposition error (BSSE ) was taken into account by the counter-poise (CP ) method. The occupied and the virtual canonical molecular orbitals (CMOS ) of the supermolecule were separately localized by the Boys' procedure. The correlation correction was calculated by the many-body perturbation theory (MBPT ) in the localized representation. Contributions of the third- and fourth-order localized diagrams were added to those of the second-order canonical diagram. This procedure gives a correction nearly equivalent to that of MP 2. The possibility to separate LMO contributions responsible for the dispersion interaction was investigated.     

Ab-Initio MRD-CI calculations for breaking a chemical bond in a molecule in a crystal or other solid environment. II.H3C—NO2 decomposition of nitromethane in a nitromethane crystal with voids

Recently we extended our strategy for MRD-CI (multireference double excitation-configuration interaction) calculations based on localized/local orbitals and an “effective” CI Hamiltonian for molecular decompositions of large molecules to breaking a chemical bond in a molecule in a crystal or other solid environment. Our technique involves solving a quantum chemical ab-initio SCF explicitly for a system of a reference molecule surrounded by a number of other molecules in the multipole environment of more distant neighbors. The resulting canonical molecular orbitals are then localized and the localized occupied and virtual orbitals in the region of interest are included explicitly in the MRD-CI with the remainder of the occupied localized orbitals being folded into an “effective” CI Hamiltonian. The MRD-CI calculations are carried out for breaking a bond in the reference molecule. This method is completely general. The space treated explicitly quantum chemically and the surrounding space can have voids, defects, deformations, dislocations, impurities, dopants, edges and surfaces, boundaries, etc. We previously applied this procedure successfully to the H₃C? NO₂ bond dissociation of nitromethane in a nitromethane crystal with extensive testing of the number of molecules that have to be included explicitly in the SCF and how many molecules have to be represented by more distant multipoles. The results indicated that it took more energy to dissociate the H₃C? NO₂ bond when the nitromethane molecule was in the crystal than it did to dissociate that bond in the free nitromethane molecule. In this present study we have investigated the effect of voids (both in the nitromethane molecules treated explicitly in the SCF and those in the environment represented by multipoles) on the calculated H₃C? NO₂ bond dissociation energies.     

A local second-order Møller-Plesset method with localized orbitals: a parallelized efficient electron correlation method

Using orthogonal localized occupied orbitals we have developed and implemented a parallelized local second-order M?ller-Plesset (MP2) method based on the idea developed by Head-Gordon and co-workers. A subset of nonorthogonal correlation functions (the orbital domain) was assigned to each of the localized occupied orbitals using a distance criterion and excitations from localized occupied orbitals that were arranged into subsets. The correlation energy was estimated using a partial diagonalization and an iterative efficient method for solving large-scale linear equations. Some illustrative calculations are provided for molecules with up to 1484 Cartesian basis sets. The orbital domain sizes were found to be independent of the molecular size, and the present local MP2 method covered about 98%-99% of the correlation energy of the conventional canonical MP2 method.     

Localization of virtual orbitals

The application of the MBPT in the localized representation requires that both the occupied and the virtual orbitals obtained by the canonical HF equation should be localized. The localization of the occupied orbitals is straightforward in general by any localization method. It is shown that by using Boys' method the localized virtual orbitals are spatially well separated and transferable not only in minimal basis sets.     

Ab initio MRD-CI calculations for breaking a chemical bond in a molecule in a crystal or other solid environment. III. Me2N?NO2 decomposition of dimethylnitramine in a large crystalline environment

Recently we extended our strategy for MRD-CI (multireference double excitation-configuration interaction) calculations, based on localized/local orbitals and an “effective” CI Hamiltonian, for molecular decompositions of large molecules to breaking a chemical bond in a molecule in a crystalline or other solid environment. Our technique begins with an explicit quantum chemical SCF calculation for a reference molecule surrounded by a number of other molecules in the multipole environment of more distant neighbors. The resulting canonical molecular orbitals are then localized, and the localized occupied and virtual orbitals in the region of interest are included explicitly in the MRD-CI with the remainder of the occupied localized orbitals being folded into an “effective” CI Hamiltonian. The MRD-CI calculations are then carried out for breaking a bond in the reference molecule. This method is completely general in that the space treated explicitly, as well as the surrounding space, may contain voids, defects, deformations, dislocations, impurities, dopants, edges and surfaces, boundaries, etc. Dimethylnitramine is the smallest prototype of the energetic R₂N—NO₂ nitramines, such as the 6-member ring RDX or the 8-member ring HMX. Decomposition of energetic compounds is initiated in the solid by a breaking of the target bond. Thus, it is crucial to know the difference in energy between breaking a bond in an isolated energetic molecule versus in the molecule in a solid. In the present study, we have carried out MRD-CI calculations for the Me₂N—NO₂ dissociation of dimethylnitramine in a dimethylnitramine crystal. The cases we investigated were one dimethylnitramine molecule (surrounded by 53 and 685 neighboring dimethylnitramine molecules represented by multipoles), three dimethylnitramine molecules, and three dimethylnitramine molecules (surrounded by 683 neighbors). All multipoles were cumulative atomic multipoles up through quadrupoles. The MRD-CI calculations on dimethylnitramine required large numbers of reference configurations from which were allowed all single and double excitations.     

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.     

Orbital Energy Order and the Through Bond Interaction in Homonuclear Diatomic Molecules

The homonuclear diatomic molecules are the simplest systems having both the σ framework and the lone pair orbitals n_a and _b for investigating their through space and through bond interaction. The striking orbital energy order n_g～ n_a+ n_b > n_n ～ n_a - n_b has been accounted for by the through bond interaction. However, when the p-content in the lone pair orbitals n_a and n_b decreases, one may have the reverse orbital energy order: n_g < n_g. A reverse orbital energy order has been found in F₂ and Cl₂, whose n_a and n_b are almost pure s-type atomic orbitals. The reverse order also occurs in molecule N₂ when the internuclear distance is larger man 1.5 Å. It is also found that the detail through space and through bond interaction and the eventual orbital energy order for n_g and n_u can be accounted for by the Fock operator within the localized molecular orbital space.     

Symmetry simplifications of space types in configuration interaction induced by orbital degeneracy

Symmetry simplifications are introduced in configuration interaction (CI ) by reducing the number of symmetry-allowed space types if there is degeneracy in some of the molecular orbitals by constructing the unique space types. A new symmetry group which we call the configuration symmetry group is defined and is shown to be expressible as a generalized wreath product group. Generating functions are derived for enumerating the equivalence classes of space types. A double coset method is expounded which constructs the representatives of all equivalence classes of space types using the cycle index of generalized wreath product and the double cosets of label subgroup with generalized wreath product in the symmetric group S_n, if n is twice the number of occupied and virtual orbitals. Method is illustrated with CI using the localized orbitals of polyenes, CI in benzene, and atomic CI for several reference states.     

Counterpoise corrected calculations at the correlated level: A simplified method using LMOs

Ab initio calculations have been performed in order to investigate the counterpoise corrections, especially at the correlated level for molecular interactions. It is pointed out that, when using a localized representation, the calculations using the MBPT/MP2 method can be simplified. The H₂O + H₂O system was studied, with the use of 6-31 G* basis set. The method allows one to determine the intramolecular correlation components in a simplified way. Boys' localization procedure was applied throughout, both in the occupied as well as in the virtual spaces.     

Finite-field many-body perturbation theory IV. Basis set optimization in MBPT calculations of molecular properties. Molecular quadrupole moments

The choice of truncated basis sets and their optimization for MBPT calculations of molecular properties are discussed. It is pointed out that computing the correlation corrections to some kth order property by using the MBPT approach requires the knowledge of accurate perturbed orbitals through the kth order. Hence, it is argued that the basis set functions can be optimized with respect to the perturbed energies calculated within the coupled Hartree-Fock method. The proposed procedure is illustrated by MBPT calculations of quadrupole moments of H₂ and FH. Additionally, also some estimates of the quadrupole polarizability tensor components for these molecules are obtained.     

A numerically stable procedure for calculating Møller-Plesset energy derivatives,derived using the theory of Lagrangians

Summary When Møller-Plesset energy derivatives are determined in the canonical Hartree-Fock basis, singularities or instabilities may arise due to degeneracies among the occupied or unoccupied orbitals. If a non-canonical basis is used these singularities disappear. Numerically stable expressions are presented for the molecular gradient and Hessian of the second-order Møller-Plesset energy, obtained by differentiating a fully variational Lagrangian of the energy constructed in a non-canonical representation. By using a non-canonical representation, singularities and instabilities are avoided, and the variational property of the Lagrangian ensures that Wigner's 2n + 1 rule is satisfied for the orbital derivatives and that the multipliers satisfy the stronger 2n + 2 rule. It is shown that the most expensive step in the calculation of the Hessian scales as Mn ₄o, where M is the number of independent Cartesian distortions, n the total number of orbitals, and o the number of occupied orbitals.     

A simple extension of the external magnasco-perico localization procedure to the virtual MO-Space

The external localization procedure of Magnasco and Perico is extended to the unoccupied molecular orbitals of the Fock-operator. The formal correspondence between bonding orbitals and localized antibonding MOs is demonstrated. Localized occupied and virtual one-electron functions are calculated within a semiempirical INDO-Hamiltonian and are analyzed; the externally localized occupied MOs are compared with energy localized orbitals computed by the Edmiston and Ruedenberg procedure. Various applications of the fully localized (occupied and virtual) MO set are discussed.     

Linear scaling local correlation approach for solving the coupled cluster equations of large systems

A linear scaling local correlation approach is proposed for approximately solving the coupled cluster doubles (CCD) equations of large systems in a basis of orthogonal localized molecular orbitals (LMOs). By restricting double excitations from spatially close occupied LMOs into their associated virtual LMOs, the number of significant excitation amplitudes scales only linearly with molecular size in large molecules. Significant amplitudes are obtained to a very good approximation by solving the CCD equations of various subsystems, each of which is made up of a cluster associated with the orbital indices of a subset of significant amplitudes and the local environmental domain of the cluster. The combined effect of these two approximations leads to a linear scaling algorithm for large systems. By using typical thresholds, which are designed to target an energy accuracy, our numerical calculations for a wide range of molecules using the 6-31G or 6-31G* basis set demonstrate that the present local correlation approach recovers more than 98.5% of the conventional CCD correlation energy.     

Energies of Association of Carbenium and Silylium Cations with Oxygen-Containing Molecules

The energies of association of Me_nH_3-n X⁺ cations (n = 0-3) with molecules of oxygen-containing bases (water, methanol, dimethyl ether) were calculated by the B3LYP/6-31G(d) method. These energies decrease with an increase in the number n of methyl groups at the X atom; this trend is more pronounced for X = C, which is due to a lower degree of charge transfer in the cation from the occupied orbitals of the methyl group to the vacant orbital of X in the case of X = Si. An increase in the association energy with an increase in the number of methyl substituents at the oxygen atom is due to an increase in the energy level of lone electron pairs of oxygen upon methyl substitution. As a result, the energy gap between the electronic levels of the unoccupied orbital of the cation and orbitals of the oxygen lone electron pairs becomes narrower, which makes the interaction between the unoccupied and occupied orbitals more efficient.__________Translated from Zhurnal Obshchei Khimii, Vol. 75, No. 5, 2005, pp. 754–756.Original Russian Text Copyright © 2005 by Ignat’ev, Kochina.     

The divide-expand-consolidate family of coupled cluster methods: numerical illustrations using second order Møller-Plesset perturbation theory

Previously, we have introduced the linear scaling coupled cluster (CC) divide-expand-consolidate (DEC) method, using an occupied space partitioning of the standard correlation energy. In this article, we show that the correlation energy may alternatively be expressed using a virtual space partitioning, and that the Lagrangian correlation energy may be partitioned using elements from both the occupied and virtual partitioning schemes. The partitionings of the correlation energy leads to atomic site and pair interaction energies which are term-wise invariant with respect to an orthogonal transformation among the occupied or the virtual orbitals. Evaluating the atomic site and pair interaction energies using local orbitals leads to a linear scaling algorithm and a distinction between Coulomb hole and dispersion energy contributions to the correlation energy. Further, a detailed error analysis is performed illustrating the error control imposed on all components of the energy by the chosen energy threshold. This error control is ultimately used to show how to reduce the computational cost for evaluating dispersion energy contributions in DEC.     

Modifications of virtual orbitals in the limited CI calculations for electron-rich molecules

The efficiency of modified virtual orbitals (MVO) of ionic type and of approximate orthogonalized natural orbitals (ONO) in the CI-SD calculations was studied for O₃ and SO₂ molecules and compared with the commonly used canonical virtual orbitals (CVOs). The systems studied represent a class of electron-rich molecules, in which the number of valence electron pairs exceeds substantially the number of formal chemical bonds. We found that the modified orbitals of the types studied appear to be less effective for these systems than in the similar calculations for the AH_n type molecules. Physical reasons for this difference were discussed. The evolution of spatial properties of virtual orbitals within the modification process was analyzed. © 1996 John Wiley & Sons, Inc.     

A modified PCILO method

The perturbative configuration interaction using strictly localized molecular orbitals, called the modified PCILO method, for which the use of the Rayleigh-Schrödinger many-body perturbation theory with the Moller-Plesset Hamiltonian partitioning is characteristic, has been proposed in this communication. On the CNDO/2 and INDO levels of Hamiltonian approximations strictly localized molecular orbitals have been constructed by solving modified Roothaan equations. From the zero and second order energy interatomic distances and harmonic force constants for some diatomic molecules have been calculated. The linear dependence of the correlation energy on the number of valence electrons in the series of the molecules CH₄, CH₃F, CH₂F₂, CHF₃ and CF₄ is perfect.     

Local orbitals by minimizing powers of the orbital variance

It is demonstrated that a set of local orthonormal Hartree-Fock (HF) molecular orbitals can be obtained for both the occupied and virtual orbital spaces by minimizing powers of the orbital variance using the trust-region algorithm. For a power exponent equal to one, the Boys localization function is obtained. For increasing power exponents, the penalty for delocalized orbitals is increased and smaller maximum orbital spreads are encountered. Calculations on superbenzene, C(60), and a fragment of the titin protein show that for a power exponent equal to one, delocalized outlier orbitals may be encountered. These disappear when the exponent is larger than one. For a small penalty, the occupied orbitals are more local than the virtual ones. When the penalty is increased, the locality of the occupied and virtual orbitals becomes similar. In fact, when increasing the cardinal number for Dunning's correlation consistent basis sets, it is seen that for larger penalties, the virtual orbitals become more local than the occupied ones. We also show that the local virtual HF orbitals are significantly more local than the redundant projected atomic orbitals, which often have been used to span the virtual orbital space in local correlated wave function calculations. Our local molecular orbitals thus appear to be a good candidate for local correlation methods.