next up previous contents
: 励起状態とエネルギーギャップ : BCS理論 : BCS理論   目次

電子対形成とBCS基底状態

電子間の相互作用がフェルミ面付近では引力になっていると仮定する.ハミルトニアンとして,
$\displaystyle H = \sum_{\mbox{\boldmath$k$}, \sigma}\epsilon_{\mbox{\boldmath$k...
...gma'}^{\dag }
c_{\mbox{\boldmath$k$}', \sigma'}c_{\mbox{\boldmath$k$}, \sigma},$     (7.11.81)

をとる. $c_{\mbox{\boldmath$k$}, \sigma}$ $c_{\mbox{\boldmath$k$}, \sigma}^{\dag }$は,それぞれ運動量$k$,スピン$\sigma $$\uparrow$$\downarrow$)の電子の消滅および生成演算子を表す.すなわち,
$\displaystyle \ [ c_{\mbox{\boldmath$k$}, \sigma}, c_{\mbox{\boldmath$k$}', \sigma'}^{\dag }\ ]$ $\textstyle =$ $\displaystyle c_{\mbox{\boldmath$k$}, \sigma}c_{\mbox{\boldmath$k$}', \sigma'}^...
...ma'} = \delta_{\mbox{\boldmath$k$}\mbox{\boldmath$k$}'}\delta_{\sigma \sigma'},$  
$\displaystyle \ [ c_{\mbox{\boldmath$k$}, \sigma}, c_{\mbox{\boldmath$k$}', \sigma'} \ ]$ $\textstyle =$ $\displaystyle \ [ c_{\mbox{\boldmath$k$}, \sigma}^{\dag }, c_{\mbox{\boldmath$k$}', \sigma'}^{\dag } \ ] = 0,$  

を満たす.

基底状態の波動関数として,

$\displaystyle \vert\Psi_{\rm G}\rangle = \prod_{\mbox{\boldmath$k$}} \left ( u_...
...row}^{\dag }c_{-{\mbox{\boldmath$k$}\downarrow}} \right ) \vert\Psi_{0}\rangle,$     (7.11.82)

をとろう.ここで $\vert\Psi_{0}\rangle$は真空状態を表す. $u_{\mbox{\boldmath$k$}}$ $v_{\mbox{\boldmath$k$}}$は,
$\displaystyle u_{\mbox{\boldmath$k$}}^{2} + v_{\mbox{\boldmath$k$}}^{2} = 1,$     (7.11.83)

の関係を満たす.(7.11.2)式では全電子数は一定でないので,その平均値がある一定値$N$に等しいとする.すなわち,
$\displaystyle \langle \Psi_{\rm G}\vert\sum_{\mbox{\boldmath$k$}, \sigma} c_{\m...
...$}, \sigma}^{\dag }c_{\mbox{\boldmath$k$}, \sigma}\vert\Psi_{\rm G}\rangle = N,$     (7.11.84)

とおく. $\vert\Psi_{\rm G}\rangle$が規格化されていることは,
$\displaystyle \langle \Psi_{\rm G}\vert\Psi_{\rm G}\rangle$ $\textstyle =$ $\displaystyle \langle\Psi_{0}\vert\prod_{\mbox{\boldmath$k$}}(u_{\mbox{\boldmat...
...parrow}^{\dag }c_{-\mbox{\boldmath$k$}'\downarrow}^{\dag })\vert\Psi_{0}\rangle$  
  $\textstyle =$ $\displaystyle \langle \Psi_{0} \vert \prod_{\mbox{\boldmath$k$}} (u_{\mbox{\bol...
...rrow}^{\dag }c_{-\mbox{\boldmath$k$}\downarrow}^{\dag }) \vert \Psi_{0} \rangle$  
  $\textstyle =$ $\displaystyle \prod_{\mbox{\boldmath$k$}} (u_{\mbox{\boldmath$k$}}^{2} + v_{\mbox{\boldmath$k$}}^{2} ) = 1,$ (7.11.85)

のように確かめられる.

(7.11.2)を用いて,エネルギーの期待値,

$\displaystyle W_{\rm G} = \langle \Psi_{\rm G} \vert H \vert \Psi_{\rm G} \rangle,$     (7.11.86)

を計算してみよう.

$\epsilon_{\mbox{\boldmath$k$}}$のかかる項は,(7.11.4)の条件から,

$\displaystyle \langle \Psi_{\rm G} \vert \sum_{\mbox{\boldmath$k$}\sigma}c_{\mb...
..._{\rm G} \rangle = 2\sum_{\mbox{\boldmath$k$}} v_{\mbox{\boldmath$k$}}^{2} = N,$      

で表される.$H$の第2項の期待値では,$k$ = $-\mbox{\boldmath$k$}'$$\sigma $ = $-\sigma'$という項と,$q$ = 0,$k$$\sigma $ $\mbox{\boldmath$k$}'\sigma'$は任意という項と, $\mbox{\boldmath$k$}'$ = $\mbox{\boldmath$k$} + \mbox{\boldmath$q$}$$\sigma $ = $\sigma'$の項だけが残る.後ろの2項は,$H$の第1項と同じ形をしているので,相互作用は小さいと仮定して無視する.こうして,
$\displaystyle W_{\rm G} = 2\sum_{\mbox{\boldmath$k$}} \epsilon_{\mbox{\boldmath...
...th$k$}}v_{\mbox{\boldmath$k$}}u_{\mbox{\boldmath$k$}'}v_{\mbox{\boldmath$k$}'}.$     (7.11.87)

(7.11.7)式を(7.11.2)式の条件下で最小にする.

$\displaystyle {\partial W_{\rm G} \over{\partial v_{\mbox{\boldmath$k$}}}}$ $\textstyle =$ $\displaystyle 4\epsilon_{\mbox{\boldmath$k$}}v_{\mbox{\boldmath$k$}} + \sum_{\m...
...$k$}}}} \sqrt{1 - v_{\mbox{\boldmath$k$}}^{2}} v_{\mbox{\boldmath$k$}} \right )$  
  $\textstyle =$ $\displaystyle 4\epsilon_{\mbox{\boldmath$k$}}v_{\mbox{\boldmath$k$}} + 2 \left ...
...$k$}\mbox{\boldmath$k$}'} u_{\mbox{\boldmath$k$}'}v_{\mbox{\boldmath$k$}'} = 0,$  

であるから,
$\displaystyle 2\epsilon_{\mbox{\boldmath$k$}}u_{\mbox{\boldmath$k$}}v_{\mbox{\b...
...{\boldmath$k$})(u_{\mbox{\boldmath$k$}}^{2} - v_{\mbox{\boldmath$k$}}^{2}) = 0,$     (7.11.88)


$\displaystyle \Delta(\mbox{\boldmath$k$}) = - \sum_{\mbox{\boldmath$k$}'}V_{\mb...
...dmath$k$}\mbox{\boldmath$k$}'}u_{\mbox{\boldmath$k$}'}v_{\mbox{\boldmath$k$}'},$     (7.11.89)

を得る.(7.11.8)式と(7.11.2)式を, $u_{\mbox{\boldmath$k$}}$ $v_{\mbox{\boldmath$k$}}$に関して解くと,
$\displaystyle u_{\mbox{\boldmath$k$}}^{2} = {1 \over{2}}\left ( 1 + {\epsilon_{...
...( 1 - {\epsilon_{\mbox{\boldmath$k$}} \over{E_{\mbox{\boldmath$k$}}}} \right ),$     (7.11.90)


$\displaystyle E_{\mbox{\boldmath$k$}} = \sqrt{ \epsilon_{\mbox{\boldmath$k$}}^{2} + \Delta(\mbox{\boldmath$k$})^{2} },$     (7.11.91)

となる.これらを用いて(7.11.9)式を書き直すと,
$\displaystyle \Delta(\mbox{\boldmath$k$}) = - \sum_{\mbox{\boldmath$k$}'} V_{\m...
...qrt{ \epsilon_{\mbox{\boldmath$k$}'}^{2} + \Delta(\mbox{\boldmath$k$}')^{2}}}},$     (7.11.92)

となる.

次に相互作用について,

$\displaystyle V_{\mbox{\boldmath$k$}\mbox{\boldmath$k$}'} = \left\{
\begin{arra...
...t\epsilon_{\mbox{\boldmath$k$}'}\vert > \hbar\omega_{\rm c}
\end{array} \right.$     (7.11.93)

の仮定をする.ここで$V > 0$とする.(7.11.13)式を仮定すれば,(7.11.9)で与えられる $\Delta(\mbox{\boldmath$k$})$$k$に依存しないので,
$\displaystyle \Delta(\mbox{\boldmath$k$}) = \left\{
\begin{array}{ll}
\displays...
...rt\epsilon_{\mbox{\boldmath$k$}}\vert > \hbar\omega_{\rm c}
\end{array} \right.$     (7.11.94)


$\displaystyle \Delta_{0} = -\sum_{\mbox{\boldmath$k$}'} V_{\mbox{\boldmath$k$}\...
...elta_{0} \over{2 \sqrt{\epsilon_{\mbox{\boldmath$k$}'}^{2} + \Delta_{0}^{2}}}},$     (7.11.95)


$\displaystyle 1 = {V \over{2}}\sum_{\mbox{\boldmath$k$}}\,' {1 \over{\sqrt{\epsilon_{\mbox{\boldmath$k$}}^{2} + \Delta_{0}^{2} }}},$     (7.11.96)

となる. $\sum_{\mbox{\boldmath$k$}}\,'$は, $\vert\epsilon_{\mbox{\boldmath$k$}}\vert \leq \hbar\omega_{\rm c}$という状態についての和をとるという意味である.$k$についての和を $\epsilon_{\mbox{\boldmath$k$}}$についての積分に置き換える. $\epsilon_{\mbox{\boldmath$k$}} = \epsilon$の状態密度を$N(\omega)$とすれば,
$\displaystyle \sum_{\mbox{\boldmath$k$}}\ \ \ \ \ \longrightarrow\ \ \ \ \ \int N(\epsilon){\rm d}\epsilon,$      

の置き換えをすればよい.$N(\epsilon)$はフェルミ面$\epsilon = 0$の近くでは緩やかな変化をしていると考えられるので,(7.11.16)を次式で近似する.
$\displaystyle 1 = {V \over{2}}N(0)V \int_{-\hbar\omega_{\rm c}}^{\hbar\omega_{\rm c}} {{\rm d}\epsilon \over{\sqrt{\epsilon^{2} + \Delta_{0}^{2} }}}.$      

したがって,
$\displaystyle \Delta_{0} = {\hbar\omega_{\rm c} \over{\sinh \left ( {1 \over{N(0)V}} \right ) }},$     (7.11.97)

を得る.

(7.11.10)式の $u_{\mbox{\boldmath$k$}}$ $v_{\mbox{\boldmath$k$}}$を用いて,(7.11.7)の基底状態を計算すると,

$\displaystyle W_{\rm G} = W_{\rm G}^{0} - 2{N(0)\hbar^{2}\omega_{\rm c}^{2} \over{\exp \left ( {2 \over{N(0)V}} \right ) - 1 }},$     (7.11.98)

となる.$W_{\rm G}^{0}$は, $\Delta_{0} = 0$のときの(常伝導状態の)エネルギーである.$N(0)$ $\approx$ 1/$E_{\rm F}$より,$N(0)V$ $\ll$ 1であるので,(7.11.17),(7.11.18)式は,
$\displaystyle \Delta_{0} = 2\hbar\omega_{\rm c}\exp \left ( {1 \over{N(0)V}} \right ),$     (7.11.99)


$\displaystyle W_{\rm G} = W_{\rm G}^{0} - 2N(0)\hbar^{2}\omega_{\rm c}^{2}\exp \left ( -{2 \over{N(0)V}} \right ),$     (7.11.100)

のように近似される.

(7.11.10)より,

$\displaystyle v_{\mbox{\boldmath$k$}}^{2} = {1 \over{2}} \left ( 1 - {\epsilon_...
...}} \over{\sqrt{\epsilon_{\mbox{\boldmath$k$}}^{2} + \Delta_{0}^{2}}}} \right ),$     (7.11.101)

である. $v_{\mbox{\boldmath$k$}}^{2}$は,フェルミエネルギーより十分低い所で1,十分高い所でゼロであり,有限温度での常伝導金属のフェルミ関数に似ている.

(7.11.19),(7.11.20)式から,

$\displaystyle W_{\rm G} - W_{\rm G}^{0} = -2N(0)\hbar^{2}\omega_{\rm c}^{2}\exp \left ( -{2 \over{N(0)V}} \right ) = -{1 \over{2}}N(0)\Delta_{0}^{2},$     (7.11.102)

となるが,これが$T = 0$における凝縮エネルギーである.このエネルギーが $\displaystyle {H_{\rm c}(0)^{2} \over{8\pi}}$に等しくなるはずである.
next up previous contents
: 励起状態とエネルギーギャップ : BCS理論 : BCS理論   目次
Masashige Onoda 平成18年4月7日