四重極応答と磁気弾性

本稿は、(i) 四重極応答理論（四重極子演算子と感受率による統一的記述）と、(ii) 磁歪・磁気弾性（磁化と格子ひずみの結合）を、「同じrank-2（2階）テンソル量が同じ対称性チャネルで結合する」という一点で接続し、式の形で整理するものである。超音波（弾性定数の軟化）、共鳴X線散乱（E2寄与や干渉散乱）、NMR/NQR（EFGと核四重極相互作用）、およびXMCD/XMLD/XLD（二色性テンソル解析）を、重複なく一つの見取り図にまとめる。

参考ドキュメント

1. 多極子・四重極子の定義と点群既約表現の整理（多極子入門）
2. 超音波と多極子（弾性定数と四重極感受率の関係、軟化の読み方）
3. E. Callen and H. Callen, “Magnetostriction, Forced Magnetostriction, and Anomalous Thermal Expansion in Ferromagnets,” Phys. Rev. 139, A455 (1965).

0. 何を「四重極」と呼ぶか

• 本稿での四重極子は、空間分布の「2階の異方性（rank-2の対称・無跡テンソル成分）」を表す量である。
• 四重極“応答”とは、外場勾配・ひずみ・結晶場変調などの「四重極子に共役な一般化外場」に対する線形応答（感受率）と、その臨界増大や動的構造である。
• 磁歪は、磁化に伴う格子ひずみであるが、その自由エネルギー記述は「ひずみ（rank-2）×磁化方向（rank-2の二次形式）」の結合になり、対称性の言葉では四重極子と同型のチャネルで整理できる。ここに統一的記述の核心がある。

1. 四重極子の定義

1.1 連続分布の四重極テンソル（電気四重極子）

電荷密度 $\rho (\mathbf{r})$ の二次の異方性（トレース除去を含む代表的定義）として

$$Q_{ij}=\int d^{3}r\rho (\mathbf{r})(3r_i{r}_j-r^2{\delta }_{ij})$$

を用いる流儀が広い。係数や無跡化の取り方は分野規約があるため、異分野（散乱・NMR・局在模型）を横断する際は「採用規約」を明示するのが必須である。

固体中の局在軌道の議論では演算子として

$${\hat{Q}}_{ij}=3{\hat{x}}_i{\hat{x}}_j-{\hat{r}}^2{\delta }_{ij},Q_{ij}=⟨{\hat{Q}}_{ij}⟩$$

と書く。$d$軌道占有の偏りや結晶場・ひずみによる混成は、$Q_{ij}$ を直接変化させる。

1.2 点群既約表現での四重極子

結晶中では $Q_{ij}$ の成分より、「点群既約表現（対称性チャネル）に沿った基底」へ写すほうが実用的である。立方晶を例にすると、代表的に

• ${\mathrm{\Gamma }}_3$（Eg型；tetragonal成分）

  $${\epsilon }_u\propto (2{\epsilon }_{zz}-{\epsilon }_{xx}-{\epsilon }_{yy}),{\epsilon }_v\propto ({\epsilon }_{xx}-{\epsilon }_{yy})$$

  に共役な四重極子（例：$O_2^0$、$O_2^2$ 型）が現れる。

• ${\mathrm{\Gamma }}_5$（T2g型；shear成分）

  $${\epsilon }_{xy},{\epsilon }_{yz},{\epsilon }_{zx}$$

  に共役な四重極子（例：$O_{xy},O_{yz},O_{zx}$）が現れる。

局在$f$電子では Stevens 演算子として

$${\hat{O}}_2^0=3J_z^2-J(J+1),{\hat{O}}_2^2=J_x^2-J_y^2,{\hat{O}}_{xy}=J_x{J}_y+J_y{J}_x$$

などが典型基底になる。遍歴$d$電子でも、対称性チャネルの整理は同様に有効である。

1.3 電気四重極子と磁気四重極子

時間反転に対する性質で、多極子は大別できる。

• 電気多極子：時間反転に偶（電荷分布起源）
• 磁気多極子：時間反転に奇（磁化分布起源）
• 磁気四重極モーメント（MQM）：磁化の空間分布の2階モーメントに相当し、磁気電気応答や温度勾配誘起磁化などの形式論で重要となる。

近年の議論では、単純なKubo項だけで書いた応答が不自然な項（境界寄与やゲージ問題に類する整理）を含む場合に、（スピン／軌道の）四重極モーメントに由来する補正（差し引き）を導入して「物理的な応答」へ整形する枠組みが議論されている。

2. 線形応答としての四重極感受率

2.1 四重極子は「外場の勾配」や「ひずみ」に結合する

双極子が電場 $\mathbf{E}$ に一次で結合するのに対し、電気四重極子は場の勾配に結合する：

$$H_{\text{int}}=-\frac{1}{6}\sum _{ij}{Q}_{ij}{\partial }_i{E}_j$$

（係数は規約）。固体物理では、場の勾配そのものよりも「結晶場の変調」や「ひずみ $\epsilon$」が、四重極子に共役な一般化外場として現れることが多い。

2.2 Kubo公式（一般化外場 $f$ に対する応答）

四重極子演算子 $Q$（または $O_{\mathrm{\Gamma }\gamma }$）に共役な外場 $f$ を導入して

$$H_{\text{int}}=-f(t)Q$$

とおくと、線形応答は

$$⟨Q⟩(\omega )={\chi }_{QQ}(\omega )f(\omega )$$

で定義され、Kubo公式は

$${\chi }_{QQ}(\omega )=\frac{i}{\hslash }{\int }_0^{\mathrm{\infty }}dt{e}^{i\omega t}⟨[Q(t),Q(0)]⟩$$

で与えられる。静的極限 $\omega \to 0$ の増大は、四重極秩序や電子ネマティックなどの「rank-2秩序」に近いゆらぎの増大を反映する。

2.3 対称性チャネルごとの感受率

実際の解析では、四重極子を $O_{\mathrm{\Gamma }\gamma }$ に分解し

$$⟨O_{\mathrm{\Gamma }\gamma }⟩={\chi }_{\mathrm{\Gamma }}(\omega )f_{\mathrm{\Gamma }\gamma }(\omega )$$

のように「どの $\mathrm{\Gamma }$ が軟いか」を追う。RPAや平均場では多極子感受率が行列になるため、基底選択（軌道・スピン・点群）が実装上の鍵となる。

---

3. 四重極子–ひずみ結合と弾性異常

3.1 四重極子–ひずみ結合

結晶中では、ひずみ ${\epsilon }_{\mathrm{\Gamma }\gamma }$ が局所結晶場を変調し、同じ対称性の四重極子と線形結合する：

$$H_{QS}=-g_{\mathrm{\Gamma }}\sum _{\gamma }{O}_{\mathrm{\Gamma }\gamma }{\epsilon }_{\mathrm{\Gamma }\gamma }$$

ここで $g_{\mathrm{\Gamma }}$ は結合定数である。これは「四重極応答」と「格子応答」を結ぶ最小結合である。

3.2 弾性定数への寄与

弾性定数は自由エネルギー $F$ の二階微分

$$C_{\mathrm{\Gamma }}=\frac{{\partial }^{2}F}{\partial {\epsilon }_{\mathrm{\Gamma }\gamma }^2}$$

で定義される。電子系の四重極応答があると、代表的に

$$C_{\mathrm{\Gamma }}(T)=C_{\mathrm{\Gamma }}^{(0)}-g_{\mathrm{\Gamma }}^2{\chi }_{\mathrm{\Gamma }}(T)$$

の形で軟化が現れる。さらに四重極子同士の相互作用（平均場的結合）を入れると

$$C_{\mathrm{\Gamma }}(T)=C_{\mathrm{\Gamma }}^0(T)-\frac{g_{\mathrm{\Gamma }}^2{\chi }_{\mathrm{\Gamma }}(T)}{1-g_{\mathrm{\Gamma }}^{\prime }{\chi }_{\mathrm{\Gamma }}(T)}$$

のような形が用いられる。超音波測定は「特定対称性の ${\chi }_{\mathrm{\Gamma }}$」を選択的に読む手段である、という位置づけが明確になる。

3.3 静的と動的の区別

超音波（MHz）は準静的に近いが、RIXSはeVスケールの動的構造因子を読む。従って「同じ四重極」という言葉でも、見ている $\chi (\omega )$ の周波数領域が異なる。議論では、$\omega$ の領域を必ず揃える（あるいは揃わないことを明示する）必要がある。

---

4. 磁歪と磁気弾性

本節は「四重極応答一般」ではなく、磁性体で磁歪・磁気弾性定数を議論するための最小式をまとめ、四重極子チャネル（${\mathrm{\Gamma }}_3,{\mathrm{\Gamma }}_5$ 等）と1対1に対応づけるのが目的である。

4.1 磁歪の定義（現象論）

磁歪は、磁化方向（あるいは磁区配置）の変化に伴って格子がひずむ現象である。ひずみテンソル ${\epsilon }_{ij}$ の変化として観測され、逆磁歪（Villari効果）では応力が磁化方向を変える。

立方晶の代表的表現として、磁化方向余弦 $\mathit{\alpha }=({\alpha }_x,{\alpha }_y,{\alpha }_z)$ を用い

$${\epsilon }_{xx}=\frac{3}{2}{\lambda }_{100}({\alpha }_x^2-\frac{1}{3}),{\epsilon }_{xy}=\frac{3}{2}{\lambda }_{111}{\alpha }_x{\alpha }_y$$

などが用いられる（座標系・方位・定義規約の明示が必須）。

4.2 立方晶の一次磁気弾性エネルギー（$b_1,b_2$）

磁化方向とひずみの結合は自由エネルギーの展開で

$$F(\mathit{\alpha },\epsilon )=F_{\mathrm{mag}}(\mathit{\alpha })+F_{\mathrm{el}}(\epsilon )+F_{\mathrm{me}}(\mathit{\alpha },\epsilon )$$

と分けられ、立方晶の一次磁気弾性項は

$$F_{\mathrm{me}}=b_1({\epsilon }_{xx}{\alpha }_x^2+{\epsilon }_{yy}{\alpha }_y^2+{\epsilon }_{zz}{\alpha }_z^2)+2b_2({\epsilon }_{xy}{\alpha }_x{\alpha }_y+{\epsilon }_{yz}{\alpha }_y{\alpha }_z+{\epsilon }_{zx}{\alpha }_z{\alpha }_x)$$

で与えられる。

ここで重要なのは、${\alpha }_i^2$ や ${\alpha }_i{\alpha }_j$ が「rank-2の基底（無跡成分）」を作るため、$F_{\mathrm{me}}$ が四重極子–ひずみ結合（第3章）の特定実現になっている点である。すなわち

• ${\mathrm{\Gamma }}_3$（tetragonal）チャネル ↔ $b_1$ 系の結合
• ${\mathrm{\Gamma }}_5$（shear）チャネル ↔ $b_2$ 系の結合 という対応が成立する。

4.3 $b_1,b_2$ と磁歪定数 ${\lambda }_{100},{\lambda }_{111}$（弾性定数を介した関係）

弾性定数 $(c_{11},c_{12},c_{44})$ を用いる代表式として

$$b_1=-\frac{3}{2}(c_{11}-c_{12}){\lambda }_{100},b_2=-3c_{44}{\lambda }_{111}$$

が広く用いられる。従って「磁歪（$\lambda$）」と「弾性（$c$）」を押さえることは、磁気弾性定数（$b_1,b_2$）を実験的に決める最短経路である。

薄膜では基板拘束、面内等方・異方、テクスチャの平均化が入るため、実効定数（面内・面外の拘束条件に沿った定義）へ置き換える必要がある。

4.4 四重極子モーメント（電子密度異方性）としての理解

磁歪の微視的起源は、スピン軌道相互作用が軌道占有（電荷分布の異方性＝四重極子成分）を再配列し、その四重極子成分がひずみ（結晶場変調）に敏感であることにある。 従って「四重極子成分を実験的に追跡できる」なら、磁歪を格子の変形ではなく電子状態の再配列として分解し、元素・軌道別に議論できる。

5. 実験プローブ

5.1 超音波（弾性定数）：${\chi }_{\mathrm{\Gamma }}$ を直接読む

第3章の関係式

$$C_{\mathrm{\Gamma }}=C_{\mathrm{\Gamma }}^{(0)}-g_{\mathrm{\Gamma }}^2{\chi }_{\mathrm{\Gamma }}$$

により、温度依存の軟化・硬化・異常から四重極感受率や相互作用の寄与を議論できる。縦波・横波・伝播方向・偏波の選択で、$\mathrm{\Gamma }$ を切り分けられるのが強みである。

5.2 共鳴X線散乱（RXS）：E2寄与・干渉散乱で四重極に感度を持たせる

吸収端近傍では散乱振幅が共鳴項で増強され、E1（双極子遷移）だけでなくE2（四重極遷移）が寄与する場合がある。禁制反射の出現、方位角依存、偏光解析、E1–E2干渉などを用いることで、四重極（テンソル）成分や多極子秩序に対する選択性を設計できる。

5.3 RIXS：動的四重極励起（対称性分解）の読み出し

RIXSは局所励起（電荷・軌道・スピン励起）をエネルギー損失として読み、選択則に多極子成分が入る。静的秩序（超音波・回折）とは別に、動的な励起対称性を押さえる位置づけである。

5.4 NMR/NQR：EFG（電場勾配）で局所四重極を読む

核スピン $I\ge 1$ は核四重極モーメント $Q_n$ を持ち、電場勾配（EFG）テンソル $V_{ij}={\partial }_i{\partial }_{j}V(\mathbf{r})$ と結合する。代表形として

$$H_Q=\frac{e{Q}_n}{4I(2I-1)}[3I_z^2-I(I+1)+\frac{\eta }{2}(I_{+}^2+I_{-}^2)]V_{zz},\eta =\frac{V_{xx}-V_{yy}}{V_{zz}}$$

で書かれる（主軸系で $|V_{zz}|\ge |V_{yy}|\ge |V_{xx}|$ を採るのが一般的）。 電子系の四重極秩序や電荷分布異方性が変わるとEFGが変化し、NQR/NMRの線分裂・線幅・周波数シフトに現れる。

6. 理論モデル：四重極感受率と秩序の増大機構（詳細版）

本章では、四重極応答（四重極感受率 ${\chi }_{\mathrm{\Gamma }}$ の増大、あるいは四重極秩序の形成）が、どのような微視的機構で生じるかを整理する。大きく分けると、(i) 局在電子の結晶場準位に由来する「単イオン（single-ion）機構」、(ii) 遍歴電子の多軌道相関に由来する「集団（itinerant/collective）機構」、(iii) 四重極に限らない高次多極子や格子自由度との結合によって「相転移の形そのものが変わる」機構、の3系統がある。

四重極応答の理論では、観測量（超音波の弾性定数、RXSの禁制反射、NMR/NQRのEFG変化など）を説明するために、

• どの対称性チャネル $\mathrm{\Gamma }$ が軟いか（${\chi }_{\mathrm{\Gamma }}$ が増大するか）
• それが「単イオンの縮退」なのか「相互作用による臨界ゆらぎ」なのか
• 静的（$\omega \to 0$）と動的（有限 $\omega$）のどちらを議論しているか を最初に切り分ける必要がある。

---

6.1 結晶場準位（単イオン）と ${\chi }_{\mathrm{\Gamma }}(T)$

6.1.1 単イオン四重極応答の基本構造

局在$f$電子（希土類・アクチノイド）では、スピン軌道相互作用により全角運動量多重項（$J$ 多重項）が形成され、その後に結晶場（crystal electric field; CEF）が縮退をさらに分裂させる。四重極子演算子 $O_{\mathrm{\Gamma }\gamma }$ はこの局在多重項空間に作用するため、四重極感受率 ${\chi }_{\mathrm{\Gamma }}(T)$ は「CEF準位構造」と「その準位間の四重極子行列要素」によってほぼ決まる。

外場（一般化外場）として、四重極子に共役な量 $f_{\mathrm{\Gamma }\gamma }$ を導入し

$$H_{\text{int}}=-\sum _{\gamma }{f}_{\mathrm{\Gamma }\gamma }{O}_{\mathrm{\Gamma }\gamma }$$

とすると、線形応答

$$⟨O_{\mathrm{\Gamma }\gamma }⟩={\chi }_{\mathrm{\Gamma }}(T)f_{\mathrm{\Gamma }\gamma }$$

で ${\chi }_{\mathrm{\Gamma }}(T)$ を定義できる。

単イオンの静的感受率は、CEF固有状態 $|n⟩$（エネルギー $E_n$）を用いて、一般に「Curie項（縮退に由来）」と「Van Vleck項（準位間混合に由来）」に分解される形で書ける。典型的には（係数や縮約は省略して）

$${\chi }_{\mathrm{\Gamma }}(T)=\underset{\text{Curie型（同一準位内のゆらぎ）}}{\underbrace{\frac{1}{k_{B}T}(⟨O_{\mathrm{\Gamma }}^2{⟩}_T-⟨O_{\mathrm{\Gamma }}{⟩}_T^2)}}+\underset{\text{Van Vleck型（準位間混合）}}{\underbrace{\sum _{n\ne m}\frac{p_n-p_m}{E_m-E_n}|⟨n|O_{\mathrm{\Gamma }}|m⟩{|}^2}}$$

のような構造になる（$p_n=e^{-\beta E_n}/Z$）。

重要なのは、

• 低温で基底状態が「四重極子に対して縮退（あるいは近縮退）」を持つと、Curie型に ${\chi }_{\mathrm{\Gamma }}\propto 1/T$ が現れやすい。
• 基底が非縮退でも、近接準位との行列要素が大きいと、Van Vleck型の“定数に近い”寄与が大きくなり得る。 という点である。

6.1.2 「$1/T$ 的増大」と「弾性定数の軟化」の関係

四重極子–ひずみ結合

$$H_{QS}=-g_{\mathrm{\Gamma }}\sum _{\gamma }{O}_{\mathrm{\Gamma }\gamma }{\epsilon }_{\mathrm{\Gamma }\gamma }$$

があると、弾性定数は

$$C_{\mathrm{\Gamma }}(T)=C_{\mathrm{\Gamma }}^{(0)}-g_{\mathrm{\Gamma }}^2{\chi }_{\mathrm{\Gamma }}(T)$$

（あるいは相互作用込みで $C_{\mathrm{\Gamma }}=C_{\mathrm{\Gamma }}^0-{\displaystyle \frac{g_{\mathrm{\Gamma }}^2{\chi }_{\mathrm{\Gamma }}}{1-g_{\mathrm{\Gamma }}^{\prime }{\chi }_{\mathrm{\Gamma }}}}$） となるため、${\chi }_{\mathrm{\Gamma }}(T)$ の増大は $C_{\mathrm{\Gamma }}(T)$ の軟化として観測される。

単イオン機構では、

• ${\chi }_{\mathrm{\Gamma }}(T)\simeq A+{\displaystyle \frac{B}{T}}$ のような「定数項＋キュリー項」がまず現れる
• その結果、$C_{\mathrm{\Gamma }}(T)$ は「なだらかな背景＋低温での急な低下」という形を取りやすい という経験則が成り立つ。したがって、弾性定数の温度依存をまずこの形でフィットし、CEFパラメータ（準位間隔、行列要素）との整合を検討するのが、単イオン四重極応答の実務的な第一歩になる。

6.1.3 反強四重極秩序への発展：平均場の見取り図

単イオン応答が大きくても、それだけでは秩序は生じない。秩序（反強四重極など）は、サイト間相互作用（四重極–四重極相互作用）により安定化される。最も単純な有効相互作用を

$$H_{QQ}=-g_{\mathrm{\Gamma }}^{\prime }\sum _{⟨ij⟩}{O}_{\mathrm{\Gamma }}(i)O_{\mathrm{\Gamma }}(j)$$

のように入れると、平均場的には

$${\chi }_{\mathrm{\Gamma }}^{\text{(MF)}}(T)=\frac{{\chi }_{\mathrm{\Gamma }}^{(0)}(T)}{1-g_{\mathrm{\Gamma }}^{\prime }{\chi }_{\mathrm{\Gamma }}^{(0)}(T)}$$

が得られ、分母がゼロになる温度が転移点 $T_Q$ を与える。ここで ${\chi }_{\mathrm{\Gamma }}^{(0)}(T)$ が単イオン感受率（前節）である。

• 単イオンの $1/T$ 増大が強いほど、比較的高温で $g_{\mathrm{\Gamma }}^{\prime }{\chi }_{\mathrm{\Gamma }}^{(0)}(T)=1$ に達しやすい
• 同じ材料でも、どの $\mathrm{\Gamma }$ チャネルの行列要素が大きいかで、実際に現れる秩序の対称性が変わる という“対称性選択”がここで起こる。

6.1.4 実験的な判別のコツ（単イオンか、臨界か）

単イオン由来かどうかを見抜くには、次の観点が有効である。

• CEF準位間隔が既知（中性子散乱・比熱・磁化など）なら、その準位構造で ${\chi }_{\mathrm{\Gamma }}(T)$ を再現できるか。
• 軟化が「特定の弾性モード（特定 $\mathrm{\Gamma }$）」に強く現れ、他のモードが比較的静かか（対称性選択が明瞭か）。
• 高温側で ${\chi }_{\mathrm{\Gamma }}(T)$ が $1/T$ に近い立ち上がりを見せるか（ただし高温では格子・電子背景も混ざるため注意）。

---

6.2 多軌道模型（遍歴性）と行列感受率

6.2.1 なぜ多軌道が本質か：四重極＝軌道自由度の組み替え

遍歴電子系（特に遷移金属$d$電子系）で四重極（ネマティック）応答が現れるとき、核心は「軌道占有や軌道混成の異方性（rank-2成分）」が自発的に増大する点にある。単一バンドの密度ゆらぎだけでは、対称性チャネルとしての四重極（例：$d_{xz}$ と $d_{yz}$ の不均衡、$t_{2g}$ 内の選択）を十分に表現できないことが多い。したがって多軌道（軌道行列）としての感受率を扱う必要がある。

6.2.2 裸感受率は「バンド構造＋軌道成分」の畳み込み

多軌道系の裸感受率は、一般に“軌道添字付き”の行列として

$${\chi }_{{\ell }_1{\ell }_2,{\ell }_3{\ell }_4}^{(0)}(\mathbf{q},i{\omega }_n)=-\frac{T}{N}\sum _{\mathbf{k},m}{G}_{{\ell }_3{\ell }_1}(\mathbf{k},i{ϵ}_m)G_{{\ell }_2{\ell }_4}(\mathbf{k}+\mathbf{q},i{ϵ}_m+i{\omega }_n)$$

のように定義される（$\ell$ は軌道添字）。ここには

• フェルミ面近傍の粒子・空孔励起（nestingやバンド交差）
• Bloch状態の軌道混成（“どの軌道成分がどの $\mathbf{k}$ に乗っているか”） が同時に入る。

この時点で、四重極チャネル（ネマティック）とは「特定の軌道組の線形結合」である。例えば

$$O_{\mathrm{\Gamma }}\sim n_{xz}-n_{yz}$$

のような“軌道占有差”は、点群の既約表現（例：$B_{2g}$ など）として四重極的に振る舞う。

6.2.3 RPA（またはそれに準じた）増大条件と「どのチャネルが勝つか」

相互作用行列 $\hat{U}$ を導入すると、RPA型には

$$\hat{\chi }(\mathbf{q},\omega )={[\hat{1}-{\hat{\chi }}^{(0)}(\mathbf{q},\omega )\hat{U}]}^{-1}{\hat{\chi }}^{(0)}(\mathbf{q},\omega )$$

となる。ここでの物理的ポイントは“行列の固有値問題”である。

• 行列 ${\hat{\chi }}^{(0)}\hat{U}$ の最大固有値が 1 に近づくと、その固有ベクトルに対応する多極子チャネルが増大する。
• 固有ベクトルは「軌道・スピン・対称性の組み合わせ」になり、これが“どの秩序（どの四重極）が主役か”を決める。

したがって理論実装では、

1. 点群の既約表現に沿って“チャネル演算子” $O_{\mathrm{\Gamma }\gamma }$ を基底として定義し
2. ${\chi }_{\mathrm{\Gamma }\gamma ,{\mathrm{\Gamma }}^{\prime }{\gamma }^{\prime }}$ を構成して
3. 対角化または最大固有値を追い
4. どの $\mathrm{\Gamma }$ が最も早く増大するかを評価する
   という流れが基本になる。

6.2.4 四重極（ネマティック）と磁性の競合・協奏

遍歴系では、四重極（軌道／電荷）チャネルとスピン（磁性）チャネルが独立ではなく、しばしば結合する。例えば

• スピンゆらぎが強いと、それが媒介する有効相互作用として四重極チャネルが強化される（“磁性ゆらぎ駆動ネマティック”のような見取り図）。
• 逆に、軌道（四重極）秩序が先に立つと、バンドの縮退が解けてスピンゆらぎの波数構造が変わる。

この種の相互作用の切り分けは、単に ${\chi }_{\mathrm{\Gamma }}(T)$ の温度依存だけでなく、$\mathbf{q}$ 依存（どの波数で増大するか）や、応答の異方性（例：${\mathrm{\Gamma }}_3$ vs ${\mathrm{\Gamma }}_5$）も含めて議論する必要がある。

6.2.5 ひずみとの結合が「見かけの転移」を作ることがある

四重極子–ひずみ結合は、電子系の秩序だけでなく格子を巻き込む。結果として

• 電子系の“裸”の臨界点より高温で、格子結合により実効転移が引き上がる
• あるいは、ひずみが外場として働き、真の相転移ではなくクロスオーバーとして観測される などが起こり得る。

したがって遍歴系では、${\chi }_{\mathrm{\Gamma }}$ の増大を議論する際に「格子との結合強度 $g_{\mathrm{\Gamma }}$」を同時に推定しないと、実験で見える軟化の形を取り違える危険がある。

---

6.3 高次多極子との結合：相転移の段階性と“見え方”を変える

6.3.1 なぜ高次多極子が効くのか

四重極秩序が現れる系では、同じ局在自由度（例：$f$電子の $J$ 多重項）が、八極子（rank-3）や十六極子（rank-4）なども同時に持ち得る。さらに、結晶対称性の下では

• ある多極子が“主秩序パラメータ”として立つ
• 別の多極子が“随伴（secondary）秩序”として誘起される
  という現象が自然に起こる。

このとき、四重極感受率だけを追っていると、

• 転移が二段階に見える
• 弾性異常の形が単純なCurie-Weiss型からずれる
• 散乱で見える禁制反射が、四重極由来か八極子由来か判別が難しい といった状況が生じ得る。

6.3.2 Landau理論の最小形：結合項が“転移の顔”を変える

四重極秩序パラメータを $Q$、高次多極子を $M$（例：八極子）と書くと、自由エネルギーは最小限でも

$$F=\frac{a_Q}{2}{Q}^2+\frac{b_Q}{4}{Q}^4+\frac{a_M}{2}{M}^2+\frac{b_M}{4}{M}^4+\lambda Q^{2}M+\cdots$$

のような結合項を持ち得る（対称性が許す場合）。このような項があると、

• $Q$ の増大が $M$ を随伴的に誘起する
• $M$ の臨界ゆらぎが $Q$ の有効係数 $a_Q$ を強く renormalize する などにより、転移温度や転移次数、観測される異常の形が変化する。

また、ひずみ $\epsilon$ を含めて

$$F=\frac{C^0}{2}{\epsilon }^2-gQ\epsilon +\frac{a_Q}{2}{Q}^2+\cdots$$

を最小化すると、$\epsilon$ を消去した実効相互作用として $-g^2{Q}^2/(2C^0)$ が生じ、格子が秩序形成を助長する（あるいは転移点をシフトさせる）という見方も得られる。

6.3.3 同定戦略：何を組み合わせると主役が見えるか

高次多極子が絡む場合、単一手法では主役同定が難しい。相補的に押さえるという意味は、具体的には次の分担である。

• 超音波（弾性定数）：どの $\mathrm{\Gamma }$ チャネルが強く軟化するか（四重極子–ひずみ結合に直結）
• RXS：禁制反射の出現条件、方位角依存、偏光解析（テンソル構造の選別）
• NMR/NQR：局所EFGの変化（電荷分布対称性の変化を局所的に確認）
• 熱力学量（比熱、熱膨張、磁化など）：転移の次数、エントロピー変化、随伴秩序の兆候
• （可能なら）外場応答：一軸応力や磁場による転移点シフト、ドメイン選択

これらを同じ対称性基底（$\mathrm{\Gamma }$）で整理して整合性を取ると、

• 四重極が主役で、八極子が随伴か
• その逆か
• あるいは両者が強く混成して“固有モード”として秩序化するか を、かなり明確に切り分けられる。

3つの機構をどう見分けるか

• 単イオン機構：CEF準位構造と行列要素で ${\chi }_{\mathrm{\Gamma }}(T)$ が決まり、$1/T$ 的増大や「定数項＋キュリー項」が見えやすい。
• 遍歴・集団機構：${\hat{\chi }}^{(0)}\hat{U}$ の固有値問題として“どのチャネルが勝つか”が決まり、$\mathbf{q}$ 依存や軌道混成が本質。
• 高次多極子結合：結合項が転移の段階性・次数・観測異常の形を変え、単一手法では主役同定が難しいため、対称性基底で複数プローブを突き合わせるのが有効。

（注）本章は「枠組みの整理」に徹しており、個別材料への適用では、点群（立方・四方など）、対象核種（NQR）、測定周波数帯（超音波 vs RIXS）などの条件を明示して、$\mathrm{\Gamma }$ の定義と外場（ひずみ・応力・場勾配）の対応を固定することが重要である。

7. 第一原理計算との接続

本章では、四重極応答（四重極感受率・秩序・ゆらぎ）や磁気弾性（磁歪・磁気弾性定数）を議論するうえで、第一原理計算（主にDFT、必要に応じてDFT+U/DMFT等）をどのように位置づけ、何を出力し、どの実験量とどう比較するかを、実務の観点で整理する。
要点は、(i) EFGのように「計算→実験量へ直結する量」をまず押さえること、(ii) ひずみ応答やSOCを含む“差分”から磁気弾性の電子論的起源を検証すること、(iii) 模型（第6章）の示す「増大機構」を、現実材料の「符号・大きさ・元素別寄与」で裏付けること、である。

---

7.1 EFG計算

7.1.1 NMR/NQRが必要とする量

核四重極相互作用で中心となるのは、核位置の電場勾配（electric field gradient; EFG）テンソル

$$V_{ij}=\frac{{\partial }^{2}V(\mathbf{r})}{\partial x_i\partial x_j}{|}_{\mathbf{r}={\mathbf{R}}_n}$$

である。これは「核位置での静電ポテンシャルの2階微分」であり、局所電荷分布の非球対称性（rank-2成分）を最も直接的に反映する。

EFGテンソルを主軸系で対角化し、主値を $V_{xx},V_{yy},V_{zz}$（慣例として $|V_{zz}|\ge |V_{yy}|\ge |V_{xx}|$）とすると、NQR/NMRで現れる主要パラメータは

• 主値 $V_{zz}$
• 非対称度

$$\eta =\frac{V_{xx}-V_{yy}}{V_{zz}}$$

である。したがって、DFTで比較すべき第一目標は $V_{zz}$ と $\eta$ を同じ規約で安定に評価することになる。

7.1.2 DFTでEFGを計算する

EFGは核近傍の電荷密度分布に敏感である。したがって実務上は、

• 全電子（all-electron）法（FLAPWなど）
• PAW法でのEFG評価（VASP等、PAW再構成により近核情報を回復） のいずれかで計算することが多い。

ここで重要なのは、EFGが「全電荷密度の角度依存（非球対称性）」に強く依存するため、次の点に注意して収束を取ることである。

• 構造の精密化：内部座標の微小差がEFGを大きく動かす場合がある
• $k$点密度：金属では特に、フェルミ面近傍の占有がEFGに影響しうる
• カットオフ（または基底サイズ）：近核の角度成分が十分表現できるか
• 電子温度（スメアリング）：金属では占有の丸めがEFGに影響する
• 磁性状態：自発磁化の向き（SOCありの場合）や磁気秩序の仮定が局所対称性を変える
• SOCの有無：重元素や強い軌道混成がある場合、EFGやその異方性がSOCで変わり得る

EFGは「絶対値の一致」を狙うより先に、(i) 計算条件を変えても値が動かない、(ii) 構造や磁性の変更に対するトレンドが物理的に説明できる、という“再現性”を確保することが大切である。

7.1.3 実験量への換算

NQR/NMR側では、核四重極モーメント $Q_n$ とEFGから四重極結合定数（あるいは四重極周波数）を組み立てる。代表的に

$${\nu }_Q\propto e{Q}_n{V}_{zz}$$

の形で周波数スケールが決まる（係数は核スピン $I$ と定義流儀に依存）。したがって比較の実務は、

1. DFTで $V_{zz}$ と $\eta$ を得る
2. 文献にある $Q_n$（核種固有）を用いて ${\nu }_Q$ に換算
3. 実験の共鳴周波数、分裂、線形状の温度依存（あるいは圧力依存）と突き合わせ
4. 一致／不一致の原因を「構造・秩序・揺らぎ」に分解する

という流れになる。

特に、電子の四重極秩序やネマティック転移がある場合には、

• $V_{zz}$ や $\eta$ の温度依存が異常を示す
• 線分裂が出る（等価サイトが不等価になる）
• 線幅が増大する（ゆらぎの増大、ドメイン形成） といった形で現れ得るため、DFTの静的EFGは「低温極限の基準」や「秩序相の候補比較」に強い。

7.1.4 静的DFTと揺らぎのギャップ

EFG計算の強みは、局所の対称性と電荷分布異方性を“元素・サイト選択的に”定量できる点にある。一方で、温度によるゆらぎや動的効果（phonon、電子揺らぎ、Kondo/valence fluctuation等）は、静的DFT単体では表現が難しい。

このギャップを埋める実務的アプローチとしては、

• 複数構造（対称性が異なる候補）でのEFG比較
• 圧力・ひずみ・原子変位のパラメトリックスキャン（擬似的にゆらぎの方向を探索）
• DFT+U/DMFT等で局在性や温度効果を取り込む（必要なら） が用いられる。

7.2 ひずみ応答

7.2.1 ひずみを入れる目的

磁歪や磁気弾性は、最終的には

• ひずみ $\epsilon$ に対して、電子状態（占有・混成・軌道異方性）がどう変わるか
• その変化がSOCを介して磁気異方性やエネルギーにどう反映されるか という“差分”の問題である。

第一原理計算は、同一のハミルトニアン（近似）に対して

• 構造（格子）
• 電子状態（DOS、占有、電荷密度）
• 磁性（スピン・軌道モーメント、磁化方向依存エネルギー） を一貫して出せるため、磁気弾性の微視的起源を追うのに適している。

7.2.2 ひずみの入れ方

四重極と磁気弾性は既約表現（対称チャネル）で整理すると見通しが良い。立方晶を例にすると、

• tetragonal型（${\mathrm{\Gamma }}_3$）：$(2{\epsilon }_{zz}-{\epsilon }_{xx}-{\epsilon }_{yy})$ や $({\epsilon }_{xx}-{\epsilon }_{yy})$
• shear型（${\mathrm{\Gamma }}_5$）：${\epsilon }_{xy},{\epsilon }_{yz},{\epsilon }_{zx}$ といった形で、四重極演算子 $O_{\mathrm{\Gamma }\gamma }$ と一対一に結び付く。

第一原理では、単に格子を伸ばす／縮めるだけでなく、

• “どの成分のひずみ”を入れたか
• 体積保存（pure shear）にするか、体積も変えるか
• 内部座標の緩和を許すか固定するか
  を明確にして、対称チャネルごとの応答を切り出すのが重要である。

7.2.3 四重極成分に相当する指標の選び方

「四重極子成分」を第一原理から抽出する方法は一つではない。目的に応じて、次のような指標を使い分ける。

• 電荷密度の異方性（局所テンソル量）
  原子周りの電荷密度を球面調和関数で分解し、$l=2$ 成分（四重極）を追う発想。
  実装はコード依存だが、EFGがその代表例である。

• 軌道占有の異方性（$d$軌道・$f$軌道の分解）
  例えば $d_{x^2-y^2}$ と $d_{3z^2-r^2}$、あるいは $t_{2g}/e_g$ の占有差など、対称性に沿った占有指標を定義する。
  これは「四重極（rank-2）＝軌道の二次異方性」という直観に対応する。

• スピン密度の非球対称性（磁気双極子項などに関連）
  XMCD和則に現れる $m_T$ のような量は、スピン密度の非球対称性を反映し、四重極的な成分と結び付く。
  第一原理側では、磁化方向を変えたときの局所スピン密度の異方性変化として追う発想になる。

実務上は、最初は「比較可能で堅牢」な指標（EFGや軌道占有差、MAE差分）を選び、必要に応じてより詳細なテンソル分解へ進むのが安全である。

7.2.4 SOCをどう扱うか

磁気異方性（MAE）や磁気弾性定数は、SOCを入れてはじめて意味を持つ量が多い。したがって、ひずみ応答で狙うべきは、概念的には

• SOCなし：ひずみで“軌道占有や結晶場”がどう変わるか（四重極の源）
• SOCあり：その変化が“エネルギー差”としてどう現れるか（磁気異方性・磁気弾性） という二段構えである。

具体的には、磁化方向を変えて全エネルギー差を取り、

$$\mathrm{MAE}(\epsilon )=E(\mathbf{M}\parallel {\hat{n}}_1;\epsilon )-E(\mathbf{M}\parallel {\hat{n}}_2;\epsilon )$$

を評価し、そのひずみ微分

$$\frac{d\mathrm{MAE}}{d\epsilon }$$

が磁気弾性（あるいは逆磁歪的な応答）に対応する、という戦略が取りやすい。

また、より直接に磁気弾性定数 $b_1,b_2$ を狙う場合は、対称ひずみを入れたときのエネルギー変化をフィットして係数を抜く（“計算版の磁歪測定”）ことになる。

7.2.5 「模型」と「第一原理」の役割

第6章の模型（結晶場モデル、多軌道Hubbard、RPA等）は、

• なぜ特定の対称チャネルの ${\chi }_{\mathrm{\Gamma }}$ が増大するのか
• 相互作用がどの形で効くと、どの秩序が優勢になるか といった“機構の見取り図”に強い。一方で第一原理は、
• 現実材料の結晶構造・元素置換・界面の効果を含めた符号と大きさ
• 元素別・サイト別寄与（どの原子のどの軌道が効くか） の裏付けに強い。

両者を接続する実務的な筋道は次の通りである。

1. 模型で「支配的チャネル（$\mathrm{\Gamma }$）」「増大の理由（相互作用の形）」を特定
2. 第一原理で、同じ対称チャネルのひずみ応答（占有異方性、MAEのひずみ微分など）を評価
3. 実験（超音波、二色性、RXS、NMR/NQR）で観測される“チャネル選択的な異常”と照合
4. 一致しない場合は、(a) 構造仮定、(b) 磁性仮定、(c) SOC/相関の不足、(d) 境界・ドメイン・不均一、に分解して仮説を更新

この往復が、四重極応答と磁気弾性を「定性的主張」ではなく「対称性＋定量」で押し切るための最短ルートになる。

小まとめ

EFGは「四重極（rank-2）情報を核位置で読む」量として、第一原理とNMR/NQRを最短距離で結ぶ。一方、ひずみ応答は「四重極成分（占有異方性）→SOC→エネルギー差」という差分構造を通じて、磁気弾性の電子状態起源を定量的に裏付ける。
模型は増大機構の見取り図、第一原理は現実材料での符号・大きさ・元素別寄与の裏付け、という役割分担を明確にし、対称チャネルを軸に実験と接続することで、議論が一段強固になる。

まとめと展望

四重極応答理論は、外場勾配・ひずみ・結晶場変調に共役な四重極子演算子と感受率で、超音波・散乱・NMR/NQRの観測量を統一して扱う枠組みである。一方、磁歪・磁気弾性は、ひずみ（rank-2）と磁化方向の二次形式（rank-2）が同じ対称性チャネルで結合するため、四重極子チャネルの言葉で自然に整理できる。二色性分光は、四重極子成分（異方性）を元素・軌道別に追跡でき、磁歪を「格子変形」ではなく「電子状態の再配列」として分解する実験設計を可能にする。

今後は、(i) 同一の対称性基底（$\mathrm{\Gamma }$）で超音波・RXS/RIXS・NMR/NQR・二色性を整合的に解析する枠組みの定着、(ii) 第一原理でEFGとひずみ応答を押さえつつ多体模型で ${\chi }_{\mathrm{\Gamma }}(\mathbf{q},\omega )$ の増大機構まで詰める統合、(iii) 磁気四重極モーメントを含む交差応答（磁気電気・熱応答）の整理と、磁気弾性との“設計変数”としての接続、が進むことで、四重極は「隠れた自由度の検出」から「材料設計の軸」へと重心を移していくと考えられる。

---

参考文献

• E. Callen and H. Callen, Phys. Rev. 139, A455 (1965).
• M. Shitade, Phys. Rev. B 99, 024404 (2019). DOI: 10.1103/PhysRevB.99.024404
• D. F. McMorrow and S. W. Lovesey, Solid State Communications 143, 367–372 (2007).
• D. Brinkmann, Prog. Nucl. Magn. Reson. Spectrosc. 48 (2006) 1–76.
• （国内資料）多極子入門、超音波と多極子の解説、KEK IMSSトピックス（2020）、SPring-8 BL25SU Outline など（各自参照）
• M. Shitade, “Theory of spin magnetic quadrupole moment and temperature-gradient-induced magnetization,” Phys. Rev. B 99, 024404 (2019). DOI: 10.1103/PhysRevB.99.024404
• KEK IMSS トピックス（2020）：四極子と磁気異方性に関する解説（題目・詳細は各自参照）
• SPring-8 BL25SU Beamline Outline（helicity switching等、運用情報）

（注）ユーザー提示の2025年文献・施設紹介等は、本稿では「位置づけの例」として扱い、書誌詳細の正確性は各自で原典確認する前提とした。