1
|
Hashem Y, Foust K, Kaledin M, Kaledin AL. Fitting Potential Energy Surfaces by Learning the Charge Density Matrix with Permutationally Invariant Polynomials. J Chem Theory Comput 2023; 19:5690-5700. [PMID: 37561135 PMCID: PMC10501011 DOI: 10.1021/acs.jctc.3c00586] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2023] [Indexed: 08/11/2023]
Abstract
The electronic energy in the Hartree-Fock (HF) theory is the trace of the product of the charge density matrix (CDM) with the one-electron and two-electron matrices represented in an atomic orbital basis, where the two-electron matrix is also a function of the same CDM. In this work, we examine a formalism of analytic representation of a generic molecular potential energy surface (PES) as a sum of a linearly parameterized HF and a correction term, the latter formally representing the electron correlation energy, also linearly parameterized, by expressing the elements of CDM using permutationally invariant polynomials (PIPs). We show on a variety of numerical examples, ranging from exemplary two-electron systems HeH+ and H3+ to the more challenging cases of methanium (CH5+) fragmentation and high-energy tautomerization of formamide to formimidic acid that such a formulation requires significantly fewer, 10-20% of PIPs, to accomplish the same accuracy of the fit as the conventional representation at practically the same computational cost.
Collapse
Affiliation(s)
- Younos Hashem
- Department
of Chemistry & Biochemistry, Kennesaw
State University, 370 Paulding Ave NW, Box # 1203, Kennesaw 30144, Georgia
| | - Katheryn Foust
- Department
of Chemistry & Biochemistry, Kennesaw
State University, 370 Paulding Ave NW, Box # 1203, Kennesaw 30144, Georgia
| | - Martina Kaledin
- Department
of Chemistry & Biochemistry, Kennesaw
State University, 370 Paulding Ave NW, Box # 1203, Kennesaw 30144, Georgia
| | - Alexey L. Kaledin
- Cherry
L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, 1515 Dickey Drive, Atlanta 30322, Georgia
| |
Collapse
|
2
|
Riera M, Knight C, Bull-Vulpe EF, Zhu X, Agnew H, Smith DGA, Simmonett AC, Paesani F. MBX: A many-body energy and force calculator for data-driven many-body simulations. J Chem Phys 2023; 159:054802. [PMID: 37526156 PMCID: PMC10550339 DOI: 10.1063/5.0156036] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 07/11/2023] [Indexed: 08/02/2023] Open
Abstract
Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.
Collapse
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Argonne National Laboratory, Computational Science Division, Lemont, Illinois 60439, USA
| | - Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Xuanyu Zhu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | | | - Andrew C. Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | | |
Collapse
|
3
|
Heindel JP, Herman KM, Xantheas SS. Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development. Annu Rev Phys Chem 2023; 74:337-360. [PMID: 37093659 DOI: 10.1146/annurev-physchem-062422-023532] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born-Oppenheimer potential energy surface. In particular, we focus on ab initio-based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through (a) parameterization of distributed multipoles, (b) explicit fitting of the MBE, (c) inclusion of many-atom features in a neural network, and (d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington, USA; ,
| |
Collapse
|
4
|
Bull-Vulpe EF, Riera M, Bore SL, Paesani F. Data-Driven Many-Body Potential Energy Functions for Generic Molecules: Linear Alkanes as a Proof-of-Concept Application. J Chem Theory Comput 2022. [PMID: 36113028 DOI: 10.1021/acs.jctc.2c00645] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a generalization of the many-body energy (MB-nrg) theoretical/computational framework that enables the development of data-driven potential energy functions (PEFs) for generic covalently bonded molecules, with arbitrary quantum mechanical accuracy. The "nearsightedness of electronic matter" is exploited to define monomers as "natural building blocks" on the basis of their distinct chemical identity. The energy of generic molecules is then expressed as a sum of individual many-body energies of incrementally larger subsystems. The MB-nrg PEFs represent the low-order n-body energies, with n = 1-4, using permutationally invariant polynomials derived from electronic structure data carried out at an arbitrary quantum mechanical level of theory, while all higher-order n-body terms (n > 4) are represented by a classical many-body polarization term. As a proof-of-concept application of the general MB-nrg framework, we present MB-nrg PEFs for linear alkanes. The MB-nrg PEFs are shown to accurately reproduce reference energies, harmonic frequencies, and potential energy scans of alkanes, independently of their length. Since, by construction, the MB-nrg framework introduced here can be applied to generic covalently bonded molecules, we envision future computer simulations of complex molecular systems using data-driven MB-nrg PEFs, with arbitrary quantum mechanical accuracy.
Collapse
Affiliation(s)
- Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Sigbjørn L. Bore
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
5
|
Schmidt M, Roy PN. On the accuracy and efficiency of different methods to calculate Raman vibrational shifts of parahydrogen clusters. J Chem Phys 2022; 156:084102. [DOI: 10.1063/5.0076403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Raman vibrational frequency shifts of pure parahydrogen and orthodeuterium clusters of sizes N = 4–9 are calculated using the Langevin equation path integral ground state method. The shifts are calculated using three different methods; the results obtained from each are compared to experiment and variance properties are assessed. The first method requires the direct calculation of energies from two simulations: one when the cluster is in the v = 0 vibrational state and one when the cluster has v = 1 total quantum of vibration. The shift is directly calculated from the difference in those two energies. The second method requires only a v = 0 simulation to be performed. The ground state energy is calculated as usual and the excited state energy is calculated by using the distribution of the v = 0 simulation and the ratio of the density matrices between the v = 1 state and the v = 0 state. The shift is calculated from the difference in those two energies. These first two are both exact methods. The final method is based on perturbation theory where the shift is calculated by averaging the pairwise difference potential over the pair distribution function. However, this is an approximate approach. It is found that for large enough system sizes, despite the approximations, the perturbation theory method has the strongest balance between accuracy and precision when weighing against computational cost.
Collapse
Affiliation(s)
- Matthew Schmidt
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
6
|
Houston PL, Qu C, Nandi A, Conte R, Yu Q, Bowman JM. Permutationally invariant polynomial regression for energies and gradients, using reverse differentiation, achieves orders of magnitude speed-up with high precision compared to other machine learning methods. J Chem Phys 2022; 156:044120. [DOI: 10.1063/5.0080506] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Paul L. Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA and Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Chen Qu
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, USA
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
7
|
Nandi A, Qu C, Houston PL, Conte R, Yu Q, Bowman JM. A CCSD(T)-Based 4-Body Potential for Water. J Phys Chem Lett 2021; 12:10318-10324. [PMID: 34662138 DOI: 10.1021/acs.jpclett.1c03152] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
High-level, ab initio calculations find that the 4-body (4-b) interaction is needed to account for near-100% of the total interaction energy for water clusters as large as the 21-mer. Motivated by this, we report a permutationally invariant polynomial potential energy surface (PES) for the 4-body interaction. This machine-learned PES is a fit to 2119 symmetry-unique, CCSD(T)-F12a/haTZ 4-b interaction energies. Configurations for these come from tetramer direct-dynamics calculations, fragments from an MD water simulation at 300 K, and tetramer fragments in a variety of water clusters. The PIP basis is purified to ensure that the PES goes rigorously to zero in monomer+trimer and dimer+dimer dissociations. The 4-b energies of isomers of the hexamer calculated with the new PES are shown to be in better agreement with benchmark CCSD(T) results than those from the MB-pol potential. Tests on larger clusters further validate the high-fidelity of the PES. The PES is shown to be fast to evaluate, taking 2.4 s for 105 evaluations on a single core of 2.4 GHz Intel Xeon processor, and significantly faster using a parallel version of the PES.
Collapse
Affiliation(s)
- Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chen Qu
- Department of Chemistry & Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
8
|
Felker PM, Lauvergnat D, Scribano Y, Benoit DM, Bačić Z. Intramolecular stretching vibrational states and frequency shifts of (H2)2 confined inside the large cage of clathrate hydrate from an eight-dimensional quantum treatment using small basis sets. J Chem Phys 2019; 151:124311. [DOI: 10.1063/1.5124051] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Peter M. Felker
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, USA
| | - David Lauvergnat
- Laboratoire de Chimie Physique, UMR-CNRS 8000, Université de Paris-Sud, Orsay F-91405, France
| | - Yohann Scribano
- Laboratoire Univers et Particule de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier Cedex, France
| | - David M. Benoit
- Department of Physics and Mathematics, E. A. Milne Centre for Astrophysics and G. W. Gray Centre for Advanced Materials, The University of Hull, Cottingham Road, Kingston upon Hull HU6 7RX, United Kingdom
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| |
Collapse
|
9
|
Lauvergnat D, Felker P, Scribano Y, Benoit DM, Bačić Z. H2, HD, and D2in the small cage of structure II clathrate hydrate: Vibrational frequency shifts from fully coupled quantum six-dimensional calculations of the vibration-translation-rotation eigenstates. J Chem Phys 2019; 150:154303. [DOI: 10.1063/1.5090573] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David Lauvergnat
- Laboratoire de Chimie Physique, UMR-CNRS 8000, Université de Paris-Sud, Orsay F-91405, France
| | - Peter Felker
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, USA
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier Cedex, France
| | - David M. Benoit
- Department of Physics and Mathematics, E.A. Milne Centre for Astrophysics and G. W. Gray Centre for Advanced Materials, The University of Hull, Cottingham Road, Kingston upon Hull HU6 7RX, United Kingdom
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
| |
Collapse
|