Magnetic anisotropy and intermediate valence in CeCo$_5$ ferromagnet

CeCo$_5$の中間価数状態と磁性を正確に記述するため、DFT+$U$にアンダーソン不純物モデルの厳密対角化を組み合わせる手法を提案し、実験値と一致する磁気異方性エネルギーや全磁気モーメントを再現することで、高性能希土類低含有永久磁石の開発指針を示しました。

要約

🧲 物語の舞台：「セリウム磁石」の謎

昔から、強力な磁石（永久磁石）を作るには、ネオジムやジスプロシウムといった「高価で貴重なレアアース」が必要でした。しかし、これらは採掘が難しく、価格も高騰しています。

そこで科学者たちは、**「セリウム（Ce）」**という、より安く手に入りやすい元素に注目しました。セリウムは地殻に豊富にあるからです。

しかし、問題が一つありました。
**「セリウムを混ぜると、なぜか磁石の力が弱まってしまう」のです。
なぜか？それはセリウムという元素が、「二面性」**を持っているからです。

🎭 主人公の正体：「二面性」を持つセリウム

セリウムは、磁石の中で**「3 価（Ce3+）」と「4 価（Ce4+）」という、2 つの異なる姿を瞬時に切り替えています。
これを「中間価数（Intermediate Valence）」**と呼びます。

• 3 価のセリウム：しっかりとした「磁石の力」を持っている。
• 4 価のセリウム：磁石の力がほとんどない。

セリウムは、この 2 つの姿を**「激しく揺れ動いている（フラクチュエーション）」状態にあります。
まるで、「静かに座っている人（3 価）」と「走り回っている人（4 価）」が、一瞬ごとに入れ替わっているような状態**です。

この「揺れ動き」が、従来の計算方法（DFT など）では捉えきれませんでした。
従来の計算は、**「セリウムは常に 3 価だ（または常に 4 価だ）」**と決めつけて計算してしまうため、実際の「揺れ動く姿」を無視してしまい、磁石の強さや方向性を間違って予測してしまっていたのです。

🔍 科学者の新しい道具：「超高速カメラ」

この論文の著者たちは、新しい計算手法**「DFT+U(ED)」という「超高速カメラ」**のようなツールを使いました。

• 従来の方法：セリウムの姿を「平均化」して見てしまうので、揺れ動きの細かさが消えてしまい、磁石の性質を正しく見られない。
• 新しい方法：セリウムの「3 価」と「4 価」の瞬間的な入れ替わり（揺れ動き）を、**量子力学の Exact Diagonalization（厳密対角化）**という技術で詳細にシミュレーションする。

これにより、セリウムが**「実は磁石の力を半分くらいしか出せていない（揺れ動いているため）」**という真実を突き止めました。

🎯 発見された 2 つの重要な事実

この新しい計算でわかったことは、以下の 2 点です。

1. 磁石の強さ（全磁気モーメント）の正解

実験で測られた磁石の強さは、**「6.70」**という値でした。

• 従来の計算：6.63〜7.23 くらい（少しズレていた）。
• 新しい計算：6.70（バッチリ一致！）

セリウムの「揺れ動き」を考慮することで、磁石の強さを正確に予測できるようになりました。

2. 磁石の「向き」を決める力（磁気異方性）の正解

磁石には「特定の方向にしか磁気が出ない」という性質（異方性）があります。これが強ければ強いほど、高性能な磁石になります。
実験では、この力が**「5.5」**でした。

• 従来の計算：2.0 くらい（弱すぎる）。
• 古い考え方の計算：7.5 くらい（強すぎる）。
• 新しい計算（セリウムだけ）：2.95（まだ弱い）。
• 新しい計算（セリウム＋コバルトの修正）：4.8（実験値 5.5 に非常に近い！）

ここが最大のポイントです。
セリウムだけでなく、磁石の骨格を作る**コバルト（Co）**の電子の動きも、セリウムの「揺れ動き」と連動して影響を受けていることがわかりました。
**「セリウムの揺れ動き」＋「コバルトの電子の相互作用」**を両方考慮することで、初めて実験に近い「強力な磁石の向き」を計算で再現できました。

💡 この研究が意味すること

この研究は、単に「CeCo5 という物質の計算ができた」というだけでなく、**「高価なレアアースを使わずに、安価なセリウムで高性能な磁石を作るための設計図」**を描き出しました。

• これまでの壁：「セリウムは磁石に向かない」と思われていた。
• 今回の突破：「セリウムの『揺れ動き』を正しく理解し、制御できれば、高性能磁石が作れる」と証明した。

🌟 まとめ

この論文は、**「セリウムという元素が、実は『二面性』を持って揺れ動いている」**という複雑な性質を、最新の計算技術で捉え直しました。

その結果、**「磁石の強さ」と「磁石の向き」**という、磁石の性能を決める 2 つの重要な要素を、実験結果とほぼ完璧に一致させることに成功しました。

これは、「高価なレアアースに頼らず、安価なセリウムで、未来の高性能磁石を作る」という夢への、大きな一歩となりました。まるで、「揺れ動くシャボン玉」の性質を理解することで、それを丈夫なボールに変える魔法を手にしたようなものです。

技術概要

以下は、Alexander B. Shick と Evgenia A. Tereshina-Chitrova による論文「Magnetic anisotropy and intermediate valence in CeCo5 ferromagnet」の技術的サマリーです。

1. 研究の背景と課題 (Problem)

• 希土類 - 遷移金属磁石の重要性: Nd, Sm, Dy などの希土類元素を含む金属間化合物は強力な永久磁石として重要ですが、資源的な制約から、より豊富で安価なセリウム (Ce) を用いた磁石の開発が急務となっています。
• CeCo5 の特異性: CeCo5 は六方晶構造を持ち、c 軸方向に容易軸を持つ強磁性体ですが、隣接する LaCo5 や PrCo5 と比較して、格子定数のランタニド収縮の傾向に従わず、磁気モーメントが異常に低く、キュリー温度も低いという特異な挙動を示します。これは Ce の混合価数（Ce3+ と Ce4+ の間の揺らぎ）に起因すると考えられています。
• 既存理論の限界:
  • 標準 DFT (LSDA/GGA): Ce の 4f 電子の強い電子間相互作用を正しく扱えず、Ce のスピン・軌道モーメントを過大評価し、磁気異方性エネルギー (MAE) を実験値（約 5.5 meV/f.u.）より大幅に過小評価（約 2.0 meV/f.u.）する傾向があります。
  • 静的 DFT+U: 動的な相関効果を無視しているため、価数揺らぎを捉えきれず、MAE の予測値が実験値と一致しないか、あるいは過大評価（HIA 近似など）するなどの不一致が生じています。
  • DFT+DMFT: 価数揺らぎやスピン・軌道モーメントの低減を成功裏に説明しましたが、計算コストが高く、MAE のような微小なエネルギー差を高精度で計算する技術的課題が残っていました。

2. 手法 (Methodology)

本研究では、CeCo5 の電子構造、磁気特性、および磁気異方性を解析するために、DFT+U(ED) 手法を採用しました。これは以下の要素を組み合わせたハイブリッド手法です。

• DFT+U とアンダーソン不純物モデル (AIM): Ce の 4f 殻に対して、密度汎関数理論 (DFT) に加え、多軌道アンダーソン不純物モデルのハミルトニアンを厳密対角化 (Exact Diagonalization: ED) して処理します。
• LDA++ 手法の適用:
  • s, p, d 殻は通常の DFT (LSDA) で記述。
  • Ce の 4f 殻は、クーロン相互作用、強いスピン軌道結合 (SOC)、結晶場分裂、および宿主バンドとのハイブリダイゼーションを明示的に含む AIM を用いて記述。
• 自己無撞着計算:
  • ED によって得られた自己エネルギー $\Sigma(z)$ を用いて、局所グリーン関数を計算し、有効ポテンシャル $V_U$ を構築します。
  • このポテンシャルを Kohn-Sham 方程式に挿入し、電荷密度の自己無撞着解を得るまで反復計算を行います。
• MAE の計算: [110] 面内方向と [001] 軸方向の磁化に対する全エネルギー差 ($E_{[110]} - E_{[001]}$) を、高精度な k 点メッシュ（3825 点）を用いて計算します。
• Co 3d 殻への相関の考慮: Ce だけでなく、Co の 3d 殻にも Coulomb 相互作用 (U) と交換相互作用 (J) を導入した DFT+U 計算を行い、二サブラティスモデルを用いて MAE を補正しました。

3. 主要な成果と結果 (Key Contributions & Results)

A. 磁気モーメントと価数揺らぎ

• Ce のスピン・軌道モーメントの低減: DFT+U(ED) により、Ce のスピンモーメント ($M_S^{4f} = -0.14 \mu_B$) と軌道モーメント ($M_L^{4f} = 0.10 \mu_B$) が、従来の DFT や静的 DFT+U に比べて著しく低減することが確認されました。これは DFT+DMFT の結果と定量的に一致します。
• 中間価数状態の再現: 計算された 4f 電子占有数 ($n_f$) は 0.64 となり、Ce3+ ($n_f \approx 1$) と Ce4+ ($n_f \approx 0$) の間の中間価数状態を示しています。基底状態は主に $f^0$ (確率 0.39) と $f^1$ (確率 0.57) の重ね合わせであり、Ce4+-Ce3+ の価数揺らぎが動的にスクリーニングされていることが示されました。
• 全磁気モーメント: 計算された単位式あたりの全磁気モーメントは 6.70 $\mu_B$ であり、実験値の範囲 (6.5 - 7.1 $\mu_B$) とよく一致します。

B. 電子状態密度 (DOS) と分光データとの一致

• 光電子分光 (PES) と Bremsstrahlung Isochromat Spectroscopy (BIS): DFT+U(ED) によって計算された 4f 状態密度は、実験的な PES および BIS スペクトルを非常に良く再現します。
  • 特に、フェルミ準位直下の占有状態や、フェルミ準位より 0.5 eV, 0.9 eV 上に現れる非占有状態のピーク構造が、静的な DFT+U(HIA) 近似では再現できず、ED によるハイブリダイゼーション効果の取り込みが不可欠であることを示しました。

C. 磁気異方性エネルギー (MAE)

• Ce 4f 殻のみの寄与: DFT+U(ED) 単独での MAE 計算値は 2.95 meV/f.u. でした。これは従来の DFT (約 2.0 meV) よりも大幅に向上していますが、実験値 (5.5 meV) よりはまだ低いです。
• Co 3d 殻の相関効果の重要性: Ce 4f 殻だけでなく、Co 3d 殻にも Coulomb 相関 (U=1.9 eV, J=0.8 eV) を考慮し、二サブラティスモデルを用いて補正を行うと、MAE は 4.84 meV/f.u. となり、実験値 (5.5 meV) と非常に良い一致を示しました。
• 局所近似 (HIA) との対比: 局所的な DFT+U(HIA) 近似では MAE が 7.5 meV と過大評価されるのに対し、動的相関を考慮した DFT+U(ED) は価数揺らぎによるモーメントの低減を正しく捉え、実験値に近い値を導き出しました。

4. 意義と結論 (Significance)

• 動的相関の重要性の確立: CeCo5 の磁気特性、特に Ce のモーメント低減と MAE の値は、静的な平均場近似ではなく、動的な電子相関とハイブリダイゼーション効果を考慮することで初めて正しく記述できることを実証しました。
• 理論手法の確立: DFT+DMFT のような高精度手法が MAE 計算にまだ技術的に困難である現状において、DFT+U(ED) 手法が、価数揺らぎを持つ希土類化合物の磁気異方性を定量的に予測する有効な代替手段であることを示しました。
• 次世代磁石への示唆: 低コストかつ高性能な希土類磁石（特に Ce 系）の開発において、Ce の価数揺らぎと Co 3d 電子の相関効果を同時に考慮した設計指針が得られました。これは、希土類含有量を削減した高品質な永久磁石の開発に向けた重要なステップとなります。

この研究は、CeCo5 の特異な磁気挙動を「中間価数」と「動的相関」の観点から統一的に説明し、理論と実験の間の長年のギャップを埋める成果と言えます。