An "ultimate" coupled cluster method based entirely on T2.

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 T2 and its products that can effectively emulate more "complete" methods that explicitly consider higher-rank, T2m, operators.

Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry.

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 with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy.

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory.

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 dynamic 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 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism.

A new "gold standard": Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing.

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 for 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 the effects of triple excitations using considerations from many-body perturbation theory (MBPT).

Perspective on Coupled-cluster Theory. The evolution toward simplicity in quantum chemistry.

Coupled-cluster theory has revolutionized quantum chemistry. It has provided the framework to effectively solve the problem of electron correlation, the main focus of the field for over 60 years. This has enabled quantum chemistry to provide predictive quality results for most quantities of interest that are obtainable from first-principle calculations.

A comparison of QTP functionals against coupled-cluster methods for EAs of small organic molecules.

EA-EOM-CCSD electron affinities and LUMO energies of various Kohn-Sham density functional theory (DFT) methods are calculated for an a priori IP benchmark set of 64 small, closed-shell molecules. The purpose of these calculations was to investigate whether the QTP KS-DFT functionals can emulate EA-EOM-CC with only a mean-field approximation. We show that the accuracy of DFT-relative to CCSD-improves significantly when elements of correlated orbital theory are introduced into the parameterization to define the QTP family of functionals.

Excited states with pair coupled cluster doubles tailored coupled cluster theory.

It is known that some non-dynamic effects of electron correlation can be included in coupled cluster theory using a tailoring technique that separates the effects of non-dynamic and dynamic correlations. Recently, the simple pCCD (pair coupled cluster doubles) wavefunction was shown to provide good results for some non-dynamic correlation problems, such as bond-breaking, in a spin-adapted way with no active space selection. In this paper, we report a study of excited states using "tailored coupled cluster singles and doubles," to attempt to use pCCD as a kernel for more complete coupled-cluster singles and doubles (CCSD) results for excited states.

On the iterative diagonalization of matrices in quantum chemistry: Reconciling preconditioner design with Brillouin-Wigner perturbation theory.

Iterative diagonalization of large matrices to search for a subset of eigenvalues that may be of interest has become routine throughout the field of quantum chemistry. Lanczos and Davidson algorithms hold a monopoly, in particular, owing to their excellent performance on diagonally dominant matrices. However, if the eigenvalues happen to be clustered inside overlapping Gershgorin disks, the convergence rate of both strategies can be noticeably degraded.

Helical Organic and Inorganic Polymers.

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set electron-correlated methods on an equal footing with molecules. This article introduces an second-order many-body Green's function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results.

Doubly Excited States with Modified Coupled Cluster Models: A Reliable Compromise between Cost and Accuracy?

To treat doubly excited states, the treatment of triple excitations is considered necessary in the framework of equation-of-motion coupled cluster (EOM-CC) methods. We investigate models without explicit triples and seek quantitative measure for the performance of EOM based on CC with singles and doubles (CCSD) or modified CCSD (Distinguishable Cluster Approximation) approaches for states with predominant double excitation character. We also test the efficacy of including triples in perturbative manner through EOM-CCSD(T) and in an iterative way through EOM-CCSDT-3 method.

Density functionals for core excitations.

The core excitation energies and related principal ionization energies are obtained for selected molecules using several density functionals and compared with benchmark equation-of-motion coupled cluster (EOM-CC) results. Both time-dependent and time-independent formulations of excitation spectra in the time-dependent density functional theory and the EOM-CC are employed to obtain excited states that are not always easily accessible with the time-independent method. Among those functionals, we find that the QTP(00) functional, which is only parameterized to reproduce the five IPs of water, provides excellent core IPs and core excitation energies, consistently yielding better excitation and ionization energies.

Examining fundamental and excitation gaps at the thermodynamic limit: A combined (QTP) DFT and coupled cluster study on trans-polyacetylene and polyacene.

Interest in ab initio property prediction of π-conjugated polymers for technological applications places significant demand on "cost-effective" and conceptual computational methods, particularly effective, one-particle theories. This is particularly relevant in the case of Kohn-Sham Density Functional Theory (KS-DFT) and its new competitors that arise from correlated orbital theory, the latter defining the QTP family of DFT functionals. This study presents large, ab initio equation of motion-coupled cluster calculations using the massively parallel ACESIII to target the fundamental bandgap of two prototypical organic polymers, trans-polyacetylene (tPA) and polyacene (Ac), and provides an assessment of the new quantum theory project (QTP) functionals for this problem.

The intermediate state approach for doubly excited dark states in EOM-coupled-cluster theory.

Solution of dark, doubly excited states using equation-of-motion coupled-cluster (EOM-CC) usually equires at least triple excitations or even quadruples beyond the standard singles and doubles (EOM-CCSD) for an appropriate treatment. A new route to obtain these doubly excited states using EOM-CCSD is demonstrated. Traditionally, EOM-CC is performed on a closed shell reference state that has a well-described single reference CC wavefunction.

Benchmarking isotropic hyperfine coupling constants using (QTP) DFT functionals and coupled cluster theory.

Significant effort has been devoted to benchmarking isotropic hyperfine coupling constants for both wavefunction and density-based approaches in recent years, as accurate theoretical predictions aid the fitting of experimental model Hamiltonians. However, literature examining the predictive quality of a Density Functional Theory (DFT) functional abiding by the Bartlett IP condition is absent. In an attempt to rectify this, we report isotropic hyperfine coupling constant predictions of 24 commonly used DFT functionals on a total of 56 radicals, with the intent of exploring the successes and failures of the Quantum Theory Project (QTP) line of DFT functionals (i.

Equation of motion coupled-cluster study of core excitation spectra II: Beyond the dipole approximation.

We present the time-independent (TI) and time-dependent (TD) equation of motion coupled-cluster (EOM-CC) oscillator strengths not limited to those obtained by the dipole approximation. For the conventional TI-EOM-CC, we implement all the terms in the multipole expansion through second order that contributes to the oscillator strength. These include contributions such as magnetic dipole, electric quadrupole, electric octupole, and magnetic quadrupole.

The Devil's Triangle of Kohn-Sham density functional theory and excited states.

Exchange-correlation (XC) functionals from Density Functional Theory (DFT) developed under the rigorous arguments of Correlated Orbital Theory (COT) address the Devil's Triangle of prominent errors in Kohn-Sham (KS) DFT. At the foundation of this triangle lie the incorrect one-particle spectrum, the lack of integer discontinuity, and the self-interaction error. At the top level, these failures manifest themselves in incorrect charge transfer and Rydberg excitation energies, along with poor activation barriers.

Index of multi-determinantal and multi-reference character in coupled-cluster theory.

A full configuration interaction calculation (FCI) ultimately defines the innate molecular orbital description of a molecule. Its density matrix and the natural orbitals obtained from it quantify the difference between having N-dominantly occupied orbitals in a reference determinant for a wavefunction to describe N-correlated electrons and how many of those N-electrons are left to the remaining virtual orbitals. The latter provides a measure of the multi-determinantal character (MDC) required to be in a wavefunction.

A route to improving RPA excitation energies through its connection to equation-of-motion coupled cluster theory.

We revisit the connection between equation-of-motion coupled cluster (EOM-CC) and random phase approximation (RPA) explored recently by Berkelbach [J. Chem. Phys.

Advanced concepts in electronic structure (ACES) software programs.

The advanced concepts in electronic structure (ACES) programs are products of the Bartlett research group at the University of Florida. They consist of ACES II, which is serial, and ACES III and Aces4, which are massively parallel. All three programs are publically available free of charge.

Adventures in DFT by a wavefunction theorist.

The attraction density functional theory (DFT) has for electronic structure theory is that it is easier to do computationally than ab initio, correlated wavefunction methods, due to its effective one-particle structure. On the contrary, ab initio theorists insist on the ability to converge to the right answer in appropriate limits, but this requires a treatment of the reduced two-particle density matrix. DFT avoids that by appealing to an "existence" theorem (not a constructive one) that all its effects are subsummed into a DFT functional of the one-particle density.

Similarity-transformed equation-of-motion coupled-cluster singles and doubles method with spin-orbit effects for excited states.

The similarity transformed equation-of-motion coupled-cluster method (STEOM-CCSD) for excited states is extended to treat spin-orbit coupling interactions (SOIs) using the spin-orbit mean field approximation of the Breit-Pauli Hamiltonian. Two possible schemes to include the spin-orbit splittings of excited states within the STEOM-CCSD model are formulated. They are identified as "diagonalize-then-perturb" and "perturb-then-diagonalize" approaches.

Vertical valence ionization potential benchmarks from equation-of-motion coupled cluster theory and QTP functionals.

The ionization potential (IP) of a molecule quantifies the energy required to remove an electron from the system. As such, it is a fundamental quantity in the context of redox chemistry, charge transfer, and molecular electronics. The accurate theoretical prediction of this property is therefore highly desirable for virtual materials design.

Perturbation Improved Natural Linear-Scaled Coupled-Cluster Method and Its Application to Conformational Analysis.

The fragment-based coupled-cluster (CC) theory utilizing the transferable functional groups through natural localized molecular orbital (NLMO), that is, the natural linear-scaled coupled-cluster (NLSCC) has been further developed to take the extra-fragment interactions into account. The correction to the interaction energies sacrificed during the fragmentation process for the previous NLSCC method is computed by a computationally efficient perturbation theory that maintains the original scaling. The new linear-scaled coupled-cluster for the singles and doubles (CCSD) method is applied to the analysis of relative energies of delicate conformational problems of polypeptides.

Low scaling EOM-CCSD and EOM-MBPT(2) method with natural transition orbitals.

A low-scaling method is presented for the equation-of-motion coupled-cluster theory with single and double (EOM-CCSD) excitations and its second-order many-body perturbation theory [EOM-MBPT(2)] approximations. For a simple description of an excited state, the particular orbitals, and , are selected from the natural transition orbitals (NTOs, ), where and refer to NTO occupied and virtual orbital indices. They are chosen based on the largest eigenvalues of the transition density matrix.