Efficient simulation of second-order phase transitions in quantum anharmonic materials

何十億もの微小な原子で構成された、巨大で複雑なダンスフロアが、音楽の変化に合わせてどのように動くかを予測しようとしていると想像してみてください。材料科学の世界では、この「ダンス」とは、結晶内部における原子の振動を指します。通常、科学者は原子を硬いバネでつながれた小さなボールと仮定することで、その動きを簡単に予測できます。それらは予測可能なリズムで跳ね回っています。これはほとんどの材料において非常にうまく機能します。しかし時として、材料がその形状を変えようとする時（例えば、ある結晶構造から別の構造へと変化する時）、原子は荒れ狂います。彼らは非常に激しく、予測不可能なほどに揺れ動きます。この時、「硬いバネ」という考え方は完全に崩壊してしまいます。これは「二次相転移」と呼ばれ、超高効率の断熱材や、電子スイッチのようにオン・オフを切り替えられる材料の背後にある「秘密のソース」なのです。

問題は、これらのような荒々しく混沌としたダンスをシミュレーションすることは極めて困難であるということです。原子は非常に敏感であるため、標準的なコンピュータ手法を使ってその動きを計算しようとすると、シミュレーションが交通渋滞に巻き込まれたようになり、前進するのに永遠の時間がかかってしまいます。それは、ハリケーンの中にある一枚の葉の軌道を予測しようとするようなものです。ズームインすればするほど、状況はより混沌としていきます。何十年もの間、科学者は最も混沌とした部分を無視するショートカット（近道）を使わざるを得ませんでした。そのため、最も興味深い瞬間を見逃してしまうことが多かったのです。しかし、もし、この交通渋滞に陥ることなく、この混沌としたダンス全体を鮮明に見ることができる新しい種類の望遠鏡を作れるとしたらどうでしょうか？ それこそが、この論文が提示しようとしていることです。

原子の混沌に対する新しい「早送り」ボタン

この論文は、科学者がこれらの荒々しい原子のダンスを、かつてないほど高速かつ正確にシミュレートすることを可能にする、全く新しい数学的なトリックを紹介しています。著者であるアンドレア・バルダンツァとロレンツォ・モナチェッリは、材料がその構造を変えようとするまさにその瞬間に、材料がどのように振る舞うかを予測するという計算上の悪夢を切り抜ける手法を開発しました。

なぜこれが大きなニュースなのかを理解するために、従来の方法を「浜辺の砂粒一つひとつを手作業で地図に書き起こそうとする試み」と考えてみてください。嵐の間にビーチ全体がどのように変化するかを知りたい場合、すべての砂粒が他のすべての砂粒とどのように相互作用しているかを計算しなければなりません。従来のコンピュータコードはこれを行っていましたが、あまりにも遅く、メモリを大量に消費したため、ごく小さなパッチ（例えば50個の原子）しか扱うことができませんでした。もし現実世界の材料の塊（数千の原子）をシミュレートしようとすれば、従来のメソッドではコンピュータを数万年間走らせる必要がありました。それは事実上、不可能だったのです。

著者たちの新しいアプローチは、ビーチの上を飛行するドローンを発明することに似ています。ドローンは、すべての砂粒を一つずつチェックすることなく、砂のパターンを即座に理解できます。彼らは、巨大で乱雑な相互作用のマップを作成する（これはメモリを使いすぎるため）代わりに、問題を一つずつ解くことができる、一連のより小さく単純なパズルとして扱うことができると気づきました。「運動量」（原子の動きをそのリズムに基づいてグループ化する方法）を用いた巧妙な数学的ショートカットを用いることで、必要な作業量を削減したのです。

その結果は驚異的です。かつてスーパーコンピュータを30万年以上走らせる必要があったシミュレーションを、わずか数時間で普通のノートパソコンで実行できるように凝縮しました。また、必要なメモリも、膨大な量の高画質映画を保存するような5ペタバイトから、写真数枚分のサイズであるわずか1ギガバイトへと削減しました。これは単なるスピードアップではありません。これまで不可能だったシミュレーションを可能にする、完全なゲームチェンジャーなのです。

実在する材料による手法の検証

新しい手法が機能することを証明するために、チームはCsSnI3と呼ばれる材料でテストを行いました。これは、太陽電池やその他のグリーンエネルギー技術のために科学者たちが注目している、鉛フリーの結晶です。この材料は「非調和的」であることで有名であり、これは、原子が非常に乱雑で、単純なルールに従わないことを意味します。この材料は、加熱されるにつれて形状を変える相転移を起こしますが、この変化はそれらの激しい原子振動によって引き起こされます。

この新しい高速手法を用いて、チームは1,080個の原子を含むこの材料の超大型セルをシミュレートしました。旧来の手法では、わずか40個の原子を含む小さな塊しか扱うことができませんでした。彼らが結果を調べたところ、驚くべきことが分かりました。小さな40個の原子のシミュレーションは、彼らに嘘をついていたのです。それは材料が実際よりもはるかに安定しているように見せかけ、約450ケルビンまで高温の形状を維持すると予測していました。しかし、新しい手法を用いて1,080個の原子の塊に対して行ったところ、景色が変わりました。材料は実際には不安定になり、約600ケルビンで形状の変化を開始したのです。

この違いは重要です。なぜなら、古いショートカットが材料の真の挙動を見逃していたことを示しているからです。新しいシミュレーションは、4つの原子間で行われる「乱雑な」相互作用（旧来の手法が無視していたもの）が、正しい答えを得るために極めて重要であることを明らかにしました。チームはまた、材料がどのように振動するかを示す指紋のようなものである「ラマン分光スペクトル」についても調査しました。その結果、転移点付近では、振動はもはや整った鋭い音ではなく、広くぼやけた混沌としたものになることが分かりました。これは、原子の「ボールとバネ」の図が、このような条件下では完全に崩壊していることを裏付けており、新しい手法こそが材料の真の混沌とした性質を見る唯一の方法であることを示しています。

なぜこれが重要なのか

この論文は単に「私たちは速くした」と言っているだけではありません。複雑な4原子間の相互作用（「バブル近似」と呼ばれる）を無視するという古い方法が、大規模なシステムに対して誤った答えを導き出すことを明確に示しています。これらのショートカットの必要性を排除することで、著者たちは以前は手の届かなかった材料のシミュレーションへの扉を開きました。

彼らは、従来のレガシーなアプローチでは不可能であったサイズである1,080個の原子のスーパーセルで、この手法が機能することを実証しました。論文は、この新しいツールが、量子効果や原子の混沌が大きな役割を果たす高圧水素や氷のような複雑な材料の真の相図を、科学者がようやく理解することを可能にすると示唆しています。この論文は手法と特定のケースであるCsSnI3に焦点を当てていますが、その含意は明白です。私たちは今、物質の「ダンスフロア」を、以前はSFレベルであった詳細さとスケールでシミュレートできるようになったのです。著者たちは、計算上の壁を管理可能なステップへと変え、研究者が材料の性質が変わる決定的な瞬間を、都市規模のスーパーコンピュータを必要とすることなく探索できるようにしました。