Exchange striction determines how fast antiferromagnetic insulators demagnetize

技術要約：反強磁性絶縁体における交換厳密性が脱磁速度を決定する

問題提起
反強磁性絶縁体における超高速脱磁は、化合物によってピコ秒からナノ秒まで極めて広いタイムスケールの広がりを示す。金属強磁性体における速度限界はエリオット・ヤフェット型のスピン反転散乱によって理解されているが、絶縁体反強磁性体における速度限界を支配するメカニズムは未解決のままである。正味の角運動量を格子に転送する必要のない補償反強磁性体において、脱磁速度は理論的にはレーザー加熱された格子からスピン系へのエネルギー流によって制限される。しかし、この広がりの大きさや、新しい化合物が具体的にどの程度の速度で脱磁するかを予測または説明できる材料パラメータは特定されていない。

手法
著者らは、酸化クロム（$\text{Cr}_2\text{O}_3$）を用いた実験的手法と、構造的に類似した $\bar{R}3c$ 空間群を持つホウ酸鉄（$\text{FeBO}_3$）との比較による理論的手法を組み合わせ、超高速脱磁を調査した。

1. 実験的アプローチ:

   • 手法: 時間分解第二高調波発生（SHG）を用い、磁気秩序のダイナミクスをプローブした。円偏光プローブパルス（1200 nm）は、ネールベクトルの配向に敏感なSHG信号（600 nm）を生成する。
   • プロトコル: 400 nmのポンプパルスが格子系を加熱し、SHGコントラストを用いて秩序パラメータ $|L|$ の減衰を追跡した。実験は、ネール温度（$T_N \approx 307.3$ K）の上下両方の最終温度に達するように、様々なポンプフルエンスで行われた。
   • 指標: 特徴的な脱磁時間 $\tau$ は、SHGコントラストが全減少量の $1/e$ まで減衰した遅延として定義された。

2. 理論的アプローチ:

   • 原子論的スピンダイナミクス（ASD）: 現象論的なギルバート減衰を用いて、スピンの急激な格子温度変化への応答をモデル化した。
   • 第一原理計算: 密度汎関数理論（DFT+U）を用い、ハイゼンベルク交換相互作用定数（$J_{ij}$）とそのイオン変位に対する微分（$\partial J/\partial u$）、すなわち交換厳密（exchange striction）を計算した。
   • 結合スピン・格子ダイナミクス（SLD）: 全ての設定パラメータなしに、DFT+Uから導出された交換テンソルと力定数を用い、スピンに対するランダウ＝リフシッツ＝ギルバート方程式とイオンに対するニュートンの運動方程式を統合した非摂動的シミュレーションを行った。
   • 運動学的解析: 一フォノンから二マグノンへの崩壊過程の位相空間を計算し、両材料におけるフォノン分散と二マグノン連続体の重なりを比較した。

主な結果

1. 実験的な格差: $\text{Cr}_2\text{O}_3$ は $T_N$ 以上に駆動された際、2 ps 以内に脱磁する。これは $\text{FeBO}_3$ よりも2桁速い速度である。$\text{FeBO}_3$ は $T_N$ 付近で臨界減速を示すが、$\text{Cr}_2\text{O}_3$ はこの領域を回避して、極めて高速な固有のエネルギー転送率を明らかにする。
2. 現象論的モデルの限界: 一定の減衰を用いる標準的なASDシミュレーションは、$T_N$ 付近の臨界減速は捉えられたものの、$T_N$ 以上で見られる $\tau$ の急激な低下を再現できなかった。これは、温度依存的なスピン・格子結合メカニズムが欠落していることを示唆していた。
3. 交換厳密の役割: 第一原理計算により、交換厳密（交換結合の原子変位依存性）が支配的な要因であることが特定された。
   • 構造的起源: $\text{Cr}_2\text{O}_3$ は、磁気相互作用が直接的な $t_{2g}$-$t_{2g}$ 重なりによって支配される短い Cr-Cr 接点（2.65 Å および 2.89 Å）を有している。これにより、交換厳密定数（$|\partial J/\partial u|$）は $\sim 100$ meV/Å に達する。対照的に、$\text{FeBO}_3$ は酸素を介した超交換相互作用に依存しており、短い陽イオン間の接触を持たないため、$|\partial J/\partial u| \lesssim 10$ meV/Å である。
   • 位相空間: $\text{Cr}_2\text{O}_3$ の強い交換相互作用は広いマグノンバンドを形成し、二マグノン連続体がフォノンスペクトル全体（最大25 THzまで）とほぼ完全に重なる。$\text{FeBO}_3$ では、連続体は 12 THz 以下に限定されており、ほとんどのフォノンに崩壊チャネルが存在しない。
4. 定量的合致: 位相空間因子（$\text{Cr}_2\text{O}_3$ の方が約3倍大きい）と行列要素因子（厳密性により約$10^2$倍大きい）の積は、約300の速度比を予測する。これは、実験的に観察された約100という比率と一致する。
5. パラメータフリーのシミュレーション: アブイニシオ（第一原理）パラメータのみを用いた結合SLDシミュレーションは、測定された比のオーダー（$\tau_{\text{Cr}_2\text{O}_3} \approx 10$ ps vs. $\tau_{\text{FeBO}_3} \approx 500$ ps）を再現し、交換積分（$J$）の変位依存性が主要な駆動力であることを裏付けた。

意義および主張
本論文は、絶縁体反強磁性の脱磁速度は経験的な謎ではなく、以下の3つの具体的な要素によって計算可能な材料特性であることを確立した：

1. 交換厳密テンソル（$\partial J/\partial u$）
2. 熱的なイオン変位の大きさ
3. 二マグノン連続体とフォノンスペクトルの重なり

著者らは、この枠組みによって、合成前に候補材料の脱磁速度を予測できると主張している。具体的には、短距離の陽イオン間結合を持つ材料（直接交換と厳密性を高める）はより速く脱磁し、酸素を介した超交換相互作用のみに依存する材料（方解石構造の $\text{MnCO}_3$ など）はより遅くなる。これは、熱アシスト型反強磁性メモリのための材料スクリーニングへの実用的な経路を提供し、絶縁体磁性系における超高速制御の理解を前進させるものである。