1
|
Eltareb A, Lopez GE, Giovambattista N. Potential energy landscape of a flexible water model: Equation of state, configurational entropy, and Adam-Gibbs relationship. J Chem Phys 2024; 160:154510. [PMID: 38639318 DOI: 10.1063/5.0200306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Accepted: 03/28/2024] [Indexed: 04/20/2024] Open
Abstract
The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid-liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam-Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully).
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, New York 11210, USA
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, New York 10016, USA
| | - Gustavo E Lopez
- Department of Chemistry, Lehman College of the City University of New York, Bronx, New York 10468, USA
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, USA
| | - Nicolas Giovambattista
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, New York 11210, USA
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, New York 10016, USA
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, USA
| |
Collapse
|
2
|
Bushuev YG, Grosu Y, Chorążewski M. Spontaneous Dipole Reorientation in Confined Water and Its Effect on Wetting/Dewetting of Hydrophobic Nanopores. ACS Appl Mater Interfaces 2024; 16:7604-7616. [PMID: 38300737 PMCID: PMC10875646 DOI: 10.1021/acsami.3c17272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Revised: 01/10/2024] [Accepted: 01/16/2024] [Indexed: 02/03/2024]
Abstract
The properties of nanoconfined fluids are important for a broad range of natural and engineering systems. In particular, wetting/dewetting of hydrophobic nanoporous materials is crucial due to their broad applicability for molecular separation and liquid purification; energy storage, conversion, recuperation, and dissipation; for catalysis, chromatography, and so on. In this work, a rapid, orchestrated, and spontaneous dipole reorientation was observed in hydrophobic nanotubes of various pore sizes d (7.9-16.5 Å) via simulations. This phenomenon leads to the fragmentation of water clusters in the narrow nanopores (d = 7.9, 10 Å) and strongly affects dewetting through cluster repulsion. The cavitation in these pores has an electrostatic origin. The dependence of hydrogen-bonded network properties on the tube aperture is obtained and is used to explain wetting (intrusion)-dewetting (extrusion) hysteresis. Computer simulations and experimental data demonstrate that d equals ca. 12.5 Å is a threshold between a nonhysteretic (spring) behavior, where intrusion-extrusion is reversible, and a hysteretic one (shock absorber), where hysteresis is prominent. This work suggests that water clustering and the electrostatic nature of cavitation are important factors that can be effectively exploited for controlling the wetting/dewetting of nanoporous materials.
Collapse
Affiliation(s)
- Yuriy G. Bushuev
- Institute
of Chemistry, University of Silesia in Katowice, Szkolna 9 Street, 40-006 Katowice, Poland
| | - Yaroslav Grosu
- Institute
of Chemistry, University of Silesia in Katowice, Szkolna 9 Street, 40-006 Katowice, Poland
- Centre
for Cooperative Research on Alternative Energies (CIC EnergiGUNE), Basque Research and Technology Alliance (BRTA), Alava Technology Park, Albert Einstein
48, Vitoria, Gasteiz 01510, Spain
| | - Mirosław Chorążewski
- Institute
of Chemistry, University of Silesia in Katowice, Szkolna 9 Street, 40-006 Katowice, Poland
| |
Collapse
|
3
|
Zhu W, Huang Y, Zhu C, Wu HH, Wang L, Bai J, Yang J, Francisco JS, Zhao J, Yuan LF, Zeng XC. Room temperature electrofreezing of water yields a missing dense ice phase in the phase diagram. Nat Commun 2019; 10:1925. [PMID: 31028288 PMCID: PMC6486617 DOI: 10.1038/s41467-019-09950-z] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2018] [Accepted: 04/05/2019] [Indexed: 11/14/2022] Open
Abstract
Water can freeze into diverse ice polymorphs depending on the external conditions such as temperature (T) and pressure (P). Herein, molecular dynamics simulations show evidence of a high-density orthorhombic phase, termed ice χ, forming spontaneously from liquid water at room temperature under high-pressure and high external electric field. Using free-energy computations based on the Einstein molecule approach, we show that ice χ is an additional phase introduced to the state-of-the-art T–P phase diagram. The χ phase is the most stable structure in the high-pressure/low-temperature region, located between ice II and ice VI, and next to ice V exhibiting two triple points at 6.06 kbar/131.23 K and 9.45 kbar/144.24 K, respectively. A possible explanation for the missing ice phase in the T–P phase diagram is that ice χ is a rare polarized ferroelectric phase, whose nucleation/growth occurs only under very high electric fields. Water can crystallize in different ice polymorphs according to temperature and pressure conditions. Here the authors predict by molecular dynamics simulations a new ice phase spontaneously forming at room temperature under high pressure and high electric field.
Collapse
Affiliation(s)
- Weiduo Zhu
- Hefei National Laboratory for Physical Sciences at Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China.,Department of Chemistry, University of Nebraska, Lincoln, NE, 68588, USA
| | - Yingying Huang
- Department of Chemistry, University of Nebraska, Lincoln, NE, 68588, USA.,Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Ministry of Education, Dalian University of Technology, Dalian, 116024, China.,Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai, 201210, China
| | - Chongqin Zhu
- Department of Chemistry, University of Nebraska, Lincoln, NE, 68588, USA
| | - Hong-Hui Wu
- Department of Chemistry, University of Nebraska, Lincoln, NE, 68588, USA
| | - Lu Wang
- Hefei National Laboratory for Physical Sciences at Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Jaeil Bai
- Department of Chemistry, University of Nebraska, Lincoln, NE, 68588, USA
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Joseph S Francisco
- Department of Chemistry, University of Nebraska, Lincoln, NE, 68588, USA
| | - Jijun Zhao
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Ministry of Education, Dalian University of Technology, Dalian, 116024, China.
| | - Lan-Feng Yuan
- Hefei National Laboratory for Physical Sciences at Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China.
| | - Xiao Cheng Zeng
- Hefei National Laboratory for Physical Sciences at Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China. .,Department of Chemistry, University of Nebraska, Lincoln, NE, 68588, USA. .,Department of Chemical & Biomolecular Engineering and Department of Mechanical and Materials Engineering, University of Nebraska, Lincoln, NE, 68588, USA.
| |
Collapse
|
4
|
Affiliation(s)
- C. Huy Pham
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Sandeep K. Reddy
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Karl Chen
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| | - Chris Knight
- Leadership
Computing Facility, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, United States
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University of California−San Diego, La Jolla, California 92093, United States
| |
Collapse
|
5
|
Huang Y, Zhu C, Wang L, Cao X, Su Y, Jiang X, Meng S, Zhao J, Zeng XC. A new phase diagram of water under negative pressure: The rise of the lowest-density clathrate s-III. Sci Adv 2016; 2:e1501010. [PMID: 26933681 PMCID: PMC4758744 DOI: 10.1126/sciadv.1501010] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2015] [Accepted: 11/30/2015] [Indexed: 05/26/2023]
Abstract
Ice and ice clathrate are not only omnipresent across polar regions of Earth or under terrestrial oceans but also ubiquitous in the solar system such as on comets, asteroids, or icy moons of the giant planets. Depending on the surrounding environment (temperature and pressure), ice alone exhibits an exceptionally rich and complicated phase diagram with 17 known crystalline polymorphs. Water molecules also form clathrate compounds with inclusion of guest molecules, such as cubic structure I (s-I), cubic structure II (s-II), hexagonal structure H (s-H), tetragonal structure T (s-T), and tetragonal structure K (s-K). Recently, guest-free clathrate structure II (s-II), also known as ice XVI located in the negative-pressure region of the phase diagram of water, is synthesized in the laboratory and motivates scientists to reexamine other ice clathrates with low density. Using extensive Monte Carlo packing algorithm and dispersion-corrected density functional theory optimization, we predict a crystalline clathrate of cubic structure III (s-III) composed of two large icosihexahedral cavities (8(6)6(8)4(12)) and six small decahedral cavities (8(2)4(8)) per unit cell, which is dynamically stable by itself and can be fully stabilized by encapsulating an appropriate guest molecule in the large cavity. A new phase diagram of water ice with TIP4P/2005 (four-point transferable intermolecular potential/2005) model potential is constructed by considering a variety of candidate phases. The guest-free s-III clathrate with ultralow density overtakes s-II and s-H phases and emerges as the most stable ice polymorph in the pressure region below -5834 bar at 0 K and below -3411 bar at 300 K.
Collapse
Affiliation(s)
- Yingying Huang
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
- Department of Chemistry and Nebraska Center for Materials and Nanoscience, University of Nebraska, Lincoln, NE 68588, USA
| | - Chongqin Zhu
- Department of Chemistry and Nebraska Center for Materials and Nanoscience, University of Nebraska, Lincoln, NE 68588, USA
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Lu Wang
- Hefei National Laboratory for Physical Sciences at the Microscale and Collaborative Innovation Center of Chemistry for Energy Materials, Department of Materials Science and Engineering, University of Science and Technology of China, Hefei 230026, China
| | - Xiaoxiao Cao
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| | - Yan Su
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| | - Xue Jiang
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| | - Sheng Meng
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Jijun Zhao
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| | - Xiao Cheng Zeng
- Department of Chemistry and Nebraska Center for Materials and Nanoscience, University of Nebraska, Lincoln, NE 68588, USA
- Hefei National Laboratory for Physical Sciences at the Microscale and Collaborative Innovation Center of Chemistry for Energy Materials, Department of Materials Science and Engineering, University of Science and Technology of China, Hefei 230026, China
| |
Collapse
|
6
|
Lin TJ, Hsing CR, Wei CM, Kuo JL. Structure prediction of the solid forms of methanol: an ab initio random structure searching approach. Phys Chem Chem Phys 2016; 18:2736-46. [DOI: 10.1039/c5cp06583f] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Liquid methanol and methanol clusters have been comprehensively studied to reveal their local structure and hydrogen bond networks.
Collapse
Affiliation(s)
- Tzu-Jen Lin
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| | - Cheng-Rong Hsing
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| | - Ching-Ming Wei
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences
- Academic Sinica
- Taipei
- Taiwan
| |
Collapse
|
7
|
Crouse J, Loock HP, Cann NM. The photoexcitation of crystalline ice and amorphous solid water: A molecular dynamics study of outcomes at 11 K and 125 K. J Chem Phys 2015. [DOI: 10.1063/1.4926666] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- J. Crouse
- Department of Chemistry, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - H.-P. Loock
- Department of Chemistry, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - N. M. Cann
- Department of Chemistry, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
8
|
Abstract
Proton ordering in ice Ih and nucleation of ice XI domains are observed for the first time in undoped H2O and D2O; the ice Ih/II/XI triple point is also identified.
Collapse
Affiliation(s)
- Fei Yen
- Key Laboratory of Materials Physics
- Institute of Solid State Physics
- Hefei Institutes of Physical Science
- Chinese Academy of Sciences
- Hefei
| | - Zhenhua Chi
- Key Laboratory of Materials Physics
- Institute of Solid State Physics
- Hefei Institutes of Physical Science
- Chinese Academy of Sciences
- Hefei
| |
Collapse
|
9
|
Abstract
In this work, we describe a procedure to evaluate the free energy of molecular solids with the GROMACS molecular dynamics package. The free energy is calculated using the Einstein molecule method that can be regarded as a small modification of the Einstein crystal method. Here, the position and orientation of the molecules is fixed by using an Einstein field that binds with harmonic springs at least three non-collinear atoms (or points of the molecule) to their reference positions. The validity of the Einstein field is tested by performing free-energy calculations of methanol, water (ice), and patchy colloids molecular solids. The free energies calculated with GROMACS show a very good agreement with those obtained using Monte Carlo and with previously published results.
Collapse
Affiliation(s)
- J L Aragones
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | | | | | | |
Collapse
|
10
|
Harsányi I, Pusztai L. Hydration structure in concentrated aqueous lithium chloride solutions: A reverse Monte Carlo based combination of molecular dynamics simulations and diffraction data. J Chem Phys 2012. [DOI: 10.1063/1.4767437] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
|
11
|
Pinnick ER, Erramilli S, Wang F. Predicting the melting temperature of ice-Ih with only electronic structure information as input. J Chem Phys 2012; 137:014510. [PMID: 22779668 DOI: 10.1063/1.4731693] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
The melting temperature of ice-Ih was calculated with only electronic structure information as input by creating a problem-specific force field. The force field, Water model by AFM for Ice and Liquid (WAIL), was developed with the adaptive force matching (AFM) method by fitting to post-Hartree-Fock quality forces obtained in quantum mechanics∕molecular mechanics calculations. WAIL predicts the ice-Ih melting temperature to be 270 K. The model also predicts the densities of ice and water, the temperature of maximum density of water, the heat of vaporizations, and the radial distribution functions for both ice and water in good agreement with experimental measurements. The non-dissociative WAIL model is very similar to a flexible version of the popular TIP4P potential and has comparable computational cost. By customizing to problem-specific configurations with the AFM approach, the resulting model is remarkably more accurate than any variants of TIP4P for simulating ice-Ih and water in the temperature range from 253 K and 293 K under ambient pressure.
Collapse
Affiliation(s)
- Eric R Pinnick
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | | | | |
Collapse
|
12
|
Abstract
A new flexible water model, TIP4P/2005f, is developed. The idea was to add intramolecular degrees of freedom to the successful rigid model TIP4P/2005 in order to try to improve the predictions for some properties, and to enable the calculation of new ones. The new model incorporates flexibility by means of a Morse potential for the bond stretching and a harmonic term for the angle bending. The parameters have been fitted to account for the peaks of the infrared spectrum of liquid water and to produce an averaged geometry close to that of TIP4P/2005. As for the intermolecular interactions, only a small change in the σ parameter of the Lennard-Jones potential has been introduced. The overall predictions are very close to those of TIP4P/2005. This ensures that the new model may be used with the same confidence as its predecessor in studies where a flexible model is advisable.
Collapse
Affiliation(s)
- Miguel A González
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | | |
Collapse
|
13
|
Pinnick ER, Calderon CE, Rusnak AJ, Wang F. Achieving fast convergence of ab initio free energy perturbation calculations with the adaptive force-matching method. Theor Chem Acc 2012; 131. [DOI: 10.1007/s00214-012-1146-6] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
14
|
|
15
|
Kambara O, Takahashi K, Hayashi M, Kuo JL. Assessment of density functional theory to calculate the phase transition pressure of ice. Phys Chem Chem Phys 2012; 14:11484-90. [DOI: 10.1039/c2cp41495c] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
16
|
|
17
|
Noya EG, Sesé LM, Ramírez R, McBride C, Conde MM, Vega C. Path integral Monte Carlo simulations for rigid rotors and their application to water. Mol Phys 2011. [DOI: 10.1080/00268976.2010.528202] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
18
|
|
19
|
|
20
|
Best RB, Mittal J. Protein Simulations with an Optimized Water Model: Cooperative Helix Formation and Temperature-Induced Unfolded State Collapse. J Phys Chem B 2010; 114:14916-23. [DOI: 10.1021/jp108618d] [Citation(s) in RCA: 210] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Robert B. Best
- Cambridge University, Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, United Kingdom, and Lehigh University, Department of Chemical Engineering, Bethlehem, Pennsylvania 18015, United States
| | - Jeetain Mittal
- Cambridge University, Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, United Kingdom, and Lehigh University, Department of Chemical Engineering, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
21
|
Affiliation(s)
- J. L. Aragones
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040, Spain
| | - L. G. MacDowell
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040, Spain
| | - C. Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040, Spain
| |
Collapse
|
22
|
Noya EG, Vega C, Doye JPK, Louis AA. The stability of a crystal with diamond structure for patchy particles with tetrahedral symmetry. J Chem Phys 2010; 132:234511. [DOI: 10.1063/1.3454907] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Zhang J, Xiao ZR, Kuo JL. Calculation of near K edge x-ray absorption spectra and hydrogen bond network in ice XIII under compression. J Chem Phys 2010. [DOI: 10.1063/1.3421650] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Conde MM, Vega C, McBride C, Noya EG, Ramírez R, Sesé LM. Can gas hydrate structures be described using classical simulations? J Chem Phys 2010; 132:114503. [DOI: 10.1063/1.3353953] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
25
|
|
26
|
|
27
|
McBride C, Vega C, Noya EG, Ramírez R, Sesé LM. Quantum contributions in the ice phases: The path to a new empirical model for water—TIP4PQ/2005. J Chem Phys 2009; 131:024506. [DOI: 10.1063/1.3175694] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
28
|
|
29
|
Noya EG, Conde MM, Vega C. Computing the free energy of molecular solids by the Einstein molecule approach: ices XIII and XIV, hard-dumbbells and a patchy model of proteins. J Chem Phys 2009; 129:104704. [PMID: 19044935 DOI: 10.1063/1.2971188] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The recently proposed Einstein molecule approach is extended to compute the free energy of molecular solids. This method is a variant of the Einstein crystal method of Frenkel and Ladd [J. Chem. Phys. 81, 3188 (1984)]. In order to show its applicability, we have computed the free energy of a hard-dumbbell solid, of two recently discovered solid phases of water, namely, ice XIII and ice XIV, where the interactions between water molecules are described by the rigid nonpolarizable TIP4P/2005 model potential, and of several solid phases that are thermodynamically stable for an anisotropic patchy model with octahedral symmetry which mimics proteins. Our calculations show that both the Einstein crystal method and the Einstein molecule approach yield the same results within statistical uncertainty. In addition, we have studied in detail some subtle issues concerning the calculation of the free energy of molecular solids. First, for solids with noncubic symmetry, we have studied the effect of the shape of the simulation box on the free energy. Our results show that the equilibrium shape of the simulation box must be used to compute the free energy in order to avoid the appearance of artificial stress in the system that will result in an increase in the free energy. In complex solids, such as the solid phases of water, another difficulty is related to the choice of the reference structure. As in some cases there is no obvious orientation of the molecules; it is not clear how to generate the reference structure. Our results will show that, as long as the structure is not too far from the equilibrium structure, the calculated free energy is invariant to the reference structure used in the free energy calculations. Finally, the strong size dependence of the free energy of solids is also studied.
Collapse
Affiliation(s)
- E G Noya
- Departamento de Quimica-Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain.
| | | | | |
Collapse
|
30
|
Vega C, Abascal JLF, Conde MM, Aragones JL. What ice can teach us about water interactions: a critical comparison of the performance of different water models. Faraday Discuss 2009; 141:251-76; discussion 309-46. [DOI: 10.1039/b805531a] [Citation(s) in RCA: 328] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
|
31
|
Aragones JL, Conde MM, Noya EG, Vega C. The phase diagram of water at high pressures as obtained by computer simulations of the TIP4P/2005 model: the appearance of a plastic crystal phase. Phys Chem Chem Phys 2009; 11:543-55. [DOI: 10.1039/b812834k] [Citation(s) in RCA: 65] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|