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

振動状態密度