1
|
Tipeev AO, Gurashkin AL, Zanotto ED. Exploring surface properties and premelting in crystals. J Chem Phys 2024; 160:224705. [PMID: 38864371 DOI: 10.1063/5.0210127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2024] [Accepted: 05/27/2024] [Indexed: 06/13/2024] Open
Abstract
Crystal surfaces play a pivotal role in governing various significant processes, such as adsorption, nucleation, wetting, friction, and wear. A fundamental property that influences these processes is the surface free energy, γ. We have directly calculated γ(T) for low-index faces of Lennard-Jones (LJ), germanium, and silicon crystals along their sublimation lines using the computational cleavage technique. Our calculations agree well with experimental values for Si(111) and Ge(111), highlighting the accuracy of the method and models used. For LJ crystals, we identified a premelting onset at Tpm = 0.75Tm, marked by a sharp increase in atom mobility within the second outermost surface layer. Notably, Tpm closely aligned with the endpoint of the LJ melting line at negative pressures, Tend = 0.76Tm. We hypothesize that the emergence and coexistence of a liquid film atop the LJ crystal at Tpm < T < Tm correspond to the metastable melting line under negative pressures experienced by stretched crystal surfaces. Furthermore, our study of thin LJ crystal slabs reveals that premelting-induced failure leads to recrystallization below the homogeneous freezing limit, offering a promising avenue to explore crystal nucleation and growth at extremely deep supercoolings. Finally, no evidence of premelting was detected in the model crystals of Ge and Si, which is consistent with the experimental observations. Overall, our findings offer valuable insights into crystal surface phenomena at the atomic scale.
Collapse
Affiliation(s)
- Azat O Tipeev
- Department of Materials Engineering, Federal University of São Carlos, 13.565-905 São Carlos, SP, Brazil
| | - Alexander L Gurashkin
- Institute of Thermal Physics, Ural Branch of the Russian Academy of Sciences, 620016 Ekaterinburg, Russia
| | - Edgar D Zanotto
- Department of Materials Engineering, Federal University of São Carlos, 13.565-905 São Carlos, SP, Brazil
| |
Collapse
|
2
|
Boström M, Li Y, Brevik I, Persson C, Carretero-Palacios S, Malyi OI. van der Waals induced ice growth on partially melted ice nuclei in mist and fog. Phys Chem Chem Phys 2023; 25:32709-32714. [PMID: 38014720 DOI: 10.1039/d3cp04157c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Ice nucleation and formation play pivotal roles across various domains, from environmental science to food engineering. However, the exact ice formation mechanisms remain incompletely understood. This study introduces a novel ice formation process, which can be either heterogeneous or homogeneous, depending on the initial conditions. The process initiates ice crystal growth from a nucleus composed of a micron-sized partially melted ice particle. We explore the role of van der Waals (Lifshitz)-free energy and its resulting stress in the accumulation of ice at the interface with water vapor. Our analysis suggests that this process could lead to thicknesses ranging from nanometers to micrometers, depending on the size and degree of initial melting of the ice nucleus. We provide evidence for the growth of thin ice layers instead of liquid water films on a partially melted ice-vapor interface, offering some insights into mist and fog formation. We also link it to potential atmospheric and astrogeophysical applications.
Collapse
Affiliation(s)
- M Boström
- Centre of Excellence ENSEMBLE3 Sp. z o. o., Wolczynska Str. 133, 01-919, Warsaw, Poland.
| | - Y Li
- School of Physics and Materials Science, Nanchang University, Nanchang 330031, China.
- Institute of Space Science and Technology, Nanchang University, Nanchang 330031, China
| | - I Brevik
- Department of Energy and Process Engineering, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| | - C Persson
- Department of Materials Science and Engineering, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
- Centre for Materials Science and Nanotechnology, Department of Physics, University of Oslo, P. O. Box 1048 Blindern, NO-0316 Oslo, Norway
| | - S Carretero-Palacios
- Departamento de Física de Materiales and Instituto de Ciencia de Materiales Nicolás Cabrera, Universidad Autónoma de Madrid, 28049 Madrid, Spain
- Instituto de Ciencia de Materiales de Madrid, ICMM-CSIC,C/Sor Juana Inés de la Cruz, 3, Madrid 28049, Spain
| | - O I Malyi
- Centre of Excellence ENSEMBLE3 Sp. z o. o., Wolczynska Str. 133, 01-919, Warsaw, Poland.
| |
Collapse
|
3
|
Li Y, Brevik I, Malyi OI, Boström M. Different pathways to anomalous stabilization of ice layers on methane hydrates. Phys Rev E 2023; 108:034801. [PMID: 37849091 DOI: 10.1103/physreve.108.034801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Accepted: 08/03/2023] [Indexed: 10/19/2023]
Abstract
We explore the Casimir-Lifshitz free-energy theory for surface freezing of methane gas hydrates near the freezing point of water. The theory enables us to explore different pathways, resulting in anomalous (stabilizing) ice layers on methane hydrate surfaces via energy minimization. Notably, we will contrast the gas hydrate material properties, under which thin ice films can form in water vapor, with those previously predicted to be required in the presence of liquid water. It is predicted that methane hydrates in water vapor near the freezing point of water nucleate ice films, instead of water films.
Collapse
Affiliation(s)
- Y Li
- School of Physics and Materials Science, Nanchang University, Nanchang 330031, China
- Institute of Space Science and Technology, Nanchang University, Nanchang 330031, China
| | - I Brevik
- Department of Energy and Process Engineering, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| | - O I Malyi
- Centre of Excellence ENSEMBLE3 Sp. z o. o., Wolczynska Strasse 133, 01-919, Warsaw, Poland
| | - M Boström
- Centre of Excellence ENSEMBLE3 Sp. z o. o., Wolczynska Strasse 133, 01-919, Warsaw, Poland
| |
Collapse
|
4
|
Kamat K, Naullage PM, Molinero V, Peters B. Oriented attachment kinetics for rod-like particles at a flat surface: Buffon's needle at the nanoscale. J Chem Phys 2022; 157:214113. [PMID: 36511557 DOI: 10.1063/5.0124531] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
The adsorption of large rod-like molecules or crystallites on a flat crystal face, similar to Buffon's needle, requires the rods to "land," with their binding sites in precise orientational alignment with matching sites on the surface. An example is provided by long, helical antifreeze proteins (AFPs), which bind at specific facets and orientations on the ice surface. The alignment constraint for adsorption, in combination with the loss in orientational freedom as the molecule diffuses toward the surface, results in an entropic barrier that hinders the adsorption. Prior kinetic models do not factor in the complete geometry of the molecule, nor explicitly enforce orientational constraints for adsorption. Here, we develop a diffusion-controlled adsorption theory for AFP molecules binding at specific orientations to flat ice surfaces. We formulate the diffusion equation with relevant boundary conditions and present analytical solutions to the attachment rate constant. The resulting rate constant is a function of the length and aspect ratio of the AFP, the distance threshold associated with binding, and solvent conditions such as temperature and viscosity. These results and methods of calculation may also be useful for predicting the kinetics of crystal growth through oriented attachment.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Pavithra M Naullage
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Baron Peters
- Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
5
|
Berrens ML, Bononi FC, Donadio D. Effect of sodium chloride adsorption on the surface premelting of ice. Phys Chem Chem Phys 2022; 24:20932-20940. [PMID: 36040383 DOI: 10.1039/d2cp02277j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We characterise the structural properties of the quasi-liquid layer (QLL) at two low-index ice surfaces in the presence of sodium chloride (Na+/Cl-) ions by molecular dynamics simulations. We find that the presence of a high surface density of Na+/Cl- pairs changes the surface melting behaviour from step-wise to gradual melting. The ions lead to an overall increase of the thickness and the disorder of the QLL, and to a low-temperature roughening transition of the air-ice interface. The local molecular structure of the QLL is similar to that of liquid water, and the differences between the basal and primary prismatic surface are attenuated by the presence of Na+/Cl- pairs. These changes modify the crystal growth rates of different facets and the solvation environment at the surface of sea-water ice with a potential impact on light scattering and environmental chemical reactions.
Collapse
Affiliation(s)
- Margaret L Berrens
- Department of Chemistry, University of California Davis, Davis, CA, 95616, USA.
| | - Fernanda C Bononi
- Department of Chemistry, University of California Davis, Davis, CA, 95616, USA.
| | - Davide Donadio
- Department of Chemistry, University of California Davis, Davis, CA, 95616, USA.
| |
Collapse
|
6
|
Luengo-Márquez J, Izquierdo-Ruiz F, MacDowell LG. Intermolecular forces at ice and water interfaces: premelting, surface freezing and regelation. J Chem Phys 2022; 157:044704. [DOI: 10.1063/5.0097378] [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
Using Lifshitz theory we assess the role of van der Waals forces at interfaces of ice and water. The results are combined with measured structural forces from computer simulations to develop a quantitative model of the surface free energy of premelting films. This input is employed within the framework of wetting theory and allows us to predict qualitatively the behavior of quasi-liquid layer thickness as a function of ambient conditions. Our results emphasize the significance of vapor pressure. The ice vapor interface is shown to exhibit only incomplete premelting, but the situation can shift to a state of complete surface melting above water saturation. The results obtained serve also to assess the role of subsurface freezing at the water-vapor interface, and we show that intermolecular forces favor subsurface ice nucleation only in conditions of water undersaturation. We show ice regelation at ambient pressure may be explained as a process of capillary freezing, without the need to invoke the action of bulk pressure melting. Our results for van der Waals forces are exploited in order to gauge dispersion interactions in empirical point charge models of water.
Collapse
Affiliation(s)
| | | | - Luis G. MacDowell
- Dpto. de Quimica Fisica, Universidad Complutense de Madrid Facultad de Ciencias Químicas, Spain
| |
Collapse
|
7
|
Cui S, Chen H, Zhao Z. Premelting layer during ice growth: role of clusters. Phys Chem Chem Phys 2022; 24:15330-15339. [PMID: 35703342 DOI: 10.1039/d2cp00412g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The premelting layer plays an important role in ice growth, but there is a significant gap in our knowledge between the atomistic premelting surface structure and the macroscopic growth mechanism. In this work, using large-scale molecular dynamics simulation, we reveal the existence of clusters on the premelting surface, as an intermediate feature bridging the gap. We show the spontaneous formation and evolution of clusters, and they form a stable distribution determined by the growth rate. We demonstrate how this stable distribution is related to the growth mode of ice, connected by the growth of clusters. We come to a bilayer-by-bilayer growth mode at simulation-reachable high growth rates, but another mechanism, namely "cluster stacking", is speculated to exist at lower growth rates. This work builds a connection between the microscopic structure of the premelting layer and the macroscopic growth of ice, making a step forward toward the full understanding of premelting and ice growth.
Collapse
Affiliation(s)
- Shifan Cui
- International Center for Quantum Materials, School of Physics, Peking University, 209 Chengfu Road, Haidian District, Beijing 100871, China.
| | - Haoxiang Chen
- School of Physics, Peking University, 209 Chengfu Road, Haidian District, Beijing 100871, China
| | - Zhengpu Zhao
- International Center for Quantum Materials, School of Physics, Peking University, 209 Chengfu Road, Haidian District, Beijing 100871, China.
| |
Collapse
|
8
|
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
|
9
|
Qiu H, Zhao W, Zhou W, Guo W. Edge premelting of two-dimensional ices. J Chem Phys 2021; 155:044706. [PMID: 34340399 DOI: 10.1063/5.0056732] [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
The surface of a three-dimensional ice crystal naturally has a quasi-liquid layer (QLL) at temperatures below its bulk melting point, due to a phenomenon called surface premelting. Here, we show that the edges of a two-dimensional (2D) bilayer hexagonal ice adsorbed on solid surfaces undergo premelting as well, resulting in the formation of quasi-liquid bands (QLBs) at the edges. Our extensive molecular dynamics simulations show that the QLB exhibits structure and dynamics indistinguishable from the bilayer liquid phase, acting as a lower-dimensional analog of the QLL on the bulk ice. We further find that at low temperatures, the width of the QLBs at armchair-type edges of the 2D ice is almost identical to that at zigzag-type edges but becomes far greater than the latter at temperatures near the melting point. The chirality-dependent edge premelting of 2D ices should add an important new ingredient to the heterogeneity of premelting.
Collapse
Affiliation(s)
- Hu Qiu
- State Key Laboratory of Mechanics and Control of Mechanical Structures and Key Laboratory for Intelligent Nano Materials and Devices of MOE, Institute of Nano Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
| | - Wen Zhao
- State Key Laboratory of Mechanics and Control of Mechanical Structures and Key Laboratory for Intelligent Nano Materials and Devices of MOE, Institute of Nano Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
| | - Wanqi Zhou
- State Key Laboratory of Mechanics and Control of Mechanical Structures and Key Laboratory for Intelligent Nano Materials and Devices of MOE, Institute of Nano Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
| | - Wanlin Guo
- State Key Laboratory of Mechanics and Control of Mechanical Structures and Key Laboratory for Intelligent Nano Materials and Devices of MOE, Institute of Nano Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
| |
Collapse
|
10
|
Niblett SP, Limmer DT. Ion Dissociation Dynamics in an Aqueous Premelting Layer. J Phys Chem B 2021; 125:2174-2181. [DOI: 10.1021/acs.jpcb.0c11286] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Samuel P. Niblett
- Materials Science Division, Lawrence Berkeley Laboratory, Berkeley, California 94720, United States
| | - David T. Limmer
- Chemistry Department, University of California, Berkeley, California 94720, United States
- Chemical Science Division, Lawrence Berkeley Laboratory, Berkeley, California 94720, United States
- Kavli Energy Nanosciences Institute, University of California, Berkeley, California 94720, United States
| |
Collapse
|
11
|
Sibley DN, Llombart P, Noya EG, Archer AJ, MacDowell LG. How ice grows from premelting films and water droplets. Nat Commun 2021; 12:239. [PMID: 33431836 PMCID: PMC7801427 DOI: 10.1038/s41467-020-20318-6] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2020] [Accepted: 11/14/2020] [Indexed: 11/08/2022] Open
Abstract
Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes. At low saturation, freezing proceeds by terrace spreading, but the motion of the buried solid is conveyed through the liquid to the outer liquid-vapor interface. At higher saturations water droplets condense, a large crater forms below, and freezing proceeds undetectably beneath the droplet. Our approach is a general framework that naturally models freezing close to three phase coexistence and provides a first principle theory of ice growth and melting which may prove useful in the geosciences.
Collapse
Affiliation(s)
- David N Sibley
- Department of Mathematical Sciences, Loughborough University, Loughborough, LE11 3TU, UK
| | - Pablo Llombart
- Instituto de Química Física Rocasolano, CSIC, Calle Serrano 119, Madrid, 28006, Spain
- Departamento de Química Física (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid, 28040, Spain
| | - Eva G Noya
- Instituto de Química Física Rocasolano, CSIC, Calle Serrano 119, Madrid, 28006, Spain
| | - Andrew J Archer
- Department of Mathematical Sciences, Loughborough University, Loughborough, LE11 3TU, UK
| | - Luis G MacDowell
- Departamento de Química Física (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid, 28040, Spain.
| |
Collapse
|
12
|
Bononi FC, Chen Z, Rocca D, Andreussi O, Hullar T, Anastasio C, Donadio D. Bathochromic Shift in the UV–Visible Absorption Spectra of Phenols at Ice Surfaces: Insights from First-Principles Calculations. J Phys Chem A 2020; 124:9288-9298. [DOI: 10.1021/acs.jpca.0c07038] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Fernanda C. Bononi
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| | - Zekun Chen
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| | - Dario Rocca
- Université de Lorraine, CNRS, LPTC, F-54000 Nancy, France
| | - Oliviero Andreussi
- Department of Physics, University of North Texas Denton, Texas 76203, United States
| | - Ted Hullar
- Department of Land, Air and Water Resources, University of California Davis Davis, California 95616-8627, United States
| | - Cort Anastasio
- Department of Land, Air and Water Resources, University of California Davis Davis, California 95616-8627, United States
| | - Davide Donadio
- Department of Chemistry, University of California Davis, Davis, California 95616-5270, United States
| |
Collapse
|
13
|
Llombart P, Noya EG, MacDowell LG. Surface phase transitions and crystal habits of ice in the atmosphere. SCIENCE ADVANCES 2020; 6:eaay9322. [PMID: 32671203 PMCID: PMC7314560 DOI: 10.1126/sciadv.aay9322] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Accepted: 03/06/2020] [Indexed: 05/28/2023]
Abstract
With climate modeling predicting a raise of at least 2°C by year 2100, the fate of ice has become a serious concern, but we still do not understand how ice grows (or melts). In the atmosphere, crystal growth rates of basal and prism facets exhibit an enigmatic temperature dependence and crossover up to three times in a range between 0° and -40°. Here, we use large-scale computer simulations to characterize the ice surface and identify a sequence of previously unidentified phase transitions on the main facets of ice crystallites. Unexpectedly, we find that as temperature is increased, the crystal surface transforms from a disordered phase with proliferation of steps to a smooth phase with small step density. This causes the anomalous increase of step free energies and provides the long sought explanation for the enigmatic crossover of snow crystal growth rates found in the atmosphere.
Collapse
Affiliation(s)
- Pablo Llombart
- Instituto de Química Física Rocasolano, Madrid, Spain
- Departamento de Química Física, Universidad Complutense de Madrid, Madrid, Spain
| | - Eva G. Noya
- Instituto de Química Física Rocasolano, Madrid, Spain
| | - Luis G. MacDowell
- Departamento de Química Física, Universidad Complutense de Madrid, Madrid, Spain
| |
Collapse
|
14
|
Llombart P, Noya EG, Sibley DN, Archer AJ, MacDowell LG. Rounded Layering Transitions on the Surface of Ice. PHYSICAL REVIEW LETTERS 2020; 124:065702. [PMID: 32109130 DOI: 10.1103/physrevlett.124.065702] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2019] [Revised: 11/29/2019] [Accepted: 01/14/2020] [Indexed: 06/10/2023]
Abstract
Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers. However, disjoining pressure curves, which closely follow expectations from a renormalized mean field liquid state theory, show that there are no layering phase transitions in the thermodynamic sense along the sublimation line. Instead, we find that transitions at mean field level are rounded due to capillary wave fluctuations. We see signatures that true first order layering transitions could arise at low temperatures, for pressures between the metastable line of water-vapor coexistence and the sublimation line. The extrapolation of the disjoining pressure curve above water-vapor saturation displays a true first order phase transition from a thin to a thick film consistent with experimental observations.
Collapse
Affiliation(s)
- Pablo Llombart
- 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, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - Eva G Noya
- Instituto de Química Física Rocasolano, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - David N Sibley
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Andrew J Archer
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Luis G MacDowell
- 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
|
15
|
Esteso V, Carretero-Palacios S, MacDowell LG, Fiedler J, Parsons DF, Spallek F, Míguez H, Persson C, Buhmann SY, Brevik I, Boström M. Premelting of ice adsorbed on a rock surface. Phys Chem Chem Phys 2020; 22:11362-11373. [DOI: 10.1039/c9cp06836h] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Considering ice-premelting on a quartz rock surface (i.e. silica) we calculate the Lifshitz excess pressures in a four layer system with rock–ice–water–air.
Collapse
|
16
|
Llombart P, Palafox MA, MacDowell LG, Noya EG. Structural transitions and bilayer formation of CTAB aggregates. Colloids Surf A Physicochem Eng Asp 2019. [DOI: 10.1016/j.colsurfa.2019.123730] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
17
|
Liang Z, Du H, Liang H, Yang Y. Grain boundary premelting of monolayer ices in 2D nano-channels. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1593532] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Zun Liang
- Physics Department, School of Physics and Material Science, East China Normal University, Shanghai, China
| | - Han Du
- Physics Department, School of Physics and Material Science, East China Normal University, Shanghai, China
| | - Hongtao Liang
- Physics Department, School of Physics and Material Science, East China Normal University, Shanghai, China
| | - Yang Yang
- Physics Department, School of Physics and Material Science, East China Normal University, Shanghai, China
| |
Collapse
|
18
|
Llombart P, Bergua RM, Noya EG, MacDowell LG. Structure and water attachment rates of ice in the atmosphere: role of nitrogen. Phys Chem Chem Phys 2019; 21:19594-19611. [PMID: 31464318 DOI: 10.1039/c9cp03728d] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work we perform computer simulations of the ice surface in order to elucidate the role of nitrogen in the crystal growth rates and crystal habits of snow in the atmosphere. In pure water vapor at temperatures typical of ice crystal formation in cirrus clouds, we find that basal and primary prismatic facets exhibit a layer of premelted ice, with thickness in the subnanometer range. For partial pressures of 1 bar, well above the expected values in the troposphere, we find that only small amounts of nitrogen are adsorbed. The adsorption takes place onto the premelted surface, and hardly any nitrogen dissolves within the premelting film. The premelting film thickness does not change either. We quantify the resulting change of the ice/vapor surface tension to be in the hundredth of mN m-1 and find that the structure of the pristine ice surface is not changed in a significant manner. We perform a trajectory analysis of colliding water molecules, and find that the attachment rates from direct ballistic collision are very close to unity irrespective of the nitrogen pressure. Nitrogen is however at sufficient density to deflect a fraction of trajectories with smaller distance than the mean free path. Our results show explicitly that the reported differences in growth rates measured in pure water vapor and a controlled nitrogen atmosphere are not related to a significant disruption of the ice surface due to nitrogen adsorption. On the contrary, we show clearly from our trajectory analysis that nitrogen slows down the crystal growth rates due to collisions between water molecules with bulk nitrogen gas. This clarifies the long standing controversy of the role of inert gases on crystal growth rates and demonstrates their influence is solely related to the diffusion limited flow of water vapor across the gas phase.
Collapse
Affiliation(s)
- Pablo Llombart
- Instituto de Química Física Rocasolano, CSIC, Calle Serrano 119, 28006 Madrid, Spain and 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.
| | - Ramon M Bergua
- 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.
| | - Eva G Noya
- Instituto de Química Física Rocasolano, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - Luis G MacDowell
- 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
|
19
|
Yin H, Sibley DN, Archer AJ. Binding potentials for vapour nanobubbles on surfaces using density functional theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:315102. [PMID: 30978706 DOI: 10.1088/1361-648x/ab18e8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We calculate density profiles of a simple model fluid in contact with a planar surface using density functional theory (DFT), in particular for the case where there is a vapour layer intruding between the wall and the bulk liquid. We apply the method of Hughes et al (2015 J. Chem. Phys. 142 074702) to calculate the density profiles for varying (specified) amounts of the vapour adsorbed at the wall. This is equivalent to varying the thickness h of the vapour at the surface. From the resulting sequence of density profiles we calculate the thermodynamic grand potential as h is varied and thereby determine the binding potential as a function of h. The binding potential obtained via this coarse-graining approach allows us to determine the disjoining pressure in the film and also to predict the shape of vapour nano-bubbles on the surface. Our microscopic DFT based approach captures information from length scales much smaller than some commonly used models in continuum mechanics.
Collapse
Affiliation(s)
- Hanyu Yin
- Department of Mathematical Sciences, Loughborough University, Loughborough, LE11 3TU, United Kingdom
| | | | | |
Collapse
|
20
|
Benet J, Llombart P, Sanz E, MacDowell LG. Structure and fluctuations of the premelted liquid film of ice at the triple point. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1583388] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Jorge Benet
- Departamento de Química-Física (Unidad Asociada de I+D+i al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid, Spain
| | - Pablo Llombart
- Departamento de Química-Física (Unidad Asociada de I+D+i al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid, Spain
- Instituto de Química Física Rocasolano, CSIC, Madrid, Spain
| | - Eduardo Sanz
- Departamento de Química-Física (Unidad Asociada de I+D+i al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid, Spain
| | - Luis G. MacDowell
- Departamento de Química-Física (Unidad Asociada de I+D+i al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid, Spain
| |
Collapse
|
21
|
|
22
|
Mitsui T, Aoki K. Fluctuation spectroscopy of surface melting of ice with and without impurities. Phys Rev E 2019; 99:010801. [PMID: 30780264 DOI: 10.1103/physreve.99.010801] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Indexed: 11/07/2022]
Abstract
Water is ubiquitous, and the surface properties of ice have been studied for some time, due to their importance. A liquidlike layer (LLL) is known to exist on ice, below the melting point. We use surface thermal fluctuation spectroscopy to study the LLL, including its thickness, for pure ice, and for ice with impurities. We find that the properties of the LLL are experimentally those of liquid water, with thickness much smaller than previous results. We also find that impurities cause the LLL to be thicker, and quite inhomogeneous, with properties depending on the dopant.
Collapse
Affiliation(s)
- Takahisa Mitsui
- Research and Education Center for Natural Sciences and Department of Physics, Hiyoshi, Keio University, Yokohama 223-8521, Japan
| | - Kenichiro Aoki
- Research and Education Center for Natural Sciences and Department of Physics, Hiyoshi, Keio University, Yokohama 223-8521, Japan
| |
Collapse
|
23
|
Murata KI, Nagashima K, Sazaki G. How Do Ice Crystals Grow inside Quasiliquid Layers? PHYSICAL REVIEW LETTERS 2019; 122:026102. [PMID: 30720327 DOI: 10.1103/physrevlett.122.026102] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2018] [Indexed: 06/09/2023]
Abstract
A microscopic understanding of crystal-melt interfaces, inseparably involved in the dynamics of crystallization, is a long-standing challenge in condensed matter physics. Here, using an advanced optical microscopy, we directly visualize growing interfaces between ice basal faces and quasiliquid layers (QLLs) during ice crystal growth. This system serves as a model for studying the molecular incorporation process of the crystal growth from a supercooled melt (the so-called melt growth), often hidden by inevitable latent heat diffusion and/or the extremely high crystal growth rate. We reveal that the growth of basal faces inside QLLs proceeds layer by layer via two-dimensional nucleation of monomolecular islands. We also find that the lateral growth rate of the islands is well described by the Wilson-Frenkel law, taking into account the slowing down of the dynamics of water molecules interfaced with ice. These results clearly indicate that, after averaging surface molecular fluctuations, the layer by layer stacking still survives even at the topmost layer on basal faces, which supports the kink-step-terrace picture even for the melt growth.
Collapse
Affiliation(s)
- Ken-Ichiro Murata
- Institute of Low Temperature Science, Hokkaido University, N19-W8, Kita-ku, Sapporo 060-0819, Japan
| | - Ken Nagashima
- Institute of Low Temperature Science, Hokkaido University, N19-W8, Kita-ku, Sapporo 060-0819, Japan
| | - Gen Sazaki
- Institute of Low Temperature Science, Hokkaido University, N19-W8, Kita-ku, Sapporo 060-0819, Japan
| |
Collapse
|
24
|
Abstract
Premelting of ice at temperatures below 0 °C is of fundamental importance for environmental processes. Various experimental techniques have been used to investigate the temperature at which liquid-like water first appears at the ice-vapor interface, reporting onset temperatures from -160 to -2 °C. The signals that identify liquid-like order at the ice-vapor interface in these studies, however, do not show a sharp initiation with temperature. That is at odds with the expected first-order nature of surface phase transitions, and consistent with recent large-scale molecular simulations that show the first premelted layer to be sparse and to develop continuously over a wide range of temperatures. Here we perform a thermodynamic analysis to elucidate the origin of the continuous formation of the first layer of liquid at the ice-vapor interface. We conclude that a negative value of the line tension of the ice-liquid-vapor three-phase contact line is responsible for the continuous character of the transition and the sparse nature of the liquid-like domains in the incomplete first layer.
Collapse
Affiliation(s)
- Yuqing Qiu
- Department of Chemistry , The University of Utah , Salt Lake City , Utah 84112-0580 , United States
| | - Valeria Molinero
- Department of Chemistry , The University of Utah , Salt Lake City , Utah 84112-0580 , United States
| |
Collapse
|
25
|
Moreira PAFP, Veiga RGDA, Ribeiro IDA, Freitas R, Helfferich J, de Koning M. Anomalous diffusion of water molecules at grain boundaries in ice I h. Phys Chem Chem Phys 2018; 20:13944-13951. [PMID: 29744498 DOI: 10.1039/c8cp00933c] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Using ab initio and classical molecular dynamics simulations, we study pre-melting phenomena in pristine coincident-site-lattice grain boundaries (GBs) in proton-disordered hexagonal ice Ih at temperatures just below the melting point Tm. Concerning pre-melt-layer thicknesses, the results are consistent with the available experimental estimates for low-disorder impurity-free GBs. With regard to molecular mobility, the simulations provide a key new insight: the translational motion of the water molecules is found to be subdiffusive for time scales from ∼10 ns up to at least 0.1 μs. Moreover, the fact that the anomalous diffusion occurs even at temperatures just below Tm where the bulk supercooled liquid still diffuses normally suggests that it is related to the confinement of the GB pre-melt layers by the surrounding crystalline environment. Furthermore, we show that this behavior can be characterized by continuous-time random walk models in which the waiting-time distributions decay according to power-laws that are very similar to those describing dynamics in glass-forming systems.
Collapse
|
26
|
Conde MM, Rovere M, Gallo P. High precision determination of the melting points of water TIP4P/2005 and water TIP4P/Ice models by the direct coexistence technique. J Chem Phys 2017; 147:244506. [DOI: 10.1063/1.5008478] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- 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
|
27
|
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
|
28
|
Demange G, Zapolsky H, Patte R, Brunel M. Growth kinetics and morphology of snowflakes in supersaturated atmosphere using a three-dimensional phase-field model. Phys Rev E 2017; 96:022803. [PMID: 28950463 DOI: 10.1103/physreve.96.022803] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2017] [Indexed: 04/26/2023]
Abstract
Simulating ice crystal growth is a major issue for meteorology and aircraft safety. Yet, very few models currently succeed in reproducing correctly the diversity of snow crystal forms, and link the model parameters to thermodynamic quantities. Here, we demonstrate that the new three-dimensional phase-field model developed in Demange et al. [npj Comput. Mater. 3, 1 (2017)2057-396010.1038/s41524-017-0015-1] is capable of reproducing properly the morphology and growth kinetics of snowflakes in supersaturated atmosphere. Aside from that, we show that the growth dynamics of snow crystals satisfies the selection theory, consistently with previous experimental observations. Finally, we link the parameters of the phase-field model to atmospheric parameters.
Collapse
Affiliation(s)
- G Demange
- GPM, UMR CNRS 6643, University of Rouen, 76575 Saint Étienne du Rouvray, France
| | - H Zapolsky
- GPM, UMR CNRS 6643, University of Rouen, 76575 Saint Étienne du Rouvray, France
| | - R Patte
- GPM, UMR CNRS 6643, University of Rouen, 76575 Saint Étienne du Rouvray, France
| | - M Brunel
- CORIA UMR 6614, University of Rouen, 76575 Saint Étienne du Rouvray, France
| |
Collapse
|
29
|
|