1
|
Romero J, Limão-Vieira P, Maihom T, Hermansson K, Probst M. A polarizable valence electron density based force field for high-energy interactions between atoms and molecules. J Chem Phys 2024; 160:235101. [PMID: 38904408 DOI: 10.1063/5.0210949] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Accepted: 05/29/2024] [Indexed: 06/22/2024] Open
Abstract
High-accuracy molecular force field models suited for hot gases and plasmas are not as abundant as those geared toward ambient pressure and temperature conditions. Here, we present an improved version of our previous electron-density based force field model that can now account for polarization effects by adjusting the atomic valence electron contributions to match ab initio calculated Mulliken partial charges. Using a slightly modified version of the Hohenberg-Kohn theorem, we also include an improved theoretical formulation of our model when applied to systems with degenerate ground states. We present two variants of our polarizable model, fitted from ab initio reference data calculated at CCSD(T)/cc-pVTZ and CCSD(T)/CEP-31G levels of theory, that both accurately model water dimer interaction energies. Further improvements include the additional interaction components with fictitious non-spherically symmetric, yet atom-centered, electron densities and fitting the exchange and correlation coefficients against analytical expressions. The latter removes all unphysical oscillations that are observed in the previous non-polarizable variant of our force field.
Collapse
Affiliation(s)
- José Romero
- Institute of Ion Physics and Applied Physics, University of Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
- Atomic and Molecular Collisions Laboratory, CEFITEC, Department of Physics, Universidade NOVA de Lisboa, 2829-516 Caparica, Portugal
| | - Paulo Limão-Vieira
- Atomic and Molecular Collisions Laboratory, CEFITEC, Department of Physics, Universidade NOVA de Lisboa, 2829-516 Caparica, Portugal
| | - Thana Maihom
- School of Molecular Science and Engineering, Vidyasirimedhi Institute of Science and Technology, Rayong 21210, Thailand
- Department of Chemistry, Faculty of Liberal Arts and Science, Kasetsart University, Kamphaeng Saen Campus, Nakhon Pathom 73140, Thailand
| | - Kersti Hermansson
- Department of Chemistry-Ångström, Uppsala University, Box 538, SE-75121 Uppsala, Sweden
| | - Michael Probst
- Institute of Ion Physics and Applied Physics, University of Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
- School of Molecular Science and Engineering, Vidyasirimedhi Institute of Science and Technology, Rayong 21210, Thailand
| |
Collapse
|
2
|
Romero J, Limão-Vieira P, Hermansson K, Probst M. A Simple Electron-Density Based Force Field Model for High-Energy Interactions between Atoms and Molecules. J Phys Chem A 2024; 128:1163-1172. [PMID: 38320398 PMCID: PMC10875675 DOI: 10.1021/acs.jpca.3c06724] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 01/11/2024] [Accepted: 01/11/2024] [Indexed: 02/08/2024]
Abstract
In high-energy molecular dynamics or Monte Carlo simulations, standard force fields optimized for simulations at ambient temperatures are inadequate. This is largely because their repulsive parts have been regarded as not very significant, even well below zero interaction energies. It is, therefore, not obvious which force fields to resort to for simulating hot gases or plasmas. A force field model that uses the electronic densities of noninteracting atoms or molecules within the pair approximation is introduced. We start by deriving a naïve model that neglects any exchange and correlation effects between the electronic clouds and then correct this model by adding a term calibrated from ab initio calculations using the CCSD(T)/cc-pVTZ level of theory. The resulting expression for this term can be regarded as a simple exchange-correlation function. We compare the results for the repulsive part of the potential energy hypersurfaces with the force fields commonly used on some dimers of small molecules.
Collapse
Affiliation(s)
- José Romero
- Institute
of Ion Physics and Applied Physics, University
of Innsbruck, Technikerstraße 25, Innsbruck 6020, Austria
- Atomic
and Molecular Collisions Laboratory, CEFITEC, Department of Physics, Universidade NOVA de Lisboa, Caparica 2829-516, Portugal
| | - Paulo Limão-Vieira
- Atomic
and Molecular Collisions Laboratory, CEFITEC, Department of Physics, Universidade NOVA de Lisboa, Caparica 2829-516, Portugal
| | - Kersti Hermansson
- Department
of Chemistry-Ångström, Uppsala
University, Uppsala SE-75121, Sweden
| | - Michael Probst
- Institute
of Ion Physics and Applied Physics, University
of Innsbruck, Technikerstraße 25, Innsbruck 6020, Austria
- School
of Molecular Science and Engineering, Vidyasirimedhi
Institute of Science and Technology, Rayong 21210, Thailand
| |
Collapse
|
3
|
Naseem-Khan S, Lagardère L, Narth C, Cisneros GA, Ren P, Gresh N, Piquemal JP. Development of the Quantum-Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed-Phase Molecular Dynamics Simulations. J Chem Theory Comput 2022; 18:3607-3621. [PMID: 35575306 PMCID: PMC10851344 DOI: 10.1021/acs.jctc.2c00029] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
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.
Collapse
Affiliation(s)
- Sehr Naseem-Khan
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Louis Lagardère
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- IP2CT, FR 2622, CNRS, Sorbonne Université, 75005 Paris, France
| | | | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Nohad Gresh
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
| | - Jean-Philip Piquemal
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
4
|
Hayakawa D, Terauchi N, Iwasaki A, Watanabe Y, Gouda H. Systematic preparation method of a molecular model explicitly describing electron distributions for halogen bonds. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139754] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|