Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Zero-point energies from bond orders and populations relationships
We report two analytical quantum mechanics (QM) models for approximating appropriately scaled harmonic zero-point energies (ZPEs) without Hessian calculations. Following our earlier bond energies from bond orders and populations model that takes a similar form as an extended Hückel model but uses well-conditioned orbital populations, this work demonstrates a proof of concept for approximating ZPEs, an important component in thermochemistry calculations, while eschewing unfavorably scaling algorithms involving Hessian matrices. The ZPE-BOP1 model uses Mulliken orbital populations from hybrid Kohn–Sham density functional theory calculations within an extended Hückel-type model that defines vibrational bond energy terms using two atom-pairwise parameters that are fit to reproduce ZPEs from B3LYP calculations. The more accurate ZPE-BOP2 model uses Mulliken orbital populations from Hartree–Fock calculations within a different extended Hückel-type model that includes a short-range anharmonic energy term and a coupled three-body oscillator energy term with seven atom-pairwise parameters. Both models predict ZPEs in molecules involving first row elements, but ZPE-BOP2 outperforms ZPE-BOP1 in strained and long-chain molecules and provides ZPEs more competitive with those from semi-empirical QM methods (e.g., AM1, PM6, PM7, and XTB-2) that compute ZPEs with Hessian calculations. This work shows progress and an outlook toward computational models that use well-conditioned orbital populations to efficiently predict useful physicochemical properties. It also shows opportunities for approximate QM models that would shift traditional computational bottlenecks away from costly algorithms such as Hessian calculations to others that focus on reliable orbital populations.  more » « less
Award ID(s):
2050944
PAR ID:
10599695
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
AIP
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
162
Issue:
8
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, The Journal of Chemical Physics)
    Bonding energies play an essential role in describing the relative stability of molecules in chemical space. Therefore, methods employed to search chemical space need to capture the bonding behavior for a wide range of molecules, including radicals. In this work, we investigate the ability of quantum alchemy to capture the bonding behavior of hypothetical chemical compounds, specifically diatomic molecules involving hydrogen with various electronic structures. We evaluate equilibrium bond lengths, ionization energies, and electron affinities of these fundamental systems. We compare and contrast how well manual quantum alchemy calculations, i.e., quantum mechanics calculations in which the nuclear charge is altered, and quantum alchemy approximations using a Taylor series expansion can predict these molecular properties. Our results suggest that while manual quantum alchemy calculations outperform Taylor series approximations, truncations of Taylor series approximations after the second order provide the most accurate Taylor series predictions. Furthermore, these results suggest that trends in quantum alchemy predictions are generally dependent on the predicted property (i.e., equilibrium bond length, ionization energy, or electron affinity). Taken together, this work provides insight into how quantum alchemy predictions using a Taylor series expansion may be applied to future studies of non-singlet systems as well as the challenges that remain open for predicting the bonding behavior of such systems. 
    more » « less
  2. ; ; (, The Journal of Chemical Physics)
    QM/MM methods have been used to study electronic structure properties and chemical reactivity in complex molecular systems where direct electronic structure calculations are not feasible. In our previous work, we showed that non-polarizable force fields, by design, describe intermolecular interactions through pairwise interactions, overlooking many-body interactions involving three or more particles. In contrast, polarizable force fields account partially for many-body effects through polarization, but still handle van der Waals and permanent electrostatic interactions pairwise. We showed that despite those limitations, polarizable and non-polarizable force fields can reproduce relative cooperativity achieved using density functional theory due to error compensation mechanisms. In this contribution, we assess the performance of QM/MM methods in reproducing these phenomena. Our study highlights the significance of the QM region size and force field choice in QM/MM calculations, emphasizing the importance of parameter validation to obtain accurate interaction energy predictions. 
    more » « less
  3. ; ; ; ; ; ; (, CrystEngComm)
    null (Ed.)
    9-Anthracene carboxylic acid derivatives comprise a family of thermally reversible photomechanical molecular crystals. The photomechanical response relies on a [4 + 4] photodimerization followed by dissociation that occurs on timescales of seconds to minutes. A combined theoretical and experimental investigation is undertaken to better understand how chemical modification of the anthracene core influences energetics of both the isolated molecule and the crystal lattice. We use both density functional theory and dispersion-corrected Moller–Plesset perturbation theory computational methods to establish orbital energies, photodimerization reaction energies, and lattice energies for a set of substituted 9-anthracene carboxylic acid molecules. The calculations reveal that steric interactions play a dominant role in the ability to form photodimers and indicate an energetic threshold of 80–90 kJ per mole for the dimerization reaction. Examination of intermolecular bonding in a subset of fluorinated 9ACs revealed the absence of H⋯F intermolecular bond formation and energy differences that can explain observed trends in the dissociation kinetics and mechanical reset times. Fluorescence recovery after photobleaching experiments shows that the photodimer dissociation kinetics depend on the amount of initial photodimer, preventing a straightforward correlation between halogen atom substitution and dissociation rates using the Bell–Evans–Polanyi principle. The results clarify how molecular structure affects intermolecular interactions and photoreactivity in this family of molecular crystals, but the origin of the complex photodimer dissociation dynamics remains an open question. 
    more » « less
  4. ; ; ; ; (, ACS Physical Chemistry Au)
    With the increasing popularity of machine-learning (ML) applications, the demand for explainable artificial intelligence (XAI) techniques to explain ML models developed for computational chemistry has also emerged. In this study, we present the development of the Boltzmann-weighted Cumulative Integrated Gradients (BCIG) approach for effective explanation of mechanistic insights into ML models trained on high-level quantum mechanical and molecular mechanical (QM/MM) minimum energy pathways (MEPs). Using the acylation reactions of the Toho-1 β-lactamases and two antibiotic molecules (ampicillin and cefalexin) as the model systems, we show that the BCIG approach could quantitatively attribute the energetic contribution in one system, and the relative reactivity of individual steps across different systems to specific chemical processes such as the bond making/breaking and proton transfers. The proposed BCIG contribution attribution method quantifies chemist-interpretable insights in terms of contributions from each elementary chemical process, which are in agreement with the validating QM/MM calculations and our intuitive mechanistic understandings of the model reactions. 
    more » « less
  5. ; (, Zenodo)
    This dataset consists of 1,000 coordinate files (in the CHARMM psf/cor format) for the QM/MM minimum energy pathways of the deacylation reactions between a Class A beta-lactamases (GES-5) and the imipenem antibiotic molecules.</p> All pathway conformations were optimized at DFTB3/3OB-f/CHARMM36 level with 36 replicas.</p> All single point calculations and charge population analysis were done at B3LYP-D3/6-31+G(d,p)/CHARMM36 level.</p> 0.paths_ges_imi_d1.tar.gz: 500 pathway conformations for GES-5/IPM-Delta1 deacylation reactions.</li>0.paths_ges_imi_d2.tar.gz: 500 pathway conformations for GES-5/IPM-Delta1 deacylation reactions.</li>1.eners.zip: The single point replica energies along all GES-5/IPM pathways.</li>1.chrgs.zip: The NBO charges of the QM region of all replica conformations along all GES-5/IPM pathways.</li>2.datasets.zip: The Python codes to postprocess the molecular data and the featurized the NumPy arrays.</li>3.gnn.zip: The Python codes that implements the edge-conditioned graph convolutional NN to predict the deacylation barriers.</li>5.representative_conf.zip: The pathway conformations of all cluster centroids and an energetic representative (pathway id 22) pathway. Note: This file also serves as a peephole of how the pathway conformations from Reaction Path with Holonomic Constrains calculations looks like.</li>6.benchmark.zip: The benchmark calculations that validates the DFTB3/3OB-f/CHARMM36 against DFTB3/3OB/CHARMM36 and B3LYP/6-31G(d,p)/CHARMM36 level of theory on the energetic representative (pathway id 22) pathway conformations. </li></ul> 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email