<100>引張荷重下タンタル単結晶の弾性挙動：分子動力学による温度、圧力、相転移知見

要旨

原子論的シミュレーションはタンタル（Ta）の力学応答を支配するメカニズムを解明できます。本研究では分子動力学（MD）を用い、温度と圧力が<100>引張荷重下Ta単結晶の弾性特性に及ぼす影響を検証。まず、Zhou-EAMとRavelo-EAMの2つの埋め込み原子法（EAM）ポテンシャルを、広い圧力範囲でのTa弾性定数再現能力で比較。Ravelo-EAMポテンシャルが優れ、140 GPaまで弾性率を精度よく再現。このポテンシャルベースのMDシミュレーションで、1 Kにおいてモデルサイズやひずみ速度に依存せず、破壊に先立って体心立方（BCC）→面心立方（FCC）相転移が生じることを発見。1-1500 Kで<100>弾性率E₁₀₀は温度とともに線形減少、降伏応力は2次関数的温度依存性を示す。圧力依存シミュレーション（0-140 GPa）でE₁₀₀は単調増加。これらの結果は、極限条件下でのTa弾性挙動を予測する包括的原子論的フレームワークを提供します。

背景

タンタルは常温常圧で体心立方（BCC）格子をとり、卓越した相安定性、高融点（3269 K）、酸化耐性で知られます。これら特性によりTaはマイクロエレクトロニクス拡散バリア、耐摩耗コーティング、高温超合金の理想的候補。ダイヤモンドアンビルセル、密度汎関数理論（DFT）、MDを用いた膨大な実験・理論研究が高圧・高温下Ta挙動を探求し、135 GPaまでBCC相持続と圧力に伴う弾性率増大を一貫して確認。これら進歩にもかかわらず、引張荷重下Ta動的応答の体系的原子論的研究は不足。本研究は、<100>引張時のTa弾性特性に対するサイズ、ひずみ速度、温度、圧力効果を評価し、引張誘起相転移の可能性を探ることでこのギャップを埋めます。

方法と実験詳細

物理モデリング

BCC単位格子を<001>、<010>、<100>方向に複製し、格子定数a=b=c=3.301 Åで立方Taセルを構築。4サイズ（12a=3.96 nm、18a=5.94 nm、24a=7.92 nm、30a=9.90 nm）で3,456-54,000原子を含む。図1は主要構造として用いた3.96 nmモデルを示す。

MDシミュレーションプロトコル

全シミュレーションをLAMMPSで実施。Zhou-EAMとRavelo-EAMの2EAMポテンシャルを0 Kおよび静水圧下での弾性定数再現能力で評価。Ravelo-EAMポテンシャル（Raveloら、Phys. Rev. B 2013, 88, 134101開発）が対象圧力範囲でZhou-EAMを上回る性能。直接法でC₁₁とC₁₂を計算し、式(1)で<100>弾性率を導出。後続シミュレーションで1 fsタイムステップ、全方向周期境界条件、50 ps NVT緩和後NPT平衡で静水圧（0-140 GPa）印加。x方向引張荷重をひずみ速度5×10⁸ s⁻¹から7.5×10⁹ s⁻¹で15%伸びまで印加。

結果と考察

伸長プロセス

図3に零圧力下<100>引張時の応力-ひずみ応答と対応する原子配置を示す。弾性変形は~7%ひずみまで線形、以降~7.4%でBCC→FCC変態開始、~9.8%で完了。FCC相は~13%で表面破壊まで持続。変態と破壊メカニズムは全モデルサイズとひずみ速度で一貫し、弾性率と降伏応力へのサイズ・速度効果は無視可能。

サイズとひずみ速度効果

表2でモデルサイズ（3.96-9.90 nm）が弾性率（~139 GPa）と降伏応力に影響しないことを実証。同様に表3でひずみ速度5×10⁸-7.5×10⁹ s⁻¹間でこれら特性に識別可能な変化なし。したがって3.96 nmセルと5×10⁸ s⁻¹ひずみ速度を後続全解析に採用。

温度依存性

図4に1 Kから1500 Kまでの応力-ひずみ曲線進化を示す。弾性率E₁₀₀と降伏応力とも温度とともに単調減少。弾性率は約42%減少（136.5→76.7 GPa）、降伏応力は約51%減少（≈8→4 GPa）。線形フィットでE₁₀₀ = 138.07 - 0.04094 T (GPa)（1500 Kまで有効）、降伏応力は2次フィットで捕捉。

圧力依存性

図6に複数温度での静水圧（0-140 GPa）に対するE₁₀₀を示す。1 Kで動的結果は式(1)ベース静的計算とよく一致。≤600 Kで20-60 GPa間に凹領域出現、先行MD研究と整合。高温では単調。全温度で~40 GPa以上でE₁₀₀増加、Ta剛性が圧縮下で増強されることを確認。

結論

Ta単結晶の<100>引張について、広範囲温度（1-1500 K）と圧力（0-140 GPa）で包括的MD調査を実施。Ravelo-EAMポテンシャルが高圧シミュレーションに信頼性高いことを実証。主な知見：(1) 低温で破壊に先立つBCC→FCC相転移；(2) 弾性率と降伏応力が温度に対しそれぞれ線形・2次関数的に減少；(3) ~40 GPa以上で圧力が弾性率を増強；(4) モデルサイズとひずみ速度は測定弾性特性に無視可能な影響。