The MP2-F12 method in the Turbomole program package

Development of a multi-step screening procedure for redox active molecules in organic radical polymer anodes and as redox flow anolytes

Benzo[d]-X-zolyl-pyridinyl (XO, S, NH) radicals represent a promising class of redox-active molecules for organic batteries. We present a multistep screening procedure to identify the most promising radical candidates. Experimental investigations and highly correlated wave function-based calculations are performed to determine benchmark redox potentials. Based on these, the accuracies of different methods (semi-empirical, density functional theory, wave function-based), solvent models, dispersion corrections, and basis sets are evaluated. The developed screening procedure consists of three steps: First, a conformer search is performed with CREST. The molecules are selected based on the redox potentials calculated using GFN2-xTB. Second, HOMO energies calculated with reparametrized B3LYP-D3(BJ) and the def2-SVP basis set are used as selection criteria. The final molecules are selected based on the redox potentials calculated from Gibbs energies using BP86-D3(BJ)/def2-TZVP. With this multistep screening approach, promising molecules can be suggested for synthesis, and structure-property relationships can be derived.

Efficient Exploitation of Numerical Quadrature with Distance-Dependent Integral Screening in Explicitly Correlated F12 Theory: Linear Scaling Evaluation of the Most Expensive RI-MP2-F12 Term

We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially. The accuracy of our numerical quadrature/CABS-RI approach and the corresponding integral screening is thoroughly assessed for interaction and isomerization energies across a variety of numerical integration grids. Our method outperforms the standard density fitting/CABS-RI approach with errors below 1 μEh even for small grid sizes and moderate screening thresholds. The choice of the grid size and screening threshold allows us to tailor our ansatz to a desired accuracy and computational efficiency. We showcase the approach's effectiveness for the chemically relevant system valinomycin, employing a triple-ζ F12 basis set combination (C54H90N6O18, 5757 AO basis functions, 10,266 CABS basis functions, 735,783 grid points). In this context, our ansatz achieves higher accuracy combined with a 135× speedup compared to the classical density fitting based variant, requiring notably less computation time than the corresponding RI-MP2 calculation. Additionally, we demonstrate near-linear scaling through calculations on linear alkanes. We achieved an 817-fold acceleration for C80H162 and an extrapolated 28,765-fold acceleration for C200H402, resulting in a substantially reduced computational time for the latter─from 229 days to just 11.5 min. Our ansatz may also be adapted to the remaining MP2-F12 terms, which will be the subject of future work.

Basis Set Limit of CCSD(T) Energies: Explicit Correlation Versus Density-Based Basis-Set Correction

A thorough comparison is carried out for explicitly correlated and density-based basis-set correction approaches, which were primarily developed to mitigate the basis-set incompleteness error of wave function methods. An efficient implementation of the density-based scheme is also presented, utilizing the density-fitting approximation. The performance of these approaches is comprehensively tested for the second-order Møller-Plesset (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] methods with respect to the corresponding complete basis set references. It is demonstrated that the density-based correction together with complementary auxiliary basis set (CABS)-corrected Hartree-Fock energies is highly robust and effectively reduces the error of the standard approaches; however, it does not outperform the corresponding explicitly correlated methods. Nevertheless, what still makes the density-corrected CCSD and CCSD(T) methods competitive is that their computational costs are roughly half of those of the corresponding explicitly correlated variants. Additionally, an incremental approach for standard CCSD and CCSD(T) is introduced. In this simple scheme, the total energies are corrected with the CABS correction and explicitly correlated MP2 contributions. As demonstrated, the resulting methods yield surprisingly good results, below 1 kcal/mol for thermochemical properties even with a double-ζ basis, while their computational expenses are practically identical to those of the density-based basis-set correction approaches.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Pyrolysis of Methyl Formate and the Reaction of Methyl Formate with H Atoms: Shock Tube Experiments and Statistical Rate Theory

Methyl formate (MF) is the smallest carboxylic ester and currently considered a promising alternative fuel. It can also serve as a model compound to study the combustion chemistry of the ester group, which is a typical structural feature in many biodiesel components. In the present work, the pyrolysis of MF was investigated behind reflected shock waves at temperatures between 1430 and 2070 K at a nominal pressure of 1.1 bar. Both time-resolved hydrogen atom resonance absorption spectroscopy (H-ARAS) and time-resolved time-of-flight mass spectrometry (TOF-MS) were used for species detection. Additionally, the reaction of MF and perdeuterated MF-d₄ with H atoms was investigated at temperatures between 1000 and 1300 K at nominal pressures of 0.4 and 1.1 bar with H-ARAS. In the latter experiments, ethyl iodide served as precursor for H atoms. Rate coefficients of seven parallel unimolecular decomposition channels of MF and five parallel reaction channels of the MF + H reaction were calculated from statistical rate theory on the basis of molecular and transition state data from quantum chemical calculations. These calculated rate coefficients were implemented into an MF pyrolysis/oxidation mechanism from the literature, and the experimental concentration-time profiles of H (from ARAS) as well as MF, CH₃OH, HCHO, and CO (from TOF-MS) were modeled. It turned out that the literature mechanism, which was originally validated against flow-reactor experiments, ignition delay times, and laminar burning velocities, was generally able to fit also the concentration-time profiles from the shock tube experiments reasonably well. The agreement could still be improved by substituting the original rate coefficients, which were estimated from structure-reactivity relationships, by the values calculated from statistical rate theory in the present work. Details of the channel branching are discussed, and the updated mechanism is given, also in machine-readable form.

Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach

Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).

Massive Assessment of the Binding Energies of Atmospheric Molecular Clusters

Quantum chemical studies of the formation and growth of atmospheric molecular clusters are important for understanding aerosol particle formation. However, the search for the lowest free-energy cluster configuration is extremely time consuming. This makes high-level benchmark data sets extremely valuable in the quest for the global minimum as it allows the identification of cost-efficient computational methodologies, as well as the development of high-level machine learning (ML) models. Herein, we present a highly versatile quantum chemical data set comprising a total of 11 749 (acid)_1-2(base)_1-2 cluster configurations, containing up to 44 atoms. Utilizing the LUMI supercomputer, we calculated highly accurate PNO-CCSD(F12*)(T)/cc-pVDZ-F12 binding energies of the full set of cluster configurations leading to an unprecedented data set both in regard to sheer size and with respect to the level of theory. We employ the constructed benchmark set to assess the performance of various semiempirical and density functional theory methods. In particular, we find that the r²-SCAN-3c method shows excellent performance across the data set related to both accuracy and CPU time, making it a promising method to employ during cluster configurational sampling. Furthermore, applying the data sets, we construct ML models based on Δ-learning and provide recommendations for future application of ML in cluster configurational sampling.

Reduced-Scaling Double Hybrid Density Functional Theory with Rapid Basis Set Convergence through Localized Pair Natural Orbital F12

Following earlier work [Mehta, N.; Martin, J. M. L. J. Chem. Theory Comput.2022, 10.1021/acs.jctc.2c00426] that showed how the slow basis set convergence of the double hybrid density functional theory can be obviated by the use of F12 explicit correlation in the GLPT2 step (second order Görling-Levy perturbation theory), we demonstrate here for the very large and chemically diverse GMTKN55 benchmark suite that the CPU time scaling of this step can be reduced (asymptotically linearized) using the localized pair natural orbital (PNO-L) approximation at negligible cost in accuracy.

Accurate computed singlet-triplet energy differences for cobalt systems: implication for two-state reactivity

Accurate singlet-triplet energy differences for cobalt and rhodium complexes were calculated by using several wave function methods, such as MRCISD, CASPT2, CCSD(T) and BCCD(T). Relaxed energy differences were obtained by considering the singlet and triplet complexes, each at the minimum of their potential energy surfaces. Active spaces for multireference calculations were carefully checked to provide accurate results. The considered systems are built by increasing progressively the first coordination sphere around the metal. We included in our set two CpCoX complexes (Cp = cyclopentadienyl, X = alkenyl ligand), which have been suggested as intermediates in cycloaddition reactions. Indeed, cobalt systems have been used for more than a decade as active species in this kind of transformations, for which a two-state reactivity has been proposed. Most of the considered systems display a triplet ground state. However, in the case of a reaction intermediate, while a triplet ground state was predicted on the basis of Density Functional Theory results, our calculations suggest a singlet ground state. This stems from the competition between the exchange term (stabilising the triplet) and the accessibility of an intramolecular coordination (stabilising the singlet). This finding has an impact on the general mechanism of the cycloaddition reaction. Analogous rhodium systems were also studied and, as expected, they have a larger tendency to electron pairing than cobalt species.

A litmus test for the balanced description of dispersion interactions and coordination chemistry of lanthanoids

The influence of long-range interactions on the structure of complexes of Eu(III) with four 9-hydroxy-phenalen-1-one ligands (HPLN) and one alkaline earth metal dication [Eu(PLN)₄AE]⁺ (AE: Mg, Ca, Sr, and Ba) is analyzed. Through the [Eu(PLN)₄Ca]⁺ complex, which is a charged complex with two metals-one of them a lanthanoid-and with four relatively fluxional π-ligands, the difficulties of describing such systems are identified. The inclusion of the D3(BJ) or D4 corrections to different density functionals introduces significant changes in the structure, which are shown to stem from the interaction between pairs of PLN ligands. This interaction is studied further with a variety of density functionals, wave-function based methods, and by means of the random phase approximation. By comparing the computed results with those from experimental evidence of gas-phase photoluminescence and ion mobility measurements it is concluded that the inclusion of dispersion corrections does not always yield structures that are in agreement with the experimental findings.

Highly Efficient and Accurate Computation of Multiple Orbital Spaces Spanning Fock Matrix Elements on Central and Graphics Processing Units for Application in F12 Theory

We employ our recently published highly efficient seminumerical exchange (sn-LinK) [Laqua, H.; Thompson, T. H.; Kussmann, J.; Ochsenfeld, C. J. Chem. Theory Comput. 2020, 16, 1456-1468] and integral-direct resolution of the identity Coulomb (RI-J) [Kussmann, J.; Laqua, H.; Ochsenfeld, C. J. Chem. Theory Comput. 2021, 17, 1512-1521] methods to significantly accelerate the computation of the demanding multiple orbital spaces spanning Fock matrix elements present in R12/F12 theory on central and graphics processing units. The errors introduced by RI-J and sn-LinK into the RI-MP2-F12 energy are thoroughly assessed for a variety of basis sets and integration grids. We find that these numerical errors are always below "chemical accuracy" (∼1 mH) even for the coarsest settings and can easily be reduced below 1 μH by employing only moderately large integration grids and RI-J basis sets. Since the number of basis functions of the multiple orbital spaces is notably larger compared with conventional Hartree-Fock theory, the efficiency gains from the superior basis scaling of RI-J and sn-LinK (O(N_bas²) instead of O(N_bas⁴) for both) are even more significant, with maximum speedup factors of 37 000 for RI-J and 4500 for sn-LinK. In total, the multiple orbital spaces spanning Fock matrix evaluation of the largest tested structure using a triple-ζ F12 basis set (5058 AO basis functions, 9267 CABS basis functions) is accelerated over 1575× using CPUs and over 4155× employing GPUs.

The unimolecular decomposition of dimethoxymethane: channel switching as a function of temperature and pressure

Branching ratios of competing unimolecular reactions often exhibit a complicated temperature and pressure dependence that makes modeling of complex reaction systems in the gas phase difficult. In particular, the competition...

Tailoring the Gibbs Free Energy of the Nitroxide Exchange Reaction: Substituent and Solvent Effects

A quantum chemical study of the nitroxide exchange reaction is presented. Inspired by the recent use of this reaction in the synthesis of dynamic covalent polymer networks, we studied the influence of substituents and solvents on the Gibbs free energy, which plays a crucial role for both the reversibility of the reaction and the extent of product formation. We provide accurate benchmark values based on CCSD(T) and COSMO-RS theory for a series of structural modifications and make suggestions for improving the molecular building blocks used so far.

Quantifying Conformational Isomerism in Chain Molecules by Linear Raman Spectroscopy: The Case of Methyl Esters

The conformational preferences of the ester group have the potential to facilitate the large amplitude folding of long alkyl chains in the gas phase. They are monitored by Raman spectroscopy in supersonic jet expansions for the model system methyl butanoate, after establishing a quantitative relationship with quantum-chemical predictions for methyl methanoate. This requires a careful analysis of experimental details, and a simulation of the rovibrational contours for near-symmetric top molecules. The technique is shown to be complementary to microwave spectroscopy in quantifying coexisting conformations. It confirms that a C-O-C(=O)-C-C chain segment can be collapsed into a single all-trans conformation by collisional cooling, whereas alkyl chain isomerism beyond this five-membered chain largely survives the jet expansion. This sets the stage for the investigation of linear alkyl alkanoates in terms of dispersion-induced stretched-chain to hairpin transitions by Raman spectroscopy.

Size-consistent explicitly correlated triple excitation correction

A new approach is proposed to reduce the basis set incompleteness error of the triple excitation correction in explicitly correlated coupled-cluster singles and doubles with perturbative triples calculations. Our method is similar to the intuitive triples correction approach of Knizia et al. [J. Chem. Phys. 130, 054104 (2009)] but, in contrast to the latter, is size-consistent. The new approximation is easy to implement, and its overhead is negligible with respect to the conventional (T) correction. The performance of the approach is assessed for atomization, reaction, and interaction energies as well as for bond lengths and harmonic vibrational frequencies. The advantages of its size consistency are also demonstrated.

Focal-point approach with pair-specific cusp correction for coupled-cluster theory

We present a basis set correction scheme for the coupled-cluster singles and doubles (CCSD) method. The scheme is based on employing frozen natural orbitals (FNOs) and diagrammatically decomposed contributions to the electronic correlation energy, which dominate the basis set incompleteness error (BSIE). As recently discussed in the work of Irmler et al. [Phys. Rev. Lett. 123, 156401 (2019)], the BSIE of the CCSD correlation energy is dominated by the second-order Møller-Plesset (MP2) perturbation energy and the particle-particle ladder term. Here, we derive a simple approximation to the BSIE of the particle-particle ladder term that effectively corresponds to a rescaled pair-specific MP2 BSIE, where the scaling factor depends on the spatially averaged correlation hole depth of the coupled-cluster and first-order pair wavefunctions. The evaluation of the derived expressions is simple to implement in any existing code. We demonstrate the effectiveness of the method for the uniform electron gas. Furthermore, we apply the method to coupled-cluster theory calculations of atoms and molecules using FNOs. Employing the proposed correction and an increasing number of FNOs per occupied orbital, we demonstrate for a test set that rapidly convergent closed and open-shell reaction energies, atomization energies, electron affinities, and ionization potentials can be obtained. Moreover, we show that a similarly excellent trade-off between required virtual orbital basis set size and remaining BSIEs can be achieved for the perturbative triples contribution to the CCSD(T) energy employing FNOs and the (T*) approximation.

Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments

Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.

A scaled explicitly correlated F12 correction to second-order Møller-Plesset perturbation theory

An empirically scaled version of the explicitly correlated F12 correction to second-order Møller-Plesset perturbation theory (MP2-F12) is introduced. The scaling eliminates the need for many of the most costly terms of the F12 correction while reproducing the unscaled explicitly correlated F12 interaction energy correction to a high degree of accuracy. The method requires a single, basis set dependent scaling factor that is determined by fitting to a set of test molecules. We present factors for the cc-pVXZ-F12 (X = D, T, Q) basis set family obtained by minimizing interaction energies of the S66 set of small- to medium-sized molecular complexes and show that our new method can be applied to accurately describe a wide range of systems. Remarkably good explicitly correlated corrections to the interaction energy are obtained for the S22 and L7 test sets, with mean percentage errors for the double-zeta basis of 0.60% for the F12 correction to the interaction energy, 0.05% for the total electron correlation interaction energy, and 0.03% for the total interaction energy, respectively. Additionally, mean interaction energy errors introduced by our new approach are below 0.01 kcal mol^-1 for each test set and are thus negligible for second-order perturbation theory based methods. The efficiency of the new method compared to the unscaled F12 correction is shown for all considered systems, with distinct speedups for medium- to large-sized structures.

Basis set extrapolation in pair natural orbital theories

We present the results of a benchmark study of the effect of Pair Natural Orbital (PNO) truncation errors on the performance of basis set extrapolation. We find that reliable conclusions from the application of Helgaker's extrapolation method are only obtained when using tight PNO thresholds of at least 10^-7. The use of looser thresholds introduces a significant risk of observing a false basis set convergence and underestimating the residual basis set errors. We propose an alternative extrapolation approach based on the PNO truncation level that only requires a single basis set and show that it is a viable alternative to hierarchical basis set extrapolation methods.

On the influence of water on urea condensation reactions: a theoretical study

The influence of water molecules on the kinetics of urea condensation reactions was studied with high-level quantum chemical methods and statistical rate theory. The study focuses on the production of biuret, triuret, and cyanuric acid from urea because of their relevance as unwanted byproducts in the urea-based selective catalytic reduction (urea-SCR) exhaust after treatment of Diesel engines. In order to characterize the potential energy surfaces and molecular reaction pathways, calculations with explicitly-correlated coupled-cluster methods were performed. It turned out that the reactions proceed via pre-reactive complexes and the inclusion of one or two water molecules into the condensation mechanisms leads to a decrease of the energy barriers. This effect is particularly pronounced in the production of biuret. Due to the pre-reactive equilibria, the rates of the overall reactions can increase or decrease by incorporating water into the mechanism, depending on the temperature and water concentration. Under the conditions of urea-SCR, the studied reactions are too slow to contribute to the observed byproduct formation.

2-Methoxyethanol: harmonic tricks, anharmonic challenges and chirality-sensitive chain aggregation

The Raman spectrum of jet-cooled, nearly monoconformational 2-methoxyethanol reveals a strong, aggregation-sensitive resonance in the backbone stretching region and pronounced chirality-dependent dimerisation effects in the OH stretching region. These effects are disentangled by a dual-detection scheme which simultaneously probes the two spectral regions. The performance of mass-scaled harmonic wavenumber predictions and of vibrational perturbation theory in reproducing these spectral features is explored.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Assessment of the DLPNO Binding Energies of Strongly Noncovalent Bonded Atmospheric Molecular Clusters

This work assesses the performance of DLPNO-CCSD(T0), DLPNO-MP2, and density functional theory methods in calculating the binding energies of a representative test set of 45 atmospheric acid-acid, acid-base, and acid-water dimer clusters. The performance of the approximate methods is compared to high level explicitly correlated CCSD(F12*)(T)/complete basis set (CBS) reference calculations. Out of the tested density functionals, ωB97X-D3(BJ) shows the best performance with a mean deviation of 0.09 kcal/mol and a maximum deviation of 0.83 kcal/mol. The RI-CC2/aug-cc-pV(T+d)Z level of theory severely overpredicts the cluster binding energies with a mean deviation of -1.31 kcal/mol and a maximum deviation up to -3.00 kcal/mol. Hence, RI-CC2/aug-cc-pV(T+d)Z should not be utilized for studying atmospheric molecular clusters. The DLPNO variants are tested both with and without the inclusion of explicit correlation (F12) in the wavefunction, with different pair natural orbital (PNO) settings (loosePNO, normalPNO, and tightPNO) and using both double and triple zeta basis sets. The performance of the DLPNO-MP2 methods is found to be independent of PNO settings and yield low mean deviations of -0.84 kcal/mol or below. However, DLPNO-MP2 requires explicitly correlated wavefunctions to yield maximum deviations below 1.40 kcal/mol. For obtaining high accuracy, with maximum deviation below ∼1.0 kcal/mol, either DLPNO-CCSD(T0)/aug-cc-pVTZ (normalPNO) calculations or DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 (normalPNO) calculations are required. The most accurate level of theory is found to be DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 using a tightPNO criterion which yields a mean deviation of 0.10 kcal/mol, with a maximum deviation of 0.20 kcal/mol, compared to the CCSD(F12*)(T)/CBS reference.

Temperature- and pressure-dependent kinetics of the competing C–O bond fission reactions of dimethoxymethane

Under typical shock tube conditions, dimethoxymethane decomposes mainly to give CH₃ + OCH₂OCH₃.

Efficient and accurate description of adsorption in zeolites

Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.

Explicitly Correlated Dispersion and Exchange Dispersion Energies in Symmetry-Adapted Perturbation Theory

The individual interaction energy terms in symmetry-adapted perturbation theory (SAPT) not only have different physical interpretations but also converge to their complete basis set (CBS) limit values at quite different rates. Dispersion energy is notoriously the slowest converging interaction energy contribution, and exchange dispersion energy, while smaller in absolute value, converges just as poorly in relative terms. To speed up the basis set convergence of the lowest-order SAPT dispersion and exchange dispersion energies, we borrow the techniques from explicitly correlated (F12) electronic structure theory and develop practical expressions for the closed-shell E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 contributions. While the latter term has been derived and implemented for the first time, the former correction was recently proposed by Przybytek [ J. Chem. Theory Comput. 2018 , 14 , 5105 - 5117 ] using an Ansatz with a full optimization of the explicitly correlated amplitudes. In addition to reimplementing the fully optimized variant of E_disp⁽²⁰⁾-F12, we propose three approximate Ansätze that substantially improve the scaling of the method and at the same time avoid the numerical instabilities of the unrestricted optimization. The performance of all four resulting flavors of E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 is first tested on helium, neon, argon, water, and methane dimers, with orbital and auxiliary basis sets up to aug-cc-pV5Z and aug-cc-pV5Z-RI, respectively. The double- and triple-ζ basis set calculations are then extended to the entire A24 database of noncovalent interaction energies and compared with CBS estimates for E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾ computed using conventional SAPT with basis sets up to aug-cc-pV6Z with midbond functions. It is shown that the F12 treatment is highly successful in improving the basis set convergence of the SAPT terms, with the F12 calculations in an X-tuple ζ basis about as accurate as conventional calculations in bases with cardinal numbers (X + 2) for E_disp⁽²⁰⁾ and either (X + 1) or (X + 2) for E_exch-disp⁽²⁰⁾. While the full amplitude optimization affords the highest accuracy for both corrections, the much simpler and numerically stable optimized diagonal Ansatz is a very close second. We have also tested the performance of the simple F12 correction based on the second-order Møller-Plesset perturbation theory, SAPT-F12(MP2) [ Frey , J. A. ; Chem. Rev. 2016 , 116 , 5614 - 5641 ] and observed that it is also quite successful in speeding up the basis set convergence of conventional E_disp⁽²⁰⁾ + E_exch-disp⁽²⁰⁾, albeit with some outliers.

High Accuracy Quantum Chemical and Thermochemical Network Data for the Heats of Formation of Fluorinated and Chlorinated Methanes and Ethanes

Reliable heats of formation are reported for numerous fluorinated and chlorinated methane and ethane derivatives by means of an accurate thermochemical protocol, which involves explicitly correlated coupled-cluster calculations augmented with anharmonic, scalar relativistic, and diagonal Born-Oppenheimer corrections. The theoretical results, along with additional experimental data, are further enhanced with the help of the thermochemical network approach. For 28 species, out of 50, this study presents the best estimates, and discrepancies with previous reports are also highlighted. Furthermore, the effects of the less accurate theoretical data on the results yielded by thermochemical networks are discussed.

Communication: Quasi-robust local density fitting

The requirement that the linear density fitting error in the integral exactly vanishes introduces unphysical long range contributions to the approximate density when the auxiliary basis is incomplete. A quasi-robust density fitting formulation is presented where spatial locality is recovered at the expense of permitting a linear error that is made small by the fitting procedure, which involves optimising the Coulomb potential of the approximate charge density. The method is shown to be stable and almost as accurate as standard robust density fitting without local approximations in practical calculations using standard density fitting basis sets.

Combining Accuracy and Efficiency: An Incremental Focal-Point Method Based on Pair Natural Orbitals

In this work, we present a new pair natural orbitals (PNO)-based incremental scheme to calculate CCSD(T) and CCSD(T0) reaction, interaction, and binding energies. We perform an extensive analysis, which shows small incremental errors similar to previous non-PNO calculations. Furthermore, slight PNO errors are obtained by using T_PNO = T_TNO with appropriate values of 10^-7 to 10^-8 for reactions and 10^-8 for interaction or binding energies. The combination with the efficient MP2 focal-point approach yields chemical accuracy relative to the complete basis-set (CBS) limit. In this method, small basis sets (cc-pVDZ, def2-TZVP) for the CCSD(T) part are sufficient in case of reactions or interactions, while some larger ones (e.g., (aug)-cc-pVTZ) are necessary for molecular clusters. For these larger basis sets, we show the very high efficiency of our scheme. We obtain not only tremendous decreases of the wall times (i.e., factors >10²) due to the parallelization of the increment calculations as well as of the total times due to the application of PNOs (i.e., compared to the normal incremental scheme) but also smaller total times with respect to the standard PNO method. That way, our new method features a perfect applicability by combining an excellent accuracy with a very high efficiency as well as the accessibility to larger systems due to the separation of the full computation into several small increments.

Coupled-Cluster in Real Space. 1. CC2 Ground State Energies Using Multiresolution Analysis

A framework to calculate CC2 approximated coupled-cluster ground state correlation energies in a multiresolution basis is derived and implemented into the MADNESS library. The CC2 working equations are formulated in first quantization which makes them suitable for real-space methods. The first quantized equations can be interpreted diagrammatically using the usual diagrams from second quantization with adjusted interpretation rules. Singularities arising from the nuclear and electronic potentials are regularized by explicitly taking the nuclear and electronic cusps into account. The regularized three- and six-dimensional cluster functions are represented directly on an adaptive grid. The resulting equations are free of singularities and virtual orbitals, which results in a low intrinsic scaling. Correlation energies close to the basis set limit are computed for small molecules. This work is the first step toward CC2 excitation energies in a multiresolution basis.

Noncovalent Interactions in Specific Recognition Motifs of Protein-DNA Complexes

In view of the importance of protein-DNA interactions in biological processes, we extracted from the Protein Data Bank several one-to-one complexes of amino acids with nucleotides that matched certain geometric and energetic specificity criteria and investigated them using quantum chemistry methods. The CCSD(T)/CBS interaction energies were used as a benchmark to compare the performance of the MP2.5, MP2-F12, DFT-D3, and PM6-D3H4 methods. All methods yielded good agreement with the reference values, with declining accuracy from MP2.5 to PM6-D3H4. Regardless of the site of interaction, the minima found after full optimization in implicit solvent with high dielectric constant were close to the structures experimentally detected in protein-DNA complexes. According to DFT-SAPT analysis, the nature of noncovalent interactions strongly depends on the type of amino acid. The negatively charged sugar-phosphate backbone of DNA heavily influences the strength of interactions and must be included in the computational model, especially in the case of interactions with charged amino acids.

Cooperative hydrogen bonds form a pseudocycle stabilizing an isolated complex of isocyanic acid with urea

The structure and properties of a gas phase complex formed between urea and isocyanic acid are characterised through microwave spectroscopy andab initiocalculations at the CCSD(T)(F12*)/aug-cc-pVTZ level.