[論文レビュー] Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation

Miguel Camarillo, J. Oller-Iscar+3 more|arXiv (Cornell University)|Jan 10, 2026

nanoparticles nucleation surface interactions被引用数 0

ひとこと要約

研究は、格子不整合、基板方位、および格子の柔軟性が水基 substrates と mW モデルを用いた不均質氷核生成に与える影響を定量化します。1% の不整合あたり約4 Kの核生成温度の低下、氷の方位間でほぼ同等の核生成能力、柔軟な基板の方が核生成効率が高いことを示します。

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 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 one per cent increase of structural mismatch leads to a decrease of approximately 4 K in the nucleation temperature. We also analyse the effect of the 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 asses the effect of lattice flexibility by comparing substrates where molecules are immobile with 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.

研究の動機と目的

• 格子不整合が氷核生成に及ぼす影響を分子間相互作用の影響から独立して分離・定量化する。
• 異なる氷-基板方位が核生成能力に与える影響を決定する。
• 基板格子の柔軟性が核生成効率に及ぼす役割を評価する。
• 不整合基板上での核生成における界面構造の機構的洞察を提供する。

提案手法

• LAMMPSのNVT系におけるmW水モデルを用い、引張られた水-基板上で氷核生成をシミュレートする。
• 共存状態で氷単位胞を伸縮して基板を構築する；構造的不整合 delta = 100*|f-1| を定義する。
• 2つの基板タイプを比較する：硬直（ immobile ）と wells（原子はポテンシャル井戸内を移動可能）。
• 核生成を定量化するために誘導時間を測定し J = 1/(2 A t_ind) を算出し、log(J/(m^2 s)) = 23.6 で T_n を定義する。
• 界面近傍のスラブに沿った径方向分布関数と q12 局所結合秩を用いて界面構造を分析する。
• 基板を曝露する3つの氷方位を検討する：basal、primary prism (pI)、secondary prism (pII)。

Figure 3: (a) Potential energy versus time for a trajectory at 272 K (1 K below the melting temperature) on a wells substrate with $\delta=0$ . (b) Potential energy versus time for nine trajectories of a liquid at 245.0 K on a wells substrate with $\delta=7$ . (c) Time evolution of the potential ene

実験結果

リサーチクエスチョン

• RQ1 氷と基板間の格子不整合は一定の面積と観測時間の下で核生成温度にどのような影響を及ぼすか？
• RQ2 3つの主要氷方位（basal、pI、pII）は、引張られた水基板上で類似の核生成能力を示すか？
• RQ3 基板の格子柔軟性は不均質氷核生成にどのような影響を与えるか？
• RQ4 等方的な引張りと圧縮が核生成温度と速度にどう影響するか？
• RQ5 不整合基板での核生成過程における界面近傍の液体構造と核の構造的特徴は何か？

主な発見

δ	A_m /nm^2	N_s	N_l	L_s	n_xn_yn_z	Plane	T_n / K
5	137.1	5400	30600	13.943	10x15x15	pII	254.6
7	142.3	5400	30600	14.208	10x15x15	pII	244.4
8	145.0	5400	30600	14.341	10x15x15	pII	241.4
5	137.1	5400	30600	9.295	10x15x15	pII	255.8
7	142.3	5400	30600	9.472	10x15x15	pII	245.1
8	145.0	5400	30600	9.561	10x15x15	pII	242.9
5	168.5	5400	30600	10.399	15x15x10	basal	255.5
5	158.3	5400	30600	10.737	15x10x15	pI	257.0
5	137.1	5400	30600	9.295	10x15x15	pII	256.0
7	174.9	5400	30600	10.597	15x15x10	basal	246.0
7	164.3	5400	30600	8.208	15x10x15	pI	247.0
7	142.3	5400	30600	9.472	10x15x15	pII	246.0
8	178.2	5400	30600	7.779	15x15x10	basal	243.0
8	167.3	5400	30600	8.285	15x10x15	pI	244.0
8	145.0	5400	30600	9.561	10x15x15	pII	243.5
8	178.2	5400	30600	7.779	15x15x10	basal	243.0

• structural mismatch が1%増えるごとに核生成温度が約4 K低下（硬直基板と wells 基板の両方で） 。
• mismatches に対する核生成速度-温度曲線は平行で、deltaと共存点（delta=0）へ外挿可能な直線状の傾向を形成する。
• 3つの氷方位（basal、pI、pII）は非常に類似した核生成能力を示し、効率には小さな差がある。
• 柔軟性のある wells 基板は硬直基板よりも高温（約1–2 K）で氷を核生成する；これは基板が氷構造に適応できるため。
• 界面近傍の液体構造は不整合が小さくなるほど氷により近づき、均一なウェットニング層よりも不均質で斑状な氷様領域が観察される；核は基板からbulk iceへと徐々に回復するような斜めの分子列を示す。

Figure 4: (a) Heterogeneous nucleation rate versus temperature for the different mismatches and the different kinds of stretched substrates studied in this work as indicated in the legend. These results correspond to the pII orientation (exposing the yz plane of the stretched and replicated unit cel