Capturing exchange-correlation spin-torque effects with a semilocal functional

技術要約：半局所的交換相関スピントルク効果の捕捉

問題提起
磁性に関する第一原理計算で使用される半局所的な交換相関（XC）汎関数は、根本的な理論的欠陥を抱えている。それは、XCスピントルクに対する無知である。局所スピン密度近似（LSDA）は、XC磁場が磁化と整列するという条件（ゼロ・スピントルクを生む）を自然に満たすが、この性質はGGA（一般化勾配近似）やメタGGAなどのLSDAを超える汎関数には成立しない。従来の非共線的密度汎関数理論（DFT）の実装では、これらの半局所的汎関数を非共線的な枠組みに強制的に適合させるため、XC磁場を局所的な磁化に整列させる回転操作が行われることが多い。この「アドホック」なゼロトルク近似は、ゼロトルク定理（XC磁場は全域的な正味のトルクは及ぼさないが、局所的なトルクは許容する）と理論的に矛盾しており、スピン電流密度を無視している。この欠陥は、非共線磁性、スピンダイナミクス、および軌道モーメントを持つ系において特に重要であり、電子収束の悪化や全エネルギーの不正確さを招く可能性がある。

手法
著者らは、Vienna ab-initio Simulation Package (VASP) において、スピン電流密度汎関数理論（SCDFT）の枠組み内で半局所的XC汎関数を実装することで、この問題に対処している。この実装では、射影原子軌道法（PAW）を利用し、$2 \times 2$ スピノル形式を導入している。

主要な手法構成要素：

• 変分原理： 全エネルギーの最小化は、スピン密度（$n_{\sigma\sigma'}$）および運動エネルギー密度（$\tau_{\sigma\sigma'}$）に加えて、$2 \times 2$ パラマグネティック・スピン電流密度（$j_{\sigma\sigma'}$）に関して拡張されている。
• XCポテンシャル： 汎関数微分により、$2 \times 2$ のXCポテンシャル演算子 $\hat{v}_{xc}$ が得られる。これには、局所成分、運動エネルギー密度に依存する項、およびスピン電流密度から導かれる非アーベル型XCベクトルポテンシャル（$A_{xc}$）が含まれる。
• 汎関数の形式： 提示された特定の汎関数は、Becke-Roussel (1989) の交換およびColle-Salvettiの相関に基づいている。
  • 交換： ラプラシアンを含まない非共線的Becke-Roussel（LFNCBR89）汎関数が構築されている。$U(1) \times SU(2)$ ゲージ不変性を確保するため、交換ホール曲率は標準的な運動エネルギーではなく、スピン電流項を含む運動エネルギー密度（$\bar{\tau}$）を用いて定義されており、これによりゲージ対称性の破れを回避している。
  • 相関： 非共線的Colle-Salvetti（NCCS）汎関数が採用されており、ゲージ不変性を維持するためにスピン電流密度を組み込んだ多体波動関数アンサンブルを利用している。
• 実装の詳細： コードは、PAW形式内で$2 \times 2$ パラマグネティック・スピン電流密度を計算し、平面波グリッドと全電子の一中心項の両方からの原子価寄与を考慮している。また、力および応力計算のために、スピン電流密度（スピノル空間におけるランク2のテンソル）の変換特性を扱うための特別な処理が行われている。

主要な貢献

1. 初のPAW実装： 本研究は、変分問題にスピン電流密度を明示的に含み、$2 \times 2$ のXCポテンシャルをもたらす、静的な非共線的SCDFT汎関数の初のフルPAW実装を報告している。
2. ゲージ不変性： 著者らは、$U(1) \times SU(2)$ ゲージ不変性を保証する導出を提供している。これは、非共線的汎関数の定式化においてしばしば見落とされる重要な要件である。
3. 新しい汎関数ファミリー： 本論文は、$2 \times 2$ スピン密度、運動エネルギー密度、およびスピン電流密度に依存する特定の「フレーバー」のSCDFT汎関数（LFNCBR89-NCCS）を提示しており、既存の他の汎関数を非共線的領域へと一般化するための道筋を示している。

結果
実装は、孤立した $Cr_3$ 分子（非共線的）と、バルクMnO（共線的反強磁性）の2つの系でテストされた。

• $Cr_3$ 分子：

  • SCDFT-LFNCBR89-NCCS 汎関数は、より計算コストの高い厳密交換（EXX）法（具体的にはKrieger-Li-Iafrate近似）によって得られるものと定性的に類似したXCスピントルクパターンを生成する。
  • スピントルクの大きさは、EXX-KLI 計算で見られるものの約90%である。
  • オンサイト磁気モーメント（$\mu \approx 3.11 \mu_B$）は、メタGGA（R2SCAN）やハイブリッド汎関数（HSE06）と同等の値であり、LSDAやGGA-PBEよりも大幅に高い。
  • ゼロ局所スピントルクを強制する（磁化をXC磁場に整列させるように制約する）と、フルSCDFT計算と比較して全エネルギー差は約0.04 eVとなり、この小さな系においてもスピントルクがエネルギー的に重要であることが浮き彫りになった。

• バルクMnO：

  • スピン軌道相互作用のないこの共線的系においては、XCスピントルクは極めて小さく（$Cr_3$ の5桁低い）、ゼロトルク近似がこのような系においては妥当な近似であることを裏付けている。
  • トルクは小さいものの、スピントルク効果の導入は磁気構成の相対的な安定性に変化を与える。SCDFT計算と従来の非共線的計算（共線的汎関数を使用）との全エネルギー差は、磁気構成や格子定数に応じて meV/f.u. のオーダーである。
  • 収束性： 特筆すべき知見として、SCDFTのスピントルクアルゴリズムは、汎関数の複雑さが増しているにもかかわらず、従来の非共線的アルゴリズムと比較して、電子最小化ステップ数が少なく（約3分の1）収束することである。
  • 構造パラメータ（格子定数、体積弾性率）および磁気モーメントは、アルゴリズムの選択（共線的 vs 非共線的 vs SCDFT）にはほとんど影響されないが、XC汎関数の選択には大きく依存する。

意義と主張
本論文は、SCDFT内でXC効果を厳密に扱うことにより、半局所的汎関数におけるXCスピントルクの欠如を解消することを主張している。著者らは以下の点を強調している：

• 共線的汎関数に基づく従来の非共線的計算は、不当なゼロトルク近似を導入しており、それが特に非共線的秩序を持つ系において、電子収束の悪化や不正確な全エネルギーを招く可能性がある。
• 新しいSCDFTフレームワークは、高レベルのEXX手法と同等の計算コストで、定性的に類似したパターンを回収しつつ、一貫したゲージ不変な方法でXCスピントルクを計算できる。
• 提示された特定の汎関数（LFNCBR89-NCCS）には、MnOの構造パラメータの記述における限界があるものの、フレームワーク自体が主要な貢献である。これは、アドホックな近似に頼ることなく、既存の確立されたXC汎関数（SCANなど）を非共線的領域へと拡張するための道を開くものである。
• 本研究は、有限の軌道角運動量または非共線磁性を持つ系において、XCエネルギーの風景を物理的に正しい形で記述するためには、スピン電流密度の包含が不可欠であることを示している。