Systematic and accurate anharmonic formation free energies of metastable defects via a constrained Bayesian Adaptive Biasing Force framework

原子という目に見えない世界において、物質が完全に静止していることは決してありません。金属の固形ブロックの中でさえ、原子は振動し、押し合い、時には場所を入れ替えます。放射線や極端な熱によって材料が損傷を受けると、「欠陥」と呼ばれる微細な不完全性が現れます。これらは単なる穴や欠落した破片ではなく、原子が集まりすぎたり、逆に隙間が生じたりといった、複雑な再配置の結果です。科学者たちは、材料が応力下でどのように振る舞うかを予測するためには、これらの欠陥のエネルギーを理解しなければならないことを古くから知っていました。しかし、このエネルギーを計算することは非常に困難です。なぜなら、原子の振動は単純で予測可能なパターンではなく、むしろ「非調和性」として知られる、混沌とした相互接続された動きを示すからです。さらに、多くのこれらの欠陥は「メタステーブル（準安定）」、つまり一時的で不安定な状態にあり、原子がわずかに大きく動くだけで容易に別の状態へと崩壊してしまう性質を持っています。数十年にわたり、この複雑さが、異なる温度におけるこれらの欠陥の挙動に関する完全なライブラリの構築を阻んできました。その結果、原子力発電所や宇宙探査のための材料設計における能力に空白が生じていました。

パリ・サクレー大学の研究チームは、この問題を解決するための実用的な方法を開発しました。彼らは、不安定な欠陥を見失うことなく、その形成自由エネルギーを正確に測定できる計算フレームワークを作成しました。研究では、原子力エネルギーに不可欠な材料である鉄に焦点を当て、4つの余剰な鉄原子がどのように集まるか、また4つの欠損原子（空孔）がどのように配置されるかという、数千通りの異なるパターンを調査しました。課題となったのは、標準的なコンピュータ・シミュレーションがここで失敗してしまうことでした。温度が上昇すると、原子が激しく振動するため、科学者が研究しようとしている特定の欠陥の形状から脱落してしまい、計算が別の無関係な状態へと崩壊してしまうのです。研究者たちは、「制約付きベイズ適応バイアス力法（constrained Bayesian Adaptive Biasing Force）」を用いることで、この問題を克服しました。この手法は、穏やかで見えないガイドのように機能し、原子を調査対象の特定の欠陥形状の境界内に留めつつ、原子が抑え込まれているという事実を数学的に補正します。これにより、コンピュータは欠陥が崩壊することなく、原子の混沌とした高温での振動をサンプリングすることが可能になります。

この研究の成果は、数百もの異なる欠陥構成に関するエネルギー値の、大規模かつ体系的なデータベースです。研究者たちは、原子間の相互作用に関する2種類の異なる数学的モデルを用いて、この手法をテストしました。一つは長年使用されてきた伝統的なモデルであり、もう一つは機械学習に基づいた新しいデータ駆動型のモデルです。その結果、伝統的なモデルはエネルギー地形（エネルギー・ランドスケープ）が粗く不規則であり、特に高温において一貫した結果を得ることが困難であることが分かりました。対照的に、機械学習モデルは非常に滑らかな地形を提供し、手法が極めて高い安定性と精度で機能することを可能にしました。チームは数千回の独立した計算を行い、幅広い温度範囲をカバーしましたが、この手法が失敗したケースは1パーセント未満でした。約1,000個の原子を含むこの規模のシステムにおいて、このレベルの信頼性は前例のないものです。彼らは、エネルギー測定の不確かさが原子あたり0.001電子ボルト未満となるレベルの精度を達成しました。これは、材料の安定性に関する現実世界の予測を行うために必要な精度です。

おそらく最も重要な点は、欠陥の形状の単純さと、そのエネルギーの複雑さとの間に深い関連性があることをこの研究が明らかにしたことです。研究者たちは、原子の混沌とした高温での挙動が、欠陥のより単純な低温での振動を見るだけで、驚くほど正確に予測できることを発見しました。これは、高温における原子運動の複雑で乱雑な実態が、決してランダムなものではなく、低温で見られるより単純なパターンのスケールアップされたバージョンであることを示唆しています。この発見は、極限条件下における欠陥の挙動を推定するために、より単純で高速な計算を用いる道を開くものです。この制約付きの手法が、数千種類の異なる欠陥形状と2種類の原子モデルに対して信頼性を持って機能することを証明することで、研究者たちは、ほぼあらゆる材料に適用できる汎用的なツールを提供しました。この成果は、材料科学における大きなボトルネックを取り除き、材料内の最小の不完全性が、いかにして過酷な環境下での強度と耐久性を決定づけるのかを理解するための明確な道筋を示しています。