1
|
Sinha VK, Das CK. Effect of confinement on water properties in super-hydrophilic pores using MD simulations with the mW model. J Mol Model 2024; 30:345. [PMID: 39316190 DOI: 10.1007/s00894-024-06145-2] [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: 06/24/2024] [Accepted: 09/12/2024] [Indexed: 09/25/2024]
Abstract
CONTEXT We explore the influence of strongly hydrophilic confinement on various properties of water, such as density, enthalpy, potential energy, radial distribution function, entropy, specific heat capacity, structural dynamics, and transition temperatures (freezing and melting temperatures), using monatomic water (mW) model. The properties of water are found to be dependent on confinement and the wall-fluid surface interaction. Hysteresis loops are observed for density, enthalpy, potential energy, and entropy around the transition temperatures, while the size of hysteresis loops varies with confinement and surface interaction. In smaller pore sizes (H ≤ 20), the solid phase displays a higher density compared to the liquid phase, which is unconventional behavior compared to bulk water systems due to the pronounced hydrophilic properties of the confinement surface. Specific heat capacity exhibits more oscillations in the confined system compared to bulk water, stemming from uneven enthalpy differences across equal temperature intervals. During phase transformation in both heating and quenching processes, there is an abrupt change observed in specific heat capacity. Confinement exerts a notable impact on entropy in the solid phase, but its influence is negligible in the liquid phase. At lower pore sizes (H < 25 Å), there is more fluctuation in freezing temperature for all wall-fluid interactions, which diminishes beyond pore sizes of H > 25 Å. Similarly, more oscillatory behavior is observed in melting temperatures at lower pore sizes (H < 40 Å), which diminishes at higher pore sizes (H > 40 Å). During the quenching process, a sudden jump in the in-plane orientational and tetrahedral order parameters indicates the formation of an ordered phase, specifically a diamond crystalline structure. The percentages of different crystalline structures (cubic diamond, hexagonal diamond, and 2D hexagonal) vary with both the confinement size and the wall-fluid interaction strength. METHODS Cooling and heating simulations are conducted with the mW water model using LAMMPS for different nanoscale confinement separation sizes ranging from 8.5 to 70 Å within the temperature range of 100-350 K. The water is modeled using two-body and three-body interaction potential (Stillinger-Weber potential) and the confinement is introduced using LJ 9-3 water-wall interaction potential. Entropy is calculated using RDF data obtained from the simulation experiments for each temperature point with increments or decrements of 2.5 K. The transition temperatures are estimated using the specific heat capacity analysis.
Collapse
Affiliation(s)
- Vikas Kumar Sinha
- Department of Chemical Engineering, National Institute of Technology Rourkela, Rourkela, Odisha, 769008, India
| | - Chandan Kumar Das
- Department of Chemical Engineering, National Institute of Technology Rourkela, Rourkela, Odisha, 769008, India.
| |
Collapse
|
2
|
Jiang J, Lai Y, Sheng D, Tang G, Zhang M, Niu D, Yu F. Two-dimensional bilayer ice in coexistence with three-dimensional ice without confinement. Nat Commun 2024; 15:5762. [PMID: 38982091 PMCID: PMC11233582 DOI: 10.1038/s41467-024-50187-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 06/26/2024] [Indexed: 07/11/2024] Open
Abstract
Icing plays an important role in various physical-chemical process. Although the formation of two-dimensional ice requires nanoscale confinement, two-dimensional bilayer ice in coexistence with three-dimensional ice without confinement remains poorly understood. Here, a critical value of a surface energy parameter is identified to characterize the liquid-solid interface interaction, above which two-dimensional and three-dimensional coexisting ice can surprisingly form on the surface. The two-dimensional ice growth mechanisms could be revealed by capturing the growth and merged of the metastable edge structures. The phase diagram about temperature and pressure vs energy parameters is predicted to distinguish liquid water, two-dimensional ice and three-dimensional ice. Furthermore, the deicing characteristics of coexisting ice demonstrate that the ice adhesion strength is linearly related to the ratio of ice-surface interaction energy to ice temperature. In addition, for gas-solid phase transition, the phase diagram about temperature and energy parameters is predicted to distinguish gas, liquid water, two-dimensional ice and three-dimensional ice. This work gives a perspective for studying the singular structure and dynamics of ice in nanoscale and provides a guide for future experimental realization of the coexisting ice.
Collapse
Affiliation(s)
- Jing Jiang
- State Key Laboratory of Frozen Soil Engineering, Northwest Institute of Eco-Environment and Resources, CAS, Lanzhou, PR China
| | - Yuanming Lai
- State Key Laboratory of Frozen Soil Engineering, Northwest Institute of Eco-Environment and Resources, CAS, Lanzhou, PR China.
- Institute of Future Civil Technology, Chongqing Jiaotong University, Chongqing, PR China.
| | - Daichao Sheng
- School of Civil and Environmental Engineering, University of Technology Sydney, Ultimo, NSW, Australia
| | - Guihua Tang
- MOE Key Laboratory of Thermo-Fluid Science and Engineering, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, PR China
| | - Mingyi Zhang
- State Key Laboratory of Frozen Soil Engineering, Northwest Institute of Eco-Environment and Resources, CAS, Lanzhou, PR China
| | - Dong Niu
- Naval Architecture and Ocean Engineering College, Dalian Maritime University, Dalian, PR China
| | - Fan Yu
- State Key Laboratory of Frozen Soil Engineering, Northwest Institute of Eco-Environment and Resources, CAS, Lanzhou, PR China
| |
Collapse
|
3
|
Wei L, Bai Q, Li X, Liu Z, Li C, Cui Y, Shen L, Zhu C, Fang W. Puckered Zigzag Monolayer Ice: Does a Confined Flat Four-Coordinated Monolayer Ice Always Have a Corresponding Puckered Phase? J Phys Chem Lett 2023; 14:8890-8895. [PMID: 37767947 DOI: 10.1021/acs.jpclett.3c02065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/29/2023]
Abstract
We note that a flat, four-coordinated monolayer ice under confinement always has a corresponding puckered phase. Recently, a monolayer ice consisting of an array of zigzag water chains (ZZMI) predicted by first-principles calculations of water under confinement is a flat four-coordinated monolayer ice. Herein, to investigate whether puckered ZZMI exists stably, we perform molecular dynamics simulations of two-dimensional (2D) ice formation for water constrained in graphene nanocapillaries. We find a novel monolayer ice structure that can be viewed as the ZZMI puckered along the direction perpendicular to the zigzag chain (pZZMI). Unlike ZZMI that does not satisfy the ice rule, each water molecule in pZZMI can form four hydrogen bonds (HBs) via forming two stable intersublayer HBs and two intrasublayer HBs. This work provides a fresh perspective on 2D confined ice, highlighting the intrinsic connections between 2D confined ices.
Collapse
Affiliation(s)
- Laiyang Wei
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Qi Bai
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Xiaojiao Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Ziyuan Liu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Chenruyuan Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Yanhong Cui
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, People's Republic of China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Chongqin Zhu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Weihai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| |
Collapse
|
4
|
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
|
5
|
Jyothirmai MV, Abraham BM, Singh JK. The pressure induced phase diagram of double-layer ice under confinement: a first-principles study. Phys Chem Chem Phys 2022; 24:16647-16654. [PMID: 35766352 DOI: 10.1039/d2cp01470j] [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/21/2022]
Abstract
Here, we present double-layer ice confined within various carbon nanotubes (CNTs) using state-of-the-art pressure induced (-5 GPa to 5 GPa) dispersion corrected density functional theory (DFT) calculations. We find that the double-layer ice exhibits remarkably rich and diverse phase behaviors as a function of pressure with varying CNT diameters. The lattice cohesive energies for various pure double layer ice phases follow the order of hexagonal > pentagonal > square tube > hexagonal-close-packed (HCP) > square > buckled-rhombic (b-RH). The confinement width was found to play a crucial role in the square and square tube phases in the intermediate pressure range of about 0-1 GPa. Unlike the phase transition in pure bilayer ice structures, the relative enthalpies demonstrate that the pentagonal phase, rather than the hexagonal structure, is the most stable ice polymorph at ambient pressure as well as in a deep negative pressure region, whereas the b-RH phase dominates under high pressure. The relatively short O⋯O distance of b-RH demonstrates the presence of a strong hydrogen bonding network, which is responsible for stabilizing the system.
Collapse
Affiliation(s)
- M V Jyothirmai
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, India.
| | - B Moses Abraham
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, India.
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, India. .,Prescience Insilico Private Limited, Bangalore 560049, India
| |
Collapse
|
6
|
|
7
|
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
|
8
|
Naullage PM, Metya AK, Molinero V. Computationally efficient approach for the identification of ice-binding surfaces and how they bind ice. J Chem Phys 2020; 153:174106. [DOI: 10.1063/5.0021631] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Pavithra M. Naullage
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Atanu K. Metya
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, USA
| |
Collapse
|
9
|
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
|
10
|
Naullage PM, Molinero V. Slow Propagation of Ice Binding Limits the Ice-Recrystallization Inhibition Efficiency of PVA and Other Flexible Polymers. J Am Chem Soc 2020; 142:4356-4366. [DOI: 10.1021/jacs.9b12943] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Pavithra M. Naullage
- 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
|
11
|
Zhu C, Gao Y, Zhu W, Jiang J, Liu J, Wang J, Francisco JS, Zeng XC. Direct observation of 2-dimensional ices on different surfaces near room temperature without confinement. Proc Natl Acad Sci U S A 2019; 116:16723-16728. [PMID: 31375634 PMCID: PMC6708332 DOI: 10.1073/pnas.1905917116] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Water-solid interfaces play important roles in a wide range of fields, including atmospheric science, geochemistry, electrochemistry, and food science. Herein, we report simulation evidence of 2-dimensional (2D) ice formation on various surfaces and the dependence of the 2D crystalline structure on the hydrophobicity and morphology of the underlying surface. Contrary to the prevailing view that nanoscale confinement is necessary for the 2D liquid-to-bilayer ice transition, we find that the liquid-to-bilayer hexagonal ice (BHI) transition can occur either on a model smooth surface or on model fcc-crystal surfaces with indices of (100), (110), and (111) near room temperature. We identify a critical parameter that characterizes the water-surface interaction, above which the BHI can form on the surface. This critical parameter increases as the temperature increases. Even at temperatures above the freezing temperature of bulk ice (Ih ), we find that BHI can also form on a superhydrophilic surface due to the strong water-surface interaction. The tendency toward the formation of BHI without confinement reflects a proper water-surface interaction that can compensate for the entropy loss during the freezing transition. Furthermore, phase diagrams of 2D ice formation are described on the plane of the adsorption energy versus the fcc lattice constant (Eads-afcc), where 4 monolayer square-like ices are also identified on the fcc model surfaces with distinct water-surface interactions.
Collapse
Affiliation(s)
- Chongqin Zhu
- Department of Earth and Environmental Science, University of Pennsylvania, Philadelphia, PA 19104
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, NE 68588
| | - Yurui Gao
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, NE 68588
| | - Weiduo Zhu
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, NE 68588
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, 230026 Anhui, China
- Department of Chemical Physics, University of Science and Technology of China, Hefei, 230026 Anhui, China
| | - Jian Jiang
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, NE 68588
| | - Jie Liu
- Key Laboratory of Green Printing, Institute of Chemistry, Chinese Academy of Sciences, 100190 Beijing, China
- School of Chemistry and Chemical Engineering, University of Chinese Academy of Sciences, 100190 Beijing, China
| | - Jianjun Wang
- Key Laboratory of Green Printing, Institute of Chemistry, Chinese Academy of Sciences, 100190 Beijing, China
- School of Chemistry and Chemical Engineering, University of Chinese Academy of Sciences, 100190 Beijing, China
| | - Joseph S Francisco
- Department of Earth and Environmental Science, University of Pennsylvania, Philadelphia, PA 19104;
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, NE 68588;
| |
Collapse
|
12
|
Yagasaki T, Yamasaki M, Matsumoto M, Tanaka H. Formation of hot ice caused by carbon nanobrushes. J Chem Phys 2019. [DOI: 10.1063/1.5111843] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masaru Yamasaki
- Graduate School of Natural Science and Technology, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
13
|
Hudait A, Qiu Y, Odendahl N, Molinero V. Hydrogen-Bonding and Hydrophobic Groups Contribute Equally to the Binding of Hyperactive Antifreeze and Ice-Nucleating Proteins to Ice. J Am Chem Soc 2019; 141:7887-7898. [DOI: 10.1021/jacs.9b02248] [Citation(s) in RCA: 56] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Arpa Hudait
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Yuqing Qiu
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Nathan Odendahl
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
14
|
Cao B, Xu E, Li T. Anomalous Stability of Two-Dimensional Ice Confined in Hydrophobic Nanopores. ACS NANO 2019; 13:4712-4719. [PMID: 30892864 DOI: 10.1021/acsnano.9b01014] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The freezing of water mostly proceeds via heterogeneous ice nucleation, a process in which an effective nucleation medium not only expedites ice crystallization but also may effectively direct the polymorph selection of ice. Here, we show that water confined within a hydrophobic slit nanopore exhibits a freezing behavior strongly distinguished from its bulk counterpart. Such a difference is reflected by a strong, non-monotonic pore-size dependence of freezing temperature but, more surprisingly, by an unexpected stacking ordering of crystallized two-dimensional ice containing just a few ice layers. In particular, confined trilayer ice is found to exclusively crystallize into a well-ordered, hexagonal stacking sequence despite the fact that nanopore exerts no explicit constraint on stacking order. The absence of cubic stacking sequence is found to be originated from the intrinsically lower thermodynamic stability of cubic ice over hexagonal ice at the interface, which contrasts sharply the nearly degenerated stability of bulk hexagonal and cubic ices. Detailed examination clearly reveals that the divergence is attributed to the inherent difference between the two ice polymorphs in their surface phonon modes, which is further found to generically occur at both hydrophobic and hydrophilic surfaces.
Collapse
Affiliation(s)
- Boxiao Cao
- Department of Civil and Environmental Engineering , George Washington University , Washington , D.C. 20052 , United States
| | - Enshi Xu
- Department of Civil and Environmental Engineering , George Washington University , Washington , D.C. 20052 , United States
| | - Tianshu Li
- Department of Civil and Environmental Engineering , George Washington University , Washington , D.C. 20052 , United States
| |
Collapse
|
15
|
Qiu Y, Hudait A, Molinero V. How Size and Aggregation of Ice-Binding Proteins Control Their Ice Nucleation Efficiency. J Am Chem Soc 2019; 141:7439-7452. [DOI: 10.1021/jacs.9b01854] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Yuqing Qiu
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0580, United States
| | - Arpa Hudait
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0580, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0580, United States
| |
Collapse
|
16
|
Franco Pinheiro Moreira PA, Gomes de Aguiar Veiga R, de Koning M. Elastic constants of ice Ihas described by semi-empirical water models. J Chem Phys 2019; 150:044503. [DOI: 10.1063/1.5082743] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
| | - Roberto Gomes de Aguiar Veiga
- Centro de Engenharia, Modelagem e Ciências Sociais Aplicadas, Universidade Federal do ABC, Santo André, 09210-580 São Paulo, Brazil
| | - Maurice de Koning
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, UNICAMP, Campinas, 13083-859 São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, Campinas, 13083-861 São Paulo, Brazil
| |
Collapse
|
17
|
Santos-Flórez PA, Ruestes CJ, de Koning M. Uniaxial-deformation behavior of ice Ih as described by the TIP4P/Ice and mW water models. J Chem Phys 2018; 149:164711. [DOI: 10.1063/1.5048517] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Pedro Antonio Santos-Flórez
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
| | - Carlos J. Ruestes
- CONICET and Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Mendoza, Argentina
| | - Maurice de Koning
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, UNICAMP, 13083-859 Campinas, São Paulo, Brazil
- Center for Computational Engineering and Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861 Campinas, São Paulo, Brazil
| |
Collapse
|
18
|
Abstract
Nanoscale confinement has a strong effect on the phase behavior of water. Studies in the last two decades have revealed a wealth of novel crystalline and quasicrystalline structures for water confined in nanoslits. Less is known, however, about the nature of ice-liquid coexistence in extremely nanoconfined systems. Here, we use molecular simulations to investigate the ice-liquid equilibrium for water confined between two nanoscopic disks. We find that the nature of ice-liquid phase coexistence in nanoconfined water is different from coexistence in both bulk water and extended nanoslits. In highly nanoconfined systems, liquid water and ice do not coexist in space because the two-phase states are unstable. The confined ice and liquid phases coexist in time, through oscillations between all-liquid and all-crystalline states. The avoidance of spatial coexistence of ice and liquid originates on the non-negligible cost of the interface between confined ice and liquid in a small system. It is the result of the small number of water molecules between the plates and has no analogue in bulk water.
Collapse
Affiliation(s)
- Noah Kastelowitz
- Department of Chemistry , The University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112-0850 , United States
| | - Valeria Molinero
- Department of Chemistry , The University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112-0850 , United States
| |
Collapse
|
19
|
Kaneko T. Elevation/depression mechanism of freezing points of liquid confined in slit nanopores. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1350785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Toshihiro Kaneko
- Department of Mechanical Engineering, Tokyo University of Science, Noda, Japan
- Research Institute for Science and Technology (RIST), Tokyo University of Science, Noda, Japan
| |
Collapse
|
20
|
Shevkunov SV. Water-vapor clustering on the surface of β-AgI crystal in the field of defects with a disordered structure. COLLOID JOURNAL 2017. [DOI: 10.1134/s1061933x1705012x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
21
|
Ferrara CG, Grigera TS. Dynamics and structural behavior of water in large confinement with planar amorphous walls. J Chem Phys 2017; 147:024705. [PMID: 28711040 DOI: 10.1063/1.4991834] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study the structure and dynamics of liquid water confined between planar amorphous walls using molecular dynamics (MD) simulations. We report MD results for systems of more than 23 000 SPC/E water molecules confined between two hydrophilic or hydrophobic walls, separated by distances of about 15 nm. We find that the walls induce ordering of the liquid and slow down the dynamics, affecting the properties of the confined water up to distances of about 8 nm at 275 K. We quantify this influence by computing dynamic and static penetration lengths and studying their temperature dependence. Our results indicate that in the temperature range considered, hydrophobic walls perturb static properties over larger lengths compared to hydrophilic walls. We also find opposite temperature trends in the dynamic penetration lengths, with hydrophobic walls increasing their range of influence on increasing the temperature.
Collapse
Affiliation(s)
- C Gastón Ferrara
- Universidad Nacional Arturo Jauretche, Florencio Varela, Argentina
| | - Tomás S Grigera
- Instituto de Física de Líquidos y Sistemas Biológicos (IFLYSIB), CONICET and Facultad de Ciencias Exactas, Universidad Nacional de La Plata, Calle 59 No. 789, B1900BTE La Plata, Argentina
| |
Collapse
|
22
|
Strength of Alkane–Fluid Attraction Determines the Interfacial Orientation of Liquid Alkanes and Their Crystallization through Heterogeneous or Homogeneous Mechanisms. CRYSTALS 2017. [DOI: 10.3390/cryst7030086] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
23
|
Qiu Y, Odendahl N, Hudait A, Mason R, Bertram AK, Paesani F, DeMott PJ, Molinero V. Ice Nucleation Efficiency of Hydroxylated Organic Surfaces Is Controlled by Their Structural Fluctuations and Mismatch to Ice. J Am Chem Soc 2017; 139:3052-3064. [PMID: 28135412 DOI: 10.1021/jacs.6b12210] [Citation(s) in RCA: 112] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Heterogeneous nucleation of ice induced by organic materials is of fundamental importance for climate, biology, and industry. Among organic ice-nucleating surfaces, monolayers of long chain alcohols are particularly effective, while monolayers of fatty acids are significantly less so. As these monolayers expose to water hydroxyl groups with an order that resembles the one in the basal plane of ice, it was proposed that lattice matching between ice and the surface controls their ice-nucleating efficiency. Organic monolayers are soft materials and display significant fluctuations. It has been conjectured that these fluctuations assist in the nucleation of ice. Here we use molecular dynamic simulations and laboratory experiments to investigate the relationship between the structure and fluctuations of hydroxylated organic surfaces and the temperature at which they nucleate ice. We find that these surfaces order interfacial water to form domains with ice-like order that are the birthplace of ice. Both mismatch and fluctuations decrease the size of the preordered domains and monotonously decrease the ice freezing temperature. The simulations indicate that fluctuations depress the freezing efficiency of monolayers of alcohols or acids to half the value predicted from lattice mismatch alone. The model captures the experimental trend in freezing efficiencies as a function of chain length and predicts that alcohols have higher freezing efficiency than acids of the same chain length. These trends are mostly controlled by the modulation of the structural mismatch to ice. We use classical nucleation theory to show that the freezing efficiencies of the monolayers are directly related to their free energy of binding to ice. This study provides a general framework to relate the equilibrium thermodynamics of ice binding to a surface and the nonequilibrium ice freezing temperature and suggests that these could be predicted from the structure of interfacial water.
Collapse
Affiliation(s)
- Yuqing Qiu
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Nathan Odendahl
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Arpa Hudait
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Ryan Mason
- Department of Chemistry, University of British Columbia , Vancouver, British Columbia V6T 1Z1, Canada
| | - Allan K Bertram
- Department of Chemistry, University of British Columbia , Vancouver, British Columbia V6T 1Z1, Canada
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92093, United States
| | - Paul J DeMott
- Department of Atmospheric Science, Colorado State University , Fort Collins, Colorado 80523-1371, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
24
|
Lupi L, Peters B, Molinero V. Pre-ordering of interfacial water in the pathway of heterogeneous ice nucleation does not lead to a two-step crystallization mechanism. J Chem Phys 2016; 145:211910. [DOI: 10.1063/1.4961652] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Laura Lupi
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| | - Baron Peters
- Department of Chemical Engineering and Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| |
Collapse
|
25
|
Lu J, Jacobson LC, Perez Sirkin YA, Molinero V. High-Resolution Coarse-Grained Model of Hydrated Anion-Exchange Membranes that Accounts for Hydrophobic and Ionic Interactions through Short-Ranged Potentials. J Chem Theory Comput 2016; 13:245-264. [PMID: 28068769 DOI: 10.1021/acs.jctc.6b00874] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Jibao Lu
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Liam C. Jacobson
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Yamila A. Perez Sirkin
- Departamento
de Química Inorgánica, Analítica y Química
Física, and INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires C1428EHA, Argentina
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
26
|
Zhu Y, Wang F, Wu H. Buckling failure of square ice-nanotube arrays constrained in graphene nanocapillaries. J Chem Phys 2016; 145:054704. [PMID: 27497569 DOI: 10.1063/1.4959902] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Graphene confinement provides a new physical and mechanical environment with ultrahigh van der Waals pressure, resulting in new quasi-two-dimensional phases of few-layer ice. Polymorphic transition can occur in bilayer constrained water/ice system. Here, we perform a comprehensive study of the phase transition of AA-stacked bilayer water constrained within a graphene nanocapillary. The compression-limit and superheating-limit (phase) diagrams are obtained, based on the extensive molecular-dynamics simulations at numerous thermodynamic states. Liquid-to-solid, solid-to-solid, and solid-to-liquid-to-solid phase transitions are observed in the compression and superheating of bilayer water. Interestingly, there is a temperature threshold (∼275 K) in the compression-limit diagram, which indicates that the first-order and continuous-like phase transitions of bilayer water depend on the temperature. Two obviously different physical processes, compression and superheating, display similar structural evolution; that is, square ice-nanotube arrays (BL-VHDI) will bend first and then transform into bilayer triangular AA stacking ice (BL-AAI). The superheating limit of BL-VHDI exhibits local maxima, while that of BL-AAI increases monotonically. More importantly, from a mechanics point of view, we propose a novel mechanism of the transformation from BL-VHDI to BL-AAI, both for the compression and superheating limits. This structural transformation can be regarded as the "buckling failure" of the square-ice-nanotube columns, which is dominated by the lateral pressure.
Collapse
Affiliation(s)
- YinBo Zhu
- CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei, Anhui 230027, China
| | - FengChao Wang
- CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei, Anhui 230027, China
| | - HengAn Wu
- CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei, Anhui 230027, China
| |
Collapse
|
27
|
Abstract
Ice crystals in the atmosphere nucleate from supercooled liquid water and grow by vapor uptake. The structure of the ice polymorph grown has strong impact on the morphology and light scattering of the ice crystals, modulates the amount of water vapor in ice clouds, and can impact the molecular uptake and reactivity of atmospheric aerosols. Experiments and molecular simulations indicate that ice nucleated and grown from deeply supercooled liquid water is metastable stacking disordered ice. The ice polymorph grown from vapor has not yet been determined. Here we use large-scale molecular simulations to determine the structure of ice that grows as a result of uptake of water vapor in the temperature range relevant to cirrus and mixed-phase clouds, elucidate the molecular mechanism of the formation of ice at the vapor interface, and compute the free energy difference between cubic and hexagonal ice interfaces with vapor. We find that vapor deposition results in growth of stacking disordered ice only under conditions of extreme supersaturation, for which a nonequilibrium liquid layer completely wets the surface of ice. Such extreme conditions have been used to produce stacking disordered frost ice in experiments and may be plausible in the summer polar mesosphere. Growth of ice from vapor at moderate supersaturations in the temperature range relevant to cirrus and mixed-phase clouds, from 200 to 260 K, produces exclusively the stable hexagonal ice polymorph. Cubic ice is disfavored with respect to hexagonal ice not only by a small penalty in the bulk free energy (3.6 ± 1.5 J mol(-1) at 260 K) but also by a large free energy penalty at the ice-vapor interface (89.7 ± 12.8 J mol(-1) at 260 K). The latter originates in higher vibrational entropy of the hexagonal-terminated ice-vapor interface. We predict that the free energy penalty against the cubic ice interface should decrease strongly with temperature, resulting in some degree of stacking disorder in ice grown from vapor in the tropical tropopause layer, and in polar stratospheric and noctilucent clouds. Our findings support and explain the evolution of the morphology of ice crystals from hexagonal to trigonal symmetry with decreasing temperature, as reported by experiments and in situ measurements in clouds. We conclude that selective growth of the elusive cubic ice polymorph by manipulation of the interfacial properties can likely be achieved at the ice-liquid interface but not at the ice-vapor interface.
Collapse
Affiliation(s)
- Arpa Hudait
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
28
|
Lu J, Chakravarty C, Molinero V. Relationship between the line of density anomaly and the lines of melting, crystallization, cavitation, and liquid spinodal in coarse-grained water models. J Chem Phys 2016; 144:234507. [DOI: 10.1063/1.4953854] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jibao Lu
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, USA
| | | | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, USA
| |
Collapse
|
29
|
Pérez-Díaz JL, Álvarez-Valenzuela MA, Rodríguez-Celis F. Surface freezing of water. SPRINGERPLUS 2016; 5:629. [PMID: 27330895 PMCID: PMC4870584 DOI: 10.1186/s40064-016-2196-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/05/2015] [Accepted: 04/20/2016] [Indexed: 11/10/2022]
Abstract
Freezing, melting, evaporation and condensation of water are essential ingredients for climate and eventually life on Earth. In the present work, we show how surface freezing of supercooled water in an open container is conditioned and triggered-exclusively-by humidity in air. Additionally, a change of phase is demonstrated to be triggered on the water surface forming surface ice crystals prior to freezing of bulk. The symmetry of the surface crystal, as well as the freezing point, depend on humidity, presenting at least three different types of surface crystals. Humidity triggers surface freezing as soon as it overpasses a defined value for a given temperature, generating a plurality of nucleation nodes. An evidence of simultaneous nucleation of surface ice crystals is also provided.
Collapse
Affiliation(s)
- J. L. Pérez-Díaz
- />Departamento de Teoría de la Señal, Universidad de Alcalá, EPS, N-II km 33,600, 28801 Alcalá de Henares, Spain
- />Division of Space Technology, Department of Computer Science, Electrical and Space Engineering, Luleå University of Technology, Kiruna, Sweden
| | | | | |
Collapse
|
30
|
Foroutan M, Fatemi SM, Shokouh F. Graphene confinement effects on melting/freezing point and structure and dynamics behavior of water. J Mol Graph Model 2016; 66:85-90. [DOI: 10.1016/j.jmgm.2016.03.011] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2015] [Revised: 02/11/2016] [Accepted: 03/24/2016] [Indexed: 10/22/2022]
|
31
|
Hudait A, Qiu S, Lupi L, Molinero V. Free energy contributions and structural characterization of stacking disordered ices. Phys Chem Chem Phys 2016; 18:9544-53. [PMID: 26983558 DOI: 10.1039/c6cp00915h] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Crystallization of ice from deeply supercooled water and amorphous ices - a process of fundamental importance in the atmosphere, interstellar space, and cryobiology - results in stacking disordered ices with a wide range of metastabilities with respect to hexagonal ice. The structural origin of this high variability, however, has not yet been elucidated. Here we use molecular dynamics simulations with the mW water model to characterize the structure of ice freshly grown from supercooled water at temperatures from 210 to 270 K, the thermodynamics of stacking faults, line defects, and interfaces, and to elucidate the interplay between kinetics and thermodynamics in determining the structure of ice. In agreement with experiments, the ice grown in the simulations is stacking disordered with a random distribution of cubic and hexagonal layers, and a cubicity that decreases with growth temperature. The former implies that the cubicity of ice is determined by processes at the ice/liquid interface, without memory of the structure of buried ice layers. The latter indicates that the probability of building a cubic layer at the interface decreases upon approaching the melting point of ice, which we attribute to a more efficient structural equilibration of ice at the liquid interface as the driving force for growth wanes. The free energy cost for creating a pair of cubic layers in ice is 8.0 J mol(-1) in experiments, and 9.7 ± 1.9 J mol(-1) for the mW water model. This not only validates the simulations, but also indicates that dispersion in cubicity is not sufficient to explain the large energetic variability of stacking disordered ices. We compute the free energy cost of stacking disorder, line defects, and interfaces in ice and conclude that a characterization of the density of these defects is required to predict the degree of metastability and vapor pressure of atmospheric ices.
Collapse
Affiliation(s)
- Arpa Hudait
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA.
| | | | | | | |
Collapse
|
32
|
Lu Q, Straub JE. Freezing Transitions of Nanoconfined Coarse-Grained Water Show Subtle Dependence on Confining Environment. J Phys Chem B 2016; 120:2517-25. [DOI: 10.1021/acs.jpcb.5b10481] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Qing Lu
- Division
of Materials Science and Engineering, Boston University, Brookline, Massachusetts 02446, United States
| | - John E. Straub
- Department
of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
33
|
Pant S, Ghorai PK. Structural anomaly of core-softened fluid confined in single walled carbon nanotube: a molecular dynamics simulation investigation. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1149242] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Shashank Pant
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| | - Pradip K. Ghorai
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| |
Collapse
|
34
|
Chen J, Schusteritsch G, Pickard CJ, Salzmann CG, Michaelides A. Two Dimensional Ice from First Principles: Structures and Phase Transitions. PHYSICAL REVIEW LETTERS 2016; 116:025501. [PMID: 26824547 DOI: 10.1103/physrevlett.116.025501] [Citation(s) in RCA: 92] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/18/2015] [Indexed: 06/05/2023]
Abstract
Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here, we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression, the pentagonal structure becomes the most stable and persists up to ∼2 GPa, at which point the square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and the width.
Collapse
Affiliation(s)
- Ji Chen
- London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Georg Schusteritsch
- Thomas Young Centre, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, United Kingdom
| | - Chris J Pickard
- Thomas Young Centre, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, United Kingdom
| | - Christoph G Salzmann
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
35
|
Zhu Y, Wang F, Bai J, Zeng XC, Wu H. AB-stacked square-like bilayer ice in graphene nanocapillaries. Phys Chem Chem Phys 2016; 18:22039-46. [DOI: 10.1039/c6cp03061k] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Water, when constrained between two graphene sheets and under ultrahigh pressure, can manifest dramatic differences from its bulk counterparts such as the van der Waals pressure induced water-to-ice transformation, known as the metastability limit of two-dimensional (2D) liquid.
Collapse
Affiliation(s)
- YinBo Zhu
- CAS Key Laboratory of Mechanical Behavior and Design of Materials
- Department of Modern Mechanics
- CAS Center for Excellence in Nanoscience
- University of Science and Technology of China
- Hefei
| | - FengChao Wang
- CAS Key Laboratory of Mechanical Behavior and Design of Materials
- Department of Modern Mechanics
- CAS Center for Excellence in Nanoscience
- University of Science and Technology of China
- Hefei
| | - Jaeil Bai
- Department of Chemistry
- University of Nebraska-Lincoln
- USA
| | - Xiao Cheng Zeng
- Department of Chemistry
- University of Nebraska-Lincoln
- USA
- Hefei National Laboratory for Physical Sciences at Microscale and Collaborative Innovation Center of Chemistry for Energy Materials
- University of Science and Technology of China
| | - HengAn Wu
- CAS Key Laboratory of Mechanical Behavior and Design of Materials
- Department of Modern Mechanics
- CAS Center for Excellence in Nanoscience
- University of Science and Technology of China
- Hefei
| |
Collapse
|
36
|
Zhu W, Zhao WH, Wang L, Yin D, Jia M, Yang J, Zeng XC, Yuan LF. Two-dimensional interlocked pentagonal bilayer ice: how do water molecules form a hydrogen bonding network? Phys Chem Chem Phys 2016; 18:14216-21. [DOI: 10.1039/c5cp07524f] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The tradeoff between the conditions of an ideal hydrogen bonding network can serve as a generic guidance to understand the rich phase behaviors of nanoconfined water.
Collapse
Affiliation(s)
- Weiduo Zhu
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Wen-Hui Zhao
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Lu Wang
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Di Yin
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Min Jia
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Xiao Cheng Zeng
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| | - Lan-Feng Yuan
- Hefei National Laboratory for Physical Sciences at Microscale
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei
- China
| |
Collapse
|
37
|
Zhu Y, Wang F, Bai J, Zeng XC, Wu H. Compression Limit of Two-Dimensional Water Constrained in Graphene Nanocapillaries. ACS NANO 2015; 9:12197-204. [PMID: 26575824 DOI: 10.1021/acsnano.5b06572] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Evaluation of the tensile/compression limit of a solid under conditions of tension or compression is often performed to provide mechanical properties that are critical for structure design and assessment. Algara-Siller et al. recently demonstrated that when water is constrained between two sheets of graphene, it becomes a two-dimensional (2D) liquid and then is turned into an intriguing monolayer solid with a square pattern under high lateral pressure [ Nature , 2015 , 519 , 443 - 445 ]. From a mechanics point of view, this liquid-to-solid transformation characterizes the compression limit (or metastability limit) of the 2D monolayer water. Here, we perform a simulation study of the compression limit of 2D monolayer, bilayer, and trilayer water constrained in graphene nanocapillaries. At 300 K, a myriad of 2D ice polymorphs (both crystalline-like and amorphous) are formed from the liquid water at different widths of the nanocapillaries, ranging from 6.0 to11.6 Å. For monolayer water, the compression limit is typically a few hundred MPa, while for the bilayer and trilayer water, the compression limit is 1.5 GPa or higher, reflecting the ultrahigh van der Waals pressure within the graphene nanocapillaries. The compression-limit (phase) diagram is obtained at the nanocapillary width versus pressure (h-P) plane, based on the comprehensive molecular dynamics simulations at numerous thermodynamic states as well as on the Clapeyron equation. Interestingly, the compression-limit curves exhibit multiple local minima.
Collapse
Affiliation(s)
- YinBo Zhu
- CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, University of Science and Technology of China , Hefei, Anhui 230027, China
| | - FengChao Wang
- CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, University of Science and Technology of China , Hefei, Anhui 230027, China
| | - Jaeil Bai
- Department of Chemistry, University of Nebraska-Lincoln , Lincoln, Nebraska 68588, United States
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska-Lincoln , Lincoln, Nebraska 68588, United States
- Hefei National Laboratory for Physical Sciences at Microscale and Collaborative Innovation Center of Chemistry for Energy Materials, University of Science and Technology of China , Hefei, Anhui 230026, China
| | - HengAn Wu
- CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, University of Science and Technology of China , Hefei, Anhui 230027, China
| |
Collapse
|
38
|
Kaneko T, Bai J, Yasuoka K, Mitsutake A, Zeng XC. Liquid-solid and solid-solid phase transition of monolayer water: high-density rhombic monolayer ice. J Chem Phys 2015; 140:184507. [PMID: 24832288 DOI: 10.1063/1.4874696] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Liquid-solid and solid-solid phase transitions of a monolayer water confined between two parallel hydrophobic surfaces are studied by molecular dynamics simulations. The solid phase considered is the high-density rhombic monolayer ice. Based on the computed free energy surface, it is found that at a certain width of the slit nanopore, the monolayer water exhibits not only a high freezing point but also a low energy barrier to crystallization. Moreover, through analyzing the oxygen-hydrogen-oxygen angle distribution and oxygen-hydrogen radial distribution, the high-density monolayer ice is classified as either a flat ice or a puckered ice. The transition between a flat ice and a puckered ice reflects a trade-off between the water-wall interactions and the electrostatic interactions among water molecules.
Collapse
Affiliation(s)
- Toshihiro Kaneko
- Department of Mechanical Engineering, Tokyo University of Science, Noda 278-8510, Japan
| | - Jaeil Bai
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, Yokohama 223-8522, Japan
| | - Ayori Mitsutake
- Department of Physics, Keio University, Yokohama 223-8522, Japan and JST, PRESTO, Yokohama 223-8522, Japan
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| |
Collapse
|
39
|
Koschke K, Limbach HJ, Kremer K, Donadio D. Freezing point depression in model Lennard-Jones solutions. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1029029] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
| | | | - Kurt Kremer
- Max Planck Institute for Polymer Research , Mainz, Germany
| | - Davide Donadio
- Max Planck Institute for Polymer Research , Mainz, Germany
- Donostia International Physics Center , San Sebastian, Spain
- IKERBASQUE, Basque Foundation for Science , Bilbao, Spain
| |
Collapse
|
40
|
Mei F, Zhou X, Kou J, Wu F, Wang C, Lu H. A transition between bistable ice when coupling electric field and nanoconfinement. J Chem Phys 2015; 142:134704. [DOI: 10.1063/1.4916521] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Feng Mei
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China
- Department of Physics, Zhejiang University, Hangzhou 310027, China
| | - Xiaoyan Zhou
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China
- Department of Physics and Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China
| | - Jianlong Kou
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China
| | - Fengmin Wu
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China
- Department of Physics and Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China
| | - Chunlei Wang
- Division of Interfacial Water and Key Laboratory of Interfacial Physics and Technology, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China
| | - Hangjun Lu
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China
| |
Collapse
|
41
|
Bullock G, Molinero V. Low-density liquid water is the mother of ice: on the relation between mesostructure, thermodynamics and ice crystallization in solutions. Faraday Discuss 2015; 167:371-88. [PMID: 24640501 DOI: 10.1039/c3fd00085k] [Citation(s) in RCA: 83] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Predicting the temperature and extent of ice freezing in aqueous solutions is crucial for areas as diverse as cryobiology and materials design. It has long been recognized that the thermodynamics of liquid water controls the temperature and kinetics of ice crystallization. Parameterizations of the freezing temperatures in terms of the water activity of the solution have been successfully established, but the fundamental origin of the thermodynamic control of the non-equilibrium crystallization of ice has remained elusive. Here we use large-scale molecular simulations to elucidate the relationship between the structure, thermodynamics, and ice crystallization temperatures for solutions of mW water and a strongly hydrophilic solute that mimics LiCI ions. Fast cooling of solutions with up to 20 mol% ions results in the formation of nanosegregated glasses with domains of low-density amorphous ice and an ion-rich vitrified solution. Slow cooling of the mixtures results in nucleation and growth of ice within the domains of four-coordinated liquid water. The temperature of crystallization Tx coincides with the temperature of appearance of nanoscopic domains of four-coordinated liquid water in the mixture, T(L). We use the insight provided by the simulations to derive a thermodynamic expression for the crystallization temperature as a function of the water activity, T(X)(a(W)), analogous to the dependence of the melting temperature, T(m)(a(W)). The simple expression derived in this work provides a good account of the experimental freezing temperatures of water and the well-known steepest dependence of Tx on solute concentration compared to that of T(m).
Collapse
|
42
|
Lu Q, Kim J, Farrell JD, Wales DJ, Straub JE. Investigating the solid-liquid phase transition of water nanofilms using the generalized replica exchange method. J Chem Phys 2014; 141:18C525. [DOI: 10.1063/1.4896513] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qing Lu
- Division of Materials Science and Engineering, Boston University, Brookline, Massachusetts 02446, USA
| | - Jaegil Kim
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - James D. Farrell
- University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - David J. Wales
- University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - John E. Straub
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| |
Collapse
|
43
|
Lupi L, Kastelowitz N, Molinero V. Vapor deposition of water on graphitic surfaces: Formation of amorphous ice, bilayer ice, ice I, and liquid water. J Chem Phys 2014; 141:18C508. [DOI: 10.1063/1.4895543] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Laura Lupi
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| | - Noah Kastelowitz
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| |
Collapse
|
44
|
Zhao WH, Wang L, Bai J, Yuan LF, Yang J, Zeng XC. Highly confined water: two-dimensional ice, amorphous ice, and clathrate hydrates. Acc Chem Res 2014; 47:2505-13. [PMID: 25088018 DOI: 10.1021/ar5001549] [Citation(s) in RCA: 76] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding phase behavior of highly confined water, ice, amorphous ice, and clathrate hydrates (or gas hydrates), not only enriches our view of phase transitions and structures of quasi-two-dimensional (Q2D) solids not seen in the bulk phases but also has important implications for diverse phenomena at the intersection between physical chemistry, cell biology, chemical engineering, and nanoscience. Relevant examples include, among others, boundary lubrication in nanofluidic and lab-on-a-chip devices, synthesis of antifreeze proteins for ice-growth inhibition, rapid cooling of biological suspensions or quenching emulsified water under high pressure, and storage of H2 and CO2 in gas hydrates. Classical molecular simulation (MD) is an indispensable tool to explore states and properties of highly confined water and ice. It also has the advantage of precisely monitoring the time and spatial domains in the sub-picosecond and sub-nanometer scales, which are difficult to control in laboratory experiments, and yet allows relatively long simulation at the 10(2) ns time scale that is impractical with ab initio molecular dynamics simulations. In this Account, we present an overview of our MD simulation studies of the structures and phase behaviors of highly confined water, ice, amorphous ice, and clathrate, in slit graphene nanopores. We survey six crystalline phases of monolayer (ML) ice revealed from MD simulations, including one low-density, one mid-density, and four high-density ML ices. We show additional supporting evidence on the structural stabilities of the four high-density ML ices in the vacuum (without the graphene confinement), for the first time, through quantum density-functional theory optimization of their free-standing structures at zero temperature. In addition, we summarize various low-density, high-density, and very-high-density Q2D bilayer (BL) ice and amorphous ice structures revealed from MD simulations. These simulations reinforce the notion that the nanoscale confinement not only can disrupt the hydrogen bonding network in bulk water but also can allow satisfaction of the ice rule for low-density and high-density Q2D crystalline structures. Highly confined water can serve as a generic model system for understanding a variety of Q2D materials science phenomena, for example, liquid-solid, solid-solid, solid-amorphous, and amorphous-amorphous transitions in real time, as well as the Ostwald staging during these transitions. Our simulations also bring new molecular insights into the formation of gas hydrate from a gas and water mixture at low temperature.
Collapse
Affiliation(s)
- Wen-Hui Zhao
- Hefei
National Laboratory for Physical Sciences at Microscale and Collaborative
Innovation Center of Chemistry for Energy Materials, Department of
Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Lu Wang
- Hefei
National Laboratory for Physical Sciences at Microscale and Collaborative
Innovation Center of Chemistry for Energy Materials, Department of
Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jaeil Bai
- Department
of Chemistry, University of Nebraska—Lincoln, Lincoln, Nebraska 68588, United States
| | - Lan-Feng Yuan
- Hefei
National Laboratory for Physical Sciences at Microscale and Collaborative
Innovation Center of Chemistry for Energy Materials, Department of
Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jinlong Yang
- Hefei
National Laboratory for Physical Sciences at Microscale and Collaborative
Innovation Center of Chemistry for Energy Materials, Department of
Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiao Cheng Zeng
- Hefei
National Laboratory for Physical Sciences at Microscale and Collaborative
Innovation Center of Chemistry for Energy Materials, Department of
Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department
of Chemistry, University of Nebraska—Lincoln, Lincoln, Nebraska 68588, United States
| |
Collapse
|
45
|
Lu J, Qiu Y, Baron R, Molinero V. Coarse-Graining of TIP4P/2005, TIP4P-Ew, SPC/E, and TIP3P to Monatomic Anisotropic Water Models Using Relative Entropy Minimization. J Chem Theory Comput 2014; 10:4104-20. [PMID: 26588552 DOI: 10.1021/ct500487h] [Citation(s) in RCA: 91] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Coarse-grained models are becoming a competitive alternative for modeling processes that occur over time and length scales beyond the reach of fully atomistic molecular simulations. Ideally, coarse-grained models should not only achieve high computational efficiency but also provide accurate predictions and fundamental insight into the role of molecular interactions, the characteristic behavior, and properties of the system they model. In this work we derive a series of monatomic coarse-grained water models mX(REM) from the most popular atomistic water models X = TIP3P, SPC/E, TIP4P-Ew, and TIP4P/2005, using the relative entropy minimization (REM) method. Each coarse-grained water molecule is represented by a single particle that interacts through short-ranged anisotropic interactions that encourage the formation of "hydrogen-bonded" structures. We systematically investigate the features of the coarse-grained models in reproducing over 20 structural, dynamic, and thermodynamic properties of the reference atomistic water models-including the existence and locus of the characteristic density anomaly. The mX(REM) coarse-grained models reproduce quite faithfully the radial and angular distribution function of water, produce a temperature of maximum density (TMD), and stabilize the ice I crystal. Moreover, the ratio between the TMD and the melting temperature of the crystal in the mX(REM) models and liquid-ice equilibrium properties show reasonable agreement with the results of the corresponding atomistic models. The mX(REM) models, however, severely underestimate the cohesive energy of the condensed water phases. We investigate which specific limitations of the coarse-grained models arise from the REM methodology, from the monatomic nature of the models, and from the Stillinger-Weber interaction potential form. Our analysis indicates that a small compromise in the accuracy of structural properties can result in a significant increase of the overall accuracy and representability of the coarse-grained water models. We evaluate the accuracy of the atomistic and the monatomic anisotropic coarse-grained water models, including the mW water model, in reproducing experimental water properties. We find that mW and mTIP4P/2005(REM) score closer to experiment than widely used atomistic water models. We conclude that monatomic models of water with short-range, anisotropic "hydrogen-bonding" three-body interactions can be competitive in accuracy with fully atomistic models for the study of a wide range of properties and phenomena at less than 1/100th of the computational cost.
Collapse
Affiliation(s)
- Jibao Lu
- Department of Chemistry, The University of Utah , Salt Lake City, Utah 84112-0850, United States.,Department of Medicinal Chemistry, The University of Utah , Salt Lake City, Utah 84112-5820, United States
| | - Yuqing Qiu
- Department of Chemistry, The University of Utah , Salt Lake City, Utah 84112-0850, United States
| | - Riccardo Baron
- Department of Medicinal Chemistry, The University of Utah , Salt Lake City, Utah 84112-5820, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah , Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
46
|
Hudait A, Molinero V. Ice Crystallization in Ultrafine Water–Salt Aerosols: Nucleation, Ice-Solution Equilibrium, and Internal Structure. J Am Chem Soc 2014; 136:8081-93. [DOI: 10.1021/ja503311r] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Arpa Hudait
- Department
of Chemistry, The University of Utah, 315 South 1400 East, Salt
Lake City, Utah 84112-0850, United States
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, 315 South 1400 East, Salt
Lake City, Utah 84112-0850, United States
| |
Collapse
|
47
|
Anick DJ. Static Density Functional Study of Graphene–Hexagonal Bilayer Ice Interaction. J Phys Chem A 2014; 118:7498-506. [PMID: 24641236 DOI: 10.1021/jp500360n] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David J. Anick
- Laboratory for Water and
Surface Studies Department of Chemistry, Tufts University, 62
Pearson Avenue, Medford, Massachusetts 02155, Unites States
| |
Collapse
|
48
|
Jia M, Zhao WH, Yuan LF. New Hexagonal-rhombic Trilayer Ice Structure Confined between Hydrophobic Plates. CHINESE J CHEM PHYS 2014. [DOI: 10.1063/1674-0068/27/01/15-19] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
49
|
Lupi L, Hudait A, Molinero V. Heterogeneous nucleation of ice on carbon surfaces. J Am Chem Soc 2014; 136:3156-64. [PMID: 24495074 DOI: 10.1021/ja411507a] [Citation(s) in RCA: 213] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Atmospheric aerosols can promote the heterogeneous nucleation of ice, impacting the radiative properties of clouds and Earth's climate. The experimental investigation of heterogeneous freezing of water droplets by carbonaceous particles reveals widespread ice freezing temperatures. It is not known which structural and chemical characteristics of soot account for the variability in ice nucleation efficiency. Here we use molecular dynamics simulations to investigate the nucleation of ice from liquid water in contact with graphitic surfaces. We find that atomically flat carbon surfaces promote heterogeneous nucleation of ice, while molecularly rough surfaces with the same hydrophobicity do not. Graphitic surfaces and other surfaces that promote ice nucleation induce layering in the interfacial water, suggesting that the order imposed by the surface on liquid water may play an important role in the heterogeneous nucleation mechanism. We investigate a large set of graphitic surfaces of various dimensions and radii of curvature and find that variations in nanostructures alone could account for the spread in the freezing temperatures of ice on soot in experiments. We conclude that a characterization of the nanostructure of soot is needed to predict its ice nucleation efficiency.
Collapse
Affiliation(s)
- Laura Lupi
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | | | | |
Collapse
|
50
|
Krott LB, Barbosa MC. Anomalies in a waterlike model confined between plates. J Chem Phys 2013; 138:084505. [DOI: 10.1063/1.4792639] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|