1
|
Chang Y, Xiao S, Yu H, Ma R, Skallerud BH, Zhang Z, He J. Unraveling Ice-Solid Interface Rupture Dynamics: Insights from Molecular Dynamics Simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:17090-17097. [PMID: 39099469 PMCID: PMC11325647 DOI: 10.1021/acs.langmuir.4c02079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/06/2024]
Abstract
Unwanted icing on exposed surfaces poses significant risks, driving the quest for effective anti-icing mechanisms. While fracture mechanics concepts have been developed for designing coatings that weaken the ice-solid interface on soft surfaces, the factors that dictate ice adhesion strength and its counterpart, ice removal force, on hard surfaces remain poorly understood. In this study, we employ molecular dynamics simulations to investigate the interface rupture between ice and a hard solid substrate. The results indicate that the ice adhesion strength is contingent on the length of the ice cube. By examining the shearing behavior, we reveal a nanoscale critical force-bearing length. The shear force required to detach the ice scales proportionally with the length of the ice cube when it is smaller than the critical length. Once the ice cube length exceeds the critical length, the shear force stabilizes at a constant maximum value, revealing the existence of a maximum ice-removal force. The results align with the so-called strength versus toughness-controlled deicing regimes and are in agreement with cohesive zone modeling at the continuum length scale and recent experimental results. Our results extend this understanding to the nanoscale, confirming consistency between macro and micro scales. This consistency suggests that the toughness of the ice-solid interface is intrinsically governed by ice-surface interactions. By unraveling key intrinsic factors and their scale-dependent effects on the interface rupture of ice on surfaces, this study lays a solid theoretical foundation for the design and fabrication of next-generation anti-icing surfaces.
Collapse
Affiliation(s)
- Yuanhao Chang
- NTNU Nanomechanical Lab, Department of Structural Engineering, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | - Senbo Xiao
- NTNU Nanomechanical Lab, Department of Structural Engineering, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | - Haiyang Yu
- Division of Applied Mechanics, Department of Materials Science and Engineering, Uppsala University, SE-75121 Uppsala, Sweden
| | - Rui Ma
- NTNU Nanomechanical Lab, Department of Structural Engineering, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | - Bjørn Helge Skallerud
- NTNU Nanomechanical Lab, Department of Structural Engineering, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | - Zhiliang Zhang
- NTNU Nanomechanical Lab, Department of Structural Engineering, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | - Jianying He
- NTNU Nanomechanical Lab, Department of Structural Engineering, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| |
Collapse
|
2
|
Camarillo M, Oller-Iscar J, M Conde M, Ramírez J, Sanz E. Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation. J Chem Phys 2024; 160:134505. [PMID: 38557847 DOI: 10.1063/5.0188929] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2023] [Accepted: 03/14/2024] [Indexed: 04/04/2024] Open
Abstract
Heterogeneous nucleation is the main path to ice formation on Earth. The ice nucleating ability of a certain substrate is mainly determined by both molecular interactions and the structural mismatch between the ice and the substrate lattices. We focus on the latter factor using molecular simulations of the mW model. Quantifying the effect of structural mismatch alone is challenging due to its coupling with molecular interactions. To disentangle both the factors, we use a substrate composed of water molecules in such a way that any variation on the nucleation temperature can be exclusively ascribed to the structural mismatch. We find that a 1% increase in structural mismatch leads to a decrease of ∼4 K in the nucleation temperature. We also analyze the effect of orientation of the substrate with respect to the liquid. The three main ice orientations (basal, primary prism, and secondary prism) have a similar ice nucleating ability. We finally assess the effect of lattice flexibility by comparing substrates where molecules are immobile to others where a certain freedom to fluctuate around the lattice positions is allowed. Interestingly, we find that the latter type of substrate is more efficient in nucleating ice because it can adapt its structure to that of ice.
Collapse
Affiliation(s)
- M Camarillo
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J Oller-Iscar
- Department of Chemical Engineering, Universidad Politécnica de Madrid, José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - M M Conde
- Department of Chemical Engineering, Universidad Politécnica de Madrid, José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - J Ramírez
- Department of Chemical Engineering, Universidad Politécnica de Madrid, José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - E Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
3
|
Marks SM, Vicars Z, Thosar AU, Patel AJ. Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling. J Phys Chem B 2023. [PMID: 37378637 DOI: 10.1021/acs.jpcb.3c01627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.
Collapse
Affiliation(s)
- Sean M Marks
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Zachariah Vicars
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Aniket U Thosar
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amish J Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
4
|
Chew PY, Reinhardt A. Phase diagrams-Why they matter and how to predict them. J Chem Phys 2023; 158:030902. [PMID: 36681642 DOI: 10.1063/5.0131028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Understanding the thermodynamic stability and metastability of materials can help us to, for example, gauge whether crystalline polymorphs in pharmaceutical formulations are likely to be durable. It can also help us to design experimental routes to novel phases with potentially interesting properties. In this Perspective, we provide an overview of how thermodynamic phase behavior can be quantified both in computer simulations and machine-learning approaches to determine phase diagrams, as well as combinations of the two. We review the basic workflow of free-energy computations for condensed phases, including some practical implementation advice, ranging from the Frenkel-Ladd approach to thermodynamic integration and to direct-coexistence simulations. We illustrate the applications of such methods on a range of systems from materials chemistry to biological phase separation. Finally, we outline some challenges, questions, and practical applications of phase-diagram determination which we believe are likely to be possible to address in the near future using such state-of-the-art free-energy calculations, which may provide fundamental insight into separation processes using multicomponent solvents.
Collapse
Affiliation(s)
- Pin Yu Chew
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Aleks Reinhardt
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
5
|
Sosso GC, Sudera P, Backes AT, Whale TF, Fröhlich-Nowoisky J, Bonn M, Michaelides A, Backus EHG. The role of structural order in heterogeneous ice nucleation. Chem Sci 2022; 13:5014-5026. [PMID: 35655890 PMCID: PMC9067566 DOI: 10.1039/d1sc06338c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Accepted: 04/07/2022] [Indexed: 01/10/2023] Open
Abstract
The freezing of water into ice is a key process that is still not fully understood. It generally requires an impurity of some description to initiate the heterogeneous nucleation of the ice crystals. The molecular structure, as well as the extent of structural order within the impurity in question, both play an essential role in determining its effectiveness. However, disentangling these two contributions is a challenge for both experiments and simulations. In this work, we have systematically investigated the ice-nucleating ability of the very same compound, cholesterol, from the crystalline (and thus ordered) form to disordered self-assembled monolayers. Leveraging a combination of experiments and simulations, we identify a “sweet spot” in terms of the surface coverage of the monolayers, whereby cholesterol maximises its ability to nucleate ice (which remains inferior to that of crystalline cholesterol) by enhancing the structural order of the interfacial water molecules. These findings have practical implications for the rational design of synthetic ice-nucleating agents. The freezing of water into ice is still not fully understood. Here, we investigate the role of structural disorder within the biologically relevant impurities that facilitate this fundamental phase transition.![]()
Collapse
Affiliation(s)
- Gabriele C Sosso
- Department of Chemistry, University of Warwick Gibbet Hill Road Coventry CV4 7AL UK
| | - Prerna Sudera
- Max Planck Institute for Polymer Research Ackermannweg 10 55128 Mainz Germany
| | - Anna T Backes
- Max Planck Institute for Chemistry Hahn-Meitner-Weg 1 55128 Mainz Germany
| | - Thomas F Whale
- Department of Chemistry, University of Warwick Gibbet Hill Road Coventry CV4 7AL UK
| | | | - Mischa Bonn
- Max Planck Institute for Polymer Research Ackermannweg 10 55128 Mainz Germany
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge Lensfield Road Cambridge CB2 1EW UK
| | - Ellen H G Backus
- Max Planck Institute for Polymer Research Ackermannweg 10 55128 Mainz Germany.,Department of Physical Chemistry, University of Vienna Währingerstrasse 42 1090 Wien Austria
| |
Collapse
|
6
|
Mukhtyar AJ, Escobedo FA. Computing free energy barriers for the nucleation of complex network mesophases. J Chem Phys 2022; 156:034502. [DOI: 10.1063/5.0079396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ankita J. Mukhtyar
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14850, USA
| | - Fernando A. Escobedo
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14850, USA
| |
Collapse
|
7
|
Kargozarfard Z, Haghtalab A, Ayatollahi S, Badizad MH. Molecular Dynamics Simulation of Calcium Sulfate Nucleation in Homogeneous and Heterogeneous Crystallization Conditions: An Application in Water Flooding. Ind Eng Chem Res 2020. [DOI: 10.1021/acs.iecr.0c04290] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Zahra Kargozarfard
- Department of Chemical Engineering, Tarbiat Modares University, PO Box 14115-143, Tehran 1411713116, Iran
| | - Ali Haghtalab
- Department of Chemical Engineering, Tarbiat Modares University, PO Box 14115-143, Tehran 1411713116, Iran
| | - Shahab Ayatollahi
- Sharif Upstream Petroleum Research Institute, Department of Chemical and Petroleum Engineering, Sharif University of Technology, Tehran 11365-9465, Iran
| | - Mohammad Hasan Badizad
- Sharif Upstream Petroleum Research Institute, Department of Chemical and Petroleum Engineering, Sharif University of Technology, Tehran 11365-9465, Iran
| |
Collapse
|
8
|
Lata NN, Zhou J, Hamilton P, Larsen M, Sarupria S, Cantrell W. Multivalent Surface Cations Enhance Heterogeneous Freezing of Water on Muscovite Mica. J Phys Chem Lett 2020; 11:8682-8689. [PMID: 32955892 DOI: 10.1021/acs.jpclett.0c02121] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Heterogeneous ice nucleation is a crucial phenomenon in various fields of fundamental and applied science. We investigate the effect of surface cations on freezing of water on muscovite mica. Mica is unique in that the exposed ion on its surface can be readily and easily exchanged without affecting other properties such as surface roughness. We investigate freezing on natural (K+) mica and mica in which we have exchanged K+ for Al3+, Mg2+, Ca2+, and Sr2+. We find that liquid water freezes at higher temperatures when ions of higher valency are present on the surface, thus exposing more of the underlying silica layer. Our data also show that the size of the ion affects the characteristic freezing temperature. Using molecular dynamics simulations, we investigate the effects that the ion valency and exposed silica layer have on the behavior of water on the surface. The results indicate that multivalent cations enhance the probability of forming large clusters of hydrogen bonded water molecules that are anchored by the hydration shells of the cations. These clusters also have a large fraction of free water that can reorient to take ice-like configurations, which are promoted by the regions on mica devoid of the ions. Thus, these clusters could serve as seedbeds for ice nuclei. The combined experimental and simulation studies shed new light on the influence of surface ions on heterogeneous ice nucleation.
Collapse
Affiliation(s)
- Nurun Nahar Lata
- Atmospheric Sciences Program, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Jiarun Zhou
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States
| | - Pearce Hamilton
- Department of Physics and Astronomy, College of Charleston, Charleston, South Carolina 29424, United States
| | - Michael Larsen
- Department of Physics and Astronomy, College of Charleston, Charleston, South Carolina 29424, United States
- Atmospheric Sciences Program and Department of Physics, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Sapna Sarupria
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States
| | - Will Cantrell
- Atmospheric Sciences Program and Department of Physics, Michigan Technological University, Houghton, Michigan 49931, United States
| |
Collapse
|
9
|
Shao M, Zhang C, Qi C, Wang C, Wang J, Ye F, Zhou X. Hydrogen polarity of interfacial water regulates heterogeneous ice nucleation. Phys Chem Chem Phys 2019; 22:258-264. [PMID: 31808477 DOI: 10.1039/c9cp04867g] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Using all-atomic molecular dynamics (MD) simulations, we show that the structure of interfacial water (IW) induced by substrates characterizes the ability of a substrate to nucleate ice. We probe the shape and structure of ice nuclei and the corresponding supercooling temperatures to measure the ability of IW with various hydrogen polarities for ice nucleation, and find that the hydrogen polarization of IW even with the ice-like oxygen lattice increases the contact angle of the ice nucleus on IW, thus lifting the free energy barrier of heterogeneous ice nucleation. The results show that not only the oxygen lattice order but the hydrogen disorder of IW on substrates are required to effectively facilitate the freezing of top water.
Collapse
Affiliation(s)
- Mingzhe Shao
- College of Light Industry Science and Engineering, Tianjin University of Science and Technology, Tianjin, China
| | | | | | | | | | | | | |
Collapse
|
10
|
Zhao W, Huang H, Bi Q, Xu Y, Lü Y. One-dimensional water nanowires induced by electric fields. Phys Chem Chem Phys 2019; 21:19414-19422. [PMID: 31460524 DOI: 10.1039/c9cp02788b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Self-aggregation of water vapour molecules under external electric fields is systemically investigated by using molecular dynamics simulations. It is found that small water clusters aggregate into one-dimensional water nanowires along the electric field direction. The electric field strength plays a crucial role in tuning the nanowire structure. Under relatively weak electric fields such as E = 0.1 V Å-1, square and pentagonal prism-like structures are preferred; when intermediate strength electric fields are applied (E = 1.0 V Å-1), water nanowires featuring a disordered mixture of four-, five- and six-membered rings are formed; and an open ordered structure which is reminiscent of two-dimensional (2D) ice is observed when the field strength becomes very high (E > 3.0 V Å-1). Bond parameter analysis based on density-functional theory calculations shows that the electric field affects anisotropically the conformation of water molecules as well as the hydrogen-bond properties. Along the electric field, the H-O bond is stretched and the hydrogen bond shrinks with field strength in contrast to the changes perpendicular to the electric field. As a result, the hydrogen bonding is enhanced along the electric field. Under very high electric fields, the anisotropic hydrogen-bond network opens up via breaking of the bonds perpendicular to the electric field and ultimately relaxes into a loose quasi-2D ordered network.
Collapse
Affiliation(s)
- Wan Zhao
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Haishen Huang
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Qingling Bi
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Yujia Xu
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Yongjun Lü
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| |
Collapse
|
11
|
Leoni F, Shi R, Tanaka H, Russo J. Crystalline clusters in mW water: Stability, growth, and grain boundaries. J Chem Phys 2019; 151:044505. [DOI: 10.1063/1.5100812] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Fabio Leoni
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | - Rui Shi
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - John Russo
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
12
|
Shi R, Tanaka H. Homogeneous nucleation of ferroelectric ice crystal driven by spontaneous dipolar ordering in supercooled TIP5P water. J Chem Phys 2019; 151:024501. [DOI: 10.1063/1.5100634] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Rui Shi
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
13
|
Sosso GC, Whale TF, Holden MA, Pedevilla P, Murray BJ, Michaelides A. Unravelling the origins of ice nucleation on organic crystals. Chem Sci 2018; 9:8077-8088. [PMID: 30542556 PMCID: PMC6238755 DOI: 10.1039/c8sc02753f] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2018] [Accepted: 08/27/2018] [Indexed: 12/01/2022] Open
Abstract
Organic molecules such as steroids or amino acids form crystals that can facilitate the formation of ice - arguably the most important phase transition on earth. However, the origin of the ice nucleating ability of organic crystals is still largely unknown. Here, we combine experiments and simulations to unravel the microscopic details of ice formation on cholesterol, a prototypical organic crystal widely used in cryopreservation. We find that cholesterol - which is also a substantial component of cell membranes - is an ice nucleating agent more potent than many inorganic substrates, including the mineral feldspar (one of the most active ice nucleating materials in the atmosphere). Scanning electron microscopy measurements reveal a variety of morphological features on the surfaces of cholesterol crystals: this suggests that the topography of the surface is key to the broad range of ice nucleating activity observed (from -4 to -20 °C). In addition, we show via molecular simulations that cholesterol crystals aid the formation of ice nuclei in a unconventional fashion. Rather than providing a template for a flat ice-like contact layer (as found in the case of many inorganic substrates), the flexibility of the cholesterol surface and its low density of hydrophilic functional groups leads to the formation of molecular cages involving both water molecules and terminal hydroxyl groups of the cholesterol surface. These cages are made of 6- and, surprisingly, 5-membered hydrogen bonded rings of water and hydroxyl groups that favour the nucleation of hexagonal as well as cubic ice (a rare occurrence). We argue that the phenomenal ice nucleating activity of steroids such as cholesterol (and potentially of many other organic crystals) is due to (i) the ability of flexible hydrophilic surfaces to form unconventional ice-templating structures and (ii) the different nucleation sites offered by the diverse topography of the crystalline surfaces. These findings clarify how exactly organic crystals promote the formation of ice, thus paving the way toward deeper understanding of ice formation in soft and biological matter - with obvious reverberations on atmospheric science and cryobiology.
Collapse
Affiliation(s)
- Gabriele C Sosso
- Department of Chemistry and Centre for Scientific Computing , University of Warwick , Gibbet Hill Road , Coventry CV4 7AL , UK .
| | - Thomas F Whale
- School of Earth and Environment , University of Leeds , Leeds LS2 9JT , UK
| | - Mark A Holden
- School of Earth and Environment , University of Leeds , Leeds LS2 9JT , UK
- Chemistry , University of Leeds , Leeds LS2 9JT , UK
| | - Philipp Pedevilla
- Thomas Young Centre , London Centre for Nanotechnology and Department of Physics and Astronomy , University College London , London WC1E 6BT , UK
| | - Benjamin J Murray
- School of Earth and Environment , University of Leeds , Leeds LS2 9JT , UK
| | - Angelos Michaelides
- Thomas Young Centre , London Centre for Nanotechnology and Department of Physics and Astronomy , University College London , London WC1E 6BT , UK
| |
Collapse
|
14
|
Pedevilla P, Fitzner M, Sosso GC, Michaelides A. Heterogeneous seeded molecular dynamics as a tool to probe the ice nucleating ability of crystalline surfaces. J Chem Phys 2018; 149:072327. [PMID: 30134662 DOI: 10.1063/1.5029336] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ice nucleation plays a significant role in a large number of natural and technological processes, but it is challenging to investigate experimentally because of the small time scales (ns) and short length scales (nm) involved. On the other hand, conventional molecular simulations struggle to cope with the relatively long time scale required for critical ice nuclei to form. One way to tackle this issue is to take advantage of free energy or path sampling techniques. Unfortunately, these are computationally costly. Seeded molecular dynamics is a much less demanding alternative that has been successfully applied already to study the homogeneous freezing of water. However, in the case of heterogeneous ice nucleation, nature's favourite route to form ice, an array of suitable interfaces between the ice seeds and the substrate of interest has to be built, and this is no trivial task. In this paper, we present a Heterogeneous SEEDing (HSEED) approach which harnesses a random structure search framework to tackle the ice-substrate challenge, thus enabling seeded molecular dynamics simulations of heterogeneous ice nucleation on crystalline surfaces. We validate the HSEED framework by investigating the nucleation of ice on (i) model crystalline surfaces, using the coarse-grained mW model, and (ii) cholesterol crystals, employing the fully atomistic TIP4P/ice water model. We show that the HSEED technique yields results in excellent agreement with both metadynamics and forward flux sampling simulations. Because of its computational efficiency, the HSEED method allows one to rapidly assess the ice nucleation ability of whole libraries of crystalline substrates-a long-awaited computational development in, e.g., atmospheric science.
Collapse
Affiliation(s)
- Philipp Pedevilla
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Martin Fitzner
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Gabriele C Sosso
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
15
|
Shevkunov S, Singh JK. Bicanonical ensemble Monte Carlo simulation of water condensation in the field of crystal lattice defects. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.05.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
16
|
Guo J, Haji-Akbari A, Palmer JC. Hybrid Monte Carlo with LAMMPS. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2018. [DOI: 10.1142/s0219633618400023] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We describe a strategy for performing canonical and isothermal-isobaric ensemble hybrid Monte Carlo (HMC) simulations with the widely-used Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) molecular dynamics (MD) software package. The overall workflow for the HMC simulations is handled using an external Python driver script, which invokes LAMMPS’ library interface to perform numerically intensive tasks such as MD integration. We document several rigorous consistency checks that have been used to validate our HMC implementation. We also demonstrate that our approach can be readily extended to implement biased HMC sampling schemes for computing free energies. Codes and input files from the documented examples are available on the web.
Collapse
Affiliation(s)
- Jingxiang Guo
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, CT 06520, USA
| | - Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, USA
| |
Collapse
|
17
|
Shevkunov SV. The Effect of Temperature on Nucleation of Condensed Water Phase on the Surface of a β-AgI Crystal. 1. Structure. COLLOID JOURNAL 2018. [DOI: 10.1134/s1061933x18020096] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
18
|
Qiu Y, Lupi L, Molinero V. Is Water at the Graphite Interface Vapor-like or Ice-like? J Phys Chem B 2018; 122:3626-3634. [PMID: 29298058 DOI: 10.1021/acs.jpcb.7b11476] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Graphitic surfaces are the main component of soot, a major constituent of atmospheric aerosols. Experiments indicate that soots of different origins display a wide range of abilities to heterogeneously nucleate ice. The ability of pure graphite to nucleate ice in experiments, however, seems to be almost negligible. Nevertheless, molecular simulations with the monatomic water model mW with water-carbon interactions parameterized to reproduce the experimental contact angle of water on graphite predict that pure graphite nucleates ice. According to classical nucleation theory, the ability of a surface to nucleate ice is controlled by the binding free energy between ice immersed in liquid water and the surface. To establish whether the discrepancy in freezing efficiencies of graphite in mW simulations and experiments arises from the coarse resolution of the model or can be fixed by reparameterization, it is important to elucidate the contributions of the water-graphite, water-ice, and ice-water interfaces to the free energy, enthalpy, and entropy of binding for both water and the model. Here we use thermodynamic analysis and free energy calculations to determine these interfacial properties. We demonstrate that liquid water at the graphite interface is not ice-like or vapor-like: it has similar free energy, entropy, and enthalpy as water in the bulk. The thermodynamics of the water-graphite interface is well reproduced by the mW model. We find that the entropy of binding between graphite and ice is positive and dominated, in both experiments and simulations, by the favorable entropy of reducing the ice-water interface. Our analysis indicates that the discrepancy in freezing efficiencies of graphite in experiments and the simulations with mW arises from the inability of the model to simultaneously reproduce the contact angle of liquid water on graphite and the free energy of the ice-graphite interface. This transferability issue is intrinsic to the resolution of the model, and arises from its lack of rotational degrees of freedom.
Collapse
Affiliation(s)
- Yuqing Qiu
- Department of Chemistry , The University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112-0850 , United States
| | - Laura Lupi
- 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
|
Glatz B, Sarupria S. Heterogeneous Ice Nucleation: Interplay of Surface Properties and Their Impact on Water Orientations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:1190-1198. [PMID: 29020452 DOI: 10.1021/acs.langmuir.7b02859] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Ice is ubiquitous in nature, and heterogeneous ice nucleation is the most common pathway of ice formation. How surface properties affect the propensity to observe ice nucleation on that surface remains an open question. We present results of molecular dynamics studies of heterogeneous ice nucleation on model surfaces. The models surfaces considered emulate the chemistry of kaolinite, an abundant component of mineral dust. We investigate the interplay of surface lattice and hydrogen bonding properties in affecting ice nucleation. We find that lattice matching and hydrogen bonding are necessary but not sufficient conditions for observing ice nucleation at these surfaces. We correlate this behavior to the orientations sampled by the metastable supercooled water in contact with the surfaces. We find that ice is observed in cases where water molecules not only sample orientations favorable for bilayer formation but also do not sample unfavorable orientations. This distribution depends on both surface-water and water-water interactions and can change with subtle modifications to the surface properties. Our results provide insights into the diverse behavior of ice nucleation observed at different surfaces and highlight the complexity in elucidating heterogeneous ice nucleation.
Collapse
Affiliation(s)
- Brittany Glatz
- Department of Chemical & Biomolecular Engineering, Clemson University , Clemson, South Carolina 29634, United States
| | - Sapna Sarupria
- Department of Chemical & Biomolecular Engineering, Clemson University , Clemson, South Carolina 29634, United States
| |
Collapse
|
20
|
Wang X, Gao Y, Xu Y, Li X, Jiang J, Hou J, Han W, Li Q, Shen C. A Prerequisite of the Poly(ε-Caprolactone) Self-Induced Nanohybrid Shish-Kebab Structure Formation: An Ordered Crystal Lamellae Orientation Morphology of Fibers. MACROMOL CHEM PHYS 2017. [DOI: 10.1002/macp.201700414] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Xiaofeng Wang
- National Center for International Research of Micro-Nano Molding Technology; Henan Key Laborotary of Micro Molding Technology; Zhengzhou University; Zhengzhou 450001 China
- School of Mechanics & Engineering Science; Zhengzhou University; Zhengzhou 450001 China
- State Key Laboratory of Molecular Engineering of Polymers (Fudan University); Shanghai 200433 China
| | - Yanhong Gao
- School of Materials Science & Engineering; Zhengzhou University; Zhengzhou 450001 China
| | - Yiyang Xu
- National Center for International Research of Micro-Nano Molding Technology; Henan Key Laborotary of Micro Molding Technology; Zhengzhou University; Zhengzhou 450001 China
- School of Mechanics & Engineering Science; Zhengzhou University; Zhengzhou 450001 China
| | - Xuyan Li
- National Center for International Research of Micro-Nano Molding Technology; Henan Key Laborotary of Micro Molding Technology; Zhengzhou University; Zhengzhou 450001 China
- School of Mechanics & Engineering Science; Zhengzhou University; Zhengzhou 450001 China
| | - Jing Jiang
- National Center for International Research of Micro-Nano Molding Technology; Henan Key Laborotary of Micro Molding Technology; Zhengzhou University; Zhengzhou 450001 China
- School of Mechanics & Engineering Science; Zhengzhou University; Zhengzhou 450001 China
| | - Jianhua Hou
- National Center for International Research of Micro-Nano Molding Technology; Henan Key Laborotary of Micro Molding Technology; Zhengzhou University; Zhengzhou 450001 China
- School of Mechanics & Engineering Science; Zhengzhou University; Zhengzhou 450001 China
| | - WenJuan Han
- School of Materials Science & Engineering; Zhengzhou University; Zhengzhou 450001 China
| | - Qian Li
- National Center for International Research of Micro-Nano Molding Technology; Henan Key Laborotary of Micro Molding Technology; Zhengzhou University; Zhengzhou 450001 China
- School of Mechanics & Engineering Science; Zhengzhou University; Zhengzhou 450001 China
| | - Changyu Shen
- National Center for International Research of Micro-Nano Molding Technology; Henan Key Laborotary of Micro Molding Technology; Zhengzhou University; Zhengzhou 450001 China
- School of Mechanics & Engineering Science; Zhengzhou University; Zhengzhou 450001 China
| |
Collapse
|
21
|
Fitzner M, Joly L, Ma M, Sosso GC, Zen A, Michaelides A. Communication: Truncated non-bonded potentials can yield unphysical behavior in molecular dynamics simulations of interfaces. J Chem Phys 2017; 147:121102. [DOI: 10.1063/1.4997698] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Martin Fitzner
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Laurent Joly
- Univ. Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne,
France
| | - Ming Ma
- Department of Mechanical Engineering, State Key Laboratory of Tribology and Center for Nano and Micro Mechanics, Tsinghua University, Beijing 100084, China
| | - Gabriele C. Sosso
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL,
United Kingdom
| | - Andrea Zen
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
22
|
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]
|
23
|
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]
|
24
|
Zhang Y, Liao Q, Wang X, Yao J, Fu H. Lattice-Matched Epitaxial Growth of Organic Heterostructures for Integrated Optoelectronic Application. Angew Chem Int Ed Engl 2017; 56:3616-3620. [PMID: 28233948 DOI: 10.1002/anie.201700447] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2017] [Indexed: 11/08/2022]
Abstract
Development of nanowire photonics requires integration of different nanowire components into highly ordered functional heterostructures. Herein, we report a sequential self-assembly of binary molecular components into branched nanowire heterostructures (BNHs) via lattice-matched epitaxial growth, in which the microribbon backbone of 2,5-Bis(5-tert-butyl-2-benzoxazolyl)thiophene (BBOT) functions as blue-emitting microlaser source to pump the nanowire branches of BODIPY. By constructing Au electrodes on both branch sides and measuring the photocurrent in them, we successfully realize the integration of an organic laser and a power meter in a single device. This work provides a new insight into the integration of 1D organic nanostructures into BNHs for realizing organic multifunctional photonic devices.
Collapse
Affiliation(s)
- Yi Zhang
- Beijing National Laboratory for Molecular Sciences (BNLMS), Institute of chemistry, Chinese Academy of Sciences, Beijing, 100190, P.R. China.,University of Chinese Academy of Sciences, Beijing, 100049, P.R. China
| | - Qing Liao
- Beijing Key Laboratory for Optical Materials and Photonic Devices, Department of Chemistry, Capital Normal University, Beijing, 100048, P.R. China
| | - Xinguo Wang
- Beijing Key Laboratory for Optical Materials and Photonic Devices, Department of Chemistry, Capital Normal University, Beijing, 100048, P.R. China
| | - Jiannian Yao
- Beijing National Laboratory for Molecular Sciences (BNLMS), Institute of chemistry, Chinese Academy of Sciences, Beijing, 100190, P.R. China
| | - Hongbing Fu
- Beijing National Laboratory for Molecular Sciences (BNLMS), Institute of chemistry, Chinese Academy of Sciences, Beijing, 100190, P.R. China.,Beijing Key Laboratory for Optical Materials and Photonic Devices, Department of Chemistry, Capital Normal University, Beijing, 100048, P.R. China.,Tianjin Key Laboratory of Molecular Optoelectronic Sciences, Department of Chemistry, Tianjin University, Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin, 300072, P.R. China
| |
Collapse
|
25
|
Zhang Y, Liao Q, Wang X, Yao J, Fu H. Lattice‐Matched Epitaxial Growth of Organic Heterostructures for Integrated Optoelectronic Application. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201700447] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Yi Zhang
- Beijing National Laboratory for Molecular Sciences (BNLMS) Institute of chemistry Chinese Academy of Sciences Beijing 100190 P.R. China
- University of Chinese Academy of Sciences Beijing 100049 P.R. China
| | - Qing Liao
- Beijing Key Laboratory for Optical Materials and Photonic Devices Department of Chemistry Capital Normal University Beijing 100048 P.R. China
| | - Xinguo Wang
- Beijing Key Laboratory for Optical Materials and Photonic Devices Department of Chemistry Capital Normal University Beijing 100048 P.R. China
| | - Jiannian Yao
- Beijing National Laboratory for Molecular Sciences (BNLMS) Institute of chemistry Chinese Academy of Sciences Beijing 100190 P.R. China
| | - Hongbing Fu
- Beijing National Laboratory for Molecular Sciences (BNLMS) Institute of chemistry Chinese Academy of Sciences Beijing 100190 P.R. China
- Beijing Key Laboratory for Optical Materials and Photonic Devices Department of Chemistry Capital Normal University Beijing 100048 P.R. China
- Tianjin Key Laboratory of Molecular Optoelectronic Sciences Department of Chemistry Tianjin University Collaborative Innovation Center of Chemical Science and Engineering (Tianjin) Tianjin 300072 P.R. China
| |
Collapse
|
26
|
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
|
27
|
Bourque AJ, Locker CR, Rutledge GC. Heterogeneous Nucleation of an n-Alkane on Tetrahedrally Coordinated Crystals. J Phys Chem B 2017; 121:904-911. [PMID: 28071905 DOI: 10.1021/acs.jpcb.6b12590] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Heterogeneous nucleation refers to the propensity for phase transformations to initiate preferentially on foreign surfaces, such as vessel walls, dust particles, or formulation additives. In crystallization, the form of the initial nucleus has ramifications for the crystallographic form, morphology, and properties of the resulting solid. Nevertheless, the discovery and design of nucleating agents remains a matter of trial and error because of the very small spatiotemporal scales over which the critical nucleus is formed and the extreme difficulty of examining such events empirically. Using molecular dynamics simulations, we demonstrate a method for the rapid screening of entire families of materials for activity as nucleating agents and for characterizing their mechanism of action. The method is applied to the crystallization of n-pentacontane, a model surrogate for polyethylene, on the family of tetrahedrally coordinated crystals, including diamond and silicon. A systematic variation of parameters in the interaction potential permits a comprehensive, physically based screening of nucleating agents in this class of materials, including both real and hypothetical candidates. The induction time for heterogeneous nucleation is shown to depend strongly on crystallographic registry between the nucleating agent and the critical nucleus, indicative of an epitaxial mechanism in this class of materials. Importantly, the severity of this registry requirement weakens with decreasing rigidity of the substrate and increasing strength of attraction to the surface of the nucleating agent. Employing this method, a high-throughput computational screening of nucleating agents becomes possible, facilitating the discovery of novel nucleating agents within a broad "materials genome" of possible additives.
Collapse
Affiliation(s)
- Alexander J Bourque
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| | - C Rebecca Locker
- ExxonMobil Research and Engineering Company , Annandale, New Jersey 08801, United States
| | - Gregory C Rutledge
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| |
Collapse
|
28
|
Desgranges C, Delhommelle J. Free energy calculations along entropic pathways. II. Droplet nucleation in binary mixtures. J Chem Phys 2016; 145:234505. [DOI: 10.1063/1.4972011] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
29
|
He X, Shen Y, Hung FR, Santiso EE. Heterogeneous nucleation from a supercooled ionic liquid on a carbon surface. J Chem Phys 2016; 145:211919. [DOI: 10.1063/1.4963336] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Xiaoxia He
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Yan Shen
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Francisco R. Hung
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
- Center for Computation & Technology, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Erik E. Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
30
|
Glatz B, Sarupria S. The surface charge distribution affects the ice nucleating efficiency of silver iodide. J Chem Phys 2016; 145:211924. [DOI: 10.1063/1.4966018] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Brittany Glatz
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, USA
| | - Sapna Sarupria
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, USA
| |
Collapse
|
31
|
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
|
32
|
Desgranges C, Delhommelle J. Free energy calculations along entropic pathways. I. Homogeneous vapor-liquid nucleation for atomic and molecular systems. J Chem Phys 2016; 145:204112. [DOI: 10.1063/1.4968231] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
33
|
Shou W, Pan H. Silicon-wall interfacial free energy via thermodynamics integration. J Chem Phys 2016; 145:184702. [PMID: 27846694 DOI: 10.1063/1.4966975] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
We compute the interfacial free energy of a silicon system in contact with flat and structured walls by molecular dynamics simulation. The thermodynamics integration method, previously applied to Lennard-Jones potentials [R. Benjamin and J. Horbach, J. Chem. Phys. 137, 044707 (2012)], has been extended and implemented in Tersoff potentials with two-body and three-body interactions taken into consideration. The thermodynamic integration scheme includes two steps. In the first step, the bulk Tersoff system is reversibly transformed to a state where it interacts with a structureless flat wall, and in a second step, the flat structureless wall is reversibly transformed into an atomistic SiO2 wall. Interfacial energies for liquid silicon-wall interfaces and crystal silicon-wall interfaces have been calculated. The calculated interfacial energies have been employed to predict the nucleation mechanisms in a slab of liquid silicon confined by two walls and compared with MD simulation results.
Collapse
Affiliation(s)
- Wan Shou
- Department of Mechanical and Aerospace Engineering, Missouri University of Science and Technology, Rolla, Missouri 65401, USA
| | - Heng Pan
- Department of Mechanical and Aerospace Engineering, Missouri University of Science and Technology, Rolla, Missouri 65401, USA
| |
Collapse
|
34
|
Sosso G, Chen J, Cox SJ, Fitzner M, Pedevilla P, Zen A, Michaelides A. Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations. Chem Rev 2016; 116:7078-116. [PMID: 27228560 PMCID: PMC4919765 DOI: 10.1021/acs.chemrev.5b00744] [Citation(s) in RCA: 379] [Impact Index Per Article: 47.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2015] [Indexed: 11/28/2022]
Abstract
The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.
Collapse
Affiliation(s)
- Gabriele
C. Sosso
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Ji Chen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | | | - Martin Fitzner
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Philipp Pedevilla
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Andrea Zen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Angelos Michaelides
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| |
Collapse
|
35
|
Pedevilla P, Cox SJ, Slater B, Michaelides A. Can Ice-Like Structures Form on Non-Ice-Like Substrates? The Example of the K-feldspar Microcline. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2016; 120:6704-6713. [PMID: 27917255 PMCID: PMC5127609 DOI: 10.1021/acs.jpcc.6b01155] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 03/02/2016] [Indexed: 05/12/2023]
Abstract
Feldspar minerals are the most common rock formers in Earth's crust. As such they play an important role in subjects ranging from geology to climate science. An atomistic understanding of the feldspar structure and its interaction with water is therefore desirable, not least because feldspar has been shown to dominate ice nucleation by mineral dusts in Earth's atmosphere. The complexity of the ice/feldspar interface arising from the numerous chemical motifs expressed on the surface makes it a challenging system. Here we report a comprehensive study of this challenging system with ab initio density functional theory calculations. We show that the distribution of Al atoms, which is crucial for the dissolution kinetics of tectosilicate minerals, differs significantly between the bulk environment and on the surface. Furthermore, we demonstrate that water does not form ice-like overlayers in the contact layer on the most easily cleaved (001) surface of K-feldspar. We do, however, identify contact layer structures of water that induce ice-like ordering in the second overlayer. This suggests that even substrates without an apparent match with the ice structure may still act as excellent ice nucleating agents.
Collapse
Affiliation(s)
- Philipp Pedevilla
- Thomas
Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- London
Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Stephen J. Cox
- Thomas
Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- London
Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Ben Slater
- Thomas
Young Centre and 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 and Department of Physics and Astronomy, University College London, Gower Street, London WC1E
6BT, United Kingdom
- E-mail: . Phone: +44 207 679
0647
| |
Collapse
|
36
|
Fitzner M, Sosso GC, Cox SJ, Michaelides A. The Many Faces of Heterogeneous Ice Nucleation: Interplay Between Surface Morphology and Hydrophobicity. J Am Chem Soc 2015; 137:13658-69. [PMID: 26434775 DOI: 10.1021/jacs.5b08748] [Citation(s) in RCA: 105] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
What makes a material a good ice nucleating agent? Despite the importance of heterogeneous ice nucleation to a variety of fields, from cloud science to microbiology, major gaps in our understanding of this ubiquitous process still prevent us from answering this question. In this work, we have examined the ability of generic crystalline substrates to promote ice nucleation as a function of the hydrophobicity and the morphology of the surface. Nucleation rates have been obtained by brute-force molecular dynamics simulations of coarse-grained water on top of different surfaces of a model fcc crystal, varying the water-surface interaction and the surface lattice parameter. It turns out that the lattice mismatch of the surface with respect to ice, customarily regarded as the most important requirement for a good ice nucleating agent, is at most desirable but not a requirement. On the other hand, the balance between the morphology of the surface and its hydrophobicity can significantly alter the ice nucleation rate and can also lead to the formation of up to three different faces of ice on the same substrate. We have pinpointed three circumstances where heterogeneous ice nucleation can be promoted by the crystalline surface: (i) the formation of a water overlayer that acts as an in-plane template; (ii) the emergence of a contact layer buckled in an ice-like manner; and (iii) nucleation on compact surfaces with very high interaction strength. We hope that this extensive systematic study will foster future experimental work aimed at testing the physiochemical understanding presented herein.
Collapse
Affiliation(s)
- Martin Fitzner
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Gabriele C Sosso
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Stephen J Cox
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
37
|
Fraux G, Doye JPK. Note: Heterogeneous ice nucleation on silver-iodide-like surfaces. J Chem Phys 2015; 141:216101. [PMID: 25481172 DOI: 10.1063/1.4902382] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Guillaume Fraux
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
38
|
Cox SJ, Kathmann SM, Slater B, Michaelides A. Molecular simulations of heterogeneous ice nucleation. II. Peeling back the layers. J Chem Phys 2015; 142:184705. [PMID: 25978903 DOI: 10.1063/1.4919715] [Citation(s) in RCA: 67] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Coarse grained molecular dynamics simulations are presented in which the sensitivity of the ice nucleation rate to the hydrophilicity of a graphene nanoflake is investigated. We find that an optimal interaction strength for promoting ice nucleation exists, which coincides with that found previously for a face centered cubic (111) surface. We further investigate the role that the layering of interfacial water plays in heterogeneous ice nucleation and demonstrate that the extent of layering is not a good indicator of ice nucleating ability for all surfaces. Our results suggest that to be an efficient ice nucleating agent, a surface should not bind water too strongly if it is able to accommodate high coverages of water.
Collapse
Affiliation(s)
- Stephen J Cox
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Ben Slater
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
39
|
Cox SJ, Kathmann SM, Slater B, Michaelides A. Molecular simulations of heterogeneous ice nucleation. I. Controlling ice nucleation through surface hydrophilicity. J Chem Phys 2015; 142:184704. [PMID: 25978902 DOI: 10.1063/1.4919714] [Citation(s) in RCA: 104] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Ice formation is one of the most common and important processes on earth and almost always occurs at the surface of a material. A basic understanding of how the physicochemical properties of a material's surface affect its ability to form ice has remained elusive. Here, we use molecular dynamics simulations to directly probe heterogeneous ice nucleation at a hexagonal surface of a nanoparticle of varying hydrophilicity. Surprisingly, we find that structurally identical surfaces can both inhibit and promote ice formation and analogous to a chemical catalyst, it is found that an optimal interaction between the surface and the water exists for promoting ice nucleation. We use our microscopic understanding of the mechanism to design a modified surface in silico with enhanced ice nucleating ability.
Collapse
Affiliation(s)
- Stephen J Cox
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Ben Slater
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
40
|
Whale TF, Rosillo-Lopez M, Murray BJ, Salzmann CG. Ice Nucleation Properties of Oxidized Carbon Nanomaterials. J Phys Chem Lett 2015; 6:3012-6. [PMID: 26267196 PMCID: PMC4550106 DOI: 10.1021/acs.jpclett.5b01096] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2015] [Accepted: 07/13/2015] [Indexed: 05/20/2023]
Abstract
Heterogeneous ice nucleation is an important process in many fields, particularly atmospheric science, but is still poorly understood. All known inorganic ice nucleating particles are relatively large in size and tend to be hydrophilic. Hence it is not obvious that carbon nanomaterials should nucleate ice. However, in this paper we show that four different readily water-dispersible carbon nanomaterials are capable of nucleating ice. The tested materials were carboxylated graphene nanoflakes, graphene oxide, oxidized single walled carbon nanotubes and oxidized multiwalled carbon nanotubes. The carboxylated graphene nanoflakes have a diameter of ∼30 nm and are among the smallest entities observed so far to nucleate ice. Overall, carbon nanotubes were found to nucleate ice more efficiently than flat graphene species, and less oxidized materials nucleated ice more efficiently than more oxidized species. These well-defined carbon nanomaterials may pave the way to bridging the gap between experimental and computational studies of ice nucleation.
Collapse
Affiliation(s)
- Thomas F. Whale
- Institute
for Climate and Atmospheric Science, School of Earth and Environment, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, United Kingdom
- E-mail:
| | - Martin Rosillo-Lopez
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Benjamin J. Murray
- Institute
for Climate and Atmospheric Science, School of Earth and Environment, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, United Kingdom
| | - Christoph G. Salzmann
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
41
|
|
42
|
Malkin TL, Murray BJ, Salzmann CG, Molinero V, Pickering SJ, Whale TF. Stacking disorder in ice I. Phys Chem Chem Phys 2015; 17:60-76. [PMID: 25380218 DOI: 10.1039/c4cp02893g] [Citation(s) in RCA: 141] [Impact Index Per Article: 15.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Traditionally, ice I was considered to exist in two well-defined crystalline forms at ambient pressure: stable hexagonal ice (ice Ih) and metastable cubic ice (ice Ic). However, it is becoming increasingly evident that what has been called cubic ice in the past does not have a structure consistent with the cubic crystal system. Instead, it is a stacking-disordered material containing cubic sequences interlaced with hexagonal sequences, which is termed stacking-disordered ice (ice Isd). In this article, we summarise previous work on ice with stacking disorder including ice that was called cubic ice in the past. We also present new experimental data which shows that ice which crystallises after heterogeneous nucleation in water droplets containing solid inclusions also contains stacking disorder even at freezing temperatures of around -15 °C. This supports the results from molecular simulations, that the structure of ice that crystallises initially from supercooled water is always stacking-disordered and that this metastable ice can transform to the stable hexagonal phase subject to the kinetics of recrystallization. We also show that stacking disorder in ice which forms from water droplets is quantitatively distinct from ice made via other routes. The emerging picture of ice I is that of a very complex material which frequently contains stacking disorder and this stacking disorder can vary in complexity depending on the route of formation and thermal history.
Collapse
Affiliation(s)
- Tamsin L Malkin
- Institute for Climate and Atmospheric Science, School of Earth and Environment, University of Leeds, Woodhouse Lane, Leeds, LS2 9JT, UK.
| | | | | | | | | | | |
Collapse
|
43
|
Cabriolu R, Li T. Ice nucleation on carbon surface supports the classical theory for heterogeneous nucleation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:052402. [PMID: 26066178 DOI: 10.1103/physreve.91.052402] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2014] [Indexed: 05/10/2023]
Abstract
The prevalence of heterogeneous nucleation in nature was explained qualitatively by the classical theory for heterogeneous nucleation established over more than 60 years ago, but the quantitative validity and the key conclusions of the theory have remained unconfirmed. Employing the forward flux sampling method and the coarse-grained water model (mW), we explicitly computed the heterogeneous ice nucleation rates in the supercooled water on a graphitic surface at various temperatures. The independently calculated ice nucleation rates were found to fit well according to the classical theory for heterogeneous nucleation. The fitting procedure further yields the estimate of the potency factor, which measures the ratio of the heterogeneous nucleation barrier to the homogeneous nucleation barrier. Remarkably, the estimated potency factor agrees quantitatively with the volumetric ratio of the critical nuclei between the heterogeneous and homogeneous nucleation. Our numerical study thus provides a strong support to the quantitative power of the theory and allows understanding ice nucleation behaviors under the most relevant freezing conditions.
Collapse
Affiliation(s)
- Raffaela Cabriolu
- Department of Civil and Environmental Engineering, George Washington University, Washington, DC 20052, USA
| | - Tianshu Li
- Department of Civil and Environmental Engineering, George Washington University, Washington, DC 20052, USA
| |
Collapse
|