1. 調和振動子と$SU(N)$代数 vol.4

1.5 高次元調和振動子

前節までの解析を3次元調和振動子に拡張するのは比較的容易である。3次元振動子のハミルトニアンは

$$H = \om \left(     a^\dagger_{1} a_1 + a^\dagger_{2} a_2 + a^\dagger_{3} a_3 + \frac{3}{2}    \right)    = \om \left(    a^\dagger_{i} a_i + \frac{3}{2}    \right)    \tag{1.47}$$

で与えられる。ただし、$i = 1,2,3$である。数演算子の和は$n = n_1 + n_2 + n_3 = a^\dagger_{i} a_i$で表される。$n$を定数とすると、縮退度は${1\over 2} (n+1)(n+2)$となる。縮退状態をつなぐ演算子には下の画像で示された$a^\dagger_2 a_1$, $a^\dagger_3 a_1$, $a^\dagger_3 a_2$とこれらに共役な$a^\dagger_1 a_2$, $a^\dagger_1 a_3$, $a^\dagger_2 a_3$が含まれる。

　演算子による閉じた代数をつくるにはこれら６つの演算子にさらに他の演算子を追加する必要がある。6つのうち2つの演算子の交換関係を計算すると

$$\begin{eqnarray}    [ a^\dagger_1 a_2 , a^\dagger_2 a_1 ] &=& a^\dagger_1 [ a_2 , a^\dagger_2 a_1 ] +    [ a^\dagger_1 , a^\dagger_2 a_1 ] a_2 \nonumber \\    &=& a^\dagger_1 a_1  - a^\dagger_2 a_2 \nonumber \\    &=& n_1  -  n_2     \tag{1.48}\end{eqnarray}$$

となる。よって、演算子のリストに$n_1 - n_2$を含めなければならない。同様に$( n_2 - n_3 )$と$( n_3 - n_1 )$も必要になる。2つ目の演算子は$( n_3 - n_1 ) = -( n_1 - n_2 ) - ( n_2 - n_3 )$と表せる。よって、閉じた演算子代数をつくるには2つの演算子を追加する必要がある。計8個の演算子が成す代数は$SU(3)$代数となる。以下では、この点を明らかにしていく。

SU(3)の表現

　$SU(2)$の場合と同様に、$SU(3)$群の要素は$g = e^{i t}$とパラメータ表示される。ただし、$t$はトレースがゼロの$(3 \times 3)$エルミート行列である。$3\times 3$行列は９つの独立な行列要素をもつがトレースゼロとなる条件があるので、8つの独立パラメータで$SU(3)$の要素を表せると予測できる。よって、行列$t$は$t = t^a \th^a$ ($a=1,2, \cdots , 8)$とパラメータ表示できる。ただし、$t^a$はトレースゼロのエルミート行列の基底を表し、$\th^a$は実パラメータである。（また、演算子$a^\dagger_i a_j$ $(i\neq j)$のパラメータとして３つの複素数、対角成分の生成子に対応するパラメータとして2つの実数が必要となるのでこれらを用いて$t$をパラメータ表示することもできる。）慣例として、この基底$t^a$は以下で与えられるゲルマン行列$\la^a$を用いて$t^a = \frac{\la^a}{2}$と表される。

$$\begin{eqnarray}    &&    \la^1 =    \left(      \begin{array}{ccc}        0 & 1 & 0 \\        1 & 0 & 0 \\        0 & 0 & 0 \\      \end{array}    \right)~~~    \la^2 =    \left(      \begin{array}{ccc}        0 & -i & 0 \\        i & 0 & 0 \\        0 & 0 & 0 \\      \end{array}    \right)~~~    \la^3 =    \left(      \begin{array}{ccc}        1 & 0 & 0 \\        0 & -1 & 0 \\        0 & 0 & 0 \\      \end{array}    \right)    \nonumber \\    &&    \la^4 =    \left(      \begin{array}{ccc}        0 & 0 & 1 \\        0 & 0 & 0 \\        1 & 0 & 0 \\      \end{array}    \right)~~~    \la^5 =    \left(      \begin{array}{ccc}        0 & 0 & -i \\        0 & 0 & 0 \\        i & 0 & 0 \\      \end{array}    \right)~~~    \la^6 =    \left(      \begin{array}{ccc}        0 & 0 & 0 \\        0 & 0 & 1 \\        0 & 1 & 0 \\      \end{array}    \right)    \nonumber \\    &&    \la^7 =    \left(      \begin{array}{ccc}        0 & 0 & 0 \\        0 & 0 & -i \\        0 & i & 0 \\      \end{array}    \right)~~~    \la^8 = \frac{1}{\sqrt{3}}    \left(      \begin{array}{ccc}        1 & 0 & 0 \\        0 & 1 & 0 \\        0 & 0 & -2 \\      \end{array}    \right)    \tag{1.49} \end{eqnarray}$$

基底$t^a = \frac{\la^a}{2}$は$\Tr (t^a t^b ) = \hf \del^{ab}$と規格化されている。交換子$[ t^a , t^b ]$は反エルミートである。よって、交換子は$i \times \mbox{(エルミート行列)}$で表せる。交換子もトレースゼロとなるので$t^a$で展開することができ、次のように書ける。

$$[ t^a , t^b ] = i f^{abc} t^c     \tag{1.50}$$

$f^{abc}$は構造定数と呼ばれる。構造定数は定数の集合であり、その具体的な値は今のところ必要ない。交換関係(1.50)は$SU(3)$のリー代数をなす。

　ここで、演算子

$$\hat{T}^a = a^\dagger_i t_{ij}^{a} a_j  = a^\dagger t^{a} a    \tag{1.51}$$

を考えよう。ただし、$t_{ij}^{a}$は$( 3\times 3)$行列$t^a$ $(i,j = 1,2,3; a= 1,2,\cdots , 8)$の行列成分を表す。以下では、これらの演算子が３次元調和振動子の縮退状態をつなぐ演算子に対応していることを見ていく。まず、$\hat{T}^a$はエルミートである。これは関係式

$$\begin{eqnarray}    ( \hat{T}^a )^\dagger &= & a^\dagger_j \left( t^a \right)_{ij}^{*} a_i    \nonumber \\    & = & a^\dagger_j  t^{a}_{ji} a_i = \hat{T}^a    \tag{1.52}\end{eqnarray}$$

から分かる。ただし、$t^a$のエルミート性$\left( t^a \right)_{ij}^{*} = t^{a}_{ji}$を用いた。つぎに、$\hat{T}^a$は$SU(3)$代数に従う。実際に

$$\begin{eqnarray}    [ \hat{T}^a , \hat{T}^b ] &=& [ a^\dagger_i t^{a}_{ij} a_j , a^\dagger_k t_{kl}^{b} a_l ]    \nonumber \\    &=& a^\dagger_i t_{ij}^{a} [ a_j ,  a^\dagger_k t_{kl}^{b} a_l ] +    [ a^\dagger_i  , a^\dagger_k t_{kl}^{b} a_l ] t^{a}_{ij} a_j    \nonumber \\    &=& a^\dagger_i t_{ij}^{a} \del_{jk} t_{kl}^{b} a_l  - \del_{il} a^\dagger_k t_{kl}^{b} t^{a}_{ij} a_j    \nonumber \\    &=& a^\dagger_i t_{ij}^{a} t_{jl}^{b} a_l -  a^\dagger_k t_{ki}^{b} t^{a}_{ij} a_j    \nonumber \\    &=& a^\dagger_i \left(  t_{ij}^{a} t_{jl}^{b} - t_{ij}^{b} t^{a}_{jl} \right) a_l    \nonumber \\    &=& a^\dagger [ t^a , t^b ] a = i f^{abc} a^\dagger t^c a = i f^{abc} \hat{T}^c     \tag{1.53}\end{eqnarray}$$

と確認できる。最後に、$[ \hat{T}^a , a^\dagger_i a_i ]$はゼロとなる。

$$[ \hat{T}^a , a^\dagger {\bf 1} a ] = a^\dagger [ t^a , {\bf 1} ] a = 0    \tag{1.54}$$

ただし、${\bf 1}$は$(3 \times 3)$単位行列を表す。式(1.54)より明らかに$[ \hat{T}^a  , H ] = 0$となる。これは総数$n = n_1 + n_2 + n_3 = a^\dagger_i a_i$は$\hat{T}^a$の作用のもとで不変であることを意味する。言い換えると、演算子$\hat{T}^a$の作用によって物理系が縮退状態の外部に持ち出されることはない。

　これらの議論から、3次元等方調和振動子の縮退状態は$SU(3)$表現を成すことが分かった。この縮退状態を$| \al \ket$で表示する。ただし、$\al = 1, 2, \cdots , N={1\over 2}(n+1)(n+2)$である。$| \al \ket$の一般形は

$${(a^\dagger_1 )^{n_1}\over \sqrt{n_1!}}\,{(a^\dagger_2 )^{n_2}\over \sqrt{n_2!}}\, {(a^\dagger_3 )^{n_3}\over \sqrt{n_3!}}~ | 0 \ket, \hskip .3in n_1 + n_2 + n_3 = n \tag{1.55}$$

と書けることに注意しよう。これより、$(N \times N)$行列$\bra \al | \hat{T}^a | \bt \ket$を考えることができ、これらは$SU(3)$代数に従う。よって、$\bra \al | \hat{T}^a | \bt \ket$は$SU(3)$代数の具体的な行列表現を与える。

　この章で紹介した形式化はより一般に$SU(N)$代数にもそのまま適用できる。その構成の概略は以下の通り。まず、$T^{a} = a^\dagger t^{a} a$の形の演算子を考える。ここで、$t^a$ ($a = 1,2, \cdots , N^2 -1$)はトレースがゼロの$(N \times N )$エルミート行列の基底を成す。また、$a^\dagger_i$と$a_i$ ($i = 1,2, \cdots N$)は$N$次元調和振動子の生成・消滅演算子にあたる。$N$次元調和振動子における縮退状態の遷移は演算子$T^{a} = a^\dagger t^{a} a$によって記述され、これらの演算子は$SU(N)$代数に従う。

　最後に、この形式化によって$SU(N)$ ($N \ge 3$) のすべてのユニタリー既約表現が得られる訳ではないことに留意する必要がある。縮退状態の構成からすぐに分かるように、ここで得られる表現は対称表現にすぎない。$SU(2)$の場合は、すべての表現が対称表現であるので調和振動子の解析により$SU(2)$代数のすべてのユニタリー既約表現を求めることができるが、これはどちらかといえば例外的なケースである。$N \ge 3$の場合、この手法によって構成できるのは$SU(N)$代数の対称表現だけである。

1.6 物理への応用例

これまで調和振動子の一般的な形式化について議論してきた。では、この結果を直接応用できる物理系はあるのだろうか。その一例として、原子核の殻模型（シェルモデル）を挙げることができる。このモデルでは原子核内の核子（陽子あるいは中性子）をそれぞれ中心力ポテンシャルの中を動く一つの粒子とみなす。このポテンシャルは平均場近似の方法で得られると理解できる。核子の間には強い相互作用が働くので核子の動きを１粒子の動力学とみなす近似は大胆なものである。しかし、強い相互作用の効果が平均場近似によってある程度取り込まれていると期待することもできる。そのようなポテンシャルの良い近似として3次元調和振動子ポテンシャル

$$V(r) = - V_0 + {1\over 2} M \omega^2 (x^2 + y^2 + z^2) \tag{1.56}$$

を採用できる。これまでの解析から核子の縮退状態は$SU(3)$表現で分類され、その縮退度は$N = {1\over 2} (n+1)(n+2)$で与えられることが分かる。最初のいくつかの縮退状態は縮退度$1, 3, 6, 10$をもつ。よって、スピンの内部自由度も考慮すると、殻（シェル）を満たす粒子の数は$2$, $2+ 6 = 8$, $2+ 6 + 12 =20$, $2 + 6+ 12+ 20 =40$などとなる。特に安定な原子核を構成するのに必要な核子の数は魔法数とよばれ、これらは$2, 8, 20, 28, 50, 82, 126$の値をとる。3次元調和振動子ポテンシャルを用いた殻模型はこれらの魔法数のうち最初の３つを再現することが分かる。$n$が大きくなるとこの一致は見られない。その一つの理由はスピン軌道相互作用の核子への影響が増大するためである。スピン軌道相互作用の効果を含めると、殻模型を用いて原子核の魔法数をうまく説明できることが知られている。