分子動力学から見た558環を持つ新規2Dカーボン同素体の熱伝導率

概要

逆非平衡分子動力学（rNEMD）と平衡分子動力学（EMD）シミュレーションを用い、新たに提案された2つの二次元カーボン同素体、オクタゴン-ペンタゴングラフェンライン（OPG-L）とオクタゴン-ペンタゴングラフェンジグザグ（OPG-Z）の熱輸送特性を調査しました。熱伝導率（TC）はサンプルサイズとともに着実に増大し、無限長への外挿でOPG-Lが310–332 W m⁻¹ K⁻¹、OPG-Zが247–228 W m⁻¹ K⁻¹となりました。グラフェン（≈3000–5000 W m⁻¹ K⁻¹）と比較して低いTCは、フォノングループ速度と平均自由行程の低下に起因します。温度と引張歪みが系統的にTCを抑制し、熱管理応用へのチューニング可能性を示唆します。

序論

ダイヤモンド、カーボンナノチューブ、グラフェンなどの炭素同素体は卓越した熱伝導率で知られ、単層ナノチューブとグラフェンは数千W m⁻¹ K⁻¹に達します。最近の理論研究で、グラファイン、オクタグラフェン、T-グラフェン、H-ネットなど、異なる環モチーフを持つ多数の2Dカーボンネットワークが明らかになり、それぞれ独自の電子的・力学的特性を示しています。Suらは八角形と五角形で構成されるエネルギー的に有利で運動論的に安定な2Dシート（558環）を2種類導入。これらOPG-LとOPG-Zの形成エネルギーはそれぞれ原子あたり0.31と0.34 eVで、グラファイン（0.76 eV/原子）より大幅に低く、異方的電子挙動を示します。しかし、熱放散特性は未解明であり、本研究の動機となっています。

モデルと方法

OPG-L（図1a）とOPG-Z（図1b）の原子モデルは繰り返し558リボンから構築、カイラリティはグラフェンに準じてアームチェアとジグザグで定義。LAMMPSを用い、修正最適化TersoffポテンシャルでC-C相互作用を精密記述。得られた平衡格子定数（OPG-L: OA=3.63 Å, OB=9.38 Å; OPG-Z: OA=6.78 Å, OB=5.04 Å）は先行第一原理計算と一致。

熱伝導率計算に逆非平衡分子動力学（rNEMD）を採用。300 KでNosé-Hoover熱浴0.25 ns平衡化後、試料を熱流方向に50スラブ分割。冷スラブ内最も高温原子と熱スラブ内最も低温原子の運動エネルギーを交換し定常熱流生成。結果の温度勾配からフーリエの法則でTC算出：κ = J / (2A ∂T/∂L) （因子2は双方向熱流、厚さはグラフェン層間隔0.34 nmを採用）。

平衡MD（EMD）計算はGreen-Kubo形式で実施。熱流自己相関関数からランニングTC取得、1–2 ns相関時間平均で無限大サイズ値抽出。

結果と考察

サイズ依存性：TCは有効サンプル長（L_eff = L/2）とともに単調増大。50 nm OPG-L（ジグザグ）とOPG-Z（ジグザグ）でTCはそれぞれ125と94 W m⁻¹ K⁻¹、1000 nmで296と236 W m⁻¹ K⁻¹へ上昇。κ⁻¹対L_eff⁻¹の線形フィットで無限大サイズTC外挿：OPG-L 310–332 W m⁻¹ K⁻¹、OPG-Z 247–228 W m⁻¹ K⁻¹、EMD結果（OPG-L 313–344、OPG-Z 261–233 W m⁻¹ K⁻¹）と優れた一致。

温度効果：200–500 KでTC減少、ウムクラップ散乱増大が原因。300 Kから500 KでTC 36–42%低下。相対伝導率傾向はグラフェンと同一傾向、2D格子におけるフォノン-フォノン散乱の普遍性を確認。

歪み効果：熱流方向に単軸引張歪みを印加すると両同素体でTC単調抑制。10%歪みでTC減少率はOPG-Lジグザグ49%、OPG-Zアームチェア31%。振動状態密度（VDOS）解析で面内モードの顕著なフォノン軟化が観察され、熱輸送減少を説明。

力学的特性：単軸引張下でOPG-LとOPG-Zのヤング率はそれぞれ538–648 GPa、492–550 GPa、グラフェンより低いがグラファインと同等。極限歪みはOPG-Lジグザグ17.2%、OPG-Lアームチェア8.7%、堅牢だが歪み感受性高い挙動。

結論

包括的なrNEMDとEMDシミュレーションで、OPG-LとOPG-ZのTCのサイズ、温度、歪み依存性を定量化。フォノン速度と平均自由行程低下でグラフェンより大幅に低いTCだが、歪み工学でチューニング可能。これらの知見は、OPG-LとOPG-Zがマイクロ・ナノスケールデバイスの調整可能な熱管理候補となり得ることを示しています。

略語

558

五-五-八員環

OPG-L

オクタゴン-ペンタゴン グラフェンライン

OPG-Z

オクタゴン-ペンタゴン グラフェンジグザグ

rNEMD

逆非平衡分子動力学

TC

熱伝導率

VDOS

振動状態密度