Assessment of Orbital-Optimized MP2.5 for Thermochemistry and Kinetics: Dramatic Failures of Standard Perturbation Theory Approaches for Aromatic Bond Dissociation Energies and Barrier Heights of Radical Reactions

Equation-of-motion orbital-optimized coupled-cluster doubles method with the density-fitting approximation: An efficient implementation

Orbital-optimized coupled-cluster methods are very helpful for theoretical predictions of the molecular properties of challenging chemical systems, such as excited states. In this research, an efficient implementation of the equation-of-motion orbital-optimized coupled-cluster doubles method with the density-fitting (DF) approach, denoted by DF-EOM-OCCD, is presented. The computational cost of the DF-EOM-OCCD method for excitation energies is compared with that of the conventional EOM-OCCD method. Our results demonstrate that DF-EOM-OCCD excitation energies are dramatically accelerated compared to EOM-OCCD. There are almost 17-fold reductions for theC 5 H 12 $$ {\mathrm{C}}_5{\mathrm{H}}_{12} $$ molecule in an aug-cc-pVTZ basis set with the RHF reference. This dramatic performance improvement comes from the reduced cost of integral transformation with the DF approach and the efficient evaluation of the particle-particle ladder (PPL) term, which is the most expensive term to evaluate. Further, our results show that the DF-EOM-OCCD approach is very helpful for the computation of excitation energies in open-shell molecular systems. Overall, we conclude that our new DF-EOM-OCCD implementation is very promising for the study of excited states in large-sized challenging chemical systems.

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter (α) was set to recover the exact dissociation limit of H2 in a minimal basis set. Alternatively α can be viewed as a regularization parameter, where each value of α represents a valid variant of BW-s2, which we denote as BW-s2(α). In this work, we semi-empirically optimize α for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal α parameter (α = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2(α) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2(α) is iterative O(N5), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

Strong External Electric Fields Reduce Explosive Sensitivity: A Theoretical Investigation into the Reaction Selectivity in NH2NO2∙∙∙NH3

Controlling the selectivity of a detonation initiation reaction of explosive is essential to reduce sensitivity, and it seems impossible to reduce it by strengthening the external electric field. To verify this, the effects of external electric fields on the initiation reactions in NH2NO2∙∙∙NH3, a model system of the nitroamine explosive with alkaline additive, were investigated at the MP2/6-311++G(2d,p) and CCSD(T)/6-311++G(2d,p) levels. The concerted effect in the intermolecular hydrogen exchange is characterized by an index of the imaginary vibrations. Due to the weakened concerted effects by the electric field along the −x-direction opposite to the “reaction axis”, the dominant reaction changes from the intermolecular hydrogen exchange to 1,3-intramolecular hydrogen transference with the increase in the field strengths. Furthermore, the stronger the field strengths, the higher the barrier heights become, indicating the lower sensitivities. Therefore, by increasing the field strength and adjusting the orientation between the field and “reaction axis”, not only can the reaction selectivity be controlled, but the sensitivity can also be reduced, in particular under a super-strong field. Thus, a traditional concept, in which the explosive is dangerous under the super-strong external electric field, is theoretically broken. Compared to the neutral medium, a low sensitivity of the explosive with alkaline can be achieved under the stronger field. Employing atoms in molecules, reduced density gradient, and surface electrostatic potentials, the origin of the reaction selectivity and sensitivity change is revealed. This work provides a new idea for the technical improvement regarding adding the external electric field into the explosive system.

Accurate Property Prediction by Second Order Perturbation Theory: The REMP and OO-REMP Hybrids

The prediction of molecular properties such as equilibrium structures or vibrationalwavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing the REMP and OO-REMP hybrid perturbation theories showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers as well as dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations.OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 pm and 0.5 pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20 cm−1 . Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3 %. Across all investigated properties it turns out that a 20 %:80 % MP:RE mixing ratio consistently provides the best results. This is in line with our previous findings featuring closed and open shell reaction energies.

UREMP, RO-REMP, and OO-REMP: Hybrid perturbation theories for open-shell electronic structure calculations

An accurate description of the electron correlation energy in closed- and open-shell molecules is shown to be obtained by a second-order perturbation theory (PT) termed REMP. REMP is a hybrid of the Retaining the Excitation degree (RE) and the Møller–Plesset (MP) PTs. It performs particularly encouragingly in an orbital-optimized variant (OO-REMP) where the reference wavefunction is given by an unrestricted Slater determinant whose spin orbitals are varied such that the total energy becomes a minimum. While the approach generally behaves less satisfactorily with unrestricted Hartree–Fock references, reasonable performance is observed for restricted Hartree–Fock and restricted open-shell Hartree–Fock references. Inclusion of single excitations to OO-REMP is investigated and found—as in similar investigations—to be dissatisfying as it deteriorates performance. For the non-multireference subset of the accurate W4-11 benchmark set of Karton et al. [Chem. Phys. Lett. 510, 165–178 (2011)], OO-REMP predicts most atomization and reaction energies with chemical accuracy (1 kcal mol−1) if complete-basis-set extrapolation with augmented and core-polarized basis sets is used. For the W4-11 related test-sets, the error estimates obtained with the OO-REMP method approach those of coupled-cluster with singles, doubles and perturbative triples [CCSD(T)] within 20%–35%. The best performance of OO-REMP is found for a mixing ratio of 20%:80% MP:RE, which is essentially independent of whether radical stabilization energies, barrier heights, or reaction energies are investigated. Orbital optimization is shown to improve the REMP approach for both closed and open shell cases and outperforms coupled-cluster theory with singles and doubles (CCSD), spin-component scaled Møller-Plesset theory at second order (SCS-MP2), and density functionals, including double hybrids in all the cases considered.

Efficient implementations of the symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method with the density-fitting approximation

Efficient implementations of the symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles (OCCD) method with the density-fitting approach, denoted by DF-OCCD(T) and DF-OCCD(T)_Λ, are presented. The computational cost of the DF-OCCD(T) method is compared with that of the conventional OCCD(T). In the conventional OCCD(T) and OCCD(T)_Λ methods, one needs to perform four-index integral transformations at each coupled-cluster doubles iterations, which limits its applications to large chemical systems. Our results demonstrate that DF-OCCD(T) provides dramatically lower computational costs compared to OCCD(T), and there are more than 68-fold reductions in the computational time for the C₅H₁₂ molecule with the cc-pVTZ basis set. Our results show that the DF-OCCD(T) and DF-OCCD(T)_Λ methods are very helpful for the study of single bond-breaking problems. Performances of the DF-OCCD(T) and DF-OCCD(T)_Λ methods are noticeably better than that of the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for the potential energy surfaces of the molecules considered. Specifically, the DF-OCCD(T)_Λ method provides dramatic improvements upon CCSD(T), and there are 8-14-fold reductions in nonparallelity errors. Overall, we conclude that the DF-OCCD(T)_Λ method is very promising for the study of challenging chemical systems, where the CCSD(T) fails.

Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals

This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N⁶) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.

OO-REMP: Approaching Chemical Accuracy with Second-Order Perturbation Theory

We present a perturbation theory (PT) providing second-order energies that reproduce main group chemistry benchmark sets for reaction energies, barrier heights, and atomization energies with mean absolute deviations below 1 kcal mol^-1. The PT is defined as a constrained mixture of the unperturbed Hamiltonians of the Retaining the Excitation degree (RE) and the Møller-Plesset (MP) PTs. The orbitals of the reference wave function, a single unrestricted Slater determinant, are iteratively optimized to minimize the total energy. For all benchmark sets, good and near optimal performance of OO-REMP was observed for an unperturbed Hamiltonian consisting of 25% MP and 75% RE contributions.

Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction

While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory (κ-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of κ-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The aug-cc-pwCVTZ basis set is used for all calculations. The use of κ-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing two multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a κ-OOMP2 orbital reference (CCSD(T):κ-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm^-1, and 8.50 cm^-1, and 8.75 cm^-1 respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets shows these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. CCSD(T) with RHF orbitals is seen to improve upon UHF for the closed-shell species, while spatial symmetry breaking in a number of restricted open-shell HF (ROHF) references leads CCSD(T) with ROHF reference orbitals to exhibit the poorest statistical performance of all methods surveyed for open-shell species. The use of κ-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T).

Energy and analytic gradients for the orbital-optimized coupled-cluster doubles method with the density-fitting approximation: An efficient implementation

Efficient implementations of the orbital-optimized coupled-cluster doubles (or simply "optimized CCD," OCCD, for short) method and its analytic energy gradients with the density-fitting (DF) approach, denoted by DF-OCCD, are presented. In addition to the DF approach, the Cholesky-decomposed variant (CD-OCCD) is also implemented for energy computations. The computational cost of the DF-OCCD method (available in a plugin version of the DFOCC module of PSI4) is compared with that of the conventional OCCD (from the Q-CHEM package). The OCCD computations were performed with the Q-CHEM package in which OCCD are denoted by OD. In the conventional OCCD method, one needs to perform four-index integral transformations at each of the CCD iterations, which limits its applications to large chemical systems. Our results demonstrate that DF-OCCD provides dramatically lower computational costs compared to OCCD, and there are almost eightfold reductions in the computational time for the C₆H₁₄ molecule with the cc-pVTZ basis set. For open-shell geometries, interaction energies, and hydrogen transfer reactions, DF-OCCD provides significant improvements upon DF-CCD. Furthermore, the performance of the DF-OCCD method is substantially better for harmonic vibrational frequencies in the case of symmetry-breaking problems. Moreover, several factors make DF-OCCD more attractive compared to CCSD: (1) for DF-OCCD, there is no need for orbital relaxation contributions in analytic gradient computations; (2) active spaces can readily be incorporated into DF-OCCD; (3) DF-OCCD provides accurate vibrational frequencies when symmetry-breaking problems are observed; (4) in its response function, DF-OCCD avoids artificial poles; hence, excited-state molecular properties can be computed via linear response theory; and (5) symmetric and asymmetric triples corrections based on DF-OCCD [DF-OCCD(T)] have a significantly better performance in near degeneracy regions.

Third-Order Møller-Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error [Bertels, L. W.; J. Phys. Chem. Lett. 2019, 10, 4170 4176]. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N⁶) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

Modified Interaction-Strength Interpolation Method as an Important Step toward Self-Consistent Calculations

The modified point charge plus continuum (mPC) model [ConstantinL. A.; Phys. Rev. B2019, 99, 085117] solves the important failures of the original counterpart, namely, the divergences when the reduced gradient of the density is large, such as in the tail of the density and in quasi-dimensional density regimes. The mPC allows us to define a modified interaction-strength interpolation (mISI) method inheriting these good features, which are important steps toward the full self-consistent treatment. Here, we provide an assessment of mISI for molecular systems (i.e., considering thermochemistry properties, correlation energies, vertical ionization potentials, and several noncovalent interactions), harmonium atoms, and functional derivatives in the strong-interaction limit. For all our tests, mISI provides a systematic improvement over the original ISI method. Semilocal approximations of the second-order Görling–Levy (GL2) perturbation theory are also considered in the mISI method, showing considerable worsening of the results. Possible further development of mISI is briefly discussed.

Atom-Based Bootstrap Embedding For Molecules

Recent developments in quantum embedding have offered an attractive approach to describing electron correlation in molecules. However, previous methods such as density matrix embedding theory (DMET) require rigid partitioning of the system into fragments, which creates significant ambiguity for molecules. Bootstrap embedding (BE) is more flexible because it allows overlapping fragments, but when done on an orbital-by-orbital basis, BE introduces ambiguity in defining the connectivity of the orbitals. In this Letter, we present an atom-based fragment definition that significantly augments BE's performance in molecules. The resulting method, which we term atom-based BE, is very effective at recovering valence electron correlation in moderate-sized bases and delivers near-chemical-accuracy results using extrapolation. We anticipate atom-based BE may lead to a low-scaling and highly accurate approach to electron correlation in large molecules.

Third-Order Møller-Plesset Perturbation Theory Made Useful? Choice of Orbitals and Scaling Greatly Improves Accuracy for Thermochemistry, Kinetics, and Intermolecular Interactions

We develop and test methods that include second- and third-order perturbation theory (MP3) using orbitals obtained from regularized orbital-optimized second-order perturbation theory, κ-OOMP2, denoted as MP3:κ-OOMP2. Testing MP3:κ-OOMP2 shows RMS errors that are 1.7-5 times smaller than those of MP3 across 7 data sets. To do still better, empirical training of the scaling factors for the second- and third-order correlation energies and the regularization parameter on one of those data sets led to an unregularized scaled (c₂ = 1.0; c₃ = 0.8) denoted as MP2.8:κ-OOMP2. MP2.8:κ-OOMP2 yields significant additional improvement over MP3:κ-OOMP2 in 4 of 6 test data sets on thermochemistry, kinetics, and noncovalent interactions. Remarkably, these two methods outperform coupled cluster with singles and doubles in 5 of the 7 data sets considered, at greatly reduced cost (no O(N6) iterations).

Analytic energy gradients for orbital-optimized MP3 and MP2.5 with the density-fitting approximation: An efficient implementation

Efficient implementations of analytic gradients for the orbital-optimized MP3 and MP2.5 and their standard versions with the density-fitting approximation, which are denoted as DF-MP3, DF-MP2.5, DF-OMP3, and DF-OMP2.5, are presented. The DF-MP3, DF-MP2.5, DF-OMP3, and DF-OMP2.5 methods are applied to a set of alkanes and noncovalent interaction complexes to compare the computational cost with the conventional MP3, MP2.5, OMP3, and OMP2.5. Our results demonstrate that density-fitted perturbation theory (DF-MP) methods considered substantially reduce the computational cost compared to conventional MP methods. The efficiency of our DF-MP methods arise from the reduced input/output (I/O) time and the acceleration of gradient related terms, such as computations of particle density and generalized Fock matrices (PDMs and GFM), solution of the Z-vector equation, back-transformations of PDMs and GFM, and evaluation of analytic gradients in the atomic orbital basis. Further, application results show that errors introduced by the DF approach are negligible. Mean absolute errors for bond lengths of a molecular set, with the cc-pCVQZ basis set, is 0.0001-0.0002 Å. © 2017 Wiley Periodicals, Inc.

Orbital-optimized linearized coupled-cluster doubles with density-fitting and Cholesky decomposition approximations: an efficient implementation

An efficient implementation of the orbital-optimized linearized coupled-cluster double method with the density-fitting (DF-OLCCD) and Cholesky decomposition (CD-OLCCD) approximations is presented.