1
|
Wang Y, Huang B, Li Z. Electric Field-Enhanced Ion Rejection Rate in Freeze Desalination. Chemphyschem 2024:e202400397. [PMID: 38960874 DOI: 10.1002/cphc.202400397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 06/11/2024] [Accepted: 06/28/2024] [Indexed: 07/05/2024]
Abstract
Freeze desalination is an appealing method for seawater desalination through freezing seawater. The percentage of ions in the liquid phase, which is termed ion rejection rate, is a critical factor affecting the performance of freeze desalination. Improving the ion rejection rate is an important topic for freeze desalination. In this work, we investigate the effects of electric fields on the ion rejection rate during the freezing of seawater through molecular dynamics simulations. It is found that the ion rejection rate increases with increasing electric field strength. The enhanced ion rejection rate is due to the reduction of the energy barrier at the ice-water interface caused by the electric field, which affects the orientation of water molecules and ion-water interactions. However, the electric field hinders the ice growth rate, which affects the productivity of freeze desalination. Nevertheless, the finding in this work offers a new idea to improve the ion rejection rate. Practically, a trade-off needs to be found to optimize the overall performance of freeze desalination.
Collapse
Affiliation(s)
- Yixiang Wang
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Baoling Huang
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Zhigang Li
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| |
Collapse
|
2
|
Chu F, Li S, Zhao C, Feng Y, Lin Y, Wu X, Yan X, Miljkovic N. Interfacial ice sprouting during salty water droplet freezing. Nat Commun 2024; 15:2249. [PMID: 38480695 PMCID: PMC10937636 DOI: 10.1038/s41467-024-46518-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2023] [Accepted: 02/26/2024] [Indexed: 03/17/2024] Open
Abstract
Icing of seawater droplets is capable of causing catastrophic damage to vessels, buildings, and human life, yet it also holds great potential for enhancing applications such as droplet-based freeze desalination and anti-icing of sea sprays. While large-scale sea ice growth has been investigated for decades, the icing features of small salty droplets remain poorly understood. Here, we demonstrate that salty droplet icing is governed by salt rejection-accompanied ice crystal growth, resulting in freezing dynamics different from pure water. Aided by the observation of brine films emerging on top of frozen salty droplets, we propose a universal definition of freezing duration to quantify the icing rate of droplets having varying salt concentrations. Furthermore, we show that the morphology of frozen salty droplets is governed by ice crystals that sprout from the bottom of the brine film. These crystals grow until they pierce the free interface, which we term ice sprouting. We reveal that ice sprouting is controlled by condensation at the brine film free interface, a mechanism validated through molecular dynamics simulations. Our findings shed light on the distinct physics that govern salty droplet icing, knowledge that is essential for the development of related technologies.
Collapse
Affiliation(s)
- Fuqiang Chu
- School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing, 100083, China
- Beijing Key Laboratory of Energy Conservation and Emission Reduction for Metallurgical Industry, University of Science and Technology Beijing, Beijing, 100083, China
| | - Shuxin Li
- School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing, 100083, China
| | - Canjun Zhao
- Department of Energy and Power Engineering, Tsinghua University, Beijing, 100084, China
| | - Yanhui Feng
- School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing, 100083, China.
- Beijing Key Laboratory of Energy Conservation and Emission Reduction for Metallurgical Industry, University of Science and Technology Beijing, Beijing, 100083, China.
| | - Yukai Lin
- Department of Energy and Power Engineering, Tsinghua University, Beijing, 100084, China
| | - Xiaomin Wu
- Department of Energy and Power Engineering, Tsinghua University, Beijing, 100084, China.
| | - Xiao Yan
- Key Laboratory of Low-grade Energy Utilization Technologies and Systems, Chongqing University, Ministry of Education, Chongqing, 400030, China.
- Institute of Engineering Thermophysics, Chongqing University, Chongqing, 400030, China.
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
| | - Nenad Miljkovic
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
- Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
- Materials Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
- International Institute for Carbon Neutral Energy Research (WPI-I2CNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka, 819-0395, Japan.
| |
Collapse
|
3
|
Wang Y, Shi D, Liu G, Huang B, Li Z. Shear-Enhanced Ion Rejection during Seawater Freezing. J Phys Chem B 2023; 127:10404-10410. [PMID: 37997846 DOI: 10.1021/acs.jpcb.3c05432] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2023]
Abstract
Ion rejection during seawater freezing is the basis for freeze desalination. A high ion rejection rate is desired for improving the performance of freeze desalination. In this work, we propose a method to enhance the ion rejection rate through external shear, which is demonstrated through molecular dynamics (MD) simulations and experiments. MD simulations show that the ion rejection rate increases with an increasing shear rate. This is attributed to the disruption of the hydration bonds between ions and water molecules in the hydration shell caused by the shear. Consequently, the mobility of ions is increased, and the energy barrier is reduced at the ice-water interface such that ions have a greater chance of diffusing into the aqueous solution, leading to an enhanced ion rejection rate. The MD results in this work are qualitatively confirmed by experiments and provide insights into the enhancement of the ion rejection rate through external parameters.
Collapse
Affiliation(s)
- Yixiang Wang
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon 999077, Hong Kong
| | - Dachuang Shi
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon 999077, Hong Kong
| | - Gongze Liu
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon 999077, Hong Kong
| | - Baoling Huang
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon 999077, Hong Kong
| | - Zhigang Li
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon 999077, Hong Kong
| |
Collapse
|
4
|
Singh Y, Santra M, Singh RS. Anomalous Vapor and Ice Nucleation in Water at Negative Pressures: A Classical Density Functional Theory Study. J Phys Chem B 2023; 127:3312-3324. [PMID: 36989467 DOI: 10.1021/acs.jpcb.2c09136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/31/2023]
Abstract
In contrast to the abundance of work on the anomalous behavior of water, the relationship between the water's thermodynamic anomalies and kinetics of phase transition from metastable water is relatively unexplored. In this work, we have employed classical density functional theory to provide a unified and coherent picture of nucleation (both vapor and ice) from metastable water at negative pressure conditions. Our results suggest a peculiar nonmonotonic temperature dependence of vapor-liquid surface tension at temperatures where vapor-liquid coexistence is metastable with respect to the ice phase. The vapor nucleation barrier on isochoric cooling also shows a nonmonotonic temperature dependence. We further report that, for low density isochores, the temperature of the minimum vapor nucleation barrier (TΔΩv/min*) does not coincide with the temperature of maximum density (TMD) where metastability is maximum. The difference between the TΔΩv/min* and the TMD, however, decreases with increasing the density of the isochore. The vapor nucleation barrier along isobars shows an interesting crossover behavior in the vicinity of the Widom line on lowering the temperature. Our results on the ice nucleation suggest an anomalous retracing behavior of the nucleation barrier along isotherms at negative pressures and theoretically validate the recent findings that the reentrant ice(Ih)-liquid coexistence line can induce a drastic change in the kinetics of ice nucleation. Thus, this study establishes a direct connection between the metastable water's thermodynamic anomalies and the (vapor and ice) nucleation kinetics. In addition, this study provides deeper insights into the origin of the isothermal compressibility maximum on isochoric cooling.
Collapse
Affiliation(s)
- Yuvraj Singh
- Department of Physics, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| | - Mantu Santra
- School of Chemical and Materials Sciences, Indian Institute of Technology Goa, Ponda, Goa 403401, India
| | - Rakesh S Singh
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati, Andhra Pradesh 517507, India
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Lin M, Cao H, Li J. Control strategies of ice nucleation, growth, and recrystallization for cryopreservation. Acta Biomater 2023; 155:35-56. [PMID: 36323355 DOI: 10.1016/j.actbio.2022.10.056] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2022] [Revised: 10/20/2022] [Accepted: 10/26/2022] [Indexed: 02/02/2023]
Abstract
The cryopreservation of biomaterials is fundamental to modern biotechnology and biomedicine, but the biggest challenge is the formation of ice, resulting in fatal cryoinjury to biomaterials. To date, abundant ice control strategies have been utilized to inhibit ice formation and thus improve cryopreservation efficiency. This review focuses on the mechanisms of existing control strategies regulating ice formation and the corresponding applications to biomaterial cryopreservation, which are of guiding significance for the development of ice control strategies. Herein, basics related to biomaterial cryopreservation are introduced first. Then, the theoretical bases of ice nucleation, growth, and recrystallization are presented, from which the key factors affecting each process are analyzed, respectively. Ice nucleation is mainly affected by melting temperature, interfacial tension, shape factor, and kinetic prefactor, and ice growth is mainly affected by solution viscosity and cooling/warming rate, while ice recrystallization is inhibited by adsorption or diffusion mechanisms. Furthermore, the corresponding research methods and specific control strategies for each process are summarized. The review ends with an outlook of the current challenges and future perspectives in cryopreservation. STATEMENT OF SIGNIFICANCE: Ice formation is the major limitation of cryopreservation, which causes fatal cryoinjury to cryopreserved biomaterials. This review focuses on the three processes related to ice formation, called nucleation, growth, and recrystallization. The theoretical models, key influencing factors, research methods and corresponding ice control strategies of each process are summarized and discussed, respectively. The systematic introduction on mechanisms and control strategies of ice formation is instructive for the cryopreservation development.
Collapse
Affiliation(s)
- Min Lin
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China; Beijing Key Laboratory for CO(2) Utilization and Reduction Technology, Tsinghua University, Beijing 100084, China
| | - Haishan Cao
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China; Beijing Key Laboratory for CO(2) Utilization and Reduction Technology, Tsinghua University, Beijing 100084, China.
| | - Junming Li
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China; Beijing Key Laboratory for CO(2) Utilization and Reduction Technology, Tsinghua University, Beijing 100084, China
| |
Collapse
|
8
|
Roy S, Bocharova V, Stack AG, Bryantsev VS. Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2SO 4 Morphology. ACS APPLIED MATERIALS & INTERFACES 2022; 14:53213-53227. [PMID: 36395432 DOI: 10.1021/acsami.2c17475] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Herein, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.
Collapse
Affiliation(s)
- Santanu Roy
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| | - Vera Bocharova
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| | - Andrew G Stack
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| | - Vyacheslav S Bryantsev
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| |
Collapse
|
9
|
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
|
10
|
Whale TF. Disordering effect of the ammonium cation accounts for anomalous enhancement of heterogeneous ice nucleation. J Chem Phys 2022; 156:144503. [DOI: 10.1063/5.0084635] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Heterogeneous nucleation of ice from supercooled water is the process responsible for triggering nearly all ice formation in the natural environment. Understanding of heterogeneous ice nucleation is particularly key for understanding the formation of ice in clouds, which impacts weather and climate. While many effective ice nucleators are known the mechanisms of their actions remain poorly understood. Some inorganic nucleators have been found to nucleate ice at warmer temperatures in dilute ammonium solution than in pure water. This is surprising, analogous to salty water melting at a warmer temperature than pure water. Here, the magnitude of this effect is rationalized as being due to thermodynamically favorable ammonium-induced disordering of the hydrogen bond network of ice critical clusters formed on inorganic ice nucleators. Theoretical calculations are shown to be consistent with new experimental measurements aimed at finding the maximum magnitude of the effect. The implication of this study is that the ice-nucleating sites and surfaces of many inorganic ice nucleators are either polar or charged and therefore tend to induce formation of hydrogen ordered ice clusters. This work corroborates various literature reports indicating that some inorganic ice nucleators are most effective when nominally neutral and implies a commonality in mechanism between a wide range of inorganic ice nucleators.
Collapse
Affiliation(s)
- Thomas F Whale
- Department of Chemistry, University of Warwick, United Kingdom
| |
Collapse
|
11
|
Bianco V, Conde MM, Lamas CP, Noya EG, Sanz E. Phase diagram of the NaCl–water system from computer simulations. J Chem Phys 2022; 156:064505. [DOI: 10.1063/5.0083371] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- V. Bianco
- Departamento de Química Física (Unidad de I+D+i asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M. M. Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C. P. Lamas
- Departamento de Química Física (Unidad de I+D+i asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Científicas, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - E. G. Noya
- Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Científicas, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - E. Sanz
- Departamento de Química Física (Unidad de I+D+i asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
12
|
P Lamas C, R Espinosa J, M Conde M, Ramírez J, Montero de Hijes P, G Noya E, Vega C, Sanz E. Homogeneous nucleation of NaCl in supersaturated solutions. Phys Chem Chem Phys 2021; 23:26843-26852. [PMID: 34817484 DOI: 10.1039/d1cp02093e] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The seeding method is an approximate approach to investigate nucleation that combines molecular dynamics simulations with classical nucleation theory. Recently, this technique has been successfully implemented in a broad range of nucleation studies. However, its accuracy is subject to the arbitrary choice of the order parameter threshold used to distinguish liquid-like from solid-like molecules. We revisit here the crystallization of NaCl from a supersaturated brine solution and show that consistency between seeding and rigorous methods, like Forward Flux Sampling (from previous work) or spontaneous crystallization (from this work), is achieved by following a mislabelling criterion to select such threshold (i.e. equaling the fraction of the mislabelled particles in the bulk parent and nucleating phases). This work supports the use of seeding to obtain fast and reasonably accurate nucleation rate estimates and the mislabelling criterion as one giving the relevant cluster size for classical nucleation theory in crystallization studies.
Collapse
Affiliation(s)
- C P Lamas
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain. .,Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Científicas, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - J R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge CB3 0H3, UK
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006, Madrid, Spain
| | - J Ramírez
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006, Madrid, Spain
| | - P Montero de Hijes
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | - E G Noya
- Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Científicas, CSIC, Calle Serrano 119, 28006 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
|
13
|
Sanchez-Burgos I, Sanz E, Vega C, Espinosa JR. Fcc vs. hcp competition in colloidal hard-sphere nucleation: on their relative stability, interfacial free energy and nucleation rate. Phys Chem Chem Phys 2021; 23:19611-19626. [PMID: 34524277 DOI: 10.1039/d1cp01784e] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hard-sphere crystallization has been widely investigated over the last six decades by means of colloidal suspensions and numerical methods. However, some aspects of its nucleation behaviour are still under debate. Here, we provide a detailed computational characterisation of the polymorphic nucleation competition between the face-centered cubic (fcc) and the hexagonal-close packed (hcp) hard-sphere crystal phases. By means of several state-of-the-art simulation techniques, we evaluate the melting pressure, chemical potential difference, interfacial free energy and nucleation rate of these two polymorphs, as well as of a random stacking mixture of both crystals. Our results highlight that, despite the fact that both polymorphs have very similar stability, the interfacial free energy of the hcp phase could be marginally higher than that of the fcc solid, which in consequence, mildly decreases its propensity to nucleate from the liquid compared to the fcc phase. Moreover, we analyse the abundance of each polymorph in grown crystals from different types of inserted nuclei: fcc, hcp and stacking disordered fcc/hcp seeds, as well as from those spontaneously emerged from brute force simulations. We find that post-critical crystals fundamentally grow maintaining the polymorphic structure of the critical nucleus, at least until moderately large sizes, since the only crystallographic orientation that allows stacking close-packed disorder is the fcc (111) plane, or equivalently the hcp (0001) one. Taken together, our results contribute with one more piece to the intricate puzzle of colloidal hard-sphere crystallization.
Collapse
Affiliation(s)
- Ignacio Sanchez-Burgos
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK.
| | - Eduardo Sanz
- 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
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK.
| |
Collapse
|
14
|
Luo S, Jin Y, Tao R, Li H, Li C, Wang J, Li Z. Molecular understanding of ion rejection in the freezing of aqueous solutions. Phys Chem Chem Phys 2021; 23:13292-13299. [PMID: 34095926 DOI: 10.1039/d1cp01733k] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, we investigate the microscopic mechanism of ion rejection phenomena during the freezing of aqueous NaCl solutions through molecular dynamics simulations. It is found that the hydration energy for the ion-water interaction is stronger than that between ions and ice, which is the fundamental reason giving rise to the phenomenon of ion rejection. The probability of ions being rejected by ice is determined by the competition between the energy barrier at the ice-water interface and the thermal effect. The ion rejection rate increases with increasing temperature. Furthermore, it is found that the rejection rate of Na+ is higher than that of Cl- because of the relatively large hydration energy difference between Na+-water and Na+-ice interactions. The role of temperature in the applications of ion rejection in freeze desalination is also discussed.
Collapse
Affiliation(s)
- Shuang Luo
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | - Yakang Jin
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | - Ran Tao
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | - Haiyang Li
- Key Laboratory of Enhanced Heat Transfer and Energy Conservation, Ministry of Education, Beijing University of Technology, Beijing 100124, P. R. China.
| | - Chu Li
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Jun Wang
- Key Laboratory of Enhanced Heat Transfer and Energy Conservation, Ministry of Education, Beijing University of Technology, Beijing 100124, P. R. China.
| | - Zhigang Li
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| |
Collapse
|
15
|
Metya AK, Molinero V. Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation. J Am Chem Soc 2021; 143:4607-4624. [PMID: 33729789 DOI: 10.1021/jacs.0c12012] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.
Collapse
Affiliation(s)
- Atanu K Metya
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
16
|
Sanchez-Burgos I, Garaizar A, Vega C, Sanz E, Espinosa JR. Parasitic crystallization of colloidal electrolytes: growing a metastable crystal from the nucleus of a stable phase. SOFT MATTER 2021; 17:489-505. [PMID: 33346291 DOI: 10.1039/d0sm01680b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Colloidal particles have been extensively used to comprehend the main principles governing liquid-crystal nucleation. Multiple mechanisms and frameworks have been proposed, through either experiments or computational approaches, to rationalise the ubiquitous formation of colloidal crystals. In this work, we elucidate the nucleation scenario behind the crystallization of oppositely charged colloids. By performing molecular dynamics simulations of colloidal electrolytes in combination with the Seeding technique, we evaluate the fundamental factors, such as the nucleation rate, free energy barrier, surface tension and kinetic pre-factor, that determine the liquid-to-solid transition of several crystalline polymorphs. Our results show that at a high packing fraction, there is a cross-over between the nucleation of the CsCl structure and that of a substitutionally disordered fcc phase, despite the CuAu crystal being the most stable phase. We demonstrate that the crucial factor in determining which phase nucleates the fastest is the free energy cost of the cluster formation rather than their kinetic ability to grow from the liquid. While at a low packing fraction, the stable phase, CsCl, is the one that nucleates and subsequently grows, we show how at moderate and high packing fractions, a disordered fcc phase subsequently grows regardless of the nature of the nucleating phase, termed parasitic crystallization. Taken together, our results provide a panoramic perspective of the complex nucleation scenario of oppositely charged colloids at moderate temperature and rationalise the different thermodynamic and kinetic aspects behind it.
Collapse
Affiliation(s)
- Ignacio Sanchez-Burgos
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK.
| | - Adiran Garaizar
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK.
| | - Carlos Vega
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK.
| |
Collapse
|
17
|
Sanchez-Burgos I, de Hijes PM, Rosales-Pelaez P, Vega C, Sanz E. Equivalence between condensation and boiling in a Lennard-Jones fluid. Phys Rev E 2020; 102:062609. [PMID: 33466022 DOI: 10.1103/physreve.102.062609] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Accepted: 11/30/2020] [Indexed: 06/12/2023]
Abstract
Condensation and boiling are phase transitions highly relevant to industry, geology, and atmospheric science. These phase transitions are initiated by the nucleation of a drop in a supersaturated vapor and of a bubble in an overstretched liquid, respectively. The surface tension between both phases, liquid and vapor, is a key parameter in the development of such nucleation stage. Whereas the surface tension can be readily measured for a flat interface, there are technical and conceptual limitations to obtain it for the curved interface of the nucleus. On the technical side, it is quite difficult to observe a critical nucleus in experiments. From a conceptual point of view, the interfacial free energy depends on the choice of the dividing surface, being the surface of tension the one relevant for nucleation. We bypass the technical limitation by performing simulations of a Lennard-Jones fluid where we equilibrate critical nuclei (both drops and bubbles). Regarding the conceptual hurdle, we find the relevant cluster size by searching the radius that correctly predicts nucleation rates and nucleation free energy barriers when combined with Classical Nucleation Theory. With such definition of the cluster size we find the same value of the surface tension for drops and bubbles of a given radius. Thus, condensation and boiling can be viewed as two sides of the same coin. Finally, we combine the data coming from drops and bubbles to obtain, via two different routes, estimates of the Tolman length, a parameter that allows describing the curvature dependence of the surface tension in a theoretical framework.
Collapse
Affiliation(s)
- I Sanchez-Burgos
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - P Montero de Hijes
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - P Rosales-Pelaez
- 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
| | - E Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
18
|
Predicting heterogeneous ice nucleation with a data-driven approach. Nat Commun 2020; 11:4777. [PMID: 32963232 PMCID: PMC7509812 DOI: 10.1038/s41467-020-18605-3] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Accepted: 08/28/2020] [Indexed: 01/05/2023] Open
Abstract
Water in nature predominantly freezes with the help of foreign materials through a process known as heterogeneous ice nucleation. Although this effect was exploited more than seven decades ago in Vonnegut's pioneering cloud seeding experiments, it remains unclear what makes a material a good ice former. Here, we show through a machine learning analysis of nucleation simulations on a database of diverse model substrates that a set of physical descriptors for heterogeneous ice nucleation can be identified. Our results reveal that, beyond Vonnegut's connection with the lattice match to ice, three new microscopic factors help to predict the ice nucleating ability. These are: local ordering induced in liquid water, density reduction of liquid water near the surface and corrugation of the adsorption energy landscape felt by water. With this we take a step towards quantitative understanding of heterogeneous ice nucleation and the in silico design of materials to control ice formation.
Collapse
|
19
|
Zhu C, Gao Y, Zhu W, Liu Y, Francisco JS, Zeng XC. Computational Prediction of Novel Ice Phases: A Perspective. J Phys Chem Lett 2020; 11:7449-7461. [PMID: 32787287 DOI: 10.1021/acs.jpclett.0c01635] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Although computational prediction of new ice phases is a niche field in water science, the scientific subject itself is representative of two important areas in physical chemistry, namely, statistical thermodynamics and molecular simulations. The prediction of a variety of novel ice phases has also attracted general public interest since the 1980s. In particular, the prediction of low-dimensional ice phases has gained momentum since the confirmation of a number of low-dimensional "computer ice" phases in the laboratory over the past decade. In this Perspective, the research advancements in computational prediction of novel ice phases over the past few years are reviewed. Particular attention is placed on new ice phases whose physical properties or dimensional structures are distinctly different from conventional bulk ices. Specific topics include the (i) formation of superionic ices, (ii) electrofreezing of water under high pressure and in a high external electric field, (iii) prediction of low-density porous ice at strongly negative pressure, (iv) ab initio computational study of two-dimensional (2D) ice under nanoscale confinement, and (v) 2D ices formed on a solid surface near ambient temperature without nanoscale confinement. Clearly, the formation of most of these novel ice phases demands certain extreme conditions. Ongoing challenges and new opportunities for predicting new ice phases from either classical molecular dynamics simulation or high-level ab initio computation are discussed.
Collapse
Affiliation(s)
- Chongqin Zhu
- Department of Earth and Environmental Science, Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Yurui Gao
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| | - Weiduo Zhu
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yuan Liu
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| | - Joseph S Francisco
- Department of Earth and Environmental Science, Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
- Department of Chemical & Biomolecular Engineering, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| |
Collapse
|
20
|
Panagiotopoulos AZ. Simulations of activities, solubilities, transport properties, and nucleation rates for aqueous electrolyte solutions. J Chem Phys 2020; 153:010903. [DOI: 10.1063/5.0012102] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
|
21
|
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
|
22
|
Schmelzer JWP, Abyzov AS, Baidakov VG. Entropy and the Tolman Parameter in Nucleation Theory. ENTROPY 2019; 21:e21070670. [PMID: 33267384 PMCID: PMC7515167 DOI: 10.3390/e21070670] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/23/2019] [Revised: 06/26/2019] [Accepted: 07/05/2019] [Indexed: 01/16/2023]
Abstract
Thermodynamic aspects of the theory of nucleation are commonly considered employing Gibbs’ theory of interfacial phenomena and its generalizations. Utilizing Gibbs’ theory, the bulk parameters of the critical clusters governing nucleation can be uniquely determined for any metastable state of the ambient phase. As a rule, they turn out in such treatment to be widely similar to the properties of the newly-evolving macroscopic phases. Consequently, the major tool to resolve problems concerning the accuracy of theoretical predictions of nucleation rates and related characteristics of the nucleation process consists of an approach with the introduction of the size or curvature dependence of the surface tension. In the description of crystallization, this quantity has been expressed frequently via changes of entropy (or enthalpy) in crystallization, i.e., via the latent heat of melting or crystallization. Such a correlation between the capillarity phenomena and entropy changes was originally advanced by Stefan considering condensation and evaporation. It is known in the application to crystal nucleation as the Skapski–Turnbull relation. This relation, by mentioned reasons more correctly denoted as the Stefan–Skapski–Turnbull rule, was expanded by some of us quite recently to the description of the surface tension not only for phase equilibrium at planar interfaces, but to the description of the surface tension of critical clusters and its size or curvature dependence. This dependence is frequently expressed by a relation derived by Tolman. As shown by us, the Tolman equation can be employed for the description of the surface tension not only for condensation and boiling in one-component systems caused by variations of pressure (analyzed by Gibbs and Tolman), but generally also for phase formation caused by variations of temperature. Beyond this particular application, it can be utilized for multi-component systems provided the composition of the ambient phase is kept constant and variations of either pressure or temperature do not result in variations of the composition of the critical clusters. The latter requirement is one of the basic assumptions of classical nucleation theory. For this reason, it is only natural to use it also for the specification of the size dependence of the surface tension. Our method, relying on the Stefan–Skapski–Turnbull rule, allows one to determine the dependence of the surface tension on pressure and temperature or, alternatively, the Tolman parameter in his equation. In the present paper, we expand this approach and compare it with alternative methods of the description of the size-dependence of the surface tension and, as far as it is possible to use the Tolman equation, of the specification of the Tolman parameter. Applying these ideas to condensation and boiling, we derive a relation for the curvature dependence of the surface tension covering the whole range of metastable initial states from the binodal curve to the spinodal curve.
Collapse
Affiliation(s)
- Jürn W. P. Schmelzer
- Institute of Physics, University of Rostock, Albert-Einstein-Strasse 23-25, 18059 Rostock, Germany
| | - Alexander S. Abyzov
- National Science Center Kharkov Institute of Physics and Technology, 61108 Kharkov, Ukraine
| | - Vladimir G. Baidakov
- Institute of Thermal Physics, Ural Branch of the Russian Academy of Sciences, Amundsen Street 107a, 620016 Yekaterinburg, Russia
- Correspondence: ; Tel.: +49-381-498-6889; Fax: +49-381-498-6882
| |
Collapse
|
23
|
Qiao Z, Zhao Y, Gao YQ. Ice Nucleation of Confined Monolayer Water Conforms to Classical Nucleation Theory. J Phys Chem Lett 2019; 10:3115-3121. [PMID: 31117689 DOI: 10.1021/acs.jpclett.9b01169] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We confirmed that monolayer water confined by parallel graphene sheets spontaneously crystallizes from a structurally and dynamically heterogeneous liquid phase under moderate supercooling via direct molecular dynamics simulation. Square-lattice-like geometric order is observed at the early stage of nucleation and is preserved during the entire nucleus growth process. The diffusion coefficient and free energy profile in the cluster space extracted from a Bayesian trajectory analysis agree well with the classical nucleation theory (CNT) prediction and yield thermodynamic quantities exhibiting linear temperature dependence. The effectiveness of maximum cluster size as the descriptor of ice nucleation dynamics in the CNT framework can be attributed to the dynamical time scale decoupling and strong structural pattern dependence of density fluctuation in the liquid phase.
Collapse
Affiliation(s)
- Zhuoran Qiao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking National Laboratory for Molecular Science , Peking University , Beijing 100871 , China
| | - Yuheng Zhao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking National Laboratory for Molecular Science , Peking University , Beijing 100871 , China
| | - Yi Qin Gao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking National Laboratory for Molecular Science , Peking University , Beijing 100871 , China
| |
Collapse
|
24
|
Supercooled Liquid Serum Physiologic Solution Instantly Crystallized on the Nurse Table Used for Cooling of Periorbital Region During Rhinoplasty. Aesthetic Plast Surg 2019; 43:453-456. [PMID: 30535556 DOI: 10.1007/s00266-018-1287-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2018] [Accepted: 11/24/2018] [Indexed: 10/27/2022]
Abstract
INTRODUCTION Formation of less periorbital ecchymosis in post-operative period of rhinoplasty is a popular trend. We present the use of instantly crystallizing supercooled serum physiologic solution for periorbital cooling. PHYSICS OF SUPERCOOLING There are circumstances in which water temperature drops below its freezing point, but no phase transition happens while water remains in the liquid phase. This is called supercooling. Pure water can be supercooled below the freezing temperature without transforming into ice. Tap water will not supercool because it contains impurities that serve as nucleation sites for crystallization. For freezer temperatures in the range of - 4 °C, - 6 °C, and - 8 °C, nucleation was not observed and pure water remained in the supercooled condition for a long time. DESCRIPTION OF THE TECHNIQUE Sterile serum physiologic solution at + 5 °C can be supercooled in the freezer at - 14 °C only between the 257 and 277 min time interval. But when it is supercooled in the freezer at - 8 °C it is possible to save it in liquid form for at least 7 days as we have observed in our trials. CLINICAL USE AND DISCUSSION It is easily possible to transform this supercooled liquid sterile serum physiologic within a few seconds into moldable snow-like ice that can be used safely and more nicely rather than solid ice for periorbital cooling in rhinoplasty operations. Its sterile inner bag is held tight and struck over the sterile nurse table and it crystallizes within a few seconds. For frozen solutions, tearing of the inner plastic bag and extracting the ice and then crushing of big masses of ice to small pieces is exhaustive and a time-consuming process. The temperature of the supercooled fluid will be zero at the moment of nucleation with no risk of frostbite. The crystallized serum physiologic solution preserves its ice-gel form for nearly 25 min. CONCLUSION The instant crystallization of supercooled liquid serum physiologic solution can be applied as a tissue cooling method in rhinoplasty and in several other surgical procedures. LEVEL OF EVIDENCE IV This journal requires that authors assign a level of evidence to each article. For a full description of these Evidence-Based Medicine ratings, please refer to the Table of Contents or the online Instructions to Authors www.springer.com/00266 .
Collapse
|
25
|
Soria GD, Espinosa JR, Ramirez J, Valeriani C, Vega C, Sanz E. A simulation study of homogeneous ice nucleation in supercooled salty water. J Chem Phys 2018; 148:222811. [DOI: 10.1063/1.5008889] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Guiomar D. Soria
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge R. Espinosa
- 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 Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- 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
|
26
|
Whale TF, Holden MA, Wilson TW, O'Sullivan D, Murray BJ. The enhancement and suppression of immersion mode heterogeneous ice-nucleation by solutes. Chem Sci 2018; 9:4142-4151. [PMID: 29780544 PMCID: PMC5941198 DOI: 10.1039/c7sc05421a] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Accepted: 03/26/2018] [Indexed: 01/05/2023] Open
Abstract
Heterogeneous nucleation of ice from aqueous solutions is an important yet poorly understood process in multiple fields, not least the atmospheric sciences where it impacts the formation and properties of clouds. In the atmosphere ice-nucleating particles are usually, if not always, mixed with soluble material. However, the impact of this soluble material on ice nucleation is poorly understood. In the atmospheric community the current paradigm for freezing under mixed phase cloud conditions is that dilute solutions will not influence heterogeneous freezing. By testing combinations of nucleators and solute molecules we have demonstrated that 0.015 M solutions (predicted melting point depression <0.1 °C) of several ammonium salts can cause suspended particles of feldspars and quartz to nucleate ice up to around 3 °C warmer than they do in pure water. In contrast, dilute solutions of certain alkali metal halides can dramatically depress freezing points for the same nucleators. At 0.015 M, solutes can enhance or deactivate the ice-nucleating ability of a microcline feldspar across a range of more than 10 °C, which corresponds to a change in active site density of more than a factor of 105. This concentration was chosen for a survey across multiple solutes-nucleant combinations since it had a minimal colligative impact on freezing and is relevant for activating cloud droplets. Other nucleators, for instance a silica gel, are unaffected by these 'solute effects', to within experimental uncertainty. This split in response to the presence of solutes indicates that different mechanisms of ice nucleation occur on the different nucleators or that surface modification of relevance to ice nucleation proceeds in different ways for different nucleators. These solute effects on immersion mode ice nucleation may be of importance in the atmosphere as sea salt and ammonium sulphate are common cloud condensation nuclei (CCN) for cloud droplets and are internally mixed with ice-nucleating particles in mixed-phase clouds. In addition, we propose a pathway dependence where activation of CCN at low temperatures might lead to enhanced ice formation relative to pathways where CCN activation occurs at higher temperatures prior to cooling to nucleation temperature.
Collapse
Affiliation(s)
- Thomas F Whale
- School of Earth and Environment , University of Leeds , Leeds , LS2 9JT , UK .
| | - Mark A Holden
- School of Earth and Environment , University of Leeds , Leeds , LS2 9JT , UK .
- School of Chemistry , University of Leeds , Leeds , LS2 9JT , UK
- School of Physics and Astronomy , University of Leeds , Leeds , LS29JT , UK
| | - Theodore W Wilson
- School of Earth and Environment , University of Leeds , Leeds , LS2 9JT , UK .
| | - Daniel O'Sullivan
- School of Earth and Environment , University of Leeds , Leeds , LS2 9JT , UK .
| | - Benjamin J Murray
- School of Earth and Environment , University of Leeds , Leeds , LS2 9JT , UK .
| |
Collapse
|
27
|
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
|
28
|
Xu Y, Petrik NG, Smith RS, Kay BD, Kimmel GA. Homogeneous Nucleation of Ice in Transiently-Heated, Supercooled Liquid Water Films. J Phys Chem Lett 2017; 8:5736-5743. [PMID: 29125304 DOI: 10.1021/acs.jpclett.7b02685] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We have investigated the nucleation and growth of crystalline ice in 0.24 μm thick, supercooled water films adsorbed on Pt(111). The films were transiently heated with ∼10 ns infrared laser pulses, which produced typical heating and cooling rates of ∼109-1010 K/s. The crystallization of these water films was monitored with infrared spectroscopy. The experimental conditions were chosen to suppress ice nucleation at both the water/metal and water/vacuum interfaces. Furthermore, internal pressure increases due to curvature effects are precluded in these flat films. Therefore, the experiments were sensitive to the homogeneous ice nucleation rate from ∼210 to 225 K. The experiments show that Jmax, the maximum for the homogeneous ice nucleation rate, J(T), needs to be ≥1026 m-3 s-1 and is likely to be ∼1029±2 m-3 s-1. We argue that such large nucleation rates are consistent with experiments on hyperquenched glassy water, which typically have crystalline fractions of ∼1% or more.
Collapse
Affiliation(s)
- Yuntao Xu
- Physical Sciences Division , Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Nikolay G Petrik
- Physical Sciences Division , Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - R Scott Smith
- Physical Sciences Division , Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Bruce D Kay
- Physical Sciences Division , Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Greg A Kimmel
- Physical Sciences Division , Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|