[论文解读] Toward a Comprehensive Model of Snow Crystal Growth: 4. Measurements of Diffusion-limited Growth at -15 C

Kenneth G. Libbrecht, Christopher W. Miller+4 more|arXiv (Cornell University)|Dec 10, 2015

nanoparticles nucleation surface interactions参考文献 20被引用 6

一句话总结

本文提出了一套全面的-15 °C下扩散限制的冰晶生长动力学模型，通过受控实验与二维元胞自动机模拟，再现了从块状柱状体到薄板状结构的形态转变。关键贡献在于证明了棱面附着系数取决于晶面宽度，当晶面宽度达到分子尺度时，成核势垒显著降低，从而能够准确预测不同过饱和度下的生长速率与形态。

ABSTRACT

We present measurements of the diffusion-limited growth of ice crystals from water vapor at different supersaturation levels in air at a temperature of -15 C. Starting with thin, c-axis ice needle crystals, the subsequent growth morphologies ranged from blocky structures on the needle tips (at low supersaturation) to thin faceted plates on the needle tips (at high supersaturation). We successfully modeled the experimental data, reproducing both growth rates and growth morphologies, using a cellular-automata method that yields faceted crystalline structures in diffusion-limited growth. From this quantitative analysis of well-controlled experimental measurements, we were able to extract information about the attachment coefficients governing ice growth under different circumstances. The results strongly support previous work indicating that the attachment coefficient on the prism surface is a function of the width of the prism facet. Including this behavior, we created a comprehensive model at -15 C that explains all the experimental data. To our knowledge, this is the first demonstration of a kinetic model that reproduces a range of diffusion-limited ice growth behaviors as a function of supersaturation.

研究动机与目标

• 开发一个从水蒸气中生长冰晶的综合性模型，能够定量再现受控条件下冰晶的生长速率与形态。
• 研究冰晶形态转变的物理起源，特别是不同过饱和度下从柱状向板状形态的转变机制。
• 确定棱面附着系数如何随晶面宽度变化，尤其是在分子尺度尺寸下的变化。
• 通过在双扩散室中精确测量冰针生长的实验数据，对模型进行验证，实验条件包括温度与过饱和度的精确控制。
• 建立一个动力学框架，解释冰生长过程中观察到的边缘锐化不稳定性，其机制源于宽度依赖的成核势垒。

提出的方法

• 实验在双扩散室中进行，用于在-15 °C下生长c轴取向的冰针晶体，并以±0.1 °C的精度控制过饱和度。
• 线性温度梯度抑制了对流，实现了对过饱和度的精确建模，生长过程通过光学显微镜以10–20秒的间隔进行监测。
• 采用二维轴对称的元胞自动机模型模拟扩散限制下的生长，通过动力学蒙特卡罗方法引入具有取向性的晶体结构。
• 模型通过参数$ R_{ESI} $引入棱面附着系数的宽度依赖性，其中$ \sigma_{0,prism} = \sigma_{0,prism,\infty}(1 - \exp(-R_c/R_{ESI})) $，且$ R_c $与晶面宽度成正比。
• 模型参数通过实验数据进行标定，包括来自先前测量的$ \sigma_{0,prism,\infty} = 3\% $，以及每分子水约$ a \approx 0.3 $ nm的晶面宽度换算系数。
• 通过调节过饱和度和$ R_{ESI} $，模型成功再现了从块状柱状体到薄板状结构的形态转变，与观测到的生长速率和形状高度一致。

实验结果

研究问题

• RQ1 在-15 °C下，随着过饱和度增加，冰晶的生长形态如何演变？
• RQ2 棱面宽度在决定水蒸气在冰表面的附着系数中起什么作用？
• RQ3 一个单一的动力学模型能否同时再现扩散限制下冰生长的速率与形态转变？
• RQ4 当棱面宽度趋近分子尺寸时，为何其成核势垒会急剧降低？
• RQ5 宽度依赖的附着系数如何解释观测到的冰晶生长中边缘锐化不稳定性？

主要发现

• 模型成功再现了当过饱和度从1.7%增加到5%时，冰针尖端从块状柱状体向薄板状结构的形态转变。
• 棱面附着系数随晶面宽度减小而降低，在分子尺度宽度下趋近于零，表明在原子台阶处附着效率接近100%。
• 参数$ R_{ESI} $（控制$ \sigma_{0,prism} $的宽度依赖性）经标定后，可准确再现三种不同过饱和度水平下的实验生长形态。
• 当棱面宽度超过约30个分子（~9 nm）时，$ \sigma_{0,prism} $趋近于3%，与[22]中先前的数据一致。
• 模型表明，边缘锐化不稳定性源于窄棱面成核势垒的急剧降低，从而解释了向板状生长的转变机制。
• 这是首个能够定量再现不同过饱和度水平下一系列扩散限制冰生长行为的动力学模型，为冰晶生长动力学提供了统一的理论框架。