[論文レビュー] Mass-density and Phonon-frequency Relaxation Dynamics of Water and Ice at Cooling

Changqing Sun|arXiv (Cornell University)|Oct 5, 2012

Material Dynamics and Properties参考文献 19被引用数 78

ひとこと要約

本稿では、水のO:H-O水素結合における結合電子対と非結合電子対の間のクーロン反発力に加え、O:HとH-O部分間の比熱差が、冷却時の異常な密度変化およびフォノン周波数の緩和ダイナミクスを支配することを提案する。液体および氷においてO:H結合は「マスター」として機能し、冷却に伴い収縮し剛性を増す一方、H-O結合は「スレーブ」として伸びて柔らかくなる。この水素結合の動的逆転が、氷の低密度および浮遊性を説明する。

ABSTRACT

Coulomb repulsion between the bonding electron pair in the H-O covalent bond (denoted by "-") and the nonbonding electron pair of O (":") and the specific-heat disparity between the O:H and the H-O segments of the entire hydrogen bond (O:H-O) are shown to determine the O:H-O bond angle-length-stiffness relaxation dynamics and the density anomalies of water and ice. The bonding part with relatively lower specific-heat is more easily activated by cooling, which serves as the "master" and contracts, while forcing the "slave" with higher specific-heat to elongate (via Coulomb repulsion) by different amounts. In the liquid and solid phases, the O:H van der Waals bond serves as the master and becomes significantly shorter and stiffer while the H-O bond becomes slightly longer and softer (phonon frequency is a measure of bond stiffness), resulting in an O:H-O cooling contraction and the seemingly "regular" process of cooling densification. In the water-ice transition phase, the master and the slave swap roles, thus resulting in an O:H-O elongation and volume expansion during freezing. In ice, the O--O distance is longer than it is in water, resulting in a lower density, so that ice floats.

研究の動機と目的

• 277 Kで密度最大を示す水の性質と、冷却にもかかわらず密度が低い氷の性質という、長年の謎を解明すること。
• 液体および固体水における通常の冷却収縮と、凍結時の異常な膨張を対比して説明すること。
• フォノン周波数のシフトと巨視的密度異常を結びつけるために、水素結合の角度・長さ・剛性の緩和ダイナミクスの役割を明確にすること。
• 電子対間のクーロン反発力とO:H-O結合における部分別比熱容量差に基づく統一的枠組みを確立すること。
• 全相態における水および氷の構造的緩和ダイナミクスと、測定されたラマンおよび赤外フォノン周波数シフトを関連付けること。

提案手法

• 水素結合を2つの部分に分割するモデルを提案：O:H（ファンデルワールス的）とH-O（共有結合的）、それぞれ比熱応答が異なる。
• デバイモデルを適用して部分別比熱容量（O:Hのη_L、H-Oのη_H）を推定し、η_Lがη_Hよりも温度上昇に伴い速やかに増加することを示した。
• 分子動力学（MD）シミュレーションを用いて、300 Kから4 Kに冷却する過程におけるO:H-O結合の緩和ダイナミクスをモデル化した。
• 結合電子対（H-O）と非結合電子対（O:）の間のクーロン反発力（f_q）を、非対称的結合緩和の駆動力として導入した。
• ラマン分光法のデータとMD予測を用いて、フォノン周波数シフト（O:Hのω_L、H-Oのω_H）と結合長および剛性の変化を相関させた。
• 温度領域を3つに分けた：T > 273 K（I）、273–258 K（II）、T < 258 K（III）、それぞれ異なる緩和行動を特定した。

実験結果

リサーチクエスチョン

• RQ1 なぜ水は277 Kで密度最大を示し、凍結するとほとんどの物質とは逆に体積が膨張するのか？
• RQ2 液体および固体水の冷却過程で観察されるH-Oフォノン周波数（ω_H）の赤方シフトと、凍結時のブルーシフトの原因は何か？
• RQ3 O:HとH-O部分間の比熱差が、結合長および剛性の緩和の方向性と大きさに与える影響は何か？
• RQ4 結合電子対と非結合電子対の間のクーロン反発力が、O:H-O水素結合における非対称的緩和を駆動する役割を果たす理由は何か？
• RQ5 水-氷転移時にO:HとH-O結合のマスター・スレーブ役割が逆転する理由は何か？これにより体積膨張が生じる仕組みは？

主な発見

• O:H結合（ファンデルワールス的）は液体および固体相において「マスター」として機能し、H-O結合の伸びよりも強く収縮し剛性を増すため、冷却による収縮が生じる。
• H-O結合（共有結合的）は「スレーブ」として機能し、O:H結合からのクーロン反発により伸びて柔らかくなり、液体相ではω_Hが約3200 cm⁻¹から約3150 cm⁻¹に赤方シフトする。
• 凍結過程（273–258 K）ではマスターとスレーブの役割が逆転する：O:H結合は伸びて柔らかくなり（ω_Lが220 cm⁻¹から215 cm⁻¹に赤方シフト）、H-O結合は収縮し剛性を増す（ω_Hが3150 cm⁻¹から3170 cm⁻¹にブルーシフト）、これが体積膨張を引き起こす。
• 氷では液体水よりもO–O距離が長くなるため、密度が9%低くなり、氷が浮遊する理由が説明される。
• 氷におけるω_Hフォノン周波数は、209 Kで3253 cm⁻¹から30 Kで3218 cm⁻¹に低下し、ナノスケールの氷クラスターの赤外分光法と整合的である。
• MDシミュレーションの動画は、水素結合の動的スイッチングと液体水中に持続するクーロン反発を確認しており、熱的フラクチュエーションを伴う四面体的配位構造の存在を支持する。