Doubly hybrid density functional for accurate descriptions of nonbond interactions, thermochemistry, and thermochemical kinetics

Chemical Potentials and the One-Electron Hamiltonian of the Second-Order Perturbation Theory from the Functional Derivative Approach

We develop a functional derivative approach to calculate the chemical potentials of second-order perturbation theory (MP2). In the functional derivative approach, the correlation part of the MP2 chemical potential, which is the derivative of the MP2 correlation energy with respect to the occupation number of frontier orbitals, is obtained from the chain rule via the noninteracting Green's function. First, the MP2 correlation energy is expressed in terms of the noninteracting Green's function, and its functional derivative to the noninteracting Green's function is the second-order self-energy. Then, the derivative of the noninteracting Green's function to the occupation number is obtained by including the orbital relaxation effect. We show that the MP2 chemical potentials obtained from the functional derivative approach agree with that obtained from the finite difference approach. The one-electron Hamiltonian, defined as the derivative of the MP2 energy with respect to the one particle density matrix, is also derived using the functional derivative approach, which can be used in the self-consistent calculations of MP2 and double-hybrid density functionals. The developed functional derivative approach is promising for calculating the chemical potentials and the one-electron Hamiltonian of approximate functionals and many-body perturbation approaches dependent explicitly on the noninteracting Green's function.

Excited-State Charge Transfer Coupling from Quasiparticle Energy Density Functional Theory

The recently developed Quasiparticle Energy (QE) scheme, based on a DFT calculation with one more (or less) electron, offers a good description of excitation energies, even with charge transfer characters. In this work, QE is further extended to calculate electron transfer (ET) couplings involving two excited states. We tested it with a donor-acceptor complex, consisting of a furan and a 1,1-dicyanoethylene (DCNE), in which two low lying charge transfer and local excitation states are involved. With generalized Mülliken-Hush and fragment charge-difference schemes, couplings from the QE approach generally agree well with those obtained from TDDFT, except that QE couplings exhibit better exponential distance dependence. Couplings from half-energy gaps with an external field are also calculated and reported. Our results show that the QE scheme is robust in calculating ET couplings with greatly reduced computational time.

Accurate Electronic and Optical Properties of Organic Doublet Radicals Using Machine Learned Range-Separated Functionals

Luminescent organic semiconducting doublet-spin radicals are unique and emergent optical materials because their fluorescent quantum yields (Φfl) are not compromised by the spin-flipping intersystem crossing (ISC) into a dark high-spin state. The multiconfigurational nature of these radicals challenges their electronic structure calculations in the framework of single-reference density functional theory (DFT) and introduces room for method improvement. In the present study, we extended our earlier development of ML-ωPBE [J. Phys. Chem. Lett., 2021, 12, 9516-9524], a range-separated hybrid (RSH) exchange-correlation (XC) functional constructed using the stacked ensemble machine learning (SEML) algorithm, from closed-shell organic semiconducting molecules to doublet-spin organic semiconducting radicals. We assessed its performance for a new test set of 64 doublet-spin radicals from five categories while placing all previously compiled 3926 closed-shell molecules in the new training set. Interestingly, ML-ωPBE agrees with the nonempirical OT-ωPBE functional regarding the prediction of the molecule-dependent range-separation parameter (ω), with a small mean absolute error (MAE) of 0.0197 a0-1, but saves the computational cost by 2.46 orders of magnitude. This result demonstrates an outstanding domain adaptation capacity of ML-ωPBE for diverse organic semiconducting species. To further assess the predictive power of ML-ωPBE in experimental observables, we also applied it to evaluate absorption and fluorescence energies (Eabs and Efl) using linear-response time-dependent DFT (TDDFT), and we compared its behavior with nine popular XC functionals. For most radicals, ML-ωPBE reproduces experimental measurements of Eabs and Efl with small MAEs of 0.299 and 0.254 eV, only marginally different from those of OT-ωPBE. Our work illustrates a successful extension of the SEML framework from closed-shell molecules to doublet-spin radicals and will open the venue for calculating optical properties for organic semiconductors using single-reference TDDFT.

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Developing correlation-consistent numeric atom-centered orbital basis sets for krypton: Applications in RPA-based correlated calculations

Localized atomic orbitals are the preferred basis set choice for large-scale explicit correlated calculations, and high-quality hierarchical correlation-consistent basis sets are a prerequisite for correlated methods to deliver numerically reliable results. At present, numeric atom-centered orbital (NAO) basis sets with valence correlation consistency (VCC), designated as NAO-VCC-nZ, are only available for light elements from hydrogen (H) to argon (Ar) [Zhang et al., New J. Phys. 15, 123033 (2013)]. In this work, we extend this series by developing NAO-VCC-nZ basis sets for krypton (Kr), a prototypical element in the fourth row of the periodic table. We demonstrate that NAO-VCC-nZ basis sets facilitate the convergence of electronic total-energy calculations using the Random Phase Approximation (RPA), which can be used together with a two-point extrapolation scheme to approach the complete basis set limit. Notably, the Basis Set Superposition Error (BSSE) associated with the newly generated NAO basis sets is minimal, making them suitable for applications where BSSE correction is either cumbersome or impractical to do. After confirming the reliability of NAO basis sets for Kr, we proceed to calculate the Helmholtz free energy for Kr crystal at the theoretical level of RPA plus renormalized single excitation correction. From this, we derive the pressure-volume (P-V) diagram, which shows excellent agreement with the latest experimental data. Our work demonstrates the capability of correlation-consistent NAO basis sets for heavy elements, paving the way toward numerically reliable correlated calculations for bulk materials.

QSAR modelling, molecular docking, molecular dynamic and ADMET prediction of pyrrolopyrimidine derivatives as novel Bruton's tyrosine kinase (BTK) inhibitors

In recent years, there has been a focus on developing and discovering novel Bruton's tyrosine kinase (BTK) inhibitors, as they offer an effective treatment strategy for B-cell malignancies. BTK plays a crucial role in B cell receptor (BCR)-mediated activation and proliferation by regulating downstream factors such as the NF-κB and MAP kinase pathways. To address this challenge and propose potential therapeutic options for B-cell lymphomas, researchers conducted 2D-QSAR and ADMET studies on pyrrolopyrimidine derivatives that act as inhibitors of the BCR site in cytochrome b. These studies aim to improve and identify new compounds that could serve as more potent potential BTK inhibitors, which would lead to the identification of new drug candidates in this field. In our study, we used 2D-QSAR (multiple linear regression, multiple nonlinear regression, and artificial neural networks), molecular docking, molecular dynamics, and ADMET properties to investigate the potential of 35 pyrrolopyrimidine derivatives as BTK inhibitors. A molecular docking study and molecular dynamics simulations of molecule 13 over 10 ns revealed that it establishes multiple hydrogen bonds with several residues and exhibits frequent stability throughout the simulation period. Based on the results obtained by molecular modeling, we proposed six new compounds (Pred1, Pred2, Pred3, Pred4, Pred5, and Pred6) with highly significant predicted activity by MLR models. A study based on the in silico evaluation of the predicted ADMET properties of the new candidate molecules is strongly recommended to classify these molecules as promising candidates for new anticancer agents specifically designed to target Bruton's tyrosine kinase (BTK) inhibition.

Accurate fundamental invariant-neural network representation of ab initio potential energy surfaces

Highly accurate potential energy surfaces are critically important for chemical reaction dynamics. The large number of degrees of freedom and the intricate symmetry adaption pose a big challenge to accurately representing potential energy surfaces (PESs) for polyatomic reactions. Recently, our group has made substantial progress in this direction by developing the fundamental invariant-neural network (FI-NN) approach. Here, we review these advances, demonstrating that the FI-NN approach can represent highly accurate, global, full-dimensional PESs for reactive systems with even more than 10 atoms. These multi-channel reactions typically involve many intermediates, transition states, and products. The complexity and ruggedness of this potential energy landscape present even greater challenges for full-dimensional PES representation. These PESs exhibit a high level of complexity, molecular size, and accuracy of fit. Dynamics simulations based on these PESs have unveiled intriguing and novel reaction mechanisms, providing deep insights into the intricate dynamics involved in combustion, atmospheric, and organic chemistry.

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost

Noncovalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller-Plesset adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields cosκos-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies. With the N4 formal scaling, cosκos-SPL2 is competitive or often outperforms more expensive dispersion-corrected double hybrids for NCIs. The accuracy of cosκos-SPL2 particularly shines for anionic halogen bonded complexes, where it surpasses standard dispersion-corrected DFT by a factor of 3 to 5.

Reproducibility of density functional approximations: How new functionals should be reported

Density functional theory is the workhorse of chemistry and materials science, and novel density functional approximations are published every year. To become available in program packages, the novel density functional approximations (DFAs) need to be (re)implemented. However, according to our experience as developers of Libxc [Lehtola et al., SoftwareX 7, 1 (2018)], a constant problem in this task is verification due to the lack of reliable reference data. As we discuss in this work, this lack has led to several non-equivalent implementations of functionals such as Becke-Perdew 1986, Perdew-Wang 1991, Perdew-Burke-Ernzerhof, and Becke's three-parameter hybrid functional with Lee-Yang-Parr correlation across various program packages, yielding different total energies. Through careful verification, we have also found many issues with incorrect functional forms in recent DFAs. The goal of this work is to ensure the reproducibility of DFAs. DFAs must be verifiable in order to prevent the reappearance of the above-mentioned errors and incompatibilities. A common framework for verification and testing is, therefore, needed. We suggest several ways in which reference energies can be produced with free and open source software, either with non-self-consistent calculations with tabulated atomic densities or via self-consistent calculations with various program packages. The employed numerical parameters-especially the quadrature grid-need to be converged to guarantee a ≲0.1 μEh precision in the total energy, which is nowadays routinely achievable in fully numerical calculations. Moreover, as such sub-μEh level agreement can only be achieved when fully equivalent implementations of the DFA are used, the source code of the reference implementation should also be made available in any publication describing a new DFA.

A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database

In two recent papers [A. D. Becke, J. Chem. Phys. 156, 214101 (2022) and A. D. Becke, J. Chem. Phys. 157, 234102 (2022)], we compared two Kohn-Sham density functionals based on physical modeling and theory with the best density-functional power-series fits in the literature. The best error statistics reported to date for a hybrid functional on the general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)] were obtained. In the present work, additional second-order perturbation-theory terms are considered. The result is a 12-parameter double-hybrid density functional with the lowest GMTKN55 WTMAD2 "weighted total mean absolute deviation" error (1.76 kcal/mol) yet seen for any hybrid or double-hybrid density-functional approximation. We call it "DH23."

Comparison of the Performance of Density Functional Methods for the Description of Spin States and Binding Energies of Porphyrins

This work analyzes the performance of 250 electronic structure theory methods (including 240 density functional approximations) for the description of spin states and the binding properties of iron, manganese, and cobalt porphyrins. The assessment employs the Por21 database of high-level computational data (CASPT2 reference energies taken from the literature). Results show that current approximations fail to achieve the "chemical accuracy" target of 1.0 kcal/mol by a long margin. The best-performing methods achieve a mean unsigned error (MUE) <15.0 kcal/mol, but the errors are at least twice as large for most methods. Semilocal functionals and global hybrid functionals with a low percentage of exact exchange are found to be the least problematic for spin states and binding energies, in agreement with the general knowledge in transition metal computational chemistry. Approximations with high percentages of exact exchange (including range-separated and double-hybrid functionals) can lead to catastrophic failures. More modern approximations usually perform better than older functionals. An accurate statistical analysis of the results also casts doubts on some of the reference energies calculated using multireference methods. Suggestions and general guidelines for users are provided in the conclusions. These results hopefully stimulate advances for both the wave function and the density functional side of electronic structure calculations.

Discriminating and understanding molecular crystal polymorphism

Polymorph discrimination for a molecular crystal has long been a challenging task, which, nonetheless, is a major concern in the pharmaceutical industry. In this work, we have investigated polymorph discrimination on three different molecular crystals, tetrolic acid, oxalic acid, and oxalyl dihydrazide, covering both packing polymorphism and conformational polymorphism. To gain more understanding, we have performed energy decomposition analysis based on many-body expansion, and have compared the results from the XO-PBC method, that is, the eXtended ONIOM method (XO) with the periodic boundary condition (PBC), with those from some commonly used dispersion corrected density functional theory (DFT-D) methods. It is shown here that, with the XYG3 doubly hybrid functional chosen as the target high level to capture the intra- and short-range intermolecular interactions, and the periodic PBE as the basic low level to take long range interactions into account, the XO-PBC(XYG3:PBE) method not only obtains the correct experimental stability orderings, but also predicts reasonable polymorph energy ranges for all three cases. Our results have demonstrated the usefulness of the present theoretical methods, in particular XO-PBC, while highlighted the importance of a better treatment of different kinds of interactions to be beneficial to polymorph control.

Accurate descriptions of molecule-surface interactions in electrocatalytic CO2 reduction on the copper surfaces

Copper-based catalysts play a pivotal role in many industrial processes and hold a great promise for electrocatalytic CO₂ reduction reaction into valuable chemicals and fuels. Towards the rational design of catalysts, the growing demand on theoretical study is seriously at odds with the low accuracy of the most widely used functionals of generalized gradient approximation. Here, we present results using a hybrid scheme that combines the doubly hybrid XYG3 functional and the periodic generalized gradient approximation, whose accuracy is validated against an experimental set on copper surfaces. A near chemical accuracy is established for this set, which, in turn, leads to a substantial improvement for the calculated equilibrium and onset potentials as against the experimental values for CO₂ reduction to CO on Cu(111) and Cu(100) electrodes. We anticipate that the easy use of the hybrid scheme will boost the predictive power for accurate descriptions of molecule-surface interactions in heterogeneous catalysis.

Statistical Equivalence of Quantum Chemical Methods for Energy Distribution in Water Clusters

Methods of quantum chemistry are instrumental in understanding molecular structures and properties. However, the results demonstrate significant variability, which is difficult to predict and rationalize. The fundamental question is whether some molecular systems exhibit properties invariant with respect to the computational method. The idea explored here is that collective properties of statistical ensembles should be more robust than characteristics of individual molecules and their arbitrary sets. This effect is demonstrated for the complete set of hydrogen-bond topologies of the dodecahedral water cluster (H₂ O)₂₀ . Non-Gaussian energy distributions produced by various methods have the same functional form despite strong differences in mean values and standard deviations. The conclusion is tested on methods of different complexity and origin employing a number of criteria. A linear mapping between the energies produced by different methods is discussed. The significance of the results is in establishing a collective equivalence property of quantum chemical methods.

Structural and dynamical studies of CH-πbonded CH4-C6H6dimer by ultrafast intermolecular Coulombic decay

The inner-valence ionization and fragmentation dynamics of CH₄-C₆H₆dimer induced by 200 eV electron impact is studied utilizing a multi-particle coincidence momentum spectroscopy. The three-dimensional momentum vectors and kinetic energy release (KER) of the CH₄⁺+C₆H₆⁺ion pairs are obtained by coincident momentum measurement. Our analysis on the absolute cross sections indicates that the intermediate dication CH₄⁺-C₆H₆⁺is preferentially produced by the removal of an inner-valence electron from CH₄or C₆H₆and subsequent relaxation of ultrafast intermolecular Coulombic decay followed by two-body Coulomb explosion. Combining withab initiomolecular dynamics (AIMD) simulations, the real-time fragmentation dynamics including translational, vibrational and rotational motions are presented as a function of propagation time. The revealed fragmentation dynamics are expected to have a potential implication for crystal structure imaging with various radiation sources.

Excitation energies of polycylic aromatic hydrocarbons by double-hybrid functionals: Assessing the PBE0-DH and PBE-QIDH models and their range-separated versions

A family of non-empirical double-hybrid (DH) density functionals, such as Perdew-Burke-Ernzerhof (PBE)0-DH, PBE-QIDH, and their range-separated exchange (RSX) versions RSX-0DH and RSX-QIDH, all using Perdew-Burke-Ernzerhof(PBE) exchange and correlationfunctionals, is applied here to calculate the excitation energies for increasingly longer linear and cyclic acenes as part of their intense benchmarking for excited states of all types. The energies for the two lowest-lying singlet ¹L_a and ¹L_b states of linear oligoacenes as well as the triplet ³L_a and ³L_b states, are calculated and compared with experimental results. These functionals clearly outperform the results obtained from hybrid functionals and favorably compare with other double-hybrid expressions also tested here, such as B2-PLYP, B2GP-PLYP, ωB2-PLYP, and ωB2GP-PLYP. The study is complemented by the computation of adiabatic S₀-T₁ singlet-triplet energy difference for linear acenes as well as the extension of the study to strained cyclic oligomers, showing how the family of non-empirical expressions robustly leads to competitive results.

Inversion Theory Leveling as a New Methodological Approach to Antioxidant Thermodynamics: A Case Study on Phenol

Antioxidants are various types of compounds that represent a link between biology and chemistry. With the development of theoretical and computational methods, antioxidants are now being studied theoretically. Here, a novel method is presented that aims to reduce the estimated wall times for DFT calculations that result in the same or higher degree of accuracy in the second derivatives over energy than is the case with the regular computational route (i.e., optimizing the reaction system at a lower model and then recalculating the energies at a higher level of theory) by applying the inversion of theory level to the universal chemical scavenger model, i.e., phenol. The resulting accuracy and wall time obtained with such a methodological setup strongly suggest that this methodology could be generally applied to antioxidant thermodynamics for some costly DFT methods with relative absolute deviation.

Building on the strengths of a double-hybrid density functional for excitation energies and inverted singlet-triplet energy gaps

It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet-triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over "safe" excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet-triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature-the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S₁ and T₁ excited states.

A Reactive Molecular Dynamics Study of Chlorinated Organic Compounds. Part II: A ChemTraYzer Study of Chlorinated Dibenzofuran Formation and Decomposition Processes

In our two-paper series, we first present the development of ReaxFF CHOCl parameters using the recently published ParAMS parametrization tool. In this second part, we update the reactive Molecular Dynamics - Quantum Mechanics coupling scheme ChemTraYzer and combine it with our new ReaxFF parameters from Part I to study formation and decomposition processes of chlorinated dibenzofurans. We introduce a self-learning method for recovering failed transition-state searches that improves the overall ChemTraYzer transition-state search success rate by 10 percentage points to a total of 48 %. With ChemTraYzer, we automatically find and quantify more than 500 reactions using transition state theory and DFT. Among the discovered chlorinated dibenzofuran reactions are numerous reactions that are new to the literature. In three case studies, we discuss the set of reactions that are most relevant to the dibenzofuran literature: (i) bimolecular reactions of the chlorinated-dibenzofuran precursors phenoxy radical and 1,3,5-trichlorobenzene, (ii) dibenzofuran chlorination and pyrolysis, and (iii) oxidation of chlorinated dibenzofurans.

Achieving vibrational energies of diatomic systems with high quality by machine learning improved DFT method

When using ab initio methods to obtain high-quality quantum behavior of molecules, it often involves a lot of trial-and-error work in algorithm design and parameter selection, which requires enormous time and computational resource costs. In the study of vibrational energies of diatomic molecules, we found that starting from a low-precision DFT model and then correcting the errors using the high-dimensional function modeling capabilities of machine learning, one can considerably reduce the computational burden and improve the prediction accuracy. Data-driven machine learning is able to capture subtle physical information that is missing from DFT approaches. The results of ¹²C¹⁶O, ²⁴MgO and Na³⁵Cl show that, compared with CCSD(T)/cc-pV5Z calculation, this work improves the prediction accuracy by more than one order of magnitude, and reduces the computation cost by more than one order of magnitude.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Assessment of advanced xDH@B3LYP methods in describing various potential energy curves driven by π-π, CH/π, and SH/π non-bonded interactions

Accurate description of potential energy curves driven by non-bonded interactions remains a great challenge for pure density functional approximations (DFAs). It is because the R−6 decay behavior of dispersion cannot be intrinsically captured by the (semi)-local ingredients and the exact-exchange used in the popular hybrid DFAs. Overemphasizing the accuracy on the equilibrium region for the functional construction would likely deteriorate the overall performance on the other regions of potential energy surfaces. In consequence, the empirical dispersion correction becomes the standard component in DFAs to treat the non-bonded interactions. In this Letter, we demonstrate that without the use of empirical dispersion correction, doubly hybrid approximations, in particular two recently proposed revXYG3 and XYG7 functionals, hold the promise to have a balanced description of non-bonded interactions on the whole potential energy curves for several prototypes of π- π, CH/ π, and SH/ π interactions. The error of revXYG3 and XYG7 for non-bonded interactions is around 0.1 kcal/mol, and their potential energy curves almost coincide with the accurate CCSD(T)/CBS curves.

Pentacoordinate Carbon Atoms in a Ferrocene Dication Derivative—[Fe(Si2-η5-C5H2)2]2+

Pentacoordinate carbon atoms are theoretically predicted here in a ferrocene dication derivative in the eclipsed-(1; C2v), gauche-(2; C2) and staggered-[Fe(Si2-η5-C5H2)2]2+(3; C2h) forms for the first time. Energetically, the relative energy gaps for 2 and 3 range from −3.06 to 16.74 and −2.78 to 40.34 kJ mol−1, respectively, when compared to the singlet electronic state of 1 at different levels. The planar tetracoordinate carbon (ptC) atom in the ligand Si2C5H2 becomes a pentacoordinate carbon upon complexation. The ligand with a ptC atom was predicted to be both a thermodynamically and kinetically stable molecule by some of us in our earlier theoretical works. Natural bond orbital and adaptive natural density partitioning analyses confirm the pentacoordinate nature of carbon in these three complexes (1–3). Although they are hypothetical at the moment, they support the idea of “hypercoordinate metallocenes” within organometallic chemistry. Moreover, ab initio molecular dynamics simulations carried out at 298 K temperature for 2000 fs suggest that these molecules are kinetically stable.

Integral-Direct Hartree-Fock and Møller-Plesset Perturbation Theory for Periodic Systems with Density Fitting: Application to the Benzene Crystal

We present an algorithm and implementation of integral-direct, density-fitted Hartree-Fock (HF) and second-order Møller-Plesset perturbation theory (MP2) for periodic systems. The new code eliminates the formerly prohibitive storage requirements and allows us to study systems 1 order of magnitude larger than before at the periodic MP2 level. We demonstrate the significance of the development by studying the benzene crystal in both the thermodynamic limit and the complete basis set limit, for which we predict an MP2 cohesive energy of -72.8 kJ/mol, which is about 10-15 kJ/mol larger in magnitude than all previously reported MP2 calculations. Compared to the best theoretical estimate from literature, several modified MP2 models approach chemical accuracy in the predicted cohesive energy of the benzene crystal and hence may be promising cost-effective choices for future applications on molecular crystals.

Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights

To better understand the thermochemical kinetics and mechanism of a specific chemical reaction, an accurate estimation of barrier heights (forward and reverse) and reaction energies is vital. Because of the large size of reactants and transition state structures involved in real-life mechanistic studies (e.g., enzymatically catalyzed reactions), density functional theory remains the workhorse for such calculations. In this paper, we have assessed the performance of 91 density functionals for modeling the reaction energies and barrier heights on a large and chemically diverse data set (BH9) composed of 449 organic chemistry reactions. We have shown that range-separated hybrid functionals perform better than the global hybrids for BH9 barrier heights and reaction energies. Except for the PBE-based range-separated nonempirical double hybrids, range separation of the exchange term helps improve the performance for barrier heights and reaction energies. The 16-parameter Berkeley double hybrid, ωB97M(2), performs remarkably well for both properties. However, our minimally empirical range-separated double hybrid functionals offer marginally better accuracy than ωB97M(2) for BH9 barrier heights and reaction energies.

Unexpected steric hindrance failure in the gas phase F- + (CH3)3CI SN2 reaction

Base-induced elimination (E2) and bimolecular nucleophilic substitution (S_N2) reactions are of significant importance in physical organic chemistry. The textbook example of the retardation of S_N2 reactivity by bulky alkyl substitution is widely accepted based on the static analysis of molecular structure and steric environment. However, the direct dynamical evidence of the steric hindrance of S_N2 from experiment or theory remains rare. Here, we report an unprecedented full-dimensional (39-dimensional) machine learning-based potential energy surface for the 15-atom F⁻ + (CH₃)₃CI reaction, facilitating the reliable and efficient reaction dynamics simulations that can reproduce well the experimental outcomes and examine associated atomic-molecular level mechanisms. Moreover, we found surprisingly high “intrinsic” reactivity of S_N2 when the E2 pathway is completely blocked, indicating the reaction that intends to proceed via E2 transits to S_N2 instead, due to a shared pre-reaction minimum. This finding indicates that the competing factor of E2 but not the steric hindrance determines the small reactivity of S_N2 for the F⁻ + (CH₃)₃CI reaction. Our study provides new insight into the dynamical origin that determines the intrinsic reactivity in gas-phase organic chemistry.

Base-induced elimination (E2) and bimolecular nucleophilic substitution (SN2) are of significant importance in physical organic chemistry. Here, the authors show that the competing factor of E2 as opposed to steric hindrance determines the low reactivity of SN2 in the F⁻ + (CH₃)₃CI reaction.

Analytic Gradients for the Long-Range-Corrected XYG3 Type of Doubly Hybrid Density Functionals: Theory, Implementation, and Assessment

An analytic gradient algorithm for the long-range-corrected (lrc-) XYG3 type of doubly hybrid functionals has been derived and implemented. Due to the introduction of a long-range second-order perturbation (lrPT2) as the correction, the Z-vector equations need to be modified, where the construction of a unique total Lagrangian of the method is required. Geometry optimizations using lrc-XYG3 and lrc-XYGJ-OS have then been applied to the A21 data set that consists of 21 noncovalently interacting systems with CCSD(T) reference data for structures, in order to testify the idea of introducing lrPT2 for a better description of intermolecular interactions. While lrc-XYG3 was observed to offer a better description in dispersion-dominant systems as compared to the XYG3 functional, lrc-XYGJ-OS stands out with an overall better balanced performance in reproducing intermolecular geometries.

Dynamical Effects of SN2 Reactivity Suppression by Microsolvation: Dynamics Simulations of the F-(H2O) + CH3I Reaction on a 21-Dimensional Potential Energy Surface

A comparison of atomistic dynamics between microsolvated and unsolvated reactions can expose the precise role of solvent molecules and thus provide deep insight into how solvation influences chemical reactions. Here we developed the first full-dimensional analytical potential energy surface of the F-(H2O) + CH3I reaction, which facilitates the efficient dynamics simulations on a quantitatively accurate level. The computed SN2 reactivity suppression ratio of the monosolvated F-(H2O) + CH3I reaction relative to the unsolvated F- + CH3I reaction as a function of collision energy first increases and then decreases steadily, forming an inverted-V shape, due to the combined dynamical effects of interaction time, steric hindrance, and collision-induced dehydration. Moreover, further analysis reveals that the steric effect of the F-(H2O) + CH3I reaction resulting from the single water molecule is manifested mainly in dragging the F- anion away from the central C atom, rather than shielding F- from C. Our study shows there is great potential in rigorously studying the role of the solvent in more complicated reactions.

Accurate Prediction of Nuclear Magnetic Resonance Parameters via the XYG3 Type of Doubly Hybrid Density Functionals

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful and versatile tools in elucidating molecular structures. To eliminate ambiguities of experimental assignments, accurate calculations of NMR spectra are of great importance. Here, a method for theoretical evaluation of the NMR shielding constants by analytic derivatives using gauge including atomic orbitals (GIAO) has been implemented for the XYG3 type of doubly hybrid density functionals (xDH), namely, the GIAO-xDH method. Benchmark calculations on shielding constants and chemical shifts demonstrate the remarkable accuracy of the GIAO-xDH method, compared to the accurate CCSD(T) references. It is shown here that the XYGJ-OS functional is able to give a mean absolute deviation (MAD) of ∼3.0 ppm in the calculated shielding constants for ¹³C, ¹⁵N, ¹⁷O, ¹⁹F, while both XYGJ-OS and xDH-PBE0 functionals are able to provide a satisfactory estimation of chemical shifts with MADs of ∼0.03 and 1.0 ppm for ¹H and ¹³C, respectively. The basis set influence upon the method has been examined and a computational scheme considering both accuracy and efficiency has been proposed and tested to predict the experimental ¹³C chemical shifts of five medium-sized natural product molecules, yielding a MAD of ∼1.0 ppm, which demonstrates the practical feasibility of the GIAO-xDH method.

Improving Results by Improving Densities: Density-Corrected Density Functional Theory

Density functional theory (DFT) calculations have become widespread in both chemistry and materials, because they usually provide useful accuracy at much lower computational cost than wavefunction-based methods. All practical DFT calculations require an approximation to the unknown exchange-correlation energy, which is then used self-consistently in the Kohn-Sham scheme to produce an approximate energy from an approximate density. Density-corrected DFT is simply the study of the relative contributions to the total energy error. In the vast majority of DFT calculations, the error due to the approximate density is negligible. But with certain classes of functionals applied to certain classes of problems, the density error is sufficiently large as to contribute to the energy noticeably, and its removal leads to much better results. These problems include reaction barriers, torsional barriers involving π-conjugation, halogen bonds, radicals and anions, most stretched bonds, etc. In all such cases, use of a more accurate density significantly improves performance, and often the simple expedient of using the Hartree-Fock density is enough. This Perspective explains what DC-DFT is, where it is likely to improve results, and how DC-DFT can produce more accurate functionals. We also outline challenges and prospects for the field.

Cov_DOX: A Method for Structure Prediction of Covalent Protein-Ligand Bindings

A handful of molecular docking tools have been extended to enable a covalent docking. However, all of them face the challenge brought by the covalent bond between proteins and ligands. Many covalent drug design scenarios still heavily rely on demanding crystallographic experiments for accurate binding structures. Aiming at filling the gap between covalent dockings and crystallographic experiments, we develop and validate a hybrid method, dubbed as Cov_DOX, in this work. Cov_DOX achieves an overall success rate of 81% with RMSD < 2 Å for the Top 1 pose prediction in the validation against a test set including 405 crystal structures for covalent protein-ligand complexes, covering various types of the warhead chemistry and receptors. Such accuracy is not far from the much more demanding crystallographic experiments, in sharp contrast to the performance of the covalent docking front runners (success rate: 40-60%).

A general justification for hybrid functionals in DFT by means of linear response theory

In the present work, resorting to linear response theory, we examine the plausibility of postulating Kohn-Sham (KS)-type equations which contain, by definition, an effective hybrid potential made up by some arbitrary mixture of local and non-local terms. In this way a general justification for the construction of hybrid functionals is provided without resorting to arguments based on the adiabatic connection, the generalized KS theory or the Levy's constrained search (or its variations). In particular, we examine the cases of single-hybrid functionals, derived from non-local exchange and of double-hybrid functionals, emerging from non-local second-order expressions obtained from the KS perturbation theory. A further generalization for higher-order hybrid functionals is also included.

Improving Density Functional Prediction of Molecular Thermochemical Properties with a Machine-Learning-Corrected Generalized Gradient Approximation

The past decade has seen an increasing interest in designing sophisticated density functional approximations (DFAs) by integrating the power of machine learning (ML) techniques. However, application of the ML-based DFAs is often confined to simple model systems. In this work, we construct an ML correction to the widely used Perdew-Burke-Ernzerhof (PBE) functional by establishing a semilocal mapping from the electron density and reduced gradient to the exchange-correlation energy density. The resulting ML-corrected PBE is immediately applicable to any real molecule and yields significantly improved heats of formation while preserving the accuracy for other thermochemical and kinetic properties. This work highlights the prospect of combining the power of data-driven ML methods with physics-inspired derivations for reaching the heaven of chemical accuracy.

Accurate Spectral Properties within Double-Hybrid Density Functional Theory: A Spin-Scaled Range-Separated Second-Order Algebraic-Diagrammatic Construction-Based Approach

Our second-order algebraic-diagrammatic construction [ADC(2)]-based double-hybrid (DH) ansatz (J. Chem. Theory Comput. 2019, 15, 4440. DOI: 10.1021/acs.jctc.9b00391) is combined with range-separation techniques. In the present scheme, both the exchange and the correlation contributions are range-separated, while spin-scaling approaches are also applied. The new methods are thoroughly tested for the most popular benchmark sets including 250 singlet and 156 triplet excitations, as well as 80 oscillator strengths. It is demonstrated that the range separation for the correlation contributions is highly recommended for both the genuine and the ADC(2)-based DH approaches. Our results show that the latter scheme slightly but consistently outperforms the former one for single excitation dominated transitions. Furthermore, states with larger fractions of double excitations are assessed as well, and challenging charge-transfer excitations are also discussed, where the recently proposed spin-scaled long-range corrected DHs fail. The suggested iterative fourth-power scaling RS-PBE-P86/SOS-ADC(2) method, using only three adjustable parameters, provides the most robust and accurate excitation energies within the DH theory. In addition, the relative error of the oscillator strengths is reduced by 65% compared to the best genuine DH functionals.

Seeking the Optimal Descriptor for SN2 Reactions through Statistical Analysis of Density Functional Theory Results

Bimolecular nucleophilic substitution is one of the fundamental reactions in organic chemistry, yet there is still knowledge to be gained on the role of the nucleophile and the substrate. A statistical treatment of over 600 density functional theory (DFT)-computed barriers for bimolecular nucleophilic substitution at methyl derivatives (S_N2@C) leads to the identification of numerical descriptors that best represent the entering and leaving ability of 26 different nucleophiles. The treatment is based on singular value decomposition (SVD) of a matrix of computed energy barriers. The current work represents the extension to a problem of reactivity of the hidden descriptor methodology that we had previously developed for the thermodynamic problem of bond dissociation energies in transition-metal complexes. The analysis of the results shows that a single descriptor is sufficient. This hidden descriptor has different values for nucleophilic and leaving abilities and, contrary to expectation, does not correlate especially well with either frontier molecular orbital descriptors or solvation descriptors. In contrast, it correlates with other thermodynamic and geometric parameters. This statistical procedure can be in principle extended to additional chemical fragments and other reactions.

Symmetry Dilemma of Doubly Hybrid Density Functionals for Equilibrium Molecular Property Calculations

In electronic structure theory, when an approximate wavefunction tends to artifactually break the symmetry of the exact Hamiltonian, the corresponding method is referred to as having a "symmetry dilemma" problem. Such types of artifacts were often reported when Hartree-Fock (HF) and the low-level post-HF methods were used, while the traditional Kohn-Sham density functional theory (KS-DFT) methods were usually found to be more resistant to this breakdown. In this work, we present a systematic study on the reliability of the doubly hybrid (DH) DFT methods for several violable cases. Almost all the commonly used B2PLYP-type (bDH) functionals are shown to have a severe "symmetry dilemma" problem and yield dramatically unreliable molecular properties, such as dipole moment, vibrational frequency, and static polarizability at the equilibrium geometry. A one-parameter bDH functional model study demonstrates that such a problem is a combined effect of the inappropriate portion of the HF exchange (over 50%) for the self-consistent field (SCF) calculation and the augmentation of the second-order perturbative contribution. It is remarkable that the XYG3-type (xDH) functionals show a good capability to resist the artifactual symmetry breaking and yield reliable molecular properties when the same critical cases are calculated. In the xDH, there are two functionals of different purposes, namely, the SCF functional and the energy functional, which have different amounts of the HF exchange and different portions of the correlation contributions. The success of the xDHs can be attributed to this flexibility in xDH construction to avoid using an improperly large portion of the HF exchange in the SCF functional. The insights gained in this work are of significance for the development of an improved DH functional.

Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data

Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)^MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol^-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH₂, Ac-Ala-NH₂, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)^MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.

A full-dimensional ab initio potential energy and dipole moment surfaces for (NH3)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) for the ammonia dimer (NH₃)₂ are reported. The database of the PES consists of 27 736 ab initio energy points and all of these points were calculated at the UCCSD(T)-F12a/AVTZ level. The PES was fitted by using the fundamental invariant neural network (FI-NN) method that satisfies the permutational symmetry of identical atoms, and the root mean square fitting error for the PES is very small as low as 0.562 meV. The geometries for the (NH₃)₂ DMS are the same as those used for the PES and are calculated at the XYG3/AVTZ level. This PES can describe a variety of internal floppy motions, including all kinds of vibrational modes no matter intermolecular or intramolecular. The CCSD(T)-PES can dissociate correctly to two NH₃ monomers, with D_e = 1135.55 cm^-1 (13.58 kJ/mol) which agrees accurately with the 13.5 ± 0.3 kJ/mol predicted by previous work.