水素様原子の波動関数
冬場と夏場は若干心臓の動きが悪いと言われているし,自身そんな気もするので,自分の体調を考えてブログの1つの記事の量も半分程度にしようと思います。
別の場所で"せいたかのっぽ"さんのコメントも頂きましたが,私は皮肉ではなく物性を定性的に捉えて説明できるような化学屋さんを尊敬しています。
私などは本質的に他人より馬鹿な(認識の閾値が低い?)ので,本来は定性的なものでさえ,なんとか定量的に捉えて数式の助けを借りなければほとんど理解できず,物理的意味を問われても数式を使用しないと答えられない人間です。
数式にでも頼ってなんとかうまく咀嚼できた後に,他人に説明できるほどの物理的イメージと語彙が整って初めて定性的説明が可能だと思うので,定性的に説明できたり,定性的に理解するのは大変な能力だと思って若干尊敬してしまいます。
さて,前記事の水,氷,水蒸気と関わる話の続きとして化学結合に関して書こうと思いましたが,まず分子以前の原子の理論まで戻り,そこから始めます。
e>0 を素電荷とするとき,中心にZe(Zは正の整数)の電荷を
持つ質量がm1の原子核があり,周りに質量がm2,電荷が(-e)の
電子が1個だけあるような系の1個の原子を想定します。
こうした原子を水素様原子と呼びますが,Z=1なら,これは水素
原子に他なりません。
さて,原子核の位置ベクトルをr1,運動量をp1,電子の位置ベクトル
をr2,その運動量をp2とすると,系の全Hamiltonian:Htotは,
Htot=p12/(2m1)+p22/(2m2)+V(|r1-r2|)となります。
ここで,V(r)=-Ze2/(4πε0r)は静電Coulombポテンシャル
による位置エネルギーです。
系の重心の位置をR≡(m1r1+m2r2)/M (M≡m1+m2)とし,
全運動量を,P≡p1+p2とします。
そして,電子の原子核に対する相対運動を重心運動から分離して
考察するために,原子核を原点と想定したときの位置ベクトルを
r≡r2-r1で定義します。
さらに,電子の換算質量mを1/m≡1/m1+1/m2で定義します。
(1/2)m(dr/dt)2=p2/(2m)が成立するように,
相対運動量:pを,
p≡m(dr/dt)=(1/m1+1/m2)-1{d(r2-r1)/dt}
=(1/m1+1/m2)-1(p2/m2-p1/m1)=(m1p2-m2p1)/M
で定義すれば,
p=p2-m2P/Mです。
そして,このとき,Htot=P2/(2M)+p2/(2m)+V(r)
=HG+H と書けます。
ここに,HG≡P2/(2M)は系の重心運動のHamiltonianです。
他方,H≡p2/(2m)+V(r)は,原子核を原点とした電子の相対
運動のHamiltoniaと考えられます。
これに対するHamiltonの正準方程式は,
dR/dt=∂Htot/∂P=∂HG/∂P=P/M=(p1+p2)/M,
dP/dt=-∂Htot/∂R=-∂HG/∂R=0,
および,
dr/dt=∂Htot/∂p=∂H/∂p=p/m=p1/m1-p2/m2,
dp/dt=-∂Htot/∂r=-∂H/∂r=-∂V(r)/∂r
です。
これらを見ると,系の重心運動と電子の相対運動は完全に分離され
ることがわかります。
特に,電子の質量は,m2=me ~ 0.51MeV,核子1個の質量は,
mN ~938 ~939Mev~1840meですから,
m1=ZmN~1840Zme です。
これは,先の水素様原子の原子核と電子の質量の定義では,
m1>>m2を意味するので,
この場合は,近似的に,
全質量:M=m1+m2 ~ m1(原子核の質量),
重心の位置:R≡(m1r1+m2r2)/M~ r1(原子核の位置),
換算質量:m≡m1m2/(m1+m2) ~ m2=me (電子質量),
相対運動量:p=p2-m2P/(m1+m2)~p2(電子の運動量)
が成立しています。
結局,原子核を不動の原点とした電子の1体問題を考察しているの
と,ほぼ同等です。
こうした一連の手続きを見ると,以前2007年10/26の記事:
「ケプラー問題(Kepler問題)」で古典力学の2体問題に関して,
次のような相対運動の定式化を与えたことを思い出します。
※"太陽,および惑星の質量をそれぞれm1,およびm2とし,位置ベク
トルをそれぞれr1,およびr2 とします。
そして,それらの間に働く力をF12(1が2に及ぼす力),および,
F21(2が1に及ぼす力)とします。
すると,作用反作用の法則(Newtonの第3法則)によって,
F21=-F12, かつ, (r1-r2)×F12=0
が成立します。
太陽と惑星の描く軌道のorderと比較して,太陽,惑星を質点で近似
すれば,この問題の運動方程式は,
m1(d2r1/dt2)=F21,および,
m2(d2r2/dt2)=F12
となります。
「物体2=惑星」の「物体1=太陽」に対する相対運動を考察する
ために,太陽を原点と想定したときの位置ベクトルをr≡r2-r1
によって定義します。
このとき,系の運動方程式は,唯1つの式:
d2r/dt2=d2(r2-r1)/dt2=F12/m1-F21/m2
となります。
すなわち,運動は方程式:d2r/dt2=(1/m1+1/m2)F12
だけで記述されます。
惑星の換算質量:mを1/m≡1/m1+1/m2で定義し,太陽1が
惑星2に及ぼす力Fを,改めてF≡F12で定義すると,
m(d2r/dt2)=Fとなり,結局,2体問題は1体問題に
帰着します。
そして,当面の問題では,太陽の質量m1が惑星の質量m2より,
はるかに大きい(m1>>m2)ため,1/m=1/m1+1/m2~1/m2
ですから,事実上m=m2として換算質量を惑星質量と同一視
してよいことになります。
さらに,後の便宜上,太陽の質量m1をMで表わすことにします。
さて,力Fが特に相対位置rだけに依存する力場,しかも保存力場:
つまり,渦なしの場(rotF=∇×F=0 )であるとすると,
謂わゆるPoincare'の補題により,あるポテンシャルV(r)が存在
して,F=-gradV(r)=-∇V(r)と書くことができます。
特に,Vがrの絶対値r≡|r|,または太陽と惑星の間の距離:
|r2-r1|だけの関数V(r)=V(r)=V(|r2-r1|)なら,
F=F12=-∇V=(-dV/dr)(r/r)
=(-dV/dr){(r2-r1)/r} です。
保存力場の内力は,F12=-∇2V,F21=-∇1Vで与えられるので,
内力のポテンシャルが2質点間の距離だけの関数で与えられる場合
には作用・反作用の法則は自動的に満たされます。"※
(再掲記事終わり)
しかし,太陽と惑星の場合と似ているようでも,水素様原子の場合は
古典力学で考えるよりも,量子力学の問題として捉えるのが適切で
あると思われます。
そこで,運動量をP=-ihc∇R,p=-ihc∇ なる座標表示の
演算子で置き換えることにより,問題を量子化ます。
(hc≡h/(2π)でhはPlanck定数)
すると,系の全Hamiltonianは,
Htot=-hc2∇R2/(2M)-hc2∇2/(2m)+V(r)
=HG+H ;
HG≡{-hc2/(2M)}∇R 2,H≡{-hc2/(2m)}∇ 2+V(r)
で与えられます。
そして系の波動関数をΨ(R,r,t)とすると,系の状態を支配
するSchrödingerの波動方程式は.
ihc{∂Ψ(R,r,t)/∂t}=HtotΨ(R,r,t)
となります。
特に,Ψ(R,r,t)が原子として閉じた1つの自由粒子の系:
つまり水素様原子として重心運動エネルギーが保存され,一定値EG
を取るHGの固有状態を表わす状態関数であるとするなら,
HGΨ(R,r,t)=EGΨ(R,r,t)です。
これは,[{-hc2/(2M)}∇R2]Ψ(R,r,t)=EGΨ(R,r,t)
と書けます。
そして,系全体:原子自体の波動関数は自由粒子としての境界条件
を満たしていて,運動方向も一定のはずですから,対応する波動ベ
クトルをhc2K2/(2M)=EGを満たすKとすれば,
Ψ(R,r,t)=Φ(r,t)exp(iKR-iEGt/hc)
と書けます。
このとき,波動方程式は,結局,
ihc{∂Φ(r,t)/∂t}=HΦ(r,t)
に帰着します。
こうして,今問題としている原子内電子の2体問題は,Φ(r,t)
を原子核に束縛された電子の波動関数と考えれば,通常の1粒子の
量子力学の問題に帰着します。
特に,Φ(r,t)が電子のエネルギーEを持つ定常状態:
つまり固有値Eに属するHの固有状態なら,
HΦ(r,t)=EΦ(r,t)なので,
Φ(r,t)=ψ(r)exp(-iEt/hc) と書けます。
そして,Φ(r,t)の時間に依存しない部分ψ(r)は
定常状態のSchrödinger方程式:
[{-hc2/(2m)}∇2+V(r)]ψ(r)=Eψ(r),or
[{-hc2/(2m)}∇2-Ze2/(4πε0r)]ψ(r)=Eψ(r)
を満たします。
右辺のLaplace演算子(Laplacian):∇2を極座標表示すると,
∇2=∂2/∂r2+(2/r)(∂/∂r)
+(1/r2)[(1/sinθ)(∂/∂θ){sinθ(∂/∂θ)}
+(1/sin2θ)(∂2/∂φ2)]です。
ψ(r)を変数分離してψ(r)≡R(r)Y(θ,φ)とし,
α=(-2mE/hc2)1/2>0,λ≡mZe2/(4πε0hc2α)
とおけば,
d2R/dr2+(2/r)(dR/dr)-βR/r2
+(-α2+2λα/r)R=0,および,
(1/sinθ)(∂/∂θ){sinθ(∂Y/∂θ)}
+(1/sin2θ)(∂2Y/∂φ2)=-βY
となります。
Y(θ,φ)をさらに変数分離しY(θ,φ)≡Θ(θ)Φ(φ)とおくと,
(1/sinθ)(d/dθ){sinθ(dΘ/dθ)}-(γ/sin2θ)Θ
=-βΘ,
d2Φ/dφ2=-γΦ
となります。
波動関数の1価性の要請から,関数Φ(φ)はφについて周期
が2πの周期関数であることが必要なので,解:Φ(φ)は,
Φm(φ)=(1/2π)1/2exp(imφ),Φ-m(φ)=(1/2π)1/2exp(-imφ)
(mは整数)の1次結合で与えられます。
それ故,適当な整数mについてγ=m2となって,
(1/sinθ)(d/dθ){sinθ(dΘ/dθ)}-(m2/sin2θ)Θ
=-βΘ
と書けます。
ここで,z=cosθと変数変換すれば,
(d/dz){(1-z2)(dΘ/dz)}+{β-m2/(1-z2)}Θ
= 0 , または,
d2Θ/dz2-{2z/(1-z2)}(dΘ/dz)+{β/(1-z2)
-m2/(1-z2)2}Θ= 0 です。
これは,z=±1,∞ に確定特異点を持つFuchs型の2階線形
常微分方程式です。
そして,この場合にFrobeniusの方法に基づいて決定方程式を解いて
得られる形式的べき級数解が収束するのは,ベキ級数が途中で有限
項で途切れて多項式の場合に限られることがわかります。
(↑ ※2007年4/28の記事「2階線形常微分方程式と確定特異点」)
つまり,β=l(l+1) (l=|m|,|m|+1,..)なる形である必要
があるため,
上の微分方程式は,
d2Θ/dz2-{2z/(1-z2)}(dΘ/dz)
+{l(l+1)/(1-z2)-m2/(1-z2)2}Θ= 0
(l=0,1,2,..,;m=-l,-l+1,..,0,1,..,l)
となります。
この方程式は,Legendreの微分方程式と呼ばれ,(l,m)に対応
する解は,Legendreの陪多項式:Plm(z)であることがわかって
います。
したがって,cを任意定数としてΘ(z)=cPlm(z)と書けます。
特に,m=0 のときのLegendreの微分方程式:
d2Θ/dz2-{2z/(1-z2)}(dΘ/dz)
+{l(l+1)/(1-z2)}Θ=0 の解に対応する多項式を,
Legendreの多項式と呼び,Pl(z)と書きます。
これは,Pl(z)≡{1/(2ll!)}{dl(z2-1)l/dzl}
で定義されます。
これから,m≠0 の陪多項式:Plm(z)は,
Plm(z)≡(1-z2)|m|/2dmPl(z)/dzm
で定義されます。
結局,(l,m)に対応する規格化された波動関数は,
Θ(θ)=Θlm(θ)
=[(2l+1)(l-|m|)!/{2(l+|m|)!}]1/2Plm(cosθ)
となり,
Y(θ,φ)=Ylm(θ,φ)=Θlm(θ)Φm(φ)
=[(2l+1)(l-|m|)!/{4π(l+|m|)!}]1/2
Plm(cosθ)exp(imφ) と書けます。
このYlm(θ,φ)は球関数,or 球面調和関数と呼ばれています。
実は,軌道角運動量をL≡r×p=(-ihc)(r×∇)と書くと,
(-hc2)[(1/sinθ)(∂/∂θ){sinθ(∂/∂θ)}
+(1/sin2θ)(∂2/∂φ2)]=L2 となるので,
∇2=∂2/∂r2+(2/r)(∂/∂r)-L2/(hc2r2)
とも表現できます。
それ故,L2Ylm(θ,φ)=l(l+1)hc2Ylm(θ,φ) です。
また,Lz=ihc(∂/∂φ)なので,
LzYlm(θ,φ)=mhcYlm(θ,φ)
とも書けます。
以上から,残る因子:動径関数の満たす動径方程式は,
d2R/dr2+(2/r)(dR/dr)-l(l+1)R/r2
+(-α2+2λα/r)R=0 となります。
ここで,改めてこの1次元の動径方程式の物理的意味を見るため,
両辺に{-hc2/(2m)}を掛けて単位を取り戻すと,
[{-hc2/(2m)}(d2/dr2)-{2hc2/(2mr)}(d/dr)
+l(l+1)hc2/(2mr2)-2λαhc2/(2mr)]R
=-{hc2/(2m)}R です。
左辺には,l≠0 のときのみ:軌道角運動量がゼロでないときのみ
出現する斥力項:l(l+1)hc2/(2mr2)=L2/(2mr2)
がありますが,
これは非線形な曲線座標への座標変換を行ったため,不可避的に
生じる仮想的な力の項です。
一見,回転運動の慣性力の1つである遠心力の表現に一致している
ので遠心力ポテンシャルと呼ばれています。
かつて,旧@nifty物理フォーラムでの"座標変換と慣性力"という
質問に関連して,2006年11/16の記事:
「非慣性座標系で現われる慣性力」でも触れましたが,
遠心力という呼称については,当時,神戸大学名誉教授の松田卓也先生
から,「こうした呼び方はふさわしくないのでは?」という内容のご
指摘を受けました。
その後,もっと前の2006年6/30の記事「慣性力の反作用」の内容とも
関連して,
松田先生ご自身が,気象庁の木下篤哉(AXION)氏
([相対論の正しい間違え方」共著者)や琉球大学の
前野(いろもの物理学者)昌弘氏らとの協議を通じ,
何らかの結論(例えば慣性力には反作用を考える必要はない。)
を得られたという連絡をもらったことを,ちょっと懐かしく
思い出しました。
余談はさておき,実際に動径方程式を解く必要があるので,
次にこれを再掲します。
まず,再び,動径方程式:d2R/dr2+(2/r)(dR/dr)
-l(l+1)R/r2+(-α2+2λα/r)R= 0 を明示します。
ここで,ρ≡2αrとおけば,これは
d2R/dρ2+(2/ρ)(dR/dρ)+
{-1/4-l(l+1)R/ρ2+λ/ρ}R= 0 となります。
そして,これの解:R(r)=Rl(r)はρ→∞で近似方程式:
d2R/dρ2-R/4= 0 を満たす必要があるため,
ρ→∞で有限なある関数をTl(ρ)として,
Rl(r)=exp(-ρ/2)Tl(ρ)と書けるはずです。
そして,Tl(ρ)に対する線形常微分方程式におけるρのベキ級数解
が有限に留まるためには,やはり有限項で途切れて多項式となるこ
とが必要です。
結局,ρ-lTl(ρ)がLaguerre多項式(ラゲール多項式)になることが
わかります。
λ=n(n=l+1,l+2,..)であることが必要で,個々の正整数n
に対し規格化を除いて,Tl(ρ)=ρlLn+l2l+1(ρ)となります。
以上から,Hのエネルギー固有値:E=-hc2α2/(2m)に属する
変数分離した固有関数は,
ψnlm(r,θ,φ)=Nnlexp(-αr)(αr)lLn+l2l+1(2αr)
Ylm(θ,φ)
(l=0,1,2,..,n-1;m=-l,-l+1,..,0,1,..,l)
で与えられることがわかります。
n=λ=mZe2/(4πε0hc2α)なので,
α=mZe2/(4πε0hc2n)ですから,理想的に電子が1個
だけの場合にはエネルギーは主量子数nだけで決まって,
E=En=-mZ2e4/{(4πε0)2(2hc2n2)}
となり,状態は縮退しています。
こうした縮退は,電子がもう1個増えたり,電場や磁場がかかるだけ
で簡単に解けます。
Σm=-ll=2l+1であり,また,Σl=0n-1(2l+1)=n2ですが,実際には
電子には,この他に上向きか下向きかというspinの自由度が2つあり
ますから,エネルギーの縮退度は2n2です。
一般,にエネルギーEnを持つ波動関数は,
ψ(r,s)
=Σl=0n-1Σm=-llΣszcnlmszψnlm(r,θ,φ)σsz(s)
となります。
ただし,σsz(s)はスピン波動関数で,これは,
sz=1/2,または-1/2に応じて,t[1.0],またはt[0,1]
と,2成分spinor形で表現されることがあります。
ここで,上添字 t は転置行列を示す添字で,今の場合は行ベクトル
を列ベクトルにすることを意味します。
そして,この縮退度の数列:{2n2}n=1,2,3..={2,8,18,..}が,
主量子数n=1,2,3,..etc.の謂わゆるK殻,L殻,M殻..etc.
の電子軌道に入ることが可能な"最大電子数=縮退した状態数"
を示しています。
原子などの系で,"状態の数と可能な占有電子の数が常に等しい。"
という性質があるのは,
電子のように,spinが半奇数であるFermi粒子(Fermion)
に対しては,
"同種のFermi粒子の系では,2つ以上の粒子が同時に同じ状態を
占めることはできない。"というPauliの禁制原理,あるいは,
Pauliの排他原理(exclusion principle)があるからです。
最後に,同位体(isotope)の話を少しします。
原子核の質量は電子質量よりはるかに大きいので,電子の換算質量
mは電子質量meにほぼ一致することは既に述べました。
例えば,通常の水素と,重陽子(deuteron)を原子核とする重水素
(deuterium)の同位体では,原子番号は同じZ=1ですから両者の
束縛状態の電子のエネルギー準位:
En=-me4/{(4πε0)2(2hc2n2)}は,
全く同一であるように見えます。
しかし,mは電子の換算質量でありm~meですがm≠meです。
ほんの僅かな違いですが,原子核の質量が異なるためmが微妙に
異なることになり,厳密に言えば同じではありません。
私は実験にはかなり疎い方なので,実際にそうした実験がされて
いるのか?,あるいは,そうした判別法は実用的ではなくてやられ
てないのか?などについて存じませんが,
上述の事実から分光学的に観測される光のスペクトルの微妙な差異を
見る方法でも,原理的には同位体を判別できると思います。
とりあえず今日はこのくらいにしておきます。
本記事を書くに当たり,量子力学のテキストをいくつか参考に
しましたが, 標準的な話を参照しただけで,これといったテキ
ストに依存するような特別な情報はなかったので,
本記事については参考文献の詳細は書かないことにします。
緊急告知!! Attention!! [広告宣伝]です。
年末、クリスマスに引き続き2008年1月20日まで期間限定で下記URLの「TRS健康ランド」でご注文の方にはSCSを除く全商品を3割引きにて提供致します。
そこのお酒のみの方,年末年始で飲食の機会の増えたあなた、体によいし特に肝臓によいウコンがいいですよ!! そして特に今回提供する沖縄原産の純粋な黒ウコンは当店が専売の新製品ですが古くから沖縄地方ではいわゆる男性機能に効果があると言われています。
おやおや、そこの静電気バチバチの人、いいものありますよ。。。
それから農薬を落とした後の皮がピカピカに光っているリンゴなど商品として販売する際の見栄えをよくするなどのために化学処理をした食品を安全に洗浄する新商品の洗浄液SCSはいかがですか。。。
http://www.rakuten.co.jp/trs-kenko-land/「TRS健康ランド」-- 黒ウコン,SCS(洗浄剤)専売などの店: 私が店長 です。
http://www.mediator.co.jp/category/pages.php?id=115「中古パソコン!メディエーター巣鴨店」
http://folomy.jp/heart/「folomy 物理フォーラム」サブマネージャーです。
人気blogランキングへ ← クリックして投票してください。(1クリック=1投票です。1人1日1投票しかできません。)
| 固定リンク
「111. 量子論」カテゴリの記事
- クライン・ゴルドン方程式(8)(2016.09.01)
- クライン・ゴルドン方程式(7)(2016.08.23)
- Dirac方程式の非相対論極限近似(2)(2016.08.14)
- Dirac方程式の非相対論極限近似(1)(2016.08.10)
- クライン・ゴルドン方程式(6)(2016.07.27)
「112. 原子・分子物理」カテゴリの記事
- 水素様原子の微細構造(補遺3-1)(2011.10.04)
- 分子と点群(3)(2009.03.01)
- 分子と点群(2)(2009.02.26)
- 分子と点群(1)(2009.02.21)
- 磁場の中の原子(ゼーマン効果)(2)(2008.04.09)
この記事へのコメントは終了しました。
コメント
ども良子さん。コメントありがとうございます。TOSHIです。
自然科学の話以外には,本人にとっては自慢話ですが,どうしようもないモットーもモラルもない変態オヤジのブログ記事をお読みいただき,その上コメントしてもらえるだけでもありがたいことです。。
よかったら,また寄ってやってください。
TOSHI
投稿: TOSHI | 2008年1月12日 (土) 19時35分
初めまして、良子です。難しいブッログは全然解りませんが、他のは楽しく読みました。多趣味ですね、女好きで、酒飲みですね、これじゃ女房が来ないよね。カレーと卵が大好きな様ですが、安上がりで良いですねー。toshitt様体を大事にして下ださいね。それでは又ね。
投稿: 良子 | 2008年1月12日 (土) 03時33分