四重極応答理論

四重極応答理論は、双極子（電荷・磁気モーメント）では捉えきれない空間分布のゆがみ（2階のモーメント）に対する物質の応答を、対称性と線形応答で定量化する理論である。結晶場・歪み・共鳴X線散乱・NMR/NQRなどで現れる観測量を、四重極子演算子と感受率の枠組みで統一して扱える点に強みがある。

参考ドキュメント

1. 日本語の多極子入門資料（多極子・四重極子の定義と対称性分類の整理）
2. 日本語の超音波と多極子に関する解説資料（弾性定数と四重極感受率の関係の議論）
3. 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

1. 四重極子の定義

1.1 多極子展開における位置づけ

電荷密度 $\rho (\mathbf{r})$ の空間分布は多極子で特徴づけられる。全電荷（単極子）、双極子 $\mathbf{p}$、四重極子テンソル $Q_{ij}$ は、分布の重心、一次の偏り、二次のゆがみを表す量である。

連続分布の電気四重極子（トレースを除いた形）を

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

で定義する流儀が広く用いられる（係数やトレース処理は分野ごとに規約がある）。

結晶中の局所自由度（軌道・スピン・格子）を扱うときは、テンソル成分 $Q_{ij}$ よりも、点群の既約表現に沿って四重極子演算子 $O_{\mathrm{\Gamma }\gamma }$ を定義するほうが実用的である。たとえば立方対称では、$O_{xy}$ や $O_2^0$ などの基底を用いて、どの対称性チャネルが軟化するかを整理する。

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

四重極子は電気的な量に限らない。時間反転に対する性質で分類すると、

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

近年は、Kubo公式だけで書いた応答が不自然な振る舞いを示す場合に、（軌道・スピンの）四重極モーメント項を補正として差し引く必要があることが系統的に議論されている。

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

2.1 外場の勾配に結合するという性質

双極子が電場 $\mathbf{E}$（あるいは磁場 $\mathbf{B}$）に一次で結合するのに対し、四重極子は場の勾配に結合する。電気四重極子なら

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

のように書ける（係数は規約で変わる）。四重極応答とは、外場の勾配や歪み（結晶場の変形）に対して四重極子がどれだけ誘起されるか、あるいは四重極子同士がどれだけ相関するかを表す。

2.2 Kubo公式

四重極子演算子 $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$ では、四重極秩序のゆらぎ（相転移に近い増大）を反映する。

ただし、磁気電気応答や重力擬ポテンシャルを介した応答などでは、単純なKubo項に加えて、磁気四重極モーメントに由来する項を差し引いて物理的な応答へ整形する議論がある。

2.3 既約表現表示

結晶場中の$f$電子や$d$電子では、四重極子はテンソル成分 $Q_{ij}$ よりも、点群の既約表現で分類した $O_{\mathrm{\Gamma }\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 }}$ は結合定数である。このとき四重極感受率は

$${\chi }_{\mathrm{\Gamma }}(T)=-\frac{1}{g_{\mathrm{\Gamma }}}\frac{\partial ⟨O_{\mathrm{\Gamma }\gamma }⟩}{\partial {\epsilon }_{\mathrm{\Gamma }\gamma }}$$

として定義できる。

3.2 弾性定数への寄与

弾性定数 $C_{\mathrm{\Gamma }}$ は自由エネルギー $F$ の二階微分であり、

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

で与えられる。$H_{QS}$ を通じて電子系が歪みに応答すると、弾性定数に

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

の形で寄与が現れる。${\chi }_{\mathrm{\Gamma }}(T)$ の増大は $C_{\mathrm{\Gamma }}$ の低下（軟化）として観測されるため、超音波測定は四重極感受率を読み出す有力な方法である。

3.3 超音波測定で得られる情報

超音波は縦波・横波により異なる対称性の歪みを選べるため、どの $\mathrm{\Gamma }$ チャネルが支配的かを分離しやすい。強相関$f$電子系では、磁気秩序ではなく四重極秩序（反強四重極など）が相転移の主役となる例があり、弾性異常が重要な手がかりになる。

4. 散乱で見る四重極

4.1 共鳴X線散乱とE2遷移

共鳴X線散乱（resonant X-ray scattering; RXS）では、吸収端近傍で散乱振幅が電子遷移行列要素により増強される。電気双極子（E1）遷移だけでなく、電気四重極子（E2）遷移が寄与する場合がある。E1とE2の干渉や、共鳴項とトムソン散乱の干渉を測ることで、四重極成分に敏感な散乱機構を切り分ける。

4.2 禁制反射と多極子秩序

結晶学的に禁制な反射が、共鳴条件下で現れることがある。これは散乱振幅が単純な電荷密度以外のテンソル成分（四重極など）を含むためである。反強四重極秩序では、方位角依存性・偏光解析・エネルギー依存性により四重極配置を決定する戦略が取られる。

4.3 国内情報源としての放射光施設発信

国内では、SPring-8などの放射光施設が、共鳴散乱や共鳴非弾性散乱（RIXS）を用いた研究成果・利用報告・プレスリリースを公開している。装置条件や測定条件に踏み込んだ資料が得られるため、実験設計の具体像を掴みやすい。

5. 局所プローブで見る四重極

5.1 核四重極相互作用

核スピン $I\ge 1$ の核は核四重極モーメント $Q_n$ を持ち、核位置の電場勾配（electric field gradient; 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 =(V_{xx}-V_{yy})/V_{zz}$ は非対称度であり、主軸系で $|V_{zz}|\ge |V_{yy}|\ge |V_{xx}|$ を採用する。

5.2 多極子秩序との接続

電子系が四重極秩序を形成すると、局所電荷分布の対称性が変わり、EFGが変化する。その結果、NQR/NMRで線分裂・線幅変化・共鳴周波数シフトとして現れる。共鳴散乱とNMR/NQRを組み合わせて、対称性と局所情報を相補的に押さえる研究も行われている。

6. 理論モデル

6.1 結晶場準位と四重極感受率

局在$f$電子では結晶場により多重項が分裂し、基底状態が四重極子に対して縮退や近縮退を持つと、低温で ${\chi }_{\mathrm{\Gamma }}(T)$ が $1/T$ 的に増大しやすい。弾性定数の温度依存が「定数項＋キュリー項」の形をとる場合、この単イオン由来の四重極応答がまず検討される。

6.2 多軌道模型と行列感受率

遍歴性を含む模型（多軌道Hubbardなど）では、裸感受率 ${\hat{\chi }}^{(0)}(\mathbf{q},\omega )$ と相互作用行列 $\hat{U}$ を用いて

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

の形で感受率を構成し、どの対称性チャネルが増大するかを評価する。多極子チャネルは、軌道自由度と対称性の組で特徴づけられるため、基底選択と群論の整理が実装上の鍵になる。

6.3 高次多極子との結合

四重極秩序が段階的に現れる系では、四重極—八極子などの結合が実効相互作用を作るという見方が重要になる。超音波・散乱・熱力学量を組み合わせて、どの多極子が主役であるかを同定する議論が進められている。

7. 近年の展開：四重極モーメントと輸送・熱応答

7.1 熱力学的四重極モーメント

自由エネルギーの外場勾配に対する応答として四重極モーメントを定義し、温度勾配下の分極（thermopolarization）などを記述する枠組みが提案されている。四重極の空間変化が境界分極に寄与するという整理は、バルクと境界の応答を同時に扱ううえで有用である。

7.2 磁気電気応答と磁気四重極モーメント

スピン磁気四重極モーメント（spin MQM）や軌道磁気四重極モーメント（orbital MQM）は、温度勾配誘起磁化や重力擬磁気電気応答の式に補正項として現れる。Kubo項のみで得られる見かけの振る舞いを、四重極モーメントを用いて物理的な応答へ整形するという論理が確立されつつある。

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

8.1 EFG計算とNMR/NQR

NMR/NQRで必要なのはEFGの主値 $V_{zz}$ と非対称度 $\eta$ である。DFTでEFGを計算し、四重極結合定数へ換算して実験と比較する手法は広く用いられてきた。基底関数や交換相関汎関数の選択で定量値が変わることがあるため、比較の前提条件を揃える必要がある。

8.2 模型と第一原理の役割分担

• 模型（結晶場＋多体相互作用）：${\chi }_{\mathrm{\Gamma }}(T)$ や秩序パラメータの対称性、相互作用機構の同定に向く
• 第一原理：局所対称性、EFG、軌道占有、歪みに対する電子状態変化の裏付けに向く
• 散乱断面積：RXS/RIXSでは遷移行列要素（E1/E2）、偏光、幾何因子が効くため、実験配置に即した理論が必要になる

9. 手法ごとの四重極情報の得られ方

手法	主に結合する量	観測される量	四重極応答理論での対応
超音波（弾性定数）	歪み ${\epsilon }_{\mathrm{\Gamma }\gamma }$	$C_{\mathrm{\Gamma }}(T)$ の軟化・異常	$C_{\mathrm{\Gamma }}=C_{\mathrm{\Gamma }}^{(0)}-g_{\mathrm{\Gamma }}^2{\chi }_{\mathrm{\Gamma }}$
RXS（共鳴X線散乱）	電子遷移（E1/E2）	禁制反射・方位角依存・偏光依存	多極子テンソル成分、E2寄与・干渉
RIXS（共鳴非弾性）	局所励起・電荷/軌道/スピン励起	エネルギー損失スペクトル	励起対称性の分解、選択則に多極子が入る
NMR/NQR	EFG $V_{ij}$ と核四重極モーメント	分裂、周波数、線幅	電荷分布対称性変化によるEFG変化

10. 注意点

10.1 規約の違い

四重極テンソル $Q_{ij}$ の定義は係数やトレース除去の扱いが分野で揺れる。散乱理論（テンソル原子因子）、局所模型（$O_{\mathrm{\Gamma }\gamma }$）、NMR（$V_{zz},\eta$）を横断するときは、採用した規約を明示する必要がある。

10.2 静的応答と動的応答

超音波はMHz帯であり静的極限に近い情報を与えることが多い一方、RIXSはeVスケールの動的構造因子を読む。$\chi (\omega )$ のどの周波数領域を議論しているかを揃えないと、同じ四重極応答という言葉でも意味がずれる。

10.3 駆動要因の切り分け

弾性異常が見えても、結晶場準位の近接（単イオン効果）、相互作用に由来する臨界ゆらぎ、格子不安定（フォノン）などが関与しうる。複数手法で対称性・局所量・相関長を相補的に押さえるのが有効である。

まとめと展望

四重極応答理論は、歪み（超音波）・共鳴散乱（E2寄与）・局所電場勾配（NMR/NQR）という異なる実験量を、四重極子演算子と感受率で統一して扱う理論である。強相関系の多極子秩序に加え、磁気四重極モーメントを介した磁気電気・熱応答の再定式化が進み、四重極は秩序パラメータとしてだけでなく輸送・熱応答の基礎量としても重要性を増している。

今後は、(i) 超音波・RXS/RIXS・NMR/NQRを同一の対称性基底で整合的に解析する枠組み、(ii) 第一原理でEFGや歪み応答を押さえつつ多体模型で ${\chi }_{\mathrm{\Gamma }}(\mathbf{q},\omega )$ の増大機構まで詰める統合、(iii) トポロジカル物質や非相反系で四重極モーメントが境界応答へ与える寄与の一般化、が進むことで、四重極応答理論は隠れた自由度の検出から設計変数としての多極子へと重心を移していくと考えられる。

その他参考にしたsources

• D. F. McMorrow and S. W. Lovesey, “Resonant X-ray scattering and the determination of multipole order,” Solid State Communications 143, 367–372 (2007).
• K. Matsumura et al., electric quadrupole order identified by RXS（arXiv:cond-mat/0403535 の系統）。
• T. Arima et al., orbital/gravitomagnetoelectric response と orbital MQM補正（Phys. Rev. B 107, 214109 (2023) などの系統）。DOI: 10.1103/PhysRevB.107.214109
• thermodynamic quadrupole moment と thermopolarization の理論提案（arXiv の関連稿の系統）。
• D. Brinkmann, “Nuclear quadrupole interactions in solids,” Progress in Nuclear Magnetic Resonance Spectroscopy 48 (2006) 1–76 の系統。
• EFGの第一原理評価と四重極結合定数（Phys. Rev. B 70, 115111 (2004) などの系統）。DOI: 10.1103/PhysRevB.70.115111