Generalized energy-based fragmentation approach and its applications to macromolecules and molecular aggregates

Automatic molecular fragmentation by evolutionary optimisation

Molecular fragmentation is an effective suite of approaches to reduce the formal computational complexity of quantum chemistry calculations while enhancing their algorithmic parallelisability. However, the practical applicability of fragmentation techniques remains hindered by a dearth of automation and effective metrics to assess the quality of a fragmentation scheme. In this article, we present the Quick Fragmentation via Automated Genetic Search (QFRAGS), a novel automated fragmentation algorithm that uses a genetic optimisation procedure to generate molecular fragments that yield low energy errors when adopted in Many Body Expansions (MBEs). Benchmark testing of QFRAGS on protein systems with less than 500 atoms, using two-body (MBE2) and three-body (MBE3) MBE calculations at the HF/6-31G* level, reveals mean absolute energy errors (MAEE) of 20.6 and 2.2 kJ  mol - 1 , respectively. For larger protein systems exceeding 500 atoms, MAEEs are 181.5 kJ  mol - 1 for MBE2 and 24.3 kJ  mol - 1 for MBE3. Furthermore, when compared to three manual fragmentation schemes on a 40-protein dataset, using both MBE and Fragment Molecular Orbital techniques, QFRAGS achieves comparable or often lower MAEEs. When applied to a 10-lipoglycan/glycolipid dataset, MAEs of 7.9 and 0.3 kJ  mol - 1 were observed at the MBE2 and MBE3 levels, respectively.Scientific Contribution This Article presents the Quick Fragmentation via Automated Genetic Search (QFRAGS), an innovative molecular fragmentation algorithm that significantly improves upon existing molecular fragmentation approaches by specifically addressing their lack of automation and effective fragmentation quality metrics. With an evolutionary optimisation strategy, QFRAGS actively pursues high quality fragments, generating fragmentation schemes that exhibit minimal energy errors on systems with hundreds to thousands of atoms. The advent of QFRAGS represents a significant advancement in molecular fragmentation, greatly improving the accessibility and computational feasibility of accurate quantum chemistry calculations.

Intermolecular interaction energies with AROFRAG-A systematic approach for fragmentation of aromatic molecules

Intermolecular interactions with polycyclic aromatic hydrocarbons (PAHs) represent an important area of physisorption studies. These investigations are often hampered by a size of interacting PAHs, which makes the calculation prohibitively expensive. Therefore, methods designed to deal with large molecules could be helpful to reduce the computational costs of such studies. Recently we have introduced a new systematic approach for the molecular fragmentation of PAHs, denoted as AROFRAG, which decomposes a large PAH molecule into a set of predefined small PAHs with a benzene ring being the smallest unbreakable unit, and which in conjunction with the Molecules-in-Molecules (MIM) approach provides an accurate description of total molecular energies. In this contribution we propose an extension of the AROFRAG, which provides a description of intermolecular interactions for complexes composed of PAH molecules. The examination of interaction energy partitioning for various test cases shows that the AROFRAG3 model connected with the MIM approach accurately reproduces all important components of the interaction energy. An additional important finding in our study is that the computationally expensive long-range electron-correlation part of the interaction energy, that is, the dispersion component, is well described at lower AROFRAG levels even without MIM, which makes the latter models interesting alternatives to existing methods for an accurate description of the electron-correlated part of the interaction energy.

The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

Fragment-Based Deep Learning for Simultaneous Prediction of Polarizabilities and NMR Shieldings of Macromolecules and Their Aggregates

Simultaneous prediction of the molecular response properties, such as polarizability and the NMR shielding constant, at a low computational cost is an unresolved issue. We propose to combine a linear-scaling generalized energy-based fragmentation (GEBF) method and deep learning (DL) with both molecular and atomic information-theoretic approach (ITA) quantities as effective descriptors. In GEBF, the total molecular polarizability can be assembled as a linear combination of the corresponding quantities calculated from a set of small embedded subsystems in GEBF. In the new GEBF-DL(ITA) protocol, one can predict subsystem polarizabilities based on the corresponding molecular wave function (thus electron density and ITA quantities) and DL model rather than calculate them from the computationally intensive coupled-perturbed Hartree-Fock or Kohn-Sham equations and finally obtain the total molecular polarizability via a linear combination equation. As a proof-of-concept application, we predict the molecular polarizabilities of large proteins and protein aggregates. GEBF-DL(ITA) is shown to be as accurate enough as GEBF, with mean absolute percentage error <1%. For the largest protein aggregate (>4000 atoms), GEBF-DL(ITA) gains a speedup ratio of 3 compared with GEBF. It is anticipated that when more advanced electronic structure methods are used, this advantage will be more appealing. Moreover, one can also predict the NMR chemical shieldings of proteins with reasonably good accuracy. Overall, the cost-efficient GEBF-DL(ITA) protocol should be a robust theoretical tool for simultaneously predicting polarizabilities and NMR shieldings of large systems.

AROFRAG─A Systematic Approach for Fragmentation of Aromatic Molecules

We present a new systematic fragmentation scheme of polycyclic aromatic hydrocarbons (PAHs), including fullerenes and nanotubes, based on an idea to treat a sextet ring as a single unbreakable unit so that the basic unit of aromaticity remains preserved upon fragmentation. In the approach, denoted as AROFRAG (from aromatic fragmentation), a set of predefined elementary subsystems, such as naphthalene and biphenyl in the first model and larger PAHs in the second and third models, is generated with appropriate weights with the aim of reproducing the structure of the original molecule. The energies of the molecules are approximated as weighted sums of the energies of these subsystems. For symmetric cases, e.g., fullerenes, the point-group symmetry is preserved during the decomposition. The AROFRAG is used in conjunction with the molecule-in-molecule (MIM) technique to obtain an accurate description of the electronic energies. The new approach has been applied for selected graphene structures and fullerene doped with boron and nitrogen atoms, for which isomerization energies were calculated, as well as for several nanotubes and regular fullerene molecules. The combination of the third AROFRAG model and the MIM approach leads to the reproduction of electronic energies with a few milli-hartree accuracy at a fraction of the computational cost of the original method.

Reformulation of All ONIOM-Type Molecular Fragmentation Approaches and Many-Body Theories Using Graph-Theory-Based Projection Operators: Applications to Dynamics, Molecular Potential Surfaces, Machine Learning, and Quantum Computing

We present a graph-theory-based reformulation of all ONIOM-based molecular fragmentation methods. We discuss applications to (a) accurate post-Hartree-Fock AIMD that can be conducted at DFT cost for medium-sized systems, (b) hybrid DFT condensed-phase studies at the cost of pure density functionals, (c) reduced cost on-the-fly large basis gas-phase AIMD and condensed-phase studies, (d) post-Hartree-Fock-level potential surfaces at DFT cost to obtain quantum nuclear effects, and (e) novel transfer machine learning protocols derived from these measures. Additionally, in previous work, the unifying strategy discussed here has been used to construct new quantum computing algorithms. Thus, we conclude that this reformulation is robust and accurate.

Capturing Weak Interactions in Surface Adsorbate Systems at Coupled Cluster Accuracy: A Graph-Theoretic Molecular Fragmentation Approach Improved through Machine Learning

The accurate and efficient study of the interactions of organic matter with the surface of water is critical to a wide range of applications. For example, environmental studies have found that acidic polyfluorinated alkyl substances, especially perfluorooctanoic acid (PFOA), have spread throughout the environment and bioaccumulate into human populations residing near contaminated watersheds, leading to many systemic maladies. Thus, the study of the interactions of PFOA with water surfaces became important for the mitigation of their activity as pollutants and threats to public health. However, theoretical study of the interactions of such organic adsorbates on the surface of water, and their bulk concerted properties, often necessitates the use of ab initio methods to properly incorporate the long-range electronic properties that govern these extended systems. Notable theoretical treatments of "on-water" reactions thus far have employed hybrid DFT and semilocal DFT, but the interactions involved are weak interactions that may be best described using post-Hartree-Fock theory. Here, we aim to demonstrate the utility of a graph-theoretic approach to molecular fragmentation that accurately captures the critical "weak" interactions while maintaining an efficient ab initio treatment of the long-range periodic interactions that underpin the physics of extended systems. We apply this graph-theoretical treatment to study PFOA on the surface of water as a model system for the study of weak interactions seen in the wide range of surface interactions and reactions. The approach divides a system into a set of vertices, that are then connected through edges, faces, and higher order graph theoretic objects known as simplexes, to represent a collection of locally interacting subsystems. These subsystems are then used to construct ab initio molecular dynamics simulations and for computing multidimensional potential energy surfaces. To further improve the computational efficiency of our graph theoretic fragmentation method, we use a recently developed transfer learning protocol to construct the full system potential energy from a family of neural networks each designed to accurately model the behavior of individual simplexes. We use a unique multidimensional clustering algorithm, based on the k-means clustering methodology, to define our training space for each separate simplex. These models are used to extrapolate the energies for molecular dynamics trajectories at PFOA water interfaces, at less than one-tenth the cost as compared to a regular molecular fragmentation-based dynamics calculation with excellent agreement with couple cluster level of full system potential energies.

Seamless Multilayer─A Novel Total Energy Partition Scheme for Embedded and Hybrid Calculations

In this paper, we provide general formulation of a multilayer approach, covering both additive and subtractive quantum mechanics/molecular mechanics (QM/MM) as special cases. After that, we suggest a novel definition of QM/MM total energy based on the consideration of a system divided into three layers. In a simplified form, it is E = E Q M ( 1 + 2 ) - E Q M ( 2 ) + E M M ( 2 + 3 ) , where layers 1, 2, and 3 represent inner QM, outer QM, and classical MM regions, respectively. The novel formulation is also not limited by only QM/MM combination of methods─in fact, any computational methods can be combined in a hybrid calculation. In this paper, we call the new approach seamless multilayer. Test calculations performed for silica and boric oxide show that the new approach requires no QM/MM interface parameterization as well as no or very simple correction terms for boundary atoms. This can greatly facilitate QM/MM studies of covalent inorganic solids. However, test calculations of α-Al2O3 show that for ionic compounds, the new method requires some additional development.

Graph-|Q⟩⟨C|: A Quantum Algorithm with Reduced Quantum Circuit Depth for Electronic Structure

The accurate determination of chemical properties is known to have a critical impact on multiple fundamental chemical problems but is deeply hindered by the steep algebraic scaling of electron correlation calculations and the exponential scaling of quantum nuclear dynamics. With the advent of new quantum computing hardware and associated developments in creating new paradigms for quantum software, this avenue has been recognized as perhaps one way to address exponentially complex challenges in quantum chemistry and molecular dynamics. In this paper, we discuss a new approach to drastically reduce the quantum circuit depth (by several orders of magnitude) and help improve the accuracy in the quantum computation of electron correlation energies for large molecular systems. The method is derived from a graph-theoretic approach to molecular fragmentation and enables us to create a family of projection operators that decompose quantum circuits into separate unitary processes. Some of these processes can be treated on quantum hardware and others on classical hardware in a completely asynchronous and parallel fashion. Numerical benchmarks are provided through the computation of unitary coupled-cluster singles and doubles (UCCSD) energies for medium-sized protonated and neutral water clusters using the new quantum algorithms presented here.

Competition between Solvent···Solvent and Solvent···Solute Interactions in the Microhydration of the Tetrafluoroborate Anion, BF4-(H2O)n=1,2,3,4

This study systematically examines the interactions of the tetrafluoroborate anion (BF4-) with up to four water molecules (BF4-(H2O)n=1,2,3,4). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using a variety of density functional theory (DFT) methods (B3LYP, B3LYP-D3BJ, and M06-2X) and the MP2 ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for B, O, and F; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the CCSD(T) ab initio method and the haTZ basis set for the mono- and dihydrate (n = 1, 2) structures. The 2-body:Many-body (2b:Mb) technique, in which CCSD(T) computations capture the 1- and 2-body contributions to the interactions and MP2 computations recover all higher-order contributions, was used to extend these demanding computations to the tri- and tetrahydrate (n = 3, 4) systems. Four, five, and eight new stationary points have been identified for the di-, tri-, and tetrahydrate systems, respectively. The global minimum of the monohydrate adopts a symmetric double ionic hydrogen bond motif with C2v symmetry and an electronic dissociation energy of 13.17 kcal mol-1 at the CCSD(T)/haTZ level of theory. This strong solvent···solute interaction, however, competes with solute···solute interactions in the lowest-energy BF4-(H2O)n=2,3,4 minima that are not seen in the other di-, tri-, or tetrahydrate minima. The latter interactions help increase the 2b:Mb dissociation energies to more than 26, 41, and 51 kcal mol-1 for n = 2, 3, and 4, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -280 cm-1 relative to an isolated H2O molecule at the 2b:Mb/haTZ level of theory.

Accurate and Efficient Prediction of Post-Hartree-Fock Polarizabilities of Condensed-Phase Systems

To accurately and efficiently predict the molecular response properties (such as polarizability) at post-Hartree-Fock levels for condensed-phase systems under periodic boundary conditions (PBC) is still an unaccomplished and ongoing task. We demonstrate that static isotropic polarizabilities can be cost-effectively predicted at post-Hartree-Fock levels by combining the linear-scaling generalized energy-based fragmentation (GEBF) and information-theoretic approach (ITA) quantities. In PBC-GEBF, the total molecular polarizability of an extended system is obtained as a linear combination of the corresponding quantities of a series of small embedded subsystems of several monomers. Here, we show that in the PBC-GEBF-ITA framework, one can obtain the molecular polarizabilities and establish linear relations to ITA quantities. Once these relations are established for smaller subsystems, one can predict the polarizabilities of larger subsystems directly from the molecular wavefunction (or electron density) via ITA quantities. Alternatively, one can determine the total molecular polarizability via a linear combination equation in PBC-GEBF. We have corroborated that this newly proposed PBC-GEBF-ITA protocol is much more efficient than the original PBC-GEBF approach but is not much less accurate and that this conclusion holds for both many-body perturbation theory and the coupled cluster calculations. Good efficiency and transferability of the PBC-GEBF-ITA protocol are demonstrated for periodic systems with several hundred atoms in a unit cell.

Toward an extreme-scale electronic structure system

Electronic structure calculations have the potential to predict key matter transformations for applications of strategic technological importance, from drug discovery to material science and catalysis. However, a predictive physicochemical characterization of these processes often requires accurate quantum chemical modeling of complex molecular systems with hundreds to thousands of atoms. Due to the computationally demanding nature of electronic structure calculations and the complexity of modern high-performance computing hardware, quantum chemistry software has historically failed to operate at such large molecular scales with accuracy and speed that are useful in practice. In this paper, novel algorithms and software are presented that enable extreme-scale quantum chemistry capabilities with particular emphasis on exascale calculations. This includes the development and application of the multi-Graphics Processing Unit (GPU) library LibCChem 2.0 as part of the General Atomic and Molecular Electronic Structure System package and of the standalone Extreme-scale Electronic Structure System (EXESS), designed from the ground up for scaling on thousands of GPUs to perform high-performance accurate quantum chemistry calculations at unprecedented speed and molecular scales. Among various results, we report that the EXESS implementation enables Hartree-Fock/cc-pVDZ plus RI-MP2/cc-pVDZ/cc-pVDZ-RIFIT calculations on an ionic liquid system with 623 016 electrons and 146 592 atoms in less than 45 min using 27 600 GPUs on the Summit supercomputer with a 94.6% parallel efficiency.

Approaching Coupled Cluster Accuracy with Density Functional Theory Using the Generalized Connectivity-Based Hierarchy

This Perspective reviews connectivity-based hierarchy (CBH), a systematic hierarchy of error-cancellation schemes developed in our group with the goal of achieving chemical accuracy using inexpensive computational techniques ("coupled cluster accuracy with DFT"). The hierarchy is a generalization of Pople's isodesmic bond separation scheme that is based only on the structure and connectivity and is applicable to any organic and biomolecule consisting of covalent bonds. It is formulated as a series of rungs involving increasing levels of error cancellation on progressively larger fragments of the parent molecule. The method and our implementation are discussed briefly. Examples are given for the applications of CBH involving (1) energies of complex organic rearrangement reactions, (2) bond energies of biofuel molecules, (3) redox potentials in solution, (4) pK_a predictions in the aqueous medium, and (5) theoretical thermochemistry combining CBH with machine learning. They clearly show that near-chemical accuracy (1-2 kcal/mol) is achieved for a variety of applications with DFT methods irrespective of the underlying density functional used. They demonstrate conclusively that seemingly disparate results, often seen with different density functionals in many chemical applications, are due to an accumulation of systematic errors in the smaller local molecular fragments that can be easily corrected with higher-level calculations on those small units. This enables the method to achieve the accuracy of the high level of theory (e.g., coupled cluster) while the cost remains that of DFT. The advantages and limitations of the method are discussed along with areas of ongoing developments.

Fragment-Based Calculations of Enzymatic Thermochemistry Require Dielectric Boundary Conditions

Electronic structure calculations on enzymes require hundreds of atoms to obtain converged results, but fragment-based approximations offer a cost-effective solution. We present calculations on enzyme models containing 500-600 atoms using the many-body expansion, comparing to benchmarks in which the entire enzyme-substrate complex is described at the same level of density functional theory. When the amino acid fragments contain ionic side chains, the many-body expansion oscillates under vacuum boundary conditions but rapid convergence is restored using low-dielectric boundary conditions. This implies that full-system calculations in the gas phase are inappropriate benchmarks for assessing errors in fragment-based approximations. A three-body protocol retains sub-kilocalorie per mole fidelity with respect to a supersystem calculation, as does a two-body calculation combined with a full-system correction at a low-cost level of theory. These protocols pave the way for application of high-level quantum chemistry to large systems via rigorous, ab initio treatment of many-body polarization.

Excited-State Polarizabilities: A Combined Density Functional Theory and Information-Theoretic Approach Study

Accurate and efficient determination of excited-state polarizabilities (α) is an open problem both experimentally and computationally. Following our previous work, (Phys. Chem. Chem. Phys. 2023, 25, 2131−2141), in which we employed simple ground-state (S0) density-related functions from the information-theoretic approach (ITA) to accurately and efficiently evaluate the macromolecular polarizabilities, in this work we aimed to predict the lowest excited-state (S1) polarizabilities. The philosophy is to use density-based functions to depict excited-state polarizabilities. As a proof-of-principle application, employing 2-(2′-hydroxyphenyl)benzimidazole (HBI), its substituents, and some other commonly used ESIPT (excited-state intramolecular proton transfer) fluorophores as model systems, we verified that either with S0 or S1 densities as an input, ITA quantities can be strongly correlated with the excited-state polarizabilities. When transition densities are considered, both S0 and S1 polarizabilities are in good relationships with some ITA quantities. The transferability of the linear regression model is further verified for a series of molecules with little or no similarity to those molecules in the training set. Furthermore, the excitation energies can be predicted based on multivariant linear regression equations of ITA quantities. This study also found that the nature of both the ground-state and excited-state polarizabilities of these species are due to the spatial delocalization of the electron density.

PESPIP: Software to fit complex molecular and many-body potential energy surfaces with permutationally invariant polynomials

We wish to describe a potential energy surface by using a basis of permutationally invariant polynomials whose coefficients will be determined by numerical regression so as to smoothly fit a dataset of electronic energies as well as, perhaps, gradients. The polynomials will be powers of transformed internuclear distances, usually either Morse variables, exp(-r_i,j/λ), where λ is a constant range hyperparameter, or reciprocals of the distances, 1/r_i,j. The question we address is how to create the most efficient basis, including (a) which polynomials to keep or discard, (b) how many polynomials will be needed, (c) how to make sure the polynomials correctly reproduce the zero interaction at a large distance, (d) how to ensure special symmetries, and (e) how to calculate gradients efficiently. This article discusses how these questions can be answered by using a set of programs to choose and manipulate the polynomials as well as to write efficient Fortran programs for the calculation of energies and gradients. A user-friendly interface for access to monomial symmetrization approach results is also described. The software for these programs is now publicly available.

Predicting the structures and vibrational spectra of molecular crystals containing large molecules with the generalized energy-based fragmentation approach

The generalized energy-based fragmentation (GEBF) approach under periodic boundary conditions (PBCs) has been developed to facilitate calculations of molecular crystals containing large molecules. The PBC-GEBF approach can help predict structures and properties of molecular crystals at different theory levels by performing molecular quantum chemistry calculations on a series of non-periodic subsystems constructed from the studied systems. A more rigorous formula of the forces on translational vectors of molecular crystals was proposed and implemented, enabling more reliable predictions of crystal structures. Our benchmark results on several typical molecular crystals show that the PBC-GEBF approach could reproduce the forces on atoms and the translational vectors and the optimized crystal structures from the corresponding conventional periodic methods. The improved PBC-GEBF approach is then applied to predict the crystal structures and vibrational spectra of two molecular crystals containing large molecules. The PBC-GEBF approach can provide a satisfactory description on the crystal structure of a molecular crystal containing 312 atoms in a unit cell at density-fitting second-order Møller-Plesset perturbation theory and density functional theory (DFT) levels and the infrared vibrational spectra of another molecular crystal containing 864 atoms in a unit cell at the DFT level. The PBC-GEBF approach is expected to be a promising theoretical tool for electronic structure calculations on molecular crystals containing large molecules.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

Efficient and accurate density-based prediction of macromolecular polarizabilities

Accurately and efficiently predicting macromolecules' polarizabilities is an open problem. In this work, we employ a few simple density-based quantities from the information-theoretic approach (ITA) to predict polarizability of proteins. We first build quantitative structure/property relationships between molecular polarizabilities and ITA quantities. We then verify the broad applicability of ITA quantities for polarizability prediction for inorganic, organic, and biological systems with both localized and delocalized electronic structure. As a proof-of-concept application, we predict the molecular polarizabilities of complex proteins. Based on the linear regression equations for 20 natural amino acid residues, 400 dipeptides, and 8000 tripeptides, one then predicts the molecular polarizability of a larger peptide or even a protein once the molecular wavefunction is obtained. Because it is extremely costly to determine the wavefunction for a macromolecule like a protein, we propose to combine the ITA with the linear-scaling generalized energy-based fragmentation (GEBF) method to predict the macromolecular polarizability. In GEBF, the total molecular polarizability is obtained as a linear combination of the corresponding quantities from a series of small subsystems. We can predict them based on the subsystem wavefunction and linear regression equations rather than compute them from the nearly-intractable coupled-perturbed Hartree-Fock or Kohn-Sham equations for the whole macromolecule. Computational results showcase that the GEBF-ITA protocol should be an inexpensive yet accurate theoretical tool for predicting macromolecular polarizabilities.

Treating Polarization Effects in Charged and Polar Bio-Molecules Through Variable Electrostatic Parameters

Polarization plays important roles in charged and hydrogen bonding containing systems. Much effort ranging from the construction of physics-based models to quantum mechanism (QM)-based and machine learning (ML)-assisted models have been devoted to incorporating the polarization effect into the conventional force fields at different levels, such as atomic and coarse grained (CG). The application of polarizable force fields or polarization models was limited by two aspects, namely, computational cost and transferability. Different from physics-based models, no predetermining parameters were required in the QM-based approaches. Taking advantage of both the accuracy of QM calculations and efficiency of molecular mechanism (MM) and ML, polarization effects could be treated more efficiently while maintaining the QM accuracy. The computational cost could be reduced with variable electrostatic parameters, such as the charge, dipole, and electronic dielectric constant with the help of linear scaling fragmentation-based QM calculations and ML models. Polarization and entropy effects on the prediction of partition coefficient of druglike molecules are demonstrated by using both explicit or implicit all-atom molecular dynamics simulations and machine learning-assisted models. Directions and challenges for future development are also envisioned.

Graph-Theoretic Molecular Fragmentation for Potential Surfaces Leads Naturally to a Tensor Network Form and Allows Accurate and Efficient Quantum Nuclear Dynamics

Molecular fragmentation methods have revolutionized quantum chemistry. Here, we use a graph-theoretically generated molecular fragmentation method, to obtain accurate and efficient representations for multidimensional potential energy surfaces and the quantum time-evolution operator, which plays a critical role in quantum chemical dynamics. In doing so, we find that the graph-theoretic fragmentation approach naturally reduces the potential portion of the time-evolution operator into a tensor network that contains a stream of coupled lower-dimensional propagation steps to potentially achieve quantum dynamics with reduced complexity. Furthermore, the fragmentation approach used here has previously been shown to allow accurate and efficient computation of post-Hartree-Fock electronic potential energy surfaces, which in many cases has been shown to be at density functional theory cost. Thus, by combining the advantages of molecular fragmentation with the tensor network formalism, the approach yields an on-the-fly quantum dynamics scheme where both the electronic potential calculation and nuclear propagation portion are enormously simplified through a single stroke. The method is demonstrated by computing approximations to the propagator and to potential surfaces for a set of coupled nuclear dimensions within a protonated water wire problem exhibiting the Grotthuss mechanism of proton transport. In all cases, our approach has been shown to reduce the complexity of representing the quantum propagator, and by extension action of the propagator on an initial wavepacket, by several orders, with minimal loss in accuracy.

Machine learning, artificial intelligence, and chemistry: how smart algorithms are reshaping simulation and the laboratory

Machine learning and artificial intelligence are increasingly gaining in prominence through image analysis, language processing, and automation, to name a few applications. Machine learning is also making profound changes in chemistry. From revisiting decades-old analytical techniques for the purpose of creating better calibration curves, to assisting and accelerating traditional in silico simulations, to automating entire scientific workflows, to being used as an approach to deduce underlying physics of unexplained chemical phenomena, machine learning and artificial intelligence are reshaping chemistry, accelerating scientific discovery, and yielding new insights. This review provides an overview of machine learning and artificial intelligence from a chemist’s perspective and focuses on a number of examples of the use of these approaches in computational chemistry and in the laboratory.

Graph Theoretic Molecular Fragmentation for Multidimensional Potential Energy Surfaces Yield an Adaptive and General Transfer Machine Learning Protocol

Over a series of publications we have introduced a graph-theoretic description for molecular fragmentation. Here, a system is divided into a set of nodes, or vertices, that are then connected through edges, faces, and higher-order simplexes to represent a collection of spatially overlapping and locally interacting subsystems. Each such subsystem is treated at two levels of electronic structure theory, and the result is used to construct many-body expansions that are then embedded within an ONIOM-scheme. These expansions converge rapidly with many-body order (or graphical rank) of subsystems and have been previously used for ab initio molecular dynamics (AIMD) calculations and for computing multidimensional potential energy surfaces. Specifically, in all these cases we have shown that CCSD and MP2 level AIMD trajectories and potential surfaces may be obtained at density functional theory cost. The approach has been demonstrated for gas-phase studies, for condensed phase electronic structure, and also for basis set extrapolation-based AIMD. Recently, this approach has also been used to derive new quantum-computing algorithms that enormously reduce the quantum circuit depth in a circuit-based computation of correlated electronic structure. In this publication, we introduce (a) a family of neural networks that act in parallel to represent, efficiently, the post-Hartree-Fock electronic structure energy contributions for all simplexes (fragments), and (b) a new k-means-based tessellation strategy to glean training data for high-dimensional molecular spaces and minimize the extent of training needed to construct this family of neural networks. The approach is particularly useful when coupled cluster accuracy is desired and when fragment sizes grow in order to capture nonlocal interactions accurately. The unique multidimensional k-means tessellation/clustering algorithm used to determine our training data for all fragments is shown to be extremely efficient and reduces the needed training to only 10% of data for all fragments to obtain accurate neural networks for each fragment. These fully connected dense neural networks are then used to extrapolate the potential energy surface for all molecular fragments, and these are then combined as per our graph-theoretic procedure to transfer the learning process to a full system energy for the entire AIMD trajectory at less than one-tenth the cost as compared to a regular fragmentation-based AIMD calculation.

Linear-Scaling Systematic Molecular Fragmentation Approach for Perturbation Theory and Coupled-Cluster Methods

The coupled-cluster (CC) singles and doubles with perturbative triples [CCSD(T)] method is frequently referred to as the “gold standard” of modern computational chemistry. However, the high computational cost of CCSD(T) [O(N⁷)], where N is the number of basis functions, limits its applications to small-sized chemical systems. To address this problem, efficient implementations of linear-scaling coupled-cluster methods, which employ the systematic molecular fragmentation (SMF) approach, are reported. In this study, we aim to do the following: (1) To achieve exact linear scaling and to obtain a pure ab initio approach, we revise the handling of nonbonded interactions in the SMF approach, denoted by LSSMF. (2) A new fragmentation algorithm, which yields smaller-sized fragments, that better fits high-level CC methods is introduced. (3) A modified nonbonded fragmentation scheme is proposed to enhance the existent algorithm. Performances of the LSSMF-CC approaches, such as LSSMF-CCSD(T), are compared with their canonical versions for a set of alkane molecules, C_nH_2n+2 (n = 6–10), which includes 142 molecules. Our results demonstrate that the LSSMF approach introduces negligible errors compared with the canonical methods; mean absolute errors (MAEs) are between 0.20 and 0.59 kcal mol^–1 for LSSMF(3,1)-CCSD(T). For a larger alkanes set (L12), C_nH_2n+2 (n = 50–70), the performance of LSSMF for the second-order perturbation theory (MP2) is investigated. For the L12 set, various bonded and nonbonded levels are considered. Our results demonstrate that the combination of bonded level 6 with nonbonded level 2, LSSMF(6,2), provides very accurate results for the MP2 method with a MAE value of 0.32 kcal mol^–1. The LSSMF(6,2) approach yields more than a 26-fold reduction in errors compared with LSSMF(3,1). Hence, we obtain substantial improvements over the original SMF approach. To illustrate the efficiency and applicability of the LSSMF-CCSD(T) approach, we consider an alkane molecule with 10,004 atoms. For this molecule, the LSSMF(3,1)-CCSD(T)/cc-pVTZ energy computation, on a Linux cluster with 100 nodes, 4 cores, and 5 GB of memory provided to each node, is performed just in ∼24 h. As a second test, we consider a biomolecular complex (PDB code: 1GLA), which includes 10,488 atoms, to assess the efficiency of the LSSMF approach. The LSSMF(3,1)-FNO–CCSD(T)/cc-pVTZ energy computation is completed in ∼7 days for the biomolecular complex. Hence, our results demonstrate that the LSSMF-CC approaches are very efficient. Overall, we conclude the following: (1) The LSSMF(m, n)-CCSD(T) methods can be reliably used for large-scale chemical systems, where the canonical methods are not computationally affordable. (2) The accuracy of bonded level 3 is not satisfactory for large chemical systems. (3) For high-accuracy studies, bonded level 5 (or higher) and nonbonded level 2 should be employed.

Combined fragment-based machine learning force field with classical force field and its application in the NMR calculations of macromolecules in solutions

We have developed a combined fragment-based machine learning (ML) force field and molecular mechanics (MM) force field for simulating the structures of macromolecules in solutions, and then compute its NMR chemical shifts with the generalized energy-based fragmentation (GEBF) approach at the level of density functional theory (DFT). In this work, we first construct Gaussian approximation potential based on GEBF subsystems of macromolecules for MD simulations and then a GEBF-based neural network (GEBF-NN) with deep potential model for the studied macromolecule. Then, we develop a GEBF-NN/MM force field for macromolecules in solutions by combining the GEBF-NN force field for the solute molecule and ff14SB force field for solvent molecules. Using the GEBF-NN/MM MD simulation to generate snapshot structures of solute/solvent clusters, we then perform the NMR calculations with the GEBF approach at the DFT level to calculate NMR chemical shifts of the solute molecule. Taking a heptamer of oligopyridine-dicarboxamides in chloroform solution as an example, our results show that the GEBF-NN force field is quite accurate for this heptamer by comparing with the reference DFT results. For this heptamer in chloroform solution, both the GEBF-NN/MM and classical MD simulations could lead to helical structures from the same initial extended structure. The GEBF-DFT NMR results indicate that the GEBF-NN/MM force field could lead to more accurate NMR chemical shifts on hydrogen atoms by comparing with the experimental NMR results. Therefore, the GEBF-NN/MM force field could be employed for predicting more accurate dynamical behaviors than the classical force field for complex systems in solutions.

A Density Functional Theory and Information-Theoretic Approach Study of Interaction Energy and Polarizability for Base Pairs and Peptides

Using density functional theory (DFT) and the information-theoretic approach (ITA) quantities to appreciate the energetics and properties of biopolymers is still an unaccomplished and ongoing task. To this end, we studied the building blocks of nucleic acid base pairs and small peptides. For base pairs, we have dissected the relative importance of energetic components by using two energy partition schemes in DFT. Our results convincingly show that the exchange-correlation effect predominantly governs the molecular stability of base pairs while the electrostatic potential plays a minor but indispensable role, and the steric effect is trivial. Furthermore, we have revealed that simple density-based ITA functions are in good relationships with molecular polarizabilities for a series of 30 hydrogen-bonded base pairs and all 20 natural α-amino acids, 400 dipeptides, and 8000 tripeptides. Based on these lines, one can easily predict the molecular polarizabilities of larger peptides, even proteins as long as the total molecular wavefunction is available, rather than solving the computationally demanding coupled-perturbed Hartree–Fock (CPHF) equation or its DFT counterpart coupled-perturbed Kohn–Sham (CPKS) equation.

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Structures and properties of ionic crystals and condensed phase ionic liquids predicted with the generalized energy-based fragmentation method

The generalized energy-based fragmentation (GEBF) approach is extended to facilitate ab initio investigations of structures, lattice energies, vibrational spectra and ¹ H NMR chemical shifts of ionic crystals and condensed-phase ionic liquids (ILs) with the periodic boundary conditions (PBC). For selected periodic systems, our results demonstrate that the so-called PBC-GEBF approach can provide satisfactory descriptions on ground-state energies, structures, and vibrational spectra of ionic crystals and IL crystals. The PBC-GEBF approach is then applied to three realistic condensed phase systems. For three ionic crystals (LiCl, NaCl, and KCl), we apply the PBC-GEBF approach with MP2 theory as well as some popular DFT methods to investigate their crystal structures and lattice energies. Our calculations indicate that the crystal structures obtained with PBC-GEBF-MP2/6-311 + G** are very close to the corresponding X-ray structures, while PBC-GEBF-ωB97X-D/6-311 + G** provides satisfactory prediction for crystal structures and lattice energies. For two polymorphs of [n-C₄ mim][Cl] crystals, we find that the PBC-GEBF approach at the M06-2X/6-311 + G** level can give a satisfactory descriptions on structures and Raman spectra of these two crystals. Furthermore, for [C₂ mim][BF₄ ] ILs, we demonstrate that their ¹ H NMR chemical shifts can be estimated from averaging over 5 typical snapshots (extracted from MD simulations) with the PBC-GEBF approach at the B97-2/pcSseg-2 level. The calculated results account for the observed experimental data quite well. Therefore, we expect that the PBC-GEBF approach, combined with various quantum chemistry methods, will become an effective tool in predicting structures and properties of ionic crystals and condensed-phase ILs.

A simple fragment-based method for van der Waals corrections over density functional theory

Modeling intermolecular noncovalent interactions between large molecules remains a challenge for the electron structure theory community due to the high cost. Fragment-based methods usually fare well in reducing the cost...

Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning

Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.

Graph-Theory-Based Molecular Fragmentation for Efficient and Accurate Potential Surface Calculations in Multiple Dimensions

We present a multitopology molecular fragmentation approach, based on graph theory, to calculate multidimensional potential energy surfaces in agreement with post-Hartree-Fock levels of theory but at the density functional theory cost. A molecular assembly is coarse-grained into a set of graph-theoretic nodes that are then connected with edges to represent a collection of locally interacting subsystems up to an arbitrary order. Each of the subsystems is treated at two levels of electronic structure theory, the result being used to construct many-body expansions that are embedded within an ONIOM scheme. These expansions converge rapidly with the many-body order (or graphical rank) of subsystems and capture many-body interactions accurately and efficiently. However, multiple graphs, and hence multiple fragmentation topologies, may be defined in molecular configuration space that may arise during conformational sampling or from reactive, bond breaking and bond formation, events. Obtaining the resultant potential surfaces is an exponential scaling proposition, given the number of electronic structure computations needed. We utilize a family of graph-theoretic representations within a variational scheme to obtain multidimensional potential surfaces at a reduced cost. The fast convergence of the graph-theoretic expansion with increasing order of many-body interactions alleviates the exponential scaling cost for computing potential surfaces, with the need to only use molecular fragments that contain a fewer number of quantum nuclear degrees of freedom compared to the full system. This is because the dimensionality of the conformational space sampled by the fragment subsystems is much smaller than the full molecular configurational space. Additionally, we also introduce a multidimensional clustering algorithm, based on physically defined criteria, to reduce the number of energy calculations by orders of magnitude. The molecular systems benchmarked include coupled proton motion in protonated water wires. The potential energy surfaces and multidimensional nuclear eigenstates obtained are shown to be in very good agreement with those from explicit post-Hartree-Fock calculations that become prohibitive as the number of quantum nuclear dimensions grows. The developments here provide a rigorous and efficient alternative to this important chemical physics problem.

Generalized energy-based fragmentation approach for calculations of solvation energies of large systems

A generalized energy-based fragmentation (GEBF) approach has been combined with a universal solvation model based on solute electron density (SMD) to compute the solvation energies of general large systems (such as protein molecules) in solutions. In the GEBF-SMD method, the solvation energy of a target system could be combined by the corresponding solvation energies of various subsystems, each of which is embedded in the background point charges and surface charges on the surface of solute cavity at the positions of its atoms and neighbouring atoms outside of the subsystem. Our results show that the GEBF-SMD model could reproduce the conventional SMD solvation energies quite well for various proteins in solutions, and could significantly reduce the computational costs for the SMD calculations of large proteins. In addition, the GEBF-SMD approach is almost independent of the basis sets and the types of solvents (including protic, polar, and nonpolar ones). Also, the GEBF-SMD approach could reproduce the relative energies of various conformers of large systems in solutions. Therefore, the GEBF-SMD method is expected to be applicable for computing the solvation energies of a broad range of large systems.

On the Accuracy of QM/MM Models: A Systematic Study of Intramolecular Proton Transfer Reactions of Amino Acids in Water

This work presents a systematic assessment of QM/QM' and QM/MM models with respect to direct QM calculations for the tautomerization (neutral to zwitterion) reactions of amino acids (glycine, alanine, valine, aspartate, and neutral and protonated histidine) solvated in a 160 water cluster. The effect of varying QM region size and choice of embedding potentials, including fixed-charge and polarizable molecular mechanics force fields (TIP3P and EFP) and various semiempirical QM methods (PM7, GFN2-xTB, DFTBA, DFTB3, HF-3c, and PBEh-3c), on the accuracy of the models was examined. A surprising finding was that molecular mechanics force fields outperformed many of the semiempirical methods. Generally, the errors in the QM/QM' and QM/MM models converge slowly with respect to the QM region size, requiring 50 or more waters to be included in the QM region before the error in the model falls below 1 kcal mol^-1 of its pure QM result. Different QM region selection schemes were also compared, and it was found that selection based on Natural Population Analysis (NPA) atomic charges significantly reduced the error in the QM/QM' and QM/MM models particularly if a low-quality embedding potential was used. It is envisaged that these results will be useful for the development of future hybrid QM models.

iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem

An iterative orbital interaction (iOI) approach is proposed to solve, in a bottom-up fashion, the self-consistent field problem in quantum chemistry. While it belongs grossly to the family of fragment-based quantum chemical methods, iOI is distinctive in that (1) it divides and conquers not only the energy but also the wave function and that (2) the subsystem sizes are automatically determined by successively merging neighboring small subsystems until they are just enough for converging the wave function to a given accuracy. Orthonormal occupied and virtual localized molecular orbitals are obtained in a natural manner, which can be used for all post-SCF purposes.

Density-Based Many-Body Expansion as an Efficient and Accurate Quantum-Chemical Fragmentation Method: Application to Water Clusters

Fragmentation methods based on the many-body expansion offer an attractive approach for the quantum-chemical treatment of large molecular systems, such as molecular clusters and crystals. Conventionally, the many-body expansion is performed for the total energy, but such an energy-based many-body expansion often suffers from a slow convergence with respect to the expansion order. For systems that show strong polarization effects such as water clusters, this can render the energy-based many-body expansion infeasible. Here, we establish a density-based many-body expansion as a promising alternative approach. By performing the many-body expansion for the electron density instead of the total energy and inserting the resulting total electron density into the total energy functional of density functional theory, one can derive a density-based energy correction, which in principle accounts for all higher-order polarization effects. Here, we systematically assess the accuracy of such a density-based many-body expansion for test sets of water clusters. We show that already a density-based two-body expansion is able to reproduce interaction energies per fragment within chemical accuracy and is able to accurately predict the energetic ordering as well as the relative interaction energies of different isomers of water clusters.

Multiscale Quantum Refinement Approaches for Metalloproteins

Biomolecules with metal ion(s) (e.g., metalloproteins) play many important biological roles. However, accurate structural determination of metalloproteins, particularly those containing transition metal ion(s), is challenging due to their complicated electronic structure, complex bonding of metal ions, and high number of conformations in biomolecules. Quantum refinement, which was proposed to combine crystallographic data with computational chemistry methods by several groups, can improve the local structures of some proteins. In this study, a quantum refinement method combining several multiscale computational schemes with experimental (X-ray diffraction) information was developed for metalloproteins. Various quantum refinement approaches using different ONIOM (our own N-layered integrated molecular orbital and molecular mechanics) combinations of quantum mechanics (QM), semiempirical (SE), and molecular mechanics (MM) methods were conducted to assess the performance and reliability on the refined local structure in two metalloproteins. The structures for two (Cu- or Zn-containing) metalloproteins were refined by combining two-layer ONIOM2(QM1/QM2) and ONIOM2(QM/MM) and three-layer ONIOM3(QM1/QM2/MM) schemes with experimental data. The accuracy of the quantum-refined metal binding sites was also examined and compared in these multiscale quantum refinement calculations. ONIOM3(QM/SE/MM) schemes were found to give good results with lower computational costs and were proposed to be a good choice for the multiscale computational scheme for quantum refinement calculations of metal binding site(s) in metalloproteins with high efficiency. Additionally, a two-center ONIOM approach was employed to speed up the quantum refinement calculations for the Zn metalloprotein with two remote active sites/ligands. Moreover, a recent quantum-embedding wavefunction-in-density functional theory (WF-in-DFT) method was also adopted as the high-level method in unprecedented ONIOM2(CCSD-in-B3LYP/MM) and ONIOM3(CCSD-in-B3LYP/SE/MM) calculations, which can be regarded as novel pseudo-three- and pseudo-four-layer ONIOM methods, respectively, to refine the key Zn binding site at the coupled-cluster singles and doubles (CCSD) level. These refined results indicate that multiscale quantum refinement schemes can be used to improve the structural accuracy obtained for local metal binding site(s) in metalloproteins with high efficiency.

Fragment-Based Ab Initio Molecular Dynamics Simulation for Combustion

We develop a fragment-based ab initio molecular dynamics (FB-AIMD) method for efficient dynamics simulation of the combustion process. In this method, the intermolecular interactions are treated by a fragment-based many-body expansion in which three- or higher body interactions are neglected, while two-body interactions are computed if the distance between the two fragments is smaller than a cutoff value. The accuracy of the method was verified by comparing FB-AIMD calculated energies and atomic forces of several different systems with those obtained by standard full system quantum calculations. The computational cost of the FB-AIMD method scales linearly with the size of the system, and the calculation is easily parallelizable. The method is applied to methane combustion as a benchmark. Detailed reaction network of methane reaction is analyzed, and important reaction species are tracked in real time. The current result of methane simulation is in excellent agreement with known experimental findings and with prior theoretical studies.

Machine learning builds full-QM precision protein force fields in seconds

Full-quantum mechanics (QM) calculations are extraordinarily precise but difficult to apply to large systems, such as biomolecules. Motivated by the massive demand for efficient calculations for large systems at the full-QM level and by the significant advances in machine learning, we have designed a neural network-based two-body molecular fractionation with conjugate caps (NN-TMFCC) approach to accelerate the energy and atomic force calculations of proteins. The results show very high precision for the proposed NN potential energy surface models of residue-based fragments, with energy root-mean-squared errors (RMSEs) less than 1.0 kcal/mol and force RMSEs less than 1.3 kcal/mol/Å for both training and testing sets. The proposed NN-TMFCC method calculates the energies and atomic forces of 15 representative proteins with full-QM precision in 10-100 s, which is thousands of times faster than the full-QM calculations. The computational complexity of the NN-TMFCC method is independent of the protein size and only depends on the number of residue species, which makes this method particularly suitable for rapid prediction of large systems with tens of thousands or even hundreds of thousands of times acceleration. This highly precise and efficient NN-TMFCC approach exhibits considerable potential for performing energy and force calculations, structure predictions and molecular dynamics simulations of proteins with full-QM precision.

Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms

We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.

Energy-based automatic determination of buffer region in the divide-and-conquer second-order Møller-Plesset perturbation theory

In the linear‐scaling divide‐and‐conquer (DC) electronic structure method, each subsystem is calculated together with the neighboring buffer region, the size of which affects the energy error introduced by the fragmentation in the DC method. The DC self‐consistent field calculation utilizes a scheme to automatically determine the appropriate buffer region that is as compact as possible for reducing the computational time while maintaining acceptable accuracy (J. Comput. Chem. 2018, 39, 909). To extend the automatic determination scheme of the buffer region to the DC second‐order Møller–Plesset perturbation (MP2) calculation, a scheme for estimating the subsystem MP2 correlation energy contribution from each atom in the buffer region is proposed. The estimation is based on the atomic orbital Laplace MP2 formalism. Based on this, an automatic buffer determination scheme for the DC‐MP2 calculation is constructed and its performance for several types of systems is assessed.

Systematic Partitioning of Proteins for Quantum-Chemical Fragmentation Methods Using Graph Algorithms

Quantum-chemical fragmentation methods offer an efficient approach for the treatment of large proteins, in particular if local target quantities such as protein-ligand interaction energies, enzymatic reaction energies, or spectroscopic properties of embedded chromophores are sought. However, the accuracy that is achievable for such local target quantities intricately depends on how the protein is partitioned into smaller fragments. While the commonly employed naı̈ve approach of using fragments with a fixed size is widely used, it can result in large and unpredictable errors when varying the fragment size. Here, we present a systematic partitioning scheme that aims at minimizing the fragmentation error of a local target quantity for a given maximum fragment size. To this end, we construct a weighted graph representation of the protein, in which the amino acids constitute the nodes. These nodes are connected by edges weighted with an estimate for the fragmentation error that is expected when cutting this edge. This allows us to employ graph partitioning algorithms provided by computer science to determine near-optimal partitions of the protein. We apply this scheme to a test set of six proteins representing various prototypical applications of quantum-chemical fragmentation methods using a simplified molecular fractionation with conjugate caps (MFCC) approach with hydrogen caps. We show that our graph-based scheme consistently improves upon the naı̈ve approach.

Transferable Multilevel Attention Neural Network for Accurate Prediction of Quantum Chemistry Properties via Multitask Learning

The development of efficient models for predicting specific properties through machine learning is of great importance for the innovation of chemistry and material science. However, predicting global electronic structure properties like Frontier molecular orbital highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy levels and their HOMO-LUMO gaps from the small-sized molecule data to larger molecules remains a challenge. Here, we develop a multilevel attention neural network, named DeepMoleNet, to enable chemical interpretable insights being fused into multitask learning through (1) weighting contributions from various atoms and (2) taking the atom-centered symmetry functions (ACSFs) as the teacher descriptor. The efficient prediction of 12 properties including dipole moment, HOMO, and Gibbs free energy within chemical accuracy is achieved by using multiple benchmarks, both at the equilibrium and nonequilibrium geometries, including up to 110,000 records of data in QM9, 400,000 records in MD17, and 280,000 records in ANI-1ccx for random split evaluation. The good transferability for predicting larger molecules outside the training set is demonstrated in both equilibrium QM9 and Alchemy data sets at the density functional theory (DFT) level. Additional tests on nonequilibrium molecular conformations from DFT-based MD17 data set and ANI-1ccx data set with coupled cluster accuracy as well as the public test sets of singlet fission molecules, biomolecules, long oligomers, and protein with up to 140 atoms show reasonable predictions for thermodynamics and electronic structure properties. The proposed multilevel attention neural network is applicable to high-throughput screening of numerous chemical species in both equilibrium and nonequilibrium molecular spaces to accelerate rational designs of drug-like molecules, material candidates, and chemical reactions.