Hydrogen-Bond Density Drives the Order-Disorder Phase Transition from Orientational to Vibrational in Polyalcohol Plastic Crystals

技術要約：ポリアルコール・プラスチック結晶における水素結合密度による秩序–無秩序相転移の駆動

問題提起
ネオペンタン誘導体のポリアルコール・プラスチック結晶は、極めて大きなエントロピー変化（$\Delta S_t$）を伴う秩序–無秩序（OD）相転移を示し、熱エネルギー貯蔵やバロカロリック冷却の有望な候補となっている。これらの転移には、配向ダイナミクスと水素結合が密接に関連していることは確立されているが、$\Delta S_t$ の微視的な起源は依然として主に定性的な理解に留まっている。重要な未解決の問いは、$\Delta S_t$ に対する配向、コンフォメーション、および振動の寄与のバランスが、個々の化合物の偶発的な特性なのか、それとも分子設計を通じて系統的に調整可能な予測可能な傾向に従うものなのか、という点である。特に、水素結合ネットワークが改変されるにつれて振動の寄与がどのように進化するか、また、この経路を圧力堅牢性を高めるためにエンジニアリングできるかどうかが不明である。

手法
本研究では、古典的分子動力学（MD）シミュレーション、情報理論、および熱力学を組み合わせたアプローチを用いて、NPG–PG–PEシリーズ（ネオペンチルグリコール：NPG（ヒドロキシル基2個）、ペンタグリセリン：PG（ヒドロキシル基3個）、ペンタエリスリトール：PE（ヒドロキシル基4個））を解析した。

• シミュレーション: バルク系（512分子）に対し、異なる圧力（0–0.5 GPa）および温度条件下で、CHARMM型力場を用いたNPT-MDシミュレーションをGROMACSを用いて実施した。相転移温度（$T_t$）は、比体積の不連続な跳躍によって特定した。
• エントロピー分解: 全転移エントロピーを、振動成分（$\Delta S_{vib}$）、配向成分（$\Delta S_{ori}$）、およびコンフォメーション成分（$\Delta S_{conf}$）の3つの成分に分解した。
• 構造検証: PGの秩序相に関するX線回折による精密化を行い、従来の正方晶対称性の仮定を修正し、三斜晶構造（空間群 $P\bar{1}$）であることを示した。
• 微視的解析: 頻度空間におけるエントロピー寄与を解明するために、振動密度状態（VDOS）を解析した。さらに、分子間のO–H···O接触の幾何学的解析を行い、秩序相および無秩序相における水素結合（HB）ネットワークの密度と連結性を定量化した。

主な結果

1. エントロピー階層の逆転: 本研究は、水素結合の連結性が増加するにつれて、相転移を駆動する主要なメカニズムが完全に逆転することを明らかにした。
   • NPG: 配向の無秩序化が $\Delta S_t$ の主要な寄与因子（~44%）であり、次いでコンフォメーションおよび振動成分が続く。
   • PG: 配向、コンフォメーション、および振動の各寄与がほぼ同等（それぞれ~32–36%）となる、バランスの取れた中間領域を占める。
   • PE: 格子振動が主要な寄与因子（47%）となり、配向の無秩序化は最小の割合（21%）に低下する。
2. 圧力堅牢性: 振動チャネルは、配向チャネルよりも静水圧に対して著しく堅牢である。NPGでは圧縮に伴い振動成分の相対的な重みが減少するのに対し、PEでは、検討された全範囲において振動成分が支配的かつ最も圧力安定な項として維持される。
3. 振動密度状態（VDOS）: この逆転現象は、秩序相と無秩序相の間の低周波領域におけるエントロピーギャップに起因する。水素結合ネットワークが緻密化する（NPGからPEへ）につれて、秩序相と無秩序相のVDOSの分離、特に低周波領域（<8 THz）での分離が拡大する。無秩序相は、HBネットワークの崩壊に伴う低エネルギー格子モードの放出により、低周波領域において振動エントロピーをより急速に蓄積する。
4. 水素結合ネットワークのダイナミクス: 分子あたりの水素結合数（$n_{HB}$）は、ヒドロキシル基の数（$n_{OH}$）とともに増加する。しかし、無秩序相への転移に伴う水素結合の「喪失」量（$\Delta n_{HB}$）も、NPGからPEにかけて大幅に増加する。PEは、水素結合の絶対量および相対量の両方で最大の喪失を示し、これが低周波振動モードの最大の放出と相関している。
5. 対称性の影響: 配向エントロピーの寄与の減少は、分子対称性の増加に起因する。メチル基がヒドロキシメチル基に置換されるにつれ、面心立方格子と適合する区別可能な分子配向の数が減少し、高温相における配向相空間が減少する。

意義および主張
本論文は、OD相転移のエントロピーの微視的な性質は偶発的なものではなく、水素結合密度によって支配される予測可能な傾向に従うことを確立している。著者らは、格子振動を単なる微小な補正としてではなく、水素結合を持つプラスチック結晶におけるOD転移の中心的かつ調整可能な駆動因子として再定義すべきであると主張している。

その主要な意義は、水素結合密度を、実用的かつ化学的にアクセス可能な設計パラメータとして特定した点にある。ヒドロキシル基の数を調整することで、特定の設計ハンドルを持つ材料を設計できる：

• 振動支配型（高いHB密度）は、外部圧力に対して堅牢なカロリック応答を持つ材料への道を開く。
• 配向支配型（低いHB密度）は、高度に圧力調整可能な応答を持つ材料への道を開く。

本研究は、局所的な化学的連結性（ヒドロキシル基の数）と集団的な格子ダイナミクスの間の微視的なリンクを提供し、熱エネルギー貯蔵およびバロカロリック材料を最適化するための新しい設計原理を提示している。これらの関係性は、他の水素結合プラスチック結晶のファミリーにも拡張可能であり、将来的な新材料探索のためのスクリーニング基準となることが示唆される。