Atomistic and coarse grained simulations reveal increased ice nucleation activity on silver iodide surfaces in slit and wedge geometries
- Institute for Atmospheric and Earth System Research/Physics, University of Helsinki, PO Box 64, FI-00014, Finland
- Institute for Atmospheric and Earth System Research/Physics, University of Helsinki, PO Box 64, FI-00014, Finland
Abstract. Ice clouds can form at low and moderate supercooling through heterogeneous ice nucleation on atmospheric particles. Typically, the nucleation requires active sites with special chemical and physical properties, including surface topology and roughness. This paper investigates microscopic mechanisms of how combinations of confinement by the surface topology and lattice match induced by the surface properties can lead to enhanced ice nucleation. We perform molecular dynamics simulations using both atomistic and coarse-grained water models, at very low supercooling, to extensively study heterogeneous ice nucleation in slit-like and concave wedge structures of silver-terminated silver iodide (0001) surfaces. We find that ice nucleation is greatly enhanced by slit-like structures when the gap width is a near-integer multiple of the thickness of an ice bilayer. For wedge systems we also do not find a simple linear dependence between ice nucleation activity and the opening angle. Instead we observe strong enhancement in concave wedge systems with angles that match the orientations of ice lattice planes, highlighting the importance of structural matching for ice nucleation in confined geometries. While in the slit systems ice cannot grow out of the slit, some wedge systems show that ice readily grows out of the wedge. In addition, some wedge systems stabilize ice structures when heating the system above the thermodynamics melting point. In the context of atmospheric ice nucleating particles, our results strongly support the experimental evidence for the importance of surface features such as cracks or pits functioning as active sites for ice nucleation at low supercooling.
Golnaz Roudsari et al.
Status: final response (author comments only)
-
RC1: 'Comment on acp-2021-1054', Anonymous Referee #1, 01 Mar 2022
General comment:
This study presents molecular dynamics simulations to investigate ice nucleation within slits and wedges of silver-terminated AgI (0001) surfaces. Ice formation depending on slit width and the opening angle of wedges was investigated. Moreover, the potential of ice to grow out of the slits and wedges was assessed. Simulations were carried out with the coarse-grained mW and the all-atom TIP4P/Ice models. It was found that slit systems promote ice nucleation when the slit width matches an integer number of ice bilayer thickness. Yet, the ice was not able to grow out of the slits. In the case of wedges, a high sensitivity to the opening angle was found with some angles inhibiting ice formation compared to the flat AgI surface and other ones enhancing it. Interestingly, the angles that enhanced or inhibited ice formation were not the same for the two water models. Moreover, ice was able to grow out of the wedges. This study exemplifies how surface geometry and templating of ice by the confining surfaces may act together to enhance ice nucleation. The differences in ice nucleation efficiency depending on the water model is well discussed and shows the potential and limitations of such simulations. Overall, it would be helpful if simulations that are described but not shown in the paper would be made accessible as supplementary information. Apart from that, this paper is well written and can be recommended for publication subject to minor revisions.
Specific comment:
Method section: information or references should be given about how cubic, hexagonal, and liquid water molecules are discriminated for the two water models.
Line 168: what is meant by “strongly hydrophilic” in terms of contact angle?
Lines 171–173: It would be interesting to see the disorder of water molecules at edges in a figure as supplementary information.
Figures 5 and 6 should correspond better with each other. Wedges with angles that are not shown in these figures should be shown as supplementary information. The meaning of "Top level" should be explained. For some wedge angles, the total number of ice molecules exceeds the number of ice molecules at the top level. This should be explained. Also the total number of water/ice molecules in each simulation should be stated.
Lines 203–205: Here, it is stated that for the W30 system ice grew to the top of the simulation cell within approximately 20 ns, but the green dashed line is already reached after less than 5 ns. This should be explained. Moreover, in the next sentence, it is stated that in the W73 system, ice grew to fill the entire cell in less than 30 ns, but the simulation only goes to 20 ns and the green dashed line is crossed already after about 12 ns. This also requires clarification.
Lines 206–207: “About seven layers of ice formed at the bottom of the W45 systems within 2 ns.” Does this statement refer to Fig. 6b? If yes, a reference to this figure could be given here.
Line 209: The systems W110 and W120 should be shown as part of SI.
Lines 212–213: Do you mean this statement in absolute terms or relative to the flat AgI surface?
Line 213: why is W45 mentioned in the bracket? There seems to be hardly any ice formation for this opening angle.
Line 226–227: the statement “none of the 15 individual simulations show ice growing beyond 3–4 ice-like layers on top of the flat surface” does not become evident based on Fig. 5g. Again, a supplementary figure would be helpful.
Line 239: what is a “compressed delay fitting parameter”? An explanation and/or a reference should be given.
Line 244: Table 1 shows an enhancement in nucleation rate of less than a factor of two for W30 and W73 compared with the flat surface. This seems to me only a minor increase in nucleation rate. What is your criterion for a “considerable” increase?
Line 273–274: again, these simulations could be shown as part of SI.
Lines 301–304: What does it mean for the predictive power of MD simulations when they are so sensitive to the specific water model? An additional comment would be helpful.
Lines 333–334: Should “Similarly to the case,” be deleted?
The content of Table 2 should be explained better.
Line 343: Which angles have been investigated at 265 K? A summarizing table with all results from both models (how many simulations produced ice and in what average time?) would be helpful.
Line 362–363: Movies as part of SI could illustrate this statement.
Line 364–365: This sentence is unclear. The formulation should be improved.
Technical comments:
Line 66: “Sects.” Instead of “Sect.”
Line 68: the point is missing at the end of the sentence.
Line 73: a comma is missing after (Abascal et al., 2005).
Figure 3: purple and blue colors are very similar and difficult to discriminate. Consider to replace e.g. purple by red.
Line 211: “times” instead of “time”.
Lines 244 and 321: “Table” should not be abbreviated by “Tab.”.
Figure caption of Fig. 7: “circles” appear to be stars.
-
RC2: 'Comment on acp-2021-1054', Anonymous Referee #2, 24 Apr 2022
This in an interesting paper that illustrate some of the complexities of heterogeneous ice nucleation on AgI by showing how the rate of nucleation in slits and wedges depends on the slit width and wedge angle. In a similar way to how AgI can enhance nucleation by providing a local matching template for ice growth, slits and wedges that provide a good match (e.g. by having a slit width that matches an integer number of ice bilayers) to the structure of the resulting ice crystal also enhance nucleation. I am happy to recommend publication once consideration has been given to the technical issues mentioned below.
1) Polar surfaces such as the AgI (0001) surface are well-known to be unstable in the absence of a polarity compensation mechanism. By keeping the ions in the AgI crystal fixed, this instability is avoided (albeit perhaps artificially), but it is important to consider how the electrostatic boundary conditions might affect the observed properties. Some discussion of these questions would be appreciated. The Sayer and Cox paper already referenced and J. Chem. Phys. 153, 164709 (2020) are interesting in this regard.
2) By keeping the ions in the AgI crystal fixed rather than allowing them to exhibit thermal vibrations around their lattice positions provides a more perfect template for ice nucleation than would be expected if vibrations were allowed. Some discussion of how this potentially affects the results would be appreciated.
Golnaz Roudsari et al.
Golnaz Roudsari et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
267 | 75 | 15 | 357 | 9 | 8 |
- HTML: 267
- PDF: 75
- XML: 15
- Total: 357
- BibTeX: 9
- EndNote: 8
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1