1
|
Wang X, Tse YLS. Flexible Polarizable Water Model Parameterized via Gaussian Process Regression. J Chem Theory Comput 2022; 18:7155-7165. [PMID: 36374554 DOI: 10.1021/acs.jctc.2c00529] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Water is one of the most common components in molecular dynamics (MD) simulations. Using Gaussian process regression for predicting the properties of a water model without the need of running a simulation whenever the parameters are changed, we obtained a flexible polarizable water model, named SWM4/Fw, that is able to reproduce many reference water properties. The added flexibility is critical for modeling chemical reactions in which chemical bonds can be stretched or even broken and for directly calculating vibrational spectra. In addition to being one of the few water models that are both flexible and polarizable, SWM4/Fw is also efficient thanks to the extended Lagrangian scheme with Drude oscillators. The overall accuracy is on par with or better than the related SWM4-NDP model.
Collapse
Affiliation(s)
- Xinyan Wang
- Department of Chemistry, The Chinese University of Hong Kong, Sha Tin, Hong Kong000000, China
| | - Ying-Lung Steve Tse
- Department of Chemistry, The Chinese University of Hong Kong, Sha Tin, Hong Kong000000, China
| |
Collapse
|
2
|
Yang H, Gladich I, Boucly A, Artiglia L, Ammann M. Orcinol and resorcinol induce local ordering of water molecules near the liquid-vapor interface. ENVIRONMENTAL SCIENCE: ATMOSPHERES 2022; 2:1277-1291. [PMID: 36561553 PMCID: PMC9648629 DOI: 10.1039/d2ea00015f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 08/17/2022] [Indexed: 12/25/2022]
Abstract
Resorcinol and orcinol are simple members of the family of phenolic compounds present in particulate matter in the atmosphere; they are amphiphilic in nature and thus surface active in aqueous solution. Here, we used X-ray photoelectron spectroscopy to probe the concentration of resorcinol (benzene-1,3-diol) and orcinol (5-methylbenzene-1,3-diol) at the liquid-vapor interface of aqueous solutions. Qualitatively consistent surface propensity and preferential orientation was obtained by molecular dynamics simulations. Auger electron yield near-edge X-ray absorption fine structure (NEXAFS) spectroscopy was used to probe the hydrogen bonding (HB) structure, indicating that the local structure of water molecules near the surface of the resorcinol and orcinol solutions tends towards a larger fraction of tetrahedrally coordinated molecules than observed at the liquid-vapor interface of pure water. The order parameter obtained from the molecular dynamics simulations confirm these observations. This effect is being discussed in terms of the formation of an ordered structure of these molecules at the surface leading to patterns of hydrated OH groups with distances among them that are relatively close to those in ice. These results suggest that the self-assembly of phenolic species at the aqueous solution-air interface could induce freezing similar to the case of fatty alcohol monolayers and, thus, be of relevance for ice nucleation in the atmosphere. We also attempted at looking at the changes of the O 1b1, 3a2 and 1b2 molecular orbitals of liquid water, which are known to be sensitive to the HB structure as well, in response to the presence of resorcinol and orcinol. However, these changes remained negligible within uncertainty for both experimentally obtained valence spectra and theoretically calculated density of states.
Collapse
Affiliation(s)
- Huanyu Yang
- Laboratory of Environmental Chemistry, Paul Scherrer Institut5232 VilligenSwitzerland,Institute of Atmospheric and Climate Science, ETH Zürich8092 ZürichSwitzerland
| | - Ivan Gladich
- Qatar Environment & Energy Research Institute, Hamad Bin Khalifa UniversityP.O. Box 34110DohaQatar
| | - Anthony Boucly
- Laboratory of Environmental Chemistry, Paul Scherrer Institut5232 VilligenSwitzerland,Electrochemistry Laboratory, Paul Scherrer Institut5232 VilligenSwitzerland
| | - Luca Artiglia
- Laboratory of Environmental Chemistry, Paul Scherrer Institut5232 VilligenSwitzerland,Laboratory for Catalysis and Sustainable Chemistry, Paul Scherrer Institut5232 VilligenSwitzerland
| | - Markus Ammann
- Laboratory of Environmental Chemistry, Paul Scherrer Institut5232 VilligenSwitzerland
| |
Collapse
|
3
|
Kwan V, Maiti SR, Saika-Voivod I, Consta S. Salt Enrichment and Dynamics in the Interface of Supercooled Aqueous Droplets. J Am Chem Soc 2022; 144:11148-11158. [PMID: 35715222 DOI: 10.1021/jacs.2c01159] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The interconversion reaction of NaCl between the contact-ion pair (CIP) and the solvent-separated ion pair (SSIP) as well as the free-ion state in cold droplets has not yet been investigated. We report direct computational evidence that the lower is the temperature, the closer to the surface the ion interconversion reaction takes place. In supercooled droplets the enrichment of the subsurface in salt becomes more evident. The stability of the SSIP relative to the CIP increases as the ion-pairing is transferred toward the droplet's outer layers. In the free-ion state, where the ions diffuse independently in the solution, the number density of Cl- shows a broad maximum in the interior in addition to the well-known maximum in the surface. In the study of the reaction dynamics, we find a weak coupling between the interionic NaCl distance reaction coordinate and the solvent degrees of freedom, which contrasts with the diffusive crossing of the free energy barrier found in bulk solution modeling. The H2O self-diffusion coefficient is found to be at least an order of magnitude larger than that in the bulk solution. We propose to exploit the enhanced surface ion concentration at low temperature to eliminate salts from droplets in native mass spectrometry ionization methods.
Collapse
Affiliation(s)
- Victor Kwan
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Shoubhik R Maiti
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada.,Department of Chemistry, The University of Sheffield, Sheffield S3 7HF, United Kingdom
| | - Ivan Saika-Voivod
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's A1B 3X7, Canada
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
4
|
Malek SMA, Kwan V, Saika-Voivod I, Consta S. Low Density Interior in Supercooled Aqueous Nanodroplets Expels Ions to the Subsurface. J Am Chem Soc 2021; 143:13113-13123. [PMID: 34375522 DOI: 10.1021/jacs.1c04142] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The interaction between water and ions within droplets plays a key role in the chemical reactivity of atmospheric and man-made aerosols. Here we report direct computational evidence that in supercooled aqueous nanodroplets a lower density core of tetrahedrally coordinated water expels the cosmotropic ions to the denser and more disordered subsurface. In contrast, at room temperature, depending on the nature of the ion, the radial distribution in the droplet core is nearly uniform or elevated toward the center. We analyze the spatial distribution of a single ion in terms of a reference electrostatic model. The energy of the system in the analytical model is expressed as the sum of the electrostatic and surface energy of a deformable droplet. The model predicts that the ion is subject to a harmonic potential centered at the droplet's center of mass. We name this effect "electrostatic confinement". The model's predictions are consistent with the simulation findings for a single ion at room temperature but not at supercooling. We anticipate this study to be the starting point for investigating the structure of supercooled (electro)sprayed droplets that are used to preserve the conformations of macromolecules originating from the bulk solution.
Collapse
Affiliation(s)
- Shahrazad M A Malek
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's A1B 3X7, Canada
| | - Victor Kwan
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Ivan Saika-Voivod
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's A1B 3X7, Canada.,Department of Applied Mathematics, Western University, London, Ontario N6A 3K7, Canada
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
5
|
Ye HF, Wang J, Zheng YG, Zhang HW, Chen Z. Machine learning for reparameterization of four-site water models: TIP4P-BG and TIP4P-BGT. Phys Chem Chem Phys 2021; 23:10164-10173. [PMID: 33951125 DOI: 10.1039/d0cp05831a] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Parameterizing an effective water model is a challenging issue because of the difficulty in maintaining a comprehensive balance among the diverse physical properties of water with a limited number of parameters. The advancement in machine learning provides a promising path to search for a reliable set of parameters. Based on the TIP4P water model, hence, about 6000 molecular dynamics (MD) simulations for pure water at 1 atm and in the range of 273-373 K are conducted here as the training data. The back-propagation (BP) neural network is then utilized to construct an efficient mapping between the model parameters and four crucial physical properties of water, including the density, vaporization enthalpy, self-diffusion coefficient and viscosity. Without additional time-consuming MD simulations, this mapping operation could result in sufficient and accurate data for high-population genetic algorithm (GA) to optimize the model parameters as much as possible. Based on the proposed parameterizing strategy, TIP4P-BG (a conventional four-site water model) and TIP4P-BGT (an advanced model with temperature-dependent parameters) are established. Both the water models exhibit excellent performance with a reasonable balance among the four crucial physical properties. The relevant mean absolute percentage errors are 3.53% and 3.08%, respectively. Further calculations on the temperature of maximum density, isothermal compressibility, thermal expansion coefficient, radial distribution function and surface tension are also performed and the resulting values are in good agreement with the experimental values. Through this water modeling example, the potential of the proposed data-driven machine learning procedure has been demonstrated for parameterizing a MD-based material model.
Collapse
Affiliation(s)
- Hong-Fei Ye
- International Research Center for Computational Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, Department of Engineering Mechanics, Faculty of Vehicle Engineering and Mechanics, Dalian University of Technology, Dalian 116024, P. R. China.
| | | | | | | | | |
Collapse
|
6
|
Effect of truncating electrostatic interactions on predicting thermodynamic properties of water–methanol systems. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1547824] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
7
|
Kong X, Toubin C, Habartova A, Pluharova E, Roeselova M, Pettersson JBC. Rapid Water Transport through Organic Layers on Ice. J Phys Chem A 2018; 122:4861-4868. [PMID: 29741896 DOI: 10.1021/acs.jpca.8b01951] [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/29/2022]
Abstract
Processes involving atmospheric aerosol and cloud particles are affected by condensation of organic compounds that are omnipresent in the atmosphere. On ice particles, organic compounds with hydrophilic functional groups form hydrogen bonds with the ice and orient their hydrophobic groups away from the surface. The organic layer has been expected to constitute a barrier to gas uptake, but recent experimental studies suggest that the accommodation of water molecules on ice is only weakly affected by condensed short-chain alcohol layers. Here, we employ molecular dynamics simulations to study the water interactions with n-butanol covered ice at 200 K and show that the small effect of the condensed layer is due to efficient diffusion of water molecules along the surface plane while seeking appropriate sites to penetrate, followed by penetration driven by the combined attractive forces from butanol OH groups and water molecules within the ice. The water molecules that penetrate through the n-butanol layer become strongly bonded by approximately three hydrogen bonds at the butanol-ice interface. The obtained accommodation coefficient (0.81 ± 0.03) is in excellent agreement with results from previous environmental molecular beam experiments, leading to a picture where an adsorbed n-butanol layer does not alter the apparent accommodation coefficient but dramatically changes the detailed molecular dynamics and kinetics.
Collapse
Affiliation(s)
- Xiangrui Kong
- Department of Chemistry and Molecular Biology, Atmospheric Science , University of Gothenburg , SE-41296 Gothenburg , Sweden
| | - Céline Toubin
- Laboratoire PhLAM, UMR 8523 , Université de Lille , F-59000 Lille , France
| | - Alena Habartova
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo nam. 2 , CZ-16610 Prague 6, Czech Republic
| | - Eva Pluharova
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo nam. 2 , CZ-16610 Prague 6, Czech Republic.,J. Heyrovsky Institute of Physical Chemistry of the Czech Academy of Sciences , Dolejskova 2155/3 , CZ-18223 Prague 8, Czech Republic
| | - Martina Roeselova
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo nam. 2 , CZ-16610 Prague 6, Czech Republic
| | - Jan B C Pettersson
- Department of Chemistry and Molecular Biology, Atmospheric Science , University of Gothenburg , SE-41296 Gothenburg , Sweden
| |
Collapse
|
8
|
Zhong J, Carignano MA, Kais S, Zeng XC, Francisco JS, Gladich I. Tuning the Stereoselectivity and Solvation Selectivity at Interfacial and Bulk Environments by Changing Solvent Polarity: Isomerization of Glyoxal in Different Solvent Environments. J Am Chem Soc 2018; 140:5535-5543. [DOI: 10.1021/jacs.8b01503] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jie Zhong
- Department of Chemistry, University of Nebraska—Lincoln, Lincoln, Nebraska 68588, United States
| | - Marcelo A. Carignano
- Qatar Environment and Energy Research Institute, Hamad Bin Khalifa University, P.O. Box 31110, Doha, Qatar
| | - Sabre Kais
- Department of Chemistry, Department of Physics and Birck Nanotechnology Center, Purdue University, West Lafayette, Indiana 47907, United States
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska—Lincoln, Lincoln, Nebraska 68588, United States
| | - Joseph S. Francisco
- Department of Chemistry, University of Nebraska—Lincoln, Lincoln, Nebraska 68588, United States
| | - Ivan Gladich
- Qatar Environment and Energy Research Institute, Hamad Bin Khalifa University, P.O. Box 31110, Doha, Qatar
| |
Collapse
|
9
|
Hoehn RD, Carignano MA, Kais S, Zhu C, Zhong J, Zeng XC, Francisco JS, Gladich I. Hydrogen bonding and orientation effects on the accommodation of methylamine at the air-water interface. J Chem Phys 2016; 144:214701. [DOI: 10.1063/1.4950951] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
|
10
|
Kar RK, Bhunia A. Will It Be Beneficial To Simulate the Antifreeze Proteins at Ice Freezing Condition or at Lower Temperature? J Phys Chem B 2015; 119:11485-95. [DOI: 10.1021/acs.jpcb.5b04919] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rajiv K. Kar
- Department
of Biophysics, Bose Institute, P-1/12, CIT Scheme VII M, Kolkata 700054, India
| | - Anirban Bhunia
- Department
of Biophysics, Bose Institute, P-1/12, CIT Scheme VII M, Kolkata 700054, India
- Biophysics
and Department of Chemistry, University of Michigan, 930 N. University
Avenue, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
11
|
Vereecken L, Glowacki DR, Pilling MJ. Theoretical Chemical Kinetics in Tropospheric Chemistry: Methodologies and Applications. Chem Rev 2015; 115:4063-114. [DOI: 10.1021/cr500488p] [Citation(s) in RCA: 146] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Luc Vereecken
- Max Planck Institute for Chemistry, 55128 Mainz, Germany
| | - David R. Glowacki
- PULSE
Institute and Department of Chemistry, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
- School
of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Department
of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
| | | |
Collapse
|
12
|
Gladich I, Oswald A, Bowens N, Naatz S, Rowe P, Roeselova M, Neshyba S. Mechanism of anisotropic surface self-diffusivity at the prismatic ice–vapor interface. Phys Chem Chem Phys 2015; 17:22947-58. [DOI: 10.1039/c5cp01330e] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Models of mesoscopic ice roughening require improved understanding of surface attachment kinetics.
Collapse
Affiliation(s)
- Ivan Gladich
- Qatar Environment and Energy Research Institute
- Hamad Bin Khalifa University
- Qatar Foundation
- Doha
- Qatar
| | - Amrei Oswald
- Chemistry Department
- University of Puget Sound
- Tacoma
- USA
| | | | - Sam Naatz
- Chemistry Department
- University of Puget Sound
- Tacoma
- USA
| | - Penny Rowe
- Department of Physics
- Universidad de Santiago de Chile
- Santiago
- Chile
- Geography Department
| | - Martina Roeselova
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- 16610 Prague 6
- Czech Republic
| | | |
Collapse
|
13
|
Kumar H, Dasgupta C, Maiti PK. Structure, dynamics and thermodynamics of single-file water under confinement: effects of polarizability of water molecules. RSC Adv 2015. [DOI: 10.1039/c4ra08730e] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Various structural, dynamic and thermodynamic properties of water molecules confined in single-wall carbon nanotubes are investigated using both polarizable and non-polarizable water models.
Collapse
Affiliation(s)
- Hemant Kumar
- Centre for Condensed Matter Theory
- Indian Institute of Science
- Bangalore-560012
- India
| | - Chandan Dasgupta
- Centre for Condensed Matter Theory
- Indian Institute of Science
- Bangalore-560012
- India
| | - Prabal K. Maiti
- Centre for Condensed Matter Theory
- Indian Institute of Science
- Bangalore-560012
- India
| |
Collapse
|
14
|
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
|
15
|
Kamath G, Deshmukh SA, Sankaranarayanan SKRS. Comparison of select polarizable and non-polarizable water models in predicting solvation dynamics of water confined between MgO slabs. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2013; 25:305003. [PMID: 23819970 DOI: 10.1088/0953-8984/25/30/305003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
We present a molecular dynamics simulation study in which we compare and contrast the performance of a polarizable shell water potential model and non-polarizable water force field-extended simple point charge (SPC/EF) model in predicting the solvation dynamics of confined water molecules sandwiched between MgO(100) slabs. Structural features based on radial distribution functions, atomic density profiles, adsorption patterns, orientational ordering and dynamical correlations such as diffusional characteristics, hydrogen bonding lifetimes and residence probabilities are used as metrics for comparison. The simulations yield significant ordering of water molecules in the two layers adjacent to the oxide interface and the extent of ordering decreases with increasing distance from the oxide-water interface. These results elucidate that the dependence of local ordering and solvation dynamics on the molecular geometry and charge distribution, observed for typical three- and four-site water models, is generally lost for confined water if polarization is explicitly included. While the interfacial water structure predicted by the polarizable and non-polarizable models are similar, the confinement and interface proximity effects on the solvation dynamics are seen to be more pronounced for polarizable water models in comparison to non-polarizable ones. The study also shows that the polarizable water model over predicts the orientational order and under predicts the transport properties of confined water. In addition, analysis of the orientational preferences and hydrogen bonding characteristics of water near oxide interfaces suggests a higher degree of tetrahedral disorder in the polarizable shell compared to the non-polarizable SPC/E flexible model. The origin of the differences in solvation behavior of confined water between oxide slabs is analyzed based on the energetic contributions of the dispersive and electrostatic terms in the two force fields. Our findings suggest some new considerations regarding the role of polarization terms in predicting confinement and interface proximity effects that may guide future development of reliable polarizable water models for confined liquids.
Collapse
Affiliation(s)
- Ganesh Kamath
- Department of Chemistry, University of Missouri-Columbia, Columbia 65211, USA
| | | | | |
Collapse
|
16
|
Baker CM, Best RB. Matching of additive and polarizable force fields for multiscale condensed phase simulations. J Chem Theory Comput 2013; 9:2826-2837. [PMID: 23997691 PMCID: PMC3752912 DOI: 10.1021/ct400116g] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Inclusion of electronic polarization effects is one of the key aspects in which the accuracy of current biomolecular force fields may be improved. The principal drawback of such approaches is the computational cost, which typically ranges from 3 - 10 times that of the equivalent additive model, and may be greater for more sophisticated treatments of polarization or other many-body effects. Here, we present a multiscale approach which may be used to enhance the sampling in simulations with polarizable models, by using the additive model as a tool to explore configuration space. We use a method based on information theory to determine the charges for an additive model that has optimal overlap with the polarizable one, and we demonstrate the feasibility of enhancing sampling via a hybrid replica exchange scheme for several model systems. An additional advantage is that, in the process, we obtain a systematic method for deriving charges for an additive model that will be the natural complement to its polarizable parent. The additive charges are found by an effective coarse-graining of the polarizable force field, rather than by ad hoc procedures.
Collapse
Affiliation(s)
- Christopher M. Baker
- University of Cambridge, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK
| | - Robert B. Best
- University of Cambridge, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892-0520, U.S.A
| |
Collapse
|
17
|
|
18
|
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
|
19
|
Reinhardt A, Doye JPK, Noya EG, Vega C. Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water. J Chem Phys 2012. [DOI: 10.1063/1.4766362] [Citation(s) in RCA: 81] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|