Search

Traditionally, excitation energies in coupled-cluster (CC) theory have been calculated by solving the CC Jacobian eigenvalue equation. However, based on our recent work [Jørgensen et al., Sci. Adv. 10, eadn3454 (2024)], we propose a reformulation of the calculation of excitation energies where excitation energies are determined as a conventional molecular property. To this end, we introduce an excitation-energy function that depends on the CC Jacobian and the right and left eigenvectors for the Jacobian eigenvalue problem. This excitation-energy function is variational with respect to the right and left eigenvectors but not with respect to the cluster amplitudes. Instead, the cluster amplitudes satisfy the cluster-amplitude equations, and we set up an excitation-energy Lagrangian by adding to the excitation-energy function the cluster-amplitude equations with an undetermined multiplier for each cluster-amplitude constraint. The excitation-energy Lagrangian is variational in all its parameters. Based on the variational property of the Lagrangian, we have determined two quadratically convergent excitation-energy series: the total-order cluster-perturbation (tCP) and variational cluster-perturbation (vCP) excitation-energy series. Calculations of the excitation energies of three small molecules have shown that the vCP series is to be preferred over the tCP series. The test calculations have been carried out for CPS(D) expansions [targeting the CC singles-and-doubles (CCSD) wave function from the CC singles wave function] and the CPSD(T) expansion [targeting the CC singles-doubles-triples (CCSDT) wave function from the CCSD wave function]. For the S(D) and SD(T) orbital excitation space calculations, we obtain in the second vCP iteration excitation energies with a mean deviation from CCSD excitation energies of about 0.04 eV for the S(D) orbital spaces, and for the SD(T) orbital space calculation, we obtain a mean deviation from the CCSDT excitation energies of 0.001 eV. [ABSTRACT FROM AUTHOR]

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.

We describe an efficient implementation of cluster perturbation and Møller–Plesset Lagrangian energy series through the fifth order that targets the coupled cluster singles and doubles energy utilizing the resolution of the identity approximation. We illustrate the computational performance of the implementation by performing ground state energy calculations on systems with up to 1200 basis functions using a single node and by comparison to conventional coupled cluster singles and doubles calculations. We further show that our hybrid message passing interface/open multiprocessing parallel implementation that also utilizes graphical processing units can be used to obtain fifth order energies on systems with almost 1200 basis functions with a 90 min "time to solution" running on Frontier at Oak Ridge National Laboratory. [ABSTRACT FROM AUTHOR]

In this paper, we develop and analyze a number of perturbation series that target the coupled cluster singles and doubles (CCSD) ground state energy. We show how classical Møller–Plesset perturbation theory series can be restructured to target the CCSD energy based on a reference CCS calculation and how the corresponding cluster perturbation series differs from the classical Møller–Plesset perturbation series. Subsequently, we reformulate these series using the coupled cluster Lagrangian framework to obtain series, where fourth and fifth order energies are determined only using parameters through second order. To test the methods, we perform a series of test calculations on molecular photoswitches of both total energies and reaction energies. We find that the fifth order reaction energies are of CCSD quality and that they are of comparable accuracy to state-of-the-art approximations to the CCSD energy based on local pair natural orbitals. The advantage of the present approach over local correlation methods is the absence of user defined threshold parameters for neglecting or approximating contributions to the correlation energy. Fixed threshold parameters lead to discontinuous energy surfaces, although this effect is often small enough to be ignored, but the present approach has a differentiable energy that will facilitate derivation and implementation of gradients and higher derivatives. A further advantage is that the calculation of the perturbation correction is non-iterative and can, therefore, be calculated in parallel, leading to a short time-to-solution. [ABSTRACT FROM AUTHOR]

A new large-scale parallel multiconfigurational self-consistent field (MCSCF) implementation in the open-source NWChem computational chemistry code is presented. The generalized active space approach is used to partition large configuration interaction (CI) vectors and generate a sufficient number of batches that can be distributed to the available cores. Massively parallel CI calculations with large active spaces can be performed. The new parallel MCSCF implementation is tested for the chromium trimer and for an active space of 20 electrons in 20 orbitals, which can now routinely be performed. Unprecedented CI calculations with an active space of 22 electrons in 22 orbitals for the pentacene systems were performed and a single CI iteration calculation with an active space of 24 electrons in 24 orbitals for the chromium tetramer was possible. The chromium tetramer corresponds to a CI expansion of one trillion Slater determinants (914 058 513 424) and is the largest conventional CI calculation attempted up to date.

A new implementation of the orbital-based two-component relativistic configuration interaction approach is reported and applied to calculations of the electronic g-shifts of three diatomic radicals: AlO, HgF, and PdH. The new implementation augments efficient routines for the calculation of nonrelativistic Hamiltonians with new vectorized routines for the calculation of the action of the one-electron spin–orbit operator and allows efficient calculations for the expansion of generalized active space type. The program makes full use of double group as well as time-reversal symmetry. Particle–hole reorganization of the operators is used to improve the efficiency for expansions with nearly fully occupied orbital spaces. The flexibility of the algorithm and program is used to investigate the convergence of electronic g-shifts for the three diatomic radicals as functions of the active space, states included in the orbital optimization, and excitation levels. It was possible to converge to the valence limits within a few percent using expansions containing up to quadruple excitations. However, when excitations from the core orbitals were added, it was not possible to demonstrate convergence to within a few percent with expansions containing at most 10 × 109 determinants. [ABSTRACT FROM AUTHOR]

We present here a massively parallel implementation of the recently developed CPS(D-3) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model. [ABSTRACT FROM AUTHOR]

The convergence of the recently developed cluster perturbation (CP) expansions [Pawlowski et al., J. Chem. Phys. 150, 134108 (2019)] is analyzed with the double purpose of developing the mathematical tools and concepts needed to describe these expansions at general order and to identify the factors that define the rate of convergence of CP series. To this end, the CP energy, amplitude, and Lagrangian multiplier equations as a function of the perturbation strength are developed. By determining the critical points, defined as the perturbation strengths for which the Jacobian becomes singular, the rate of convergence and the intruder and critical states are determined for five small molecules: BH, CO, H2O, NH3, and HF. To describe the patterns of convergence for these expansions at orders lower than the high-order asymptotic limit, a model is developed where the perturbation corrections arise from two critical points. It is shown that this model allows for rationalization of the behavior of the perturbation corrections at much lower order than required for the onset of the asymptotic convergence. For the H2O, CO, and HF molecules, the pattern and rate of convergence are defined by critical states where the Fock-operator underestimates the excitation energies, whereas the pattern and rate of convergence for BH are defined by critical states where the Fock-operator overestimates the excitation energy. For the NH3 molecule, both forms of critical points are required to describe the convergence behavior up to at least order 25. [ABSTRACT FROM AUTHOR]

We have extended cluster perturbation (CP) theory to comprehend the Lagrangian framework of coupled cluster (CC) theory and derived the CP Lagrangian energy series (LCP) where the 2n + 1/2n + 2 rules for the cluster amplitudes and multipliers are used to get the energy corrections. We have also developed the variational CP ( L C P ) series, where the total cluster amplitudes and multipliers are determined through the same orders as in the LCP series, but the energy is obtained by inserting the total cluster amplitudes and multipliers in the Lagrangian. The energies of the L C P series have errors that are bilinear in the errors of the total cluster amplitudes and multipliers. Test calculations have been performed for S(D) and SD(T) orbital excitation spaces. With the exception of molecular systems that have a low lying doubly excited state compared to the electronic ground state configuration, we find that the fourth order models L C P S (D−4), L C P S D (T−4), and LCPSD(T−4) give energies of CC target state quality. For the L C P S (D−4) model, CC target state quality is obtained as the L C P S (D−4) calculation determines more than 99.7% of the coupled cluster singles and doubles (CCSD) correlation energy as the numerical deviations of the L C P S (D−4) energy from the CCSD energy were more than an order of magnitude smaller than the triples correlation contribution. For the L C P S D (T−4) and LCPSD(T−4) models, CC target state quality was obtained, given that the L C P S D (T−4) and LCPSD(T−4) calculations recover more than 99% of the coupled cluster singles doubles and triples (CCSDT) correlation contribution and as the numerical deviations of the L C P S D (T−4) and LCPSD(T−4) energies from the CCSDT energy were nearly and order of magnitude smaller than the quadruples correlation contribution. We, thus, suggest that the fourth order models may replace the full target CC models with no or very limited loss of accuracy. [ABSTRACT FROM AUTHOR]

We have extended cluster perturbation (CP) theory to comprehend the calculation of first order properties (FOPs). We have determined CP FOP series where FOPs are determined as a first energy derivative and also where the FOPs are determined as a generalized expectation value of the external perturbation operator over the coupled cluster state and its biorthonormal multiplier state. For S(D) orbital excitation spaces, we find that the CP series for FOPs that are determined as a first derivative, in general, in second order have errors of a few percent in the singles and doubles correlation contribution relative to the targeted coupled cluster (CC) results. For a SD(T) orbital excitation space, we find that the CP series for FOPs determined as a generalized expectation value in second order have errors of about ten percent in the triples correlation contribution relative to the targeted CC results. These second order models, therefore, constitute viable alternatives for determining high quality FOPs. [ABSTRACT FROM AUTHOR]

A simple two-state model has previously been shown to be able to describe and rationalize the convergence of the most common perturbation method for including electron correlation, the Møller-Plesset expansion. In particular, this simple model has been able to predict the convergence rate and the form of the higher-order corrections for typical Møller-Plesset expansions of the correlation energy. In this paper, the convergence of nondegenerate perturbation expansions in the two-state model is analyzed in detail for a general form of two-state perturbation expansion by examining the analytic expressions of the corrections and series of the values of the corrections for various choices of the perturbation. The previous analysis that covered only a single form of the perturbation is thereby generalized to arbitrary forms of the perturbation. It is shown that the convergence may be described in terms of four characteristics: archetype, rate of convergence, length of recurring period, and sign pattern. The archetype defines the overall form of a plot of the energy-corrections, and the remaining characteristics specify details of the archetype. For symmetric (Hermitian) perturbations, five archetypes are observed: zigzag, interspersed zigzag, triadic, ripples, and geometric. Two additional archetypes are obtained for an asymmetric perturbation: zigzag-geometric and convex-geometric. For symmetric perturbations, each archetype has a distinctive pattern that recurs with a period which depends on the perturbation parameters, whereas no such recurrence exists for asymmetric perturbations from a series of numerical corrections. The obtained relations between the form of a two-state perturbation and the energy corrections allow us to obtain additional insights into the convergence behavior of the Møller-Plesset and other forms of perturbation expansions. This is demonstrated by analyzing several diverging or slowly converging perturbation expansions of ground state and excitation energies. It is demonstrated that the higher-order corrections of these expansions can be described using the two-state model and each expansion can therefore be described in terms of an archetype and the other three characteristics. Examples of all archetypes except the zigzag and convex-geometric archetypes are given. For each example, it is shown how the characteristics may be extracted from the higher-order corrections and used to identify the term in the perturbation that is the cause of the observed slow convergence or divergence. [ABSTRACT FROM AUTHOR]

Cluster perturbation (CP) theory was developed in Paper I [F. Pawłowski et al., J. Chem. Phys. 150, 134108 (2019)] for a coupled cluster (CC) target state and is extended in this paper to comprehend a cluster linear (CL) target state, for which the embedding of a CC parent state in the target excitation space is described using a linear parametrization. The theory is developed for determining the energy and molecular properties for a CL state. When CP theory is applied to a CL target state, a series of corrections is determined in orders of the CC parent-state similarity-transformed fluctuation potential, where the zeroth-order term is the energy or molecular property of the CC parent state and where the series formally converges to the energy or molecular property of the CL target state. The determination of energies and molecular properties is simpler for a CL state than for a CC state because the CL state is linearly parametrized. The amplitude equations are quadratic for a CL target state, while quartic for a CC target state, and molecular property expressions for a CL target state have the same simple structure as for a configuration interaction state. The linear parametrization introduces non-size-extensive contributions in the energy and molecular property expressions. However, since the linear parametrization describes the embedding of the CC parent state in the target excitation space, the energy and molecular properties for a CL state are weakly size-extensive. For the energy, weak size-extensivity means that non-size-extensive contributions enter in sixth and higher orders in the CP energy series, whereas for molecular properties, weak size-extensivity means that non-size-extensive contributions enter in second and higher orders. Weak size-extensivity therefore has a little or vanishing effect on calculated energies or molecular properties. The determination of the CP energy and molecular property corrections does not require that amplitude or response equations are solved explicitly for the target state and it becomes computationally tractable to use low-order corrections from these series to obtain energies and molecular properties of CL target state quality. For three simple molecules, HF, N2, and CH2, the accuracy of the CL approach for ground-state energies is tested using a parent state including single and double excitations (i.e., the CC singles-and-doubles state, CCSD) and a target state that includes triple excitations. It is found that the size-extensive fifth-order CL energies deviate by less than 0.0001 hartree from the energies of a target CC that includes triple excitations (i.e., the CC singles-doubles-and-triples state, CCSDT). CP theory with a CL target state therefore becomes a very attractive replacement of standard CC theory for high-accuracy energy and molecular property calculations, in which triple and higher excitation levels are considered. [ABSTRACT FROM AUTHOR]