前の単元 <<<
量子力学Ⅰ
>>> 次の単元
復習:x,y,z 座標で解いた3次元調和振動子†
バネ定数 K=mω2 (ω は系の固有角振動数)
の3次元調和振動子に対する時間を含まないシュレーディンガー方程式を
x,y,z について変数分離して解けば、その解は φ(x,y,z)=X(x)Y(y)Z(z)
の形となり、X(x),Y(y),Z(z) はそれぞれ1次元調和振動子に対する解となることを
量子力学Ⅰ/調和振動子#y596d643 で学んだ。
すなわち r0=h−/mω およびエルミート多項式 Hn を使って、
φnxnynz(x,y,z)=Hnx(x/r0)Hny(y/r0)Hnz(z/r0)e−2r02r2
対応するエネルギーは N=nx+ny+nz を用いて、
εnxnynz=h−ω(nx+ny+nz+3/2)=h−ω(N+3/2)
したがって、エネルギー量子数 N に対する縮退度は次のようになる。
N=nx+ny+nz | (nxnynz) | 縮退度 |
0 | (0 0 0) | 1 |
1 | (1 0 0), (0 1 0), (0 0 1) | 3 |
2 | (2 0 0), (0 2 0), (0 0 2) | 6 |
(1 1 0), (1 0 1), (0 1 1) |
3 | (3 0 0), (0 3 0), (0 0 3) | 10 |
(2 1 0), (0 2 1), (1 0 2) |
(1 2 0), (0 1 2), (2 0 1) |
(1 1 1) |
球座標で解いた3次元調和振動子†
一方、この問題は球対称なポテンシャル
V(r,θ,ϕ)=21Kr2
の中での運動であるから、球座標で展開して解けば
φnlm′(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ)
の形の解が得られる。このとき、
ε=εnl
と表せ、エネルギーは m の値によらない。
状態 | l | m | 縮退度 |
s | 0 | 0 | 1 |
p | 1 | −1,0,+1 | 3 |
d | 2 | −2,−1,0,1,2 | 5 |
f | 3 | −3,−2,−1,0,1,2,3 | 7 |
ただし、水素原子のところでも見たとおり、ポテンシャルの形状によっては m
に対する縮退以外にも、異なる l を持つ状態が縮退することがある。
両者の関係は?†
次に見るとおり、
- N=0 の縮退度は 1 → l=0 (縮退度 1 の s 状態)
- N=1 の縮退度は 3 → l=1 (縮退度 3 の p 状態)
- N=2 の縮退度は 6 → l=2 (縮退度 5 の d 状態) と l=0 (縮退度 1 の s 状態) が縮退
- N=3 の縮退度は 10 → l=3 (縮退度 7 の f 状態) と l=1 (縮退度 3 の p 状態) が縮退
にそれぞれ対応している。
演習:球座標による解の縮退†
x,y,z 座標で変数分離して得られる解と、
r,θ,ϕ 座標で変数分離して得られる解との関係を明らかにしたい。
(1) 3次元調和振動子に対する動径方向の方程式は、
r/r0=ξ と書き換え、rR(r)=X(ξ) と置けば、
X′′(ξ)=[ξ2−h−ω2ε+ξ2l(l+1)]X(ξ)
となることを示せ。(l=0 のとき、この式は1次元調和振動子の式と一致する)
(2) dξ2d2(ξke−ξ2/2)=(ξ2−(2k+1)+ξ2(k−1)k)ξke−ξ2/2 を確かめよ。
(3) X(ξ)=ξl+1e−2ξ2 が
ε=(l+23)h−ω に対する固有関数となることを示せ。
(4) X(ξ)=(ξl+1−2l+32ξl+3)e−ξ2/2
が ε=(l+27)h−ω
に対する固有関数となることを示せ。
(5) (発展)Xnl(ξ)=(∑k=0nckξ2k)ξl+1e−ξ2/2
が ε=(l+2n+23)h−ω
に対する固有関数となるために要求される ck に対する漸化式を求めよ。
● 解答はこちら
(1) において l=0 とすれば1次元調和振動子と同じ方程式となる。
(3) で l=0 とすれば R(r)∝e−ξ2/2 となって、N=0 の解に相当する。
(3) で l=1 とすれば R(r)∝ξe−ξ2/2 となって、N=1 の解に相当する。
(3) で l=2 とした R(r)∝ξ2e−ξ2/2 と、
(4) で l=0 とした R(r)∝(ξ−32ξ3)e−ξ2/2 とは同じエネルギーを持ち(縮退しており)、N=2 の解に相当する。
(3) で l=3 とした R(r)∝ξ3e−ξ2/2 と、
(4) で l=1 とした R(r)∝(ξ2−52ξ4)e−ξ2/2 とは同じエネルギーを持ち(縮退しており)、N=3 の解に相当する。
(5) のように、一般には
Xnl(ξ)=(ξ の 2n 次多項式) ξl+1e−ξ2/2
の形で ε=(l+2n+3/2)h−ω に対する固有関数を作れるから(n=0,1,2,…)、
N=l+2n となるように l,n を選べば N と同じエネルギーを与える。
たとえば N=4 なら以下の通り15重に縮退する(図中の赤の点線に相当する)。
- n=2, l=0 の 1 個
- n=1, l=2 の 5 個
- n=0, l=4 の 9 個
グラフ形状†
距離を r0 を単位として測った動径分布関数 ∣Xnl(ξ)∣2 をプロットした。
上でも述べたように l=0 の解は1次元調和振動子の解にもなっているが、
εnl=h−ω(2n+3/2) となっていることからも分かるとおり、
(1次元調和振動子の n)=2(動径運動量の n)+1
という対応になっている。例えば上の n=0 のグラフは、一次元調和振動子の
n=1 のグラフの右半分、n=1 のグラフは一次元調和振動子の n=3 のグラフの右半分と重なる。
球対称な箱形ポテンシャルと同様に、
- n が増加すると振動回数が増える → r 方向の運動量が大きくなり振幅も増える
- l が増加すると角運動量が増える → 原点付近の存在確率密度が低下する
l,n が増加すると運動エネルギーが増えるほか、分布が原点から遠ざかるためポテンシャルエネルギーが増加する。当然、l が増えれば遠心力ポテンシャルのエネルギーも増加する。
エネルギーは (l+2n+3/2)h−ω なので、
n の増加は l の増加の2倍、エネルギーを増やす。
N=0 の解†
φnxnynz000(x,y,z)=(πr01)3/2e−(x2+y2+z2)/2r02=R00(r)4π(πr01)3/2e−2r02r2⋅Y00(r,θ,ϕ)4π1=φn,l,m0,0,0(r,θ,ϕ)
すなわち n=0 の解は 1s 状態そのものであることが分かる。s 状態の特徴として、波動関数自体が球対称になっている。
N=1 の解†
N=l+2n=1 は n=0,l=1 を表すから、球座標における解は
φ0,1,m(r,θ,ϕ)∝∝R01(r)r0re−2r02r2Y1m(θ,ϕ)
と表される。ただし m=−1,0,1 である。
一方、x,y,z 座標で解いた解は N=nx+ny+nz=1 より、(nx,ny,nz)=(100),(010),(001) のいずれかとなり、それぞれ px,py,pz 軌道に相当する。
両者の関係は、
px: φ100(x,y,z)∝r0xe−2r02r2=r0rsinθcosϕe−2r02r2∝=R01(r)r0re−2r02r2⋅{∝Y11(−sinθeiϕ)−∝Y1−1sinθe−iϕ}/2=2−1{φ0,1,1(r,θ,ϕ)−φ0,1,−1(r,θ,ϕ)}
同様に、
py: φ010(x,y,z)=2i{φ0,1,1(r,θ,ϕ)+φ0,1,−1(r,θ,ϕ)}
pz: φ001(x,y,z)=φ0,1,0(r,θ,ϕ)
の関係があることが分かる。
同じ固有値に属する固有関数の線形結合はやはり固有関数であるため、
これらの px,py,pz はすべて l=1 の固有関数、すなわち p 状態である。
一方、pz に関しては m=0 の固有関数となるものの、px,py に関しては m=1 と m=−1 の状態の線形結合となるため、l^z の固有関数とはならない。
これら2つについて lz を測定すれば、1/2 の確率で h− が、
1/2 の確率で −h− が、観測されることになる(それぞれ m=+1,−1 に対応する)。
数学的にはエネルギー固有値 ε=(1+23)h−ω の固有空間は3次元であり、
φ100(x,y,z),φ010(x,y,z),φ001(x,y,z)
と、
φ1,0,−1(r,θ,ϕ),φ1,0,0(r,θ,ϕ),φ1,0,+1(r,θ,ϕ)
とはどちらもこの空間を張る正規直交基底となっている。
両者を結ぶ上記の方程式は基底の変換に対応し、行列を使えば
(φ100φ010φ001)=(φ0,1,1φ0,1,0φ0,1,−1)−1/201/2i/20i/2010
と書ける。これは正規直交基底から正規直交基底への基底の変換行列であるからユニタリー行列となっている。
N=2 および N=3 の解†
N=2 の解のうち、
φ200, φ020, φ002 は l=0 の解と l=2 の解の線形結合で、
φ110, φ011, φ101 は l=2 の解のみで
それぞれ表せる。すなわち後者は d 状態に対応するが、前者は s 状態と d 状態の混合、
すなわち l^2 の固有関数ではなく、
全角運動量を測定すれば 1/3 の確率で l=0 を、2/3 の確率で l=2 を取る。
N=3 の解のうち、
φ111 は純粋な f 状態であるが、
他は、p 状態と f 状態の混合となっている。
詳しい計算はこちら
正規直交系同士の変換†
ここで見たように、エネルギー量子数 N=nx+ny+nz=l+2n で指定されるエネルギー固有値 εN は N/=0 に対して縮退しており、その固有空間に取った2つの正規直交完全系が
{φnxnynz} および {φn,l,m} である。
このような正規直交系同士の間の変換はユニタリー行列で表される。
ここで2つの正規直交系の変換行列 U は 線形代数IIで学んだ基底の変換行列 のように、
(φnxnynz φnx′ny′nz′ … φnx′′ny′′nz′′)=(φn,l,m φn′,l′,m′ … φn′′,l′′,m′′)U
として定義される。
例えば、N=0 の
( φnxnynz000 )=( φn,l,m0,0,0 )( 1 )
は自明だとしても、N=1 の
( φnxnynz100 φ010 φ001 )=( φ0,1,−1 φ0,1,0 φn,l,m0,1,1 ) U−1/201/2i/20i/2010
に対しては、
U†U=−1/2−i/200011/2−i/20−1/201/2i/20i/2010=100010001
を確かめられる。
2つの正規直交系 {ϕi} と {ψi} の間の変換行列 U がユニタリーになることは以下のように導ける。
U の要素を uij として、
U=( uij )
と置くと、
(ϕ1ϕ1…ϕn)=(ψ1ψ1…ψn)U
は、
ϕi=k∑ψkuki
の関係を表す。これを使って {ϕi} の正規直交性を変形すると、
(ϕi, ϕj)=k∑l∑(ψk, ψl)uki∗ulj=k∑l∑δkluki∗ulj=k∑uki∗ukj=δij
を得るが、これはそのままユニタリーの条件
U†U=E
を表している。
前の単元 <<<
量子力学Ⅰ
>>> 次の単元