[논문 리뷰] Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation

Miguel Camarillo, J. Oller-Iscar외 3명|arXiv (Cornell University)|2026. 01. 10.

nanoparticles nucleation surface interactions인용 수 0

한 줄 요약

연구는 격자 불일치, 기판 방향성, 격자 유연성이 물 기반 기판과 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.

연구 동기 및 목표

• 격자 불일치가 분자 간 상호작용과 무관하게 아이스 핵생성에 미치는 효과를 분리하고 정량화한다.
• 확산된 물-기판 방향에서 핵생성 능력에 미치는 아이스-기판 방향성을 결정한다.
• 핵생성 효율에 대한 기판 격자 유연성의 역할을 평가한다.
• 불일치 기판에서의 핵생성 동안 계면 구조에 대한 기계적 통찰력을 제공한다.

제안 방법

• NVT 엔밴에서 LAMMPS의 mW 물 모델을 사용하여 확장된 물-기판에서 아이스 핵생성을 시뮬레이션한다.
• 공존에서 얼음 단위 셀을 늘리거나 압축하여 기판을 구성하고 구조적 불일치 delta = 100*|f-1|를 정의한다.
• 두 가지 기판 타입을 비교한다: 강체(움직이지 않음)와 wells(원자가 포텐셜 웰 내에서 움직일 수 있음).
• induced 시간 측정과 J = 1/(2 A t_ind)을 계산하여 핵생성을 정량화하고 log(J/(m^2 s)) = 23.6에서 T_n을 정의한다.
• 계면 근처의 슬래브를 따라 반경 분포 함수와 q12 지역 결합 순서를 이용해 계면 구조를 분석한다.
• 세 가지 아이스 방향성(기판 노출 방향): 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 세 가지 주요 아이스 방향성(basal, pI, pII)이 늘어난 물 기판에서 유사한 핵생성 능력을 보이는가?
• RQ3 기판 격자 유연성이 비균질 아이스 핵생성에 어떤 영향을 주는가?
• RQ4 기판의 등방성 늘리기 vs 압축이 핵생성 온도와 속도에 어떤 영향을 미치는가?
• 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

• 구조적 불일치가 1% 증가할 때 두 경우의 기판(강체 및 wells)에서 핵생성 온도가 대략 4 K 감소한다.
• 핵생성 속도 대비 온도 곡선은 불일치 간에 평행하고, delta에 대한 T_n의 선형 유사 경향이 공존점(delta=0)으로 외삽되게 한다.
• 세 가지 아이스 방향성(basal, pI, pII)은 매우 유사한 핵생성 능력을 보이며 효율성에 작은 차이가 있다.
• 유연한(wells) 기판은 강체 기판보다 더 높은 온도에서 얼음을 핵생성하는데, 이는 기판이 얼음 구조에 적응할 수 있기 때문이며 약 1–2 K의 차이가 있다.
• 계면 근처의 액상 구조가 불일치가 감소할수록 더 아이스와 유사해지고, 균일한 젖음층이 아닌 이질적이고 패치된 아이스 유사 영역이 나타나며 핵은 기판에서 벌크 아이스까지 점진적 구조 회복을 나타내는 기울어진 분자 열으로 보인다.

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