🔬 condensed matter

The origin of complex behavior of liquid carbon: an insight from computer simulation

本論文では、機械学習ポテンシャル GAP-20 を用いた分子動力学シミュレーションにより、液体炭素の複雑な挙動の起源が解明され、そのモデルにおける比較的低い臨界温度がグラファイトの融点の実験測定結果に影響を与える可能性が示されました。

原著者： Yu. D. Fomin

公開日 2026-03-02

📖 1 分で読めます☕ さくっと読める

原著者： Yu. D. Fomin

原論文は CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

1. 問題点：なぜ温度の答えがバラバラなの？

黒鉛（鉛筆の芯やダイヤモンドの原料）を熱すると、いつ溶けるのか？
これまでの実験では、**「4000 度で溶ける」という人もいれば、「6700 度で溶ける」**という人もいて、答えが全く一致しません。

これまでは、実験のやり方の違いや測定ミスのせいだと思われていました。しかし、この論文の著者（ドミトリー・フォミン氏）は、**「もしかして、黒鉛が溶ける瞬間に『ある特別な状態』に陥っているから、実験が難しくなっているのではないか？」**と考えました。

2. 発見：「临界点（きょうかいてん）」という魔法の境界線

液体と気体の間には、**「临界点（クリティカルポイント）」**という不思議な境界線があります。

水の場合： 100 度で沸騰しますが、それ以上圧力をかけると、液体と気体の区別がつかない「超臨界流体」という状態になります。
黒鉛の場合： この論文では、**「黒鉛の液体と気体の境界（临界点）が、実は黒鉛が溶ける温度（融点）のすぐそばにある」**という可能性を突き止めました。

【イメージ】
Imagine 氷が溶けて水になり、さらに蒸気になるところを想像してください。
通常、氷→水→蒸気と明確に分かれます。
しかし、黒鉛の場合は、**「氷が溶け始めた瞬間、すでに蒸気と液体が混ざり合った『どろどろの霧』のような状態」**になっている可能性があります。

3. シミュレーションの結果：「揺らぎ」が実験を狂わせる

著者は、AI（機械学習）を使って原子レベルの動きをシミュレーションしました。その結果、以下のようなことがわかりました。

临界点が低い： 黒鉛の液体と気体の境界は、予想よりずっと低い温度（約 4200 度）にあり、黒鉛が溶ける温度（約 4000〜6000 度）と重なり合っている。
激しい揺らぎ： この境界付近では、物質の性質（熱の吸収具合や体積の膨らみ方）が**「激しく揺れ動いている」**。

【アナロジー：揺れる橋】
実験室で黒鉛を測ろうとするのは、**「激しく揺れている吊り橋の上で、静かに定規で長さを測ろうとする」**ようなものです。

橋が揺れれば、測るたびに長さが変わります。
熱の吸収具合も、揺れによって「すごく熱を吸う瞬間」と「あまり吸わない瞬間」が混ざり合います。

過去の研究者たちは、「熱の吸収具合は一定だ」という前提で計算していましたが、実は**「揺れている最中」**だったので、前提が崩れてしまい、実験結果がバラバラになったのです。

4. なぜこれが重要なのか？

この研究は、**「実験結果がバラバラなのは、実験者のミスではなく、物質そのものが『超臨界状態』という不安定な領域にいるからかもしれない」**と示唆しています。

実験の難しさ： 温度を少し変えるだけで、液体と気体の割合が激しく変わり、密度や熱の性質が予測不能になります。
新しい視点： これまでの「一定の値」という考え方を捨て、**「揺らぎ（フラクチュエーション）」**を考慮して実験データを解釈する必要があると提言しています。

5. 結論：まだ完全な答えではないが、大きな一歩

著者は正直に言っています。「今回のシミュレーションに使った AI モデルは、固体（黒鉛）の性質を学ぶために作られたもので、気体や液体の極限状態を完璧に再現しているわけではない」と。

しかし、**「临界点が融点のすぐそばにある可能性」**というアイデアは、長年の謎を解くための重要な鍵になるでしょう。

まとめると：
黒鉛が溶ける温度がなぜわからないのか？
それは、「溶けかけの黒鉛が、液体と気体の『境界線』で激しく揺れているから」。
まるで**「沸騰したお湯の泡が弾ける瞬間」**のように、安定した状態ではなく、常に揺れ動いているため、正確な温度を測るのが極めて難しい、というのがこの論文のメッセージです。

一言で言うと：
「黒鉛の融点の謎は、**『溶ける瞬間に物質が激しく揺れ動いているから、測るのが難しい』**という新しい視点で解決できるかもしれない！」という、物理学界への挑戦状です。

以下は、Yu. D. Fomin 氏による論文「The origin of complex behavior of liquid carbon: an insight from computer simulation（液体炭素の複雑な挙動の起源：コンピュータシミュレーションからの洞察）」の技術的な要約です。

1. 研究の背景と課題 (Problem)

黒鉛（グラファイト）の融点は物理学における長年の未解決問題の一つです。

実験結果の不一致: 20 世紀初頭からの実験以来、黒鉛の融点に関する推定値は $T_m = 4000\,\text{K}$ から $6700\,\text{K}$ まで広範に散らばっており、合意が得られていません。
既存モデルの限界: 従来の臨界点（液体 - ガス転移）の推定（Leider らによる 1973 年の研究など）は、融解時の体積変化や熱膨張係数に関する仮定に依存しており、これらの仮定が実験結果と矛盾している可能性があります。
未考慮の要因: 多くの実験解析において、高温領域における気相（昇華した炭素）の影響や、臨界点近傍での熱力学的な揺らぎ（比熱や熱膨張係数の急増など）が十分に考慮されていません。特に、融解直後の状態が「液体」なのか「超臨界流体（SCF）」なのかの区別が不明確です。

2. 手法 (Methodology)

本研究では、機械学習ポテンシャル GAP-20 を用いた分子動力学（MD）シミュレーションを実施しました。

モデル: 4096 個の炭素原子を含む立方体ボックス（周期的境界条件）。
ポテンシャル: 密度汎関数理論（DFT）計算に基づいて学習された GAP-20 ポテンシャルを使用。
シミュレーション条件:
- 正準集団（NVT）で実施。
- 密度範囲： $0.05 \sim 2.0\,\text{g/ml}$ （ステップ $0.05\,\text{g/ml}$ ）。
- 温度範囲： $4000 \sim 5000\,\text{K}$ （ステップ $200\,\text{K}$ ）に加え、 $4100, 6000, 7000\,\text{K}$ の等温線も計算。
- 平衡化時間：10 ps、生産計算時間：20 ps。
解析手法:
- 圧力と内部エネルギーのデータを多項式関数（ $X = \sum a_{ij} T^i \rho^j$ ）にフィットさせ、状態方程式（EoS）を構築。
- 等温線から臨界点（変曲点）を特定。
- 熱力学的応答関数（定圧比熱 $c_P$ 、熱膨張係数 $\alpha_P$ 、等温圧縮率 $\beta_T$ ）を計算し、その極大値（ウィドム線：Widom lines）の位置を特定。

3. 主要な結果 (Results)

臨界点の推定:
- GAP-20 モデルに基づく計算結果、臨界点は $T_c \approx 4235\,\text{K}$ 、 $\rho_c \approx 0.189\,\text{g/ml}$ 、 $P_c \approx 107\,\text{bar}$ と求められました。
- これは従来の推定値（ $T_c \approx 6810\,\text{K}$ ）よりも著しく低い値ですが、黒鉛の融点（実験値で $4000\sim6000\,\text{K}$ 程度）と非常に近接しています。
臨界点近傍の挙動:
- 臨界点の近くで、定圧比熱 $c_P$ 、熱膨張係数 $\alpha_P$ 、等温圧縮率 $\beta_T$ が極大値を示すことが確認されました（ウィドム線）。
- 特に $c_P$ は臨界点近傍で粒子あたり $25 k_B$ に達し、急激に減少します。
- $T=4100\,\text{K}$ 、 $\rho=1.55\,\text{g/ml}$ 付近では、液体から黒鉛相への自発的な結晶化が観測され、臨界点が融解線に極めて近いことが示唆されました。
実験的不整合のメカニズム:
- 黒鉛の融解直後、系は超臨界流体（SCF）状態に入り、臨界点近傍の強い熱力学的揺らぎの影響を受けます。
- これにより、比熱や密度などの物理量が大きく変動し、実験測定（高温計測や密度測定）に大きな誤差をもたらす可能性があります。
- 従来の実験解析で「比熱は一定」「熱膨張係数は黒鉛と同等」といった仮定が用いられていますが、臨界点近傍ではこれらが完全に破綻していることが示されました。

4. 考察と限界 (Discussion & Limitations)

GAP-20 モデルの限界: GAP-20 は凝縮相（固体・液体）の再現に焦点を当てて学習されたため、希薄な気相や超臨界流体の記述には限界がある可能性があります。したがって、算出された絶対的な臨界点の値（ $4235\,\text{K}$ ）は実物質の値と一致しないかもしれませんが、「臨界現象が融解実験に重大な影響を与える」という定性的な結論は有効です。
DFT との課題: 臨界点の計算は DFT 自体の課題（高温での精度、分散力の扱い、臨界揺らぎのサンプリング難易度）とも深く関連しており、高精度な計算は依然として困難です。
実験との対比: 最近の実験（Ivanov et al., 2026）では、融点付近の液体炭素の定積比熱が Dulong-Petit の法則を大きく上回る値を示しており、これは臨界揺らぎまたは構造のクロスオーバーによるものと考えられます。

5. 意義と結論 (Significance & Conclusion)

実験結果の散らばりの説明: 黒鉛の融点に関する実験結果の大きなばらつきは、融解直後の系が臨界点（またはその近傍）に位置し、熱力学的な揺らぎが極大化していることによる影響が、実験データの解析で無視されていることに起因する可能性が高いと提唱しました。
今後の展望: 固体、液体、気相のすべての性質を同時に記述できる新たな炭素モデルの構築が必要であり、それによってより精密なシミュレーションと実験データの整合性が図られると結論付けています。

総括:
本論文は、機械学習ポテンシャルを用いたシミュレーションを通じて、黒鉛の融解過程が単なる「固体から液体への転移」ではなく、臨界点近傍の超臨界流体領域を含む複雑な過程であることを示唆しました。臨界点近傍での熱力学的応答関数の急激な変化が、従来の実験測定値の不一致や誤差の主要原因である可能性を指摘した点が、この研究の最大の貢献です。