Crystallisation kinetics of supercooled liquid palladium

本研究では、古典的分子動力学シミュレーションを用いて過冷却液体パラジウムの結晶化速度論を特性化し、拡散律速成長と、時間分解X線回折実験と一致する融点付近$0.5 T_{\mathrm{m}}$における均質核生成の極大を明らかにしており、これは急速にクエンチされたPd薄膜における達成可能な過冷却が均質核生成によって支配されていることを示している。

要約

パラジウム（光沢のある銀白色の金属）の溶融した金属のポットを想像してみてください。もし、これをゆっくりと冷却させると、水が氷に変わるように、自然に結晶へと戻ろうとします。しかし、もし、結晶化する暇を与えないほど信じられないほど速く冷却することができたらどうなるでしょうか？結晶になる代わりに、乱れた無秩序な状態で「凍結」され、金属ガラスになります。

この論文は、パラジウムが過冷却状態において、いかに速く結晶に戻ろうとするのか、そして、結晶化を防ぐために十分に速く冷却できるのかを探る、探偵物語です。

二つの対立する力

液体金属を、混み合ったダンスフロアだと考えてみてください。

1. 組織化への欲求（熱力学）: 金属が冷えるにつれて、原子はより安定した状態である「整然とした列（結晶化）」に並ぼうとします。温度が下がれば下がるほど、この衝動は強くなります。
2. エネルギーの欠如（動力学）: しかし、温度が下がると、原子は動きが鈍くなります。まるで、粘り気のあるシロップの中を動いているかのように、彼らの動きはどんどん遅くなります。彼らは、整然とした列にたどり着くための方法を見つけることができません。

「組織化したい」という欲求と、「動くのが遅すぎる」という事実の間の戦いが、その金属が結晶になるか、あるいはガラスになるかを決定します。

実験：デジタル・タイムマシン

研究者たちは、顕微鏡を使って個々の原子の動きをリアルタイムで観察することはできませんでした。なぜなら、それはあまりにも速い（10億分の1秒単位）からです。その代わりに、彼らは大規模なデジタル・シミュレーション（137万個の原子による「映画」）を構築し、何が起きているのかを観察しました。

また、超強力なX線レーザー（高速カメラのようなもの）を使用して、パラジウムの薄膜を照射し、溶かし、そして冷却される様子を観察するという、現実世界での実験も行いました。

彼らが発見したこと

1. 原子における「速度制限」
金属が冷却されるにつれて、原子は非常に予測可能な方法で、どんどんゆっくりと動いていくことがわかりました。これは、車が坂道を登る際に減速していくようなものです。坂が急であればあるほど（温度が低ければ低いほど）、車の速度は落ちます。彼らは、原子が一段階進むためにどれだけのエネルギーが必要かを正確に算出しました。

2. 「種」の問題（核生成）
結晶になるためには、液体の中に成長を開始するための「種（シード）」が必要です。

• 発見: この金属は、これらの種を作るのが驚くほど得意です。非常に低温であっても、小さな結晶の種が至る所で同時に自然発生します。
• 比喩: 人々がコンガ・ライン（一列になって踊る列）を作るのを阻止しようとする場面を想像してください。多くの物質では、阻止できるかもしれません。しかし、パラジウムの場合、音楽が止まった瞬間（冷却が始まった瞬間）、人々は即座に腕を組み始めます。研究者たちは、これらの種が形成される「完璧な」温度が、この金属の融点の約半分であることを突き止めました。この温度では、組織化への欲求が強く、かつ原子が結合するのに十分な速さで動いています。

3. 成長速度
一度種が形成されると、それは急速に成長します。

• 発見: 結晶の最前線は、秒速数メートルという速度で移動します。
• メカニズム: 研究者たちは、どのように成長するかについて二つの理論を検証しました。
  • 理論A（衝突制限）: 原子が結晶に衝突すると、フロントガラスに当たる雨のように、瞬時に付着する。
  • 理論B（拡散制限）: 原子が、混雑した劇場の席を探す人々のように、付着する場所を見つけるために液体の中をうごめき、移動しなければならない。
  • 結論: データは、理論Bが正しいことを示しました。原子は自分の居場所を見つけるために、うごめいて移動する必要があります。「衝突して付着する」という理論では、金属は実際よりも100倍速く成長すると予測されていました。

4. 「ガラス」という目標
この研究の究極の目標は、結晶化する前に、パラジウムをガラス化（ガラス転移）させるほど速く冷却できるかどうかを見ることでした。

• 結果: 結晶の形成を止めるには、毎秒10兆度（10¹³ K/s）という速度で金属を冷却する必要があります。
• 現実的な検証: 彼らが行った現実世界の実験では、毎秒5000億度（5×10¹¹ K/s）で冷却されました。
• 結論: 現実世界の冷却は遅すぎました。金属がガラスへと凍結する前に、結晶化を避けるための時間が足りなかったのです。結晶の「種」が形成され、成長する前に、金属は固まってしまいました。

総括

この論文は、純粋なパラジウムは、金属ガラスを作るという点において「質の悪い市民」であることを教えてくれます。それは結晶に戻ろうとする意欲が強すぎます。現在利用可能な最も速い冷却技術を用いたとしても、原子はあまりにも速く組織化されてしまいます。

研究者たちは、スーパーコンピュータによるシミュレーションを用いて、いつ、どこで結晶が形成され始めるかを正確に予測しましたが、その予測は現実世界のX線レーザー実験の結果と完璧に一致しました。これは、これらの薄膜において、結晶が汚れや容器の壁から始まるのではなく（不均質核生成）、ゼロから直接形成されている（均質核生成）ことを裏付けています。

要約すると: 純粋なパラジウムをガラスに変えることは容易ではありません。なぜなら、それはあまりにも組織化が得意だからです。それを実現するには、現在の彼らの実験が許容している自然の摂理よりも、もっと速く冷却する必要があるのです。

技術概要

技術要約：過冷却液体パラジウムの結晶化動力学

問題提起
過冷却液体（SCL）の結晶化動力学を理解することは、材料科学における根本的な課題である。特に純金属は、核生成および成長速度が極めて高いため、ガラス化が非常に困難である。液体の結晶化とガラス形成を決定づけるのは、熱力学的駆動力と原子移動度の間の競争であるが、純金属における急速な結晶成長を支配する具体的なメカニズムは未だ解明されていない。中心となる未解決の問いは、成長が原子拡散（熱活性化プロセス）によって制限されているのか、それとも原子衝突（音速による制約）によって制限されているのかという点である。さらに、実験設定において、不純物などの不均一なサイトが支配的になることが多いため、均質核生成と不均質な核生成を区別することは困難である。本研究では、多くの金属ガラスの主要成分である代表的な面心立方（FCC）金属である純パラジウム（Pd）に焦点を当て、これらの動力学を定量化し、超高速実験データとの検証を行う。

手法
著者らは、LAMPFコードを用い、Pd用に特別に開発された埋め込み原子法（EAM）ポテンシャルを使用した古典的分子動力学（MD）シミュレーションを採用した。統計的な堅牢性を確保し、有限サイズ効果を最小限に抑えるため、主要なシミュレーションには1,372,000個の原子からなる系を使用し、アンサンブル平均による時間－温度－変態（TTT）図の構築には、256,000個の原子からなる補完的な系を用いた。

シミュレーション・プロトコルは以下の通りである：

1. 熱履歴を消去するため、初期のFCC構成を2500 Kで溶融させる。
2. 0.38から0.65 $T_m$（700–1150 K）の温度範囲へ瞬時に急冷する。
3. OVITOソフトウェアを用い、多面体テンプレートマッチング（PTM）および粒界セグメンテーションアルゴリズムを使用して、微細構造の進化を監視し、秩序化した結晶粒と無秩序な液体を区別する。
4. 平均二乗変位（MSD）から自己拡散係数を算出し、結晶粒密度およびサイズの経時変化から核生成速度（$J$）および成長速度（$U$）を決定する。

これらのシミュレーションは、欧州XFELにおける光学的に溶融したPd薄膜から得られた、時間分解X線回折（XRD）データと直接比較された。実験セットアップは、フェムト秒レーザー励起とそれに続くX線プローブを用いており、約 $5 \times 10^{11}$ K/s の冷却速度を実現している。

主な結果

• 自己拡散: 調査した温度範囲における液相Pdの自己拡散係数（$D$）はアレニウス挙動に従い、活性化エネルギーは467(6) meV/atom、前指数因子は $5.8(4) \times 10^{-8}$ m²/s であった。液体－液体相転移を示すクロスオーバーは観察されなかった。
• 核生成動力学: 定常状態の均質核生成速度は、$0.5 T_m$ 付近で約 $4 \times 10^{35}$ m⁻³ s⁻¹ の極大値を示す。データは古典的核生成理論（CNT）の枠組み内でうまく適合し、液体－結晶界面エネルギー $\gamma \approx 0.17(1)$ J/m² を得た。
• 成長メカニズム: 結晶成長速度は、秒速数メートル（1–8 m/s）のオーダーであった。温度依存性と速度の大きさは、拡散制限型のウィルソン–フレンケルモデルと密接に一致した。対照的に、衝突制限型の成長モデルは速度を2桁近く過大評価しており、成長は衝突制限ではなく熱活性化プロセスであることを示している。
• TTT図および臨界冷却速度: シミュレーションから構築されたTTT図は、$0.5 T_m$ 付近の約100 psに「ノーズ」を示している。これは、$\sim 10^{13}$ K/s のガラス化のための臨界冷却速度に対応する。
• 実験的検証: $5 \times 10^{11}$ K/s の冷却速度におけるMD予測の結晶化開始時刻および温度は、実験によるXRDの観測結果と一致した。実験的に観察された過冷却は、およそ $0.65 T_m$ に達しており、これはシミュレーションの予測と一致している。

意義および主張
本論文は、急速に急冷された純Pd薄膜において、達成可能な過冷却は不均質な核生成ではなく、均質な核生成によって支配されることを確立した。MDシミュレーションとXFELポンプ・プローブ実験との一致は、大規模なMDが超高速結晶化動力学の研究に有効であることを検証している。

本研究は以下の結論を得ている：

1. 純PdのようなFCC金属は、その臨界冷却速度（$\sim 10^{13}$ K/s）が、研究された薄膜実験で達成可能な速度（$5 \times 10^{11}$ K/s）を大きく上回るため、ガラス化が極めて困難である。
2. 過冷却Pdにおける結晶化前線の伝播は拡散制限型であり、衝突制限領域よりもウィルソン–フレンケル機構を支持している。
3. これらの結果は、実験データに対する一貫した物理的枠組みを提供し、観察された限界が実験的なアーティファクトや不十分な急冷技術のみによるものではなく、材料固有の均質核生成特性に起因することを裏付けている。