1
|
Shang R, Wu T, Meguid SA. Molecular dynamics simulations of the effect of static electric field on progressive ice formation. J Chem Phys 2024; 161:094504. [PMID: 39230380 DOI: 10.1063/5.0226624] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2024] [Accepted: 08/20/2024] [Indexed: 09/05/2024] Open
Abstract
Ice accumulation under static electric fields presents a significant hazard to transmission lines and power grids. Contemporary computational studies of electrofreezing predominantly probed excessive electric fields (109 V/m) that are significantly higher than those typically encountered in proximity to transmission lines. To elucidate the influence of realistic electric fields (105 V/m) on ice crystallization, we run extensive molecular dynamics (MD) simulations across dual ice-water coexistence systems. Three aspects of work were accordingly examined. First, we investigated the influence of the effect of static electric fields, with a strength of 105 V/m, along three orthogonal axes on the phase transition during the encountered freezing and melting processes. Second, we established the mechanism of how the direction of an electric field, the initial ice crystallography, and the adjacent crystal planes influence the solidification process. Third, the results of our MD simulations were further post-processed to determine the dipole moment, radial distribution, and angle distribution resulting from the static electric field. Our results indicate that while weak electric fields do not cause complete polarization of liquid water molecules, they can induce a transition to a more structured ice-like geometry of the water molecules at the ice-water interphase region, particularly when applied perpendicular to the ice-water interphase. Notably, the interface adjacent to cubic ice exhibits a greater response to the electric fields than that adjacent to hexagonal ice. This is attributable to the intrinsic differences in their original hydrogen bonding networks.
Collapse
Affiliation(s)
- Ruiqi Shang
- Mechanics and Aerospace Design Laboratory, University of Toronto, Toronto, Ontario M5S 3G8, Canada
- The Institute of Advanced Technologies in Energy and Electrical Engineering, Shenzhen International Graduate School, Tsinghua University, Shenzhen, Guangdong, People's Republic of China
| | - Tongyu Wu
- Mechanics and Aerospace Design Laboratory, University of Toronto, Toronto, Ontario M5S 3G8, Canada
| | - S A Meguid
- Mechanics and Aerospace Design Laboratory, University of Toronto, Toronto, Ontario M5S 3G8, Canada
| |
Collapse
|
2
|
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]
|
3
|
Matsumoto RA, Thompson MW, Vuong VQ, Zhang W, Shinohara Y, van Duin ACT, Kent PRC, Irle S, Egami T, Cummings PT. Investigating the Accuracy of Water Models through the Van Hove Correlation Function. J Chem Theory Comput 2021; 17:5992-6005. [PMID: 34516134 DOI: 10.1021/acs.jctc.1c00637] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present molecular-simulation-based calculations of the Van Hove correlation function (VHF) of water using multiple modeling approaches: classical molecular dynamics with simple three-site nonpolarizable models, with a polarizable model, and with a reactive force field; density functional tight-binding molecular dynamics; and ab initio molecular dynamics. Due to the many orders of magnitude difference in the computational cost of these approaches, we investigate how small and short the simulations can be while still yielding sufficiently accurate and interpretable results for the VHF. We investigate the accuracy of the different models by comparing them to recently published inelastic X-ray scattering measurements of the VHF. We find that all of the models exhibit qualitative agreement with the experiments, and in some models and for some properties, the agreement is quantitative. This work lays the foundation for future simulation approaches to calculating the VHF for aqueous solutions in bulk and under nanoconfinement.
Collapse
Affiliation(s)
- Ray A Matsumoto
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37212, United States.,Multiscale Modeling and Simulation Center, Vanderbilt University, Nashville, Tennessee 37212, United States
| | - Matthew W Thompson
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37212, United States.,Multiscale Modeling and Simulation Center, Vanderbilt University, Nashville, Tennessee 37212, United States
| | - Van Quan Vuong
- Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Weiwei Zhang
- Department of Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, State College, Pennsylvania 16802, United States
| | - Yuya Shinohara
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Adri C T van Duin
- Department of Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, State College, Pennsylvania 16802, United States
| | - Paul R C Kent
- Center for Nanophase Materials Sciences and Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Stephan Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Takeshi Egami
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States.,Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37916, United States.,Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Peter T Cummings
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37212, United States.,Multiscale Modeling and Simulation Center, Vanderbilt University, Nashville, Tennessee 37212, United States
| |
Collapse
|
4
|
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
|
5
|
Yue S, Panagiotopoulos AZ. Dynamic properties of aqueous electrolyte solutions from non-polarisable, polarisable, and scaled-charge models. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1645901] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Shuwen Yue
- Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ, 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
|
|
8
|
Fundamental interfacial mechanisms underlying electrofreezing. Adv Colloid Interface Sci 2018; 251:26-43. [PMID: 29289337 DOI: 10.1016/j.cis.2017.12.003] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2017] [Revised: 11/19/2017] [Accepted: 12/03/2017] [Indexed: 11/24/2022]
Abstract
This article reviews the fundamental interfacial mechanisms underlying electrofreezing (promotion of ice nucleation via the application of an electric field). Electrofreezing has been an active research topic for many decades, with applications in food preservation, cryopreservation, cryogenics and ice formation. There is substantial literature detailing experimental and simulations-based studies, which aim to understand the complex mechanisms underlying accelerated ice nucleation in the presence of electric fields and electrical charge. This work provides a critical review of all such studies. It is noted that application-focused studies of electrofreezing are excluded from this review; such studies have been previously reviewed in literature. This review focuses only on fundamental studies, which analyze the physical mechanisms underlying electrofreezing. Topics reviewed include experimental studies on electrofreezing (DC and AC electric fields), pyroelectricity-based control of freezing, molecular dynamics simulations of electrofreezing, and thermodynamics-based explanations of electrofreezing. Overall, it is seen that electrofreezing can enable disruptive advancements in the control of liquid-to-solid phase change, and that our current understanding of the underlying mechanisms can be significantly improved through further studies of various interfacial effects coming into play.
Collapse
|
9
|
Székely E, Baranyai A. Local order in water: The fifth neighbor. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.02.087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
10
|
Székely E, Varga IK, Baranyai A. Tetrahedrality and hydrogen bonds in water. J Chem Phys 2016; 144:224502. [PMID: 27306013 DOI: 10.1063/1.4953555] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
We carried out extensive calculations of liquid water at different temperatures and pressures using the BK3 model suggested recently [P. T. Kiss and A. Baranyai, J. Chem. Phys. 138, 204507 (2013)]. In particular, we were interested in undercooled regions to observe the propensity of water to form tetrahedral coordination of closest neighbors around a central molecule. We compared the found tetrahedral order with the number of hydrogen bonds and with the partial pair correlation functions unfolded as distributions of the closest, the second closest, etc. neighbors. We found that contrary to the number of hydrogen bonds, tetrahedrality changes substantially with state variables. Not only the number of tetrahedral arrangements increases with lowering the pressure, the density, and the temperature but the domain size of connecting tetrahedral structures as well. The difference in tetrahedrality is very pronounced between the two sides of the Widom line and even more so between the low density amorphous (LDA) and high density amorphous (HDA) phases. We observed that in liquid water and in HDA, the 5th water molecule, contrary to ice and LDA, is positioned between the first and the second coordination shell. We found no convincing evidence of structural heterogeneity or regions referring to structural transition.
Collapse
Affiliation(s)
- Eszter Székely
- Institute of Chemistry, Eötvös University, P.O. Box 32, Budapest 112 1518, Hungary
| | - Imre K Varga
- Institute of Chemistry, Eötvös University, P.O. Box 32, Budapest 112 1518, Hungary
| | - András Baranyai
- Institute of Chemistry, Eötvös University, P.O. Box 32, Budapest 112 1518, Hungary
| |
Collapse
|
11
|
Lemkul J, Huang J, Roux B, MacKerell AD. An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications. Chem Rev 2016; 116:4983-5013. [PMID: 26815602 PMCID: PMC4865892 DOI: 10.1021/acs.chemrev.5b00505] [Citation(s) in RCA: 379] [Impact Index Per Article: 47.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2015] [Indexed: 11/28/2022]
Abstract
Molecular mechanics force fields that explicitly account for induced polarization represent the next generation of physical models for molecular dynamics simulations. Several methods exist for modeling induced polarization, and here we review the classical Drude oscillator model, in which electronic degrees of freedom are modeled by charged particles attached to the nuclei of their core atoms by harmonic springs. We describe the latest developments in Drude force field parametrization and application, primarily in the last 15 years. Emphasis is placed on the Drude-2013 polarizable force field for proteins, DNA, lipids, and carbohydrates. We discuss its parametrization protocol, development history, and recent simulations of biologically interesting systems, highlighting specific studies in which induced polarization plays a critical role in reproducing experimental observables and understanding physical behavior. As the Drude oscillator model is computationally tractable and available in a wide range of simulation packages, it is anticipated that use of these more complex physical models will lead to new and important discoveries of the physical forces driving a range of chemical and biological phenomena.
Collapse
Affiliation(s)
- Justin
A. Lemkul
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Baltimore, Maryland 21201, United States
| | - Jing Huang
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Baltimore, Maryland 21201, United States
| | - Benoît Roux
- Department
of Biochemistry and Molecular Biology, University
of Chicago, Chicago, Illinois 60637, United
States
| | - Alexander D. MacKerell
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Baltimore, Maryland 21201, United States
| |
Collapse
|
12
|
Jiang H, Moultos OA, Economou IG, Panagiotopoulos AZ. Gaussian-Charge Polarizable and Nonpolarizable Models for CO2. J Phys Chem B 2016; 120:984-94. [DOI: 10.1021/acs.jpcb.5b11701] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/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
|
13
|
Jiang H, Mester Z, Moultos OA, Economou IG, Panagiotopoulos AZ. Thermodynamic and Transport Properties of H2O + NaCl from Polarizable Force Fields. J Chem Theory Comput 2015; 11:3802-10. [DOI: 10.1021/acs.jctc.5b00421] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Zoltan Mester
- 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
|
14
|
Chialvo AA, Moucka F, Vlcek L, Nezbeda I. Vapor–Liquid Equilibrium and Polarization Behavior of the GCP Water Model: Gaussian Charge-on-Spring versus Dipole Self-Consistent Field Approaches to Induced Polarization. J Phys Chem B 2015; 119:5010-9. [DOI: 10.1021/acs.jpcb.5b00595] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ariel A. Chialvo
- Chemical
Sciences Division, Geochemistry and Interfacial Sciences Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, United States
| | - Filip Moucka
- Faculty
of Science, J. E. Purkinje University, 40096 Usti nad
Labem, Czech Republic
| | - Lukas Vlcek
- Chemical
Sciences Division, Geochemistry and Interfacial Sciences Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, United States
- Joint
Institute for Computational Sciences, Oak Ridge National Laboratory, Oak
Ridge, Tennessee 37831-6173, United States
| | - Ivo Nezbeda
- Faculty
of Science, J. E. Purkinje University, 40096 Usti nad
Labem, Czech Republic
- E.
Hala Laboratory of Thermodynamics, Institute of Chemical Process Fundamentals Academy of Sciences, 16502 Prague 6, Czech Republic
| |
Collapse
|
15
|
Vlcek L, Uhlik F, Moucka F, Nezbeda I, Chialvo AA. Thermodynamics of Small Alkali Metal Halide Cluster Ions: Comparison of Classical Molecular Simulations with Experiment and Quantum Chemistry. J Phys Chem A 2015; 119:488-500. [DOI: 10.1021/jp509401d] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lukas Vlcek
- Chemical Sciences Division, Geochemistry & Interfacial Sciences Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, United States
- Joint
Institute for Computational Sciences, Oak Ridge National Laboratory, Oak
Ridge, Tennessee 37831-6173, United States
| | - Filip Uhlik
- Department
of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, 128 43 Prague 2, Czech Republic
| | - Filip Moucka
- Faculty
of Science, J. E. Purkinje University, 400 96 Usti nad
Labem, Czech Republic
| | - Ivo Nezbeda
- Faculty
of Science, J. E. Purkinje University, 400 96 Usti nad
Labem, Czech Republic
- E.
Hala Laboratory of Thermodynamics, Institute of Chemical Process Fundamentals, Academy of Sciences, 165 02 Prague 6, Czech Republic
| | - Ariel A. Chialvo
- Chemical Sciences Division, Geochemistry & Interfacial Sciences Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, United States
| |
Collapse
|
16
|
Steele-MacInnis M, Reimer J, Bachmann S. Hydrothermal properties of the COS/D2 water model: a polarizable charge-on-spring water model, at elevated temperatures and pressures. RSC Adv 2015. [DOI: 10.1039/c5ra13495a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
Molecular simulations have been conducted to assess the pVT properties and static permittivity of the charge-on-spring polarizable water model COS/D2 at hydrothermal conditions from 300 to 450 °C and bulk densities of 0.001 to 1.0 g cm−3.
Collapse
Affiliation(s)
- Matthew Steele-MacInnis
- Department of Geosciences
- The University of Arizona
- Tucson
- USA
- Institute for Geochemistry and Petrology
| | - Joachim Reimer
- Laboratory for Bioenergy and Catalysis
- Paul Scherrer Institute
- 5232 Villigen PSI
- Switzerland
| | - Stephan Bachmann
- Laboratory of Physical Chemistry
- ETH Zürich
- CH-8093 Zürich
- Switzerland
| |
Collapse
|
17
|
Kiss PT, Sega M, Baranyai A. Efficient Handling of Gaussian Charge Distributions: An Application to Polarizable Molecular Models. J Chem Theory Comput 2014; 10:5513-9. [DOI: 10.1021/ct5009069] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Péter T. Kiss
- Institute
of Chemistry, Eötvös University, 1518 Budapest 112, P.O. Box 32, Hungary
| | - Marcello Sega
- Department
of Computational Biological Chemistry, Univiersity of Vienna, Währingerstr.
17, 1090 Vienna, Austria
| | - András Baranyai
- Institute
of Chemistry, Eötvös University, 1518 Budapest 112, P.O. Box 32, Hungary
| |
Collapse
|
18
|
Chialvo AA, Vlcek L. Ewald Summation Approach to Potential Models of Aqueous Electrolytes Involving Gaussian Charges and Induced Dipoles: Formal and Simulation Results. J Phys Chem B 2014; 118:13658-70. [DOI: 10.1021/jp509074p] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ariel A. Chialvo
- Chemical Sciences Division, Geochemistry & Interfacial Sciences Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, United States
| | - Lukas Vlcek
- Chemical Sciences Division, Geochemistry & Interfacial Sciences Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, United States
- Joint
Institute for Computational Sciences, Oak Ridge National Laboratory, Oak
Ridge, Tennessee 37831-6173, United States
| |
Collapse
|
19
|
Kiss PT, Baranyai A. A new polarizable force field for alkali and halide ions. J Chem Phys 2014; 141:114501. [DOI: 10.1063/1.4895129] [Citation(s) in RCA: 65] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Péter T. Kiss
- Institute of Chemistry, Eötvös University, 1518 Budapest 112, PO BOX 32, Hungary
| | - András Baranyai
- Institute of Chemistry, Eötvös University, 1518 Budapest 112, PO BOX 32, Hungary
| |
Collapse
|
20
|
|
21
|
Kolafa J, Viererblová L. Static Dielectric Constant from Simulations Revisited: Fluctuations or External Field? J Chem Theory Comput 2014; 10:1468-76. [DOI: 10.1021/ct500025m] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Jiří Kolafa
- Department of Physical
Chemistry, Institute of Chemical Technology, Prague, Technická
5, 166 28 Praha 6, Czech Republic
| | - Linda Viererblová
- Department of Physical
Chemistry, Institute of Chemical Technology, Prague, Technická
5, 166 28 Praha 6, Czech Republic
| |
Collapse
|
22
|
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
|
23
|
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
|
24
|
Moučka F, Nezbeda I, Smith WR. Computationally efficient Monte Carlo simulations for polarisable models: multi-particle move method for water and aqueous electrolytes. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.804183] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
25
|
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
|
26
|
Finkelstein Y, Moreh R. Proton dynamics in ice VII at high pressures. J Chem Phys 2013; 139:044716. [DOI: 10.1063/1.4816630] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
|
27
|
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
|
28
|
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
|
29
|
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
|
30
|
Kiss P, Darvas M, Baranyai A, Jedlovszky P. Surface properties of the polarizable Baranyai-Kiss water model. J Chem Phys 2012; 136:114706. [DOI: 10.1063/1.3692602] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Péter Kiss
- Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/a, H-1117 Budapest, Hungary
| | | | | | | |
Collapse
|
31
|
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
|
32
|
Baranyai A, Kiss PT. Polarizable model of water with field-dependent polarization. J Chem Phys 2011; 135:234110. [DOI: 10.1063/1.3670962] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Kiss PT, Baranyai A. Transferable model of water with variable molecular size. J Chem Phys 2011; 134:214111. [DOI: 10.1063/1.3596727] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
34
|
Kiss PT, Baranyai A. Sources of the deficiencies in the popular SPC/E and TIP3P models of water. J Chem Phys 2011; 134:054106. [DOI: 10.1063/1.3548869] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Viererblová L, Kolafa J. A classical polarizable model for simulations of water and ice. Phys Chem Chem Phys 2011; 13:19925-35. [DOI: 10.1039/c1cp22198a] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
36
|
Vega C, Abascal JLF. Simulating water with rigid non-polarizable models: a general perspective. Phys Chem Chem Phys 2011; 13:19663-88. [DOI: 10.1039/c1cp22168j] [Citation(s) in RCA: 658] [Impact Index Per Article: 50.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|