An Efficient Fragment-Based Approach for Predicting the Ground-State Energies and Structures of Large Molecules

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.

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.

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.

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.

Graph-based quantum response theory and shadow Born-Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born-Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. The proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

Constructing Potential Energy Surface with Correlated Theory for Dipeptides Using Molecular Tailoring Approach

We demonstrate a cost-effective alternative employing the fragment-based molecular tailoring approach (MTA) for building the potential energy surface (PES) for two dipeptides viz. alanine-alanine and alanine-proline employing correlated theory, with augmented Dunning basis sets. About 1369 geometries are generated for each test dipeptide by systematically varying the dihedral angles Φ ${{\rm{\Phi }}}$ and Ψ ${{{\Psi }}}$ . These conformational geometries are partially optimized by relaxing all the other Z-matrix parameters, fixing the values of Φ ${{\rm{\Phi }}}$ and Ψ ${{{\Psi }}}$ . The MP2 level PES is constructed from the MTA-energies of chemically intact geometries using minimal hardware. The fidelity of MP2/aug-cc-pVDZ level PES is brought out by comparing it with its full calculation counterpart. Further, we bring out the power of the method by reporting the MTA-based CCSD/aug-cc-pVDZ level PES for these two dipeptides containing 498 and 562 basis functions respectively.

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.

A deep transfer learning-based protocol accelerates full quantum mechanics calculation of protein

Effective full quantum mechanics (FQM) calculation of protein remains a grand challenge and of great interest in computational biology with substantial applications in drug discovery, protein dynamic simulation and protein folding. However, the huge computational complexity of the existing QM methods impends their applications in large systems. Here, we design a transfer-learning-based deep learning (TDL) protocol for effective FQM calculations (TDL-FQM) on proteins. By incorporating a transfer-learning algorithm into deep neural network (DNN), the TDL-FQM protocol is capable of performing calculations at any given accuracy using models trained from small datasets with high-precision and knowledge learned from large amount of low-level calculations. The high-level double-hybrid DFT functional and high-level quality of basis set is used in this work as a case study to evaluate the performance of TDL-FQM, where the selected 15 proteins are predicted to have a mean absolute error of 0.01 kcal/mol/atom for potential energy and an average root mean square error of 1.47 kcal/mol/$ {\rm A^{^{ \!\!\!o}}} $ for atomic forces. The proposed TDL-FQM approach accelerates the FQM calculation more than thirty thousand times faster in average and presents more significant benefits in efficiency as the size of protein increases. The ability to learn knowledge from one task to solve related problems demonstrates that the proposed TDL-FQM overcomes the limitation of standard DNN and has a strong power to predict proteins with high precision, which solves the challenge of high precision prediction in large chemical and biological systems.

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.

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.

Implementation of Occupied and Virtual Edmiston-Ruedenberg Orbitals Using Cholesky Decomposed Integrals

We present a trust-region optimization of the Edmiston–Ruedenberg orbital localization function. The approach is used to localize both the occupied and the virtual orbitals and is the first demonstration of general virtual orbital localization using the Edmiston–Ruedenberg localization function. In the Edmiston–Ruedenberg approach, the sum of the orbital self-repulsion energies is maximized to obtain the localized orbitals. The Cholesky decomposition reduces the cost of transforming the electron repulsion integrals, and the overall scaling of our implementation is . The optimization is performed with all quantities in the molecular orbital basis, and the localization of the occupied orbitals is often less expensive than the corresponding self-consistent field (SCF) optimization. Furthermore, the occupied orbital localization scales linearly with the basis set. For the virtual space, the cost is significantly higher than the SCF optimization. The orbital spreads of the resulting virtual Edmiston–Ruedenberg orbitals are larger than for other, less expensive, orbital localization functions. This indicates that other localization procedures are more suitable for applications such as local post-Hartree–Fock calculations.

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.

All-Atom Nonadiabatic Dynamics Simulation of Hybrid Graphene Nanoribbons Based on Wannier Analysis and Machine Learning

Trajectory surface hopping combined with ab initio electronic structure calculations is a popular and powerful approach for on-the-fly nonadiabatic dynamics simulations. For large systems, however, this remains a significant challenge because of the unaffordable computational cost of large-scale electronic structure calculations. Here, we present an efficient divide-and-conquer approach to construct the system Hamiltonian based on Wannier analysis and machine learning. In detail, the large system under investigation is first decomposed into small building blocks, and then all possible segments formed by building blocks within a cutoff distance are found out. Ab initio molecular dynamics is carried out to generate a sequence of geometries for each equivalent segment with periodicity. The Hamiltonian matrices in the maximum localized Wannier function (MLWF) basis are obtained for all geometries and utilized to train artificial neural networks (ANNs) for the structure-dependent Hamiltonian elements. Taking advantage of the orthogonality and spatial locality of MLWFs, the one-electron Hamiltonian of a large system at arbitrary geometry can be directly constructed by the trained ANNs. As demonstrations, we study charge transport in a zigzag graphene nanoribbon (GNR), a coved GNR, and a series of hybrid GNRs with a state-of-the-art surface hopping method. The interplay between delocalized and localized states is found to determine the electron dynamics in hybrid GNRs. Our approach has successfully studied GNRs with >10 000 atoms, paving the way for efficient and reliable all-atom nonadiabatic dynamics simulation of general systems.

Terahertz Spectral Properties of 5-Substituted Uracils

Applications of terahertz time-domain spectroscopy (THz-TDS) in the fields of chemistry and biomedicine have recently received increased attention. Specifically, THz-TDS is particularly effective for the identification of alkaloid molecules, because it can distinguish the vibration types of base molecules in the THz band and provide a direct characteristic spectrum for the configuration and conformation of biomolecules. However, when THz-TDS technology is used to identify alkaloid molecules, most of them are concentrated in the 0.1-3.0 THz band, limiting the amount of information that can be obtained. In this work, a wide-spectrum THz-TDS system was independently built to explore the absorption spectra of uracil and its 5-substituents in the range of 1.3-6.0 THz. We found that, in the THz band, uracil and its 5-substituents have similar absorption peaks near 4.9 and 3.3 THz, while the 5-substituents have an additional absorption peak in the range of 1.5-2.5 THz. This absorption peak is red-shifted as the relative atomic mass of the 5-substituted atoms increases. Gaussian software was adopted to calculate the absorption spectra of the samples. The simulation conclusions were in good agreement with the experimental results, revealing that the vibration of the base molecule at low frequencies can be attributed to the inter-molecular vibration. This work demonstrates that THz-TDS technology can be used to accurately identify biomolecules with similar molecular structures, reflecting the importance of molecular structure in biological activity.

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.

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.

Three-Layer Multiscale Approach Based on Extremely Localized Molecular Orbitals to Investigate Enzyme Reactions

Quantum mechanics/molecular mechanics (QM/MM) calculations are widely used embedding techniques to computationally investigate enzyme reactions. In most QM/MM computations, the quantum mechanical region is treated through density functional theory (DFT), which offers the best compromise between chemical accuracy and computational cost. Nevertheless, to obtain more accurate results, one should resort to wave function-based methods, which however lead to a much larger computational cost already for relatively small QM subsystems. To overcome this drawback, we propose the coupling of our QM/ELMO (quantum mechanics/extremely localized molecular orbital) approach with molecular mechanics, thus introducing the three-layer QM/ELMO/MM technique. The QM/ELMO strategy is an embedding method in which the chemically relevant part of the system is treated at the quantum mechanical level, while the rest is described through frozen ELMOs. Since the QM/ELMO method reproduces results of fully QM computations within chemical accuracy and with a much lower computational effort, it can be considered a suitable strategy to extend the range of applicability and accuracy of the QM/MM scheme. In this paper, other than briefly presenting the theoretical bases of the QM/ELMO/MM technique, we will also discuss its validation on the well-tested deprotonation of acetyl coenzyme A by aspartate in citrate synthase.

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.

Structures and Spectroscopic Properties of Large Molecules and Condensed-Phase Systems Predicted by Generalized Energy-Based Fragmentation Approach

ConspectusThe structures and spectroscopic properties of molecules and condensed-phase systems are usually experimentally characterized by X-ray, infrared (IR), Raman, nuclear magnetic resonance (NMR), and electronic absorption/emission spectra. Quantum mechanics (QM) calculations are critical in quantitatively understanding the relationship between the structure and physicochemical properties of various chemical systems. However, it is very challenging to apply traditional QM methods to large molecules and condensed-phase systems with large unit cells due to their steep computational scaling with the system size. To overcome this difficulty, theoretical chemists have developed various linear (or low) scaling QM methods, among which energy-based fragmentation methods have achieved great success for large molecules or clusters. One of the most popular energy-based fragmentation methods is the generalized energy-based fragmentation (GEBF) approach developed by us.In this approach, the ground-state energy of a large molecule can be evaluated from the ground-state energies of a series of embedded subsystems. In this Account, we focus on the recent developments and applicability of the GEBF approach for the structures and spectroscopic properties of complicated large molecules and condensed-phase systems. With new fragmentation schemes, the GEBF approach can now describe ionic liquid clusters and metal-containing supramolecular systems accurately and can provide accurate binding energies for host-guest complexes. In addition, the GEBF approach is now available for describing the localized excited states of large systems including a chromophore. More importantly, the GEBF approach under periodic boundary conditions (PBC-GEBF) has been developed to deal with periodic molecular crystals and liquids. Then, the ground-state energy (or property) per unit cell of a periodic condensed phase system can be predicted with QM calculations on nonperiodic embedded subsystems. This feature enables accurate electron correlation calculations on molecular crystals and liquids to be feasible on ordinary workstations. The PBC-GEBF approach has been applied to predict the crystal structures, lattice energies, and spectroscopic properties of some typical molecular crystals and solutions. By combining the GEBF method and machine learning (ML) method, a GEBF-ML force field has been developed for long normal alkanes, and the IR spectra of long alkanes can be obtained from the GEBF-ML molecular dynamics (MD) simulations. The GEBF and its periodic variant are expected to play increasingly important roles in investigating real-life chemical systems of broad interests at the ab initio levels.

Adaptive density-guided approach to double incremental potential energy surface construction

We present a combination of the recently developed double incremental expansion of potential energy surfaces with the well-established adaptive density-guided approach to grid construction. This unique methodology is based on the use of an incremental expansion for potential energy surfaces, known as n-mode expansion; an incremental many-body representation of the electronic energy; and an efficient vibrational density-guided approach to automated determination of grid dimensions and granularity. The reliability of the method is validated calculating potential energy surfaces and obtaining fundamental excitation energies for three moderate-size chain-like molecular systems. The use of our methodology leads to considerable computational savings for potential energy surface construction compared to standard approaches while maintaining a high level of accuracy in the resulting potential energy surfaces. Additional investigations indicate that our method can be applied to covalently bound and strongly interacting molecular systems, even though these cases are known to be very unfavorable for fragmentation schemes. We therefore conclude that the presented methodology is a robust and flexible approach to potential energy surface construction, which introduces considerable computational savings without compromising the accuracy of vibrational spectra calculations.

eT 1.0: An open source electronic structure program with emphasis on coupled cluster and multilevel methods

The e^T program is an open source electronic structure package with emphasis on coupled cluster and multilevel methods. It includes efficient spin adapted implementations of ground and excited singlet states, as well as equation of motion oscillator strengths, for CCS, CC2, CCSD, and CC3. Furthermore, e^T provides unique capabilities such as multilevel Hartree-Fock and multilevel CC2, real-time propagation for CCS and CCSD, and efficient CC3 oscillator strengths. With a coupled cluster code based on an efficient Cholesky decomposition algorithm for the electronic repulsion integrals, e^T has similar advantages as codes using density fitting, but with strict error control. Here, we present the main features of the program and demonstrate its performance through example calculations. Because of its availability, performance, and unique capabilities, we expect e^T to become a valuable resource to the electronic structure community.

Kinetics of Hydrogen Abstraction Reactions of Methyl Palmitate and Octadecane by Hydrogen Atoms

Hydrogen abstractions play a crucial role in the consumption of fuel molecules during fuel pyrolysis and combustion processes. In this study, a generalized energy-based fragmentation approach was used to obtain CCSD(T)-F12a/cc-pVTZ energy barriers of hydrogen abstraction reactions by hydrogen atoms from methyl palmitate (C₁₅H₃₁COOCH₃), a key component of biodiesel. The accuracy of M06-2X/6-311++G(d,p) for obtaining the energy barriers was evaluated against the CCSD(T) results. Based on the quantum chemical results, the high-pressure-limit rate constants for C₁₅H₃₁COOCH₃ + H were calculated and compared with those of octadecane ( n-C₁₈H₃₈) reacting with H. The treatment of hindered internal rotations for such long-chain molecules was discussed and the rate rules for different abstraction sites were summarized. The results show that in the C₁₅H₃₁COOCH₃ + H system, the α hydrogen abstraction no longer plays a dominant role as in small methyl esters, and the hydrogen atoms of CH₂ groups far away from the ester group are more easily abstracted than those near the ester group.

THz vibrational spectroscopy for RNA basepair cocrystals and oligonucleotide sequences

Fourier infrared spectrometer and X-ray diffractometer were used to detect the spectra of lab-made U:A (uracil and adenine hydrate dried at room temperature), lab-made C:G (cytosine and guanine hydrate dried at room temperature), U + A (mixture of uracil and adenine), and C + G (mixture of cytosine and guanine). The results of our testing showed that U:A did not form a eutectic, but C:G did. In order to further characterize the vibrational modes of RNA base pair crystals, the absorption spectra of 1‑methyl‑5‑bromouracil‑9‑methyladenine (abbreviated as 1M5BU:9MA) and 1‑methylcytosine‑9‑ethylguanine (abbreviated as 1MC:9EG) were calculated based on the PBC-GEBF (generalized energy-based fragmentation approach under periodic boundary conditions) method. To further study the effect of substituents on the vibrational mode of the crystal structure, the substituents of 1M5BU:9MA and 1MC:9EG were artificially removed. The results of calculation brought out that methyl and ethyl as substituents have little effect on the vibrational spectrum, but halogen atoms such as Br atom in 1M5BU:9MA has a certain influence on the spectrum. Furthermore, THz (terahertz) spectra of the RNA nucleotide sequence 5'-AUCG-3' was analyzed. In the perspective of the spectra with DNA signatures, their spectra show a great deal of similarity regardless of RNA or DNA, or the base sequence difference. This study will provide a very important information for revealing the role of RNA in the formation of biological macromolecules and its mechanism.

A Fragment Quantum Mechanical Method for Metalloproteins

An accurate energy calculation of metalloprotein is of crucial importance and also a theoretical challenge. In this work, a metal molecular fractionation with conjugate caps (metal-MFCC) approach is developed for efficient linear-scaling quantum calculation of potential energy and atomic forces of metalloprotein. In this approach, the potential energy of a given protein is calculated by a linear combination of potential energies of the neighboring residues, two-body interaction energy between non-neighboring residues that are spatially in close contact and the potential energy of the metal binding group. The calculation of each fragment is embedded in a field of point charges representing the remaining protein environment. Numerical studies were carried out to check the performance of this method, and the calculated potential energies and atomic forces all show excellent agreement with the full system calculations at the M06-2X/6-31G(d) level. By combining the energy calculation with molecular dynamic simulation, we performed an ab initio structural optimization for a zinc finger protein with high efficiency. The present metal-MFCC approach is linear-scaling with a low prefactor and trivially parallelizable. The individual fragment typically contains about 50 atoms, and it is thus possible to be calculated at higher levels of the quantum chemistry method. This fragment method can be routinely applied to perform structural optimization and ab initio molecular dynamic simulation for metalloproteins of any size.

Libraries of Extremely Localized Molecular Orbitals. 3. Construction and Preliminary Assessment of the New Databanks

The fast and reliable determination of wave functions and electron densities of macromolecules has been one of the goals of theoretical chemistry for a long time, and in this context, several linear scaling techniques have been successfully devised over the years. Different approaches have been adopted to tackle this problem, and one of them exploits the fact that, according to the traditional chemical perception, molecules can be seen as constituted of recurring units (e.g., functional groups) with well-defined chemical features. This has led to the development of methods in which the global wave functions or electron densities of macromolecules are obtained by simply transferring density matrices or fuzzy electron densities associated with molecular fragments. In this context, we propose an alternative strategy that aims at quickly reconstructing wave functions and electron densities of proteins through the transfer of extremely localized molecular orbitals (ELMOs), which are orbitals strictly localized on small molecular units and, for this reason, easily transferable from molecule to molecule. To accomplish this task we have constructed original libraries of ELMOs that cover all the possible elementary fragments of the 20 natural amino acids in all their possible protonation states and forms. Our preliminary test calculations have shown that, compared to more traditional methods of quantum chemistry, the transfers from the novel ELMO databanks allow to obtain wave function and electron densities of large polypeptides and proteins at a significantly reduced computational cost. Furthermore, notwithstanding expected discrepancies, the obtained electron distributions and electrostatic potentials are in very good agreement with those obtained at Hartree-Fock and density functional theory (DFT) levels. Therefore, the results encourage to use the new libraries as alternatives to the popular pseudoatom-databases of crystallography in the refinement of crystallographic structures of macromolecules. In particular, in this context, we have already envisaged the coupling of the ELMO databanks with the promising Hirshfeld atom refinement technique to extend the applicability of the latter to very large systems.

Harmonizing accuracy and efficiency: A pragmatic approach to fragmentation of large molecules

Fragmentation methods offer an attractive alternative for ab initio treatment of large molecules and molecular clusters. However, balancing the accuracy and efficiency of these methods is a tight-rope-act. With this in view, we present an algorithm for automatic molecular fragmentation within Molecular Tailoring Approach (MTA) achieving this delicate balance. The automated code is tested out on a variety of molecules and clusters at the Hartree-Fock (HF)- and Møller-Plesset second order perturbation theory as well as density functional theory employing augmented Dunning basis sets. The results show remarkable accuracy and efficiency vis-à-vis the respective full calculations. Thus the present work forms an important step toward the development of an MTA-based black box code for implementation of HF as well as correlated quantum chemical calculations on large molecular systems.

THz spectra and corresponding vibrational modes of DNA base pair cocrystals and polynucleotides

The generalized energy-based fragmentation (GEBF) approach has been applied to study the THz spectra and vibrational modes of base pair cocrystals under periodic boundary conditions (denoted as PBC-GEBF). Results of vibrational mode reveal that hydrogen bonds play a pivotal role in the pairing process of base crystals, where most NH and CH bonds stretch to some extent. We also found that hydrogen bonds of a self-made A:T cocrystal completely break in a transition from liquid to the solid state, while self-made C:G cocrystal is different and easier to form a cocrystal, as confirmed by X-ray diffraction (XRD) and terahertz (THz) spectra. Furthermore, we have studied DNA polynucleotides (in both A and B forms) found that the vibrational modes changed a lot during the process of their forming double strand. Despite the key role played by hydrogen bonds, the key contribution originates from collective motions of the main skeleton. A comparative study of the spectra of some stranded fragments suggests that different sequences or forms have similar spectra in THz band. They distinguish from each other mainly in the low-frequency regions, especially below 1 THz. This study would make great contributions to the molecular dynamics model based DNA long-chain structure simulation in the future study.

Aggregation-induced visible light absorption makes reactant 1,2-diisocyanoarenes act as photosensitizers in double radical isocyanide insertions

The joint computational and experimental efforts reveal that the organic molecule 1,2-diisocyano-4,5-dimethylbenzene (1) acts as both a reactant and a photosensitizer (PS) in a metal-free reaction with perfluoroalkylhalide (2) to produce 2-perfluoroalkyl quinoxalines (3) under visible light. Both the π-π stacking aggregation in crystals and the solvation in various solvents of PS 1 exhibited visible-light absorption at 466 nm in spite of its smaller coefficient than that of the ultraviolet-light absorption. Such an aggregation-assisted visible-light absorption phenomenon is rationalized by theoretical calculations of the condensed-phase properties with the consideration of the explicit polarization effect from the neighboring molecules. Upon irradiation with different wavelengths, the emission colors changed from navy to bright yellow. Fluorescence lifetime measurements show that the emission of 1 comes from its singlet excited state. The aggregation induced emission when excited at 420 nm has a shorter lifetime (0.45 ns) than that of the emission from isolated molecules (2.71 ns) when excited at 381 nm. It is conceived that the aggregation assisted visible light absorption properties may be general in other photo-reactive molecules, such as 1,4-diisocyano-2,5-dimethylbenzene (4), 1,4-dicyanobenzene (5), and 1,4-diisocyanobenzene (6), which are widely used in many photochemical reactions in the absence of any external photosensitizer.

Accurate prediction of the structure and vibrational spectra of ionic liquid clusters with the generalized energy-based fragmentation approach: critical role of ion-pair-based fragmentation

The GEBF method with the ion-pair-based fragmentation has been developed to facilitate ab initio calculations of general ionic liquid clusters.

Molecular Mechanism of Self-Assembly of Aromatic Oligoamides into Interlocked Double-Helix Foldamers

Foldamer, inspired by the structures and functions of biopolymers, is defined as an artificial molecular architecture that can fold into a three-dimensional structure in solution and has been a growing and active field in supramolecular chemistry. The central issue in foldamer science is to understand how the primary sequence of oligomer folds into conformationally ordered structures as well as how individual subunits self-associate into assembly. For duplex structures, these two issues are always interrelated and inseparable with each other. Although the emergence of new foldamer keeps growing, the detailed mechanism remains elusive. On the basis of an artificially synthesized arylamide oligoamide foldamer with its crystal structure available, we constructed a set of four foldamers with a similar backbone but different substituents and aimed at dissecting the folding and self-association mechanisms of a double-helical foldamer with computations. Using molecular simulations at a microsecond time scale, we observed very consistent processes of the spontaneous self-assembly of two single-helical motifs into an entwined complex. Our results reveal that aggregation of two single-helical motifs driven by extensive π-π interactions is energetically favorable and that this spontaneous self-assembly proceeds through an "unwinding-threading-rewinding" mechanism. The detailed mechanisms about the folding and self-assembly in an aromatic oligoamide foldamer we present here disclose how the sequence is associated with a well-ordered three-dimensional structure at atomic level and therefore may have implications for designing new foldamers with versatile functions.

Terahertz spectra of DNA nucleobase crystals: A joint experimental and computational study

Terahertz (THz) spectra of DNA nucleobase crystals were experimentally studied by terahertz time domain spectroscopy (THz-TDS), Fourier transform infrared spectroscopy (FTIR), and computationally studied by the generalized energy-based fragmentation approach under periodic boundary conditions (denoted as PBC-GEBF). We analyzed the vibrational spectra of solid-state DNA nucleobases and assigned the corresponding vibrational modes to the main peaks in the experimental spectra with the PBC-GEBF results. The computational results were verified to be in good accordance with the experimental data. Harmonic vibrational frequency results revealed that all the vibrational modes belong to collective vibrational modes, which involve complicated mixtures of inter- and intramolecular displacements, somewhere in the vicinity of 0.5-9THz.

Quantum mechanical calculation of nanomaterial-ligand interaction energies by molecular fractionation with conjugated caps method

Molecular fractionation with conjugate caps (MFCC) method is introduced for the efficient estimation of quantum mechanical (QM) interaction energies between nanomaterial (carbon nanotube, fullerene, and graphene surface) and ligand (charged and neutral). In the calculations, nanomaterials are partitioned into small fragments and conjugated caps that are properly capped, and the interaction energies can be obtained through the summation of QM calculations of the fragments from which the contribution of the conjugated caps is removed. All the calculations were performed by density functional theory (DFT) and dispersion contributions for the attractive interactions were investigated by dispersion corrected DFT method. The predicted interaction energies by MFCC at each computational level are found to give excellent agreement with full system (FS) calculations with the mean energy deviation just a fractional kcal/mol. The accurate determination of nanomaterial-ligand interaction energies by MFCC suggests that it is an effective method for performing QM calculations on nanomaterial-ligand systems.

The MOD-QM/MM Method: Applications to Studies of Photosystem II and DNA G-Quadruplexes

Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and catalytic sites embedded in macrobiomolecules (eg, proteins and nucleic acids). QM/MM methodologies are highly efficient since they implement quantum chemistry methods for modeling only the portion of the system involving bond-breaking/forming processes (QM layer), as influenced by the surrounding molecular environment described in terms of molecular mechanics force fields (MM layer). Some of the limitations of QM/MM methods when polarization effects are not explicitly considered include the approximate treatment of electrostatic interactions between QM and MM layers. Here, we review recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in biomacromolecules and structural refinement, beyond the common limitations of QM/MM hybrid methods. We focus on photosystem II (PSII) with emphasis on the description of the oxygen-evolving complex (OEC) and its high-resolution extended X-ray absorption fine structure spectra (EXAFS) in conjunction with Monte Carlo structural refinement. Furthermore, we review QM/MM structural refinement studies of DNA G4 quadruplexes with embedded monovalent cations and direct comparisons to NMR data.