Showing papers on "Coupled cluster published in 2016"

The theory of variational hybrid quantum-classical algorithms

Abstract: Many quantum algorithms have daunting resource requirements when compared to what is available today. To address this discrepancy, a quantum-classical hybrid optimization scheme known as 'the quantum variational eigensolver' was developed (Peruzzo et al 2014 Nat. Commun. 5 4213) with the philosophy that even minimal quantum resources could be made useful when used in conjunction with classical routines. In this work we extend the general theory of this algorithm and suggest algorithmic improvements for practical implementations. Specifically, we develop a variational adiabatic ansatz and explore unitary coupled cluster where we establish a connection from second order unitary coupled cluster to universal gate sets through a relaxation of exponential operator splitting. We introduce the concept of quantum variational error suppression that allows some errors to be suppressed naturally in this algorithm on a pre-threshold quantum device. Additionally, we analyze truncation and correlated sampling in Hamiltonian averaging as ways to reduce the cost of this procedure. Finally, we show how the use of modern derivative free optimization techniques can offer dramatic computational savings of up to three orders of magnitude over previously used optimization techniques.

Scalable Quantum Simulation of Molecular Energies

Abstract: We report the first electronic structure calculation performed on a quantum computer without exponentially costly precompilation. We use a programmable array of superconducting qubits to compute the energy surface of molecular hydrogen using two distinct quantum algorithms. First, we experimentally execute the unitary coupled cluster method using the variational quantum eigensolver. Our efficient implementation predicts the correct dissociation energy to within chemical accuracy of the numerically exact result. Second, we experimentally demonstrate the canonical quantum algorithm for chemistry, which consists of Trotterization and quantum phase estimation. We compare the experimental performance of these approaches to show clear evidence that the variational quantum eigensolver is robust to certain errors. This error tolerance inspires hope that variational quantum simulations of classically intractable molecules may be viable in the near future.

Sparse maps--A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory.

Abstract: Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

Decomposition of Intermolecular Interaction Energies within the Local Pair Natural Orbital Coupled Cluster Framework

Abstract: The local coupled cluster method DLPNO-CCSD(T) allows calculations on systems containing hundreds of atoms to be performed while typically reproducing canonical CCSD(T) energies with chemical accuracy. In this work, we present a scheme for decomposing the DLPNO-CCSD(T) interaction energy between two molecules into physical meaningful contributions, providing a quantification of the most important components of the chemical interaction. The method, called Local Energy Decomposition (LED), is straightforward and requires negligible additional computing time. Both the Hartree-Fock and the correlation energy are decomposed into contributions from localized or pairs of localized occupied orbitals. Assigning these localized orbitals to fragments allows one to differentiate between intra- and intermolecular contributions to the interaction energy. Accordingly, the interaction energy can be decomposed into electronic promotion, electrostatic, exchange, dynamic charge polarization, and dispersion contributions. The LED scheme is applied to a number of test cases ranging from weakly, dispersively bound complexes to systems with strong ionic interactions. The dependence of the results on the one-particle basis set and various technical aspects, such as the localization scheme, are carefully studied in order to ensure that the results do not suffer from technical artifacts. A numerical comparison between the DLPNO-CCSD(T)/LED and the popular symmetry adapted perturbation theory (DFT-SAPT) is made, and the limitations of the proposed scheme are discussed.

On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice

Abstract: The MB-pol many-body potential has recently emerged as an accurate molecular model for water simulations from the gas to the condensed phase. In this study, the accuracy of MB-pol is systematically assessed across the three phases of water through extensive comparisons with experimental data and high-level ab initio calculations. Individual many-body contributions to the interaction energies as well as vibrational spectra of water clusters calculated with MB-pol are in excellent agreement with reference data obtained at the coupled cluster level. Several structural, thermodynamic, and dynamical properties of the liquid phase at atmospheric pressure are investigated through classical molecular dynamics simulations as a function of temperature. The structural properties of the liquid phase are in nearly quantitative agreement with X-ray diffraction data available over the temperature range from 268 to 368 K. The analysis of other thermodynamic and dynamical quantities emphasizes the importance of explicitly including nuclear quantum effects in the simulations, especially at low temperature, for a physically correct description of the properties of liquid water. Furthermore, both densities and lattice energies of several ice phases are also correctly reproduced by MB-pol. Following a recent study of DFT models for water, a score is assigned to each computed property, which demonstrates the high and, in many respects, unprecedented accuracy of MB-pol in representing all three phases of water.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods.

Abstract: The performance of different GW methods is assessed for a set of 24 organic acceptors. Errors are evaluated with respect to coupled cluster singles, doubles, and perturbative triples [CCSD(T)] reference data for the vertical ionization potentials (IPs) and electron affinities (EAs), extrapolated to the complete basis set limit. Additional comparisons are made to experimental data, where available. We consider fully self-consistent GW (scGW), partial self-consistency in the Green's function (scGW0), non-self-consistent G0W0 based on several mean-field starting points, and a "beyond GW" second-order screened exchange (SOSEX) correction to G0W0. We also describe the implementation of the self-consistent Coulomb hole with screened exchange method (COHSEX), which serves as one of the mean-field starting points. The best performers overall are G0W0+SOSEX and G0W0 based on an IP-tuned long-range corrected hybrid functional with the former being more accurate for EAs and the latter for IPs. Both provide a balanced treatment of localized vs delocalized states and valence spectra in good agreement with photoemission spectroscopy (PES) experiments.

Towards a pair natural orbital coupled cluster method for excited states

Abstract: The use of back-transformed pair natural orbitals in the calculation of excited state energies, ionization potentials, and electron affinities is investigated within the framework of equation of motion coupled cluster theory and its similarity transformed variant. Possible approaches to a more optimal use of pair natural orbitals in these methods are indicated.

Benchmark ab Initio Conformational Energies for the Proteinogenic Amino Acids through Explicitly Correlated Methods. Assessment of Density Functional Methods

Abstract: The relative energies of the YMPJ conformer database of the 20 proteinogenic amino acids, with N- and C-termination, have been re-evaluated using explicitly correlated coupled cluster methods. Lower-cost ab initio methods such as MP2-F12 and CCSD-F12b actually are outperformed by double-hybrid DFT functionals; in particular, the DSD-PBEP86-NL double hybrid performs well enough to serve as a secondary standard. Among range-separated hybrids, ωB97X-V performs well, while B3LYP-D3BJ does surprisingly well among traditional DFT functionals. Treatment of dispersion is important for the DFT functionals; for the YMPJ set, D3BJ generally works as well as the NL nonlocal dispersion functional. Basis set sensitivity for DFT calculations on these conformers is weak enough that def2-TZVP is generally adequate. For conformer corrections to heats of formation, B3LYP-D3BJ and especially DSD-PBEP86-D3BJ or DSD-PBEP86-NL are adequate for all but the most exacting applications. The revised geometries and energetics for the...

Coupled Cluster Method with Single and Double Excitations Tailored by Matrix Product State Wave Functions

Abstract: In the past decade, the quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as the method of choice for calculations of strongly correlated molecular systems. Despite its favorable scaling, it is in practice not suitable for computations of dynamic correlation. We present a novel method for accurate “post-DMRG” treatment of dynamic correlation based on the tailored coupled cluster (CC) theory in which the DMRG method is responsible for the proper description of nondynamic correlation, whereas dynamic correlation is incorporated through the framework of the CC theory. We illustrate the potential of this method on prominent multireference systems, in particular, N2 and Cr2 molecules and also oxo-Mn(Salen), for which we have performed the first post-DMRG computations in order to shed light on the energy ordering of the lowest spin states.

Combined 4-component and relativistic pseudopotential study of ThO for the electron electric dipole moment search

Abstract: A precise theoretical study of the electronic structure of heavy atom diatomic molecules is of key importance to interpret the experiments in the search for violation of time-reversal (T) and spatial-parity (P) symmetries of fundamental interactions in terms of the electron electric dipole moment, eEDM, and dimensionless constant, kT,P, characterizing the strength of the T,P-odd pseudoscalar-scalar electron-nucleus neutral current interaction. The ACME collaboration has recently improved limits on these quantities using a beam of ThO molecules in the electronic H3Δ1 state [J. Baron et al., Science 343, 269 (2014)]. We apply the combined direct relativistic 4-component and two-step relativistic pseudopotential/restoration approaches to a benchmark calculation of the effective electric field, Eeff, parameter of the T,P-odd pseudoscalar-scalar interaction, WT,P, and hyperfine structure constant in Δ13 state of the ThO molecule. The first two parameters are required to interpret the experimental data in terms of the eEDM and kT,P constant. We have investigated the electron correlation for all of the 98 electrons of ThO simultaneously up to the level of the coupled cluster with single, double, and noniterative triple amplitudes, CCSD(T), theory. Contributions from iterative triple and noniterative quadruple cluster amplitudes for the valence electrons have been also treated. The obtained values are Eeff = 79.9 GV/cm, WT,P = 113.1 kHz. The theoretical uncertainty of these values is estimated to be about two times smaller than that of our previous study [L. V. Skripnikov and A. V. Titov, J. Chem. Phys., 142, 024301 (2015)]. It was found that the correlation of the inner- and outer-core electrons contributes 9% to the effective electric field. The values of the molecule frame dipole moment of the Δ13 state and the H3Δ1→X1Σ+ transition energy of ThO calculated within the same methods are in a very good agreement with the experiment.

Correlation consistent basis sets for lanthanides: The atoms La–Lu

Abstract: Using the 3rd-order Douglas-Kroll-Hess (DKH3) Hamiltonian, all-electron correlation consistent basis sets of double-, triple-, and quadruple-zeta quality have been developed for the lanthanide elements La through Lu. Basis sets designed for the recovery of valence correlation (defined here as 4f5s5p5d6s), cc-pVnZ-DK3, and outer-core correlation (valence + 4s4p4d), cc-pwCVnZ-DK3, are reported (n = D, T, and Q). Systematic convergence of both Hartree-Fock and correlation energies towards their respective complete basis set (CBS) limits are observed. Benchmark calculations of the first three ionization potentials (IPs) of La through Lu are reported at the DKH3 coupled cluster singles and doubles with perturbative triples, CCSD(T), level of theory, including effects of correlation down through the 4s electrons. Spin-orbit coupling is treated at the 2-component HF level. After extrapolation to the CBS limit, the average errors with respect to experiment were just 0.52, 1.14, and 4.24 kcal/mol for the 1st, 2nd, and 3rd IPs, respectively, compared to the average experimental uncertainties of 0.03, 1.78, and 2.65 kcal/mol, respectively. The new basis sets are also used in CCSD(T) benchmark calculations of the equilibrium geometries, atomization energies, and heats of formation for Gd2, GdF, and GdF3. Except for the equilibrium geometry and harmonic frequency of GdF, which are accurately known from experiment, all other calculated quantities represent significant improvements compared to the existing experimental quantities. With estimated uncertainties of about ±3 kcal/mol, the 0 K atomization energies (298 K heats of formation) are calculated to be (all in kcal/mol): 33.2 (160.1) for Gd2, 151.7 (-36.6) for GdF, and 447.1 (-295.2) for GdF3.

Spectral functions of the uniform electron gas via coupled-cluster theory and comparison to the G W and related approximations

Abstract: We use ab initio coupled-cluster theory to compute the spectral function of the uniform electron gas at a Wigner-Seitz radius of r_s=4. The coupled-cluster approximations we employ go significantly beyond the diagrammatic content of state-of-the-art GW theory. We compare our calculations extensively to GW and GW-plus-cumulant theory, illustrating the strengths and weaknesses of these methods in capturing the quasiparticle and satellite features of the electron gas. Our accurate calculations further allow us to address the long-standing debate over the occupied bandwidth of metallic sodium. Our findings indicate that the future application of coupled-cluster theory to condensed phase material spectra is highly promising.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules II: Non-Empirically Tuned Long-Range Corrected Hybrid Functionals.

Abstract: The performance of non-empirically tuned long-range corrected hybrid functionals for the prediction of vertical ionization potentials (IPs) and electron affinities (EAs) is assessed for a set of 24 organic acceptor molecules. Basis set-extrapolated coupled cluster singles, doubles, and perturbative triples [CCSD(T)] calculations serve as a reference for this study. Compared to standard exchange-correlation functionals, tuned long-range corrected hybrid functionals produce highly reliable results for vertical IPs and EAs, yielding mean absolute errors on par with computationally more demanding GW calculations. In particular, it is demonstrated that long-range corrected hybrid functionals serve as ideal starting points for non-self-consistent GW calculations.

Ionization Energies and Aqueous Redox Potentials of Organic Molecules: Comparison of DFT, Correlated ab Initio Theory and Pair Natural Orbital Approaches.

Abstract: The calculation of redox potentials involves large energetic terms arising from gas phase ionization energies, thermodynamic contributions, and solvation energies of the reduced and oxidized species. In this work we study the performance of a wide range of wave function and density functional theory methods for the prediction of ionization energies and aqueous one-electron oxidation potentials of a set of 19 organic molecules. Emphasis is placed on evaluating methods that employ the computationally efficient local pair natural orbital (LPNO) approach, as well as several implementations of coupled cluster theory and explicitly correlated F12 methods. The electronic energies are combined with implicit solvation models for the solvation energies. With the exception of MP2 and its variants, which suffer from enormous errors arising at least partially from the poor Hartree–Fock reference, ionization energies can be systematically predicted with average errors below 0.1 eV for most of the correlated wave functi...

Toward a W4-F12 approach: Can explicitly correlated and orbital-based ab initio CCSD(T) limits be reconciled?

Abstract: In the context of high-accuracy computational thermochemistry, the valence coupled cluster with all singles and doubles (CCSD) correlation component of molecular atomization energies presents the most severe basis set convergence problem, followed by the (T) component. In the present paper, we make a detailed comparison, for an expanded version of the W4-11 thermochemistry benchmark, between, on the one hand, orbital-based CCSD/AV{5,6}Z + d and CCSD/ACV{5,6}Z extrapolation, and on the other hand CCSD-F12b calculations with cc-pVQZ-F12 and cc-pV5Z-F12 basis sets. This latter basis set, now available for H–He, B–Ne, and Al–Ar, is shown to be very close to the basis set limit. Apparent differences (which can reach 0.35 kcal/mol for systems like CCl4) between orbital-based and CCSD-F12b basis set limits disappear if basis sets with additional radial flexibility, such as ACV{5,6}Z, are used for the orbital calculation. Counterpoise calculations reveal that, while total atomization energies with V5Z-F12 basis sets are nearly free of BSSE, orbital calculations have significant BSSE even with AV(6 + d)Z basis sets, leading to non-negligible differences between raw and counterpoise-corrected extrapolated limits. This latter problem is greatly reduced by switching to ACV{5,6}Z core-valence basis sets, or simply adding an additional zeta to just the valence orbitals. Previous reports that all-electron approaches like HEAT (high-accuracy extrapolated ab-initio thermochemistry) lead to different CCSD(T) limits than “valence limit + CV correction” approaches like Feller-Peterson-Dixon and Weizmann-4 (W4) theory can be rationalized in terms of the greater radial flexibility of core-valence basis sets. For (T) corrections, conventional CCSD(T)/AV{Q,5}Z + d calculations are found to be superior to scaled or extrapolated CCSD(T)-F12b calculations of similar cost. For a W4-F12 protocol, we recommend obtaining the Hartree-Fock and valence CCSD components from CCSD-F12b/cc-pV{Q,5}Z-F12 calculations, but the (T) component from conventional CCSD(T)/aug’-cc-pV{Q,5}Z + d calculations using Schwenke’s extrapolation; post-CCSD(T), core-valence, and relativistic corrections are to be obtained as in the original W4 theory. W4-F12 is found to agree slightly better than W4 with ATcT (active thermochemical tables) data, at a substantial saving in computation time and especially I/O overhead. A W4-F12 calculation on benzene is presented as a proof of concept.

The Delicate Balance of Static and Dynamic Electron Correlation.

Abstract: Multi-configurational approaches yield universal wave function parametrizations that can qualitatively well describe electronic structures along reaction pathways. For quantitative results, multi-reference perturbation theory is required to capture dynamic electron correlation from the otherwise neglected virtual orbitals. Still, the overall accuracy suffers from the finite size and choice of the active orbital space and peculiarities of the perturbation theory. Fortunately, the electronic wave functions at equilibrium structures of reactants and products can often be well described by single-reference methods and hence are accessible to accurate coupled cluster calculations. Here, we calculate the heterolytic double dissociation energy of four 3d-metallocenes with the complete active space self-consistent field method and compare to highly accurate coupled cluster data. Our coupled cluster data are well within the experimental error bars. This accuracy can also be approached by complete active space calc...

How important is thermal expansion for predicting molecular crystal structures and thermochemistry at finite temperatures

Abstract: Molecular crystals expand appreciably upon heating due to both zero-point and thermal vibrational motion, yet this expansion is often neglected in molecular crystal modeling studies. Here, a quasi-harmonic approximation is coupled with fragment-based hybrid many-body interaction calculations to predict thermal expansion and finite-temperature thermochemical properties in crystalline carbon dioxide, ice Ih, acetic acid and imidazole. Fragment-based second-order Moller-Plesset perturbation theory (MP2) and coupled cluster theory with singles, doubles and perturbative triples [CCSD(T)] predict the thermal expansion and the temperature dependence of the enthalpies, entropies and Gibbs free energies of sublimation in good agreement with experiment. The errors introduced by neglecting thermal expansion in the enthalpy and entropy cancel somewhat in the Gibbs free energy. The resulting ∼ 1-2 kJ mol(-1) errors in the free energy near room temperature are comparable to or smaller than the errors expected from the electronic structure treatment, but they may be sufficiently large to affect free-energy rankings among energetically close polymorphs.

Linear Absorption Spectra from Explicitly Time-Dependent Equation-of-Motion Coupled-Cluster Theory

Abstract: We report an explicitly time-dependent approach to the generation of linear absorption spectra for molecular systems within the framework of equation-of-motion (EOM) coupled-cluster (CC) theory. While most time-dependent CC approaches consider the perturbation and time-evolution of a CC wave function, the present work considers the time-evolution of a CC dipole function. The dipole function formalism introduces no approximations and requires the evolution of only one time-dependent quantity, either the left or right dipole function. This time-dependent framework can be used to compute linear absorption spectra for molecules with a high density of states over a broad spectral range, a case for which conventional frequency-domain computations may become impractical. We validate the approach by comparing absorption spectra for small molecules computed at EOM second-order approximate CC (CC2) and time-dependent EOM-CC2 (TD-EOM-CC2) levels of theory. TD-EOM-CC2 computations are also used to predict extreme ult...

Coupled Cluster Evaluation of the Stability of Atmospheric Acid-Base Clusters with up to 10 Molecules.

Abstract: We investigate the utilization of the domain local pair natural orbital coupled cluster (DLPNO-CCSD(T)) method for calculating binding energies of atmospherical molecular clusters. Applied to small complexes of atmospherical relevance we find that the DLPNO method significantly reduces the scatter in the binding energy, which is commonly present in DFT calculations. For medium sized clusters consisting of sulfuric acid and bases the DLPNO method yields a systematic underestimation of the binding energy compared to canonical coupled cluster results. The errors in the DFT binding energies appear to be more random, while the systematic nature of the DLPNO results allows the establishment of a scaling factor, to better mimic the canonical coupled cluster calculations. Based on the trends identified for the small and medium sized systems, we further extend the application of the DLPNO method to large acid - base clusters consisting of up to 10 molecules, which have previously been out of reach with accurate co...

A new approach to approximate equation-of-motion coupled cluster with triple excitations.

Abstract: Accurate methods for excited, ionized, and electron-attached states are critical to the study of many chemical species such as reactive intermediates, radicals, and ionized systems. The equation-of-motion coupled cluster singles, doubles, and triples (EOM-CCSDT) family of methods is very accurate (roughly similar in accuracy as for CCSDT calculations of the ground state), but the computational cost scales iteratively as the eighth power of the system size. Many approximations already exist, although most either correct only the excited state or require an iterative 𝒪(n7) procedure which can also be prohibitively expensive. In this paper, new methods, termed EOM-CCSD(T)(a) and EOM-CCSD(T)(a)*, are proposed which correct both the ground and excited states based on a shared effective Hamiltonian, and the latter of which includes only non-iterative corrections to both the CCSD and EOM-CCSD energies. These methods are found to significantly improve the description of excited and ionized potential energy surfac...

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules I. Reference Data at the CCSD(T) Complete Basis Set Limit.

Abstract: In designing organic materials for electronics applications, particularly for organic photovoltaics (OPV), the ionization potential (IP) of the donor and the electron affinity (EA) of the acceptor play key roles. This makes OPV design an appealing application for computational chemistry since IPs and EAs are readily calculable from most electronic structure methods. Unfortunately reliable, high-accuracy wave function methods, such as coupled cluster theory with single, double, and perturbative triples [CCSD(T)] in the complete basis set (CBS) limit are too expensive for routine applications to this problem for any but the smallest of systems. One solution is to calibrate approximate, less computationally expensive methods against a database of high-accuracy IP/EA values; however, to our knowledge, no such database exists for systems related to OPV design. The present work is the first of a multipart study whose overarching goal is to determine which computational methods can be used to reliably compute IP...

Polynomial Similarity Transformation Theory: A smooth interpolation between coupled cluster doubles and projected BCS applied to the reduced BCS Hamiltonian

Abstract: We present a similarity transformation theory based on a polynomial form of a particle-hole pair excitation operator. In the weakly correlated limit, this polynomial becomes an exponential, leading to coupled cluster doubles. In the opposite strongly correlated limit, the polynomial becomes an extended Bessel expansion and yields the projected BCS wave function. In between, we interpolate using a single parameter. The effective Hamiltonian is non-Hermitian and this polynomial similarity transformation theory follows the philosophy of traditional coupled cluster, left projecting the transformed Hamiltonian onto subspaces of the Hilbert space in which the wave function variance is forced to be zero. Similarly, the interpolation parameter is obtained through minimizing the next residual in the projective hierarchy. We rationalize and demonstrate how and why coupled cluster doubles is ill suited to the strongly correlated limit, whereas the Bessel expansion remains well behaved. The model provides accurate wave functions with energy errors that in its best variant are smaller than 1% across all interaction strengths. The numerical cost is polynomial in system size and the theory can be straightforwardly applied to any realistic Hamiltonian.

Speeding up equation of motion coupled cluster theory with the chain of spheres approximation

Abstract: In the present paper, the chain of spheres exchange (COSX) approximation is applied to the highest scaling terms in the equation of motion (EOM) coupled cluster equations with single and double excitations, in particular, the terms involving integrals with four virtual labels. It is found that even the acceleration of this single term yields significant computational gains without compromising the desired accuracy of the method. For an excitation energy calculation on a cluster of five water molecules using 585 basis functions, the four virtual term is 9.4 times faster using COSX with a loose grid than using the canonical implementation, which yields a 2.6 fold acceleration for the whole of the EOM calculation. For electron attachment calculations, the four virtual term is 15 times and the total EOM calculation is 10 times faster than the canonical calculation for the same system. The accuracy of the new method was tested using Thiel’s test set for excited states using the same settings and the maximum absolute deviation over the whole test set was found to be 12.945 cm−1 (59 μHartree) for excitation energies and 6.799 cm−1 (31 μHartree) for electron attachments. Using MP2 amplitudes for the ground state in combination with the parallel evaluation of the full EOM equations in the manner discussed in this paper enabled us to perform calculations for large systems. Electron affinity values for the two lowest states of a Zn protoporphyrine model compound (224 correlated electrons and 1120 basis functions) were obtained in 3 days 19 h using 4 cores of a Xeon E5-2670 processor allocating 10 GB memory per core. Calculating the lowest two excitation energies for trans-retinal (114 correlated electrons and 539 basis functions) took 1 day 21 h using eight cores of the same processor and identical memory allocation per core.

Communication: Finite size correction in periodic coupled cluster theory calculations of solids.

Abstract: We present a method to correct for finite size errors in coupled cluster theory calculations of solids. The outlined technique shares similarities with electronic structure factor interpolation methods used in quantum Monte Carlo calculations. However, our approach does not require the calculation of density matrices. Furthermore we show that the proposed finite size corrections achieve chemical accuracy in the convergence of second-order Moller-Plesset perturbation and coupled cluster singles and doubles correlation energies per atom for insulating solids with two atomic unit cells using 2 × 2 × 2 and 3 × 3 × 3 k-point meshes only.

Predicting finite-temperature properties of crystalline carbon dioxide from first principles with quantitative accuracy

Abstract: Molecular crystal structures, thermodynamics, and mechanical properties can vary substantially with temperature, and predicting these temperature-dependencies correctly is important for many practical applications in the pharmaceutical industry and other fields. However, most electronic structure predictions of molecular crystal properties neglect temperature and/or thermal expansion, leading to potentially erroneous results. Here, we demonstrate that by combining large basis set second-order Moller–Plesset (MP2) or even coupled cluster singles, doubles, and perturbative triples (CCSD(T)) electronic structure calculations with a quasiharmonic treatment of thermal expansion, experimentally observable properties such as the unit cell volume, heat capacity, enthalpy, entropy, sublimation point and bulk modulus of phase I crystalline carbon dioxide can be predicted in excellent agreement with experiment over a broad range of temperatures. These results point toward a promising future for ab initio prediction of molecular crystal properties at real-world temperatures and pressures.

A Projector-Embedding Approach for Multiscale Coupled-Cluster Calculations Applied to Citrate Synthase

Abstract: Projector-based embedding has recently emerged as a robust multiscale method for the calculation of various electronic molecular properties. We present the coupling of projector embedding with quantum mechanics/molecular mechanics modeling and apply it for the first time to an enzyme-catalyzed reaction. Using projector-based embedding, we combine coupled-cluster theory, density-functional theory (DFT), and molecular mechanics to compute energies for the proton abstraction from acetyl-coenzyme A by citrate synthase. By embedding correlated ab initio methods in DFT we eliminate functional sensitivity and obtain high-accuracy profiles in a procedure that is straightforward to apply.

Analytic energy derivatives for the calculation of the first-order molecular properties using the domain-based local pair-natural orbital coupled-cluster theory

Abstract: The domain-based local pair-natural orbital coupled-cluster (DLPNO-CC) theory has recently emerged as an efficient and powerful quantum-chemical method for the calculation of energies of molecules comprised of several hundred atoms. It has been demonstrated that the DLPNO-CC approach attains the accuracy of a standard canonical coupled-cluster calculation to about 99.9% of the basis set correlation energy while realizing linear scaling of the computational cost with respect to system size. This is achieved by combining (a) localized occupied orbitals, (b) large virtual orbital correlation domains spanned by the projected atomic orbitals (PAOs), and (c) compaction of the virtual space through a truncated pair natural orbital (PNO) basis. In this paper, we report on the implementation of an analytic scheme for the calculation of the first derivatives of the DLPNO-CC energy for basis set independent perturbations within the singles and doubles approximation (DLPNO-CCSD) for closed-shell molecules. Perturbati...

π-π stacking between polyaromatic hydrocarbon sheets beyond dispersion interactions.

Abstract: High level ab initio calculations ranging from coupled cluster methods including explicitly correlated approaches to standard second order Moller–Plesset theory using spin scaling (SOS-MP2) have been performed on sandwich and slipped parallel dimer structures of a series of quasi one-dimensional acenes and on two-dimensional sheets containing the pyrene to coronene series encircled with two layers of benzene rings. Sandwich (graphitic AA type) and slipped parallel (AB type) structures were considered and, within the given symmetry restrictions, full geometry optimizations were performed. Basis set superposition effects have been considered. The computed geometries show a significant biconcave deviation of the two-dimensional sheets from planarity with the central intersheet C⋯C distances considerably smaller that van der Waals distances. The computed intersheet binding energy per carbon atom extrapolated for N → ∞ of −74.3 meV (1.713 kcal mol−1) per atom agrees quite well with an experimental defoliation energy of −52 meV (1.199 kcal mol−1) per atom (−67 meV (1.545 kcal mol−1) per carbon atom without corrections for H binding contributions) for polyaromatic hydrocarbons (PAHs) from graphite. A limited investigation of density functional theory (DFT) calculations using empirical dispersion contributions has been performed also showing a significant underbinding character of the D3 method. For most of the DFT variants investigated the graphene sheet models retain a quasi-planar structure in strong contrast to the aforementioned SOS-MP2 results.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules IV: Electron-Propagator Methods.

Abstract: Comparison of ab initio electron-propagator predictions of vertical ionization potentials and electron affinities of organic, acceptor molecules with benchmark calculations based on the basis set-extrapolated, coupled cluster single, double, and perturbative triple substitution method has enabled identification of self-energy approximations with mean, unsigned errors between 0.1 and 0.2 eV. Among the self-energy approximations that neglect off-diagonal elements in the canonical, Hartree–Fock orbital basis, the P3 method for electron affinities, and the P3+ method for ionization potentials provide the best combination of accuracy and computational efficiency. For approximations that consider the full self-energy matrix, the NR2 methods offer the best performance. The P3+ and NR2 methods successfully identify the correct symmetry label of the lowest cationic state in two cases, naphthalenedione and benzoquinone, where some other methods fail.