1
|
Holzapfel WB, Klotz S. Thermophysical properties of H2O and D2O ice Ih with contributions from proton disorder, quenching, relaxation, and extended defects: A model case for solids with quenching and relaxation. J Chem Phys 2024; 160:154508. [PMID: 38634494 DOI: 10.1063/5.0203614] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 04/01/2024] [Indexed: 04/19/2024] Open
Abstract
Application of the coherent thermodynamic model [W. Holzapfel and S. Klotz, J. Chem. Phys. 155, 024506 (2021)] for H2O ice Ih to the more detailed data for D2O ice Ih provides better insight into the contributions from quenched proton disorder and offers a new basis for understanding the apparent differences between the data for thermal expansion measured with neutron diffraction on polycrystalline samples [A. Fortes, Acta Crystallogr., Sect. B: Struct. Sci., Cryst. Eng. Mater. 74, 196 (2018) and A. Fortes, Phys. Chem. Chem. Phys 21, 8264 (2019)] and macroscopic dilatation measurements on single crystals [D. Buckingham et al., Phys. Rev. Lett. 121, 185505 (2018)]. The comparison points to contributions from defects effecting the two techniques in different ways. The uncertainties in thermodynamic data due to the contributions from proton disorder and additional defects are compared with the "reference data" [R. Feistel and W. Wagner, J. Phys. Chem. Ref. Data 35, 1021 (2006)] for H2O ice Ih.
Collapse
Affiliation(s)
- W B Holzapfel
- Department Physik, Universität Paderborn, Paderborn D-33095, Germany
| | - S Klotz
- IMPMC, CNRS UMR 7590, Sorbonne Université, 4 Place Jussieu, Paris F-75252, France
| |
Collapse
|
2
|
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 PMCID: PMC11184974 DOI: 10.1063/5.0200306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [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
- Authors to whom correspondence should be addressed: ; ; and
| | | | | |
Collapse
|
3
|
Eltareb A, Lopez GE, Giovambattista N. A continuum of amorphous ices between low-density and high-density amorphous ice. Commun Chem 2024; 7:36. [PMID: 38378859 PMCID: PMC10879119 DOI: 10.1038/s42004-024-01117-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 02/01/2024] [Indexed: 02/22/2024] Open
Abstract
Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY, 11210, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| | - Gustavo E Lopez
- Department of Chemistry, Lehman College of the City University of New York, Bronx, NY, 10468, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| | - Nicolas Giovambattista
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY, 11210, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| |
Collapse
|
4
|
Hai Y, Jiang M, Tian H, Zhong G, Li W, Yang C, Chen X, Lin H. Superconductivity Above 100 K Predicted in Carbon-Cage Network. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2023; 10:e2303639. [PMID: 37807820 PMCID: PMC10667821 DOI: 10.1002/advs.202303639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2023] [Revised: 08/22/2023] [Indexed: 10/10/2023]
Abstract
To explore carbide superconductors with higher transition temperature, two novel carbon structures of cage-network are designed and their superconductivity is studied by doping metals. MC6 and MC10 are respectively identified as C24 and C32 cage-network structures. This study finds that both carbon structures drive strong electron-phonon interaction and can exhibit superconductivity above liquid nitrogen temperature. Importantly, the superconducting transition temperatures above 100 K are predicted to be achieved in C24 -cage-network systems doped by Na, Mg, Al, In, and Tl at ambient pressure, which is far higher than those in graphite, fullerene, and other carbides. Meanwhile, the superconductivity of cage-network carbides is also found to be sensitive to the electronegativity and concentration of dopant M. The result indicates that the higher transition temperatures can be obtained by optimizing the carbon-cage-network structures and the doping conditions. The study suggests that the carbon-cage-network structure is a direction to explore high-temperature superconducting carbides.
Collapse
Affiliation(s)
- Yu‐Long Hai
- Shenzhen Institute of Advanced TechnologyChinese Academy of SciencesShenzhen518055China
- Nano Science and Technology InstituteUniversity of Science and Technology of ChinaSuzhou215123China
| | - Meng‐Jing Jiang
- Shenzhen Institute of Advanced TechnologyChinese Academy of SciencesShenzhen518055China
- Nano Science and Technology InstituteUniversity of Science and Technology of ChinaSuzhou215123China
| | - Hui‐Li Tian
- Shenzhen Institute of Advanced TechnologyChinese Academy of SciencesShenzhen518055China
- Nano Science and Technology InstituteUniversity of Science and Technology of ChinaSuzhou215123China
| | - Guo‐Hua Zhong
- Shenzhen Institute of Advanced TechnologyChinese Academy of SciencesShenzhen518055China
- University of Chinese Academy of SciencesBeijing100049China
| | - Wen‐Jie Li
- Shenzhen Institute of Advanced TechnologyChinese Academy of SciencesShenzhen518055China
- University of Chinese Academy of SciencesBeijing100049China
| | - Chun‐Lei Yang
- Shenzhen Institute of Advanced TechnologyChinese Academy of SciencesShenzhen518055China
- University of Chinese Academy of SciencesBeijing100049China
| | - Xiao‐Jia Chen
- School of ScienceHarbin Institute of TechnologyShenzhen518055China
- Center for High Pressure Science and Technology Advanced ResearchShanghai201203China
| | - Hai‐Qing Lin
- School of PhysicsZhejiang UniversityHangzhou310058China
| |
Collapse
|
5
|
Wang J, Ghosh DB, Zhang Z. Computational Materials Design for Ceramic Nuclear Waste Forms Using Machine Learning, First-Principles Calculations, and Kinetics Rate Theory. MATERIALS (BASEL, SWITZERLAND) 2023; 16:4985. [PMID: 37512262 PMCID: PMC10383080 DOI: 10.3390/ma16144985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 07/05/2023] [Accepted: 07/10/2023] [Indexed: 07/30/2023]
Abstract
Ceramic waste forms are designed to immobilize radionuclides for permanent disposal in geological repositories. One of the principal criteria for the effective incorporation of waste elements is their compatibility with the host material. In terms of performance under environmental conditions, the resistance of the waste forms to degradation over long periods of time is a critical concern when they are exposed to natural environments. Due to their unique crystallographic features and behavior in nature environment as exemplified by their natural analogues, ceramic waste forms are capable of incorporating problematic nuclear waste elements while showing promising chemical durability in aqueous environments. Recent studies of apatite- and hollandite-structured waste forms demonstrated an approach that can predict the compositions of ceramic waste forms and their long-term dissolution rate by a combination of computational techniques including machine learning, first-principles thermodynamics calculations, and modeling using kinetic rate equations based on critical laboratory experiments. By integrating the predictions of elemental incorporation and degradation kinetics in a holistic framework, the approach could be promising for the design of advanced ceramic waste forms with optimized incorporation capacity and environmental degradation performance. Such an approach could provide a path for accelerated ceramic waste form development and performance prediction for problematic nuclear waste elements.
Collapse
Affiliation(s)
- Jianwei Wang
- Department of Geology and Geophysics, Center for Computation and Technology, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Dipta B Ghosh
- Department of Geology and Geophysics, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Zelong Zhang
- Department of Geology and Geophysics, Louisiana State University, Baton Rouge, LA 70803, USA
| |
Collapse
|
6
|
Eltareb A, Lopez GE, Giovambattista N. The Importance of Nuclear Quantum Effects on the Thermodynamic and Structural Properties of Low-Density Amorphous Ice: A Comparison with Hexagonal Ice. J Phys Chem B 2023; 127:4633-4645. [PMID: 37178124 DOI: 10.1021/acs.jpcb.3c01025] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We study the nuclear quantum effects (NQE) on the thermodynamic properties of low-density amorphous ice (LDA) and hexagonal ice (Ih) at P = 0.1 MPa and T ≥ 25 K. Our results are based on path-integral molecular dynamics (PIMD) and classical MD simulations of H2O and D2O using the q-TIP4P/F water model. We show that the inclusion of NQE is necessary to reproduce the experimental properties of LDA and ice Ih. While MD simulations (no NQE) predict that the density ρ(T) of LDA and ice Ih increases monotonically upon cooling, PIMD simulations indicate the presence of a density maximum in LDA and ice Ih. MD and PIMD simulations also predict a qualitatively different T-dependence for the thermal expansion coefficient αP(T) and bulk modulus B(T) of both LDA and ice Ih. Remarkably, the ρ(T), αP(T), and B(T) of LDA are practically identical to those of ice Ih. The origin of the observed NQE is due to the delocalization of the H atoms, which is identical in LDA and ice Ih. H atoms delocalize considerably (over a distance ≈ 20-25% of the OH covalent-bond length) and anisotropically (preferentially perpendicular to the OH covalent bond), leading to less linear hydrogen bonds HB (larger HOO angles and longer OO separations) than observed in classical MD simulations.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, New York 11210, United States
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Gustavo E Lopez
- Department of Chemistry, Lehman College of the City University of New York, Bronx, New York 10468, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Nicolas Giovambattista
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, New York 11210, United States
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, New York 10016, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| |
Collapse
|
7
|
Rasti S, Jónsson EÖ, Jónsson H, Meyer J. New Insights into the Volume Isotope Effect of Ice Ih from Polarizable Many-Body Potentials. J Phys Chem Lett 2022; 13:11831-11836. [PMID: 36520035 PMCID: PMC9791686 DOI: 10.1021/acs.jpclett.2c03212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2022] [Accepted: 12/06/2022] [Indexed: 06/17/2023]
Abstract
The anomalous volume isotope effect (VIE) of ice Ih is calculated and analyzed based on the quasi-harmonic approximation to account for nuclear quantum effects in the Helmholtz free energy. While a lot of recently developed polarizable many-body potential functions give a normal VIE contrary to experimental results, we find that one of them, MB-pol, yields the anomalous VIE in good agreement with the most recent high-resolution neutron diffraction measurements─better than DFT calculations. The short-range three-body terms in the MB-pol function, which are fitted to CCSD(T) calculations, are found to have a surprisingly large influence. A vibrational mode group decomposition of the zero-point pressure together with a hitherto unconsidered benchmark value for the intramolecular stretching modes of H2O ice Ih obtained from Raman spectroscopy data unveils the reason for the VIE: a delicate competition between the latter and the librations.
Collapse
Affiliation(s)
- Soroush Rasti
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands
| | - Elvar Örn Jónsson
- Science
Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| | - Hannes Jónsson
- Science
Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107Reykjavík, Iceland
| | - Jörg Meyer
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands
| |
Collapse
|
8
|
Molecular dynamics analysis of elastic properties and new phase formation during amorphous ices transformations. Sci Rep 2022; 12:13325. [PMID: 35922440 PMCID: PMC9349219 DOI: 10.1038/s41598-022-17666-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 07/28/2022] [Indexed: 12/02/2022] Open
Abstract
Unlike conventional first-order phase transitions, the kinetics of amorphous-amorphous transitions has been much less studied. The ultrasonic experiments on the transformations between low-density and high-density amorphous ice induced by pressure or heating provided the pressure and temperature dependencies of elastic moduli. In this article, we make an attempt to build a microscopic picture of these experimentally studied transformations using the molecular dynamics method with the TIP4P/Ice water model. We study carefully the dependence of the results of elastic constants calculations on the deformation rates. The system size effects are considered as well. The comparison with the experimental data enriches our understanding of the transitions observed. Our modeling gives new information about the formation mechanisms of new phase clusters during the transition between low-density and high-density amorphous ices. We analyse the applicability of the term “nucleation” for these processes.
Collapse
|
9
|
Tanaka H, Yagasaki T, Matsumoto M. On the role of intermolecular vibrational motions for ice polymorphs. III. Mode characteristics associated with negative thermal expansion. J Chem Phys 2021; 155:214502. [PMID: 34879657 DOI: 10.1063/5.0068560] [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/14/2022] Open
Abstract
Low-pressure ice forms, such as hexagonal and cubic ice, expand on cooling below temperature 60 K. This negative thermal expansivity has been explored in terms of phonon frequency modulation with varying volume and attributed to the negative Grüneisen parameters unique mostly to tetrahedrally coordinated substances. However, an underlying mechanism for the negative Grüneisen parameters has not been known except some schematic analyses. We investigate in this study the characteristics of the intermolecular vibrational modes whose Grüneisen parameters are negative by examining the individual vibrational modes rigorously. It is found that the low frequency modes below 100 cm-1, which we explicitly show are mostly bending motions of three hydrogen-bonded molecules, necessarily accompany elongation of the hydrogen bond length at peak amplitudes compared with that at the equilibrium position in executing the vibrational motions. The elongation gives rise to a decrease in the repulsive interaction while an increase in the Coulombic one. The decrease in the repulsive interaction is relaxed substantially by expansion due to its steep slope against molecular separation compared with the sluggish increase in the Coulombic one, and therefore, the negative Grüneisen parameters are obtainable. This scenario is tested against some variants of cubic ice with various water potential models. It is demonstrated that four interaction-site models are suitable to describe the intermolecular vibrations and the thermal expansivity because of the moderate tendency to favor the tetrahedral coordination.
Collapse
Affiliation(s)
- Hideki Tanaka
- Toyota Physical and Chemical Research Institute, Nagakute 480-1192, Japan
| | - Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
10
|
Cherubini M, Monacelli L, Mauri F. The microscopic origin of the anomalous isotopic properties of ice relies on the strong quantum anharmonic regime of atomic vibration. J Chem Phys 2021; 155:184502. [PMID: 34773945 DOI: 10.1063/5.0062689] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Water ice is a unique material presenting intriguing physical properties, such as negative thermal expansion and anomalous volume isotope effect (VIE). They arise from the interplay between weak hydrogen bonds and nuclear quantum fluctuations, making theoretical calculations challenging. Here, we employ the stochastic self-consistent harmonic approximation to investigate how thermal and quantum fluctuations affect the physical properties of ice XI with ab initio accuracy. Regarding the anomalous VIE, our work reveals that quantum effects on hydrogen are so strong to be in a nonlinear regime: When progressively increasing the mass of hydrogen from protium to infinity (classical limit), the volume first expands and then contracts, with a maximum slightly above the mass of tritium. We observe an anharmonic renormalization of about 10% in the bending and stretching phonon frequencies probed in IR and Raman experiments. For the first time, we report an accurate comparison of the low-energy phonon dispersion with the experimental data, possible only thanks to high-level accuracy in the electronic correlation and nuclear quantum and thermal fluctuations, paving the way for the study of thermal transport in ice from first-principles and the simulation of ice under pressure.
Collapse
Affiliation(s)
- Marco Cherubini
- Dipartimento di Fisica, Università di Roma Sapienza, Piazzale Aldo Moro 5, I-00185 Roma, Italy
| | - Lorenzo Monacelli
- Dipartimento di Fisica, Università di Roma Sapienza, Piazzale Aldo Moro 5, I-00185 Roma, Italy
| | - Francesco Mauri
- Dipartimento di Fisica, Università di Roma Sapienza, Piazzale Aldo Moro 5, I-00185 Roma, Italy
| |
Collapse
|
11
|
Holzapfel WB, Klotz S. Coherent thermodynamic model for ice Ih-A model case for complex behavior. J Chem Phys 2021; 155:024506. [PMID: 34266244 DOI: 10.1063/5.0049215] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
New data on the variation of the thermal expansion of ice Ih with temperature at ambient pressure together with new evaluations of the bulk modulus and earlier data for the heat capacity provide the basis for a coherent thermodynamic modeling of the main thermophysical properties of ice Ih over its whole range of stability. The quasi-harmonic approximation with one Debye term and seven Einstein terms, together with explicit anharmonicity, represents the dominant contribution next to minor "anomalies" from hydrogen ordering and lattice defects. The model accurately fits the main features of all experimental data and provides a basis for the comparison with earlier determinations of the phonon density of states and the Grüneisen parameters.
Collapse
Affiliation(s)
| | - Stefan Klotz
- IMPMC, CNRS UMR 7590, Sorbonne Université, 4 Place Jussieu, F-75252 Paris, France
| |
Collapse
|
12
|
Rossi M. Progress and challenges in ab initio simulations of quantum nuclei in weakly bonded systems. J Chem Phys 2021; 154:170902. [PMID: 34241065 DOI: 10.1063/5.0042572] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Atomistic simulations based on the first-principles of quantum mechanics are reaching unprecedented length scales. This progress is due to the growth in computational power allied with the development of new methodologies that allow the treatment of electrons and nuclei as quantum particles. In the realm of materials science, where the quest for desirable emergent properties relies increasingly on soft weakly bonded materials, such methods have become indispensable. In this Perspective, an overview of simulation methods that are applicable for large system sizes and that can capture the quantum nature of electrons and nuclei in the adiabatic approximation is given. In addition, the remaining challenges are discussed, especially regarding the inclusion of nuclear quantum effects (NQEs) beyond a harmonic or perturbative treatment, the impact of NQEs on electronic properties of weakly bonded systems, and how different first-principles potential energy surfaces can change the impact of NQEs on the atomic structure and dynamics of weakly bonded systems.
Collapse
Affiliation(s)
- Mariana Rossi
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
13
|
Eltareb A, Lopez GE, Giovambattista N. Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water. Phys Chem Chem Phys 2021; 23:6914-6928. [PMID: 33729222 DOI: 10.1039/d0cp04325g] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We perform path-integral molecular dynamics (PIMD) simulations of H2O and D2O using the q-TIP4P/F model. Simulations are performed at P = 1 bar and over a wide range of temperatures that include the equilibrium (T≥ 273 K) and supercooled (210 ≤T < 273 K) liquid states of water. The densities of both H2O and D2O calculated from PIMD simulations are in excellent agreement with experiments in the equilibrium and supercooled regimes. We also evaluate important thermodynamic response functions, specifically, the thermal expansion coefficient αP(T), isothermal compressibility κT(T), isobaric heat capacity CP(T), and static dielectric constant ε(T). While these properties are in excellent [αP(T) and κT(T)] or semi-quantitative agreement [CP(T) and ε(T)] with experiments in the equilibrium regime, they are increasingly underestimated upon further cooling. It follows that the inclusion of nuclear quantum effects in PIMD simulations of (q-TIP4P/F) water is not sufficient to reproduce the anomalous large fluctuations in density, entropy, and electric dipole moment characteristic of supercooled water. It has been hypothesized that water may exhibit a liquid-liquid critical point (LLCP) in the supercooled regime at P > 1 bar and that such a LLCP generates a maximum in CP(T) and κT(T) at 1 bar. Consistent with this hypothesis and in particular, with experiments, we find a maximum in the κT(T) of q-TIP4P/F light and heavy water at T≈ 230-235 K. No maximum in CP(T) could be detected down to T≥ 210 K. We also calculate the diffusion coefficient D(T) of H2O and D2O using the ring-polymer molecular dynamics (RPMD) technique and find that computer simulations are in remarkable good agreement with experiments at all temperatures studied. The results from RPMD/PIMD simulations are also compared with the corresponding results obtained from classical MD simulations of q-TIP4P/F water where atoms are represented by single interacting sites. Surprisingly, we find minor differences in most of the properties studied, with CP(T), D(T), and structural properties being the only (expected) exceptions.
Collapse
Affiliation(s)
- Ali Eltareb
- Department of Physics, Brooklyn College of the City University of New York, Brooklyn, NY 11210, USA.
| | | | | |
Collapse
|
14
|
Engel EA. Identification of synthesisable crystalline phases of water – a prototype for the challenges of computational materials design. CrystEngComm 2021. [DOI: 10.1039/d0ce01260b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We discuss the identification of experimentally realisable crystalline phases of water to outline and contextualise some of the diverse building blocks of a computational materials design process.
Collapse
Affiliation(s)
- Edgar A. Engel
- TCM Group
- Cavendish Laboratory
- University of Cambridge
- Cambridge CB3 0HE
- UK
| |
Collapse
|
15
|
Fuentes-Landete V, Rasti S, Schlögl R, Meyer J, Loerting T. Calorimetric Signature of Deuterated Ice II: Turning an Endotherm to an Exotherm. J Phys Chem Lett 2020; 11:8268-8274. [PMID: 32902994 PMCID: PMC7528406 DOI: 10.1021/acs.jpclett.0c02368] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2020] [Accepted: 09/09/2020] [Indexed: 06/11/2023]
Abstract
Calorimetric studies on ice II reveal a surprising H2O/D2O isotope effect. While the ice II to ice Ic transition is endothermic for H2O, it is exothermic for D2O samples. The transition enthalpies are +40 and -140 J/mol, respectively, where such a sign change upon isotope substitution is unprecedented in ice research. To understand the observations we employ force field calculations using two water models known to perform well for H2O ice phases and their vibrational properties. These simulations reveal that the isotope effect can be traced back to zero-point energy. q-TIP4P/F fares better and is able to account for approximately three-fourths of the isotope effect, while MB-pol only catches approximately one-third. Phonon and configurational entropy contributions are necessary to predict reasonable transition enthalpies, but they do not have an impact on the isotope effect. We suggest to use these calorimetric isotope data as a benchmark for water models.
Collapse
Affiliation(s)
- Violeta Fuentes-Landete
- Institute
of Physical Chemistry, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria
- Max
Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Soroush Rasti
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Robert Schlögl
- Max
Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Jörg Meyer
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Thomas Loerting
- Institute
of Physical Chemistry, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria
| |
Collapse
|
16
|
Tanaka H, Yagasaki T, Matsumoto M. On the role of intermolecular vibrational motions for ice polymorphs I: Volumetric properties of crystalline and amorphous ices. J Chem Phys 2019; 151:114501. [PMID: 31542026 DOI: 10.1063/1.5119748] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Intermolecular vibrations and volumetric properties are investigated using the quasiharmonic approximation with the TIP4P/2005, TIP4P/Ice, and SPC/E potential models for most of the known crystalline and amorphous ice forms that have hydrogen-disordering. The ice forms examined here cover low pressure ices (hexagonal and cubic ice I, XVI, and hypothetical dtc ice), medium pressure ices (III, IV, V, VI, XII, hydrogen-disordered variant of ice II), and high pressure ice (VII) as well as the low density and the high density amorphous forms. We focus on the thermal expansivities and the isothermal compressibilities in the low temperature regime over a wide range of pressures calculated via the intermolecular vibrational free energies. Negative thermal expansivity appears only in the low pressure ice forms. The sign of the thermal expansivity is elucidated in terms of the mode Grüneisen parameters of the low frequency intermolecular vibrational motions. Although the band structure for the low frequency region of the vibrational density of state in the medium pressure ice has a close resemblance to that in the low pressure ice, its response against volume variation is opposite. We reveal that the mixing of translational and rotational motions in the low frequency modes plays a crucial role in the appearance of the negative thermal expansivity in the low pressure ice forms. The medium pressure ices can be further divided into two groups in terms of the hydrogen-bond network flexibility, which is manifested in the properties on the molecular rearrangement against volume variation, notably the isothermal compressibility.
Collapse
Affiliation(s)
- Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
17
|
Rasti S, Meyer J. Importance of zero-point energy for crystalline ice phases: A comparison of force fields and density functional theory. J Chem Phys 2019; 150:234504. [DOI: 10.1063/1.5097021] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Soroush Rasti
- Gorlaeus Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Gorlaeus Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
18
|
Herrero CP, Ramírez R. Nuclear quantum effects in graphene bilayers. J Chem Phys 2019; 150:204707. [DOI: 10.1063/1.5096602] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Affiliation(s)
- Carlos P. Herrero
- Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain
| | - Rafael Ramírez
- Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
19
|
Cheng B, Engel EA, Behler J, Dellago C, Ceriotti M. Ab initio thermodynamics of liquid and solid water. Proc Natl Acad Sci U S A 2019; 116:1110-1115. [PMID: 30610171 PMCID: PMC6347673 DOI: 10.1073/pnas.1815117116] [Citation(s) in RCA: 152] [Impact Index Per Article: 30.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Thermodynamic properties of liquid water as well as hexagonal (Ih) and cubic (Ic) ice are predicted based on density functional theory at the hybrid-functional level, rigorously taking into account quantum nuclear motion, anharmonic fluctuations, and proton disorder. This is made possible by combining advanced free-energy methods and state-of-the-art machine-learning techniques. The ab initio description leads to structural properties in excellent agreement with experiments and reliable estimates of the melting points of light and heavy water. We observe that nuclear-quantum effects contribute a crucial [Formula: see text] to the stability of ice Ih, making it more stable than ice Ic. Our computational approach is general and transferable, providing a comprehensive framework for quantitative predictions of ab initio thermodynamic properties using machine-learning potentials as an intermediate step.
Collapse
Affiliation(s)
- Bingqing Cheng
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland;
| | - Edgar A Engel
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Universität Göttingen, 37073 Göttingen, Germany
| | | | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
20
|
Shaykomalova ES, Zhuravlyov YN. Structure and Properties of Ice Phase States. J STRUCT CHEM+ 2019. [DOI: 10.1134/s0022476619010086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
21
|
Abraham NS, Shirts MR. Thermal Gradient Approach for the Quasi-harmonic Approximation and Its Application to Improved Treatment of Anisotropic Expansion. J Chem Theory Comput 2018; 14:5904-5919. [PMID: 30281302 DOI: 10.1021/acs.jctc.8b00460] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a novel approach to efficiently implement thermal expansion in the quasi-harmonic approximation (QHA) for both isotropic and more importantly, anisotropic expansion. In this approach, we rapidly determine a crystal's equilibrium volume and shape at a given temperature by integrating along the gradient of expansion from 0 Kelvin up to the desired temperature. We compare our approach to previous isotropic methods that rely on a brute-force grid search to determine the free energy minimum, which is infeasible to carry out for anisotropic expansion, as well as quasi-anisotropic approaches that take into account the contributions to anisotropic expansion from the lattice energy. We compare these methods for experimentally known polymorphs of piracetam and resorcinol and show that both isotropic methods agree to within error up to 300 K. Using the Grüneisen parameter causes up to 0.04 kcal/mol deviation in the Gibbs free energy, but for polymorph free energy differences there is a cancellation in error with all isotropic methods within 0.025 kcal/mol at 300 K. Anisotropic expansion allows the crystals to relax into lattice geometries 0.01-0.23 kcal/mol lower in energy at 300 K relative to isotropic expansion. For polymorph free energy differences all QHA methods produced results within 0.02 kcal/mol of each other for resorcinol and 0.12 kcal/mol for piracetam, the two molecules tested here, demonstrating a cancellation of error for isotropic methods. We also find that with expansion in more than a single volume variable, there is a non-negligible rate of failure of the basic approximations of QHA. Specifically, while expanding into new harmonic modes as the box vectors are increased, the system often falls into alternate, structurally distinct harmonic modes unrelated by continuous deformation from the original harmonic mode.
Collapse
Affiliation(s)
- Nathan S Abraham
- Department of Chemical and Biological Engineering , University of Colorado Boulder , Boulder , Colorado 80309 , United States
| | - Michael R Shirts
- Department of Chemical and Biological Engineering , University of Colorado Boulder , Boulder , Colorado 80309 , United States
| |
Collapse
|
22
|
Herrero CP, Ramírez R. Thermal properties of graphene from path-integral simulations. J Chem Phys 2018; 148:102302. [PMID: 29544269 DOI: 10.1063/1.4997178] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Carlos P. Herrero
- Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain
| | - Rafael Ramírez
- Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
23
|
Pamuk B, Allen PB, Fernández-Serra MV. Insights into the Structure of Liquid Water from Nuclear Quantum Effects on the Density and Compressibility of Ice Polymorphs. J Phys Chem B 2018; 122:5694-5706. [DOI: 10.1021/acs.jpcb.8b00110] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Betül Pamuk
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
- Physics and Astronomy Department, Stony Brook University, Stony Brook, New York 11794-3800, United States
| | - P. B. Allen
- Physics and Astronomy Department, Stony Brook University, Stony Brook, New York 11794-3800, United States
| | - M.-V. Fernández-Serra
- Physics and Astronomy Department, Stony Brook University, Stony Brook, New York 11794-3800, United States
| |
Collapse
|
24
|
Herrero CP, Ramírez R. Path-integral simulation of graphene monolayers under tensile stress. Phys Chem Chem Phys 2017; 19:31898-31909. [DOI: 10.1039/c7cp06821b] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Finite-temperature properties of graphene monolayers under tensile stress have been studied by path-integral molecular dynamics (PIMD) simulations.
Collapse
Affiliation(s)
- Carlos P. Herrero
- Instituto de Ciencia de Materiales de Madrid
- Consejo Superior de Investigaciones Científicas (CSIC)
- 28049 Madrid
- Spain
| | - Rafael Ramírez
- Instituto de Ciencia de Materiales de Madrid
- Consejo Superior de Investigaciones Científicas (CSIC)
- 28049 Madrid
- Spain
| |
Collapse
|
25
|
Herrero CP, Ramírez R. Quantum effects in graphene monolayers: Path-integral simulations. J Chem Phys 2016; 145:224701. [DOI: 10.1063/1.4971453] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Carlos P. Herrero
- Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain
| | - Rafael Ramírez
- Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
26
|
Heit YN, Beran GJO. How important is thermal expansion for predicting molecular crystal structures and thermochemistry at finite temperatures? ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2016; 72:514-529. [PMID: 27484373 DOI: 10.1107/s2052520616005382] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/16/2016] [Accepted: 03/30/2016] [Indexed: 06/06/2023]
Abstract
Molecular crystals expand appreciably upon heating due to both zero-point and thermal vibrational motion, yet this expansion is often neglected in molecular crystal modeling studies. Here, a quasi-harmonic approximation is coupled with fragment-based hybrid many-body interaction calculations to predict thermal expansion and finite-temperature thermochemical properties in crystalline carbon dioxide, ice Ih, acetic acid and imidazole. Fragment-based second-order Möller-Plesset perturbation theory (MP2) and coupled cluster theory with singles, doubles and perturbative triples [CCSD(T)] predict the thermal expansion and the temperature dependence of the enthalpies, entropies and Gibbs free energies of sublimation in good agreement with experiment. The errors introduced by neglecting thermal expansion in the enthalpy and entropy cancel somewhat in the Gibbs free energy. The resulting ∼ 1-2 kJ mol(-1) errors in the free energy near room temperature are comparable to or smaller than the errors expected from the electronic structure treatment, but they may be sufficiently large to affect free-energy rankings among energetically close polymorphs.
Collapse
Affiliation(s)
- Yonaton N Heit
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
27
|
Salim MA, Willow SY, Hirata S. Ice Ih anomalies: Thermal contraction, anomalous volume isotope effect, and pressure-induced amorphization. J Chem Phys 2016; 144:204503. [DOI: 10.1063/1.4951687] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael A. Salim
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - Soohaeng Yoo Willow
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| |
Collapse
|
28
|
|
29
|
Parmentier A, Shephard JJ, Romanelli G, Senesi R, Salzmann CG, Andreani C. Evolution of Hydrogen Dynamics in Amorphous Ice with Density. J Phys Chem Lett 2015; 6:2038-2042. [PMID: 26266499 DOI: 10.1021/acs.jpclett.5b00711] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The single-particle dynamics of hydrogen atoms in several of the amorphous ices are reported using a combination of deep inelastic neutron scattering (DINS) and inelastic neutron scattering (INS). The mean kinetic energies of the hydrogen nuclei are found to increase with increasing density, indicating the weakening of hydrogen bonds as well as a trend toward steeper and more harmonic hydrogen vibrational potential energy surfaces. DINS shows much more pronounced changes in the O-H stretching component of the mean kinetic energy going from low- to high-density amorphous ices than indicated by INS and Raman spectroscopy. This highlights the power of the DINS technique to retrieve accurate ground-state kinetic energies beyond the harmonic approximation. In a novel approach, we use information from DINS and INS to determine the anharmonicity constants of the O-H stretching modes. Furthermore, our experimental kinetic energies will serve as important benchmark values for path-integral Monte Carlo simulations.
Collapse
Affiliation(s)
- A Parmentier
- †Dipartimento di Fisica and NAST Centre, Università degli Studi di Roma Tor Vergata, Via della Ricerca Scientifica 1, 00133 Roma, Italy
| | - J J Shephard
- ‡Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- §Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - G Romanelli
- †Dipartimento di Fisica and NAST Centre, Università degli Studi di Roma Tor Vergata, Via della Ricerca Scientifica 1, 00133 Roma, Italy
| | - R Senesi
- †Dipartimento di Fisica and NAST Centre, Università degli Studi di Roma Tor Vergata, Via della Ricerca Scientifica 1, 00133 Roma, Italy
- ∥CNR-IPCF Sezione di Messina, Viale F. Stagno D'Alcontres 37, 98158 Messina, Italy
| | - C G Salzmann
- ‡Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - C Andreani
- †Dipartimento di Fisica and NAST Centre, Università degli Studi di Roma Tor Vergata, Via della Ricerca Scientifica 1, 00133 Roma, Italy
- ∥CNR-IPCF Sezione di Messina, Viale F. Stagno D'Alcontres 37, 98158 Messina, Italy
| |
Collapse
|
30
|
Moustafa SG, Schultz AJ, Kofke DA. Effects of Finite Size and Proton Disorder on Lattice-Dynamics Estimates of the Free Energy of Clathrate Hydrates. Ind Eng Chem Res 2014. [DOI: 10.1021/ie504008h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sabry G. Moustafa
- Department
of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| | - Andrew J. Schultz
- Department
of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| | - David A. Kofke
- Department
of Chemical and
Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, United States
| |
Collapse
|
31
|
Conde MM, Gonzalez MA, Abascal JLF, Vega C. Determining the phase diagram of water from direct coexistence simulations: The phase diagram of the TIP4P/2005 model revisited. J Chem Phys 2013; 139:154505. [DOI: 10.1063/1.4824627] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
32
|
Ramírez R, Neuerburg N, Herrero CP. The phase diagram of ice: A quasi-harmonic study based on a flexible water model. J Chem Phys 2013; 139:084503. [DOI: 10.1063/1.4818875] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Herrero CP, Ramírez R. Topological characterization of crystalline ice structures from coordination sequences. Phys Chem Chem Phys 2013; 15:16676-85. [PMID: 23986009 DOI: 10.1039/c3cp52167b] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Affiliation(s)
- Carlos P Herrero
- Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain.
| | | |
Collapse
|
34
|
Ramírez R, Neuerburg N, Herrero CP. The phase diagram of ice Ih, II, and III: A quasi-harmonic study. J Chem Phys 2012; 137:134503. [DOI: 10.1063/1.4757064] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Herrero CP, Ramírez R. High-density amorphous ice: A path-integral simulation. J Chem Phys 2012; 137:104505. [DOI: 10.1063/1.4750027] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
|