1
|
Wang W, Yan D, Cai Y, Xu D, Ma J, Wang Q. General Charge Transfer Dipole Model for AMOEBA-Like Force Fields. J Chem Theory Comput 2023; 19:2518-2534. [PMID: 37125725 DOI: 10.1021/acs.jctc.2c01084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
The development of highly accurate force fields is always an importance aspect in molecular modeling. In this work, we introduce a general damping-based charge transfer dipole (D-CTD) model to describe the charge transfer energy and the corresponding charge flow for H, C, N, O, P, S, F, Cl, and Br elements in common bio-organic systems. Then, two effective schemes to evaluate the charge flow from the corresponding induced dipole moment between the interacting molecules were also proposed and discussed. The potential applicability of the D-CTD model in ion-containing systems was also demonstrated in a series of ion-water complexes including Li+, Na+, K+, Mg2+, Ca2+, Fe2+, Zn2+, Pt2+, F-, Cl-, Br-, and I- ions. In general, the D-CTD model demonstrated good accuracy and good transferability in both charge transfer energy and the corresponding charge flow for a wide range of model systems. By distinguishing the intermolecular charge redistribution (charge transfer) under the influence of an external electric field from the accompanying intramolecular charge redistribution (polarization), the D-CTD model is theoretically consistent with current induced dipole-based polarizable dipole models and hence can be easily implemented and parameterized. Along with our previous work in charge penetration-corrected electrostatics, a bottom-up approach constructed water model was also proposed and demonstrated. The structure-maker and structure-breaker roles of cations and anions were also correctly reproduced using Na+, K+, Cl-, and I- ions in the new water model, respectively. This work demonstrates a cost-effective approach to describe the charge transfer phenomena. The water and ion models also show the feasibility of a modulated development approach for future force fields.
Collapse
Affiliation(s)
- Wei Wang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Dengjie Yan
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu 610041, China
| | - Yao Cai
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Dingguo Xu
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Jianyi Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Qiantao Wang
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu 610041, China
| |
Collapse
|
2
|
Aldossary A, Gimferrer M, Mao Y, Hao H, Das AK, Salvador P, Head-Gordon T, Head-Gordon M. Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions. J Phys Chem A 2023; 127:1760-1774. [PMID: 36753558 DOI: 10.1021/acs.jpca.2c08061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.
Collapse
Affiliation(s)
- Abdulrahman Aldossary
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Martí Gimferrer
- Institut de Química Computacional i Catàlsi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia Spain
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Hongxia Hao
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Akshaya K Das
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Pedro Salvador
- Institut de Química Computacional i Catàlsi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia Spain
| | - Teresa Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Martin Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| |
Collapse
|
3
|
Unveiling the reactivity of truxillic and truxinic acids (TXAs): deprotonation, anion…H–O, cation…O and cation…$${\varvec{\pi}}$$ interactions in TXA0…Y+ and TXA0…Z− complexes (Y = Li, Na, K; Z = F, Cl, Br). Struct Chem 2022. [DOI: 10.1007/s11224-022-01965-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
4
|
Intramolecular resonance-assisted hydrogen bonds: Insights from symmetry adapted perturbation theory. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111474] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
5
|
Das AK, Liu M, Head-Gordon T. Development of a Many-Body Force Field for Aqueous Alkali Metal and Halogen Ions: An Energy Decomposition Analysis Guided Approach. J Chem Theory Comput 2022; 18:953-967. [DOI: 10.1021/acs.jctc.1c00537] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Akshaya Kumar Das
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Meili Liu
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
- Department of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Teresa Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
- Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
6
|
Herbert JM, Paul SK. Interaction Energy Analysis of Monovalent Inorganic Anions in Bulk Water Versus Air/Water Interface. Molecules 2021; 26:6719. [PMID: 34771133 PMCID: PMC8588468 DOI: 10.3390/molecules26216719] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2021] [Revised: 10/29/2021] [Accepted: 11/03/2021] [Indexed: 11/16/2022] Open
Abstract
Soft anions exhibit surface activity at the air/water interface that can be probed using surface-sensitive vibrational spectroscopy, but the structural implications of this surface activity remain a matter of debate. Here, we examine the nature of anion-water interactions at the air/water interface using a combination of molecular dynamics simulations and quantum-mechanical energy decomposition analysis based on symmetry-adapted perturbation theory. Results are presented for a set of monovalent anions, including Cl-, Br-, I-, CN-, OCN-, SCN-, NO2-, NO3-, and ClOn- (n=1,2,3,4), several of which are archetypal examples of surface-active species. In all cases, we find that average anion-water interaction energies are systematically larger in bulk water although the difference (with respect to the same quantity computed in the interfacial environment) is well within the magnitude of the instantaneous fluctuations. Specifically for the surface-active species Br-(aq), I-(aq), ClO4-(aq), and SCN-(aq), and also for ClO-(aq), the charge-transfer (CT) energy is found to be larger at the interface than it is in bulk water, by an amount that is greater than the standard deviation of the fluctuations. The Cl-(aq) ion has a slightly larger CT energy at the interface, but NO3-(aq) does not; these two species are borderline cases where consensus is lacking regarding their surface activity. However, CT stabilization amounts to <20% of the total induction energy for each of the ions considered here, and CT-free polarization energies are systematically larger in bulk water in all cases. As such, the role of these effects in the surface activity of soft anions remains unclear. This analysis complements our recent work suggesting that the short-range solvation structure around these ions is scarcely different at the air/water interface from what it is in bulk water. Together, these observations suggest that changes in first-shell hydration structure around soft anions cannot explain observed surface activities.
Collapse
Affiliation(s)
- John M. Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, OH 43210, USA;
| | | |
Collapse
|
7
|
Herbert JM. Neat, Simple, and Wrong: Debunking Electrostatic Fallacies Regarding Noncovalent Interactions. J Phys Chem A 2021; 125:7125-7137. [PMID: 34388340 DOI: 10.1021/acs.jpca.1c05962] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Multipole moments such as charge, dipole, and quadrupole are often invoked to rationalize intermolecular phenomena, but a low-order multipole expansion is rarely a valid description of electrostatics at the length scales that characterize nonbonded interactions. This is illustrated by examining several common misunderstandings rooted in erroneous electrostatic arguments. First, the notion that steric repulsion originates in Coulomb interactions is easily disproved by dissecting the interaction potential for Ar2. Second, the Hunter-Sanders model of π-π interactions, which is based on quadrupolar electrostatics, is shown to have no basis in accurate calculations. Third, curved "buckybowls" exhibit unusually large dipole moments, but these are ancillary to the forces that control their intermolecular interactions, as illustrated by two examples involving corannulene. Finally, the assumption that interactions between water and small anions are dictated by the dipole moment of H2O is shown to be false in the case of binary halide-water complexes. These examples present a compelling case that electrostatic explanations based on low-order multipole moments are very often counterfactual for nonbonded interactions at close range and should not be taken seriously in the absence of additional justification.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
8
|
Dong X, Yuan X, Song Z, Wang Q. The development of an Amber-compatible organosilane force field for drug-like small molecules. Phys Chem Chem Phys 2021; 23:12582-12591. [PMID: 34037028 DOI: 10.1039/d1cp01169c] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
As members of the group IVA elements, silicon and carbon have long been thought of as isosteres of each other in drug design. However, the lack of silicon parameters in current main stream force fields hinders the computational study of this important element in drug discovery. Thus, in this study, we attempted to supplement the parameters of organosilanes in the General Amber Force Field (GAFF2). The parameters have been designed following the principles of GAFF2 to make it compatible with the Amber force field family. The accuracy of the parameters was discussed by comparing the pair interaction energy, the liquid properties, and the structures and alchemical binding free energy differences for a set of protein-ligand complexes.
Collapse
Affiliation(s)
- Xue Dong
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu, 610041, P. R. China.
| | - Xinghang Yuan
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu, 610041, P. R. China.
| | - Zhenlei Song
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu, 610041, P. R. China.
| | - Qiantao Wang
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu, 610041, P. R. China.
| |
Collapse
|
9
|
Elgengehi SM, El-Taher S, Ibrahim MA, El-Kelany KE. Unexpected favourable noncovalent interaction between chlorine oxyanions (ClO−; x = 1–4) and benzene: Benchmarking DFT and SAPT methods with respect to CCSD(T). COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113214] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
10
|
Laplaza R, Peccati F, A. Boto R, Quan C, Carbone A, Piquemal J, Maday Y, Contreras‐García J. NCIPLOT
and the analysis of noncovalent interactions using the reduced density gradient. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1497] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Rubén Laplaza
- CNRS, Laboratoire de Chimie Théorique, LCT Sorbonne Université Paris France
- Departamento de Química Física Universidad de Zaragoza Zaragoza Spain
| | - Francesca Peccati
- CNRS, Laboratoire de Chimie Théorique, LCT Sorbonne Université Paris France
- Institut des Sciences du Calcul et des Données, ISCD, Sorbonne Université Paris France
| | - Roberto A. Boto
- CNRS, Laboratoire de Chimie Théorique, LCT Sorbonne Université Paris France
- Centro de Física de Materiales CFM‐MPC (CSIC‐UPV/EHU) Donostia Spain
| | - Chaoyu Quan
- Institut des Sciences du Calcul et des Données, ISCD, Sorbonne Université Paris France
- SUSTech International Center for Mathematics, and Department of Mathematics Southern University of Science and Technology Shenzhen China
| | - Alessandra Carbone
- CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative (LCQB) Sorbonne Université Paris France
- Institut Universitaire de France Paris France
| | - Jean‐Philip Piquemal
- CNRS, Laboratoire de Chimie Théorique, LCT Sorbonne Université Paris France
- Institut Universitaire de France Paris France
| | - Yvon Maday
- SUSTech International Center for Mathematics, and Department of Mathematics Southern University of Science and Technology Shenzhen China
- Institut Universitaire de France Paris France
| | | |
Collapse
|
11
|
de Oliveira BG, Zabardasti A, do Rego DG, Pour MM. The formation of H···X hydrogen bond, C···X carbon-halide or Si···X tetrel bonds on the silylene-halogen dimers (X = F or Cl): intermolecular strength, molecular orbital interactions and prediction of covalency. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-02644-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
12
|
Garcia J, Podeszwa R, Szalewicz K. SAPT codes for calculations of intermolecular interaction energies. J Chem Phys 2020; 152:184109. [PMID: 32414261 DOI: 10.1063/5.0005093] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.
Collapse
Affiliation(s)
- Javier Garcia
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia at Katowice, Szkolna 9, Katowice, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
13
|
Inakollu VS, Geerke DP, Rowley CN, Yu H. Polarisable force fields: what do they add in biomolecular simulations? Curr Opin Struct Biol 2020; 61:182-190. [PMID: 32044671 DOI: 10.1016/j.sbi.2019.12.012] [Citation(s) in RCA: 56] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 12/18/2019] [Accepted: 12/19/2019] [Indexed: 12/11/2022]
Abstract
The quality of biomolecular simulations critically depends on the accuracy of the force field used to calculate the potential energy of the molecular configurations. Currently, most simulations employ non-polarisable force fields, which describe electrostatic interactions as the sum of Coulombic interactions between fixed atomic charges. Polarisation of these charge distributions is incorporated only in a mean-field manner. In the past decade, extensive efforts have been devoted to developing simple, efficient, and yet generally applicable polarisable force fields for biomolecular simulations. In this review, we summarise the latest developments in accounting for key biomolecular interactions with polarisable force fields and applications to address challenging biological questions. In the end, we provide an outlook for future development in polarisable force fields.
Collapse
Affiliation(s)
- Vs Sandeep Inakollu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong NSW 2522, Australia; Molecular Horizons, University of Wollongong, Wollongong NSW 2522 Australia; Illawarra Health and Medical Research Institute, Wollongong NSW 2522, Australia
| | - Daan P Geerke
- AIMMS Division of Molecular and Computational Toxicology, Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, the Netherlands.
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador, Canada.
| | - Haibo Yu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong NSW 2522, Australia; Molecular Horizons, University of Wollongong, Wollongong NSW 2522 Australia; Illawarra Health and Medical Research Institute, Wollongong NSW 2522, Australia.
| |
Collapse
|
14
|
Kruse H, Šponer J. Revisiting the Potential Energy Surface of the Stacked Cytosine Dimer: FNO-CCSD(T) Interaction Energies, SAPT Decompositions, and Benchmarking. J Phys Chem A 2019; 123:9209-9222. [PMID: 31560201 DOI: 10.1021/acs.jpca.9b05940] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Nucleobase stacking interactions are crucial for the stability of nucleic acids. This study investigates base stacking energies of the cytosine homodimer in different configurations, including intermolecular separation plots, detailed twist dependence, and displaced structures. Highly accurate ab initio quantum chemical single point energies using an energy function based on MP2 complete basis set extrapolation ([6 → 7]ZaPa-NR) and a CCSD(T)/cc-pVTZ-F12 high-level correction are presented as new reference data, providing the most accurate stacking energies of nucleobase dimers currently available. Accurate SAPT2+(3)δMP2 energy decomposition is used to obtain detailed insights into the nature of base stacking interactions at varying vertical distances and twist values. The ab initio symmetry adapted perturbation theory (SAPT) energy decomposition suggests that the base stacking originates from an intricate interplay between dispersion attraction, short-range exchange-repulsion, and Coulomb interaction. The interpretation of the SAPT data is a complex issue as key energy terms vary substantially in the region of optimal (low energy) base stacking geometries. Thus, attempts to highlight one leading stabilizing SAPT base stacking term may be misleading and the outcome strongly depends on the used geometries within the range of geometries sampled in nucleic acids upon thermal fluctuations. Modern dispersion-corrected density functional theory (among them DSD-BLYP-D3, ωB97M-V, and ωB97M-D3BJ) is benchmarked and often reaches up to spectroscopic accuracy (below 1 kJ/mol). The classical AMBER force field is benchmarked with multiple different sets of point-charges (e.g. HF, DFT, and MP2-based) and is found to produce reasonable agreement with the benchmark data.
Collapse
Affiliation(s)
- Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , CZ-61265 Brno , Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , CZ-61265 Brno , Czech Republic.,Central European Institute of Technology , Masaryk University , Kamenice 753/5 , 62500 Brno , Czech Republic
| |
Collapse
|
15
|
Das AK, Urban L, Leven I, Loipersberger M, Aldossary A, Head-Gordon M, Head-Gordon T. Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis. J Chem Theory Comput 2019; 15:5001-5013. [DOI: 10.1021/acs.jctc.9b00478] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
16
|
Leduc T, Aubert E, Espinosa E, Jelsch C, Iordache C, Guillot B. Polarization of Electron Density Databases of Transferable Multipolar Atoms. J Phys Chem A 2019; 123:7156-7170. [PMID: 31294565 DOI: 10.1021/acs.jpca.9b05051] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Polarizability is a key molecular property involved in either macroscopic (i.e., dielectric constant) and microscopic properties (i.e., interaction energies). In rigid molecules, this property only depends on the ability of the electron density (ED) to acquire electrostatic moments in response to applied electric fields. Databases of transferable electron density fragments are a cheap and efficient way to access molecular EDs. This approach is rooted in the relative conservation of the atomic ED between different molecules, termed transferability principle. The present work discusses the application of this transferability principle to the polarizability, an electron density-derived property, partitioned in atomic contributions using the Quantum Theory of Atoms In Molecules topology. The energetic consequences of accounting for in situ deformation (polarization) of database multipolar atoms are investigated in detail by using a high-quality quantum chemical benchmark.
Collapse
Affiliation(s)
- Théo Leduc
- Université de Lorraine, CNRS, CRM2 , F-54000 Nancy , France
| | | | | | | | | | - Benoît Guillot
- Université de Lorraine, CNRS, CRM2 , F-54000 Nancy , France
| |
Collapse
|
17
|
Liu C, Piquemal JP, Ren P. AMOEBA+ Classical Potential for Modeling Molecular Interactions. J Chem Theory Comput 2019; 15:4122-4139. [PMID: 31136175 DOI: 10.1021/acs.jctc.9b00261] [Citation(s) in RCA: 76] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Classical potentials based on isotropic and additive atomic charges have been widely used to model molecules in computers for the past few decades. The crude approximations in the underlying physics are hindering both their accuracy and transferability across chemical and physical environments. Here we present a new classical potential, AMOEBA+, to capture essential intermolecular forces, including permanent electrostatics, repulsion, dispersion, many-body polarization, short-range charge penetration, and charge transfer, by extending the polarizable multipole-based AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) model. For a set of common organic molecules, we show that AMOEBA+ with general parameters can reproduce both quantum mechanical interactions and energy decompositions according to Symmetry-Adapted Perturbation Theory (SAPT). Additionally, a new water model based on the AMOEBA+ framework captures various liquid-phase properties in molecular dynamics simulations while remaining consistent with SAPT energy decompositions, utilizing both ab initio data and experimental liquid properties. Our results demonstrate that it is possible to improve the physical basis of classical force fields to advance their accuracy and general applicability.
Collapse
Affiliation(s)
- Chengwen Liu
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States.,Laboratoire de Chimie Théorique , Sorbonne Université, UMR7616 CNRS , Paris 75252 , France.,Institut Universitaire de France , Paris Cedex 05, 75005 , France
| | - Pengyu Ren
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
18
|
Mao Y, Ge Q, Horn PR, Head-Gordon M. On the Computational Characterization of Charge-Transfer Effects in Noncovalently Bound Molecular Complexes. J Chem Theory Comput 2018; 14:2401-2417. [PMID: 29614855 DOI: 10.1021/acs.jctc.7b01256] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a "CT-free" state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the "CT-free" state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal-carbonyl complexes, and complexes formed by water and metal cations. For most of these systems, CDFT yields a much smaller equilibrium CT energy compared to that given by the ALMO-based definition. This is mainly because the CDFT population constraint does not fully inhibit CT, which means that the CDFT "CT-free" state is in fact CT-contaminated. Examples of this contamination include (i) matching forward and backward donation (e.g., formic acid dimer) and (ii) unidirectional CT without changing fragment populations. The magnitude of the latter effect is quantified in systems such as the water dimer by employing a 3-space density constraint in addition to the orbital constraint. Furthermore, by means of the adiabatic EDA, it is shown that several observable effects of CT, such as the "pyramidalization" of the planar BH3 molecule upon the complexation with Lewis bases, already appear on the "CT-free" CDFT surface. These results reveal the essential distinctions between the ALMO and CDFT definitions of CT and suggest that the former is more consistent with accepted understanding of the role of CT in intermolecular binding.
Collapse
Affiliation(s)
- Yuezhi Mao
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California at Berkeley , Berkeley , California 94720 , United States
| | - Qinghui Ge
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California at Berkeley , Berkeley , California 94720 , United States.,Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Paul R Horn
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California at Berkeley , Berkeley , California 94720 , United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California at Berkeley , Berkeley , California 94720 , United States.,Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
19
|
Thaunay F, Jana C, Clavaguéra C, Ohanessian G. Strategy for Modeling the Infrared Spectra of Ion-Containing Water Drops. J Phys Chem A 2018; 122:832-842. [PMID: 29266957 DOI: 10.1021/acs.jpca.7b10554] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Hydrated ions are ubiquitous in environmental and biological media. Understanding the perturbation exerted by an ion on the water hydrogen bond network is possible in the nanodrop regime by recording vibrational spectra in the O-H bond stretching region. This has been achieved experimentally in recent years by forming gaseous ions containing tens to hundreds of water molecules and recording their infrared photodissociation spectra. In this paper, we demonstrate the capabilities of a modeling strategy based on an extension of the AMOEBA polarizable force field to implement water atomic charge fluctuations along with those of intramolecular structure along the dynamics. This supplementary flexibility of nonbonded interactions improves the description of the hydrogen bond network and, therefore, the spectroscopic response. Finite temperature IR spectra are obtained from molecular dynamics simulations by computing the Fourier transform of the dipole moment autocorrelation function. Simulations of 1-2 ns are required for extensive sampling in order to reproduce the experimental spectra. Furthermore, bands are assigned with the driven molecular dynamics approach. This method package is shown to compare successfully with experimental spectra for 11 ions in water drops containing 36-100 water molecules. In particular, band frequency shifts of the free O-H stretching modes at the cluster surface are well reproduced as a function of both ion charge and drop size.
Collapse
Affiliation(s)
- Florian Thaunay
- LCM, CNRS, Ecole Polytechnique, Université Paris Saclay , 91128 Palaiseau, France
| | - Chandramohan Jana
- LCM, CNRS, Ecole Polytechnique, Université Paris Saclay , 91128 Palaiseau, France
| | - Carine Clavaguéra
- Laboratoire de Chimie Physique, Université Paris Sud - CNRS, Université Paris Saclay , 15, avenue Jean Perrin, 91405 Orsay Cedex, France
| | - Gilles Ohanessian
- LCM, CNRS, Ecole Polytechnique, Université Paris Saclay , 91128 Palaiseau, France
| |
Collapse
|