《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

NMR で観測する重い電子系の電子状態
兵庫県立大学水戸毅*1

1 はじめに重い電子系とは，一般には，局在性の強い f 電子が伝導電子と混成し，電子の有効質量が大きく増大した状態が観測される物質系を指す．一見，f 電子を有するランタノイド系やアクチノイド系化合物に特化した狭い領域の物理現象のように思われるかも知れないが，この分野に多くの研究者が魅了される理由の一つは，重い電子系やその周辺領域に多様な物性が観測されるからであろう．それらの物性が現れる機構の底流にある伝導電子と f 電子の混成の度合いを，温度・圧力・磁場・元素置換などの外的要因によって比較的容易に制御することができることも，この物質系の大きな魅力である．一般に，f 電子の局在性が強い領域では磁気秩序状態が生じるが，多くの場合は金属のままである．混成が増して局在性が弱まると有効質量が増大した非磁性のフェルミ液体状態が現れるが，この磁性−非磁性転移近傍では非フェルミ
液体的な振舞いが観測され，非従来的な超伝導が出現することもある．また，遍歴性が強い領域では，f電子の混成バンドへの出入りによってランタノイドイオン，或いはアクチノイドイオンの平均価数が整数からずれる価数揺動状態になるが，価数揺動を示す物質の中には低温で半導体的に振舞い，近藤絶縁体と呼ばれる物質群もある，その他，f 電子特有の大きな軌道角運動量に起因する多極子秩序を示す物質も知られる．こうした多彩な物性の本質を解明するために，様々な実験研究が行われるわけであるが，その中で核磁気共鳴（NMR）や核四重極共鳴（NQR）測定による研究は独自の重要な情報を与える．しかも，微視的な実験手法とみなされるものの中で，NMR や NQR はそれほど大掛かりな実験装置を要しないた

めに手軽に研究がなされ，自ら実験を行わなくとも，それらのデータを目にすることはよくあるだろう．
本講義では，NMR や NQR の実験データそのものから，或いはそれらを他の実験データと比較しながら，重い電子系の特性を紹介していきたい．しかし，本講義を受ける方々の全てが NMR に詳しいわけではないだろうから，講義では NMR や NQR の解説も並行して行い，このテキストはむしろそのために使っていただきたい．ただし，NMR の原理と実験データから得られる静的な情報について記している

うちに既定のページ数に達してしまったために，動的な情報（つまり，T1 などの緩和現象）について紙面を割くことができなくなってしまった．これらについては別途講義で解説する．先に記した通り，ひと口に重い電子系と言っても，その周辺の現象も含めると相当に広い領域である．そこで，特に超伝導に関してはこれまでにも多く優れた解説があり、超伝導出現によって重い電子系としての性質がマスクされてしまう側面もあるので、本講義ではそれ以外の物性に重点を置こうと思う。

*1 E-mail: mito@sci.u-hyogo.ac.jp

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）)'$!%*#$)&'

µ!

µ"
+!

−!

(%&$&'!"!#$%&'
図1

1

H における電子と原子核の模式図．

2 共鳴の条件
電子が持つ磁気モーメント µe と磁場H0 との相互作用は

HZeeman,e = −µe · H0であるが，この µe はスピン角運動量 !s （! はプランク定数 h の 1/2π）を用いて，µe = −gµB s= −γe !sと表わされる．ここで，g は電子の g 因子，µB はボーア磁子，γe は電子の磁気回転比である．同様に，原子核も角運動量 !I を持つことから磁気モーメント µN を有し，H0 との間に

HZeeman,N = −µN · H 0= −gN µN I · H 0= −γN !I · H 0(1)の相互作用がはたらく．ここで，gN は核の g 因子，µN は核磁子，γN は核スピンの磁気回転比である．

µB =

e!
e!
M
，µN =
で，核と電子の質量比∼ 1800 より，HZeeman,N は HZeeman,e に比べて 3 桁
2mc
2M c
m

ほど小さいことがわかる．

式 (1) において，磁場の方向を z に選び，H 0 = (0, 0, H0 ) とすると

HZeeman,N = −γN !H0 Iz ．
このエネルギー固有値 EZ,m = ⟨m′ |HZeeman,n |m⟩ は，m = I, I − 1, ..., −I として，m′ = m のときのみ値をとり，EZ,m = −γN !H0 ⟨m|Iz |m⟩= −γN !H0 m．

(2)

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

m=−

3
2

m=−

1
2

m=+

1
2

m=+

3
2

∆ E m = γhH 0

!!"!"

!!"≠"

3
E−3 / 2 = γhH 021E−1/ 2 = γhH 02
1
E1/ 2 = − γhH 023E3 / 2 = − γhH 02

図 2 ゼーマン相互作用による核スピンのエネルギー準位（I = 3/2 の場合）．左：H0 = 0 のとき．
右：H0 ̸= 0 のとき．

共鳴は準位間の遷移を観測するので，隣接するエネルギー準位間の差を求めると，∆EZ,m = EZ,m − EZ,m+1= γN !H0
γN
=
hH0 ．
2π
図 2 に，I = 3/2 の場合の例（7 Li，8 Be，11 B，23 Na，35 Cl，63,65 Cu，75 As，等がこれにあたる）を示
す．エネルギー準位は均等のエネルギー差で 4 つに分裂し，∆EZ,m = hν0 を満たす周波数 ν0 を持つ電磁波を照射することによって遷移を起こすことができる．つまり，共鳴の条件は，

ν0 =

γN
H0
2π

(3)

で与えられる．この式は，NMR において磁場 ⇔ 周波数の変換にも頻繁に使われる．

3 NMR/NQR スペクトル
3.1 シフト式 (3) によって得られる共鳴条件は，原子核が孤立して存在しているときに相当する．実際の物質中

では，核の周囲にある電子系が外部磁場によって分極し，余分な局所的内部磁場 ∆Hint を核位置につくるために，共鳴は式 (3) の条件からずれて生じる．外部磁場 H によって誘起された ∆Hint の割り合い∆Hint /H をシフトと呼ぶ．∆Hint は様々な機構によって生じ，ここではそれらについて詳しくは述べな

いが，一般的には，閉殻電子が外部磁場を遮蔽する効果や他の核スピンによる双極子磁場などが考えられる．一方，ランタノイドやアクチノイド元素など磁性が強い元素を含む化合物では，上記の機構よりも，伝導電子や局在磁気モーメントが原因となるシフト（ナイトシフトと呼ばれる）が支配的になる．

図 3(a) と (b) に，観測核位置に ∆Hint が生じる機構の模式図を示す．重い電子系化合物の研究においては，図 3(a) のように，磁性原子つまりランタノイドやアクチノイド原子を直接観測した例は，これらのサイトが非磁性状態に近い場合や磁気秩序状態にある場合などを除いて，非常に少ない．これは，f 電子が有する大きな角運動量に由来してナイトシフトが非常に大きくなり，共鳴磁場、或いは共鳴周波数を見出しにくかったり，局在性の強い f 電子の揺らぎによって緩和が非常に短くなるためであると思われる．従って，f 電子系の NMR 測定は，図 3(b) にあるように，磁性原子近傍の非磁性サイトを観測する方が圧倒的に多いのが現状である．どちらの場合であっても，磁場によって分極してできた磁気モーメン
ト ⟨M ⟩（ここでは ⟨ ⟩ は時間平均を表す）と，それが原因で観測核位置に生じる内部磁場 ⟨∆Hint ⟩ との

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

!"#

!$#

〈"〉
〈"〉
〈∆!$%#〉
!!"#

!0

!!"#

〈∆!$%#〉!"#$%"&'()*+%,'+-'!./!"#$%"&'()*+%,'+-'!./

図 3 磁気モーメント ⟨M ⟩ が観測核位置に内部磁場 ⟨∆Hint ⟩ を生じる様子を表す模式図．(a) 磁性
元素を直接観測している場合．(b) 磁性元素に隣接する非磁性元素を観測している場合．

間には，⟨∆Hint ⟩ = A⟨M ⟩の関係があることが容易に予想される．このとき，両辺を印加している磁場 Hext で割れば，χ を帯磁率として，⟨∆Hint ⟩
Hext
⟨M ⟩
=A
Hext
= Aχ

K≡

(4)

(5)

となり，ナイトシフトは基本的に帯磁率に比例する量であるということが分かる．ただし，ここで χ は

局所的な帯磁率に相当すること，また，帯磁率には Curie-Weiss 項，van Vleck 項，Pauli 常磁性項などがあるように，ナイトシフトも複数の帯磁率成分を含み，それぞれの項に対応する結合定数 A を持つことに言及しておく．

以下では，先の述べた核位置における局所内部磁場 ⟨∆Hint ⟩ に対して，シフト（ナイトシフト）の符号について述べる．核スピンは，核位置での磁場の総量 |⟨Htotal ⟩| に対して式 (3) に従う．つまり，γN|⟨Htotal ⟩|2πγN|Hres + ⟨∆Hint ⟩|．
=2π

νres =(6)

νres とH res はそれぞれ共鳴周波数と共鳴磁場である．実際に NMR スペクトルを測定するには，測定周波数を固定して磁場を掃引する場合と，磁場を固定して周波数を掃引する場合がある．

まず，磁場掃引の場合の模式図を図 4 に示す．固定周波数は ν0 であり，縦軸は磁場を表す．今，

簡単のために Hres と ⟨∆Hint ⟩ が平行であるとし，Hres の向きを正にとる．図 4(a) はシフトがゼロ（⟨∆Hint ⟩ = 0）のときを表し，このとき Hres = H0 とする（式 (3) 参照）．図 4(b) の場合は，⟨∆Hint ⟩ > 0である．このとき，式 (6) を満たして共鳴を起こすためには Hres < H0 となり，式 (4) よりK=

H0 − Hres
Hres

(7)

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

(a)

(b)

(c)
K>0

K<0

H
〈∆Hint〉〈∆Hint〉
H0

Hres
Hres

ν0 =

γN
H0
2π

ν0 =

γN
H res + ∆H int
2π() ν = γ2π (H − ∆H )
N

0

res

∆!

∆!

H

H
Hres H0=2πν!/γN
図4

(a)

int

H0

Hres

磁場掃引の測定によるシフトの符号．

(b)

(c)
K>0

K<0

H
〈∆Hint〉〈∆Hint〉(H + ∆H )

H0

0

ν0 =

int

= 2πν res / γ N

H0H0(H − ∆H )0int= 2πν res / γ N

γN
H0
2π

∆ν

∆ν

ν
ν!

ν"#$

ν
ν"#$ ν!

図 5 周波数掃引の測定によるシフトの符号．

は正のシフトを与える．図 4(c) の場合は，Hres と ⟨∆Hint ⟩ の向きが反平行であるので，共鳴を起こす
ためには Hres > H0 となり，同様に式 (7) を用いてシフトは負の値になる．

次に，磁場を固定して周波数を掃引する場合の模式図を図 5 に示す．同様に縦軸は磁場で，固定磁場H 0 の向きが正である．図 5(a) はシフトがゼロのときに相当しており，このときは式 (3) が成り立つ．

図 5(b) は，一定磁場 H0 の印加によって核位置に正の向きに ⟨∆Hint ⟩ が誘起される状態で，⟨H total ⟩ = H 0 + ⟨∆Hint ⟩．

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

式 (6) より

νres =

γN|H 0 + ⟨∆Hint ⟩|2πであるから，νres > ν0 となる．シフトは， 式 (4) と (3) より

∆Hint
H0
2πνres /γN − H0=
H0
νres − ν0
=
．
ν0

K=

(8)

図 5(c) は，核位置に負の向きに ⟨∆Hint ⟩ が誘起された状態で，νres < ν0 であるから式 (8) よりシフトは負になる．

以下に，ナイトシフト測定の具体例を示す．多くの重い電子系化合物は，高温で f 電子が局在化し，帯磁率は Curie-Weiss 則に従う．一方，低温では近藤効果によって f 電子の局在モーメントが非磁性化され，帯磁率は温度変化を示さなくなる（Fermi 液体状態）．通常，これら二つの状態間はクロスオーバー

で移り変わるが，ここで示す YbInCu4 は，TV = 42 K での一次相転移によって，高温側の局在状態と低温側の Fermi 液体状態の間を不連続に変化するという特徴がある．TV は価数転移温度で，高温相でのYb 価数はほぼ +3 価，低温相では約 2.84 価と見積もられており [1]，上記の両極の f 電子状態を一度に見ることができる都合のよい物質である．

図 6(a) に，この物質の温度 − 圧力相図を示す．価数転移は Pc ≈ 2.4 GPa の圧力で消失し，それより高い圧力下での基底状態は強磁性秩序状態である [2, 3]．図 6(b) に，常圧下高温相における 115 In-NMR
スペクトルの温度依存性を示す．115 In は核スピン I = 9/2 であるが，YbInCu4 の In（4c）サイトは対称性が高いために（内挿図），後述の核四重極相互作用によるスペクトルの分裂は生じない．この測定は磁場掃引で行っているから図 4 の場合に相当し，測定周波数が 65.307 MHz，115 In の γN /2π の値が9.3295 MHz/T であるので，図 4(a) に示すようにもし K = 0 であれば，共鳴磁場は H0 ≈ 7.000 T である．実際は，スペクトルが H0 よりも高磁場側で観測されているから，図 4(c) より K < 0 であり，低温になるに従ってスペクトルがより高磁場側にシフトしていることは，シフトが負の方向に増大していることを意味している．式 (7) を用いて見積もったシフトの温度依存性を図 4(c) に示す．高温相における Curie-Weiss 的な温度変化に対し，TV 以下では絶対値が急激に減少し温度依存性が殆ど無い．これは，先に記した通り，f 電子の局在性（高温）と遍歴性（低温）に起因した振舞いである．高温相における局在性の程度については以下の様な考察ができる．式 (5) より，K と帯磁率 χ には線形の関係が期待されるので，図 4(d) のように温度を陰関数として K 対 χ のプロットをすると，両者間の結合定数 Ahf（超微細結合定数）が得られる．この結合定数が圧力によって変化しないと仮定すると，高圧下におけるK から χ を見積もることができる．図 4(e) には，そうして見積もった χNMR の逆数を温度に対してプ
ロットしてある．この結果から見積もった，YbInCu4 の高温相における µeﬀ は，Yb3+ 状態に期待される µeﬀ = 4.53µB に近く，Yb モーメントは良く局在した状況，つまり Yb 価数は 3+ に近い値にあると結論される．つまり，YbInCu4 では Yb3+ に近い状態を有する高温相が，そのまま圧力によって低温領域まで拡大し，Pc より高圧では，高温相の局在モーメントが強磁性的に秩序する状態の方が遍歴性の強い低温相よりも安定になる描像が明らかになった．

実験技術的な観点からは，汎用の SQUID（superconducting quantum interference device）装置では，試料スペースの制限等から 2 GPa 以上の高圧下で常磁性帯磁率を測定することは簡単ではない現状がある．一方，NMR 測定に関しては，測定核種等に制限はあるものの，最近は 2 GPa 以上の高圧下で

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）!

ac susceptibility
present dc magnetization

TM

TM

(c)

-8

YbInCu4115

In-NMR-6
0.1 MPa
2.0 GPa
2.6 GPaTV (0.1MPa)-4!

Knight Shift (%)

Resistivity
NQR
present dc magnetiztion

TV

10

(a)

!

Temperature [K]

YbInCu4TV!

TV (2.0GPa)

100

-2
0
0

50

100

150

200

Temperature (K)
1
2

3

5

10

15

20

25

0.0

30

Pressure [GPa]!

YbInCu4
In NMR

211 K

115

NMR Intensity (arb. units)

198 K

65.307 MHz

Yb

180 K
160 K

-2.0
0.00

X
(mol/emu)!

68 K

YbXCu458 K

0.03

YbInCu4115

In-NMR

(e)
50

NMR

88 K

45 K
7.10

0.02

100

109 K

7.05

0.01

χ (emu/mol)
!

127 K

7.00

-1.5

Cu

146 K

(b)

-1.0

0.1 MPa
2.0 GPa
2.6 GPa

1/χ

K=0

(d)
-0.5

!

1

Knight shift (%)

0

0
7.15

Magnetic field (T)

7.20

0

50

100

150

200

Temperature (K)

図 6 YbInCu4 における 115 In-NMR 測定の結果と相図 [3]．(a) 電気抵抗，NQR，交流／直流帯磁率測定によって求められた温度 − 圧力相図．図中の実線（価数転移）と点線（強磁性転移）はそれ
ぞれ一次と二次の相転移を，斜線の領域は二つの相が共存していることをそれぞれ意味する．(b) 常
圧下における 115 In-NMR スペクトルの温度変化．内挿図は YbXCu4 の結晶構造．(c) 115 In-NMR
スペクトルから見積もったナイトシフトの温度依存性（常圧，2.0 GPa，2.6 GPa）．(d) 帯磁率 vsナイトシフトのプロット（常圧）．(e) (d) から見積もられた超微細結合定数（−3.4 kOe/µB ）を用いて求めた帯磁率 χNMR の逆数（常圧，2.0 GPa，2.6 GPa）．直線と点線は，それぞれ µeﬀ = 4.5µBと 4.0µB を仮定した Curie-Weiss 則を表す．

あっても多くの測定がなされており，上記の様に高圧下における帯磁率の情報を得ることができるようになってきた．

YbInCu4 の帯磁率について説明したので，Yb 系化合物における局在−非局在の様子を知るために，
YbXCu4 (X = Cu, Ag) の帯磁率との比較をすることは参考になるであろう．これらの YbXCu4 物質
は全て同じ結晶構造（AuB5 型（C15b）構造．図 6(b) 内挿図参照）を持つが，X 元素によって遍歴性の強い物質から局在性の強いものまで幅広く変化する [5]．先述の通り，X = In の高温相は局在性の強い性質を示すが，X = Cu と Ag はそれぞれ電子比熱係数が γ = 550 mJ/molK2 [6] と 240 mJ/molK2 [7]の重い電子系化合物である．図 7(a) はこれら 3 つの物質の帯磁率を比較したもので，X = In と Cu の帯
磁率 χNMR は，先に記したのと同様，ナイトシフト測定から見積もったものである．X = Ag の χ の温度依存性は，重い電子系化合物，或いは中間価数物質と呼ばれるものに良く見られる振舞いである．温度

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

!

!

0.15

5

0.04
!

(emu/mol)
NMR

0.05 YbCu

(b)

0.03

YbCu 5
2.4 GPa
0.1 MPaχ(emu/mol)

0.10

0.05

(a)

2.5GPa2.0GPa0.1MPaχ, χ

NMR

YbInCu4

0.02

YbAgCu4
0.00
0

50

100

150

Temperature (K)

0

20

40

60

80

100

Temperature (K)

図 7 (a) YbXCu4 (X = In, Cu, Ag) の帯磁率の温度依存性．X = In と Cu の帯磁率 χNMR はそれぞれ 115 In-NMR，63 Cu-NMR によって測定されたナイトシフトから見積もられた．X = Ag の
χ は直流帯磁率測定による．(b) X = Cu の χNMR における圧力依存性．文献 [4] より引用．

減少と共に，高温の Curie-Weiss 則からずれて最大値を示すのは近藤効果の発達によるものと考えられ，
低温でも大きな状態密度を反映して通常金属よりも高い値で一定値となる．X = Ag → Cu と変化する

ことによって χ が極大を示す温度が低温側にシフトし，低温の χ の値が増大するのは，明らかに f 電子の局在性が強まったことを意味しており，そのさらに延長上に X = In があると考えられる．また，Yb系化合物への圧力印加が一般に f 電子を局在化する方向にはたらくことは，図 7(b) に示された X = Cu
における圧力依存性を見ても明らかである．同様な圧力効果は Sm 系でも観測され，Ce 系とは逆である
ことは良く知られる．

3.2 シフトの異方性
式 (4) では，ナイトシフトを単なるスカラーとして表式しているが，実際の固体中ではシフトは次式のようにテンソルとして表される．

!

⟨∆H⟩x
⟨∆H⟩y
⟨∆H⟩z

"

=

!

K11
K12
K13

K12
K22
K23

K13
K23
K33

"! "
Hx
Hy
Hz

(9)

ここで，シフトテンソルについて相反定理（Kij = Kji ）を用いた．シフトテンソルは，観測している核
位置での点群によって多くの場合より簡単に表現される．

例えば，重い電子系化合物の中でも最も有名な結晶構造の一つである ThCr2 Si2 型結晶構造（No. 139,

空間群 I4/mmm）の場合，2a サイト（位置対称 4/mmm），4d サイト（4̄m2），4e サイト（4mm）には，4 回（4d サイトは回反）軸とその軸を含む鏡映面が存在する．4 回軸を主軸（z 軸）とし，対称操作

xy → −xy ，yz → −yz ，xz → −xz や x2 → y 2 の考察から，K12 = K23 = K13 = 0 や K11 = K22 が#は得られるので，式 (9) 中のシフトテンソル K
# =
K

!

K11
0
0

0
K11
0

0
0
K33

"

(10)

と書けることがわかる．このような軸対称性がある場合，図 9 のように z 軸とのなす角が θ となる外部

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）!a!"#$%&'$%()"d!",$-.$-,)"e!*+)a=b∫c
図 8 ThCr2 Si2 型結晶構造（No. 139, 空間群 I4/mmm）．CeCu2 Si2 ，YbRh2 Si2 ，URu2 Si2 などがこの結晶構造を成す．

図 9 4 回対称性があるときの外部磁場

ex 方向の内部磁場 ⟨∆int ⟩(θ) と主軸方向成分 ⟨∆int ⟩x,y,z ．

磁場 Hex = (Hex sin θ, 0, Hex cos θ) を印加すると，⟨∆Hint ⟩x = K11 Hex sin θ⟨∆Hint ⟩z = K33 Hex cos θ⟨∆Hint ⟩(θ)
Hex⟨∆Hint ⟩(θ) = ⟨∆Hint ⟩z cos θ + ⟨∆Hint ⟩x sin θK(θ) =などを使って，磁場の方向のナイトシフト K(θ) は，

K(θ) = K11 sin2 θ + K33 cos2 θ(11)となる．従って，単結晶試料を測定する場合は，θ を 0◦ から 90◦ まで変化させると，式 (7) と (8) に従っ《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）(a)

K (0°) = K11

K (θ )
H ex

(b)
°
9080
°
70°

c

60°
50°
40°

θ

30°
20°
10°

K (90°) = K 33
H 0 /(1 + K11 )H01°

H 0 /(1 + K 33 )図 10 (a) 外部磁場方向に対する多結晶試料，または粉末試料中の微結晶（図中の四角柱）の向きを示す模式図．微結晶が軸対称性を有する場合は，磁場と対称軸とのなす角 θ が同一の微結晶からの
NMR 信号は同じ共鳴磁場（周波数）で観測される．(b) θ = 0◦ ∼ 90◦ で様々な値を成す微結晶群からの NMR 信号．信号強度は sin θ に比例し，K11 > 0 と K33 < 0 と仮定．

て K(θ) は

θ = 0◦
ωres = ω0 (1 + K33 )Hres = H0 /(1 + K33 )から

θ = 90◦ωres = ω0 (1 + K11 )Hres = H0 /(1 + K11 )まで連続的に変化する．
次に，式 (11) を多結晶試料，または粉末試料の測定に適用する．このときシフトは θ にしか依らないので，試料を形成する微結晶の中で，主軸と磁場のなす角が θ であるものは全て同じ共鳴磁場（或いは共鳴周波数）で信号が観測される．従って，微結晶が完全にランダムな方向を向いているとみなせる状況

では，図 10(a) に示すような単位球を考えて，θ 方向の信号強度は相対的に sin θ に比例することになる．
図 10(b) に示すように，θ を 0 から π/2 まで変化させたときの信号を全て足し合わせると粉末パターンと呼ばれる形状を持つ信号が観測される（図 11）．このように，NMR 測定ではたとえ粉末試料であっても，その粉末パターンを解析することによって，K11 と K33 成分を抽出することが可能である．

図 12 は，重い電子系として最初の超伝導体として有名な CeCu2 Si2 [8] の 29 Si-NMR スペクトルである [9]．測定には粉末試料を用いており，それをパラフィンワックスで固めて得た図 12(a) のスペクトルは，先に図 11 で示したものと同様の粉末パターンを示している．図 12(b) にはパラフィンワックスで固めていない試料で得られたスペクトルが示されているが，c 軸に平行なシフト成分 K∥ での信号強度が増大している．これは帯磁率の c 軸に平行な成分が垂直成分よりも大きいために，印加磁場によって自然に配向した（微結晶が c 軸 ∥ 磁場となるように向きを変えた）ためである．さらに人為的に小さな機械的振動を与えると，粉末試料はほぼ完全に配向し，図 12(c) に示したほぼ K∥ 成分のみのスペクトルにな《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

π /2

2∑ 2π sin θ !(θ )θ =0! (θ )

H 0 /(1 + K11 )H0

H 0 /(1 + K 33 )図 11 図 10(b) の各信号を足し合わせて得た，シフトの異方性によって生じる粉末パターン．

図 12 粉末 CeCu2 Si2 試料の 29 Si-NMR スペクトル [9]．29 Si の核スピンは I = 1/2．(a) 粉末試料が配向しないようにパラフィンワックスで固めて 1.3 K で測定．K⊥ ，K∥ はそれぞれ対称軸である c 軸に垂直なナイトシフト成分（= K11 ）と平行（= K33 ）成分．(b) と (c) はそれぞれ 1.0 K と
0.39 K で測定されたスペクトルで，共にパラフィンワックスは用いていないため，一部または全体が配向している．

る．このような磁場中配向は，粉末試料の NMR 測定ではしばしば用いられる手法で，c 軸に垂直な成分（今の場合は ab 成分）にランダムネスは残るが，シフトの一成分を得るには有効な手法である．配向した状態をパラフィンワックス等で固めて維持すれば，それに垂直な成分を測定することも可能である．

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

3.3 核四重極相互作用
3.3.1

ゼロ磁場下におけるスペクトル

ここまでは，主に核と周囲の環境との磁気的な相互作用について述べた．この節では，電気的な相互作用（核四重極相互作用）について述べるが，核四重極共鳴（NQR）の原理を理解するためだけではなく，多くの NMR スペクトルがこの相互作用の影響を受けた形状を示すことからも，避けるわけにはいかな

い内容である．核四重極相互作用とは，一般にはあまり耳慣れない用語かも知れないが，物質中に存在する結晶電場の原理を思い浮かべると理解しやすいのではないだろうか．結晶電場が，着目している電子と周囲の電荷とのクーロン相互作用に基づいたものであるのと同様に，核と周囲の電荷とのクーロン相互作
用を考え，それを原点（核位置）の周辺で展開（多重極展開）したときの四重極項がそれに当たる．

核外の電荷が核の電荷が分布する領域につくる静電ポテンシャルを φ(r ) とする．核の大きさは周囲の電荷分布の広がりに比べて十分小さいので φ(r ) を展開して

%
%
3
1 $
∂φ %%
∂ 2 φ %%
φ(r) = φ(0) +
xi
+
xi xj
+ ···
∂xi %r=0 2! i,j=1∂xi ∂xj %r=0
i=13$

ここで，x1 , x2 , x3 = x, y, z である．従って，核内の電荷密度を ρ(r ) とすると，静電エネルギーはE = φ(0)

&

%
3
$
∂φ %%
ρ(r)dτ +
∂xi %

r=0

i=1

&

%
&
3
1 $ ∂ 2 φ %%xi ρ(r)dτ +xi xj ρ(r)dτ + · · ·2! i,j=1 ∂xi ∂xj %r=0と表わされる．この内，実質的に核のエネルギー準位に影響を与えるのは第 3 項の核四重極相互作用で，ウィグナー-エッカルトの定理を使うなどして，最終的に次式で表わされる．

HQ =

e2 qQ
2
2
+ I−
)}．
{(3Iz2 − I 2 ) + η(I+
4I(2I − 1)(12)%
∂ 2 φ %%
ここで，Q は核の四極子モーメント，eq =は電場勾配である．η は電場勾配の非対称パラメタ∂z 2 %r=0
φxx − φyy∂2φでη =(ただし，φxi xi ≡），z 軸は電場勾配の主値の中で最大値をとる方向である．式
φzz∂x2i(12) の導出は例えば文献 [26] 等に詳しい．
以下，簡単のために η = 0 の場合，つまり軸対称性を有している系を考える．このとき，式 (12) は

HQ =

e2 qQ
(3I 2 − I 2 )．
4I(2I − 1) z

(13)

式 (13) のエネルギー固有値を求めると，
QEm= ⟨m|HQ |m⟩= A{3m − I(I + 1)}
2

'

e2 qQ
A≡
4I(2I − 1)

(

(14)

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

I= 5/2
m=!5/2
12A ≡ 2hνQ
m=!3/2
6A ≡ hνQ
m=!1/2
HQ=0

HQ≠0

図 13 核四重極相互作用によるエネルギー準位分裂（I = 5/2 の場合）．

となり，エネルギー準位間の差は
Q
Q
Q∆Em↔m−1= Em− Em−1= 3A(2m − 1)．

(15)

従って，外部磁場を印加しなくても核スピンのエネルギー状態が分裂し，そのエネルギー差に等しい電磁波を照射することによって共鳴を観測する手法が核四重極共鳴（NQR）である．例えば I = 5/2 のとき，
図 13 に示すように二重縮退した 3 つのエネルギー準位が生じる．式 (15) より，
Q∆E3/2↔1/2
= 6A

≡ hνQ
と νQ を定義すると，∆E5/2↔3/2 = 2hνQ となる．ただし，η ̸= 0 のときは，∆E3/2↔1/2 : ∆E5/2↔3/2
Q

Q

Qは整数比 1 : 2 からずれる．

ところで，A と同様 νQ も電場勾配に比例する量である．電磁気学の基礎で学ぶように，核位置を原点とし，点 (x, y, z) にある電荷が原点につくる静電ポテンシャルは φ ∼であるから，∂2φ∂z 2
3z 2 − r 2∼
r5
3 cos2 θ − 1=
r3

1（ここで，r 2 = x2 + y 2 + z 2 ）
r

φzz =

(16)

最後の式では，座標 (x, y, z) を極座標 (r, θ,ϕ ) を用いて表し直した．従って，νQ は電荷との距離 r に関
して 1/r3 の依存性（体積の次元に反比例）を持つ．

図 14(a) は，最近新たに報告された重い電子系化合物 YbCo2 Zn20 の 59 Co-NQR スペクトルである．

この物質は，γ = 7.9 J/molK2 と極めて大きな電子比熱係数を示すことで注目を集めており，磁気秩

序は 20 mK までの低温で確認されていない [11]．図 14 中に結晶構造が示されているが，この物質はCeCr2 Al20 型の立方晶（No. 227-2，空間群 F d3̄m）構造で，Yb（8a）サイトはダイヤモンド構造，Co（16d）サイトはパイロクロア構造を成す．Co 核には式 (12) の核四重極相互作用がはたらくが，Co サイ
トは 3 回回反軸を有する（位置対称 .3̄m）ので，η = 0 となるはずである．59 Co の核スピンは I = 7/2
で，式 (14) に従って 4 つのエネルギー準位に分裂し，3 つの共鳴が観測されるであろう．実験結果はま《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

!

!

(a)

Co
!

Yb
!

!

!
!

LuCo2Zn20
!

!
2.006

59

Co NQR Intensity (arb. units)!

YbCo2Zn20(b)Frequency (MHz)

2.004

!

!

νQ (MHz)1.6 1.8 2.0 2.2 2.4 3.6 3.8 4.0 4.2 4.4 5.6 5.8 6.0 6.2 6.4

2.002

2.000
0.0

0.2

0.4

0.6

Temperature (K)図 14 (a) YbCo2 Zn20 と非磁性参照物質 LuCo2 Zn20 の 59 Co-NQR スペクトル（T = 4.2 K）．右上図は YbCo2 Zn20 の結晶構造（簡単のため，Zn 元素は省略）．(b) YbCo2 Zn20 における 59 Co-NQR
周波数 νQ の温度依存性（T < 0.6 K）．文献 [10] より引用．

さにそれと合致しており，図 14(a) の上段のスペクトルから，共鳴周波数は νQ = 2.005 MHz である．
また，νQ の丁度 2 倍と 3 倍の周波数にも共鳴が観測されていることから，結晶対称性から η = 0 となる
先述の予想も正しかったことが分かる．

図 14(a) には，非磁性参照物質である LuCo2 Zn20 の測定結果も示されている（Lu3+ イオンは 4f 14 電子配置で S = L = 0）．LuCo2 Zn20 の νQ が 2.039 MHz と YbCo2 Zn20 に比べて僅かに大きいのは，式(16) で示した通り，νQ がおおよそ 1/r3 （r は電荷までの距離）の依存性を示す量であることに関係する．つまり，ランタノイド収縮のために LuCo2 Zn20 の単位胞体積が YbCo2 Zn20 に比べて小さくなり，
それによって電場勾配が増大したためだと考えられる．このように，体積変化が主因となって νQ が変化する例として，図 6 で示した YbInCu4 の価数転移が挙げられる [12]．Yb 系化合物は，Yb2+ 成分が増大すると単位胞体積が大きくなる傾向がある．YbInCu4 の低温相もそれに違わず，TV 以下で体積が増大し，同時に νQ が減少する．図 15(a) に，63 Cu-NQR スペクトルの温度依存性が示されているが，相転移温度をまたぐと，離れた周波数に別の相からの信号が不連続に出現し，ある温度幅で片方から他方に信号強度が移り変わっていく様子は，一次相転移の典型である [13]．この νQ 変化の不連続性が Pc ≈ 2.4

GPa まで残るので（図 15(b)），相図 6(a) における価数転移の境界線は常圧から Pc まで全て一次相転移であることがわかる．

YbCo2 Zn20 に戻ると，図 14(b) に示すように，約 200 mK 以下で νQ が減少する振舞いが観測されている [10]．この物質では、電気抵抗や比熱の測定からほぼ同じ温度以下で Fermi 液体状態に典型な振舞
いが報告されている [11, 14, 15]． f 電子の遍歴性が増して Yb2+ 成分が増えると，先述の様に格子の膨張が期待されるので，これが νQ 減少を引き起こすと考えられる．或いは，f が伝導バンドに寄与することによって系の電荷分布が変わり，このために電場勾配が変化することも考えられるが，どちらにしてもこの νQ の減少が重い電子状態の発達に深く関係していることに変わりはない．このような極低温下で重い電子状態の発達が観測される例は他に殆ど無く，YbCo2 Zn20 の近藤温度が非常に低いこと，またフェルミ準位における状態密度が非常に大きいというこの物質の特徴を裏付ける結果である．

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

(a)

(b)
!

(Just above TV)

LT Phase!

15.0(T=4.2K)

14.5

63

νQ (MHz)15.5 HT Phase

63

14.0

Cu-NQR0.0 0.5 1.0 1.5 2.0 2.5
P (GPa)
図 15 (a) YbInCu4 の 2.3 GPa で測定された 63 Cu-NQR スペクトルの温度依存性．高周波数側が高温相（T > TV ），低周波数側が低温相（T < TV ）からの信号．文献 [13] より引用．(b) 高温相（TV 直上）と低温相（4.2 K）で測定された 63 Cu-NQR の νQ の圧力依存性．磁気秩序が現れる臨界圧力近傍（∼ 2.5 GPa）においても両相の νQ が不連続に変化している．

3.3.2

磁場中におけるスペクトル

核四重極相互作用がはたらく状況に磁場を印加した場合は，核四重極相互作用とゼーマン相互作用の両方を考えることになる．そのときのハミルトニアンはH = HZ + Hhf + HQ

(17)となるであろう．第一項は外部磁場によるゼーマン項で，式 (1) に等しい．第二項は節 3.1 で出た磁気的な超微細相互作用で，第三項は節 3.3 で出た式 (12) である．このような状況下で NMR のスペクトルがどのような形状になるかを以下に記す．簡単のために，式 (17) の第二項 Hhf が小さいとして無視する．
このとき，式 (1) と (12) より，

H = −γn !I · Hex +hνQ
2
2{(3Iz2′ − I 2 ) + η(I+′ + I−′ )}6(18)であるが，HZ 項と HQ 項はそれぞれの主軸に基づいた表式となっているので，同一の座標軸上での表現に変換する必要がある．

今，HQ 項が HZ 項に比べて十分に小さい場合を考える．また，結晶に軸対称性があるとする（つまり η = 0 で，HQ 項として式 (13) を採用）．このとき，HZ 項を無摂動系として磁場の方向を z 軸にとり（(0, 0, Hex ) ∥ z ），これと HQ の主軸 z ′ とのなす角を θ として，二つの項は以下のようになる（図 16参照）．

HZ,0 = −γn !Iz Hex
)
hνQ 1
3
HQ =
(3 cos2 θ − 1){3Iz2 − I(I + 1)} + sin θ cos θ{Iz (I+ + I− ) + (I+ + I− )Iz }
6 2
2
*
3
2
2
+ sin2 θ(I+
+ I−
) ．
(19)
4
HQ を摂動として一次の近似でエネルギー固有値を求めると，式 (2) で与えられる無摂動系のエネルギー

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

z’
z’

z
z
θ
θ
x’

y , y’
y’ x
x’
図 16

HZ の主軸を x, y, z ，HQ の主軸を x′ , y ′ , z ′ とするとき，HQ 系に軸対称性があれば，z 軸が
x − z 面内にあるように x′ 軸を決めても一般性は失われない．このとき，Iz ′ = Iz cos θ + Ix sin θの関係を式 (13) に代入することにより，式 (19) が得られる．
′′

m=−5/2
−3/2
−1/2
hν0
1/2
3/2
5/2
HQ=0
図 17

HQ≠0

核四重極相互作用を摂動として一次の近似まで求めたエネルギー準位の模式図（I = 5/2 の場合）．

固有値 EZ,m に対して，(1)

EZ,m = ⟨m|HQ |m⟩=hνQ 3 cos2 θ − 1{3m2 − I(I + 1)}
6
2

(20)

だけシフトを引き起こす．その様子を図 17 に示す．m = ±1/2 の準位は同じだけシフトするためにエ

ネルギー差は変わらないが，それ以外の準位間はこのエネルギーシフトのために元の無摂動系とは異なるエネルギー差になり，図 18 に示すように m = 1/2 ↔ −1/2 の遷移とはずれて共鳴線が観測される．

m = 1/2 ↔ −1/2 の共鳴信号を一般にセンターラインと呼ぶのに対し，それ以外はサテライトラインと呼ばれる．式 (20) から分かるように，サテライトラインの共鳴位置は θ に依存し，1st サテライト間の分裂幅は θ = 0 のときは 2νQ に，θ = 90◦ のときは νQ に相当する（図 18(a)）．式 (20) において，3 cos2 θ − 1 = 0 になる角度（θ ≈ 54.7◦ ）では，サテライトラインが全てセンターラインの位置に集まり，この角度はマジックアングルと呼ばれる．ここまでの話は，単結晶試料（或いは配向した粉末試料）
を，磁場の方向に対して θ を制御して測定した場合に相当する．

次に，上述の内容を多結晶試料，または粉末試料に適用する．これまでと同様に，核四重極相互作用を摂動項として取り扱い，一次までの近似で，ナイトシフトの影響は無視する．節 3.2 でシフトの異方性に

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

z // H0(a)-10

Intensity (arb. unit)

(b)

!

ν0

2 νQ
Intensity (arb. unit)

0
10

30
40

ν0

!

50
Intensity (arb. unit)

Angle θ (Degrees)

20

z'
θ

60
70
80
90
νQ

ν0

100

Frequency

Intensity (arb. unit)

ν0

θ =π

ν0

Frequency

図 18(a) 共鳴周波数の角度 θ 依存性（I = 5/2 の場合．一次摂動で計算）．θ = 0◦ のときと
θ = 90 のときの 1st サテライトの分裂幅がそれぞれ 2νQ と νQ に相当する．本文参照．(b) θ を変化させたときのスペクトルの形状と信号強度（I = 5/2 の場合．一次摂動で計算）．右の単位円は，磁場の方向（z 軸）と四重極相互作用の主軸（z ′ 軸）とのなす角 θ を示す．
◦

依る粉末パターンを考えたのと同様に，今の場試料粉末中の微結晶が等しい θ を成していれば，同一の周波数（或いは磁場）にて信号が観測されるはずである．図 18(b) に示すように単位円を考えると，同一の θ を成す微結晶を見出す確率は sin θ に比例する．求める粉末パターンは全ての角度のスペクトルを足し合わせたものであり，図 19(a) におけるライン A（理想的にはライン B）のようになる．図 19(a) を見てもわかるように，粉末パターンの 1st サテライトは θ = 90◦ に一致するところでピークを持つので（図中縦の点線），1st サテライト間の分裂は νQ に相当する．ゼロ磁場にて NQR 測定を行う際に共鳴周波数が未知の場合には，まず NMR スペクトルを測定して，このように νQ の値を見積もる場合が多い．

粉末パターンにおけるセンターラインは，二次の摂動まで考慮に入れると図 19(b) に示すように二つの山を持つような形状になることが分かる．このような形状は，νQ が ν0 (= γN Hex /2π) に対して比較的大きな場合や，スペクトル幅が狭い場合などに観測される．サテライトラインが現れる周波数（磁場）も一次の近似に対してシフトするが，近似が用いられる範囲内では，1st サテライト間の分裂幅が νQ に相当するということに変わりはない．

図 20 に，図 14 で既出の YbCo2 Zn20 について得られた 59 Co-NMR の粉末パターンを示す．上述の様に，センターラインは二つの山を持つ形状を成しているが，図 19(b) とは，横軸が磁場 ⇔ 周波数の違い

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）!

!

1st satellite lines

60°
70°
80°
90°

B
A

Center line

!

θ=

10°
20°
30°
40°
50°

Intensity (arb. units)

!$#
θ=

!

Intensity (arb. units)

!"#

ν0

ν0

Frequency

Frequency

図 19 (a) 一次摂動で計算した θ = 10◦ ∼ 90◦ のスペクトルとそれらの和を取ったものがライン A.
ライン B は理想的なスペクトル．(b) 二次摂動で計算した粉末パターン．(a) と (b) 共に横軸は周波数で，シフトはゼロとしている．

!

59

!

Intensity (arb. unit)

YbCo2Zn20 Co-NMR

2.6

2.7

2.8

2.9

3.0

3.1

3.2

3.3

Mangetic field (T)図 20 YbCo2 Zn20 の 59 Co-NMR 測定と計算による粉末パターンとスペクトル．計算では，νQ =2.004 MHz，K11 = 1.1%，K33 = −0.4% を用いた．横軸は磁場で，図 19 とは異なることに注意．
測定周波数は ν0 = 30.09 MHz．約 2.71 T に見られる信号は，NMR 用コイルに用いた銅線（63 Cu）の信号．

があるために，山の高さが左右反対になっている．1st サテライト間の分裂幅が約 0.2 T であり，59 Co
は γN /2π = 10.03 MHz/T であるから，νQ ∼ 10.03 × 0.2 ∼ 2 MHz となって，図 14 の 59 Co-NQR 測定で得られた値によく一致していることが分かるであろう．図 20 には，計算によって求めたスペクトルも示してあるが，近年はパソコンの発達のおかげで，式 (18) を摂動や対角化によって解くことは難しくないので，NMR スペクトルも比較的良く再現することができる．

3.4 磁気秩序状態
最後に，磁気秩序状態における NMR スペクトルについて紹介する．重い電子系化合物は局在性の強い

f 電子の磁気モーメントを持ち，重い電子状態は磁気秩序相の近傍に出現する現象であるから，磁気秩序

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

(b)
(a)

Si
Hex
U
Hex+∆H ord
inHex−∆H ordin(c)
図 21 URu2 Si2 について，8.3 kbar の圧力下で測定された 29 Si-NMR スペクトルの温度変化．測定に用いた試料は，粉末試料を大きな磁気異方性（χc > χab ）を利用して磁場中で配向させ，スタイ
キャストで固めたもの．測定磁場は Hex (= 4.3 T) ∥ c．内挿図は中性子回折実験によって決められた反強磁性磁気秩序構造 [16, 17]．(b) U の磁気モーメントが Si サイトに生じる内部磁場の模式図．
ord(c) 内部磁場 Hin（図中では Hin と示されている）と反強磁性体積分率の温度依存性．文献 [19] より引用．

状態における NMR スペクトルから得られる情報について知っておくことは役に立つであろう．良く知られるように，NMR や NQR 測定は磁気秩序に敏感な実験手段である．以下に，URu2 Si2 の 29 Si-NMR
測定について述べる．

URu2 Si2 （結晶構造は図 8）は，Tc = 1.4 K の重い電子系超伝導体として良く知られる．加えて，
TO = 17.5 K にて秩序機構の明らかになっていない，いわゆる Hidden Order（隠れた秩序）と呼ばれる相転移（以降，HO と記す）を示すことでも有名である [18]．その研究の歴史をみると，HO 相内では0.03 µB 程度の非常に微小な反強磁性秩序モーメントが観測されていたが，この微小モーメントの秩序がHO の機構であるとすると，TO での比熱の大きな飛びが説明できない．従って，この物質の発見からほぼ 15 年間，その微小モーメントの起源について議論されてきたのであるが，以下に記すように NMR 測

定から重要な知見が得られた．

図 21(a) は，8.3 kbar の圧力下における 29 Si-NMR スペクトルの温度依存性である [19]．TO 以上の常磁性状態では 1 本の信号が観測されているのに対して，TO 以下では常磁性状態の信号が最低温度まで残存し，加えてその両脇に異なる信号が現れている．両脇に 2 本の信号が現れる機構は，次の様に解釈される．今，磁気モーメントが内挿図に示された type-I と呼ばれる反強磁性構造（図 21(a) 内挿図）で秩序しているところに，外部磁場 Hex が c 軸に平行に印加されると，図 21(b) に示したように，Si サイトにord
は Hex と U の磁気秩序モーメントが物質内につくる磁場 ∆Hin
が共に c 軸に平行にかかることになる．

ord
ord
従って，type-I の秩序構造に従って，Hex + ∆Hinと Hex − ∆Hinの磁場を受ける Si サイトが交互に並ぶので，反強磁性領域からの信号は 2 等分に分裂することになる．図 21(a) で両脇の反強磁性信号間の分

ord
裂幅は約 1.6 MHz で，29 Si は γN /2π = 8.45577 MHz/T であるから，∆Hin∼ 1.6/(γN /2π) ∼ 0.09 Tである．一方，内部磁場を受けていない信号も中心に残っているので，これが HO 相からの信号に対応

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

(a)

θ =0º

Hex
(b)

∆Hinord!

180°∑ 2π! (θ )

θ =90ºθ = 0°!"(θ )!

θ

Intensity

ideal lineν 0−∆νν0ν0+∆νFrequency
図 22 多結晶試料，或いは粉末試料が磁気秩序状態にあるときのスペクトル形状．(a) 多結晶試料，
或いは粉末試料中の微結晶が磁場に対してランダムな方向を向いている様子．θ は外部磁場 ex とord
磁気秩序による ∆ inの成す角度．(b) 各微結晶からの信号の共鳴周波数が θ だけに依るものと◦◦し，θ = 0 から 180 までを 15◦ 間隔で加え合わせて得られたスペクトルと理想的な粉末パターン
ord（破線）．∆ν = (γN /2π) ∆Hin．

する．ここで，HO 相の高圧側には type-I 構造の反強磁性相が存在することは分かっているので，得られた結果は，8.3 kbar の圧力下で既に試料の一部が反強磁性秩序を起こしていることを意味する．NMRの信号強度は原理的に観測している原子数に比例するので，図 21(a) の常磁性相と反強磁性相の信号強度比を取れば，両者の体積分率が見積もられる．その結果は図 21(c) に示されており，2 kbar から 10 kbarに圧力を増加させていくにつれて，反強磁性相の領域が徐々に増えていくことが分かる．ここで重要なのordは，この圧力領域で ∆Hinの値が殆ど変わらないことである．こうして，先述の微小磁気モーメントの根拠となった中性子回折の非常に弱いブラッグピークの原因は，磁気モーメントの小ささにあるのではな
く，反強磁性領域が常圧下では僅かしかないことに起因していたことが明らかになったのである．

上述のように，磁気秩序状態では物質内部に自発的な内部磁場が生じ，NMR の共鳴は外部磁場とこordの内部磁場の和 Hex + ∆Hinに対して起こる．測定している試料が多結晶試料や粉末試料の場合，これまでと同様に試料は微結晶の集まったものであると考えると，図 22(a) に示すように，微結晶は磁場中でランダムな方向を向き，それぞれに固定された磁気秩序モーメントによる内部磁場が生じている．

ord
ord
従って、核位置における磁場の和の大きさは Hex − ∆Hinから Hex + ∆Hinまで分布するであろう．

ord
各々の微結晶に対する共鳴周波数（或いは共鳴磁場）が Hex と ∆Hinの成す角 θ だけに依るとみなすことができるときは，スペクトルの形状を計算するのは簡単で，シフトの異方性に起因する粉末パターンを考えた時と同様，各 θ に対応するスペクトル g(θ) を θ = 0◦ から 180◦ まで加え合わせる．もし，ord
Hex ≫ ∆Hinの場合は，図 22(b) に示すように，台形型の粉末パターンが得られる．ただし，このスペクトルの形状のみからは強磁性秩序と反強磁性秩序の区別はつかない．

図 23 は，CePd2 Si2 粉末試料を 4.2 K で測定して得た 29 Si-NMR スペクトルである [20]．CePd2 Si2の結晶構造も図 8 で示した ThCr2 Si2 型である．この物質は TN ∼ 10 K の反強磁性体 [21] であり，約27 kbar の圧力を印加すると，反強磁性相が消失して超伝導相（Tc ∼ 0.4 K）が出現することが知られている [22]．常圧下で測定された図 23 のスペクトルは，上で説明した典型的な磁気秩序状態における粉末パターンである．ところで，CePd2 Si2 における磁気構造は，次の 105 Pd-NQR 測定によって，ある程度推測することができる．105 Pd の核スピンは I = 5/2 であるので，常磁性状態では図 13 のエネルギー準位分裂に従って，2 本の 105 Pd-NQR 信号が観測される．これに相当するのが，図 24 の上段のスペクト《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

図 23

CePd2 Si2 の 反 強 磁 性 秩 序 相 に お け る 29 Si-NMR 粉 末 パ タ ー ン ．測 定 周 波 数 は ν0 =11.1 MHz．文献 [20] より引用．左図は，中性子回折実験によって決められた反強磁性磁気秩序構造 [21]．

NQR intensity (arb. units)

CePd2Si2105

Pd-NQRT = 11 K

T = 4.2 K

4

5

8

9

Frequency (MHz)図 24 CePd2 Si2 の常磁性状態（T = 11 K> TN ）と反強磁性状態（T = 4.2 K< TN ）における
105
Pd-NQR スペクトル [23]．

ルである．一方，磁気秩序状態では自発的な内部磁場が生じるために NQR スペクトルが変化すると予想されるが，図 24 の上段と下段のスペクトルを比較する限り変化は生じていない．このことは，周囲の磁気モーメントが生じる内部磁場が，対称性の良い Pd（4d）サイトでは（つまり形状因子によって）丁度相殺していることを示唆している．このことから，磁気秩序構造が反強磁性的であるということがわかる．
これまで述べてきた NMR や NQR 測定の例は，全て磁性元素そのものではなく，磁性サイトの周囲を占める非磁性原子を観測たものであった（図 3(b) に示す状況に相当）．以下では，磁性原子を直接観測した例について述べる．節 3.1 で記した通り，磁性を示すランタノイドやアクチノイド元素の直接観測は容易ではないが，ある幾つかの条件下ではそれが可能になる．その一つが，磁気秩序状態での NMR 測定である．これは，磁気モーメントが秩序するために磁気的な揺らぎ（核スピンの緩和に大きく寄与する）が常磁性状態に比べて大幅に抑えられるため，緩和時間が比較的長くなって信号が観測しやすくなるためである．ただし，オンサイトの磁気モーメントが秩序することによって核位置には大きな内部磁場が生じるので，信号検出を成功させるためにはその内部磁場の大きさを探し当てなくてはならない．

図 25 に，Eu 化合物である EuGa4 について Eu-NMR 測定に成功した例を示す [24]．Eu イオンは，+3 価のときは S = L = 3 で J = 0 となり非磁性であるが，+2 価の状態では J = 7/2（S = 7/2，
L = 0）で gJ JµB = 7 µB の比較的大きな磁気モーメントを有する．この物質は，図 25(a) に示す通り，
図 8 と同じ空間群 I4/mmm の正方晶系の構造を持ち，TN = 16 K の反強磁性体である．図 25(b) には，

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

(a)

(b)

(c)

図 25 EuGa4 の (a) 結晶構造（BaAl4 型正方晶構造），(b) 中性子回折実験で決められた反強磁性
磁気構造 (c) 153 Eu のゼロ磁場 NMR スペクトル（上段）．下段は Eu 核位置での内部磁場が核四重極相互作用の主軸（z 軸 ∥ 結晶 c 軸）に垂直で，内挿図に示した νQ の分布を仮定して求めた計算結果．文献 [24] より引用．

中性子回折実験によって決められた磁気構造を示す [25]．この反強磁性状態では，Eu 核位置に大きな内部磁場（後に記すように約 27 T）が生じるため，核スピンとの間にゼーマン相互作用が生じる．加えて，153 Eu の核スピンは I = 5/2 で，核四重極相互作用もはたらくため，状況としては式 (18) において，
ord
Hex を内部磁場Hinに置き換えたものに等しい．実際に観測された 153 Eu-NMR スペクトルが図 25(c)の上段に示されている．先にも記した通り，この NMR 信号観測には外部磁場を用いず，物質内の自発的な内部磁場を利用したものであるから，このような共鳴はゼロ磁場 NMR と呼ばれる．図 25(c) の下段には，実験結果を再現するシミュレーションが示されているが，ここでは結晶には c 軸方向に 4 回対称ord性があることから，HQ 項に η = 0 の式 (13) を用い，νQ = 30.5 MHz，Eu 核位置で Hint
= 27.08 T，
ord
Hin ⊥ c としている．この結果は，図 25(b) で示した磁気構造に合致するものであるが，このように磁気秩序状態でのゼロ磁場 NMR 測定，特に四重極相互作用も共存する状況下では，磁気構造に関する情報が得られる場合が多い．

4 終わりにNMR や NQR 測定から得られる静的な情報を中心に，f 電子系で報告された例を紹介しながら解説を行った．講義の前にこれだけ読めば NMR については大体分かる，というものにしたいと思い執筆したが，NMR スペクトルの形状やそこから得られる情報については大方お分かりいただけるものと思う．
一方，動的な情報について記すスペースが無かったために片手落ちになってしまったのが残念である．

NMR 測定から得られる情報量の多さや，様々な分野への適応能力の高さには，学生時代からこの測定を行っている身に取っても，いつも新鮮で驚きさえ感じるときもある．本稿が NMR や重い電子系の理解の一助となれば幸いである．

5 謝辞
世話人の高野俊氏，推敲を手伝っていただいた小山岳秀氏に感謝致します．本テキスト執筆に当たって，[26–28] の書を参考にしました．

《講義ノート》

物性研究・電子版 Vol. 4, No. 1, 041206 （2015年2月号）

参考文献[1] S. Suga et al., J. Phys. Soc. Jpn. 78, 074704 (2009).
[2] T. Mito et al., Phys. Rev. B 75, 224409 (2003).
[3] T. Mito et al., Phys. Rev. B 75, 134401 (2007).
[4] T. Mito et al., J. Phys. Soc. Jpn. 81, 033706 (2012).
[5] J. L. Sarrao et al., Phys. Rev. B 59, 6855 (1999).
[6] N. Tsujii et al., Phys. Rev. B 56, 8103 (1997).
[7] C. Rossel et al., Phys. Rev. B 35, 1914 (1987).
[8] F. Steglich, et al., Phys. Rev. Lett. 43, (1979) 1892.
[9] K. Ueda et al., J. Phys. Soc. Jpn. 56, 867 (1987).
[10] T. Mito et al., J. Phys. Soc. Jpn. 82, 103704 (2013).
[11] M. S. Torikachvili et al., Proc. Natl. Acad. Sci. U.S.A. 104, 9960 (2007).
[12] H. Nakamura et al., J. Phys. Soc. Jpn. 59, 28 (1990).
[13] T. Koyama et al., Phys. Rev. B 71, 184437 (2005).
[14] Y. Saiga: Dr. Thesis, Saitama University, (2008).
[15] M. Ohya et al., J. Phys. Soc. Jpn. 79, 083601 (2010).
[16] C. Broholm et al., Phys. Rev. Lett. 58, 1467 (1987).
[17] C. Broholm et al., Phys. Rev. B 43, 12809 (1991).
[18] レビューとして，J. A. Maydosh et al., Rev. Mod. Phys. 83, 1301 (2011).
H. Amitsuka et al., J. Magn. Magn. Mater. 310, 214 (2007).
[19] K. Matsuda, et al., Phys. Rev. Lett. 87, 087203 (2001).
[20] Y. Kawasaki et al., Phys. Rev. B 58, 8634 (1998).
[21] B. H. Grier et al., Phys. Rev. B 29, 2664 (1984).
[22] F. M. Grosche et al., Physica B 237-238, 197 (1997).
[23] S. Kawasaki, unpublished.
[24] M. Yogi et al., J. Phys. Soc. Jpn. 82, 103701 (2013).
[25] T. Kawasaki et al., JPS Conf. Proc. 1, 014009 (2014).
[26] A. Abragam 著，富田和久・田中基之 共訳，「核の磁性」，（吉岡書店，1964）.
[27] 朝山邦輔 著，「遍歴電子系の核磁気共鳴」，（裳華房，2001）.
[28] C. P. Slichter 著，益田義賀・雑賀亜幌 共訳，「磁気共鳴の原理」，（岩波書店，1966）.