On the limits of the energetic coupling between field dislocation mechanics and phase field crystal

技術要約：電場転位力学（FDM）とフェーズフィールド結晶（PFC）間のエネルギー結合の限界

問題提起
本論文は、転位を介した塑性に関する二つの異なるモデリング・フレームワーク、すなわち電場転位力学（FDM）とフェーズフィールド結晶（PFC）モデルを統合しようとする理論的および計算的な課題に取り組んでいる。

• FDMは、転位場を持つ連続体の初期値問題（IBVP）を解くことに成功しており、正しい弾性応力の平衡と輸送ダイナミクスを捉えることができる。しかし、固有の結晶構造を欠いており、転位コアの局在化は、経路依存的な塑性歪み変数に依存する非凸エネルギー項によって規定されなければならず、これは熱力学的な不整合をもたらす。
• PFCは、基礎となる原子スケールの結晶構造を捉え、非凸なスウィフト・ホーヘンバーグ自由エネルギー汎関数によってコンパクトな転位コアを自然に維持する。しかし、その支配方程式は散逸的な考慮のみから導出されているため、弾性波の伝播ではなく、非物理的な拡散的機械的緩和を引き起こし、正しい弾性挙動を捉えることができない。

AcharyaとViñals（2020）によって提案された結合フレームワーク（PFCFDM）は、エネルギーペナルティ項を導入することで、これらのモデルを統一しようと試みた。この項は、FDMからの弾性歪み（$U^e$）と、PFCの秩序パラメータから導かれる構成歪み（$Q$）の$L^2$距離にペナルティを課すものである。その目的は、経路依存的な変数を自由エネルギーに導入することなく、PFCに正しい弾性を与え、FDMに固有の結晶構造を与えることにある。本論文は、この結合が意図した成果を達成しているかどうかを調査するものである。

手法
著者らは、PFCFDMの結合を解体するために、変分解析と数値シミュレーションを組み合わせて用いている。

1. 変分解析： 著者らは、秩序パラメータ $\psi$ および弾性歪み $U^e$ に関する結合自由エネルギーの変分微分を導出する。そして、結果として得られるフェーズフィールド発展方程式における強制項と、修正された応力テンソルの構造を分析する。
2. 数学的分解： 著者らは、ストークス・ヘルムホルツ分解を用いて、歪み場を適合的な部分（回転ゼロ）と不適合な部分（発散ゼロ）に分離する。そして、ペナルティ項がこれらの成分とどのように相互作用するかを検討し、特に、転位のトポロジー（バーガース・ベクトル）を担う不適合成分に焦点を当てる。
3. 数値シミュレーション： 2次元六方格子に対してスペクトル法を用いて結合モデルを実装する。シミュレーションは以下の3つのシナリオを対象とする：
   • 静止したエッジ転位双極子による、$U^e$ と $Q$ の間の平衡状態の検討。
   • 外部荷重がない状態でのコアの進化による、コアの広がりを防ぐ能力のテスト。
   • マクロな剪断応力下での転位の滑りと消滅。ここでは、FDM駆動型（転位密度の輸送方程式を解く）と、PFC駆動型（塑性歪み速度をPFCの転位フラックスに等置する）という二つの駆動パラダイムを比較する。

主要な貢献と理論的知見
本論文は、$L^2$ エネルギー結合が根本的な構造的限界を抱えていることを実証している：

• 不適合歪みに対する無感度： 秩序パラメータ $\psi$ に関するペナルティ項の変分微分は、ミスマッチのダイバージェンス（$\nabla \cdot (Q - U^e)$）に依存する。その結果、結合力は歪み場の適合的な部分にのみ作用する。この結合は、転位のトポロジー（歪みの回転）に関するすべての情報を含む不適合な（ダイバージェンス・フリーな）部分に対しては盲目である。したがって、この結合は、極転位密度（$\alpha = \nabla \times U^e$ および $\tilde{\alpha} = \nabla \times Q$）を一致させることはできない。
• 力学の拡散的伝達： 結合は、弾性的緩和ではなく、フェーズ拡散を通じてFDMドメインからPFCドメインへと機械的な境界条件を伝達する。ペナルティ項は複素振幅の位相に対する拡散演算子として作用する。転位コア付近では、振幅が消失するため、この拡散係数は発散し、非物理的な挙動を引き起こし、関連する時間スケールでの弾性情報の伝達を妨げる。
• コアの広がりとピンニング： この結合は、FDMにおける転位コアの不自然な広がりを防ぐことができない。さらに、ペナルティ項はPFCの転位を固定する「復元スプリング」として作用する。著者らは、結合パラメータ $c_{pen}$ における矛盾を特定している。すなわち、$c_{pen}$ はスウィフト・ホーヘンバーグの緩和を克服できるほど大きく（$c_{pen} \gg 1$）なければならない一方で、転位をロックしない程度に小さくなければならない（$c_{pen} \ll 1$）。これにより、モデルが高速なリラクサーとしても、移動する転位の駆動源としても機能しない「空の動作ウィンドウ」が生じている。
• PFC駆動型アプローチにおける熱力学的不整合： FDM駆動型アプローチが二つの異なる解決不能な転位密度に苦しむ一方で、PFC駆動型アプローチ（$\dot{U}^p = J^\psi$ を強制する）は、塑性散逸が非負であることを保証できない。数値結果によれば、塑性散逸が負になることがあり、これはPFCによる運動が純粋に応力駆動ではなく、特定の領域において熱力学的整合性を損なっていることを示している。

数値結果

• 平衡： 静的なシミュレーションにおいて、システムは $U^e$ と $Q$ のダイバージェンスが一致する（$\nabla \cdot U^e = \nabla \cdot Q$）局所的な最小状態に達するが、不適合な部分（したがって転位コア）は依然として不一致のままである。
• コアの進化： 外部荷重がない場合、FDM密度 $\alpha$ は結合強度に関わらず拡散的に広がる。PFC密度 $\tilde{\alpha}$ はコンパクトなままであるが、これら二つの密度は独立して進化する（$\alpha \neq \tilde{\alpha}$）。
• 消滅： 負荷が印加されると、転位の運動は著しく阻害される。FDM駆動型の場合、コアは広がり、効率的な消滅に失敗する。PFC駆動型の場合、運動はペナルティによる「スプリング」によってロックされるか、あるいは純粋な拡散的PFC進化よりも遅い速度で発生し、弾性緩和の回復は見られない。

意義と主張
本論文は、$L^2$ エネルギー結合が、結晶学と転位力学を統合することには成功していないと結論付けている。

• 結合項はダイバージェンス駆動の強制力として作用し、歪み場の適合的な部分のみを一致させるが、トポロジカルな情報（不適合な部分）は結合されない。
• 機械的な情報は弾性的にではなく拡散的に伝達されるため、非物理的なコア挙動が生じ、弾性平衡を満たすことができない。
• モデルは、コアの局在化と転位の移動能を同時に実現することを妨げる、結合パラメータにおける根本的なトレードオフに直面している。
• FDM駆動型とPFC駆動型の両方の実装において決定的な欠陥が存在する。前者は転位密度の不一致に、後者は熱力学第二法則の違反（負の塑性散逸）に直面している。

著者らは、この結合の定式化自体は熱力学的に一貫しているものの、その変分構造が、FDMの連続体力学とPFCの結晶学的解像度を調和させる能力を本質的に制限していると主張している。本研究は、この特定の結合の適用領域を定義しており、現在の形態では、二つの理論間の堅牢な架け橋として機能し得ないことを事実上示している。