-
An "ultimate" coupled cluster method based entirely on $T_2$
Authors:
Zachary W. Windom,
Ajith Perera,
Rodney J. Bartlett
Abstract:
Electronic structure methods built around double-electron excitations have a rich history in quantum chemistry. However, it seems to be the case that such methods are only suitable in particular situations and are not naturally equipped to simultaneously handle the variety of electron correlations that might be present in chemical systems. To this end, the current work seeks a computationally effi…
▽ More
Electronic structure methods built around double-electron excitations have a rich history in quantum chemistry. However, it seems to be the case that such methods are only suitable in particular situations and are not naturally equipped to simultaneously handle the variety of electron correlations that might be present in chemical systems. To this end, the current work seeks a computationally efficient, low-rank, "ultimate" coupled cluster method based exclusively on $T_2$ and its products which can effectively emulate more "complete" methods that explicitly consider higher-rank, $T_{2m}$ operators. We introduce a hierarchy of methods designed to systematically account for higher, even order cluster operators - like $T_4, T_6, \cdots, T_{2m}$ - by invoking tenets of the factorization theorem of perturbation theory and expectation-value coupled cluster theory. It is shown that each member within this methodological hierarchy is defined such that both the wavefunction and energy are correct through some order in many-body perturbation theory (MBPT), and can be extended up to arbitrarily high orders in $T_2$. The efficacy of such approximations are determined by studying the potential energy surface of several prototypical systems that are chosen to represent both non-dynamic, static, and dynamic correlation regimes. We find that the proposed hierarchy of augmented $T_2$ methods essentially reduce to standard CCD for problems where dynamic electron correlations dominate, but offer improvements in situations where non-dynamic and static correlations become relevant. A notable highlight of this work is that the cheapest methods in this hierarchy - which are correct through fifth-order in MBPT - consistently emulate the behavior of the $\mathcal{O}(N^{10})$ CCDQ method, yet only require a $\mathcal{O}(N^{6})$ algorithm by virtue of factorized intermediates.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry
Authors:
James H. Thorpe,
Zachary W. Windom,
Rodney J. Bartlett,
Devin A. Matthews
Abstract:
Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as $\mathcal{O}(N^9)$ with the number of basis functions. It very difficult to predict, a pri…
▽ More
Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as $\mathcal{O}(N^9)$ with the number of basis functions. It very difficult to predict, a priori, the effect correlation past CCSD(T) has on a give reaction energies. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of many-body perturbation theory that may address these challenges. We show that the $\mathcal{O}(N^7)$, factorized CCSD(TQ${}_\text{f}$) method introduces minimal error to predicted correlation and reaction energies as compared to the $\mathcal{O}(N^9)$ CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)${}_Λ$ contributions to reaction energies, as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQ${}_\text{f}$)/cc-pVDZ is capable of predicting CCSDT(Q)${}_Λ$/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol${}^{-1}$ and a RMST of 0.52 kcal mol${}^{-1}$ when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably ballpark how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
An attractive way to correct for missing singles excitations in unitary coupled cluster doubles theory
Authors:
Zachary W. Windom,
Daniel Claudino,
Rodney J. Bartlett
Abstract:
Coupled cluster methods based exclusively on double excitations are comparatively "cheap" and interesting model chemistries, as they are typically able to capture the bulk of the dynamical electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron pair-restricted…
▽ More
Coupled cluster methods based exclusively on double excitations are comparatively "cheap" and interesting model chemistries, as they are typically able to capture the bulk of the dynamical electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron pair-restricted $T_2$ operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent in which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Our analysis includes the derivations of finite-order, UCC energy functionals which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing singles excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
An assessment of frozen natural orbitals and band gaps using equation of motion coupled cluster theory: a case study on polyacene and trans-polyacetylene
Authors:
Zachary W. Windom,
AV Lam,
Ajith Perera,
Rodney J. Bartlett
Abstract:
Frozen natural orbitals (FNOs) are used to augment IP/EA-EOM-CCSD calculations targeting the band gap of trans-polyacetylene and polyacene. We show the resulting electron affinities (EAs), ionization potentials (IPs), and extrapolated band gaps incur errors that are largely tunable to a desired accuracy, yet require many orders of magnitude fewer core-hours as compared to the corresponding full ca…
▽ More
Frozen natural orbitals (FNOs) are used to augment IP/EA-EOM-CCSD calculations targeting the band gap of trans-polyacetylene and polyacene. We show the resulting electron affinities (EAs), ionization potentials (IPs), and extrapolated band gaps incur errors that are largely tunable to a desired accuracy, yet require many orders of magnitude fewer core-hours as compared to the corresponding full calculation. The relationship between various FNO truncation schemes and (cc-pV$n$Z) basis set is also examined.
△ Less
Submitted 13 February, 2024;
originally announced February 2024.
-
A new "gold standard": perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing
Authors:
Zachary W. Windom,
Daniel Claudino,
Rodney J. Bartlett
Abstract:
A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework in dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitatio…
▽ More
A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework in dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the "gold standard" CCSD(T) method, which bolsters the baseline CCSD with effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes, which can be obtained on a quantum computer, are leveraged by a classical computer to evaluate energy corrections associated with triple excitations - leading to the UCCSD[T] and UCCSD(T*) methods. The rationale behind these choices is shown to be rigorous by studying the properties of finite-order UCC energy functionals. Although our efforts do not support the addition of the fifth-order contribution as in the (T) correction, comparisons are nevertheless made using a hybrid UCCSD(T) approach. We assess the performance of these approaches on a collection of small molecules, and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.
△ Less
Submitted 11 January, 2024;
originally announced January 2024.
-
Nonconvergence of the Feynman-Dyson diagrammatic perturbation expansion of propagators
Authors:
So Hirata,
Ireneusz Grabowski,
J. V. Ortiz,
Rodney J. Bartlett
Abstract:
Using a general-order ab initio many-body Green's function method, we numerically illustrate several pathological behaviors of the Feynman-Dyson diagrammatic perturbation expansion of one-particle many-body Green's functions as electron Feynman propagators. (i) The perturbation expansion of the frequency-dependent self-energy is not convergent at the exact self-energy in many frequency domains. (i…
▽ More
Using a general-order ab initio many-body Green's function method, we numerically illustrate several pathological behaviors of the Feynman-Dyson diagrammatic perturbation expansion of one-particle many-body Green's functions as electron Feynman propagators. (i) The perturbation expansion of the frequency-dependent self-energy is not convergent at the exact self-energy in many frequency domains. (ii) An odd-perturbation-order self-energy has a qualitatively wrong shape and, as a result, many roots of the corresponding Dyson equation are nonphysical in that the poles may be complex or residues can exceed unity or be negative. (iii) A higher even-order self-energy consists of vertical lines at many frequencies, predicting numerous phantom poles with zero residues. (iv) Infinite partial resummations of diagrams by vertex or edge renormalization tend to exacerbate these pathologies. (v) The nonconvergence is caused by the nonanalyticity of the rational-function form of the exact Green's function at many frequencies, where the radius of convergence of its Taylor expansion is zero. This is consistent with the fact that (vi) Padé approximants (power-series expansions of a rational function) can largely restore the correct shape and poles of the Green's function. Nevertheless, not only does the nonconvergence render higher-order Feynman-Dyson diagrammatic perturbation theory useless for many lower-lying ionization or higher-lying electron-attachment states, but it also calls into question the validity of its combined use with the ansätze requiring the knowledge of all poles and residues. Such ansätze include the Galitskii-Migdal identity, the self-consistent Green's function methods, and some models of the algebraic diagrammatic construction.
△ Less
Submitted 25 April, 2024; v1 submitted 5 December, 2023;
originally announced December 2023.
-
A route to improving RPA excitation energies through its connection to equation-of-motion coupled cluster theory
Authors:
Varun Rishi,
Ajith Perera,
Rodney J. Bartlett
Abstract:
We revisit the connection between equation-of-motion coupled cluster (EOM-CC) and random phase approximation (RPA) explored recently by Berkelbach [J. Chem. Phys. 149, 041103 (2018)] and unify various methodological aspects of these diverse treatment of ground and excited states. The identity of RPA and EOM-CC based on the ring coupled cluster doubles is established with numerical results which wa…
▽ More
We revisit the connection between equation-of-motion coupled cluster (EOM-CC) and random phase approximation (RPA) explored recently by Berkelbach [J. Chem. Phys. 149, 041103 (2018)] and unify various methodological aspects of these diverse treatment of ground and excited states. The identity of RPA and EOM-CC based on the ring coupled cluster doubles is established with numerical results which was proved previously on theoretical grounds. We then introduce new approximations in EOM-CC and RPA family of methods, assess their numerical performance and explore a way to reap the benefits of such a connection to improve on excitation energies. Our results suggest that addition of perturbative corrections to account for double excitations and missing exchange effects could result in significantly improved estimates.
△ Less
Submitted 6 November, 2020; v1 submitted 1 August, 2020;
originally announced August 2020.
-
Computing solvated excited states using fragment-effective-field coupled-cluster perturbation theory with application to the electronic spectra of nucleobases in water
Authors:
Jason N. Byrd,
Robert W. Molt, Jr.,
Rodney J. Bartlett,
Beverly A. Sanders
Abstract:
Molecular fragment or embedding methods are powerful techniques for overcoming scalability limitations in electronic structure theory by dividing large molecular systems into individual units that are small enough to be treated using standard techniques. Our group has developed a linear-scaling coupled-cluster based fragment effective field method, molecular cluster perturbation theory (MCPT), tha…
▽ More
Molecular fragment or embedding methods are powerful techniques for overcoming scalability limitations in electronic structure theory by dividing large molecular systems into individual units that are small enough to be treated using standard techniques. Our group has developed a linear-scaling coupled-cluster based fragment effective field method, molecular cluster perturbation theory (MCPT), that provides a systematically improvable description of the electronic wavefunction and is scalable to thousands of atoms using the Aces4 massively parallel software package. We will present a systematic computational benchmarking study of small organic molecule solvation free energies obtained from MCPT using various perturbation corrections using experimental and large-scale ab initio reference values. We will also present recent work on the relative crystal lattice polymorph energies of several molecular crystals where the roles of many-body dispersion and intermol ecular electron exchange will be investigated. The molecular solvent and crystal systems presented here are carefully chosen so as to provide a robust test set for the fragment electronic structure community.
△ Less
Submitted 10 August, 2017;
originally announced August 2017.
-
Predictive coupled-cluster isomer orderings for some Si${}_n$C${}_m$ ($m, n\le 12$) clusters; A pragmatic comparison between DFT and complete basis limit coupled-cluster benchmarks
Authors:
Jason N. Byrd,
Jesse J. Lutz,
Yifan Jin,
Duminda S. Ranasinghe,
John A. Montgomery Jr.,
Ajith Perera,
Xiaofeng F. Duan,
Larry W. Burggraf,
Beverly A. Sanders,
Rodney J. Bartlett
Abstract:
The accurate determination of the preferred ${\rm Si}_{12}{\rm C}_{12}$ isomer is important to guide experimental efforts directed towards synthesizing SiC nano-wires and related polymer structures which are anticipated to be highly efficient exciton materials for opto-electronic devices. In order to definitively identify preferred isomeric structures for silicon carbon nano-clusters, highly accur…
▽ More
The accurate determination of the preferred ${\rm Si}_{12}{\rm C}_{12}$ isomer is important to guide experimental efforts directed towards synthesizing SiC nano-wires and related polymer structures which are anticipated to be highly efficient exciton materials for opto-electronic devices. In order to definitively identify preferred isomeric structures for silicon carbon nano-clusters, highly accurate geometries, energies and harmonic zero point energies have been computed using coupled-cluster theory with systematic extrapolation to the complete basis limit for set of silicon carbon clusters ranging in size from SiC$_3$ to ${\rm Si}_{12}{\rm C}_{12}$. It is found that post-MBPT(2) correlation energy plays a significant role in obtaining converged relative isomer energies, suggesting that predictions using low rung density functional methods will not have adequate accuracy. Utilizing the best composite coupled-cluster energy that is still computationally feasible, entailing a 3-4 SCF and CCSD extrapolation with triple-$ζ$ (T) correlation, the {\it closo} ${\rm Si}_{12}{\rm C}_{12}$ isomer is identified to be the preferred isomer in support of previous calculations [J. Chem. Phys. 2015, 142, 034303]. Additionally we have investigated more pragmatic approaches to obtaining accurate silicon carbide isomer energies, including the use of frozen natural orbital coupled-cluster theory and several rungs of standard and double-hybrid density functional theory. Frozen natural orbitals as a way to compute post MBPT(2) correlation energy is found to be an excellent balance between efficiency and accuracy.
△ Less
Submitted 21 June, 2016;
originally announced June 2016.
-
Approximating electronically excited states with equation-of-motion linear coupled-cluster theory
Authors:
Jason N. Byrd,
Varun Rishi,
Ajith Perera,
Rodney J. Bartlett
Abstract:
A new perturbative approach to canonical equation-of-motion coupled-cluster theory is presented using coupled-cluster perturbation theory. A second-order Møller-Plesset partitioning of the Hamiltonian is used to obtain the well known equation-of-motion many-body perturbation theory (EOM-MBPT(2)) equations and two new equation-of-motion methods based on the linear coupled-cluster doubles (EOM-LCCD)…
▽ More
A new perturbative approach to canonical equation-of-motion coupled-cluster theory is presented using coupled-cluster perturbation theory. A second-order Møller-Plesset partitioning of the Hamiltonian is used to obtain the well known equation-of-motion many-body perturbation theory (EOM-MBPT(2)) equations and two new equation-of-motion methods based on the linear coupled-cluster doubles (EOM-LCCD) and linear coupled-cluster singles and doubles (EOM-LCCSD) wavefunctions. This is achieved by performing a short-circuiting procedure on the MBPT(2) similarity transformed Hamiltonian. These new methods are benchmarked against very accurate theoretical and experimental spectra from 25 small organic molecules. It is found that the proposed methods have excellent agreement with canonical EOM-CCSD state for state orderings and relative excited state energies as well as acceptable quantitative agreement for absolute excitation energies compared with the best estimate theory and experimental spectra.
△ Less
Submitted 7 July, 2015;
originally announced July 2015.
-
Molecular Cluster Perturbation Theory. I. Formalism
Authors:
Jason N. Byrd,
Nakul Jindal,
Robert W. Molt, Jr.,
Rodney J. Bartlett,
Beverly A. Sanders,
Victor F. Lotrich
Abstract:
We present second-order molecular cluster perturbation theory (MCPT(2)), a linear scaling methodology to calculate arbitrarily large systems with explicit calculation of individual wavefunctions in a coupled-cluster framework. This new MCPT(2) framework uses coupled-cluster perturbation theory and an expansion in terms of molecular dimer interactions to obtain molecular wavefunctions that are infi…
▽ More
We present second-order molecular cluster perturbation theory (MCPT(2)), a linear scaling methodology to calculate arbitrarily large systems with explicit calculation of individual wavefunctions in a coupled-cluster framework. This new MCPT(2) framework uses coupled-cluster perturbation theory and an expansion in terms of molecular dimer interactions to obtain molecular wavefunctions that are infinite-order in both the electronic fluctuation operator and all possible dimer (and products of dimers) interactions. The MCPT(2) framework has been implemented in the new SIA/Aces4 parallel architecture, making use of the advanced dynamic memory control and fine grained parallelism to perform very large explicit molecular cluster calculations. To illustrate the power of this method, we have computed energy shifts, lattice site dipole moments, and harmonic vibrational frequencies via explicit calculation of the bulk system for the polar and non-polar polymorphs of solid hydrogen fluoride. The explicit lattice size (without using any periodic boundary conditions) was expanded up to 1,000 HF molecules, with 32,000 basis functions and 10,000 electrons. Our obtained HF lattice site dipole moments and harmonic vibrational frequencies agree well with the existing literature.
△ Less
Submitted 23 March, 2015; v1 submitted 6 September, 2014;
originally announced September 2014.
-
Correlation correction to configuration interaction singles from coupled cluster perturbation theory
Authors:
Jason N. Byrd,
Victor F. Lotrich,
Rodney J. Bartlett
Abstract:
A new state specific correlation correction to configuration interaction singles (CIS) excitation energies is preseted using coupled cluster perturbation theory (CCPT). General expressions for CIS-CCPT are derived and expanded explicitly to first order in the wavefunction and second order in the energy. By virtue of the nature of CCPT this method is {\it a priori} size extensive and incorporates i…
▽ More
A new state specific correlation correction to configuration interaction singles (CIS) excitation energies is preseted using coupled cluster perturbation theory (CCPT). General expressions for CIS-CCPT are derived and expanded explicitly to first order in the wavefunction and second order in the energy. By virtue of the nature of CCPT this method is {\it a priori} size extensive and incorporates infinite order effects into the wavefunction. This results in a balanced singles space excited state theory that at second order is an improvement over the ubiquitous CIS(D) method and comparable in quality to equation of motion coupled cluster (EOM-CC). A modest test set composed of the first four excited states from nine small organic molecules was used to quantify the accuracy and consistency of the CIS-CCPT2 excitation energies and density of states. We find that CIS-CCPT2 has a standard deviation error of 0.18 eV for excitation energies and 0.14 eV for density of states compared to EOM-CC, a factor of two better than CIS(D) with a significant reduction in the maximum deviation as well
△ Less
Submitted 23 April, 2014;
originally announced April 2014.