Chapter 4 On the Selection of Domains and Orbital Pairs in Local Correlation Treatments
Section snippets
1. Introduction
Local correlation treatments, as originally proposed by Pulay [1] and first implemented by Pulay and Saebø[2], [3], [4], [5], [6] employ local orbital spaces to restrict the number of excited configurations in the wave function. The occupied orbitals are localized using standard procedures (e.g. Boys [7] or Pipek-Mezey [8] localization), while the virtual space is spanned by non-orthogonal projected atomic orbitals (PAOs). By fully exploiting the locality, our group has achieved linear scaling
Method
In the following sections, we summarize the criteria used to define the local orbital spaces and domains as implemented in the MOLPRO package of ab initio programs [35]. The keywords for thresholds are the same as used in this program.
Standard domains
For each correlated occupied orbital an orbital domain [i] is generated according to the procedure proposed by Boughton and Pulay [17]. Some modifications were made in our program, which are outlined in the following. The purpose of the method is to include all PAOs in the domain, which arise from AOs (basis functions), that significantly contribute to the considered LMO. The PAOs arising from AOs at a given atom are considered as a group. Therefore, the first step is to select a set of
Pair classes
The orbital pairs (ij) are classified according to the closest distance R(ij) between atoms in the primary domains [i] and [j]. This classification is independent of domain extensions. Furthermore, only atoms in the primary domains are considered for the pair classification if the atomic Löwdin charge is larger than CHGMIN_PAIRS (default value 0.2). This criterion was introduced in order to reduce the dependence of the pair selection on localization tails. The strong pairs (0≤R(ij)<RCLOSE)
Triple excitations
As for the double excitations, the triple excitations are restricted to domains (cf. Section 3.3) and by a triples list (ijk) of LMOs, so that the total number of triple excitations scales linearly with molecular size. As discussed in detail in Ref. [14], the triples list (ijk) is defined by the condition that the pairs (ij), (ik), or (jk) must be either strong or close pairs. Additionally, at least one of these pairs must be strong. Second, the close pair amplitudes as determined in the
Dependence of the correlation energy on the domain approximation
In order to demonstrate the effect of the domain approximation on correlation and reaction energies, we have studied 52 chemical reactions involving 59 molecules [38]. Here we present only a representative subset of the results. The geometries of all molecules have been optimized at the MP2/aug-cc-pV(T+d)Z level. Table 1, Table 2, Table 3 show the computed correlation energies of 22 molecules for LMP2, LCCSD, and LCCSD(T0). In each case, the results are compared to the full non-local
Conclusions
We have reviewed the approximations made in the linear scaling local correlation methods developed in our group. The accuracy of local coupled-cluster calculations depends on (i) the domain sizes and (ii) the definition of the strong-pair list included in the LCCSD. These approximations can be controlled using connectivity or distance criteria. The domain sizes can be determined by a single parameter (IEXT or REXT), and we have shown that domain extensions lead to rapid and systematic
Acknowledgements
This work has been supported by the Deutsche Forschungsgemeinschaft and the Fonds der Chemischen Industrie. The authors thank Dr. A. Schäfer for providing the experimental data as well as zero-point and thermodynamic corrections for the reaction enthalpies.
References (45)
Localizability of dynamic electron correlation
Chem. Phys. Lett.
(1983)- et al.
Local configuration-interaction: an efficient approach for larger molecules
Chem. Phys. Lett.
(1985) - et al.
Multipole approximation of distant pair energies in local MP2 calculations
Chem. Phys. Lett.
(1998) - et al.
Local perturbative triples correction (T) with linear cost scaling
Chem. Phys. Lett.
(2000) - et al.
Impact of local approximations on MP2 vibrational frequencies
Spectrochim. Acta A
(1999) - et al.
Local correlation in coupled cluster calculations of molecular response properties
Chem. Phys. Lett.
(2004) - et al.
Orbital-invariant formulation and second-order gradient evaluation in Møller–Plesset perturbation theory
Theor. Chim. Acta
(1986) - et al.
Fourth-order Møller–Plesset perturbation theory in the local correlation treatment. I. Method
J. Chem. Phys.
(1987) - et al.
The local correlation treatment. II. Implementation and tests
J. Chem. Phys.
(1988) - et al.
Local treatment of electron correlation
Annu. Rev. Phys. Chem.
(1993)
A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions
J. Chem. Phys.
Local treatment of electron correlation in coupled cluster theory
J. Chem. Phys.
Low-order scaling local electron correlation methods. I. Linear scaling local MP2
J. Chem. Phys.
Low-order scaling local correlation methods. II: Splitting the Coulomb operator in linear scaling local second-order Møller–Plesset perturbation theory
J. Chem. Phys.
Low-order scaling local electron correlation methods. III. Linear scaling local perturbative triples correction (T)
J. Chem. Phys.
Low-order scaling local electron correlation methods. IV. Linear scaling local coupled-cluster (LCCSD)
J. Chem. Phys.
Linear scaling local connected triples beyond local (T): Local CCSDT-1b with O(N) scaling
J. Chem. Phys.
Comparison of the Boys and Pipek–Mezey localizations in the local correlation approach and automatic virtual basis selection
J. Comput. Chem.
Analytical energy gradients for local second-order Møller–Plesset perturbation theory
J. Chem. Phys.
Analytical energy gradients for local coupled-cluster methods
Phys. Chem. Chem. Phys.
Analytical energy gradients for local second-order Møller–Plesset perturbation theory using density fitting approximations
J. Chem. Phys.
Cited by (88)
Shape resonance induced electron attachment to cytosine: The effect of aqueous media
2023, Journal of Chemical PhysicsBimolecular sinks of Criegee intermediates derived from hydrofluoroolefins - a computational analysis
2023, Environmental Science: AtmospheresProjector-free capped-fragment scheme within density functional embedding theory for covalent and ionic compounds
2021, Journal of Chemical Theory and ComputationAn efficient Fock space multi-reference coupled cluster method based on natural orbitals: Theory, implementation, and benchmark
2021, Journal of Chemical Physics