1
|
Lin M, Xiong Z, Cao H. Bridging classical nucleation theory and molecular dynamics simulation for homogeneous ice nucleation. J Chem Phys 2024; 161:084504. [PMID: 39206829 DOI: 10.1063/5.0216645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2024] [Accepted: 08/14/2024] [Indexed: 09/04/2024] Open
Abstract
Water freezing, initiated by ice nucleation, occurs widely in nature, ranging from cellular to global phenomena. Ice nucleation has been experimentally proven to require the formation of a critical ice nucleus, consistent with classical nucleation theory (CNT). However, the accuracy of CNT quantitative predictions of critical cluster sizes and nucleation rates has never been verified experimentally. In this study, we circumvent this difficulty by using molecular dynamics (MD) simulation. The physical properties of water/ice for CNT predictions, including density, chemical potential difference, and diffusion coefficient, are independently obtained using MD simulation, whereas the calculation of interfacial free energy is based on thermodynamic assumptions of CNT, including capillarity approximation among others. The CNT predictions are compared to the MD evaluations of brute-force simulations and forward flux sampling methods. We find that the CNT and MD predicted critical cluster sizes are consistent, and the CNT predicted nucleation rates are higher than the MD predicted values within three orders of magnitude. We also find that the ice crystallized from supercooled water is stacking-disordered ice with a stacking of cubic and hexagonal ices in four representative types of stacking. The prediction discrepancies in nucleation rate mainly arise from the stacking-disordered ice structure, the asphericity of ice cluster, the uncertainty of ice-water interfacial free energy, and the kinetic attachment rate. Our study establishes a relation between CNT and MD to predict homogeneous ice nucleation.
Collapse
Affiliation(s)
- Min Lin
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhewen Xiong
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| | - Haishan Cao
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
2
|
Wang J, Hei H, Zheng Y, Zhang H, Ye H. Five-Site Water Models for Ice and Liquid Water Generated by a Series-Parallel Machine Learning Strategy. J Chem Theory Comput 2024. [PMID: 39133036 DOI: 10.1021/acs.jctc.4c00440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/13/2024]
Abstract
Icing, a common natural phenomenon, always originates from a molecule. Molecular simulation is crucial for understanding the relevant process but still faces a great challenge in obtaining a uniform and accurate description of ice and liquid water with limited model parameters. Here, we propose a series-parallel machine learning (ML) approach consisting of a classification back-propagation neural network (BPNN), parallel regression BPNNs, and a genetic algorithm to establish conventional TIP5P-BG and temperature-dependent TIP5P-BGT models. The established water models exhibit a comprehensive balance among the crucial physical properties (melting point, density, vaporization enthalpy, self-diffusion coefficient, and viscosity) with mean absolute percentage errors of 2.65 and 2.40%, respectively, and excellent predictive performance on the related properties of liquid water. For ice, the simulation results on the critical nucleus size and growth rate are in good accordance with experiments. This work offers a powerful molecular model for phase transition and icing in nanoconfinement and a construction strategy for a complex molecular model in the extreme case.
Collapse
Affiliation(s)
- Jian Wang
- International Research Center for Computational Mechanics, State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, Department of Engineering Mechanics, Faculty of Vehicle Engineering and Mechanics, Dalian University of Technology, Dalian 116024, P. R. China
| | - Haitao Hei
- International Research Center for Computational Mechanics, State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, Department of Engineering Mechanics, Faculty of Vehicle Engineering and Mechanics, Dalian University of Technology, Dalian 116024, P. R. China
| | - Yonggang Zheng
- International Research Center for Computational Mechanics, State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, Department of Engineering Mechanics, Faculty of Vehicle Engineering and Mechanics, Dalian University of Technology, Dalian 116024, P. R. China
- DUT-BSU Joint Institute, Dalian University of Technology, Dalian 116024, P. R. China
| | - Hongwu Zhang
- International Research Center for Computational Mechanics, State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, Department of Engineering Mechanics, Faculty of Vehicle Engineering and Mechanics, Dalian University of Technology, Dalian 116024, P. R. China
| | - Hongfei Ye
- International Research Center for Computational Mechanics, State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, Department of Engineering Mechanics, Faculty of Vehicle Engineering and Mechanics, Dalian University of Technology, Dalian 116024, P. R. China
| |
Collapse
|
3
|
Hayton JA, Davies MB, Whale TF, Michaelides A, Cox SJ. The limit of macroscopic homogeneous ice nucleation at the nanoscale. Faraday Discuss 2024; 249:210-228. [PMID: 37791990 DOI: 10.1039/d3fd00099k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
Nucleation in small volumes of water has garnered renewed interest due to the relevance of pore condensation and freezing under conditions of low partial pressures of water, such as in the upper troposphere. Molecular simulations can in principle provide insight on this process at the molecular scale that is challenging to achieve experimentally. However, there are discrepancies in the literature as to whether the rate in confined systems is enhanced or suppressed relative to bulk water at the same temperature and pressure. In this study, we investigate the extent to which the size of the critical nucleus and the rate at which it grows in thin films of water are affected by the thickness of the film. Our results suggest that nucleation remains bulk-like in films that are barely large enough accommodate a critical nucleus. This conclusion seems robust to the presence of physical confining boundaries. We also discuss the difficulties in unambiguously determining homogeneous nucleation rates in nanoscale systems, owing to the challenges in defining the volume. Our results suggest any impact on a film's thickness on the rate is largely inconsequential for present day experiments.
Collapse
Affiliation(s)
- John A Hayton
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | - Michael B Davies
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK
| | - Thomas F Whale
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK
- School of Earth and Environment, University of Leeds, Leeds, UK
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| |
Collapse
|
4
|
Zhu X, Vandamme M, Jiang Z, Brochard L. Molecular simulation of the confined crystallization of ice in cement nanopore. J Chem Phys 2023; 159:154704. [PMID: 37850696 DOI: 10.1063/5.0169783] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 09/26/2023] [Indexed: 10/19/2023] Open
Abstract
Freezing of water under nanoconfinement exhibits physical peculiarities with respect to the bulk water. However, experimental observations are extremely challenging at this scale, which limits our understanding of the effect of confinement on water properties upon freezing. In this study, we use molecular dynamic simulations to investigate how confinement affects the kinetics of growth of ice and the thermodynamic equilibrium of ice-liquid coexistence. TIP4P/Ice water model and CSH-FF model were applied to simulate ice crystallization in a confined cement system at temperatures down to 220 K. We adapted an interface detection algorithm and reparameterized the CHILL/CHILL+ algorithm to capture ice growth. The confinement leads to a shift of the maximum growth rate of ice to a higher temperature than for bulk water. Both the confinement and surface impurities contribute to slowing down the ice growth. For the ice-liquid coexistence at equilibrium, we derive a formulation of Thomson's equation adapted to statistical physics quantities accessible by molecular simulation, and we show that this adapted equation predicts accurately the melting line of bulk and confined ice Ih as a function of pressure. The confinement decreases systematically the melting temperature of ice of about 5 K compared with bulk ice Ih. A premelted water film about 1 nm thick is observed between the solid wall and ice, and its thickness is found to decrease continuously as temperature is lowered. We note that the surface impurities are key to the formation of the premelted water nanofilm when the temperature is lower than 250 K.
Collapse
Affiliation(s)
- Xinping Zhu
- Key Laboratory of Advanced Civil Engineering Materials of Ministry of Education, School of Materials Science and Engineering, Tongji University, Shanghai 201804, China
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| | - Matthieu Vandamme
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| | - Zhengwu Jiang
- Key Laboratory of Advanced Civil Engineering Materials of Ministry of Education, School of Materials Science and Engineering, Tongji University, Shanghai 201804, China
| | - Laurent Brochard
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| |
Collapse
|
5
|
Spesyvyi A, Žabka J, Polášek M, Charvat A, Schmidt J, Postberg F, Abel B. Charged Ice Particle Beams with Selected Narrow Mass and Kinetic Energy Distributions. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2023; 34:878-892. [PMID: 37018538 DOI: 10.1021/jasms.2c00357] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Small ice particles play an important role in atmospheric and extraterrestrial chemistry. Circumplanetary ice particles that are encountered by space probes at hypervelocities play a critical role in the determination of surface and subsurface properties of their source bodies. Here we present an apparatus for the generation of low-intensity beams of single mass-selected charged ice particles under vacuum. They are produced via electrospray ionization of water at atmospheric pressure and undergo evaporative cooling when transferred to vacuum through an atmospheric vacuum interface. m/z selection is achieved through two subsequent quadrupole mass filters operated in the variable-frequency mode within a range of m/z values between 8 × 104 and 3 × 107. Velocity and charge of the selected particles are measured using a nondestructive single-pass image charge detector. From the known electrostatic acceleration potentials and settings of the quadrupoles the particle masses could be obtained and be accurately controlled. It has been shown that the droplets are frozen within the transit time of the apparatus such that ice particles are present after the quadrupole stages and finally detected. The demonstrated correspondence between particle mass and specific quadrupole potentials in this device allows preparation of beams of single particles with a repetition rate between 0.1 and 1 Hz with various diameter distributions from 50 to 1000 nm at 30-250 eV of kinetic energy per charge. This corresponds to velocities and particle masses quickly available between 600 m/s (80 nm) and 50 m/s (900 nm) and particle charge numbers (positive) between 103 and 104[e], depending upon size.
Collapse
Affiliation(s)
- Anatolii Spesyvyi
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences, Dolejškova 2155/3, 18223 Prague 8, Czech Republic
| | - Ján Žabka
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences, Dolejškova 2155/3, 18223 Prague 8, Czech Republic
| | - Miroslav Polášek
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences, Dolejškova 2155/3, 18223 Prague 8, Czech Republic
| | - Ales Charvat
- Institute of Chemical Technology and Wilhelm-Ostwald-Institute of Physical and Theoretical Chemistry, Linnestrasse 3, D-04103 Leipzig, Germany
- Leibniz Institute of Surface Engineering, Permoserstrasse 15, D-04318 Leipzig, Germany
| | - Jürgen Schmidt
- Institute of Geological Sciences, Freie Universität Berlin, Malteserstraße 74-100, D-12249 Berlin, Germany
| | - Frank Postberg
- Institute of Geological Sciences, Freie Universität Berlin, Malteserstraße 74-100, D-12249 Berlin, Germany
| | - Bernd Abel
- Institute of Chemical Technology and Wilhelm-Ostwald-Institute of Physical and Theoretical Chemistry, Linnestrasse 3, D-04103 Leipzig, Germany
- Leibniz Institute of Surface Engineering, Permoserstrasse 15, D-04318 Leipzig, Germany
| |
Collapse
|
6
|
Montero de Hijes P, R Espinosa J, Vega C, Dellago C. Minimum in the pressure dependence of the interfacial free energy between ice Ih and water. J Chem Phys 2023; 158:124503. [PMID: 37003785 DOI: 10.1063/5.0140814] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023] Open
Abstract
Despite the importance of ice nucleation, this process has been barely explored at negative pressures. Here, we study homogeneous ice nucleation in stretched water by means of molecular dynamics seeding simulations using the TIP4P/Ice model. We observe that the critical nucleus size, interfacial free energy, free energy barrier, and nucleation rate barely change between isobars from -2600 to 500 bars when they are represented as a function of supercooling. This allows us to identify universal empirical expressions for homogeneous ice nucleation in the pressure range from -2600 to 500 bars. We show that this universal behavior arises from the pressure dependence of the interfacial free energy, which we compute by means of the mold integration technique, finding a shallow minimum around -2000 bars. Likewise, we show that the change in the interfacial free energy with pressure is proportional to the excess entropy and the slope of the melting line, exhibiting in the latter a reentrant behavior also at the same negative pressure. Finally, we estimate the excess internal energy and the excess entropy of the ice Ih-water interface.
Collapse
Affiliation(s)
| | - J R Espinosa
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Dellago
- Faculty of Physics, University of Vienna, A-1090 Vienna, Austria
| |
Collapse
|
7
|
Grabowska J, Blazquez S, Sanz E, Noya EG, Zeron IM, Algaba J, Miguez JM, Blas FJ, Vega C. Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations. J Chem Phys 2023; 158:114505. [PMID: 36948790 DOI: 10.1063/5.0132681] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023] Open
Abstract
In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log10 (J/(m3 s)) = -7(5).
Collapse
Affiliation(s)
- J Grabowska
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E G Noya
- Instituto de Química Física Rocasolano, CSIC, C/ Serrano 119, 28006 Madrid, Spain
| | - I M Zeron
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - J M Miguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
8
|
Goy C, Caupin F, Caresana M, Cremonesi L, Kalinin A, Grübel G, Potenza MAC, Grisenti RE. Refractive Index of Supercooled Water Down to 230.3 K in the Wavelength Range between 534 and 675 nm. J Phys Chem Lett 2022; 13:11872-11877. [PMID: 36520590 DOI: 10.1021/acs.jpclett.2c02803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Knowledge of the refractive index of water in the deeply supercooled metastable liquid state is important, for example, for an accurate description of optical reflection and refraction processes occurring in clouds. However, a measurement of both the temperature and wavelength dependence of the refractive index under such extreme conditions is challenging. Here, we employ Raman spectroscopy in combination with microscopic water jets in vacuum to obtain the refractive index of supercooled water to a lowest temperature of 230.3 K. While our approach is based on the analysis of Mie resonances in Raman spectra measured by using a single excitation wavelength at 532 nm, it allows us to obtain the refractive index in a wide visible wavelength range from 534 to 675 nm. Because of a direct link between the refractive index and density of water, our results provide a promising approach to help improve our understanding of water's anomalous behavior.
Collapse
Affiliation(s)
- Claudia Goy
- Centre for Molecular Water Science (CMWS), Deutsches Elektronen-Synchrotron (DESY), 22607Hamburg, Germany
| | - Frédéric Caupin
- Institut Lumiére Matiére, Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Universitaire de France, 69622Villeurbanne, France
| | - Michele Caresana
- Institut für Kernphysik, J. W. Goethe-Universität Frankfurt am Main, 60438Frankfurt am Main, Germany
| | - Llorenç Cremonesi
- Dipartimento di Fisica, Università degli Studi di Milano, 20133Milano, Italy
- DISAT-Dipartimento di Scienze dell'Ambiente e della Terra, Università di Milano-Bicocca, 20126Milano, Italy
| | - Anton Kalinin
- GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291Darmstadt, Germany
| | - Gerhard Grübel
- Centre for Molecular Water Science (CMWS), Deutsches Elektronen-Synchrotron (DESY), 22607Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, 22761Hamburg, Germany
| | - Marco A C Potenza
- Dipartimento di Fisica, Università degli Studi di Milano, 20133Milano, Italy
| | - Robert E Grisenti
- Institut für Kernphysik, J. W. Goethe-Universität Frankfurt am Main, 60438Frankfurt am Main, Germany
- GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291Darmstadt, Germany
| |
Collapse
|
9
|
Zhang W, Liu H, Fu H, Shao X, Cai W. Revealing the Mechanism of Irreversible Binding of Antifreeze Glycoproteins to Ice. J Phys Chem B 2022; 126:10637-10645. [PMID: 36513495 DOI: 10.1021/acs.jpcb.2c06183] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Antifreeze glycoproteins (AFGPs) are a special kind of antifreeze proteins with strong flexibility. Whether their antifreeze activity is achieved by reversibly or irreversibly binding to ice is widely debated, and the molecular mechanism of irreversible binding remains unclear. In this work, the antifreeze mechanism of the smallest AFGP isoform, AFGP8, is investigated at the atomic level. The results indicate that AFGP8 can bind to ice both reversibly through its hydrophobic methyl groups (peptide binding) and irreversibly through its hydrophilic disaccharide moieties (saccharide binding). Although peptide binding occurs faster than saccharide binding, free-energy calculations indicate that the latter is energetically more favorable. In saccharide binding, at least one disaccharide moiety is frozen in the grown ice, resulting in irreversible binding, while the other moieties significantly perturb the water hydrogen-bonding network, thus inhibiting ice growth more effectively. The present study reveals the coexistence of reversible and irreversible bindings of AFGP8, both contributing to the inhibition of ice growth and further provides molecular mechanism of irreversible binding.
Collapse
Affiliation(s)
- Weijia Zhang
- Research Center for Analytical Sciences, Haihe Laboratory of Sustainable Chemical Transformations, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin300071, China
| | - Han Liu
- Research Center for Analytical Sciences, Haihe Laboratory of Sustainable Chemical Transformations, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin300071, China
| | - Haohao Fu
- Research Center for Analytical Sciences, Haihe Laboratory of Sustainable Chemical Transformations, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin300071, China
| | - Xueguang Shao
- Research Center for Analytical Sciences, Haihe Laboratory of Sustainable Chemical Transformations, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin300071, China
| | - Wensheng Cai
- Research Center for Analytical Sciences, Haihe Laboratory of Sustainable Chemical Transformations, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, College of Chemistry, Nankai University, Tianjin300071, China
| |
Collapse
|
10
|
Grabowska J, Blazquez S, Sanz E, Zerón IM, Algaba J, Míguez JM, Blas FJ, Vega C. Solubility of Methane in Water: Some Useful Results for Hydrate Nucleation. J Phys Chem B 2022; 126:8553-8570. [PMID: 36222501 DOI: 10.1021/acs.jpcb.2c04867] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this paper, the solubility of methane in water along the 400 bar isobar is determined by computer simulations using the TIP4P/Ice force field for water and a simple LJ model for methane. In particular, the solubility of methane in water when in contact with the gas phase and the solubility of methane in water when in contact with the hydrate has been determined. The solubility of methane in a gas-liquid system decreases as temperature increases. The solubility of methane in a hydrate-liquid system increases with temperature. The two curves intersect at a certain temperature that determines the triple point T3 at a certain pressure. We also determined T3 by the three-phase direct coexistence method. The results of both methods agree, and we suggest 295(2) K as the value of T3 for this system. We also analyzed the impact of curvature on the solubility of methane in water. We found that the presence of curvature increases the solubility in both the gas-liquid and hydrate-liquid systems. The change in chemical potential for the formation of hydrate is evaluated along the isobar using two different thermodynamic routes, obtaining good agreement between them. It is shown that the driving force for hydrate nucleation under experimental conditions is higher than that for the formation of pure ice when compared at the same supercooling. We also show that supersaturation (i.e., concentrations above those of the planar interface) increases the driving force for nucleation dramatically. The effect of bubbles can be equivalent to that of an additional supercooling of about 20 K. Having highly supersaturated homogeneous solutions makes possible the spontaneous formation of the hydrate at temperatures as high as 285 K (i.e., 10K below T3). The crucial role of the concentration of methane for hydrate formation is clearly revealed. Nucleation of the hydrate can be either impossible or easy and fast depending on the concentration of methane which seems to play the leading role in the understanding of the kinetics of hydrate formation.
Collapse
Affiliation(s)
- Joanna Grabowska
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain.,Department of Physical Chemistry, Faculty of Chemistry and BioTechMed Center, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Samuel Blazquez
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Iván M Zerón
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Jesús Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - José Manuel Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Felipe J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Carlos Vega
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
11
|
Phan AD. Screening and collective effects in randomly pinned fluids: a new theoretical framework. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:435101. [PMID: 35985315 DOI: 10.1088/1361-648x/ac8b51] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Accepted: 08/19/2022] [Indexed: 06/15/2023]
Abstract
We propose a theoretical framework for the dynamics of bulk isotropic hard-sphere systems in the presence of randomly pinned particles and apply this theory to supercooled water to validate it. Structural relaxation is mainly governed by local and non-local activated process. As the pinned fraction grows, a local caging constraint becomes stronger and the long range collective aspect of relaxation is screened by immobile obstacles. Different responses of the local and cooperative motions results in subtle predictions for how the alpha relaxation time varies with pinning and density. Our theoretical analysis for the relaxation time of water with pinned molecules quantitatively well describe previous simulations. In addition, the thermal dependence of relaxation for unpinned bulk water is also consistent with prior computational and experimental data.
Collapse
Affiliation(s)
- Anh D Phan
- Faculty of Materials Science and Engineering, Phenikaa Institute for Advanced Study, Phenikaa University, Hanoi 12116, Vietnam
| |
Collapse
|
12
|
Abstract
Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.
Collapse
|
13
|
Sanchez-Burgos I, Tejedor AR, Vega C, Conde MM, Sanz E, Ramirez J, Espinosa JR. Homogeneous ice nucleation rates for mW and TIP4P/ICE models through Lattice Mold calculations. J Chem Phys 2022; 157:094503. [DOI: 10.1063/5.0101383] [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
Water freezing is the most common liquid-to-crystal phase transition on Earth, however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240K. We employ two different water models, mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic non-polarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih, Ic and a stacking mixture of ice Ih/Ic; reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner which contributes to benchmarking the freezing behaviour of two popular water models.
Collapse
Affiliation(s)
| | | | - Carlos Vega
- Departamento de Quimica Fisica, Universidad Complutense de Madrid Facultad de Ciencias Químicas, Spain
| | - Maria M. Conde
- Universidad Politécnica de Madrid Escuela Técnica Superior de Ingenieros Industriales, Spain
| | | | - Jorge Ramirez
- Chemical Engineering, Universidad Politécnica de Madrid Escuela Técnica Superior de Ingenieros Industriales, Spain
| | | |
Collapse
|
14
|
Pereyra RG, Sebastianelli P, Ávila EE. Homogeneous nucleation in supercooled liquid water. Determination of ice germ size and activation energy barrier in Molecular Dynamics simulations. MOLECULAR SIMULATION 2022. [DOI: 10.1080/08927022.2022.2068801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Rodolfo G. Pereyra
- Fa.M.A.F., Universidad Nacional de Córdoba, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
- IFEG-CONICET, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
| | - Paolo Sebastianelli
- Fa.M.A.F., Universidad Nacional de Córdoba, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
- School of Chemistry, University of New South Wales, Sydney, Australia
| | - Eldo E. Ávila
- Fa.M.A.F., Universidad Nacional de Córdoba, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
- IFEG-CONICET, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
| |
Collapse
|
15
|
de Almeida Ribeiro I, de Koning M, Molinero V. Is It Possible to Follow the Structural Evolution of Water in "No-Man's Land" Using a Pulsed-Heating Procedure? J Phys Chem Lett 2022; 13:1085-1089. [PMID: 35080178 DOI: 10.1021/acs.jpclett.1c04106] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The anomalous increase in compressibility and heat capacity of supercooled water has been attributed to its structural transformation of into a four-coordinated liquid. Experiments revealed that κT and Cp peak at TWthermo ≈ 229 K [Kim et al. Science 2017, 358, 1589; Pathak et al. Proc. Natl. Acad. Sci. 2021, 118, e2018379118]. Recently, a pulsed heating procedure (PHP) was employed to interrogate the structure of water, reporting a steep increase in tetrahedrality around TWPHP = 210 ± 3 K [Kringle et al. Science 2020, 369, 1490]. This discrepancy questions whether water structure and thermodynamics are decoupled, or if the shift in TW is an artifact of PHP. Here we implement PHP in molecular simulations. We find that the stationary states captured at the bottom of the pulse are not representative of the thermalized liquid or its inherent structure. Our analysis reveals a temperature-dependent distortion that shifts TWPHP to ∼20 K below TWthermo. We conclude that 2 orders of magnitude faster rates are required to sample water's inherent structure with PHP.
Collapse
Affiliation(s)
- Ingrid de Almeida Ribeiro
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Maurice de Koning
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861, Campinas, São Paulo, Brazil
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
16
|
Nikiforidis VM, Datta S, Borg MK, Pillai R. Impact of surface nanostructure and wettability on interfacial ice physics. J Chem Phys 2021; 155:234307. [PMID: 34937379 DOI: 10.1063/5.0069896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface-ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations.
Collapse
Affiliation(s)
- Vasileios-Martin Nikiforidis
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Saikat Datta
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Matthew K Borg
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Rohit Pillai
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| |
Collapse
|
17
|
Moore FJ, Royall CP, Liverpool TB, Russo J. Crystallisation and polymorph selection in active Brownian particles. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:121. [PMID: 34580776 PMCID: PMC8476478 DOI: 10.1140/epje/s10189-021-00108-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Accepted: 07/20/2021] [Indexed: 06/13/2023]
Abstract
We explore crystallisation and polymorph selection in active Brownian particles with numerical simulation. In agreement with previous work (Wysocki et al. in Europhys Lett 105:48004, 2014), we find that crystallisation is suppressed by activity and occurs at higher densities with increasing Péclet number ([Formula: see text]). While the nucleation rate decreases with increasing activity, the crystal growth rate increases due to the accelerated dynamics in the melt. As a result of this competition, we observe the transition from a nucleation and growth regime at high [Formula: see text] to "spinodal nucleation" at low [Formula: see text]. Unlike the case of passive hard spheres, where preference for FCC over HCP polymorphs is weak, activity causes the annealing of HCP stacking faults, thus strongly favouring the FCC symmetry at high [Formula: see text]. When freezing occurs more slowly, in the nucleation and growth regime, this tendency is much reduced and we see a trend towards the passive case of little preference for either polymorph.
Collapse
Affiliation(s)
- Fergus J. Moore
- Bristol Centre for Functional Nanomaterials, University of Bristol, Bristol, BS8 1FD UK
- H.H. Wills Physics Laboratory, Tyndall Ave., Bristol, BS8 1TL UK
| | - C. Patrick Royall
- H.H. Wills Physics Laboratory, Tyndall Ave., Bristol, BS8 1TL UK
- Gulliver UMR CNRS 7083, ESPCI Paris, Université PSL, 75005 Paris, France
- School of Chemistry, Cantock’s Close, University of Bristol, Bristol, BS8 1TS UK
| | | | - John Russo
- School of Mathematics, University of Bristol, Bristol, BS8 1UG UK
- Department of Physics, Sapienza University of Rome, P.le Aldo Moro 5, 00185 Rome, Italy
| |
Collapse
|
18
|
Moritz C, Geissler PL, Dellago C. The microscopic mechanism of bulk melting of ice. J Chem Phys 2021; 155:124501. [PMID: 34598556 DOI: 10.1063/5.0064380] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
We study the initial stages of homogeneous melting of a hexagonal ice crystal at coexistence and at moderate superheating. Our trajectory-based computer simulation approach provides a comprehensive picture of the events that lead to melting, from the initial accumulation of 5+7 defects, via the formation of L-D and interstitial-vacancy pairs, to the formation of a liquid nucleus. Of the different types of defects that we observe to be involved in melting, a particular kind of 5+7 type defect (type 5) plays a prominent role as it often forms prior to the formation of the initial liquid nucleus and close to the site where the nucleus forms. Hence, like other solids, ice homogeneously melts via the prior accumulation of defects.
Collapse
Affiliation(s)
- Clemens Moritz
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Phillip L Geissler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| |
Collapse
|
19
|
Smith RS, Tylinski M, Kimmel GA, Kay BD. Crystallization kinetics of amorphous acetonitrile nanoscale films. J Chem Phys 2021; 154:144703. [PMID: 33858151 DOI: 10.1063/5.0045461] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We measure the isothermal crystallization kinetics of amorphous acetonitrile films using molecular beam dosing and reflection adsorption infrared spectroscopy techniques. Experiments on a graphene covered Pt(111) substrate revealed that the crystallization rate slows dramatically during long time periods and that the overall kinetics cannot be described by a simple application of the Avrami equation. The crystallization kinetics also have a thickness dependence with the thinner films crystallizing much slower than the thicker ones. Additional experiments showed that decane layers at both the substrate and vacuum interfaces can also affect the crystallization rates. A comparison of the crystallization rates for CH3CN and CD3CN films showed only an isotope effect of ∼1.09. When amorphous films were deposited on a crystalline film, the crystalline layer did not act as a template for the formation of a crystalline growth front. These overall results suggest that the crystallization kinetics are complicated, indicating the possibility of multiple nucleation and growth mechanisms.
Collapse
Affiliation(s)
- R Scott Smith
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M Tylinski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Greg A Kimmel
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Bruce D Kay
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
20
|
Structural relaxation and crystallization in supercooled water from 170 to 260 K. Proc Natl Acad Sci U S A 2021; 118:2022884118. [PMID: 33790015 DOI: 10.1073/pnas.2022884118] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The origin of water's anomalous properties has been debated for decades. Resolution of the problem is hindered by a lack of experimental data in a crucial region of temperatures, T, and pressures where supercooled water rapidly crystallizes-a region often referred to as "no man's land." A recently developed technique where water is heated and cooled at rates greater than 109 K/s now enables experiments in this region. Here, it is used to investigate the structural relaxation and crystallization of deeply supercooled water for 170 K < T < 260 K. Water's relaxation toward a new equilibrium structure depends on its initial structure with hyperquenched glassy water (HQW) typically relaxing more quickly than low-density amorphous solid water (LDA). For HQW and T > 230 K, simple exponential relaxation kinetics is observed. For HQW at lower temperatures, increasingly nonexponential relaxation is observed, which is consistent with the dynamics expected on a rough potential energy landscape. For LDA, approximately exponential relaxation is observed for T > 230 K and T < 200 K, with nonexponential relaxation only at intermediate temperatures. At all temperatures, water's structure can be reproduced by a linear combination of two, local structural motifs, and we show that a simple model accounts for the complex kinetics within this context. The relaxation time, τ rel , is always shorter than the crystallization time, τ xtal For HQW, the ratio, τ xtal /τ rel , goes through a minimum at ∼198 K where the ratio is about 60.
Collapse
|
21
|
Tarn MD, Sikora SNF, Porter GCE, Shim JU, Murray BJ. Homogeneous Freezing of Water Using Microfluidics. MICROMACHINES 2021; 12:223. [PMID: 33672200 PMCID: PMC7926757 DOI: 10.3390/mi12020223] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Revised: 02/16/2021] [Accepted: 02/17/2021] [Indexed: 01/17/2023]
Abstract
The homogeneous freezing of water is important in the formation of ice in clouds, but there remains a great deal of variability in the representation of the homogeneous freezing of water in the literature. The development of new instrumentation, such as droplet microfluidic platforms, may help to constrain our understanding of the kinetics of homogeneous freezing via the analysis of monodisperse, size-selected water droplets in temporally and spatially controlled environments. Here, we evaluate droplet freezing data obtained using the Lab-on-a-Chip Nucleation by Immersed Particle Instrument (LOC-NIPI), in which droplets are generated and frozen in continuous flow. This high-throughput method was used to analyse over 16,000 water droplets (86 μm diameter) across three experimental runs, generating data with high precision and reproducibility that has largely been unrepresented in the microfluidic literature. Using this data, a new LOC-NIPI parameterisation of the volume nucleation rate coefficient (JV(T)) was determined in the temperature region of -35.1 to -36.9 °C, covering a greater JV(T) compared to most other microfluidic techniques thanks to the number of droplets analysed. Comparison to recent theory suggests inconsistencies in the theoretical representation, further implying that microfluidics could be used to inform on changes to parameterisations. By applying classical nucleation theory (CNT) to our JV(T) data, we have gone a step further than other microfluidic homogeneous freezing examples by calculating the stacking-disordered ice-supercooled water interfacial energy, estimated to be 22.5 ± 0.7 mJ m-2, again finding inconsistencies when compared to theoretical predictions. Further, we briefly review and compile all available microfluidic homogeneous freezing data in the literature, finding that the LOC-NIPI and other microfluidically generated data compare well with commonly used non-microfluidic datasets, but have generally been obtained with greater ease and with higher numbers of monodisperse droplets.
Collapse
Affiliation(s)
- Mark D. Tarn
- School of Earth and Environment, University of Leeds, Leeds LS2 9JT, UK; (S.N.F.S.); (G.C.E.P.)
- School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK;
| | - Sebastien N. F. Sikora
- School of Earth and Environment, University of Leeds, Leeds LS2 9JT, UK; (S.N.F.S.); (G.C.E.P.)
| | - Grace C. E. Porter
- School of Earth and Environment, University of Leeds, Leeds LS2 9JT, UK; (S.N.F.S.); (G.C.E.P.)
- School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK;
| | - Jung-uk Shim
- School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK;
| | - Benjamin J. Murray
- School of Earth and Environment, University of Leeds, Leeds LS2 9JT, UK; (S.N.F.S.); (G.C.E.P.)
| |
Collapse
|
22
|
Bianco V, de Hijes PM, Lamas CP, Sanz E, Vega C. Anomalous Behavior in the Nucleation of Ice at Negative Pressures. PHYSICAL REVIEW LETTERS 2021; 126:015704. [PMID: 33480790 DOI: 10.1103/physrevlett.126.015704] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 10/14/2020] [Accepted: 12/08/2020] [Indexed: 06/12/2023]
Abstract
Ice nucleation is a phenomenon that, despite the relevant implications for life, atmospheric sciences, and technological applications, is far from being completely understood, especially under extreme thermodynamic conditions. In this work we present a computational investigation of the homogeneous ice nucleation at negative pressures. By means of the seeding technique we estimate the size of the ice critical nucleus N_{c} for the TIP4P/Ice water model. This is done along the isotherms 230, 240, and 250 K, from positive to negative pressures until reaching the liquid-gas kinetic stability limit (where cavitation cannot be avoided). We find that N_{c} is nonmonotonic upon depressurization, reaching a minimum at negative pressures in the doubly metastable region of water. According to classical nucleation theory we establish the nucleation rate J and the surface tension γ, revealing a retracing behavior of both when the liquid-gas kinetic stability limit is approached. We also predict a reentrant behavior of the homogeneous nucleation line. The reentrance of these properties is related to the reentrance of the coexistence line at negative pressure, revealing new anomalies of water. The results of this work suggest the possibility of having metastable samples of liquid water for long times at negative pressure provided that heterogeneous nucleation is suppressed.
Collapse
Affiliation(s)
- Valentino Bianco
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - P Montero de Hijes
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - Cintia P Lamas
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| |
Collapse
|
23
|
Bianco V, Espinosa JR, Vega C. Antifreeze proteins and homogeneous nucleation: On the physical determinants impeding ice crystal growth. J Chem Phys 2020; 153:091102. [PMID: 32891082 DOI: 10.1063/5.0023211] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Antifreeze proteins (AFPs) are biopolymers capable of interfering with ice growth. Their antifreeze action is commonly understood considering that the AFPs, by pinning the ice surface, force the crystal-liquid interface to bend forming an ice meniscus, causing an increase in the surface free energy and resulting in a decrease in the freezing point ΔTmax. Here, we present an extensive computational study for a model protein adsorbed on a TIP4P/Ice crystal, computing ΔTmax as a function of the average distance d between AFPs, with simulations spanning over 1 µs. First, we show that the lower the d, the larger the ΔTmax. Then, we find that the water-ice-protein contact angle along the line ΔTmax(d) is always larger than 0°, and we provide a theoretical interpretation. We compute the curvature radius of the stable solid-liquid interface at a given supercooling ΔT ≤ ΔTmax, connecting it with the critical ice nucleus at ΔT. Finally, we discuss the antifreeze capability of AFPs in terms of the protein-water and protein-ice interactions. Our findings establish a unified description of the AFPs in the contest of homogeneous ice nucleation, elucidating key aspects of the antifreeze mechanisms and paving the way for the design of novel ice-controlling materials.
Collapse
Affiliation(s)
- Valentino Bianco
- Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid 28040, Spain
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge CB3 0H3, United Kingdom
| | - Carlos Vega
- Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid 28040, Spain
| |
Collapse
|
24
|
Rosales-Pelaez P, Sanchez-Burgos I, Valeriani C, Vega C, Sanz E. Seeding approach to nucleation in the NVT ensemble: The case of bubble cavitation in overstretched Lennard Jones fluids. Phys Rev E 2020; 101:022611. [PMID: 32168559 DOI: 10.1103/physreve.101.022611] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Accepted: 01/24/2020] [Indexed: 11/07/2022]
Abstract
Simulations are widely used to study nucleation in first order phase transitions due to the fact that they have access to the relevant length and time scales. However, simulations face the problem that nucleation is an activated process. Therefore, rare event simulation techniques are needed to promote the formation of the critical nucleus. The Seeding method, where the simulations are started with the nucleus already formed, has proven quite useful in efficiently providing estimates of the nucleation rate for a wide range of orders of magnitude. So far, Seeding has been employed in the NPT ensemble, where the nucleus either grows or redissolves. Thus, several trajectories have to be run in order to find the thermodynamic conditions that make the seeded nucleus critical. Moreover, the nucleus lifetime is short and the statistics for obtaining its properties is consequently poor. To deal with these shortcomings we extend the Seeding method to the NVT ensemble. We focus on the problem of bubble nucleation in a metastable Lennard Jones fluid. We show that, in the NVT ensemble, it is possible to equilibrate and stabilise critical bubbles for a long time. The nucleation rate inferred from NVT-Seeding is fully consistent with that coming from NPT-Seeding. The former is quite suitable to obtain the nucleation rate along isotherms, whereas the latter is preferable if the dependence of the rate with temperature at constant pressure is required. Care should be taken with finite size effects when using NVT-Seeding. Further work is required to extend NVT seeding to other sorts of phase transitions.
Collapse
Affiliation(s)
- P Rosales-Pelaez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - I Sanchez-Burgos
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Valeriani
- Departamento de Estructura de la Materia, Física Termica y Electronica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
25
|
Karmakar T, Piaggi PM, Parrinello M. Molecular Dynamics Simulations of Crystal Nucleation from Solution at Constant Chemical Potential. J Chem Theory Comput 2019; 15:6923-6930. [DOI: 10.1021/acs.jctc.9b00795] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Tarak Karmakar
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computationali, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - Pablo M. Piaggi
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computationali, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computationali, Università della Svizzera Italiana, Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| |
Collapse
|
26
|
Montero de Hijes P, Espinosa JR, Sanz E, Vega C. Interfacial free energy of a liquid-solid interface: Its change with curvature. J Chem Phys 2019; 151:144501. [PMID: 31615240 DOI: 10.1063/1.5121026] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
We analyze the changes in the interfacial free energy between a spherical solid cluster and a fluid due to the change of the radius of the solid. Interfacial free energies from nucleation studies using the seeding technique for four different systems, being hard spheres, Lennard-Jones, and two models of water (mW and TIP4P/ICE), were plotted as a function of the inverse of the radius of the solid cluster. In all cases, the interfacial free energy was a linear function of the inverse of the radius of the solid cluster and this is consistent with Tolman's equation. This linear behavior is shown not only in isotherms but also along isobars. The effect of curvature on the interfacial free energy is more pronounced in water, followed by hard spheres, and smaller for Lennard-Jones particles. We show that it is possible to estimate nucleation rates of Lennard-Jones particles at different pressures by using information from simple NpT simulations and taking into account the variation of the interfacial free energy with the radius of the solid cluster. Neglecting the effects of the radius on the interfacial free energy (capillarity approximation) leads to incorrect values of the nucleation rate. For the Lennard-Jones system, the homogeneous nucleation curve is not parallel to the melting curve as was found for water in previous work. This is due to the increase in the interfacial free energy along the coexistence curve as the pressure increases. This work presents a simple and relatively straightforward way to approximately estimate nucleation rates.
Collapse
Affiliation(s)
- P Montero de Hijes
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge CB3 0H3, United Kingdom
| | - Eduardo Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
27
|
Leoni F, Shi R, Tanaka H, Russo J. Crystalline clusters in mW water: Stability, growth, and grain boundaries. J Chem Phys 2019; 151:044505. [DOI: 10.1063/1.5100812] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Fabio Leoni
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | - Rui Shi
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - John Russo
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
28
|
Montero de Hijes P, Espinosa JR, Vega C, Sanz E. Ice growth rate: Temperature dependence and effect of heat dissipation. J Chem Phys 2019; 151:044509. [DOI: 10.1063/1.5103273] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- P. Montero de Hijes
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. R. Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, United Kingdom
| | - C. Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E. Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
29
|
Niu H, Yang YI, Parrinello M. Temperature Dependence of Homogeneous Nucleation in Ice. PHYSICAL REVIEW LETTERS 2019; 122:245501. [PMID: 31322390 DOI: 10.1103/physrevlett.122.245501] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2019] [Revised: 05/06/2019] [Indexed: 06/10/2023]
Abstract
Ice nucleation is a process of great relevance in physics, chemistry, technology, and environmental sciences; much theoretical effort has been devoted to its understanding, but it still remains a topic of intense research. We shed light on this phenomenon by performing atomistic based simulations. Using metadynamics and a carefully designed set of collective variables, reversible transitions between water and ice are able to be simulated. We find that water freezes into a stacking disordered structure with the all-atom transferable intermolecular potential with 4 points/ice (TIP4P/ice) model, and the features of the critical nucleus of nucleation at the microscopic level are revealed. We have also estimated the ice nucleation rates along with other nucleation parameters at different undercoolings. Our results are in agreement with recent experimental and other theoretical works, and they confirm that nucleation is preceded by a large increase in tetrahedrally coordinated water molecules.
Collapse
Affiliation(s)
- Haiyang Niu
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, and National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
| | - Yi Isaac Yang
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, and National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, and National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera Italiana, Via Giuseppe Buffi 13, 6900 Lugano, Switzerland
- Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
30
|
Liang S, Hall KW, Laaksonen A, Zhang Z, Kusalik PG. Characterizing key features in the formation of ice and gas hydrate systems. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2019; 377:20180167. [PMID: 30982452 PMCID: PMC6501917 DOI: 10.1098/rsta.2018.0167] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 02/26/2019] [Indexed: 05/16/2023]
Abstract
Crystallization in liquids is critical to a range of important processes occurring in physics, chemistry and life sciences. In this article, we review our efforts towards understanding the crystallization mechanisms, where we focus on theoretical modelling and molecular simulations applied to ice and gas hydrate systems. We discuss the order parameters used to characterize molecular ordering processes and how different order parameters offer different perspectives of the underlying mechanisms of crystallization. With extensive simulations of water and gas hydrate systems, we have revealed unexpected defective structures and demonstrated their important roles in crystallization processes. Nucleation of gas hydrates can in most cases be characterized to take place in a two-step mechanism where the nucleation occurs via intermediate metastable precursors, which gradually reorganizes to a stable crystalline phase. We have examined the potential energy landscapes explored by systems during nucleation, and have shown that these landscapes are rugged and funnel-shaped. These insights provide a new framework for understanding nucleation phenomena that has not been addressed in classical nucleation theory. This article is part of the theme issue 'The physics and chemistry of ice: scaffolding across scales, from the viability of life to the formation of planets'.
Collapse
Affiliation(s)
- Shuai Liang
- Key Laboratory of Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangdong Key Laboratory of New and Renewable Energy Research and Development, Guangzhou, People's Republic of China
| | - Kyle Wm. Hall
- Department of Chemistry, Temple University, Philadelphia, PA, USA
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
- Department of Chemistry-Ångström Laboratory, Uppsala University, 75121 Uppsala, Sweden
- Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Aleea Grigore Ghica-Voda, 41A, 700487 Iasi, Romania
| | - Zhengcai Zhang
- Department of Chemistry, University of Calgary, Calgary, Canada
| | | |
Collapse
|
31
|
Moreau DW, Atakisi H, Thorne RE. Ice formation and solvent nanoconfinement in protein crystals. IUCRJ 2019; 6:346-356. [PMID: 31098016 PMCID: PMC6503922 DOI: 10.1107/s2052252519001878] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Accepted: 01/31/2019] [Indexed: 05/06/2023]
Abstract
Ice formation within protein crystals is a major obstacle to the cryocrystallographic study of protein structure, and has limited studies of how the structural ensemble of a protein evolves with temperature in the biophysically interesting range from ∼260 K to the protein-solvent glass transition near 200 K. Using protein crystals with solvent cavities as large as ∼70 Å, time-resolved X-ray diffraction was used to study the response of protein and internal solvent during rapid cooling. Solvent nanoconfinement suppresses freezing temperatures and ice-nucleation rates so that ice-free, low-mosaicity diffraction data can be reliably collected down to 200 K without the use of cryoprotectants. Hexagonal ice (Ih) forms in external solvent, but internal crystal solvent forms stacking-disordered ice (Isd) with a near-random stacking of cubic and hexagonal planes. Analysis of powder diffraction from internal ice and single-crystal diffraction from the host protein structure shows that the maximum crystallizable solvent fraction decreases with decreasing crystal solvent-cavity size, and that an ∼6 Å thick layer of solvent adjacent to the protein surface cannot crystallize. These results establish protein crystals as excellent model systems for the study of nanoconfined solvent. By combining fast cooling, intense X-ray beams and fast X-ray detectors, complete structural data sets for high-value targets, including membrane proteins and large complexes, may be collected at ∼220-240 K that have much lower mosaicities and comparable B factors, and that may allow more confident identification of ligand binding than in current cryocrystallographic practice.
Collapse
Affiliation(s)
- David W. Moreau
- Physics Department, Cornell University, Ithaca, NY 14853, USA
| | - Hakan Atakisi
- Physics Department, Cornell University, Ithaca, NY 14853, USA
| | | |
Collapse
|
32
|
Franco Pinheiro Moreira PA, Gomes de Aguiar Veiga R, de Koning M. Elastic constants of ice Ihas described by semi-empirical water models. J Chem Phys 2019; 150:044503. [DOI: 10.1063/1.5082743] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
| | - Roberto Gomes de Aguiar Veiga
- Centro de Engenharia, Modelagem e Ciências Sociais Aplicadas, Universidade Federal do ABC, Santo André, 09210-580 São Paulo, Brazil
| | - Maurice de Koning
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, UNICAMP, Campinas, 13083-859 São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, Campinas, 13083-861 São Paulo, Brazil
| |
Collapse
|
33
|
Espinosa JR, Diez AL, Vega C, Valeriani C, Ramirez J, Sanz E. Ice Ih vs. ice III along the homogeneous nucleation line. Phys Chem Chem Phys 2019; 21:5655-5660. [DOI: 10.1039/c8cp07432a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The melting and the homogeneous nucleation lines of the TIP4P/Ice water model closely resemble the experimental ones. Computer simulations show that a change in the nucleating ice polymorph from ice Ih to ice III justifies the slope sign inversion of the homogeneous nucleation line.
Collapse
Affiliation(s)
- Jorge R. Espinosa
- Departamento de Quimica Fisica
- Facultad de Ciencias Quimicas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| | - Angel L. Diez
- Departamento de Quimica Fisica
- Facultad de Ciencias Quimicas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| | - Carlos Vega
- Departamento de Quimica Fisica
- Facultad de Ciencias Quimicas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| | - Chantal Valeriani
- Departamento de Estructura de la Materia
- Fisica Termica y Electronica
- Facultad de Ciencias Fisicas
- Universidad Complutense de Madrid
- 28040 Madrid
| | - Jorge Ramirez
- Departamento de Ingenieria Quimica Industrial y Medio Ambiente
- Escuela Tecnica Superior de Ingenieros Industriales
- Universidad Politecnica de Madrid
- 28006 Madrid
- Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica
- Facultad de Ciencias Quimicas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| |
Collapse
|
34
|
Zubieta I, Vázquez del Saz M, Llombart P, Vega C, Noya EG. Nucleation of pseudo hard-spheres and dumbbells at moderate metastability: appearance of A15 Frank–Kasper phase at intermediate elongations. Phys Chem Chem Phys 2019; 21:1656-1670. [PMID: 30383878 DOI: 10.1039/c8cp04964e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Crystal nucleation of repulsive hard-dumbbells from the sphere to the two tangent spheres limit is investigated at moderately high metastability by brute-force molecular dynamics simulations.
Collapse
Affiliation(s)
- Itziar Zubieta
- Instituto de Química Física Rocasolano
- Consejo Superior de Investigaciones Científicas
- CSIC
- 28006 Madrid
- Spain
| | - Miguel Vázquez del Saz
- Instituto de Química Física Rocasolano
- Consejo Superior de Investigaciones Científicas
- CSIC
- 28006 Madrid
- Spain
| | - Pablo Llombart
- Instituto de Química Física Rocasolano
- Consejo Superior de Investigaciones Científicas
- CSIC
- 28006 Madrid
- Spain
| | - Carlos Vega
- Departamento de Química Física (Unidad Asociada de I+D+i al CSIC)
- Facultad de Ciencias Químicas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| | - Eva G. Noya
- Instituto de Química Física Rocasolano
- Consejo Superior de Investigaciones Científicas
- CSIC
- 28006 Madrid
- Spain
| |
Collapse
|
35
|
Santos-Flórez PA, Ruestes CJ, de Koning M. Uniaxial-deformation behavior of ice Ih as described by the TIP4P/Ice and mW water models. J Chem Phys 2018; 149:164711. [DOI: 10.1063/1.5048517] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Pedro Antonio Santos-Flórez
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
| | - Carlos J. Ruestes
- CONICET and Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Mendoza, Argentina
| | - Maurice de Koning
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Center for Computational Engineering and Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861 Campinas, São Paulo, Brazil
| |
Collapse
|
36
|
Zhang Z, Kusalik PG, Guo GJ. Bridging solution properties to gas hydrate nucleation through guest dynamics. Phys Chem Chem Phys 2018; 20:24535-24538. [PMID: 30247495 DOI: 10.1039/c8cp04466j] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
By investigating the aqueous solution properties of several hydrate guests with molecular simulations, we find that with increasing guest concentration, the guest's hydration shell becomes more ordered and the system entropy decreases. A common critical value of the self-diffusion coefficient of different guest molecules is identified, below which hydrates will nucleate very readily.
Collapse
Affiliation(s)
- Zhengcai Zhang
- Key Laboratory of Petroleum Resource Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
| | | | | |
Collapse
|
37
|
Saito S, Bagchi B, Ohmine I. Crucial role of fragmented and isolated defects in persistent relaxation of deeply supercooled water. J Chem Phys 2018; 149:124504. [DOI: 10.1063/1.5044458] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shinji Saito
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan and The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Biman Bagchi
- Indian Institute of Science, Bangalore 560012, India
| | - Iwao Ohmine
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan and The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
38
|
Svoboda M, Lísal M. Concentrated aqueous sodium chloride solution in clays at thermodynamic conditions of hydraulic fracturing: Insight from molecular dynamics simulations. J Chem Phys 2018; 148:222806. [DOI: 10.1063/1.5017166] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Martin Svoboda
- Department of Molecular and Mesoscopic Modelling, Institute of Chemical Process Fundamentals of the CAS, v. v. i., Prague, Czech Republic
- Department of Physics, Faculty of Science, J. E. Purkinje University, Ústí nad Labem, Czech Republic
| | - Martin Lísal
- Department of Molecular and Mesoscopic Modelling, Institute of Chemical Process Fundamentals of the CAS, v. v. i., Prague, Czech Republic
- Department of Physics, Faculty of Science, J. E. Purkinje University, Ústí nad Labem, Czech Republic
| |
Collapse
|
39
|
Zaragoza A, Espinosa JR, Ramos R, Antonio Cobos J, Luis Aragones J, Vega C, Sanz E, Ramírez J, Valeriani C. Phase boundaries, nucleation rates and speed of crystal growth of the water-to-ice transition under an electric field: a simulation study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:174002. [PMID: 29508769 DOI: 10.1088/1361-648x/aab464] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We investigate with computer simulations the effect of applying an electric field on the water-to-ice transition. We use a combination of state-of-the-art simulation techniques to obtain phase boundaries and crystal growth rates (direct coexistence), nucleation rates (seeding) and interfacial free energies (seeding and mold integration). First, we consider ice Ih, the most stable polymorph in the absence of a field. Its normal melting temperature, speed of crystal growth and nucleation rate (for a given supercooling) diminish as the intensity of the field goes up. Then, we study polarised cubic ice, or ice Icf, the most stable solid phase under a strong electric field. Its normal melting point goes up with the field and, for a given supercooling, under the studied field (0.3 V nm-1) ice Icf nucleates and grows at a similar rate as Ih with no field. The net effect of the field would then be that ice nucleates at warmer temperatures, but in the form of ice Icf. The main conclusion of this work is that reasonable electric fields (not strong enough to break water molecules apart) are not relevant in the context of homogeneous ice nucleation at 1 bar.
Collapse
Affiliation(s)
- Alberto Zaragoza
- Departamento de Estructura de la Materia, Fisica Termica y Electronica, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain. Departamento de Ingenieria Fisica, Division de Ciencias e Ingenierias, Universidad de Guanajuato, Loma del Bosque 103, Col. Lomas del Campestre, CP 37150 Leon, Mexico
| | | | | | | | | | | | | | | | | |
Collapse
|
40
|
Affiliation(s)
- Andrew J. Amaya
- William G. Lowrie Department of Chemical and Biomolecular Engineering, The Ohio State University, Columbus, Ohio 43210, USA
| | - Barbara E. Wyslouzil
- William G. Lowrie Department of Chemical and Biomolecular Engineering, The Ohio State University, Columbus, Ohio 43210, USA
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
41
|
Cheng B, Dellago C, Ceriotti M. Theoretical prediction of the homogeneous ice nucleation rate: disentangling thermodynamics and kinetics. Phys Chem Chem Phys 2018; 20:28732-28740. [DOI: 10.1039/c8cp04561e] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We computed the homogeneous ice nucleation rate, and disentangled and investigated all the relevant physical quantities.
Collapse
Affiliation(s)
- Bingqing Cheng
- Laboratory of Computational Science and Modeling
- Institute of Materials
- École Polytechnique Fédérale de Lausanne
- 1015 Lausanne
- Switzerland
| | | | - Michele Ceriotti
- Laboratory of Computational Science and Modeling
- Institute of Materials
- École Polytechnique Fédérale de Lausanne
- 1015 Lausanne
- Switzerland
| |
Collapse
|
42
|
Pipolo S, Salanne M, Ferlat G, Klotz S, Saitta AM, Pietrucci F. Navigating at Will on the Water Phase Diagram. PHYSICAL REVIEW LETTERS 2017; 119:245701. [PMID: 29286747 DOI: 10.1103/physrevlett.119.245701] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Indexed: 06/07/2023]
Abstract
Despite the simplicity of its molecular unit, water is a challenging system because of its uniquely rich polymorphism and predicted but yet unconfirmed features. Introducing a novel space of generalized coordinates that capture changes in the topology of the interatomic network, we are able to systematically track transitions among liquid, amorphous, and crystalline forms throughout the whole phase diagram of water, including the nucleation of crystals above and below the melting point. Our approach, based on molecular dynamics and enhanced sampling or free energy calculation techniques, is not specific to water and could be applied to very different structural phase transitions, paving the way towards the prediction of kinetic routes connecting polymorphic structures in a range of materials.
Collapse
Affiliation(s)
- S Pipolo
- Sorbonne Universités, UPMC Université Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, IMPMC, F-75005 Paris, France
| | - M Salanne
- Sorbonne Universités, UPMC Université Paris 06, CNRS, Laboratoire PHENIX, F-75005 Paris, France
| | - G Ferlat
- Sorbonne Universités, UPMC Université Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, IMPMC, F-75005 Paris, France
| | - S Klotz
- Sorbonne Universités, UPMC Université Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, IMPMC, F-75005 Paris, France
| | - A M Saitta
- Sorbonne Universités, UPMC Université Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, IMPMC, F-75005 Paris, France
| | - F Pietrucci
- Sorbonne Universités, UPMC Université Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, IMPMC, F-75005 Paris, France
| |
Collapse
|
43
|
Pugliese P, Conde MM, Rovere M, Gallo P. Freezing Temperatures, Ice Nanotubes Structures, and Proton Ordering of TIP4P/ICE Water inside Single Wall Carbon Nanotubes. J Phys Chem B 2017; 121:10371-10381. [DOI: 10.1021/acs.jpcb.7b06306] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- P. Pugliese
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. M. Conde
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. Rovere
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - P. Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| |
Collapse
|
44
|
Brini E, Fennell CJ, Fernandez-Serra M, Hribar-Lee B, Lukšič M, Dill KA. How Water's Properties Are Encoded in Its Molecular Structure and Energies. Chem Rev 2017; 117:12385-12414. [PMID: 28949513 PMCID: PMC5639468 DOI: 10.1021/acs.chemrev.7b00259] [Citation(s) in RCA: 199] [Impact Index Per Article: 28.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2017] [Indexed: 11/29/2022]
Abstract
How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties.
Collapse
Affiliation(s)
- Emiliano Brini
- Laufer
Center for Physical and Quantitative Biology, Department of Physics and Astronomy, and Department of
Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Christopher J. Fennell
- Department
of Chemistry, Oklahoma State University, Stillwater, Oklahoma 74078, United States
| | - Marivi Fernandez-Serra
- Laufer
Center for Physical and Quantitative Biology, Department of Physics and Astronomy, and Department of
Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Barbara Hribar-Lee
- Faculty
of Chemistry and Chemical Technology, University
of Ljubljana, Večna
pot 113, SI-1000 Ljubljana, Slovenia
| | - Miha Lukšič
- Faculty
of Chemistry and Chemical Technology, University
of Ljubljana, Večna
pot 113, SI-1000 Ljubljana, Slovenia
| | - Ken A. Dill
- Laufer
Center for Physical and Quantitative Biology, Department of Physics and Astronomy, and Department of
Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
45
|
Espinosa JR, Soria GD, Ramirez J, Valeriani C, Vega C, Sanz E. Role of Salt, Pressure, and Water Activity on Homogeneous Ice Nucleation. J Phys Chem Lett 2017; 8:4486-4491. [PMID: 28876070 DOI: 10.1021/acs.jpclett.7b01551] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Pure water can be substantially supercooled below the melting temperature without transforming into ice. The achievable supercooling can be enhanced by adding solutes or by applying hydrostatic pressure. Avoiding ice formation is of great importance in the cryopreservation of food or biological samples. In this Letter, we investigate the similarity between the effects of pressure and salt on ice formation using a combination of state-of-the-art simulation techniques. We find that both hinder ice formation by increasing the energetic cost of creating the ice-fluid interface. Moreover, we examine the widely accepted proposal that the ice nucleation rate for different pressures and solute concentrations can be mapped through the activity of water [ Koop , L. ; Tsias , P. Nature , 2000 , 406 , 611 ]. We show that such a proposal is not consistent with the nucleation rates predicted in our simulations because it does not include all parameters affecting ice nucleation. Therefore, even though salt and pressure have a qualitatively similar effect on ice formation, they cannot be quantitatively mapped onto one another.
Collapse
Affiliation(s)
- Jorge R Espinosa
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid , 28040 Madrid, Spain
| | - Guiomar D Soria
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid , 28040 Madrid, Spain
| | - Jorge Ramirez
- Departamento de Ingenieria Quimica Industrial y Medio Ambiente, Escuela Tecnica Superior de Ingenieros Industriales, Universidad Politecnica de Madrid , 28006 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Fisica Aplicada I, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid , 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid , 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid , 28040 Madrid, Spain
| |
Collapse
|