Tensor hypercontracted ppRPA: Reducing the cost of the particle-particle random phase approximation from O(r 6) to O(r 4)

Tensor hypercontraction for fully self-consistent imaginary-time GF2 and GWSOX methods: Theory, implementation, and role of the Green's function second-order exchange for intermolecular interactions

We present an efficient MPI-parallel algorithm and its implementation for evaluating the self-consistent correlated second-order exchange term (SOX), which is employed as a correction to the fully self-consistent GW scheme called scGWSOX (GW plus the SOX term iterated to achieve full Green's function self-consistency). Due to the application of the tensor hypercontraction (THC) in our computational procedure, the scaling of the evaluation of scGWSOX is reduced from O(nτnAO5) to O(nτN2nAO2). This fully MPI-parallel and THC-adapted approach enabled us to conduct the largest fully self-consistent scGWSOX calculations with over 1100 atomic orbitals with only negligible errors attributed to THC fitting. Utilizing our THC implementation for scGW, scGF2, and scGWSOX, we evaluated energies of intermolecular interactions. This approach allowed us to circumvent issues related to reference dependence and ambiguity in energy evaluation, which are common challenges in non-self-consistent calculations. We demonstrate that scGW exhibits a slight overbinding tendency for large systems, contrary to the underbinding observed with non-self-consistent RPA. Conversely, scGWSOX exhibits a slight underbinding tendency for such systems. This behavior is both physical and systematic and is caused by exclusion-principle violating diagrams or corresponding corrections. Our analysis elucidates the role played by these different diagrams, which is crucial for the construction of rigorous, accurate, and systematic methods. Finally, we explicitly show that all perturbative fully self-consistent Green's function methods are size-extensive and size-consistent.

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (Nk) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, NISDF) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of NISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. We perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

When Gold Is Not Enough: Platinum Standard of Quantum Chemistry with N7 Cost

In this paper, we extend the rank-reduced coupled-cluster formalism to the calculation of non-iterative energy corrections due to quadruple excitations. There are two major components of the proposed formalism. The first is an approximate compression of the quadruple excitation amplitudes using the Tucker format. The second is a modified functional used for the evaluation of the corrections which gives exactly the same results for the exact amplitudes, but is less susceptible to errors resulting from the aforementioned compression. We show, both theoretically and numerically, that the computational cost of the proposed method scales as the seventh power of the system size. Using reference results for a set of small molecules, the method is calibrated to deliver relative accuracy of a few percent in energy corrections. To illustrate the potential of the theory, we calculate the isomerization energy of ortho/meta benzyne (C₆H₄) and the barrier height for the Cope rearrangement in bullvalene (C₁₀H₁₀). The method retains a near-black-box nature of the conventional coupled-cluster formalism and depends on only one additional parameter that controls the accuracy.

Tensor-Hypercontracted MP2 First Derivatives: Runtime and Memory Efficient Computation of Hyperfine Coupling Constants

We employ our recently introduced tensor-hypercontracted (THC) second-order Møller-Plesset perturbation theory (MP2) method [Bangerter, F. H., Glasbrenner, M., Ochsenfeld, C. J. Chem. Theory Comput. 2021, 17, 211-221] for the computation of hyperfine coupling constants (HFCCs). The implementation leverages the tensor structure of the THC factorized electron repulsion integrals for an efficient formation of the integral-based intermediates. The computational complexity of the most expensive and formally quintic scaling exchange-like contribution is reduced to effectively subquadratic, by making use of the intrinsic, exponentially decaying coupling between tensor indices through screening based on natural blocking. Overall, this yields an effective subquadratic scaling with a low prefactor for the presented THC-based AO-MP2 method for the computation of isotropic HFCCs on DNA fragments with up to 500 atoms and 5000 basis functions. Furthermore, the implementation achieves considerable speedups with up to a factor of roughly 600-1000 compared to previous implementations [Vogler, S., Ludwig, M., Maurer, M., Ochsenfeld, C. J. Chem. Phys. 2017, 147, 024101] for medium-sized organic radicals, while also significantly reducing storage requirements.

Static and dynamic Bethe-Salpeter equations in the T-matrix approximation

While the well-established GW approximation corresponds to a resummation of the direct ring diagrams and is particularly well suited for weakly correlated systems, the T-matrix approximation does sum ladder diagrams up to infinity and is supposedly more appropriate in the presence of strong correlation. Here, we derive and implement, for the first time, the static and dynamic Bethe-Salpeter equations when one considers T-matrix quasiparticle energies and a T-matrix-based kernel. The performance of the static scheme and its perturbative dynamical correction are assessed by computing the neutral excited states of molecular systems. A comparison with more conventional schemes as well as other wave function methods is also reported. Our results suggest that the T-matrix-based formalism performs best in few-electron systems where the electron density remains low.

Near-Exact CCSDT Energetics from Rank-Reduced Formalism Supplemented by Non-iterative Corrections

We introduce a non-iterative energy correction, added on top of the rank-reduced coupled-cluster method with single, double, and triple substitutions, that accounts for excitations excluded from the parent triple excitation subspace. The formula for the correction is derived by employing the coupled-cluster Lagrangian formalism, with an additional assumption that the parent excitation subspace is closed under the action of the Fock operator. Owing to the rank-reduced form of the triple excitation amplitudes tensor, the computational cost of evaluating the correction scales as N⁷, where N is the system size. The accuracy and computational efficiency of the proposed method is assessed for both total and relative correlation energies. We show that the non-iterative correction can fulfill two separate roles. If the accuracy level of a fraction of kJ/mol is sufficient for a given system, the correction significantly reduces the dimension of the parent triple excitation subspace needed in the iterative part of the calculations. Simultaneously, it enables reproducing the exact CCSDT results to an accuracy level below 0.1 kJ/mol, with a larger, yet still reasonable, dimension of the parent excitation subspace. This typically can be achieved at a computational cost only several times larger than required for the CCSD(T) method. The proposed method retains the black-box features of the single-reference coupled-cluster theory; the dimension of the parent excitation subspace remains the only additional parameter that has to be specified.

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.

Hole-hole Tamm-Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole-hole Tamm-Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole-hole counterpart to the more established particle-particle TDA (pp-TDA) method, both of which are derived from the particle-particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. The treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. We show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states-particularly those with both low-lying ππ^* and nπ^* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

Accurate Quasiparticle Spectra from the T-Matrix Self-Energy and the Particle-Particle Random Phase Approximation

The GW self-energy, especially G₀W₀ based on the particle-hole random phase approximation (phRPA), is widely used to study quasiparticle (QP) energies. Motivated by the desirable features of the particle-particle (pp) RPA compared to the conventional phRPA, we explore the pp counterpart of GW, that is, the T-matrix self-energy, formulated with the eigenvectors and eigenvalues of the ppRPA matrix. We demonstrate the accuracy of the T-matrix method for molecular QP energies, highlighting the importance of the pp channel for calculating QP spectra.

Nature of ground and electronic excited states of higher acenes

Higher acenes have drawn much attention as promising organic semiconductors with versatile electronic properties. However, the nature of their ground state and electronic excited states is still not fully clear. Their unusual chemical reactivity and instability are the main obstacles for experimental studies, and the potentially prominent diradical character, which might require a multireference description in such large systems, hinders theoretical investigations. Here, we provide a detailed answer with the particle-particle random-phase approximation calculation. The (1)Ag ground states of acenes up to decacene are on the closed-shell side of the diradical continuum, whereas the ground state of undecacene and dodecacene tilts more to the open-shell side with a growing polyradical character. The ground state of all acenes has covalent nature with respect to both short and long axes. The lowest triplet state (3)B2u is always above the singlet ground state even though the energy gap could be vanishingly small in the polyacene limit. The bright singlet excited state (1)B2u is a zwitterionic state to the short axis. The excited (1)Ag state gradually switches from a double-excitation state to another zwitterionic state to the short axis, but always keeps its covalent nature to the long axis. An energy crossing between the (1)B2u and excited (1)Ag states happens between hexacene and heptacene. Further energetic consideration suggests that higher acenes are likely to undergo singlet fission with a low photovoltaic efficiency; however, the efficiency might be improved if a singlet fission into multiple triplets could be achieved.