Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Jun 14;18(6):3607-3621.
doi: 10.1021/acs.jctc.2c00029. Epub 2022 May 16.

Development of the Quantum-Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed-Phase Molecular Dynamics Simulations

Sehr Naseem-Khan  1   2 , Louis Lagardère  1   3 , Christophe Narth  1 , G Andrés Cisneros  2 , Pengyu Ren  4 , Nohad Gresh  1 , Jean-Philip Piquemal  1   4   5
Affiliations

Development of the Quantum-Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed-Phase Molecular Dynamics Simulations

Sehr Naseem-Khan et al. J Chem Theory Comput. .

Abstract

We present the extension of the Sum of Interactions Between Fragments Ab initio Computed (SIBFA) many-body polarizable force field to condensed-phase molecular dynamics (MD) simulations. The quantum-inspired SIBFA procedure is grounded on simplified integrals obtained from localized molecular orbital theory and achieves full separability of its intermolecular potential. It embodies long-range multipolar electrostatics (up to quadrupole) coupled to a short-range penetration correction (up to charge-quadrupole), exchange repulsion, many-body polarization, many-body charge transfer/delocalization, exchange dispersion, and dispersion (up to C10). This enables the reproduction of all energy contributions of ab initio symmetry-adapted perturbation theory (SAPT(DFT)) gas-phase reference computations. The SIBFA approach has been integrated within the Tinker-HP massively parallel MD package. To do so, all SIBFA energy gradients have been derived and the approach has been extended to enable periodic boundary conditions simulations using smooth particle mesh Ewald. This novel implementation also notably includes a computationally tractable simplification of the many-body charge transfer/delocalization contribution. As a proof of concept, we perform a first computational experiment defining a water model fitted on a limited set of SAPT(DFT) data. SIBFA is shown to enable a satisfactory reproduction of both gas-phase energetic contributions and condensed-phase properties highlighting the importance of its physically motivated functional form.

PubMed Disclaimer

Figures

Figure 1:
Figure 1:
Comparison of SIBFA21 and SAPT(DFT) energy components for the linear water dimer. The PBE0 functional and the aug-cc-pVTZ basis were used in SAPT(DFT) calculations.
Figure 2:
Figure 2:
Comparison of the radial distribution functions (O-O), (O-H) and (H-H) computed with SIBFA21 to the experiment at 298K.
Figure 3:
Figure 3:
Comparison of condensed phase properties (from 261K to 369K) computed with SIBFA21, AMOEBA+, AMOEBA14, MB-POL and MB-UCB-MDQ to the experiment.
See this image and copyright information in PMC

References

    1. Stone A The theory of intermolecular forces; Oxford University Press, Oxford, UK., 2013.
    1. Gresh N; Cisneros GA; Darden TA; Piquemal J-P Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy. J. Chem. Theory Comput 2007, 3, 1960–1986. - PMC - PubMed
    1. Warshel A; Kato M; Pisliakov AV Polarizable Force Fields: History, Test Cases, and Prospects. Journal of Chemical Theory and Computation 2007, 3, 2034–2045, - PubMed
    1. Shi Y; Ren P; Schnieders M; Piquemal J-P Reviews in Computational Chemistry Volume 28; John Wiley & Sons, Ltd, 2015; Chapter 2, pp 51–86.
    1. Melcr J; Piquemal J-P Accurate Biomolecular Simulations Account for Electronic Polarization. Frontiers in Molecular Biosciences 2019, 6, 143. - PMC - PubMed

LinkOut - more resources