TIP4P/Iceモデルを用いて250Kでの多結晶氷の古典的分子動力学シミュレーションを行った。多結晶氷の構造は、過冷却水中で氷粒子を成長させることで作成した。最近開発された次数パラメータを用いて、液-液相転移シナリオの観点から局所構造を特徴付けることに成功した。その結果、氷の結晶粒界と三重接合は、ほとんどの水分子が4つの水素結合を形成し、O-O-O角が109.47°の四面体角からずれている低密度の液体水と構造的に類似していることが示された。結晶粒界の厚さは∼1nmである。本研究で計算された粒界に沿った水分子の拡散係数は5.0×10-13m2 s-1であり、実験データとよく一致している。また、トリプルジャンクションに沿った拡散は、粒界に沿った拡散よりも3.4倍速いことがわかった。シミュレーション結果を用いて多結晶氷中の水分子の拡散率の粒径依存性をモデル化したところ、ミリオーダーの粒径を持つ典型的な多結晶氷試料では、粒界と三重接合が拡散率に与える影響は無視できることがわかった。また、結晶粒界の特性は、同じ温度でも氷/蒸気界面の特性とは大きく異なることも明らかにした。 (DeepLによる機械翻訳)

Takuma Yagasaki, Masakazu Matsumoto, and Hideki Tanaka, Molecular dynamics study of grain boundaries and triple junctions in ice, J. Chem. Phys. 153, 124502 (2020). https://doi.org/10.1063/5.0021635]