Communication: The distinguishable cluster approximation

Rank-Reduced Equation-of-Motion Coupled Cluster Triples: an Accurate and Affordable Way of Calculating Electronic Excitation Energies

In the present work, we report an implementation of the rank-reduced equation-of-motion coupled cluster method with approximate triple excitations (RR-EOM-CC3). The proposed variant relies on tensor decomposition techniques in order to alleviate the high cost of computing and manipulating the triply excited amplitudes. In the RR-EOM-CC3 method, both ground-state and excited-state triple-excitation amplitudes are compressed according to the Tucker-3 format. This enables factorization of the working equations such that the formal scaling of the method is reduced to N6, where N is the system size. An additional advantage of our method is the fact that the accuracy can be strictly controlled by proper choice of two parameters defining sizes of triple-excitation subspaces in the Tucker decomposition for the ground and excited states. Optimal strategies of selecting these parameters are discussed. The developed method has been tested in a series of calculations of electronic excitation energies and compared to its canonical EOM-CC3 counterpart. Errors several times smaller than the inherent error of the canonical EOM-CC3 method (in comparison to FCI) are straightforward to achieve. This conclusion holds both for valence states dominated by single excitations and for states with pronounced doubly excited character. Taking advantage of the decreased scaling, we demonstrate substantial computational costs reductions (in comparison with the canonical EOM-CC3) in the case of two large molecules - l-proline and heptazine. This illustrates the usefulness of the RR-EOM-CC3 method for accurate determination of excitation energies of large molecules.

Orbital optimisation in xTC transcorrelated methods

We present a combination of the bi-orthogonal orbital optimisation framework with the recently introduced xTC version of transcorrelation. This allows us to implement non-iterative perturbation based methods on top of the transcorrelated Hamiltonian. Additionally, the orbital optimisation influences results of other truncated methods, such as the distinguishable cluster with singles and doubles. The accuracy of these methods in comparison to standard xTC methods is demonstrated, and the advantages and disadvantages of the orbital optimisation are discussed.

Quantum Information Orbitals (QIO): Unveiling Intrinsic Many-Body Complexity by Compressing Single-Body Triviality

The simultaneous treatment of static and dynamic correlations in strongly correlated electron systems is a critical challenge. In particular, finding a universal scheme for identifying a single-particle orbital basis that minimizes the representational complexity of the many-body wave function is a formidable and longstanding problem. As a contribution toward its solution, we show that the total orbital correlation actually reveals and quantifies the intrinsic complexity of the wave function, once it is minimized via orbital rotations. To demonstrate the power of this concept in practice, an iterative scheme is proposed to optimize the orbitals by minimizing the total orbital correlation calculated by the tailored coupled cluster singles and doubles (TCCSD) ansatz. The optimized orbitals enable the limited TCCSD ansatz to capture more nontrivial information on the many-body wave function, indicated by the improved wave function and energy. An initial application of this scheme shows great improvement of TCCSD in predicting the singlet ground state potential energy curves of the strongly correlated C2 and Cr2 molecule.

Toward Reliable Dipole Moments without Single Excitations: The Role of Orbital Rotations and Dynamical Correlation

The dipole moment is a crucial molecular property linked to a molecular system's bond polarity and overall electronic structure. To that end, the electronic dipole moment, which results from the electron density of a system, is often used to assess the accuracy and reliability of new electronic structure methods. This work analyses electronic dipole moments computed with the pair coupled cluster doubles (pCCD) ansätze and its linearized coupled cluster (pCCD-LCC) corrections using the canonical Hartree-Fock and pCCD-optimized (localized) orbital bases. The accuracy of pCCD-based dipole moments is assessed against experimental and CCSD(T) reference values using relaxed and unrelaxed density matrices and different basis set sizes. Our test set comprises molecules of various bonding patterns and electronic structures, exposing pCCD-based methods to a wide range of electron correlation effects. Additionally, we investigate the performance of pCCD-in-DFT dipole moments of some model complexes. Finally, our work indicates the importance of orbital relaxation in the pCCD model and shows the limitations of the linearized couple cluster corrections in predicting electronic dipole moments of multiple-bonded systems. Most importantly, pCCD with a linearized CCD correction can reproduce the dipole moment surfaces in singly bonded molecules, which are comparable to the multireference ones.

State-Specific Coupled-Cluster Methods for Excited States

We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.

Two determinant distinguishable cluster

A two reference determinant version of the distinguishable cluster with singles and doubles (DCSD) has been developed. We have implemented the two determinant distinguishable cluster (2D-DCSD) and the corresponding traditional 2D-CCSD method in a new open-source package written in Julia called ElemCo.jl. The methods were benchmarked on singlet and triplet excited states of valence and Rydberg character, as well as for singlet-triplet gaps of diradicals. It is demonstrated that the distinguishable cluster approximation improves the accuracy of 2D-CCSD.

Comprehensive quantum chemical analysis of the (ro)vibrational spectrum of thiirane and its deuterated isotopologue

The (ro)vibrational spectra of thiirane, c-C2H4S, and its fully deuterated isotopologue, c-C2D4S, have been studied by means of vibrational configuration interaction theory, VCI, its incremental variant, iVCI, and subsequent variational rovibrational calculations, RVCI, which rely on multidimensional potential energy surfaces of coupled-cluster quality including up to four-mode coupling terms. Accurate geometrical parameters, fundamental vibrational transitions and first overtones, rovibrational spectra and rotational spectroscopic constants have been determined from these calculations and were compared with experimental results whenever available. A number of tentative misassignments in the vibrational spectra could be resolved and most results for the deuterated thiirane are high-level predictions, which may guide experiments to come. Besides this, a new implementation of infrared intensities within the iVCI framework has been tested for the transitions of the title compounds and are compared with results obtained from standard VCI calculations.

Highly Accurate Electronic Structure of Metallic Solids from Coupled-Cluster Theory with Nonperturbative Triple Excitations

Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T))-often considered the "gold standard" of main-group quantum chemistry-is inapplicable to three-dimensional metals due to an infrared divergence, preventing its application to many important problems in materials science. We study the full, nonperturbative inclusion of triple excitations (CCSDT) and propose a new, iterative method, which we call ring-CCSDT, that resums the essential triple excitations with the same N^{7} run-time scaling as CCSD(T). CCSDT and ring-CCSDT are used to calculate the correlation energy of the uniform electron gas at metallic densities and the structural properties of solid lithium. Inclusion of connected triple excitations is shown to be essential to achieving high accuracy. We also investigate semiempirical CC methods based on spin-component scaling and the distinguishable cluster approximation and find that they enhance the accuracy of their parent ab initio methods.

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Mid-infrared Laser Spectroscopy of Jet-Cooled Formic Acid Trimer: Mode-Dependent Line Broadening in the C-O Stretching Region

Building on recent progress in the vibrational spectroscopy of the formic acid trimer, we present the first high-resolution measurements of the jet-cooled laser absorption spectrum of (HCOOH)3. The spectra of the lowest- and highest-frequency C-O stretching fundamentals are analyzed whereas the third band is not observed, complicated by monomer and dimer absorptions at 1219 cm-1 (8.2 μm). Vibration-rotation parameters are obtained for the band at 1172.31512(68) cm-1 whereas the C-O stretch at 1246.33(5) cm-1 exhibits a significantly larger breadth, allowing only resolution of the coarse PQR structure. Vibrational predissociation can be ruled out, and intramolecular vibrational redistribution mechanisms are discussed, particularly coupling to the concerted proton exchange within the cyclic dimer subunit. Ultimately, the question remains open. The prospects of high-resolution measurements of other trimer bands or isotope substitution experiments, which might assist in revealing the mode-specificity of the underlying broadening mechanisms, are discussed.

Quantum chemical rovibrational analysis of aminoborane and its isotopologues

Aminoborane, H₂ NBH₂ and its isotopologues, H₂ N¹⁰ BH₂ , D₂ NBD₂ , and D₂ N¹⁰ BD₂ , have been studied by high-level ab initio methods. All calculations rely on multidimensional potential energy surfaces and dipole moment surfaces including high-order mode coupling terms, which have been obtained from electronic structure calculations at the level of explicitly correlated coupled-cluster theory, CCSD(T)-F12, or the distinguishable cluster approximation, DCSD. Subsequent vibrational structure calculations based on second-order vibrational perturbation theory, VPT2, and vibrational configuration interaction theory, VCI, were used to determine rotational constants, centrifugal distortion constants, vibrationally averaged geometrical parameters and (an)harmonic vibrational frequencies. The impact of core-correlation effects is discussed in detail. Rovibrational VCI calculations were used to simulate the gas phase spectra of these species and an in-depth analysis of the ν₇ band of aminoborane is provided. Color-coding is used to reveal the identity of the individual progressions of the rovibrational transitions for this particular mode.

Dark 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. Extended similarity transformed EOM-CCSD(EXT-STEOM-CCSD) method is also tested and provides superior quality results at comparatively low cost. We use the QUEST2 benchmark set of double excitations proposed by Loos et al. [ J. Chem. Theory Comput.2019, 15, 1939] to investigate the performance of methods such as EOM-CCSD, EOM-DCSD, EXT-STEOM-CCSD, ΔCCSD, and ΔDCSD. We also test a tailored CC approach, ΔpairCCD-TCCSD.

Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Benchmarking Approximate Methods for Adiabatic Spin-State Energy Differences in Fe(II) Complexes

The CASPT2+δMRCI composite approach reported in a companion paper has been extended and used to provide high-quality reference data for a series of adiabatic spin gaps (defined as ΔE = E_quintet - E_singlet) of [Fe^IIL₆]²⁺ complexes (L = CNH, CO, NCH, NH₃, H₂O), either at nonrelativistic level or including scalar relativistic effects. These highly accurate data have been used to evaluate the performance of various more approximate methods. Coupled-cluster theory with singles, doubles, and perturbative triples, CCSD(T), is found to agree well with the new reference data for Werner-type complexes but exhibits larger underestimates by up to 70 kJ/mol for the π-acceptor ligands, due to appreciable static correlation in the low-spin states of these systems. Widely used domain-based local CCSD(T) calculations, DLPNO-CCSD(T), are shown to depend very sensitively on the cutoff values used to construct the localized domains, and standard values are not sufficient. A large number of density functional approximations have been evaluated against the new reference data. The B2PLYP double hybrid gives the smallest deviations, but several functionals from different rungs of the usual ladder hierarchy give mean absolute deviations below 20 kJ/mol. This includes the B97-D semilocal functional, the PBE0* global hybrid with 15% exact-exchange admixture, as well as the local hybrids LH07s-SVWN and LH07t-SVWN. Several further functionals achieve mean absolute errors below 30 kJ/mol (M06L-D4, SSB-D, B97-1-D4, LC-ωPBE-D4, LH12ct-SsirPW92-D4, LH12ct-SsifPW92-D4, LH14t-calPBE-D4, LHJ-HFcal-D4, and several further double hybrids) and thereby also still overall outperform CCSD(T) or uncorrected CASPT2. While exact-exchange admixture is a crucial factor in favoring high-spin states, the present evaluations confirm that other aspects can be important as well. A number of the better-performing functionals underestimate the spin gaps for the π-acceptor ligands but overestimate them for L = NH₃, H₂O. In contrast to a previous suggestion, non-self-consistent density functional theory (DFT) computations on top of Hartree-Fock orbitals are not a promising path to produce accurate spin gaps in such complexes.

Increase of Radiative Forcing through Midinfrared Absorption by Stable CO2 Dimers?

We performed matrix-isolation infrared (MI-IR) spectroscopy of carbon dioxide monomers, CO₂, and dimers, (CO₂)₂, trapped in neon and in air. On the basis of vibration configuration interaction (VCI) calculations accounting for mode coupling and anharmonicity, we identify additional infrared-active bands in the MI-IR spectra due to the (CO₂)₂ dimer. These bands are satellite bands next to the established CO₂ monomer bands, which appear in the infrared window of Earth’s atmosphere at around 4 and 15 μm. In a systematic carbon dioxide mixing ratio study using neon matrixes, we observe a significant fraction of the dimer at mixing ratios above 300 ppm, with a steep increase up to 1000 ppm. In neon matrix, the dimer increases the IR absorbance by about 15% at 400 ppm compared to the monomer absorbance alone. This suggests a high fraction of the (CO₂)₂ dimer in our matrix experiments. In atmospheric conditions, such increased absorbance would significantly amplify radiative forcings and, thus, the greenhouse warming. To enable a comparison of our laboratory experiment with various atmospheric conditions (Earth, Mars, Venus), we compute the thermodynamics of the dimerization accordingly. The dimerization is favored at low temperatures and/or high carbon dioxide partial pressures. Thus, we argue that matrix isolation does not trap the gas composition “as is”. Instead, the gas is precooled to 40 K, where CO₂ dimerizes before being trapped in the matrix, already at very low carbon dioxide partial pressures. In the context of planetary atmospheres, our results improve understanding of the greenhouse effect for planets of rather thick CO₂ atmospheres such as Venus, where a significant fraction of the (CO₂)₂ dimer can be expected. There, the necessity of including the mid-IR absorption by stable (CO₂)₂ dimers in databases used for modeling radiative forcing, such as HITRAN, arises.

Comparison of body definitions for incremental vibrational configuration interaction theory (iVCI)

Within incremental vibrational configuration interaction theory (iVCI), the vibrational state energy is determined by means of a many-body expansion, i.e., it is a sum of terms of increasing order, which allow for an embarrassingly parallel evaluation. The convergence of this expansion depends strongly on the definition of the underlying bodies, which essentially decompose the correlation space into fragments. The different definitions considered here comprise mode-based bodies, excitation level-based bodies, and energy-based bodies. An analysis of the convergence behavior revealed that accounting for resonances within these definitions is mandatory and leads to a substantial improvement of the convergence, that is, the expansions can be truncated at lower orders. Benchmark calculations and systematic comparisons of the different body definitions for a small set of molecules, i.e., ketene, ethene, and diborane, have been conducted to study the overall performance of these iVCI implementations with respect to accuracy and central processing unit time.

FCIQMC-Tailored Distinguishable Cluster Approach: Open-Shell Systems

A recently proposed tailored approach based on the distinguishable cluster method and the stochastic FCI solver, FCIQMC [J. Chem. Theory Comput. 2020, 16, 5621], is extended to open-shell molecular systems. The method is employed to calculate spin gaps of various Fe(II) complexes, including a Fe(II) porphyrin model system. Both distinguishable cluster and fully relaxed CASSCF natural orbitals were used in this work as reference for the subsequent tailored distinguishable cluster calculations. The distinguishable cluster natural orbitals occupation numbers were also used as an aid to the selection of the active space. The effect of the active space sizes and of the explicit correlation correction (F12) onto the predicted spin gaps is investigated. The tailored distinguishable cluster with singles and doubles yields consistently more accurate results compared to the tailored coupled cluster with singles and doubles.

A shortcut to the thermodynamic limit for quantum many-body calculations of metals

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely dense k-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors.

Transcorrelated coupled cluster methods

Transcorrelated coupled cluster and distinguishable cluster methods are presented. The Hamiltonian is similarity transformed with a Jastrow factor in the first quantization, which results in up to three-body integrals. The coupled cluster with singles and doubles equations on this transformed Hamiltonian are formulated and implemented. It is demonstrated that the resulting methods have a superior basis set convergence and accuracy to the corresponding conventional and explicitly correlated methods. Additionally, approximations for three-body integrals are suggested and tested.

Dressing of Vertices by Cumulants in Multi-Reference Coupled Cluster

A new scheme is introduced in Multi-Reference (MR) Coupled Cluster (CC) based on the MR Generalized Normal Ordering (MRGNO) and the corresponding MR Generalized Wick Theorem (MRGWT) of Kutzelnigg and Mukherjee. The key element is the identification of a structure in MRGWT generated terms, facilitated by Goldstone diagram techniques. This allows for bundling the many terms of the MRGWT expansion and introduces a hierarchy in the equations that can be harnessed in devising approximations. One- and two-particle interaction vertices are found to be uniformly substituted for their counterpart dressed by density cumulants. This allows for a straightforward rewriting of the ordinary energy expression of CC with interaction dressed (id) one- and two-particle terms and reveals the presence of three- and higher-rank dressed interaction vertices too. Cumulants appearing out of dressed interaction vertices contribute to the amplitude equations and can be interpreted to have an amplitude dressing role. Dressing of one- and two-particle interaction vertices is the most straightforward to implement and does not hinder computational feasibility, provided that the reference function involves strictly limited active space sizes. The Generalized Valence Bond wave function, underlying pilot numerical tests, fulfills this criterion. Results on multiple bond breaking scenarios point to the need of stepping beyond one- and two-particle id. An extremely simple version of incorporating amplitude dressing in addition to one- and two-particle id is seen to cure the potential energy curves remarkably, stimulating further investigations along this line.

Variational coupled cluster for ground and excited states

In single-reference coupled-cluster (CC) methods, one has to solve a set of non-linear polynomial equations in order to determine the so-called amplitudes that are then used to compute the energy and other properties. Although it is of common practice to converge to the (lowest-energy) ground-state solution, it is also possible, thanks to tailored algorithms, to access higher-energy roots of these equations that may or may not correspond to genuine excited states. Here, we explore the structure of the energy landscape of variational CC and we compare it with its (projected) traditional version in the case where the excitation operator is restricted to paired double excitations (pCCD). By investigating two model systems (the symmetric stretching of the linear H₄ molecule and the continuous deformation of the square H₄ molecule into a rectangular arrangement) in the presence of weak and strong correlations, the performance of variational pCCD (VpCCD) and traditional pCCD is gauged against their configuration interaction (CI) equivalent, known as doubly occupied CI, for reference Slater determinants made of ground- or excited-state Hartree-Fock orbitals or state-specific orbitals optimized directly at the VpCCD level. The influence of spatial symmetry breaking is also investigated.

Accuracy of the distinguishable cluster approximation for triple excitations for open-shell molecules and excited states

The distinguishable cluster approximation for triple excitations has been applied to calculate thermochemical properties and excited states involving closed-shell and open-shell species, such as small molecules, 3d transition metal atoms, ozone, and an iron-porphyrin model. Excitation energies have been computed using the ΔCC approach by directly optimizing the excited states. A fixed-reference technique has been introduced to target selected spin-states for open-shell molecular systems. The distinguishable cluster approximation consistently improves coupled cluster with singles doubles and triples results for absolute and relative energies. For excited states dominated by a single configuration state function, the fixed-reference approach combined with high-level coupled-cluster methods has a comparable accuracy to the corresponding equation-of-motion coupled-cluster methods with a negligible amount of spin contamination.

High-Level Rovibrational Calculations on Ketenimine

From an astrochemical point of view ketenimine (CH₂CNH) is a complex organic molecule (COM) and therefore likely to be a building block for biologically relevant molecules. Since it has been detected in the star-forming region Sagittarius B2(N), it is of high relevance in this field. Although experimental data are available for certain bands, for some energy ranges such as above 1200 cm^-1 reliable data virtually do not exist. In addition, high-level ab initio calculations are neither reported for ketenimine nor for one of its deuterated isotopologues. In this paper, we provide for the first time data from accurate quantum chemical calculations and a thorough analysis of the full rovibrational spectrum. Based on high-level potential energy surfaces obtained from explicitly correlated coupled-cluster calculations including up to 4-mode coupling terms, the (ro)vibrational spectrum of ketenimine has been studied in detail by variational calculations relying on rovibrational configuration interaction (RVCI) theory. Strong Fermi resonances were found for all isotopologues. Rovibrational infrared intensities have been obtained from dipole moment surfaces determined from the distinguishable cluster approximation. A comparison of the spectra of the CH₂CNH molecule with experimental data validates our results, but also reveals new insight about the system, which shows very strong Coriolis coupling effects.

Tailored coupled cluster theory in varying correlation regimes

The tailored coupled cluster (TCC) approach is a promising ansatz that preserves the simplicity of single-reference coupled cluster theory while incorporating a multi-reference wave function through amplitudes obtained from a preceding multi-configurational calculation. Here, we present a detailed analysis of the TCC wave function based on model systems, which require an accurate description of both static and dynamic correlation. We investigate the reliability of the TCC approach with respect to the exact wave function. In addition to the error in the electronic energy and standard coupled cluster diagnostics, we exploit the overlap of TCC and full configuration interaction wave functions as a quality measure. We critically review issues, such as the required size of the active space, size-consistency, symmetry breaking in the wave function, and the dependence of TCC on the reference wave function. We observe that possible errors caused by symmetry breaking can be mitigated by employing the determinant with the largest weight in the active space as reference for the TCC calculation. We find the TCC model to be promising in calculations with active orbital spaces which include all orbitals with a large single-orbital entropy, even if the active spaces become very large and then may require modern active-space approaches that are not restricted to comparatively small numbers of orbitals. Furthermore, utilizing large active spaces can improve on the TCC wave function approximation and reduce the size-consistency error because the presence of highly excited determinants affects the accuracy of the coefficients of low-excited determinants in the active space.

Doubly Occupied Pair Coupled Cluster F12 Approach

Inspired by reports of the good performance of the doubly occupied pair coupled cluster (pCCD) theory in describing static electron correlation, we have introduced and implemented a variant thereof that includes single excitations and explicitly treats the dynamic electron correlation using the F12 methodology (pCCSD-F12). This drastically reduces the computation scaling with respect to the standard method using the full double-excitation operator (CCSD-F12). Slater-type geminals as a correlation factor, together with fixed cusp conditions, were used, which is known as the SP-ansatz. For sample model systems, we have investigated the performance of reference states constructed from either canonical or localized molecular orbitals. Finaly, the employment of Brueckner orbitals has been tested, which causes the single excitations to naturally vanish from the wave function expansion (B-pCCD-F12). Our test systems include different-sized rings of hydrogen atoms and dissociation curves for small molecules such as HF, N₂, and CO₂; and comparison with CCSD-F12 is presented for a series of reaction enthalpies.

On the synergy of matrix-isolation infrared spectroscopy and vibrational configuration interaction computations

The key feature of matrix-isolation infrared (MI-IR) spectroscopy is the isolation of single guest molecules in a host system at cryogenic conditions. The matrix mostly hinders rotation of the guest molecule, providing access to pure vibrational features. Vibrational self-consistent field (VSCF) and configuration interaction computations (VCI) on ab initio multimode potential energy surfaces (PES) give rise to anharmonic vibrational spectra. In a single-sourced combination of these experimental and computational approaches, we have established an iterative spectroscopic characterization procedure. The present article reviews the scope of this procedure by highlighting the strengths and limitations based on the examples of water, carbon dioxide, methane, methanol, and fluoroethane. An assessment of setups for the construction of the multimode PES on the example of methanol demonstrates that CCSD(T)-F12 level of theory is preferable to compute (a) accurate vibrational frequencies and (b) equilibrium or vibrationally averaged structural parameters. Our procedure has allowed us to uniquely assign unknown or disputed bands and enabled us to clarify problematic spectral regions that are crowded with combination bands and overtones. Besides spectroscopic assignment, the excellent agreement between theory and experiment paves the way to tackle questions of rather fundamental nature as to whether or not matrix effects are systematic, and it shows the limits of conventional notations used by spectroscopists.

FCIQMC-Tailored Distinguishable Cluster Approach

The tailored approach is applied to the distinguishable cluster method together with a stochastic FCI solver (FCIQMC). It is demonstrated that the new method is more accurate than the corresponding tailored coupled cluster and the pure distinguishable cluster methods. An F12 correction for tailored methods and FCIQMC is introduced, which drastically improves the basis set convergence. A new black-box approach to define the active space using the natural orbitals from the distinguishable cluster is evaluated and found to be a convenient alternative to the usual CASSCF approach.

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Duality of Ring and Ladder Diagrams and Its Importance for Many-Electron Perturbation Theories

We present a diagrammatic decomposition of the transition pair correlation function for the uniform electron gas. We demonstrate explicitly that ring and ladder diagrams are dual counterparts that capture significant long- and short-ranged interelectronic correlation effects, respectively. Our findings help to guide the further development of approximate many-electron theories and reveal that the contribution of the ladder diagrams to the electronic correlation energy can be approximated in an effective manner using second-order perturbation theory. We employ the latter approximation to reduce the computational cost of coupled cluster theory calculations for insulators and semiconductors by 2 orders of magnitude without compromising accuracy.

Improving the Efficiency of the Multireference Driven Similarity Renormalization Group via Sequential Transformation, Density Fitting, and the Noninteracting Virtual Orbital Approximation

This study examines several techniques to improve the efficiency of the linearized multireference driven similarity renormalization group truncated to one- and two-body operators [MR-LDSRG(2)]. We propose a sequential MR-LDSRG(2) [sq-MR-LDSRG(2)] scheme, in which one-body substitutions are folded exactly into the Hamiltonian. This new approach is combined with density fitting (DF) to reduce the storage cost of two-electron integrals. To further avoid storage of large four-index intermediates, we propose a noninteracting virtual orbital (NIVO) approximation of the Baker-Campbell-Hausdorff series that neglects commutators terms with three and four virtual indices. The NIVO approximation reduces the computational prefactor of the MR-LDSRG(2), bringing it closer to that of coupled cluster with singles and doubles (CCSD). We test the effect of the DF and NIVO approximations on the MR-LDSRG(2) and sq-MR-LDSRG(2) methods by computing properties of eight diatomic molecules. The diatomic constants obtained by DF-sq-MR-LDSRG(2)+NIVO are found to be as accurate as those from the original MR-LDSRG(2) and coupled cluster theory with singles, doubles, and perturbative triples. Finally, we demonstrate that the DF-sq-MR-LDSRG(2)+NIVO scheme can be applied to chemical systems with more than 550 basis functions by computing the automerization energy of cyclobutadiene with a quintuple-ζ basis set. The predicted automerization energy is found to be similar to the value computed with Mukherjee's state-specific multireference coupled cluster theory with singles and doubles.

Role of Valence and Semicore Electron Correlation on Spin Gaps in Fe(II)-Porphyrins

The role of valence and semicore correlation in differentially stabilizing the intermediate spin state of Fe(II)-porphyrins is analyzed. For CASSCF treatments of valence correlation, a (32,34) active space containing metal 3d, d′ orbitals and the entire π system of the porphyrin is necessary to stabilize the intermediate spin state. Semicore correlation provides a minor (−1.6 kcal/mol) but quantitatively significant correction. Accounting for valence, semicore, and correlation beyond the active space enlarges the (³E_g–⁵A_1g) spin gap to −5.7 kcal/mol.

Systematic High-Accuracy Prediction of Electron Affinities for Biological Quinones

Quinones play vital roles as electron carriers in fundamental biological processes; therefore, the ability to accurately predict their electron affinities is crucial for understanding their properties and function. The increasing availability of cost-effective implementations of correlated wave function methods for both closed-shell and open-shell systems offers an alternative to density functional theory approaches that have traditionally dominated the field despite their shortcomings. Here, we define a benchmark set of quinones with experimentally available electron affinities and evaluate a range of electronic structure methods, setting a target accuracy of 0.1 eV. Among wave function methods, we test various implementations of coupled cluster (CC) theory, including local pair natural orbital (LPNO) approaches to canonical and parameterized CCSD, the domain-based DLPNO approximation, and the equations-of-motion approach for electron affinities, EA-EOM-CCSD. In addition, several variants of canonical, spin-component-scaled, orbital-optimized, and explicitly correlated (F12) Møller-Plesset perturbation theory are benchmarked. Achieving systematically the target level of accuracy is challenging and a composite scheme that combines canonical CCSD(T) with large basis set LPNO-based extrapolation of correlation energy proves to be the most accurate approach. Methods that offer comparable performance are the parameterized LPNO-pCCSD, the DLPNO-CCSD(T₀ ), and the orbital optimized OO-SCS-MP2. Among DFT methods, viable practical alternatives are only the M06 and the double hybrids, but the latter should be employed with caution because of significant basis set sensitivity. A highly accurate yet cost-effective DLPNO-based coupled cluster approach is used to investigate the methoxy conformation effect on the electron affinities of ubiquinones found in photosynthetic bacterial reaction centers. © 2018 Wiley Periodicals, Inc.