Diagrammatic Simplification of Linearized Coupled Cluster Theory

TL;DR

This paper introduces a simplified and more accurate method for calculating energies in chemically complex systems using a modified version of coupled cluster theory.

Contribution

The novel contribution is the development of a size-consistent LinLCCD method that avoids divergences by retaining only linear ladder diagrams.

Findings

LinLCCD is robust in statically correlated systems and achieves O(nocc4nvir2) scaling with a hole–hole approximation.

LinLCCD outperforms full CCD and CCSD for noncovalent interactions in small-to-medium molecules.

The method provides accurate energy differences and can be adapted to include more dynamical correlation.

Abstract

Linearized Coupled Cluster Doubles (LinCCD) often provides
near-singular
energies in small-gap systems that exhibit static correlation. This
has been attributed to the lack of quadratic T̂
2
2 terms that
typically balance out small energy denominators in the CCD amplitude
equations. Herein, I show that exchange contributions to ring and
crossed-ring contractions (not small denominators per se) cause the divergent behavior of LinCC­(S)­D approaches. Rather than
omitting exchange terms, I recommend a regular and size-consistent
method that retains only linear ladder diagrams. As LinCCD and configuration
interaction doubles (CID) equations are isomorphic, this also implies
that simplification (rather than quadratic extensions) of CID amplitude
equations can lead to a size-consistent theory. Linearized ladder
CCD (LinLCCD) is robust in statically correlated systems and can be
made
O(nocc4nvir2)
with a hole–hole approximation.
The results presented here show that LinLCCD and its hole–hole
approximation can accurately capture energy differences, even outperforming
full CCD and CCSD for noncovalent interactions in small-to-medium
sized molecules, setting the stage for further adaptations of these
approaches that incorporate more dynamical correlation.

Figures

Diagrammatic representation of the linearized CCD equations. The “driver” terms are displayed in black, “ladder” terms in blue, and “ring and crossed-ring” in red. When permutation operations are shown, the associated diagram corresponds to the contraction displayed in parentheses.

Dissociation curves for (a) hydrogen molecule in the minimal STO-3G basis set, (b) Hydrogen fluoride molecule in the aug-cc-pVQZ basis, , and (c) H6 in the cc-pVDZ basis. Dashed lines of like-color indicate the dissociation limit for a particular method. Dissociation limits were estimated at R = 106 Å, except in the case of H6 where R = 103 Å was used instead due to convergence difficulties.

Error statistics for noncovalent interaction energies extrapolated to the complete basis set limit for (a) the A24 set of small dimers and (b) the S22 data set of small to medium sized dimers. The inset numbers indicate the mean absolute error in kcal/mol.

Systems that comprise the L7 data set along with their commonly employed acronyms.

Mean absolute errors (shown as colored bars) and maximum errors (shown as gray bars) for the lowest energy singlet/triplet gaps in Quest #8 as computed by various ΔCC and ΔMP2 methods. The bars are colored according to the category of approximation being applied, where Møller–Plesset perturbation theory is red, linearized CC methods are orange, and nonlinear CC approaches are green. The systems in question are shown in the inset.

Dissociation of CF3COCl along the C–Cl bond using ΔLinLCCD from reference configurations that represent the ground electronic state as well as n → π*, p → π*, and p → σ* transitions.

Tables

1. Interaction Energies for L7 Data Set (kcal/mol)

system	LinLCCD(hh)	BW-s2,	MP2,	CCSD(T0),
C2C2PD	–34.54	–33.32	–38.08	–20.93 ± 0.44
C3A	–25.06	–24.11	–27.09	–17.49 ± 0.31
C3GC	–41.54	–40.40	–45.37	–29.24 ± 0.91
GCGC	–16.50	–17.00	–18.99	–13.54 ± 0.27
GGG	–3.31	–3.63	–4.54	–2.08 ± 0.09
CBH	–9.40	–10.90	–11.83	–11.00 ± 0.17
PHE	–25.94	–25.73	–26.32	–25.46 ± 0.01
MAE	5.36	4.78	7.18

Full text

Introduction

1

Coupled cluster (CC) theory with double substitutions (CCD) is the simplest form of CC that captures electron correlation.? There are a host of advantages to linearized CCD methods (LinCCD) over full CCD, ?,? including reductions in memory demands, ease of spin-adapting the LinCCD wave function (albeit there is no rigorous wave function in linearized CC approximations),? and simpler physical interpretation of the equations. My research group is particularly interested in the Hermitian formulation that is offered by linearized CC methods, as this can be useful in developing excited-state theories.? Hermitian approaches are also quite powerful in the sense that they satisfy the generalized Hellmann–Feynman theorem, permitting simpler evaluation of forces (whereas left-eigenvectors are required in non-Hermitian CC approaches).?

However, linearized CC approaches often encounter near singularities in small-gap systems, affecting the performance of LinCCD away from equilibrium.? Small orbital-energy gaps are often a qualitative indicator of static correlation,? where the near-singular behavior of LinCCD can be further understood as a deficiency resulting from the lack of quadratic T̂ 2 ^2^ terms that fold in higher-order correlation effects necessary to describe static correlation. In other words, LinCCD lacks the implicit account for quadruple excitations that is found in the CCD amplitude equations, making it unable to counteract small energy denominators. Consequently, LinCCSD has been combined with Tikhonov regularization? as a means of sidestepping divergences.? Multireference LinCC approaches have also been developed to avoid divergences in systems that exhibit static correlation, albeit at increased cost. ?−? ? Furthermore, LinCC methods have been applied to capture additional correlation effects atop geminal reference states that naturally incorporate static correlation. ?−? ? ?

Beyond LinCCD, there are classes of CC approaches that attempt to correct errors within single-reference CC that arise due to static correlation. These “addition-by-subtraction” (ABS) CC methods take the seemingly paradoxical approach of removing components of the T̂ operator that are found to be particularly ill-behaved in the face of static correlation.? Perhaps the most well-known ABS-CC approach is pair CCD,? where only pair double substitution clusters are retained, leading to an approach that can describe single and double-bond dissociation. ?−? ? ? Alternatively, it is also possible to decouple the singlet- and triplet-paired amplitudes in CCD to achieve similarly well-behaved bond dissociation curves. ?,? Though, such singlet/triplet-pair couplings occur through a quadratic term in the CCD equations so the divergence of LinCCD must be attributable to other factors.

Another flavor of ABS-CC approach restricts the CCD equations to certain classes of diagrams. For instance, CCD with only ring diagrams is equivalent to the particle-hole random-phase approximation (ph-RPA) ?−? ? and is especially applicable in the case of the high-density homogeneous electron gas. Beyond single-reference approaches, diagrammatic resummations of the ring diagrams have been recently applied to extend ph-RPA to the multireference case.?

At the other end of the spectrum, CCD restricted to ladder diagrams (ladder-CCD) is formally equivalent to particle–particle RPA (pp-RPA) ?,? and is especially suitable in the limit of the low-density electron gas due to its explicit account of particle–particle correlations. Diagrammatic analysis of the ring and ladder CC equations has revealed that ring and ladder diagrams mainly describe long and short-ranged correlation effects, respectively.? These naturally imposed length scales have been leveraged in combinations of ladder- and ring-CCD via range-separation techniques, leading to promising methods for describing systems that do not fall into either extreme.?

In this work, I present an ABS linearized CCD approach that linearizes the ladder CCD amplitude equations. By removing the terms associated with ring and crossed-ring diagrams from LinCCD, the resultant linearized ladder CCD (LinLCCD) approach avoids the near-singularities encountered in these diagrams. Furthermore, the isomorphism between LinCCD and configuration interaction with double substitutions (CID) suggests that LinCCD equations are not size-consistent.? A lack of size-consistency implies that, for well-separated molecular fragments A and B, E _ A∪B _ = E _ A _ + E _ B _ is not satisfied by LinCCD/CID. While quadratic corrections have been added to CID to obtain CCD equations, ?−? ? ? ? ? revealing the role of T̂ 2 ^2^ terms in size-consistent approaches, I propose removing ring/crossed-ring terms from the LinCCD (or equivalently CID) equations as an alternative route to obtain a size-consistent, size-extensive, orbital invariant, and naturally regular method.

Theory

2

Throughout this work, I will denote occupied orbitals as i, j, k, l, ···, virtual orbitals as a, b, c, d, ···, and general unspecified orbitals as p, q, r, s, ···. The abbreviations n _ v _ and n _ o _ will be used for the number of virtual orbitals and the number of occupied orbitals, respectively. Einstein summation notation is used except in limited cases where the summation is explicitly written out.

Linearized Coupled Cluster Theory

2.1

LinCCD invokes the approximation that the usual exponential parameterization of the wave function

$$[eqn]$$

is Taylor-expanded through first order such that

$$[eqn]$$

By retaining only strongly connected diagrams, the LinCCD energy can be cast in terms of a Hermitian Hamiltonian

$$[eqn]$$

where

$$[eqn]$$

is the usual double-substitution operator and the subscript SC indicates that only the strongly connected diagrams are retained. Connected diagrams are defined as those whose components are all connected via directed lines. Strongly connected diagrams are a subclass of connected diagrams in which – for operators with dual-space components such as T̂ ^†^ ĤT̂ in eq – the removal of one T̂ ^†^ or T̂ still results in a connected diagram. On the other hand, if such a removal results in a disconnected diagram the term “weakly connected” is used. ?,? Restricting Hermitian linearized CC equations to the subclass of strongly connected diagrams ensures size-extensivity.

The LinCCD doubles amplitude equations are

$$[eqn]$$

where v _ pq _ ^ rs ^ are antisymmetrized 2-electron integrals ⟨rs∥pq⟩ and $[eqn]$ are index permutation operators. The first 3 terms (first line) of eq are often referred to as the driver terms (the latter two of which are responsible for the energy denominator in perturbation theory). The fourth and fifth terms (middle line of eq) are associated with ladder diagrams, and the final term (last line) emerges from ring and crossed-ring diagrams. The LinCCD diagrams are explicitly shown alongside their corresponding mathematical incarnation in Figure. For more detail pertaining to the CCD (and LinCCD) diagrams (including the basis set convergence of each term) the reader is referred to ref ?. Finally, the LinCCD energy expression in the spatial-orbital basis is explicitly:

$$[eqn]$$

Diagrammatic representation of the linearized CCD equations. The “driver” terms are displayed in black, “ladder” terms in blue, and “ring and crossed-ring” in red. When permutation operations are shown, the associated diagram corresponds to the contraction displayed in parentheses.

Linearized Ladder Coupled Cluster Doubles

2.2

Interestingly, to my knowledge it has yet to be observed that retaining only driver terms and ladder-type diagrams within LinCCD leads to naturally regular equations that are strongly resistant to divergence in molecular systems. Notably, linearized ladder approximations in various combinations with other RPA terms have been explored in the context of the homogeneous electron gas, ?−? ? ? ? but – to the best of my knowledge – were essentially abandoned before being applied in the context of chemistry. Given the tremendous volume of interest in pp-RPA/ladder-CCD approaches in quantum chemistry ?−? ? ? ? ? ? ? it appears pertinent to explore linearized ladder approximations.

Applying the same precedent of diagrammatic simplification as pp-RPA/ladder-CCD, I restrict LinCCD to ladder diagrams yielding

$$[eqn]$$

These LinLCCD equations are naturally regular. In principle, the inclusion of only ladder diagrams incorporates the most important contributions for describing strong correlation. However, to understand precisely how this is the case, it is helpful to notice that the contractions in terms 2 and 5 and terms 3 and 4 can be grouped together

$$[eqn]$$

By choosing a clever basis for eq, such as the one that diagonalizes the n _ v _ ^2^ × n _ v _ ^2^ matrix in term 1 and the n _ o _ ^2^ × n _ o _ ^2^ matrix in term 2, the amplitude equation can be reduced to a highly revealing linear form. Specifically, one can define a particle–particle (pp)- and hole–hole (hh)-blocked super Fock matrix with elements

$$[eqn]$$ $$[eqn]$$

where the eigenvalues of eq, the pp-Fock matrix (F pp), and eq, the hh-Fock matrix (F hh), are the hh- and pp-pair energies as estimated by sums of canonical one-particle orbital energies. Taking V pp and V hh as matrix representations of the pp- and hh-integrals from eq, one can solve the eigenvalue problems

$$[eqn]$$ $$[eqn]$$

where the eigenvalues are pp-pair (λ) and hh-pair (η) orbital energies in the dressed orbital basis, respectively. They can be decomposed into a one-particle pair contribution and a two-particle “dressing” supplied by the integrals

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ are the contributions from one-particle pair energies and the tilde designates quantities that are in the dressed orbital basis.

Rotation of the supermatrices defined above by eigenvectors X in the pp space and Y in the hh space leads to a diagonal representation of eq

$$[eqn]$$

The quantity in parentheses is diagonal and is thus invertible, leading to the amplitudes

$$[eqn]$$

Here, the elements of λ and η are expanded to emphasize that the denominator is constructed of inseparable dressed pair energies. This distinguishes LinLCCD as a coupled pair theory as opposed to the independent electron pair theory that constitutes second-order Møller–Plesset perturbation theory (MP2), to which these equations bear a striking resemblance. Finally, note that this basis transformation was accomplished via independent occupied-occupied and virtual–virtual orbital rotations, so the orbital-invariant LinLCCD energy does not change.

Such an isomorphism between MP2 and CCD equations has been noticed before, but only in the context of mosaic CCD.? In the linearized ladder context, the hh- and pp-pair energies have been shifted by the hh and pp integrals, respectively. This essentially results in a set of screened first-order amplitudes wherein the energy gap is widened by adding hh correlations to the one-particle occupied pair energies and pp correlations to one-particle virtual pair energies, making LinLCCD robust against divergence in small-gap systems. Unlike mosaic CCD and myriad other renormalized MP2 theories, ?−? ? ? ? the LinLCCD gap is widened in an amplitude-independent way, suggesting that the LinLCCD equations could be solved noniteratively, albeit such an approach would be quite impractical as it would require the diagonalization of a n _ v _ ^2^ × n _ v _ ^2^ matrix.

Incredibly, removing the ring and crossed-ring contractions from LinCCD also corrects the size-consistency errors in the parent method. I note that beyond size-consistency, LinLCCD is also size-extensive and orbital invariant. As the LinLCCD equations can be recast in the form of eq, the proof for size-consistency in this basis is trivial. Consider a system wherein molecules A and B are sufficiently far apart that their respective molecular orbitals may be trivially fragment-ascribed. After a rotation into the aforementioned dressed orbital basis eq becomes zero for all sums over disjoint orbitals by means of the electron repulsion integrals in the numerator, ensuring that the resultant energy satisfies E _ A∪B _ = E _ A _ + E _ B _.

My group is particularly interested in low-scaling approximations that describe static correlation qualitatively. While I will demonstrate the utility of LinLCCD, it does retain the most expensive $[eqn]$ particle–particle ladder term that is responsible for the $[eqn]$ cost of CCD. In the spirit of exploring low-scaling variants of linear ladder theories, I introduce one further approximation by completely removing the costly particle–particle ladder term to achieve

$$[eqn]$$

thus shifting only the occupied-pair energies by the hole–hole ladder term. Unlike the LinLCCD case, there is no diagrammatic justification for removing the particle–particle ladder term, but the resultant LinLCCD­(hh) approach is a potentially fruitful approximation that scales much more favorably as $[eqn]$ .

It should be noted that the LinLCCD­(hh) approximation does not sacrifice orbital invariance, size-consistency, or size-extensivity. Such claims can be verified by applying analogous occupied-occupied orbital rotations to eq such that term 2 becomes diagonal. Once again, this gives way to a set of dressed amplitude equations that can be written as

$$[eqn]$$

where now only the occupied pair orbital energies have been dressed with hh correlation. While these basis transformations are revealing as to the underlying nature of various linearized CC approximations, in practice my implementation solves eq self-consistently.

Results and Discussion

3

Bond Breaking

3.1

I first show the relative robustness of LinLCCD methods through a few simple bond dissociation potential energy surfaces. The simplest case of H_2_ in the STO-3G basis is shown in Figurea. As expected, LinCCD diverges rapidly as the H–H bond is stretched. However, both LinLCCD and the further pruned LinLCCD­(hh) approaches smoothly dissociate H_2_ toward some limit, which is exact in the case of LinLCCD­(hh). Of course, I note that these results also show that LinLCCD is no longer exact for all two electron systems, but neither is its (chemically very useful) parent ladder approximation, pp-RPA.

Dissociation curves for (a) hydrogen molecule in the minimal STO-3G basis set, (b) Hydrogen fluoride molecule in the aug-cc-pVQZ basis, , and (c) H6 in the cc-pVDZ basis. Dashed lines of like-color indicate the dissociation limit for a particular method. Dissociation limits were estimated at R = 106 Å, except in the case of H6 where R = 103 Å was used instead due to convergence difficulties.

Interestingly, if the linear direct ring and crossed-ring terms (i.e., without antisymmetrizing the two-electron repulsion integrals) are retained alongside the fully antisymmetric ladder terms, the resultant LinLdRxRCCD approach also does not diverge. LinLdRxRCCD differs from full LinCCD only by the exchange component of the linear ring and crossed-ring contractions, implying that the term most responsible for the instability of LinCCD is not the small orbital-energy denominator per se, but the exchange ring/crossed-ring terms. In the case of HF molecule dissociation in Figureb, this finding helps to explain why regularization by means of eliminating near-zero denominator terms in LinCCSD is only somewhat effective,? eventually breaking down at large R H–F, whereas LinLCCD, LinLCCD­(hh), and LinLdRxRCCD are all stable out to 10^6^ Å.

I note that others have put forth that the divergence of CC equations alongside other deficiencies in statically correlated systems manifest due to various exchange terms. ?,?,?,? Similar instabilities were also recently reported for renormalized propagator methods that include ring and crossed-ring exchange terms.? Of course, removing exchange terms can lead to undesirable self-interaction artifacts, as is well-known in the case of direct ring CCD (otherwise known as the direct ph-RPA), which overbinds significantly at equilibrium geometries. ?−? ? ? Though less common, removing exchange terms from ladder CCD has also been explored but appears to lead to less satisfactory results for bond dissociation energies.? Thus, if a self-interaction-free theory is desired that can smoothly dissociate bonds to a clear asymptotic limit, LinLCCD and LinLCCD­(hh) represent suitable options.

Lastly in the series of bond dissociation curves, I investigate hexagonal H_6_ dissociation in Figurec. The hexagonal H_6_ system is prototypical of strongly correlated systems in chemistry and is reminiscent of the Hubbard model Hamiltonian. As the H_6_ ring is expanded LinCCD rapidly diverges while all methods that exclude ring/crossed-ring exchange diagrams remain stable. The estimated asymptotic limit of LinLCCD­(hh) appears somewhat deceptive as plotted because it very slightly overestimates the correlation energy at large R, dipping below the full configuration interaction (FCI) reference curve by about 4 mE_h_. While the LinLCCD­(hh) potential curve is far too repulsive at intermediate R, the asymptotic limit remains very impressive for such an approximate scheme. Overall, all methods that remove the exchange ring/crossed-ring terms are robust for the bond stretching coordinates investigated in Figure while LinCCD fails in all cases, suggesting that ring/crossed-ring exchange contractions are to blame for the near-singular behavior of LinCCD.

Non-Bonded Interactions

3.2

Recall that LinLCCD is an approximation to pp-RPA, so it should not be expected to be a quantitatively accurate approach for total energies as pp-RPAand therefore LinLCCDdoes not offer ideal coverage of dynamical correlation effects. ?,?,? These deficiencies can be accounted for in pp-RPA by combination with density functionals, but I will not explore this here. ?,?,?−? ? ? The lack of dynamical correlation can be immediately seen in the potential energy curves of Figure, as LinLCCD and LinLCCD­(hh) both underestimate the magnitude of the correlation energy relative to FCI and LinCCD near equilibrium. However, much of the utility in RPA methods comes from how accurately they predict energy differences rather than the total energies themselves.? In this section I will explore the accuracy of LinLCCD for noncovalent interactions computed via

$$[eqn]$$

where E _ AB _ is the energy of the complex and E _ X _ is the monomer energy of fragment X. All interaction energy calculations have been counterpoise corrected.?

First, I examine the performance of LinLCCD, LinLCCD­(hh), LinCCD, CCD, and CCSD on the A24 data set of small dimers in Figurea.? These results are largely unremarkable as all methods perform statistically about the same with relatively low errors ranging between 0.1–0.2 kcal/mol. Some notable differences are seen in Figureb for the S22 data set which features several medium-sized π-stacked systems. ?,? Interestingly, while LinLCCD and LinLCCD­(hh) perform similarly with quite low mean absolute errors (MAE) of 0.4 kcal/mol, LinCCD performs poorly with a MAE of 1.4 kcal/mol. This poorer performance is in line with the full CCD and CCSD methods, which also under-perform on S22. While it is sensible that LinCCD, CCD, and CCSD perform similarly at equilibrium geometries, it is reasonable to wonder whether any particular interaction type (H-bonding, dispersion-bound, π-stacking, or mixed interactions) is responsible for the uniformly poor performance of these methods. However, Figure S1 suggests that the performance of all approaches remains consistent across interaction types. This is especially interesting in the case of LinLCCD methods, which are akin to renormalized MP2. Whereas MP2 tends to have a significant propensity to overbind π-stacked complexes by upward of 100%, ?,? no such bias is noted in LinLCCD approaches for the dimers in S22.

Error statistics for noncovalent interaction energies extrapolated to the complete basis set limit for (a) the A24 set of small dimers and (b) the S22 data set of small to medium sized dimers. The inset numbers indicate the mean absolute error in kcal/mol.

Both LinLCCD and LinLCCD­(hh) methods appear to perform at least as well as regularized perturbation theory approaches. ?,? These results suggest that, despite a lack of dynamical correlation, the energy differences obtained with LinLCCD approaches are quite accurate for nonbonded interactions, surpassing those of full CCD and CCSD.

In an effort to emphasize the relative affordability of the $[eqn]$ LinLCCD­(hh) approximation, I also present complete basis set limit extrapolated interaction energies for the L7 data set of large dimers in Figure.? I compare to the complete basis set limit domain-localized pair natural orbital (DLPNO) ?−? ? CCSD­(T_0_) benchmark data of Lao and co-workers.? I note that DLPNO–CCSD­(T_0_) results are quite sensitive to the particular thresholds chosen, so the benchmark data have some unknown error in addition to those imposed by the (T 0) correction. This error is roughly 2–6 kcal/mol for the π-stacked complexes in L7.?

Systems that comprise the L7 data set along with their commonly employed acronyms.

The results in Table show that MP2 dramatically overestimates the interaction energies in L7,? which features several large π-stacked systems. As LinLCCD­(hh) can be viewed as a renormalized MP2, it is of interest to contrast its performance with MP2 and with the renormalized MP2 approach known as size-consistent Brillouin–Wigner perturbation theory (BW-s2).? The results in Table suggest that including hole–hole relaxation in the one-particle energies can temper the overestimated interaction energies of conventional MP2, but not dramatically so. The MAE is reduced relative to MP2 by nearly 2 kcal/mol, which is an improvement but suggests that linear ladder correlation is insufficient to achieve quantitative accuracy. The renormalization supplied by BW-s2 is somewhat more effective at suppressing overcorrelation in the largest π-stacked systems than LinLCCD­(hh), but less aggressive in systems like GCGC, GGG, and CBH, leading to an overall MAE that is about 0.6 kcal/mol lower than LinLCCD­(hh). That said, the results between BW-s2 and LinLCCD­(hh) are comparable to within 12% of one another. While I have shown that LinLCCD­(hh) is affordable enough to be applied to such large systems and that the results are somewhat improved relative to MP2, some empiricism,? or other means of incorporating additional many-body screening effects could beget improvements.

1: Interaction Energies for L7 Data Set (kcal/mol)

Quest #8 Singlet/Triplet Gaps

3.3

Spin state energetics of transition metal complexes are a key quantity in the predictive modeling of energy relevant systems such as photoredox catalysts that are used to generate solar fuels such as H_2_ and to facilitate organic syntheses. ?−? ? ? While there remains a paucity of benchmark-quality spin-state energetics data for transition metal complexes in the literature, the small Quest #8 data set of Loos and co-workers does provide excellent data for comparisons with nonrelativistic quantum chemistry methods. ?,? The Quest #8 set contains 11 diatomic, monometallic transition metal molecules with nonrelativistic theoretical best estimate (TBE) values computed in the gas phase, making comparison with other nonrelativistic quantum chemistry methods more straightforward than (still useful) back-corrected experimental values.? Herein, I evaluate the performance of various correlated wave function theoretic approaches on the singlet/triplet gaps of the 7 molecule subset in Quest #8 that has a singlet ground state.

MAEs and maximum errors for the Quest #8 singlet/triplet gaps are shown in Figure. None of the methods that truncate at double substitutions are quantitatively accurate relative to the TBE benchmarks, but the inclusion of perturbative triple excitations in CCSD­(T) clearly has a large effect on the accuracy of the predicted gaps, reducing the MAE from 0.3 eV with CCSD to 0.07 eV with the inclusion of triples. As one might expect, LinCCD performs similarly well, albeit with a larger maximum error of 1.2 eV. Interestingly, while LinLCCD features a slightly larger MAE than LinCCD at 0.5 eV, it has a smaller maximum error than both LinCCD and CCSD at just 0.7 eV. Of all of the linearized CC approaches, LinLCCD­(hh) performs the worst with a MAE of 0.6 eV, but it still outperforms MP2. As the linearized ladder CCD approaches can be conceptualized as intermediate theories between MP2 and LinCCD, it makes sense that their MAEs fall in a hierarchical order MP2 > LinLCCD­(hh) > LinLCCD > LinCCD. While the results presented here are in line with expectations and are reasonably accurate, they once again point toward a need for incorporating more dynamical correlation into the LinLCCD approximation. Such studies are currently underway in my group.

Mean absolute errors (shown as colored bars) and maximum errors (shown as gray bars) for the lowest energy singlet/triplet gaps in Quest #8 as computed by various ΔCC and ΔMP2 methods. The bars are colored according to the category of approximation being applied, where Møller–Plesset perturbation theory is red, linearized CC methods are orange, and nonlinear CC approaches are green. The systems in question are shown in the inset.

Photolysis of Volatile Organic Compounds

3.4

The previous results focus mainly on energy differences between the electronic ground state and a high-spin triplet state of the system. This is one example of a ΔCC calculation, where the energy of each electronic configuration is optimized at the self-consistent field level and used as the reference for a subsequent non-Aufbau CC calculation. Excitation energies are then obtained by taking the difference between the CC energies of each respective calculation. In this section, I will explore the calculation of excited state bond dissociation curves via the ΔCC approach to assess LinLCCD for its recovery of potential surfaces of open-shell systems.

Volatile organic compounds are of intense interest in the atmospheric chemistry community and understanding their photochemistry can inform on human health ?,? and global climate modeling. ?−? ? One downstream product of CHCl_2_F, a recently phased-out refrigerant, is CF_3_COCl, which is known to decompose under UV irradiation.? The results in Figure show ΔLinLCCD calculations on several of the lowest-energy excited states in CF_3_COCl along the C–Cl bond stretching coordinate. Typically I would include single substitutions within the CC ansatz to model open-shell systems, but the presence of singles clusters threatens the stability of the excited-state configuration. By Thouless’ theorem,? single substitutions are equivalent to orbital rotations that could push the desired excited-state solutions toward the ground state. Therefore, a caveat to bear in mind in the following analysis is that the reference determinant for each ΔLinLCCD excited state is spin contaminated and without single substitutions to aid in spin purification the calculated excitation energies are likely underestimated. ?−? ? While approximate spin projection could be employed,? I do not expect this to impact the qualitative validity of the resultsespecially for the purposes of modeling the topography of each potential surface at large C–Cl distances.

Dissociation of CF3COCl along the C–Cl bond using ΔLinLCCD from reference configurations that represent the ground electronic state as well as n → π, p → π*, and p → σ* transitions.*

The dissociation curves for various states of CF_3_COCl are shown in Figure. All of the LinLCCD results are qualitatively consistent with the multireference calculations of ref ?. Namely, the S 1 state corresponds to a bound n → π* transition at 4.7 eV, which is in good agreement with the experimental band maximum of about 4.9 eV.? Furthermore, the S 2 and S 4 states are unbound and correspond to p → π* and p → σ* transitions, respectively. Both states lead to free dissociation to two different limits. This is expected, because dissociation along the ground state potential surface should lead to the homolytic cleavage of the C–Cl σ bond, populating the σ* orbital and resulting in the same dissociation limit as the S 4 state. In the case of the S 2 state, the occupied π* orbital corresponds to a qualitatively different configuration at dissociation.

The dissociation limits for the ground state and S 2 states are 3.5 and 6.5 eV, respectively. The former is very close to the extended multistate complete active space second-order perturbation theory (XMS-CASPT2)? dissociation energy for S 0 predicted in ref ?. of 3.44 eV. The energy difference between S 0 and S 2 states at dissociation is 3 eV and is also in excellent agreement with XMS-CASPT2 results. Despite spin contamination, the qualitative curvature of each surface is reasonable and LinLCCD with open-shell singlet references provides excellent bond dissociation energies in the case of CF_3_COCl.

Conclusions

4

In summary, I have introduced linearized ladder CCD equations alongside a linear hole–hole ladder approximation and examined their applicability to a range of chemistry contexts. LinLCCD and LinLCCD­(hh) are both robust when static correlation becomes important, and LinLCCD­(hh) serves as an affordable $[eqn]$ approximation that can be applied to large systems. I have also shown that the most problematic terms in LinCCD that lead to near-singular correlation energies are the ring and crossed-ring contractionsparticularly the exchange contributions therein. It is especially notable that LinLCCD is a size-consistent CID method that is obtained by removing linear terms rather than adding quadratic ones.? With future adaptations of LinLCCD and LinLCCD­(hh) to incorporate more dynamical correlation, these approaches could prove to be quite useful in the design of new computational methodologies for modeling strongly correlated systems reasonably well within single-reference approximations.

5. Computational Details

All calculations make use of the resolution-of-the-identity approximation and were performed using a developer version of Q-Chem v6.2.? The A24 and S22 complete basis set limit results were obtained using two-point aug-cc-pVDZ/aug-cc-pVTZ extrapolation using β = 2.51 and α = 4.3 as per Neese and Valeev.? The S22 calculations made use of frozen natural orbitals (FNOs), retaining 99.6% of the natural orbital occupation. ?−? ? The use of FNOs has been shown to be a robust approximation that provides benchmark-quality noncovalent interaction energies on the S22 set.? The L7 calculations were extrapolated to the complete basis set limit using the somewhat smaller Def2-ma-SVP and Def2-ma-TZVP basis sets for heavy atoms and the corresponding Def2-SVP/Def2-TZVP for H atoms. ?−? ? As Neese and Valeev do not provide parameters to minimize extrapolation errors for Karlsruhe basis sets with diffuse functions,? I computed the Hartree–Fock energies for L7 with the Def2-ma-QZVPP basis set and extrapolated the correlation energy with the more typical β = 3 parameter. The L7 calculations are the only ones that use the frozen core approximation. The Quest#8 calculations employ the aug-cc-pVTZ basis for even-handed comparison with the aug-cc-pVTZ reference data in ref ?. (due to limitations in Q-Chem, all I angular-momentum functions were removed from the auxiliary basis). Quest#8 systems of triplet multiplicity and the CF_3_COCl potential energy surfaces were calculated using unrestricted Hartree–Fock reference orbitals. Non-Aufbau reference configurations for CF_3_COCl were stabilized using a combination of state-targeted energy projection and initial maximum overlap method algorithms. ?−? ?

Supplementary Material