Bachelorthesis: Calculation of the magnetic properties of quarternary ThMn$_{12}$-type compounds with Zr as a substitution for Nd

本学士論文は、Zr と Ti を含む Nd 置換 ThMn$_{12}$型化合物の固有特性および磁性を密度汎関数理論計算を用いて調査し、希少希土類磁石の有望な代替品となる可能性を秘めた Nd 低含有量の四元材料を同定するものである。

要約

以下は、論文「Zr を Nd の置換体とする四元 ThMn12 型化合物の磁気特性の計算」を平易な言葉と日常的な比喩を用いて解説したものです。

全体像：より安価で環境に優しい磁石の発見

現代社会が電気によって動く巨大な機械だと想像してください。この機械（電気自動車、風力タービン、ハードディスクなど）を効率的に動かすためには、強力な磁石が必要です。現在、これらの磁石における「ゴールドスタンダード」は、Nd-Fe-B（ネオジム - 鉄 - ボロン）と呼ばれる材料です。

ネオジム（Nd）を、あるパーティのVIP 客だと考えてみてください。ネオジムは磁石を驚くほど強力にしますが、高価で入手が難しく、ほぼすべての国が供給を一つの国（中国）に依存しています。これは、誰もが渡ろうとする単一の狭い橋のような供給網のボトルネックを生み出しています。

この研究の目的は、その VIP 客に代わる代替品を見つけることです。著者のニコ・ヤニック・メルトは、ジルコニウム（Zr）を使用することを提案しています。ジルコニウムは、豊富で安価、かつ入手が容易な頼もしい近所の人のような存在です。問題は、VIP を近所の人と入れ替えた場合、パーティ（磁石）は依然として同じように機能するかどうかです。

問題点：「不安定な家」

研究対象となっている特定の磁石構造は、ThMn12（または「1:12 相」）と呼ばれます。

• 設計図：1 つの希土類原子（VIP）と 12 個の鉄原子で構成される、特定の設計図で建てられた家だと想像してください。
• 問題点：もし VIP（ネオジム）と 12 個の鉄原子だけでこの家を建てようとすると、家は不安定になります。砂の基礎の上に高層ビルを建てようとするようなもので、崩壊してしまいます。
• 解決策：家が立つようにするためには、「安定化剤」が必要です。この場合、研究者たちはチタン（Ti）を使用します。チタンは、家が崩壊しないように枠組みに追加する鋼鉄製の梁のようなものです。

実験：仮想建設現場

実際の研究所でこれらの磁石を製造するのは高価で時間がかかるため、著者はスーパーコンピュータを用いて建設をシミュレーションしました。これを密度汎関数理論（DFT）と呼びます。

• シミュレーション：ビーカーで化学薬品を混ぜる代わりに、コンピュータは原子同士がどのように「感じ合っているか」を計算します。「もしここにジルコニウムを入れたら、家は立ち続けるか？磁力はどれほど強くなるか？」という問いに答えるのです。

コンピュータが見つけたこと

この論文では、ネオジムをジルコニウムに置き換えることが磁石の性能にどのような影響を与えるかを見るために、いくつかの「もしも」のシナリオを実行しました。以下が主要な発見です。

1. 安定性（家は立ち続けるか？）

• 純粋なネオジム：助けなしでは、家は不安定です。
• 純粋なジルコニウム：驚くべきことに、ジルコニウムと鉄だけで作られた家は安定しています。
• 混合（50/50）：ネオジムとジルコニウムを半々で混ぜた場合、家は少しぐらつきます。安定させるために、より多くの「鋼鉄製の梁」（チタン）を追加する必要があります。
• 結論：ネオジムをジルコニウムに置き換えることは可能ですが、構造を安定させるためにレシピに注意を払う必要があります。

2. 強度（磁石はどれほど強力か？）

• トレードオフ：VIP（ネオジム）は本質的に非常に強い磁性を持っています。一方、近所の人（ジルコニウム）は磁性が弱いです。
• 結果：ネオジムをジルコニウムに置き換えると、磁石はわずかに弱くなります。超強力なロープを、わずかに細いロープに置き換えるようなものです。しかし、磁石は依然として非常に強力であり、実用に耐える十分な強さがあります。
• 「エネルギー積」：これは磁石が蓄えることのできるエネルギーの尺度です。新しいジルコニウムベースの磁石は非常に高いスコアを記録し、いくつかの古い種類の磁石を凌駕し、現在のチャンピオン（Nd-Fe-B）に迫る結果となりました。

3. 耐熱性（キュリー温度）

• 磁石は熱くなりすぎると力を失います。「キュリー温度」とは、磁石が諦めて機能を停止するポイントです。
• 発見：新しいジルコニウム磁石は、ネオジム磁石とほぼ同様に熱に耐えることができます。高温の電気モーターの中で溶けたり、力を失ったりすることはありません。

4. 方向性（「一方通行」）

• 良い永久磁石は、消磁されにくい「硬い」ものである必要があります。一方通行の道路のように、その方向をしっかりと保持する必要があります。
• 発見：ジルコニウム磁石は方向を保持するのが非常に得意です。実際、いくつかの計算では、ジルコニウム磁石はネオジム磁石よりもさらに優れて方向を保持していました。

結論：勝者か？

この論文は、ジルコニウムがネオジムの非常に有望な代替品であると結論付けています。

• 長所：安価で豊富、かつ供給網にとって重要度が低いです。結果として得られる磁石は安定しており、優れた磁気特性を持っています。
• 短所：磁石は純粋なネオジム製のものよりもわずかに弱く、現時点では「完全に硬い」磁石（「半硬質」です）になるにはわずかに足りていません。
• 未来：著者は、少しの調整（窒素の添加やレシピの調整など）を行えば、これらのジルコニウム磁石が、現在使用されている高価なネオジム磁石に代わる現実的な代替品になりうると提案しています。

要約すると：著者はコンピュータを用いて、高価な VIP（ネオジム）の代わりに、安価で豊富な近所の人（ジルコニウム）を使って、強力で安定した磁石を構築できることを証明しました。それはまだ VIP 版ほど強力ではありませんが、電気自動車やグリーンエネルギー向けの磁石の製造方法を革命的に変えるのに十分なほど近いものです。

技術概要

ニコ・ヤニック・メルクトによる学士論文「Zr を Nd の置換元素とする四元 ThMn12 型化合物の磁性計算」の詳細な技術的概要を以下に示す。

1. 問題提起と動機

高性能永久磁石の世界的需要は、再生可能エネルギー技術と電気自動車によって牽引されている。現在、産業界はNd$_2$Fe$_{14}$B磁石に大きく依存している。しかし、ネオジム（Nd）やその他の希土類（RE）元素は、供給網の独占（主に中国）と埋蔵量の減少により、重要資源と見なされている。

• 課題: 重要元素への依存なしに高い磁気性能を維持する「希土類低含有」または「希土類フリー」磁石の開発。
• 具体的なギャップ: **ThMn$_{12}$型（1:12）**構造（RFe$_{12}$）は Nd$_2$Fe$_{14}$B よりも希土類含有量が低いものの、純粋な RFe$_{12}$相は熱力学的に不安定である。これらは安定化（通常はチタン、Ti による）を必要とし、Nd$_2$Fe$_{14}$B に比べてキュリー温度や異方性が低いという欠点にしばしば直面する。
• 目的: この研究は、1:12 相における Nd の置換元素として**ジルコニウム（Zr）**を調査する。Zr は豊富で安価であり、重要資源ではない。本研究の目的は、Zr 置換四元 (Zr,Nd)Fe$_{12-y}$Ti$_y$および五元 (Zr,Nd)Fe$_{10}$CoTi 化合物が、Nd-Fe-B 磁石の実用的で高性能な代替品となり得るかどうかを決定することである。

2. 手法

本研究は、固有の磁気特性を予測するために密度汎関数理論（DFT）に基づく第一原理計算手法を採用している。

• ソフトウェアとフレームワーク:
  • VASP（Vienna Ab-initio Simulation Package）: 構造緩和、全エネルギー、磁気モーメント、結晶異方性の計算に使用。
  • AkaiKKR（MACHIKANEYAMA）: コリング・コーン・ロストッカー（KKR）グリーン関数法を介してキュリー温度（$T_C$）を計算するために使用。
• 理論的アプローチ:
  • 交換相関汎関数: 主に**一般化勾配近似（GGA-PBE）**を使用。結晶異方性エネルギー（MAE）については、これらの系においてより正確な異方性定数をもたらす傾向があるため、**局所スピン密度近似（LSDA）**も採用された。
  • f 電子の扱い: Nd の 4f 電子は、標準的な DFT が正確に捉えることができない自己相互作用誤差と強い電子相関を修正するために、DFT+U法（Dudarev 手法、$U = 6$ eV）を用いて処理された。
  • 擬ポテンシャル: プロジェクター増幅波（PAW）法。
• 実施された計算:
  • 相安定性: Fe サイトへの Ti と Co の置換、および RE サイトへの Zr の置換の安定な構成を特定するための生成エネルギーと溶解エンタルピーの計算。
  • 磁気特性: 全磁気モーメント（$m_{tot}$）、飽和磁化（$\mu_0 M_S$）、および最大エネルギー積（$|BH|_{max}$）。
  • 異方性: 磁化方向を回転させる（[001] から [010] へ）ことによる結晶異方性エネルギー（MAE）の計算から、異方性定数（$K_1$）と異方性場（$H_a$）を導出。
  • キュリー温度（$T_C$）: 熱揺らぎを考慮するために平均場近似（MFA）と局所モーメント乱雑（LMD）モデルを用いて計算。
  • 硬さ係数（$\kappa$）: 材料が硬磁石（$\kappa > 1$）に該当するかどうかを判断するために計算。

3. 主要な貢献

1. 体系的な Zr 置換分析: この研究は、Y や Ce に焦点を当てた以前の研究を超えて、1:12 ThMn$_{12}$構造における Nd の Zr による置換に関する最初の包括的な理論的調査を提供する。
2. 安定化メカニズム: 純粋な (Zr,Nd)Fe$_{12}$は不安定であるが、Ti（8i サイトへ）と Co（8f/8j サイトへ）の添加が相を安定化することを確認した。また、Zr 置換により、純粋な Nd 系アナログと比較して Ti 濃度を低減できることを確立した。
3. 手法の検証: これらの系における異方性の計算にLSDAを使用することを検証した。GGA+U は Nd を含む化合物に対して $K_1$ を過小評価する傾向があるのに対し、LSDA は実験データに近い結果を提供することを指摘した。
4. 有望な候補の特定: コスト（Nd 削減）と性能のバランスが取れた有望な候補として、(Zr,Nd)Fe$_{11}$Tiおよび**(Zr,Nd)Fe$_{10}$CoTi**を特定した。

4. 主要な結果

A. 安定性と構造

• Ti による安定化: Ti は優先的に8i結晶学サイトを占有する。約 7.7 at.% の Ti 濃度（26 原子超格子あたり 2 原子）は、(Zr,Nd)Fe$_{12}$相を安定化するのに十分である。
• Co 置換: Co は8fおよび8jサイトを好む。Co は磁気モーメントをわずかに改善するが、Ti に比べて安定化効果は小さい。
• 格子定数: Zr の原子半径が小さいため、Zr 置換は純粋な Nd 化合物に比べて格子体積の減少をもたらす。傾向は以下の通り：$Nd_2Fe_{24} > ZrNdFe_{24} > Zr_2Fe_{24}$。

B. 磁気モーメントと飽和磁化（$\mu_0 M_S$）

• Zr の影響: Nd を Zr に置換すると、Zr の磁気モーメント（-0.47 $\mu_B$）が Nd（GGA では -0.27 $\mu_B$、ただし DFT+U を通して f 電子を含めると -3.3 $\mu_B$）よりも低いため、全磁気モーメントが減少する。
• 計算値:
  • NdFe$_{11}$Ti: $\mu_0 M_S \approx 1.69$ T（DFT+U）、実験範囲（1.38–1.70 T）と一致。
  • (Zr,Nd)Fe$_{11}$Ti: $\mu_0 M_S \approx 1.62$ T。
  • ZrFe$_{11}$Ti: $\mu_0 M_S \approx 1.51$ T。
• 結論: Zr 置換は $M_S$ をわずかに低下させるが、値は実用的な応用にとって十分に高いままである。

C. 最大エネルギー積（$|BH|_{max}$）

• 性能: 四元化合物**(Zr,Nd)Fe$_{11.5}$Ti$_{0.5}$は、理論的な $|BH|_{max}$ が約 600 kJ/m³**を示し、Nd$_2$Fe$_{14}$B の理論限界（約 512 kJ/m³）を超えている。
• 安定相: 安定な**(Zr,Nd)Fe$_{11}$Ti化合物は約 510–525 kJ/m³**を生み出し、Nd$_2$Fe$_{14}$B と同等であり、Sm$_2$Co$_{17}$（294 kJ/m³）よりも著しく高い。
• Co の効果: Co 置換は $M_S$ をわずかに増加させるが、体積変化により $|BH|_{max}$ をわずかに減少させる。

D. キュリー温度（$T_C$）

• 傾向: $T_C$は Ti 含有量の増加とともに低下するが、Co 置換とともに増加する。
• 値:
  • NdFe$_{11}$Ti: $T_C \approx 652$ K（LMD）/ 792 K（FM）。
  • (Zr,Nd)Fe$_{11}$Ti: $T_C \approx 652$ K（LMD）/ 783 K（FM）。
  • 比較: これらの値は Nd$_2$Fe$_{14}$B（$T_C \approx 588$ K）よりも高く、高温用途に適している。
• 注記: FM 近似は $T_C$を過大評価する。LMD 値の方がより現実的と見なされるが、依然として実用的な動作温度を示している。

E. 結晶異方性（$K_1$）と硬さ係数（$\kappa$）

• 異方性定数（$K_1$）: Zr 置換は $K_1$を著しく増加させる。
  • NdFe$_{11}$Ti: $K_1 \approx 1.68$ MJ/m³（LSDA）。
  • (Zr,Nd)Fe$_{11}$Ti: $K_1 \approx 2.41$ MJ/m³（LSDA）。
  • ZrFe$_{11}$Ti: $K_1 \approx 2.47$ MJ/m³（LSDA）。
  • 意義: Zr は保磁力にとって重要な特性である異方性を増強する。
• 硬さ係数（$\kappa$）:
  • (Zr,Nd)Fe$_{11}$Ti は $\kappa \approx 0.97$（LSDA）を生み出す。
  • これは材料を半硬（$\kappa < 1$）に分類するが、Nd$_2$Fe$_{14}$B は硬（$\kappa \approx 1.54$）である。
  • この値は硬磁石の閾値に非常に近いため、わずかな修正（例えば、窒素化など）により硬磁石領域に押し上げられる可能性があることを示唆している。

5. 意義と展望

• 戦略的重要性: この研究は、ジルコニウムが、飽和磁化やキュリー温度などの重要な磁気特性を犠牲にすることなく、ThMn$_{12}$型磁石におけるネオジムの代替として、低コストかつ豊富に存在する viable な選択肢であることを実証している。
• 性能の可能性: 計算された特性は、(Zr,Nd)Fe$_{11}$Ti 基磁石がエネルギー積と熱的安定性の面で Nd-Fe-B と競合でき、より安全な供給網を提供する可能性を示唆している。
• 将来の方向性:
  • 実験的検証: 理論的予測は、合成（例えば、選択的レーザー溶融による）と実験的検証を必要とする。
  • 最適化: 硬さ係数が 1.0 をわずかに下回っているため、将来の研究では、$\kappa > 1.0$に押し上げるための窒素化（R-Fe-N 相の形成）またはさらなる Co/Ti の調整を探求すべきである。
  • 拡大: この手法は、他の希土類低含有構造（2:14、2:17 相）や他の豊富元素（例えば、La、Ce）に適用され、重要希土類への依存をさらに軽減できる。

結論として、本論文は**(Zr,Nd)Fe$_{12-y}$Ti$_y$**化合物を、次世代の持続可能で高性能な永久磁石の有望な候補として特定することに成功した。