球面調和関数

• 媒介変数表示
• 球面調和関数
• 種々の量子数
• 変形した原子核
  [ver.4]

種々の軌道角運動量と磁気量子数

l=1,m=-1,+1

球面調和関数Y[lm]のl,mは，量子力学の軌道角運動量と磁気量子数です．
ここでは，小さいlに対するグラフを幾つか表示してみます．やり方は，前項
と同様，媒介変数を用いてx,y,zを表す方法です．

mが奇数の場合は球面調和関数はexp(-imφ)の項の為に虚部を持ちます．
m=-1,1ののこの関数は次のようになります．

 Y[1,-1](t,p) =  sqrt(3/8pi) sin(t) exp(-i p) Y[1, 1](t,p) = -sqrt(3/8pi) sin(t) exp( i p) 

|Y|^2はYとその複素共役をかければよく，上の2つの式はいづれも同じY
(θ)=3/8π sin^2θ となります．なお，球面調和関数の複素共役は簡単に求
めることができ， (Y[l,m])^* = (-1)^m Y[l,-m] で与えられます．
つまり，mの符号を変え，さらにmが奇数の場合は全体の符号を変えます．

 gnuplot> set parametric  dummy variable is t for curves, u/v for surfaces gnuplot> set angle degree gnuplot> set urange [0:360] gnuplot> set vrange [-90:90] gnuplot> set isosample 36,18 gnuplot> set ticslevel 0 gnuplot> set size 0.65,1.0 gnuplot> a=3.0/(8*pi) gnuplot> fx(u,v)=cos(u)*cos(v) gnuplot> fy(u,v)=sin(u)*cos(v) gnuplot> fz(v)=sin(v) gnuplot> g(v)=cos(v)*cos(v) gnuplot> splot a*g(v)*fx(u,v),a*g(v)*fy(u,v),a*g(v)*fz(v) 

左がY[11]を自乗したグラフで，z軸を主軸にしたトーラスとなります．但し
内側に穴はありません．関数はsin(θ)を含みますが，角度θをvで
表すとcos(v)になるので， g(v)=cos(v)*cos(v)と関数定義して
います．

l=2,m=-2,-1,0,+1,+2

l=2のときは，磁気量子数mが -2〜2 までの5つの関数があります．それぞれの
関数は以下で与えられます．

 Y[2,-2](t,p) =  sqrt(15/32pi) sin(t)sin(t) exp(-2i p) Y[2,-1](t,p) =  sqrt(15/8pi) sin(t)cos(t) exp(-i p) Y[2, 0](t,p) =  sqrt(5/16pi) (3cos(t)cos(t)-1) Y[2, 1](t,p) = -sqrt(15/8pi) sin(t)cos(t) exp( i p) Y[2, 2](t,p) =  sqrt(15/32pi) sin(t)sin(t) exp( 2i p) 

自乗にして表示すると|m|が同じものは同じグラフになりますので，m=0,1,2の
3つの場合について表示します．

m=0

 gnuplot> a= 5.0/(16*pi) gnuplot> g(v)= (3*sin(v)*sin(v)-1)**2 gnuplot> splot a*g(v)*fx(u,v),a*g(v)*fy(u,v),a*g(v)*fz(v) 

m=1

 gnuplot> a=15.0/( 8*pi) gnuplot> g(v)= (sin(v)*cos(v))**2 gnuplot> splot a*g(v)*fx(u,v),a*g(v)*fy(u,v),a*g(v)*fz(v) 

m=2

 gnuplot> a=15.0/(32*pi) gnuplot> g(v)= (cos(v)*cos(v))**2 gnuplot> splot a*g(v)*fx(u,v),a*g(v)*fy(u,v),a*g(v)*fz(v)