1
|
Canales M, Guàrdia E. Computer simulation study of ion-water and water-water hydrogen bonds in methanesulfonic acid solutions at room temperature. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121518] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
|
2
|
Yue S, Muniz MC, Calegari Andrade MF, Zhang L, Car R, Panagiotopoulos AZ. When do short-range atomistic machine-learning models fall short? J Chem Phys 2021; 154:034111. [DOI: 10.1063/5.0031215] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Shuwen Yue
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | | |
Collapse
|
3
|
Zhao CL, Zhao DX, Jiang QY, Zhang HX, Li S, Yang ZZ. Polarizable TIP7P Water Model with Perturbation Charges Evaluated from ABEEM. J Phys Chem B 2020; 124:2450-2464. [PMID: 32141292 DOI: 10.1021/acs.jpcb.9b11775] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A polarizable version of the rigid seven-site (TIP7P) water model with the atom-bond electronegativity equalization method (ABEEM) is proposed. The model uses direct polarization, where an isolated water monomer in the equilibrium geometry is assumed as a reference state and the polarization of the monomer arises from interacting with other molecules as a perturbation of the reference state. The charge on each site of the monomer splits into reference charge and perturbation charge. The perturbation charge arises only because of other reference charges. The interaction of the perturbation charge with other perturbation charges is replaced using polarization scaling to enhance the interaction of perturbation charge with the reference charges of the sites from other molecules. The perturbation charges are updated by evaluating explicit expressions once. This direct polarization is time-reversible because the charge update is independent of the charges in previous simulation steps. A Slater-type damping function moderates the short-range electrostatics to treat charge diffusion. The Ewald method corrects the long-range electrostatics both in the nuclei movement and in electronegativity equalization to diminish the size effect. The water model is parameterized by fitting the ab initio results of water clusters and the experimental results of water monomers and thermodynamic properties for liquid water. Owing to polarizability, the model performs better than the TIP7P model in terms of vaporization enthalpy, isothermal compressibility, and shear viscosity of the liquid phase. It performs better at the melting point of ice but slightly worse under critical conditions than the TIP7P model. Direct polarization has a low time complexity of O(N) and is a good choice for ABEEM to improve its computational efficiency.
Collapse
Affiliation(s)
- Chong-Li Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Qian-Ying Jiang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Hai-Xia Zhang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Shenmin Li
- College of Environment and Chemical Engineering, Dalian University, Dalian 116622, People's Republic of China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| |
Collapse
|
4
|
Moučka F, Zamfir S, Bratko D, Luzar A. Molecular polarizability in open ensemble simulations of aqueous nanoconfinements under electric field. J Chem Phys 2019; 150:164702. [PMID: 31042910 DOI: 10.1063/1.5094170] [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/15/2022] Open
Abstract
Molecular polarization at aqueous interfaces involves fast degrees of freedom that are often averaged-out in atomistic-modeling approaches. The resulting effective interactions depend on a specific environment, making explicit account of molecular polarizability particularly important in solutions with pronounced anisotropic perturbations, including solid/liquid interfaces and external fields. Our work concerns polarizability effects in nanoscale confinements under electric field, open to an unperturbed bulk environment. We model aqueous molecules and ions in hydrophobic pores using the Gaussian-charge-on-spring BK3-AH representation. This involves nontrivial methodology developments in expanded ensemble Monte Carlo simulations for open systems with long-ranged multibody interactions and necessitates further improvements for efficient modeling of polarizable ions. Structural differences between fixed-charge and polarizable models were captured in molecular dynamics simulations for a set of closed systems. Our open ensemble results with the BK3 model in neat-aqueous systems capture the ∼10% reduction of molecular dipoles within the surface layer near the hydrophobic pore walls in analogy to reported quantum mechanical calculations at water/vapor interfaces. The polarizability affects the interfacial dielectric behavior and weakens the electric-field dependence of water absorption at pragmatically relevant porosities. We observe moderate changes in thermodynamic properties and atom and charged-site spatial distributions; the Gaussian distribution of mobile charges on water and ions in the polarizable model shifts the density amplitudes and blurs the charge-layering effects associated with increased ion absorption. The use of polarizable force field indicates an enhanced response of interfacial ion distributions to applied electric field, a feature potentially important for in silico modeling of electric double layer capacitors.
Collapse
Affiliation(s)
- F Moučka
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23221, USA
| | - S Zamfir
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23221, USA
| | - D Bratko
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23221, USA
| | - A Luzar
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23221, USA
| |
Collapse
|
5
|
Shafiei M, von Domaros M, Bratko D, Luzar A. Anisotropic structure and dynamics of water under static electric fields. J Chem Phys 2019; 150:074505. [DOI: 10.1063/1.5079393] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Mahdi Shafiei
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284-2006, USA
| | - Michael von Domaros
- Department of Chemistry, University of California, Irvine, Irvine, California 92697, USA
| | - Dusan Bratko
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284-2006, USA
| | - Alenka Luzar
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284-2006, USA
| |
Collapse
|
6
|
Ghahremanpour MM, van Maaren PJ, Caleman C, Hutchison GR, van der Spoel D. Polarizable Drude Model with s-Type Gaussian or Slater Charge Density for General Molecular Mechanics Force Fields. J Chem Theory Comput 2018; 14:5553-5566. [PMID: 30281307 DOI: 10.1021/acs.jctc.8b00430] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Gas-phase electric properties of molecules can be computed routinely using wave function methods or density functional theory (DFT). However, these methods remain computationally expensive for high-throughput screening of the vast chemical space of virtual compounds. Therefore, empirical force fields are a more practical choice in many cases, particularly since force field methods allow one to routinely predict the physicochemical properties in the condensed phases. This work presents Drude polarizable models, to increase the physical realism in empirical force fields, where the core particle is treated as a point charge and the Drude particle is treated either as a 1 s-Gaussian or a ns-Slater ( n = 1, 2, 3) charge density. Systematic parametrization to large high-quality quantum chemistry data obtained from the open access Alexandria Library ( https://doi.org/10.5281/zenodo.1004711 ) ensures the transferability of these parameters. The dipole moments and isotropic polarizabilities of the isolated molecules predicted by the proposed Drude models are in agreement with experiment with accuracy similar to DFT calculations at the B3LYP/aug-cc-pVTZ level of theory. The results show that the inclusion of explicit polarization into the models reduces the root-mean-square deviation with respect to DFT calculations of the predicted dipole moments of 152 dimers and clusters by more than 50%. Finally, we show that the accuracy of the electrostatic interaction energy of the water dimers can be improved systematically by the introduction of polarizable smeared charges as a model for charge penetration.
Collapse
Affiliation(s)
- Mohammad Mehdi Ghahremanpour
- Uppsala Center for Computational Chemistry, Department of Cell and Molecular Biology , Uppsala University , Husargatan 3 , Box 596, SE-75124 Uppsala , Sweden
| | - Paul J van Maaren
- Uppsala Center for Computational Chemistry, Department of Cell and Molecular Biology , Uppsala University , Husargatan 3 , Box 596, SE-75124 Uppsala , Sweden
| | - Carl Caleman
- Department of Physics and Astronomy , Uppsala University , Box 516, SE-75120 Uppsala , Sweden.,Center for Free-Electron Laser Science , Deutsches Elektronen-Synchrotron , DE-22607 Hamburg , Germany
| | - Geoffrey R Hutchison
- Department of Chemistry , University of Pittsburgh , Pittsburgh , Pennsylvania 15260 , United States
| | - David van der Spoel
- Uppsala Center for Computational Chemistry, Department of Cell and Molecular Biology , Uppsala University , Husargatan 3 , Box 596, SE-75124 Uppsala , Sweden
| |
Collapse
|
7
|
Jiang H, Moultos OA, Economou IG, Panagiotopoulos AZ. Hydrogen-Bonding Polarizable Intermolecular Potential Model for Water. J Phys Chem B 2016; 120:12358-12370. [DOI: 10.1021/acs.jpcb.6b08205] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Othonas A. Moultos
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
8
|
Tröster P, Lorenzen K, Tavan P. Polarizable six-point water models from computational and empirical optimization. J Phys Chem B 2014; 118:1589-602. [PMID: 24437570 DOI: 10.1021/jp4125765] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Tröster et al. (J. Phys. Chem B 2013, 117, 9486-9500) recently suggested a mixed computational and empirical approach to the optimization of polarizable molecular mechanics (PMM) water models. In the empirical part the parameters of Buckingham potentials are optimized by PMM molecular dynamics (MD) simulations. The computational part applies hybrid calculations, which combine the quantum mechanical description of a H2O molecule by density functional theory (DFT) with a PMM model of its liquid phase environment generated by MD. While the static dipole moments and polarizabilities of the PMM water models are fixed at the experimental gas phase values, the DFT/PMM calculations are employed to optimize the remaining electrostatic properties. These properties cover the width of a Gaussian inducible dipole positioned at the oxygen and the locations of massless negative charge points within the molecule (the positive charges are attached to the hydrogens). The authors considered the cases of one and two negative charges rendering the PMM four- and five-point models TL4P and TL5P. Here we extend their approach to three negative charges, thus suggesting the PMM six-point model TL6P. As compared to the predecessors and to other PMM models, which also exhibit partial charges at fixed positions, TL6P turned out to predict all studied properties of liquid water at p0 = 1 bar and T0 = 300 K with a remarkable accuracy. These properties cover, for instance, the diffusion constant, viscosity, isobaric heat capacity, isothermal compressibility, dielectric constant, density, and the isobaric thermal expansion coefficient. This success concurrently provides a microscopic physical explanation of corresponding shortcomings of previous models. It uniquely assigns the failures of previous models to substantial inaccuracies in the description of the higher electrostatic multipole moments of liquid phase water molecules. Resulting favorable properties concerning the transferability to other temperatures and conditions like the melting of ice are also discussed.
Collapse
Affiliation(s)
- Philipp Tröster
- Lehrstuhl für Biomolekulare Optik, Fakultät für Physik, Ludwig-Maximilians-Universität München , Oettingenstrasse 67, D-80538 Müunchen, Germany
| | | | | |
Collapse
|
9
|
Kiss PT, Baranyai A. A systematic development of a polarizable potential of water. J Chem Phys 2014; 138:204507. [PMID: 23742493 DOI: 10.1063/1.4807600] [Citation(s) in RCA: 110] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Based on extensive studies of existing potentials we propose a new molecular model for water. The new model is rigid and contains three Gaussian charges. Contrary to other models, all charges take part in the polarization of the molecule. They are connected by harmonic springs to their gas-phase positions: the negative charge to a prescribed point on the main axis of the molecule; the positive charges to the hydrogens. The mechanical equilibrium between the electrostatic forces and the spring forces determines the polarization of the molecule which is established by iteration at every timestep. The model gives excellent estimates for ambient liquid properties and reasonably good results from high-pressure solids to gas-phase clusters. We present a detailed description of the development of this model and a large number of calculated properties compared to the estimates of the nonpolarizable TIP4P∕2005 [J. L. F. Abascal and C. Vega, J. Chem. Phys. 123, 234505 (2005)], the polarizable GCPM [P. Paricaud, M. Predota, A. A. Chialvo, and P. T. Cummings, J. Chem. Phys. 122, 244511 (2005)], and our earlier BKd3 model [P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)]. The best overall performance is shown by the new model.
Collapse
Affiliation(s)
- Péter T Kiss
- Institute of Chemistry, Eötvös University, P.O. Box 32, 1518 Budapest 112, Hungary
| | | |
Collapse
|
10
|
Tröster P, Lorenzen K, Schwörer M, Tavan P. Polarizable water models from mixed computational and empirical optimization. J Phys Chem B 2013; 117:9486-500. [PMID: 23844727 DOI: 10.1021/jp404548k] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Here we suggest a mixed computational and empirical approach serving to optimize the parameters of complex and polarizable molecular mechanics (PMM) models for complicated liquids. The computational part of the parameter optimization relies on hybrid calculations combining density functional theory (DFT) for a solute molecule with a PMM treatment of its solvent environment at well-defined thermodynamic conditions. As an application we have developed PMM models for water featuring ν = 3, 4, and 5 points of force action, a Gaussian inducible dipole and a Buckingham potential at the oxygen, the experimental liquid phase geometry, the experimental gas phase polarizability α(exp)(g) = 1.47 ų, and, for ν = 4 and 5, the gas phase value μ(exp)(g) = 1.855 D for the static dipole moment. The widths of the Gaussian dipoles and, for ν = 4 and 5, also the electrostatic geometries of these so-called TLνP models are derived from self-consistent DFT/PMM calculations, and the parameters of the Buckingham potentials (and the static TL3P dipole moment) are estimated from molecular dynamics (MD) simulations. The high quality of the resulting models is demonstrated for the observables targeted during optimization (potential energy per molecule, pressure, radial distribution functions) and a series of predicted properties (quadrupole moments, density at constant pressure, dielectric constant, diffusivity, viscosity, compressibility, heat capacity) at certain standard conditions. Remaining deficiencies and possible ways for their removal are discussed.
Collapse
Affiliation(s)
- Philipp Tröster
- Lehrstuhl für Biomolekulare Optik, Fakultät für Physik, Ludwig-Maximilians-Universität München, Oettingenstrasse 67, D-80538 München, Germany
| | | | | | | |
Collapse
|
11
|
Kiss PT, Baranyai A. Testing the recent charge-on-spring type polarizable water models. II. Vapor-liquid equilibrium. J Chem Phys 2012. [DOI: 10.1063/1.4767064] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
12
|
Kiss PT, Bertsyk P, Baranyai A. Testing recent charge-on-spring type polarizable water models. I. Melting temperature and ice properties. J Chem Phys 2012. [DOI: 10.1063/1.4767063] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
|
13
|
Kiss PT, Baranyai A. Density maximum and polarizable models of water. J Chem Phys 2012; 137:084506. [DOI: 10.1063/1.4746419] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
|
14
|
Kiss PT, Baranyai A. On the pressure calculation for polarizable models in computer simulation. J Chem Phys 2012; 136:104109. [DOI: 10.1063/1.3692733] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|