Robust Tensor Hypercontraction of the Particle-Particle Ladder Term in Equation-of-Motion Coupled Cluster Theory.

One method of representing a high-rank tensor as a (hyper-)product of lower-rank tensors is the tensor hypercontraction (THC) method of Hohenstein et al. This strategy has been found to be useful for reducing the polynomial scaling of coupled-cluster methods by representation of a four-dimensional tensor of electron-repulsion integrals in terms of five two-dimensional matrices. Pierce et al. have already shown that the application of a robust form of THC to the particle-particle ladder (PPL) term reduces the cost of this term in couple-cluster singles and doubles (CCSD) from O(N6) to O(N5) with negligible errors in energy with respect to the density-fitted variant. In this work, we have implemented the least-squares variant of THC (LS-THC) which does not require a nonlinear tensor factorization, including the robust form (R-LS-THC), for the calculation of the excitation and electron attachment energies using equation-of-motion coupled cluster methods EOMEE-CCSD and EOMEA-CCSD, respectively. We have benchmarked the effect of the R-LS-THC-PPL approximation on excitation energies using the comprehensive QUEST database and the accuracy of electron attachment energies using the NAB22 database. We find that errors on the order of 1 meV are achievable with a reduction in total calculation time of approximately 5 ×.

Transition Moments for STEOM-CCSD with Core Triples.

Similarity-transformed equation-of-motion coupled-cluster theory (STEOM-CC) is an alternative approach to equation-of-motion coupled-cluster theory for excited states (EOMEE-CC), which uses a second similarity transformation of the Hamiltonian, followed by diagonalization in a small (CI singles-like) excitation space, even when single and double excitations are included in the transformation. In addition to vertical excitation energies, transition moments measure the strength of the interactions between states determining absorption, emission, and other processes. In STEOM-CCSD, transition moments are calculated in a straightforward manner as biorthogonal expectation values using both the left- and right-hand solutions, with the main difference from EOMEE-CC being the inclusion of the transformation operator. We recently developed an extension of STEOM-CCSD to core excitations, CVS-STEOM-CCSD+cT, which includes triple excitations and the well-known core-valence separation for the core ionization potential calculations. In this work, we derived transition moments for core-excited states with core triple excitations, including both ground-to-core-excited and valence-to-core-excited transitions. The improvement of the computed transition moments of the CVS-STEOM-CCSD+cT method is compared to standard CVS-STEOMEE-CCSD and CVS-EOMEE-CCSD for our previously published small-molecule benchmark set.

Open-Shell Tensor Hypercontraction.

The extension of least-squares tensor hypercontracted second- and third-order Møller-Plesset perturbation theory (LS-THC-MP2 and LS-THC-MP3) to open-shell systems is an important development due to the scaling reduction afforded by THC and the ubiquity of molecular ions, radicals, and other open-shell reactive species. The complexity of wavefunction-based quantum-chemical methods such as Møller-Plesset and coupled cluster theory is reflected in the steep scaling of the computational costs with the molecular size. The least-squares tensor hypercontraction (LS-THC) method is an efficient, single-step factorization for the two-electron integral tensor but can also be used to factorize the double excitation amplitudes, leading to significant scaling reduction. Here we extend this promising method to open-shell variants of LS-THC-MP2 and -MP3 by using diagrammatic techniques and explicit spin summation. The accuracy of the resulting methods for open-shell species is benchmarked on standard test systems such as regular alkanes as well as realistic systems involving bond breaking, radical stabilization, and other effects. We find that open-shell LS-THC-MPn methods exhibit errors highly comparable to those produced by closed-shell LS-THC-MPn and are highly insensitive to particular chemical interactions, geometries, or even moderate spin contamination.

Ab Initio Investigation of Intramolecular Charge Transfer States in DMABN by Calculation of Excited State X-ray Absorption Spectra.

Dual fluorescence in 4-(dimethylamino)benzonitrile (DMABN) and its derivatives in polar solvents has been studied extensively for the past several decades. An intramolecular charge transfer (ICT) minimum on the excited state potential energy surface, in addition to the localized low-energy (LE) minimum, has been proposed as a mechanism for this dual fluorescence, with large geometric relaxation and molecular orbital reorganization a key feature of the ICT pathway. Herein, we have used both equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) and time-dependent density functional (TDDFT) methods to investigate the landscape of excited state potential energy surfaces across a number of geometric conformations proposed as ICT structures. In order to correlate these geometries and valence excited states in terms of potential experimental observables, we have calculated the nitrogen K-edge ground and excited state absorption spectra for each of the predicted "signpost" structures and identified several key spectral features that could be used to interpret a future time-resolved X-ray absorption experiment.

Theoretical Investigation of the X-ray Stark Effect in Small Molecules.

We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.

Allosteric regulation in STAT3 interdomains is mediated by a rigid core: SH2 domain regulation by CCD in D170A variant.

Signal Transducer and Activator of Transcription 3 (STAT3) plays a crucial role in cancer development and thus is a viable target for cancer treatment. STAT3 functions as a dimer mediated by phosphorylation of the SRC-homology 2 (SH2) domain, a key target for therapeutic drugs. While great efforts have been employed towards the development of compounds that directly target the SH2 domain, no compound has yet been approved by the FDA due to a lack of specificity and pharmacologic efficacy. Studies have shown that allosteric regulation of SH2 via the coiled-coil domain (CCD) is an alternative drug design strategy. Several CCD effectors have been shown to modulate SH2 binding and affinity, and at the time of writing at least one drug candidate has entered phase I clinical trials. However, the mechanism for SH2 regulation via CCD is poorly understood. Here, we investigate structural and dynamic features of STAT3 and compare the wild type to the reduced function variant D170A in order to delineate mechanistic differences and propose allosteric pathways. Molecular dynamics simulations were employed to explore conformational space of STAT3 and the variant, followed by structural, conformation, and dynamic analysis. The trajectories explored show distinctive conformational changes in the SH2 domain for the D170A variant, indicating long range allosteric effects. Multiple analyses provide evidence for long range communication pathways between the two STAT3 domains, which seem to be mediated by a rigid core which connects the CCD and SH2 domains via the linker domain (LD) and transmits conformational changes through a network of short-range interactions. The proposed allosteric mechanism provides new insight into the understanding of intramolecular signaling in STAT3 and potential pharmaceutical control of STAT3 specificity and activity.

A Mountaineering Strategy to Excited States: Revising Reference Values with EOM-CC4.

In the framework of the computational determination of highly accurate vertical excitation energies in small organic compounds, we explore the possibilities offered by the equation-of-motion formalism relying on the approximate fourth-order coupled-cluster (CC) method, CC4. We demonstrate, using an extended set of more than 200 reference values based on CC including up to quadruples excitations (CCSDTQ), that CC4 is an excellent approximation to CCSDTQ for excited states with a dominant contribution from single excitations with an average deviation as small as 0.003 eV. We next assess the accuracy of several additive basis set correction schemes, in which vertical excitation energies obtained with a compact basis set and a high-order CC method are corrected with lower-order CC calculations performed in a larger basis set. Such strategies are found to be overall very beneficial, though their accuracy depends significantly on the actual scheme. Finally, CC4 is employed to improve several theoretical best estimates of the QUEST database for molecules containing between four and six (nonhydrogen) atoms, for which previous estimates were computed at the CCSDT level.

Accurate Core-Excited States via Inclusion of Core Triple Excitations in Similarity-Transformed Equation-of-Motion Theory.

The phenomenon of orbital relaxation upon excitation of core electrons is a major problem in the linear-response treatment of core-hole spectroscopies. Rather than addressing relaxation through direct dynamical correlation of the excited state via equation-of-motion coupled cluster theory (EOMEE-CC), we extend the alternative similarity-transformed equation-of-motion coupled cluster theory (STEOMEE-CC) by including the core-valence separation (CVS) and correlation of triple excitations only within the calculation of core ionization energies. This new method, CVS-STEOMEE-CCSD+cT, significantly improves on CVS-EOMEE-CCSD and unmodified CVS-STEOMEE-CCSD when compared to full CVS-EOM-CCSDT for K-edge core-excitation energies of a set of small molecules. The improvement in both absolute and relative (shifted) peak positions is nearly as good as that for transition-potential coupled cluster (TP-CC), which includes an explicit treatment of orbital relaxation, and CVS-EOMEE-CCSD*, which includes a perturbative treatment of triple excitations.

Quadratic Unitary Coupled-Cluster Singles and Doubles Scheme: Efficient Implementation, Benchmark Study, and Formulation of an Extended Version.

An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin-orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S0 → S1 vertical excitation energies for cis- and trans-4a,4b-dihydrotriphenylene. The qUCCSD scheme has been shown to provide improved excitation energies compared with the UCC3 scheme formulated based on perturbation theory. A UCC truncation scheme that can provide excitation energies correct through the fourth order is also presented to further improve the accuracy of the qUCCSD scheme.

Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs.

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

How accurate are EOM-CC4 vertical excitation energies?

We report the first investigation of the performance of EOM-CC4-an approximate equation-of-motion coupled-cluster model, which includes iterative quadruple excitations-for vertical excitation energies in molecular systems. By considering a set of 28 excited states in 10 small molecules for which we have computed CC with singles, doubles, triples, quadruples, and pentuples and full configuration interaction reference energies, we show that, in the case of excited states with a dominant contribution from the single excitations, CC4 yields excitation energies with sub-kJ mol-1 accuracy (i.e., error below 0.01 eV), in very close agreement with its more expensive CC with singles, doubles, triples, and quadruples parent. Therefore, if one aims at high accuracy, CC4 stands as a highly competitive approximate method to model molecular excited states, with a significant improvement over both CC3 and CC with singles, doubles, and triples. Our results also evidence that, although the same qualitative conclusions hold, one cannot reach the same level of accuracy for transitions with a dominant contribution from the double excitations.

A critical analysis of least-squares tensor hypercontraction applied to MP3.

The least-squares tensor hypercontraction (LS-THC) approach is a promising method of reducing the high polynomial scaling of wavefunction methods, for example, those based on many-body perturbation theory or coupled cluster. Here, we focus on LS-THC-MP3 and identify four variants with differing errors and efficiency characteristics. The performance of LS-THC-MP3 is analyzed for regular test systems with up to 40 first-row atoms. We also analyze the size-extensivity/size-consistency and grid- and basis set dependence of LS-THC-MP3. Overall, the errors observed are favorably small in comparison with standard density fitting, and a more streamlined method of generating grids via pruning is suggested. A practical crossover (the point at which LS-THC-MP3 is cheaper than the canonical method) is achieved around 240 correlated electrons. Despite several drawbacks of LS-THC that have been identified: an initial non-linearity of error when increasing system size, poor description of angular correlation, and a potentially large increase in error with the basis set size, the results show that LS-THC has significant potential for practical application to MP3 and other wavefunction methods.

Transition-potential coupled cluster.

The problem of orbital relaxation in computational core-hole spectroscopies, including x-ray absorption and x-ray photoionization, has long plagued linear response approaches, including equation-of-motion coupled cluster with singles and doubles (EOM-CCSD). Instead of addressing this problem by including additional electron correlation, we propose an explicit treatment of orbital relaxation via the use of "transition potential" reference orbitals, leading to a transition-potential coupled cluster (TP-CC) family of methods. One member of this family, in particular, TP-CCSD(12), is found to essentially eliminate the orbital relaxation error and achieve the same level of accuracy for the core-hole spectra as is typically expected of EOM-CCSD in the valence region. These results show that very accurate x-ray absorption spectra for molecules with first-row atoms can be computed at a cost essentially the same as that for EOM-CCSD.

Analytic Gradients of Approximate Coupled Cluster Methods with Quadruple Excitations.

The analytic gradient theory for both iterative and noniterative coupled-cluster approximations that include connected quadruple excitations is presented. These methods include, in particular, CCSDT(Q), which is an analog of the well-known CCSD(T) method which starts from the full CCSDT method rather than CCSD. The resulting methods are implemented in the CFOUR program suite, and pilot applications are presented for the equilibrium geometries and harmonic vibrational frequencies of the simplest Criegee intermediate, CH2OO, as well as to the isomerization pathway between dimethylcarbene and propene. While all methods are seen to approximate the full CCSDTQ results well for "well-behaved" systems, the more difficult case of the Criegee intermediate shows that CCSDT(Q), as well as certain iterative approximations, display problematic behavior.

Coupled-cluster techniques for computational chemistry: The CFOUR program package.

An up-to-date overview of the CFOUR program system is given. After providing a brief outline of the evolution of the program since its inception in 1989, a comprehensive presentation is given of its well-known capabilities for high-level coupled-cluster theory and its application to molecular properties. Subsequent to this generally well-known background information, much of the remaining content focuses on lesser-known capabilities of CFOUR, most of which have become available to the public only recently or will become available in the near future. Each of these new features is illustrated by a representative example, with additional discussion targeted to educating users as to classes of applications that are now enabled by these capabilities. Finally, some speculation about future directions is given, and the mode of distribution and support for CFOUR are outlined.

Improved Grid Optimization and Fitting in Least Squares Tensor Hypercontraction.

A new method for generating fitting grids for least-squares tensor hypercontraction (LS-THC) is presented. This method draws inspiration from the related interpolative separable density fitting (ISDF) technique but uses only a pivoted Cholesky decomposition of the metric matrix, S, already computed as a matter of course in LS-THC. The size and quality of the resulting grid are controlled by a user-defined cutoff parameter and the size of the starting grid. Additionally, the Cholesky-based method provides an alternative and possible more numerically stable method for performing the least-squares fit. The quality of the grids produced is evaluated for LS-DF-THC-MP2 calculations on retinal and benzene, the former with a large starting grid and small cc-pVDZ basis set, and the latter with a wide range of grids and basis sets. The error and grid size are found to be well-controlled by either the cutoff parameter (with a large starting grid) or the starting grid size (with a tight cutoff) and highly predictable. The Cholesky-based method is also able to generate unique grids tailored to different charge distributions, for example the (ab|, (ai|, and (ij| distributions that arise in the molecular orbital integrals. While only the (ai| grid directly affects the MP2 energy, the relative sizes of the other grids are examined.

Accuracy of Coupled Cluster Excited State Potential Energy Surfaces.

The validation of the quality of the description of excited electronic states is of special importance in quantum chemistry as the general reliability of ab initio methods shows a much larger variation for these states than for the ground state. In this study, we investigate the quality of excited state energy gradients and potential energy surfaces on selected systems, as provided by the single reference coupled cluster variants CC2, CCSD, CCSD(T)(a)*, and CC3. Gradients and surface plots that follow the Franck-Condon forces are compared to the respective CCSDT reference values, thereby establishing a useful strategy for judging each variant's accuracy. The results reveal serious flaws of lower order methods - in particular, CC2 - in several situations where they otherwise give accurate vertical excitation energies, as well as excellent accuracy and consistency of the recently proposed CCSD(T)(a)* method.

Fourth-order vibrational perturbation theory with the Watson Hamiltonian: Report of working equations and preliminary results.

A derivation of fourth-order vibrational perturbation theory (VPT4) based on the Watson Hamiltonian in dimensionless rectilinear normal coordinates is presented. Terms that are linear and cubic in the (nk + 12 ), with nk being the zeroth-order harmonic oscillator quantum numbers, appear in fourth order and extend the much simpler second-order vibrational perturbation theory model. The rather involved expressions for the fourth-order terms are derived with Rayleigh-Schrödinger perturbation theory, the process of verifying their correctness is described, and a computer code to generate the VPT4 constants from the potential energy surface derivatives is provided. The paper concludes with numerical examples featuring the H2O, Si2C, and cyclic-C3H2 molecules.

Geometric Energy Derivatives at the Complete Basis Set Limit: Application to the Equilibrium Structure and Molecular Force Field of Formaldehyde.

Geometric energy derivatives which rely on core-corrected focal-point energies extrapolated to the complete basis set (CBS) limit of coupled cluster theory with iterative and noniterative quadruple excitations, CCSDTQ and CCSDT(Q), are used as elements of molecular gradients and, in the case of CCSDT(Q), expansion coefficients of an anharmonic force field. These gradients are used to determine the CCSDTQ/CBS and CCSDT(Q)/CBS equilibrium structure of the S0 ground state of H2CO where excellent agreement is observed with previous work and experimentally derived results. A fourth-order expansion about this CCSDT(Q)/CBS reference geometry using the same level of theory produces an exceptional level of agreement to spectroscopically observed vibrational band origins with a MAE of 0.57 cm-1. Second-order vibrational perturbation theory (VPT2) and variational discrete variable representation (DVR) results are contrasted and discussed. Vibration-rotation, anharmonicity, and centrifugal distortion constants from the VPT2 analysis are reported and compared to previous work. Additionally, an initial application of a sum-over-states fourth-order vibrational perturbation theory (VPT4) formalism is employed herein, utilizing quintic and sextic derivatives obtained with a recursive algorithmic approach for response theory.

Convergence of coupled cluster perturbation theory.

The convergence of a recently proposed coupled cluster (CC) family of perturbation series [J. J. Eriksen et al., J. Chem. Phys. 140, 064108 (2014)], in which the energetic difference between two CC models-a low-level parent and a high-level target model-is expanded in orders of the Møller-Plesset (MP) fluctuation potential, is investigated for four prototypical closed-shell systems (Ne, singlet CH2, distorted HF, and F-) in standard and augmented basis sets. In these investigations, energy corrections of the various series have been calculated to high orders and their convergence radii have been determined by probing for possible front- and back-door intruder states, the existence of which would make the series divergent. In summary, we conclude how it is primarily the choice of the target state, and not the choice of the parent state, which ultimately governs the convergence behavior of a given series. For example, restricting the target state to, say, triple or quadruple excitations might remove intruders present in series which target the full configuration interaction limit, such as the standard MP series. Furthermore, we find that whereas a CC perturbation series might converge within standard correlation consistent basis sets, it may start to diverge whenever these become augmented by diffuse functions, similar to the MP case. However, unlike for the MP case, such potential divergences are not found to invalidate the practical use of the low-order corrections of the CC perturbation series.