Processing math: 100%

Shironetsu Blog

@shironetsuのブログ

球面調和関数のレシピ - 角運動量の合成から

 球面上では、自分の方程式の基本解が球面調和関数になることにヤルダは気付いた。以前、地震学の講義でいちどだけお目にかかったことのある種類の波形だ。球面全体で成り立つどんな複雑な解でも、それぞれの調和振動にその寄与をあらわす適切な係数をかけて足しあわせることで表現できる。
グレッグ・イーガン『クロックワーク・ロケット』)

 物理学では水素原子のSchrödinger方程式を解くとき, 角度方向と動径方向の変数分離をして球面上の微分作用素so(3)su(2)代数に従うことを利用することで球面調和関数に出会うようになっている(このあたりの歴史を自分は知らない). 物理っぽい(2次元)球面調和関数の導出法として, 他にWignerのD行列の次元を落とす方法もある. ここではスピンの合成を利用して導いてみる.

個のスピン1系の合成によってスピン系を得る方法は1通りしかない.

スピン1系はこのように表記する.

|11=|+,   |10=|+,   |11=|

まずスピン1系を個繋げてm=状態を用意する.

|=|+|+|+

下降演算子を作用させて適当に規格化定数をかければ全てのスピン状態が得られる

|m=(+m)!(2)!(m)! Jm|

スピン1の状態たちは同じ空間に住んでいるものとして順序を考慮しなければ

Jm|=p+q+r=, pr=mA(p,q,r)|+p(2|0q(2|)r

と書ける. 0と-にかかっている2,2はスピン1系への下降演算子の作用で出てくるお釣り, Aはスピン1系へ下降演算子を作用させることでp,q,rに辿り着く道の総数.

Aはこのように考える. 1からまでの|+状態をここでは区別して, まず下降演算子を1回だけ作用させるq個を選ぶ. 残った中から2回作用させるr個を選ぶ. そのr個については重複を許して作用させる順番を決める. 最後にそのr個の順序をキャンセルするために2rで割る. つまり,

A(p,q,r)=CqqCr(q+2r)!2r=!(q+2r)!p!q!r!2r

よって,

Jm|=q!(m)!2q(+mq2)!q!(mq2)! |++mq2|0q|mq2

変数qの範囲は冪が非負整数となるようなすべての整数. 規格化定数をかけて,

|m=!(+m)!(m)!2(21)!!q2q(+mq2)!q!(mq2)! |++mq2|0q|mq2

となる. さて, スピン1系は

|+=|xi|y2 ,  |0=|z ,  |=|xi|y2

であった. 球面上の関数として

xiy2=sinθeiφ2 ,  z=cosθ ,  xiy2=sinθeiφ2

も同じ代数に従い, 3つに共通の定数倍の違いを除いて=1の球面調和関数に等しい. 一般のについて, その規格化定数は球面上の積分内積として定義されるべきであり,|に対して

S2dS(x2+y22)=πθ=02πφ=0sinθdθdφ(sin2θ2)=4π!(2+1)!!

の逆数の平方根. 結局,

Ym=(2+1)(+m)!(m)!4πq(xiy)+mq2zq(xiy)mq22q(+mq2)!q!(mq2)!=(2+1)(+m)!(m)!4π(+m)/2k=max{0,m}(1)k(x+iy)kz+m2k(xiy)km22kmk!(+m2k)!(km)!=(2+1)(+m)!(m)!4π(+m)/2k=max{0,m}(1)ksin2kmθcos+m2kθ eimφ22kmk!(+m2k)!(km)!

という表式を得る(ほんとうは4πの因子は面積要素のほうに押し付けたい). x,y,z座標での表示が先に出てきて便利. これを使って正20面体対称性を持った多項式の形を見ていきたい(続).

(追記)
 確認のため利用したMaximaで使われている定義とはmが奇数の時だけ-1倍違った. この-1倍というのは下降演算子の作用に由来している. これをCondon - Shortley phaseというらしい.
Condon-Shortley Phase -- from Wolfram MathWorld
 2人は角運動量固有状態間の位相についての取り決めを定めた有名な本の著者. Clebsch - Gordan係数表はだいたいCondon - Shortley conventionに従っている. それを思うと球面調和関数の定義がいまいちよく統一されていないのはよく分からない.