ハミルトニアンの選び方で固体を考える 第5部 第2章 §5.2.4 ★★★ 2つのサイトの模型を解くと、直接交換 − 2 K -2K − 2 K (強磁性)と運動交換 4 t 2 / ( U − V ) 4t^2/(U-V) 4 t 2 / ( U − V ) (反強磁性)の和が J J J になる。酸素を挟むと4次の超交換になり、結合の角度で符号が変わる(Goodenough–Kanamori 則)。スピン軌道相互作用を戻すと、反対称な DM 相互作用 D ⋅ ( S ^ i × S ^ j ) \bm{D}\cdot(\hat{\bm{S}}_i\times\hat{\bm{S}}_j) D ⋅ ( S ^ i × S ^ j ) が加わる
2つのサイト:直接交換と運動交換
式 (5.2.4) の模型に電子を2個入れる。例題 5.2.1 の分け方で、3重項は空間の波動関数が共有結合的な反対称のもの1つだけで、飛び移り t t t では2個が同じサイトに入る状態へ移れない(Pauli の排他律)。したがって3重項のエネルギーは t t t によらず、式 (5.2.7) から
E T = V − K E_T=V-K E T = V − K
である。1重項は3つあり、共有結合的な状態(エネルギー V + K V+K V + K )と、イオン的な状態の対称な組 1 2 ( ∣ 1 ↑ ↓ ⟩ + ∣ 2 ↑ ↓ ⟩ ) \frac{1}{\sqrt2}(| 1\!\uparrow\downarrow \rangle +| 2\!\uparrow\downarrow \rangle ) 2 1 ( ∣1 ↑↓ ⟩ + ∣2 ↑↓ ⟩) (U + K U+K U + K 。対で移る項が + K +K + K を足す)と反対称な組(U − K U-K U − K )である。飛び移りは共有結合的な状態と対称なイオン的な状態だけを行列要素 − 2 t -2t − 2 t で結ぶので、1重項の最低のエネルギーは 2 × 2 2\times2 2 × 2 行列( V + K − 2 t − 2 t U + K ) \begin{pmatrix}V+K&-2t\\-2t&U+K\end{pmatrix} ( V + K − 2 t − 2 t U + K ) の小さいほうの固有値で、
J = E T − E S = − 2 K − U − V 2 + ( U − V 2 ) 2 + 4 t 2 ≈ − 2 K + 4 t 2 U − V ( t ≪ U − V ) (5.2.14) J=E_T-E_S=-2K-\frac{U-V}{2}+\sqrt{\Bigl(\frac{U-V}{2}\Bigr)^2+4t^2}\;\approx\;-2K+\frac{4t^2}{U-V}\qquad(t\ll U-V)\tag{5.2.14} J = E T − E S = − 2 K − 2 U − V + ( 2 U − V ) 2 + 4 t 2 ≈ − 2 K + U − V 4 t 2 ( t ≪ U − V ) ( 5.2.14 )
を得る。これが2サイトの模型の厳密な J J J で、右の近似は式 (5.2.8) の2次の摂動と一致する(Fock 空間の 6 × 6 6\times6 6 × 6 の行列の対角化と合わせた。検算 scripts/verify_ch5-02.py)。
式 (5.2.14) の2つの項は符号が逆である。第1項 − 2 K -2K − 2 K は直接交換 で、軌道の重なりの電荷の Coulomb エネルギー(命題 5.2.2 )のために3重項を下げ、強磁性の向きに働く。第2項 4 t 2 / ( U − V ) 4t^2/(U-V) 4 t 2 / ( U − V ) は運動交換 で、1重項だけが隣のサイトへ一度飛び移って戻る仮想的な過程を持てる(3重項は Pauli の排他律で飛べない)ために1重項を下げ、反強磁性の向きに働く。分母の U − V U-V U − V は、電子を隣へ移してサイトに2個詰めるのにかかるエネルギーである。V = K = 0 V=K=0 V = K = 0 とすると第1部第1章 の J = 4 t 2 / U J=4t^2/U J = 4 t 2 / U に戻る。
軌道が1つずつで、電子がサイトに1個ずつの場合、遷移金属の化合物では普通 4 t 2 / U 4t^2/U 4 t 2 / U が 2 K 2K 2 K より大きく、反強磁性になる。Hubbard 模型(V = K = 0 V=K=0 V = K = 0 )については、連結な2部格子(隣どうしが別の副格子に属する格子)で半分詰めのとき、U > 0 U>0 U > 0 ならどんな大きさでも基底状態の全スピンが ∣ N A − N B ∣ / 2 \lvert N_A-N_B\rvert /2 ∣ N A − N B ∣ /2 (N A N_A N A 、N B N_B N B は副格子のサイト数)になることが証明されている(Lieb(1989)、Physical Review Letters 62、1201)。運動交換の反強磁性は摂動の極限だけの結果でなく、この意味で U U U の大きさによらない。強磁性を作るには、軌道が直交したまま重なる(t = 0 t=0 t = 0 で K > 0 K>0 K > 0 )ような配置や、同じサイトに複数の軌道があって Hund 則が働く場合が必要になる。次の超交換の 90° の結合がその例である。水素分子では、例題 5.2.1 のイオン的な状態の重みが t / U t/U t / U で決まり、第2部第1章 で分子軌道の近似が「半分イオン的」になって解離を誤ったのは、U U U を無視して重みを 1 2 \frac12 2 1 に固定したためである。
例題 5.2.4 直接交換と運動交換のどちらが勝つか ★★★
オリジナル 式 (5.2.4) で U = 4 e V U=4\,\mathrm{eV} U = 4 eV 、V = 1 e V V=1\,\mathrm{eV} V = 1 eV 、K = 0.02 e V K=0.02\,\mathrm{eV} K = 0.02 eV とする。
(1) t = 0.3 e V t=0.3\,\mathrm{eV} t = 0.3 eV のとき、J J J を2次の摂動と式 (5.2.14) の厳密な形で求め、比べよ。
(2) 2次の摂動で J = 0 J=0 J = 0 になる t c t_{\mathrm c} t c を求めよ。t < t c t<t_{\mathrm c} t < t c でどちらの磁性になるか。厳密な形での t c t_{\mathrm c} t c と比べよ。
解答 (1) 2次の摂動は J = − 2 K + 4 t 2 / ( U − V ) = − 0.04 + 0.36 / 3 = 0.080 e V J=-2K+4t^2/(U-V)=-0.04+0.36/3=0.080\,\mathrm{eV} J = − 2 K + 4 t 2 / ( U − V ) = − 0.04 + 0.36/3 = 0.080 eV 。厳密な形は − 0.04 − 1.5 + 2.25 + 0.36 = 0.0756 e V -0.04-1.5+\sqrt{2.25+0.36}=0.0756\,\mathrm{eV} − 0.04 − 1.5 + 2.25 + 0.36 = 0.0756 eV で、摂動より 5.6 % 小さい。展開の次の項は − 16 t 4 / ( U − V ) 3 -16t^4/(U-V)^3 − 16 t 4 / ( U − V ) 3 の大きさで、t / ( U − V ) = 0.1 t/(U-V)=0.1 t / ( U − V ) = 0.1 ではその程度のずれになる。
(2) 4 t c 2 / ( U − V ) = 2 K 4t_{\mathrm c}^2/(U-V)=2K 4 t c 2 / ( U − V ) = 2 K から t c = K ( U − V ) / 2 = 0.03 = 0.173 e V t_{\mathrm c}=\sqrt{K(U-V)/2}=\sqrt{0.03}=0.173\,\mathrm{eV} t c = K ( U − V ) /2 = 0.03 = 0.173 eV 。t < t c t<t_{\mathrm c} t < t c では直接交換が勝って J < 0 J<0 J < 0 、強磁性になる。厳密な形では ( ( U − V ) / 2 ) 2 + 4 t 2 = ( U − V ) / 2 + 2 K \sqrt{((U-V)/2)^2+4t^2}=(U-V)/2+2K (( U − V ) /2 ) 2 + 4 t 2 = ( U − V ) /2 + 2 K から t c 2 = K ( U − V ) / 2 + K 2 t_{\mathrm c}^2=K(U-V)/2+K^2 t c 2 = K ( U − V ) /2 + K 2 、t c = 0.174 e V t_{\mathrm c}=0.174\,\mathrm{eV} t c = 0.174 eV で、差は1 % 以内である。飛び移り t t t と重なりの電荷による K K K は、どちらも2つの軌道の重なりで決まり、原子を離すとどちらも小さくなる。そのため t t t と K K K を独立に動かせるわけではなく、どちらが勝つかは軌道の形と距離による。(検算 scripts/verify_ch5-02.py)
酸素を挟んだ超交換
遷移金属の酸化物では、磁性を担う金属イオンどうしは直接に接しておらず、間に酸素のイオン O2 − ^{2-} 2 − がある。金属の d d d 軌道どうしの重なりは小さく、t t t も K K K もほとんど 0 だが、それでも金属のスピンの間に大きな反強磁性の相互作用が働く。これを超交換 と呼ぶ。Anderson は、酸素の p p p 軌道を仲立ちにした飛び移りの高次の過程として超交換を定式化した(Anderson(1959)、Physical Review 115、2)。
銅の酸化物を例に、Cu2 + ^{2+} 2 + (3 d 9 3d^9 3 d 9 )を「d d d 殻に正孔が1個」、O2 − ^{2-} 2 − (2 p 6 2p^6 2 p 6 )を「正孔が無い」とみなし、正孔の言葉で模型を作る。Cu–O–Cu が一直線に並び、d d d 軌道1つと酸素の p p p 軌道1つだけを残すと、
H ^ d p d = Δ c t n ^ p − t p d ∑ j = 1 , 2 ∑ σ ( d ^ j σ † p ^ σ + p ^ σ † d ^ j σ ) + U d ∑ j = 1 , 2 n ^ j ↑ d n ^ j ↓ d + U p n ^ ↑ p n ^ ↓ p (5.2.15) \hat H_{\mathrm{dpd}}=\Delta_{\mathrm{ct}}\,\hat n_{\mathrm p}-t_{\mathrm{pd}}\sum_{j=1,2}\sum_\sigma\bigl(\hat d^\dagger_{j\sigma}\hat p_\sigma+\hat p^\dagger_\sigma\hat d_{j\sigma}\bigr)
+U_{\mathrm d}\sum_{j=1,2}\hat n^{\mathrm d}_{j\uparrow}\hat n^{\mathrm d}_{j\downarrow}+U_{\mathrm p}\hat n^{\mathrm p}_\uparrow\hat n^{\mathrm p}_\downarrow\tag{5.2.15} H ^ dpd = Δ ct n ^ p − t pd j = 1 , 2 ∑ σ ∑ ( d ^ j σ † p ^ σ + p ^ σ † d ^ j σ ) + U d j = 1 , 2 ∑ n ^ j ↑ d n ^ j ↓ d + U p n ^ ↑ p n ^ ↓ p ( 5.2.15 )
である。d ^ j σ \hat d_{j\sigma} d ^ j σ は銅 j j j の、p ^ σ \hat p_\sigma p ^ σ は酸素の正孔の消滅演算子、Δ c t \Delta_{\mathrm{ct}} Δ ct は正孔を銅から酸素へ移すエネルギー(電荷移動エネルギー)、U d U_{\mathrm d} U d ・U p U_{\mathrm p} U p は同じ軌道に正孔が2個入ったときの斥力である。正孔は2個で、低い準位の組は各銅に1個ずつの4状態である。
2つの銅は直接つながっていないので、2次の摂動では正孔が酸素へ出て同じ銅へ戻る過程しかなく、スピンによらない。スピンに依る最初の項は4次で、2種類の過程がある。1つは、正孔が銅1から酸素を経て銅2へ移り(銅2に正孔2個、エネルギー U d U_{\mathrm d} U d )、同じ道を戻る過程で、運動交換の t t t を t e f f = t p d 2 / Δ c t t_{\mathrm{eff}}=t_{\mathrm{pd}}^2/\Delta_{\mathrm{ct}} t eff = t pd 2 / Δ ct に置き換えたものにあたる。もう1つは、両方の銅の正孔が酸素へ出て(酸素に正孔2個、エネルギー 2 Δ c t + U p 2\Delta_{\mathrm{ct}}+U_{\mathrm p} 2 Δ ct + U p )から戻る過程である。どちらも、2個の正孔が同じ軌道に入る途中の状態を通るので、Pauli の排他律のために、1重項だけがこの途中の状態を経由できる。経路の数と分母を数えると、
J 180 ∘ = 4 t p d 4 Δ c t 2 ( 1 U d + 2 2 Δ c t + U p ) > 0 (5.2.16) J_{180^\circ}=\frac{4t_{\mathrm{pd}}^4}{\Delta_{\mathrm{ct}}^2}\Bigl(\frac{1}{U_{\mathrm d}}+\frac{2}{2\Delta_{\mathrm{ct}}+U_{\mathrm p}}\Bigr)>0\tag{5.2.16} J 18 0 ∘ = Δ ct 2 4 t pd 4 ( U d 1 + 2 Δ ct + U p 2 ) > 0 ( 5.2.16 )
となり、反強磁性である(Δ c t ≫ U d \Delta_{\mathrm{ct}}\gg U_{\mathrm d} Δ ct ≫ U d では第1項 4 t e f f 2 / U d 4t_{\mathrm{eff}}^2/U_{\mathrm d} 4 t eff 2 / U d が残る)。図 5.2.2 は、3つのサイトの模型を厳密に対角化した J J J と式 (5.2.16) を比べたものである。t p d / Δ c t = 0.1 t_{\mathrm{pd}}/\Delta_{\mathrm{ct}}=0.1 t pd / Δ ct = 0.1 では差は 5 % だが、0.3 では厳密な値が4次の式の 0.70 倍で、t p d / Δ c t t_{\mathrm{pd}}/\Delta_{\mathrm{ct}} t pd / Δ ct が 0.3 を超える範囲では4次の式は大きめの見積もりになる。
図 5.2.2 一直線の Cu–O–Cu の超交換。Δ c t = 3 e V \Delta_{\mathrm{ct}}=3\,\mathrm{eV} Δ ct = 3 eV 、U d = 8 e V U_{\mathrm d}=8\,\mathrm{eV} U d = 8 eV 、U p = 4 e V U_{\mathrm p}=4\,\mathrm{eV} U p = 4 eV として、正孔2個の3サイトの模型 (5.2.15) を厳密に対角化した J = E T − E S J=E_T-E_S J = E T − E S (点線)と、4次の摂動の式(実線)。厳密な値と式の比は t p d / Δ c t = 0.1 , 0.2 , 0.3 , 0.4 t_{\mathrm{pd}}/\Delta_{\mathrm{ct}}=0.1,\,0.2,\,0.3,\,0.4 t pd / Δ ct = 0.1 , 0.2 , 0.3 , 0.4 で 0.955 , 0.841 , 0.700 , 0.563 0.955,\,0.841,\,0.700,\,0.563 0.955 , 0.841 , 0.700 , 0.563 である。
結合の角度と Goodenough–Kanamori 則
金属–酸素–金属の角度が 90° のときは、2つの金属の d d d 軌道がつながる相手が、酸素の互いに直交する p x p_x p x と p y p_y p y になる(図 5.2.3 )。4次の過程で両方の正孔が酸素へ出ても、入る軌道が別々なので、Pauli の排他律による3重項の制限は無い。1重項と3重項の差は、酸素の上で2個の正孔が p x p_x p x ・p y p_y p y に1個ずつ入った途中の状態のエネルギーだけに残り、その差は酸素の中の交換積分 K p K_{\mathrm p} K p (p x p_x p x と p y p_y p y の間。5.2.3 節の Hund の第1の規則と同じもの)による。U ′ = U p − 2 K p U'=U_{\mathrm p}-2K_{\mathrm p} U ′ = U p − 2 K p を p x p_x p x と p y p_y p y の正孔の Coulomb 積分とすると、途中の状態は3重項で 2 Δ c t + U ′ − K p 2\Delta_{\mathrm{ct}}+U'-K_{\mathrm p} 2 Δ ct + U ′ − K p 、1重項で 2 Δ c t + U ′ + K p 2\Delta_{\mathrm{ct}}+U'+K_{\mathrm p} 2 Δ ct + U ′ + K p であり、
J 90 ∘ = − 4 t p d 4 Δ c t 2 ⋅ 2 K p ( 2 Δ c t + U ′ ) 2 − K p 2 < 0 (5.2.17) J_{90^\circ}=-\frac{4t_{\mathrm{pd}}^4}{\Delta_{\mathrm{ct}}^2}\cdot\frac{2K_{\mathrm p}}{(2\Delta_{\mathrm{ct}}+U')^2-K_{\mathrm p}^2}<0\tag{5.2.17} J 9 0 ∘ = − Δ ct 2 4 t pd 4 ⋅ ( 2 Δ ct + U ′ ) 2 − K p 2 2 K p < 0 ( 5.2.17 )
となる。90° の結合の J J J は負で、弱い強磁性になる。K p = 0 K_{\mathrm p}=0 K p = 0 なら J 90 ∘ = 0 J_{90^\circ}=0 J 9 0 ∘ = 0 である(どちらも4つの軌道の模型の厳密な対角化で確かめた)。
図 5.2.3 金属–酸素–金属の結合の角度と超交換の符号。180° では2つの金属の正孔が酸素の同じ p p p 軌道を経由するので、Pauli の排他律で1重項だけが仮想的な過程を持つ。90° では経由する軌道が直交する p x p_x p x と p y p_y p y に分かれ、酸素の中の Hund の交換 K p K_{\mathrm p} K p の分だけ3重項が下がる。
このように、超交換の符号と大きさは、金属の d d d 軌道の占め方と結合の角度から決まる。180° の結合で半分詰まった軌道どうしなら強い反強磁性、90° なら弱い強磁性、というこの経験則をGoodenough–Kanamori 則 と呼ぶ(Goodenough(1955)、Physical Review 100、564。Kanamori(1959)、Journal of Physics and Chemistry of Solids 10、87)。式 (5.2.15) の模型の2つのエネルギー U d U_{\mathrm d} U d と Δ c t \Delta_{\mathrm{ct}} Δ ct のどちらが小さいかで、電荷を動かすのに要る最小のエネルギー(絶縁体のギャップ)の中身が変わる。U d < Δ c t U_{\mathrm d}<\Delta_{\mathrm{ct}} U d < Δ ct の物質を Mott–Hubbard 型、Δ c t < U d \Delta_{\mathrm{ct}}<U_{\mathrm d} Δ ct < U d の物質を電荷移動型と呼んで分ける(Zaanen, Sawatzky, Allen(1985)、Physical Review Letters 55、418)。銅の酸化物は後者に入る。
例題 5.2.5 180° と 90° の超交換を比べる ★★
オリジナル t p d = 1.0 e V t_{\mathrm{pd}}=1.0\,\mathrm{eV} t pd = 1.0 eV 、Δ c t = 3.0 e V \Delta_{\mathrm{ct}}=3.0\,\mathrm{eV} Δ ct = 3.0 eV 、U d = 8 e V U_{\mathrm d}=8\,\mathrm{eV} U d = 8 eV 、U p = 4 e V U_{\mathrm p}=4\,\mathrm{eV} U p = 4 eV 、K p = 0.6 e V K_{\mathrm p}=0.6\,\mathrm{eV} K p = 0.6 eV とする。
(1) 式 (5.2.16) で J 180 ∘ J_{180^\circ} J 18 0 ∘ を求め、2つの項の割合を述べよ。
(2) 式 (5.2.17) で J 90 ∘ J_{90^\circ} J 9 0 ∘ を求め、J 180 ∘ J_{180^\circ} J 18 0 ∘ との比を求めよ。
(3) 同じパラメータで模型を厳密に対角化すると、J 180 ∘ = 0.094 e V J_{180^\circ}=0.094\,\mathrm{eV} J 18 0 ∘ = 0.094 eV 、J 90 ∘ = − 4.1 m e V J_{90^\circ}=-4.1\,\mathrm{meV} J 9 0 ∘ = − 4.1 meV になる。4次の式とのずれの向きと大きさを述べよ。
解答 (1) 4 t p d 4 / Δ c t 2 = 4 / 9 = 0.444 e V 4t_{\mathrm{pd}}^4/\Delta_{\mathrm{ct}}^2=4/9=0.444\,\mathrm{eV} 4 t pd 4 / Δ ct 2 = 4/9 = 0.444 eV 、括弧の中は 1 / 8 + 2 / 10 = 0.125 + 0.200 = 0.325 e V − 1 1/8+2/10=0.125+0.200=0.325\,\mathrm{eV^{-1}} 1/8 + 2/10 = 0.125 + 0.200 = 0.325 e V − 1 で、J 180 ∘ = 0.144 e V J_{180^\circ}=0.144\,\mathrm{eV} J 18 0 ∘ = 0.144 eV 。酸素に正孔が2個入る過程(第2項)が 0.200/0.325 = 62 % を占め、銅に2個入る過程より大きい。これは 2 Δ c t + U p = 10 e V 2\Delta_{\mathrm{ct}}+U_{\mathrm p}=10\,\mathrm{eV} 2 Δ ct + U p = 10 eV でも、経路の数が2倍あるためである。
(2) U ′ = 4 − 1.2 = 2.8 e V U'=4-1.2=2.8\,\mathrm{eV} U ′ = 4 − 1.2 = 2.8 eV 、( 2 Δ c t + U ′ ) 2 − K p 2 = 8.8 2 − 0.36 = 77.08 e V 2 (2\Delta_{\mathrm{ct}}+U')^2-K_{\mathrm p}^2=8.8^2-0.36=77.08\,\mathrm{eV^2} ( 2 Δ ct + U ′ ) 2 − K p 2 = 8. 8 2 − 0.36 = 77.08 e V 2 で、J 90 ∘ = − 0.444 × 1.2 / 77.08 = − 6.92 × 10 − 3 e V J_{90^\circ}=-0.444\times1.2/77.08=-6.92\times10^{-3}\,\mathrm{eV} J 9 0 ∘ = − 0.444 × 1.2/77.08 = − 6.92 × 1 0 − 3 eV 。比は J 90 ∘ / J 180 ∘ = − 0.048 J_{90^\circ}/J_{180^\circ}=-0.048 J 9 0 ∘ / J 18 0 ∘ = − 0.048 で、符号が逆で大きさは約20分の1である。
(3) どちらも厳密な値のほうが小さく、180° で 0.65 倍、90° で 0.59 倍である。t p d / Δ c t = 0.33 t_{\mathrm{pd}}/\Delta_{\mathrm{ct}}=0.33 t pd / Δ ct = 0.33 は図 5.2.2 で4次の式が大きめに出る範囲にあり、6次以降の過程(正孔が酸素へ出たまま往復を重ねる過程)の寄与で J J J が小さくなっている。符号と、90° が 180° よりずっと小さいという大小の関係は、厳密な値でも変わらない。(検算 scripts/verify_ch5-02.py)
伝導電子を介した交換:RKKY 相互作用
金属の中に局在したスピン(希土類の 4f 電子や、薄く混ぜた磁性の不純物)があると、その間の交換は伝導電子が仲立ちする。局在したスピン S ^ j \hat{\bm{S}}_j S ^ j と、その位置の伝導電子のスピン密度 s ^ ( r ) \hat{\bm{s}}(\bm{r}) s ^ ( r ) がJ s d ∑ j S ^ j ⋅ s ^ ( R j ) J_{\mathrm{sd}}\sum_j\hat{\bm{S}}_j\cdot\hat{\bm{s}}(\bm{R}_j) J sd ∑ j S ^ j ⋅ s ^ ( R j ) で結合するとする(J s d J_{\mathrm{sd}} J sd はエネルギー×体積の次元を持つ)。スピン1は伝導電子にとってスピンに依るポテンシャル ± 1 2 J s d S 1 z δ ( 3 ) ( r − R 1 ) \pm\frac12J_{\mathrm{sd}}S_{1z}\delta^{(3)}(\bm{r}-\bm{R}_1) ± 2 1 J sd S 1 z δ ( 3 ) ( r − R 1 ) になり、上向きと下向きの電子の密度を逆向きにずらす。そのずれは第5部第1章 の静的な Lindhard 関数の実空間の形 χ 0 ( r ) \chi_0(r) χ 0 ( r ) (両スピンの和。1つのスピンあたりはその半分)で決まり、位置 R 2 \bm{R}_2 R 2 に誘起されるスピン密度は s z ( R 2 ) = 1 4 J s d χ 0 ( R ) S 1 z s_z(\bm{R}_2)=\frac14J_{\mathrm{sd}}\chi_0(R)S_{1z} s z ( R 2 ) = 4 1 J sd χ 0 ( R ) S 1 z (R = ∣ R 1 − R 2 ∣ R=\lvert \bm{R}_1-\bm{R}_2\rvert R = ∣ R 1 − R 2 ∣ )である。スピン2がこれを感じるエネルギーを回転で不変な形に直すと
H ^ R K K Y = J ( R ) S ^ 1 ⋅ S ^ 2 , J ( R ) = J s d 2 4 χ 0 ( R ) = − J s d 2 m k F 4 2 π 3 ℏ 2 sin x − x cos x x 4 , x = 2 k F R (5.2.18) \hat H_{\mathrm{RKKY}}=J(R)\,\hat{\bm{S}}_1\cdot\hat{\bm{S}}_2,\qquad
J(R)=\frac{J_{\mathrm{sd}}^2}{4}\chi_0(R)=-\frac{J_{\mathrm{sd}}^2mk_{\mathrm F}^4}{2\pi^3\hbar^2}\,\frac{\sin x-x\cos x}{x^4},\qquad x=2k_{\mathrm F}R\tag{5.2.18} H ^ RKKY = J ( R ) S ^ 1 ⋅ S ^ 2 , J ( R ) = 4 J sd 2 χ 0 ( R ) = − 2 π 3 ℏ 2 J sd 2 m k F 4 x 4 sin x − x cos x , x = 2 k F R ( 5.2.18 )
を得る。これを RKKY 相互作用 と呼ぶ(Ruderman, Kittel(1954)、Physical Review 96、99。Kasuya(1956)、Progress of Theoretical Physics 16、45。Yosida(1957)、Physical Review 106、893)。近い距離(x ≪ 1 x\ll1 x ≪ 1 )では ( sin x − x cos x ) / x 4 ≈ 1 / ( 3 x ) > 0 (\sin x-x\cos x)/x^4\approx1/(3x)>0 ( sin x − x cos x ) / x 4 ≈ 1/ ( 3 x ) > 0 で J < 0 J<0 J < 0 (強磁性)、最初に符号が変わるのは tan x = x \tan x=x tan x = x の根 2 k F R = 4.49 2k_{\mathrm F}R=4.49 2 k F R = 4.49 で、遠方では J ∝ cos ( 2 k F R ) / R 3 J\propto\cos(2k_{\mathrm F}R)/R^3 J ∝ cos ( 2 k F R ) / R 3 と符号を変えながら減る。5.2.2 節の縮退摂動の言葉では、これは伝導電子の励起を「高い準位」とした2次の過程で、Friedel 振動と同じく Fermi 面の角(q = 2 k F q=2k_{\mathrm F} q = 2 k F )で振動の波数が決まる。2つのスピンの距離によって強磁性にも反強磁性にもなるので、磁性の不純物がでたらめに散らばった金属ではスピンの向きがそろわずに凍りつく状態(スピングラス)が現れ、磁性層を非磁性の金属の層で挟んだ多層膜では、層の厚さで上下の層の向きが平行と反平行に交互に変わる。
スピン軌道相互作用を戻す:DM 相互作用
ここまで式 (5.2.1) の第3項を 0 にしてきた。これを戻すと、飛び移りの行列要素がスピンに依るようになる。最も簡単な形は、2つのサイトの間の飛び移りを− ∑ σ σ ′ c ^ 1 σ † ( t δ σ σ ′ + i λ ⋅ σ σ σ ′ ) c ^ 2 σ ′ + h . c . -\sum_{\sigma\sigma'}\hat c^\dagger_{1\sigma}\bigl(t\,\delta_{\sigma\sigma'}+\mathrm{i}\bm{\lambda}\cdot\bm{\sigma}_{\sigma\sigma'}\bigr)\hat c_{2\sigma'}+\mathrm{h.c.} − ∑ σ σ ′ c ^ 1 σ † ( t δ σ σ ′ + i λ ⋅ σ σ σ ′ ) c ^ 2 σ ′ + h.c. とすることで、実数のベクトル λ \bm{\lambda} λ がスピン軌道相互作用で生じたスピンの回転を表す(∣ λ ∣ ≪ t \lvert \bm{\lambda}\rvert \ll t ∣ λ ∣ ≪ t )。式 (5.2.8) の2次の摂動を Hubbard の U U U について計算すると
H ^ e f f = 4 U [ ( t 2 − λ 2 ) S ^ 1 ⋅ S ^ 2 + 2 t λ ⋅ ( S ^ 1 × S ^ 2 ) + 2 ( λ ⋅ S ^ 1 ) ( λ ⋅ S ^ 2 ) ] + 定数 (5.2.19) \hat H_{\mathrm{eff}}=\frac{4}{U}\Bigl[(t^2-\lambda^2)\,\hat{\bm{S}}_1\cdot\hat{\bm{S}}_2+2t\,\bm{\lambda}\cdot(\hat{\bm{S}}_1\times\hat{\bm{S}}_2)+2(\bm{\lambda}\cdot\hat{\bm{S}}_1)(\bm{\lambda}\cdot\hat{\bm{S}}_2)\Bigr]+\text{定数}\tag{5.2.19} H ^ eff = U 4 [ ( t 2 − λ 2 ) S ^ 1 ⋅ S ^ 2 + 2 t λ ⋅ ( S ^ 1 × S ^ 2 ) + 2 ( λ ⋅ S ^ 1 ) ( λ ⋅ S ^ 2 ) ] + 定数 ( 5.2.19 )
を得る(λ = ∣ λ ∣ \lambda=\lvert \bm{\lambda}\rvert λ = ∣ λ ∣ 。4 × 4 4\times4 4 × 4 の行列として検算した)。第1項は等方的な交換 J = 4 ( t 2 − λ 2 ) / U J=4(t^2-\lambda^2)/U J = 4 ( t 2 − λ 2 ) / U 、第3項は λ 2 \lambda^2 λ 2 の大きさの異方的な交換で、第2項が新しく現れた反対称な項
H ^ D M = D ⋅ ( S ^ 1 × S ^ 2 ) , D = 8 t λ U , ∣ D ∣ J ≈ 2 λ t (5.2.20) \hat H_{\mathrm{DM}}=\bm{D}\cdot\bigl(\hat{\bm{S}}_1\times\hat{\bm{S}}_2\bigr),\qquad\bm{D}=\frac{8t\bm{\lambda}}{U},\qquad\frac{\lvert \bm{D}\rvert }{J}\approx\frac{2\lambda}{t}\tag{5.2.20} H ^ DM = D ⋅ ( S ^ 1 × S ^ 2 ) , D = U 8 t λ , J ∣ D ∣ ≈ t 2 λ ( 5.2.20 )
である。これを Dzyaloshinskii–Moriya(DM)相互作用 と呼ぶ。Dzyaloshinsky は対称性からこの項の存在を指摘し(Dzyaloshinsky(1958)、Journal of Physics and Chemistry of Solids 4、241)、Moriya は超交換にスピン軌道相互作用を1次で入れてこれを導き、大きさを ∣ D ∣ ∼ ( Δ g / g ) J \lvert \bm{D}\rvert \sim(\Delta g/g)J ∣ D ∣ ∼ ( Δ g / g ) J と見積もった(Moriya(1960)、Physical Review 120、91)。Δ g = g − 2 \Delta g=g-2 Δ g = g − 2 は g g g 因子の2からのずれで、スピン軌道相互作用の強さを測る量である。
D \bm{D} D の向きは、2つのサイトを結ぶ結合の対称性で制限される。スピンは軸性ベクトルなので、空間の反転で向きを変えない。結合の中点に反転の中心があれば、反転で2つのサイトが入れ替わり、S ^ 1 × S ^ 2 → S ^ 2 × S ^ 1 = − S ^ 1 × S ^ 2 \hat{\bm{S}}_1\times\hat{\bm{S}}_2\to\hat{\bm{S}}_2\times\hat{\bm{S}}_1=-\hat{\bm{S}}_1\times\hat{\bm{S}}_2 S ^ 1 × S ^ 2 → S ^ 2 × S ^ 1 = − S ^ 1 × S ^ 2 となる。H ^ \hat H H ^ が反転で不変なら D = − D \bm{D}=-\bm{D} D = − D で、D = 0 \bm{D}=0 D = 0 である。2つのサイトを含む鏡映面があれば、鏡映で軸性ベクトルの面内の成分だけが符号を変えるので、D \bm{D} D はその面に垂直でなければならない。DM 相互作用は、反転の中心を持たない結合にだけ現れる。
図 5.2.4 DM 相互作用による反強磁性のスピンの傾き(傾きの角は誇張して描いた)。等方的な交換 J > 0 J>0 J > 0 のエネルギーは2つのスピンが反平行のとき、D ⋅ ( S ^ 1 × S ^ 2 ) \bm{D}\cdot(\hat{\bm{S}}_1\times\hat{\bm{S}}_2) D ⋅ ( S ^ 1 × S ^ 2 ) のエネルギーは直交するとき最も低い。2つの和を最小にする向きでは各スピンが φ \varphi φ だけ傾き、2つの磁気モーメントの和が 0 でなくなる(弱い強磁性)。
DM 相互作用のエネルギーは2つのスピンが直交するとき最も低く、等方的な交換のエネルギーは平行か反平行のとき最も低い。反強磁性の副格子が2つある結晶で、隣どうしの D \bm{D} D が打ち消し合わずにそろっていると、各スピンが反平行の向きから少し傾き、副格子の磁化の和が残る(図 5.2.4 )。これを弱い強磁性 と呼び、Dzyaloshinsky が赤鉄鉱(α \alpha α -Fe2 _2 2 O3 _3 3 )の弱い自発磁化をこの機構で説明した。
例題 5.2.6 DM 相互作用によるスピンの傾きと、1本の結合では消せること ★★
オリジナル
(1) 長さ S S S の古典的なスピン2つが E = J S 1 ⋅ S 2 + D ⋅ ( S 1 × S 2 ) E=J\bm{S}_1\cdot\bm{S}_2+\bm{D}\cdot(\bm{S}_1\times\bm{S}_2) E = J S 1 ⋅ S 2 + D ⋅ ( S 1 × S 2 ) (J > 0 J>0 J > 0 )で相互作用し、スピンは D \bm{D} D に垂直な面の中にあるとする。各スピンが反平行の向きから φ \varphi φ だけ逆向きに傾くとして、E E E を最小にする φ \varphi φ が tan 2 φ = D / J \tan2\varphi=D/J tan 2 φ = D / J を満たすことを示せ。D / J = 0.02 D/J=0.02 D / J = 0.02 、S = 1 2 S=\frac12 S = 2 1 、g = 2 g=2 g = 2 のとき、φ \varphi φ とスピン1個あたりの正味の磁気モーメントを求めよ。
(2) 量子の2つのスピン 1 2 \frac12 2 1 で H ^ = J S ^ 1 ⋅ S ^ 2 + D ( S ^ 1 x S ^ 2 y − S ^ 1 y S ^ 2 x ) \hat H=J\hat{\bm{S}}_1\cdot\hat{\bm{S}}_2+D(\hat S_{1x}\hat S_{2y}-\hat S_{1y}\hat S_{2x}) H ^ = J S ^ 1 ⋅ S ^ 2 + D ( S ^ 1 x S ^ 2 y − S ^ 1 y S ^ 2 x ) の準位を求めよ。
(3) 式 (5.2.19) の3つの項をすべて残すと、1本の結合の準位は等方的な J ′ = 4 ( t 2 + λ 2 ) / U J'=4(t^2+\lambda^2)/U J ′ = 4 ( t 2 + λ 2 ) / U の Heisenberg 模型と同じになる(検算で確かめた)。このことが、DM 相互作用の効果を見るのに結晶の中で結合が複数要る理由になっていることを説明せよ。
解答 (1) 2つのスピンのなす角は π − 2 φ \pi-2\varphi π − 2 φ である。D \bm{D} D に垂直な面の中で S 1 × S 2 \bm{S}_1\times\bm{S}_2 S 1 × S 2 は D \bm{D} D と平行か反平行で、大きさは S 2 sin ( π − 2 φ ) = S 2 sin 2 φ S^2\sin(\pi-2\varphi)=S^2\sin2\varphi S 2 sin ( π − 2 φ ) = S 2 sin 2 φ 。E E E が下がる向きに傾くとして
E ( φ ) = J S 2 cos ( π − 2 φ ) − D S 2 sin 2 φ = − S 2 ( J cos 2 φ + D sin 2 φ ) . E(\varphi)=JS^2\cos(\pi-2\varphi)-DS^2\sin2\varphi=-S^2\bigl(J\cos2\varphi+D\sin2\varphi\bigr). E ( φ ) = J S 2 cos ( π − 2 φ ) − D S 2 sin 2 φ = − S 2 ( J cos 2 φ + D sin 2 φ ) .
d E d φ = 2 S 2 ( J sin 2 φ − D cos 2 φ ) = 0 \frac{\mathrm{d} E}{\mathrm{d} \varphi}=2S^2(J\sin2\varphi-D\cos2\varphi)=0 d φ d E = 2 S 2 ( J sin 2 φ − D cos 2 φ ) = 0 から tan 2 φ = D / J \tan2\varphi=D/J tan 2 φ = D / J 。D / J = 0.02 D/J=0.02 D / J = 0.02 で φ = 1 2 arctan 0.02 = 0.0100 r a d = 0.573 ∘ \varphi=\frac12\arctan0.02=0.0100\,\mathrm{rad}=0.573^\circ φ = 2 1 arctan 0.02 = 0.0100 rad = 0.57 3 ∘ 。副格子の磁気モーメントは傾いた方向に g μ B S sin φ g\mu_{\mathrm B}S\sin\varphi g μ B S sin φ ずつ残るので、スピン1個あたり 2 × 1 2 × 0.0100 μ B = 0.0100 μ B 2\times\frac12\times0.0100\,\mu_{\mathrm B}=0.0100\,\mu_{\mathrm B} 2 × 2 1 × 0.0100 μ B = 0.0100 μ B で、反強磁性の副格子の磁化の約1 % である。
(2) S ^ 1 x S ^ 2 y − S ^ 1 y S ^ 2 x \hat S_{1x}\hat S_{2y}-\hat S_{1y}\hat S_{2x} S ^ 1 x S ^ 2 y − S ^ 1 y S ^ 2 x は全スピンの z z z 成分 S ^ 1 z + S ^ 2 z \hat S_{1z}+\hat S_{2z} S ^ 1 z + S ^ 2 z と交換するので、S z = ± 1 S_z=\pm1 S z = ± 1 の2状態 ∣ ↑ ↑ ⟩ | \uparrow\uparrow \rangle ∣ ↑↑ ⟩ 、∣ ↓ ↓ ⟩ | \downarrow\downarrow \rangle ∣ ↓↓ ⟩ には作用せず、エネルギーは J / 4 J/4 J /4 のままである。S z = 0 S_z=0 S z = 0 の1重項 ∣ s ⟩ | s \rangle ∣ s ⟩ と3重項 ∣ t 0 ⟩ | t_0 \rangle ∣ t 0 ⟩ の間には ± i D / 2 \pm\mathrm{i} D/2 ± i D /2 の行列要素があり、( − 3 J / 4 i D / 2 − i D / 2 J / 4 ) \begin{pmatrix}-3J/4&\mathrm{i} D/2\\-\mathrm{i} D/2&J/4\end{pmatrix} ( − 3 J /4 − i D /2 i D /2 J /4 ) を対角化してE = − J 4 ± 1 2 J 2 + D 2 E=-\frac J4\pm\frac12\sqrt{J^2+D^2} E = − 4 J ± 2 1 J 2 + D 2 を得る。DM 相互作用は1重項と3重項を混ぜ、D 2 / ( 2 J ) D^2/(2J) D 2 / ( 2 J ) の程度だけ最低の準位を下げる。
(3) 2つのサイトのスピンの座標軸を、λ \bm{\lambda} λ の周りに互いに逆向きに少し回すと、飛び移り t + i λ ⋅ σ t+\mathrm{i}\bm{\lambda}\cdot\bm{\sigma} t + i λ ⋅ σ のスピンの回転の部分を消すことができ、模型は等方的な Hubbard 模型と同じ準位を持つ。1本の結合だけなら DM 相互作用と第3項の異方性は、この座標の取り直しで消える組み合わせになっている。結晶の中ではスピンは複数の結合で隣とつながり、すべての結合の回転を同時に消す座標の取り直しは一般に無い。そのため、DM 相互作用の効果(弱い強磁性やらせん状の秩序)は、結合が並んだ格子の上で初めて現れる。(検算 scripts/verify_ch5-02.py)
ほかの巻では
近似の数理 第13章 :Hubbard の2量体の厳密解と、t / U t/U t / U の展開を超えて基底状態のスピンを決める Lieb の定理が読める
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com