« 甲斐智美女流二段がマイナビ新女王に!! | トップページ | 原子核のγ崩壊とメスバウアー効果(3) »

2010年4月22日 (木)

原子核のγ崩壊とメスバウアー効果(2)

原子核のγ崩壊とメスバウアー効果の続きです。

 

いきなり本題に入ります。

電磁波を与える真空中のマクウェル方程式(Maxwell eq.)は線型なので,電場と磁場の任意の解,は電気的波(E波=TM波;transverse magnetic wave)と磁気的波(M波=TE波;transverse electric wave)に分解できます。

すなわち,,EM,EM,(ⅰ)E波=TM波:EEr=Er,BEr=0 (or rEErE,rBE=0 )(ⅱ)M波=TE波:EMr=0,BMr=Br (or rEM=0,rBMrB)の和に分解されます。

 

(例えば2009年11/7の記事「光(電磁波)の散乱(4)」参照)

一方,前記事で述べたように電磁ポテンシャルφ(,t)=φ^()exp(-iωt),(,t)=^()exp(-iωt)の空間部分:φ^(),^()の任意の1成分ψ()はヘルムホルツ方程式(Helmholtz eq.):(∇2+k2)ψ()=0 の解です。

 

すなわち,ψ()=ψ(r,θ,φ)は[r-2(∂/∂r){r2(∂/∂r)}+k2-(2/hc2)/r2]ψ()=0 を満たします。

そこで,動径関数(radial fuction)因子をflm(kr)とすると,ψ()は独立な変数分離解の和としてψ()=Σl=0Σm=-lllm(kr)lm(Ω)と多重極(multi-pole)に展開できます。

lm(Ω)=Ylm(θ,φ)は球面調和関数(spherical harmonics)です。これは軌道角運動量の固有値方程式:2lm(Ω) =hc2l(l+1),Lzlm(Ω)=mYlm(Ω)を満たします。

そこで,φ(,t),A(,t)の空間微分,時間微分の線形和である=-∇φ-∂/∂t,=∇×も多重極展開できるはずです。

先出の2009年11/7の記事光(電磁波)の散乱(4)でも,M.Born,E.Wolf著(草川徹)「光学の原理」(東海大学出版会)などを参考にして球状散乱体によるレーリー(Raileigh)散乱,ミイ(Mie)散乱の境界条件を満たすE,E,M,Mの具体的な展開形を与えましたが,ここでは別の一般的方法で多重極展開の具体形を得たいと思います。

まず,マクスウェル方程式(Maxwell eq.):∇=ρe,0,∇×/∂t,∇×/∂t+e,および連続方程式:e=-∂ρe/∂tにおいて,,,,の時間依存性がexp(-iωt)のみである角振動数ωの単色波の形式を考えます。

ただし,以下では混乱は生じないと思われるので,例えば電場(,t)=^()exp(-iωt)の空間部分^()も,(,t)と同じく(),あるいは単にと表わすことにします。

そして,振動数がωの単色波では上記のマクスウェルの方程式系と連続方程式は=ρe,0,∇×,∇×=-iωe,およびe=-∂ρe/∂tに帰着します。

特に"ρee0 のとき=真空中"では,この方程式系は0,0,∇×,∇×=-0ε0ω=-iωc-2です。

これらの式から,真空中の,を具体的に得るには(∇2+k2)=0,および∇=0 を解いて(ic/k)∇×,=--1∇×とすればよいことがわかります。

そこで,まず後の便宜のためベクトルヘルムホルツ方程式:(∇2+k2)=0 の一般解()を,動径関数を球面ハンケル(Hankel)関数hl(ν)(x)(ν=1,2)の線形結合とする展開式:()=Σl=0Σm=-ll[lm(1)l(1)(kr)+lm(2)l(1)(kr)]Ylm(Ω)で表現します。

球面ハンケル関数hl(1)(x),hl(2)(x)というのはヘルムホルツ動径方程式のjl(x),nl(x)とは異なる選択の独立解:hl(1)(x)=jl(x)+inl(x),hl(2)(x)=jl(x)-inl(x)です。

さて,()の多重極展開の表現式を∇0 に代入するとΣl,m{lm(ν)l(ν)(kr)Ylm(Ω)}0 (ν=1,2)を得ます。

ここで,右辺の∇に公式:∇=(/r)(∂/∂r)-{i/(c2)}(×)を用います。

 一応,この公式を証明しておきます。

(証明):=(-ic)(×∇)より,(×)i/(-ic)ijkjk)/(-ic)=εijkεklmjlmilδjm-δimδjl)xjlm=xijj-xjji=xi(∇)-r2iです。

 

 つまり,(×)/(-ic)=()-r2∇なので,∇=(/r)(∂/∂r)-{i/(c2)}(×)を得ます。(証明終わり)

あるいは,別の極座標を用いた証明もあります。

(別証明):∇=r(∂/∂r)+θ-1(∂/∂θ)+φ-1sin-1θ(∂/∂φ)であり,r/r,r×r=0,r×θφ,r×φ=-θですから,=(-ic)(×∇)=(-ic){-θsin-1θ(∂/∂φ)r-1(∂/∂θ)+φ(∂/∂θ)}です。

さらに,×=(-ic){θ(∂/∂θ) +φsin-1θ(∂/∂φ)}が得られます。

そこで,この極座標表現からも公式:∇=(/r)(∂/∂r)-{i/(c2)}(×)を導くことができました。(証明終わり)

さらに,F(r)をrだけの微分可能な任意関数とすれば,{/(-ihc)}F(r)=(-i)(×∇)F(r)=(-i/r)(×)(dF/dr)=0,つまり(r)=0 です。そこではrだけの関数を素通りします。

以上から,式:∇0 or Σl,m{lm(ν)l(ν)(kr)Ylm(Ω)}0 はΣl[{dhl(ν)(kr)/dr}Σm{lm(ν)lm(Ω)}-c-1-1l(ν)(kr)×{Σmlm(ν)lm(Ω)}]=0 を意味します。

もしも,()=Σl=0Σm=-ll[lm(1)l(1)(kr)+lm(2)l(1)(kr)]Ylm(Ω)がTM波(E波):EであればrBE=0より,Σm=-ll[lm(ν)l(ν)(kr)Ylm(Ω)=0 ですから,[×{Σmlm(ν)lm(Ω)}]=0 を得ます。

Σm=-ll[lm(ν)l(ν)(kr)Ylm(Ω)=0,かつ[×{Σmlm(ν)lm(Ω)}]=0 が常に成立するためにはElm(ν)を定係数として,Σmlm(ν)lm(Ω)=c-1ΣmElm(ν)lm(Ω)と書ければ十分です。

この表現を採用してE()=Σl=0Σm=-llΣνlm(ν)l(ν)(kr)Ylm(Ω)に代入すると,E()=Σl=0Σm=-llΣνElm(ν)l(ν)(kr)c-1lm(Ω)を得ます。

そして,(ic/k)(∇×)より,E()=(ic/k)Σl=0Σm=-llΣνElm(ν)l(ν)(kr){∇×c-1lm(Ω)}です。

同様にrEM0 よりM()=Σl=0Σm=-llΣνMlm(ν)l(ν)(kr)c-1lm(Ω)です。

 

そこで,=--1∇×からM()=-i(ck)-1∇×M()={i/(ck)}Σl=0Σm=-llΣνΣm=-llMlm(ν)l(ν)(kr){∇×c-1lm(Ω)}を得ます。

ここで,∫(lm(Ω)*l'm'(Ω))dΩ=δll'δmm'と直交規格化されたベクトル球面調和関数:lm(Ω)をlm(Ω)≡hc-1lm(Ω)/{l(l+1)}1/2(i){×∇lm(Ω)}/{l(l+1)}1/2で定義して導入します。

また,fE,Mlm(kr)≡E,Mlm(1)l(1)(kr)+CE,Mlm(2)l(2)(kr),E,Mlm(ν)≡{l(l+1)}1/2E,Mlm(ν)と置きます。

すると場は()=E()+M()=Σl=1Σm=-ll[AElmElm(kr)lm(Ω)-AMlm{i/(ck)}{∇×{fMlm()lm(Ω)}},()=E()+M()=Σl=1Σm=-ll[AElm(i/k){∇×{fElm()lm(Ω)}+AMlmMlm(kr)lm(Ω)}と表わすことができます。

これで,真空中の磁場()と電場()の多重極展開の具体的な形が得られました。

 

右辺の級数和:Σl=0をl=1から始まるΣl=1に書き換えたのはl=0 (球対称なs波)なら,(r)=0 のためlm(Ω)=hc-1lm(Ω)}/{l(l+1)}1/2が存在しないからです。

次に,原子核の中心=0 付近の限られた領域を考えると,ここは真空ではなくρe(,t)=ρe()exp(-iωt),e(,t)=e()exp(-iωt)なる振動電荷,振動電流が存在します。

 

そして,物質の電磁気学の現象論からεを核の誘電率,μを核の透磁率として=ε=ε0,μ-1μ0-1と書けます。

ただし,実際には非負電荷のみから成る原子核では電気分極:はゼロですからε=ε0,=ε0です。核磁気モーメントの方は存在して(,t)=()exp(-iωt)と書けます。

そこで,先に与えた単色波のマクスウェル方程式:=ρe,0,∇×,∇×=-iωe,連続方程式:eρeに,=ε0,μ0-1を代入してρee0 の真空のときと同じく,を消去してみます。

=ρe/ε0,0,∇×,∇×(μ0)=-iε0μ0ω+μ0e=-ic-2ω+μ0e,eiωρeですね。

そして∇=ρe/ε0eiωρeから,∇{i/(ε0ω)e}=0 となりρeを消去できます。

 

つまり,'を'≡i/(ε0ω)eで定義すれば∇'=0 です。さらに,∇×(μ0)=-ic-2ω+μ0e'を用いると∇×(μ0)=-ic-2ω'と書き直せます。

そして,発散(divergence)がゼロの,'については,∇×(∇×)=∇(∇)-∇2=-∇2,∇×(∇×)=∇(∇')-∇2'=-∇2'が成り立ちます。

 故に,∇×(μ0)=-ic-2ω'より∇×{∇×(μ0)}=-ic-2ω∇×'=ic-2ω∇×+μ0∇×eですが,∇×なので(2+k2)=-μ0{∇×e+∇×(∇×)}を得ます。

また,∇×より∇×{'-i/(ε0ω)e}ですから∇×{∇×{'-i/(ε0ω)e}}=iω∇×です。

 

故に,∇×(μ0)=-ic-2ω'+μ0∇×を用いて-∇2'-i/(ε0ω)∇×(∇×e)=k2'+iμ0ω∇×となります。

 以上から,(2+k2)'=-iμ0-1{∇×(∇×e)+2∇×}が得られます。

並べて書くと,(2+k2)=-μ0{∇×e+∇×(∇×)},(2+k2)'=-iμ0-1{k2∇×∇×(∇×e)}です。

 

式が過剰となりますから整合性が必要ですが,これらと∇0,∇'0,∇×=-ic-2ω'μ0∇×,および∇×'=iω+i/(ε0ω)∇×eを組み合わせます。

真空のときと同様,EM,''E'Mと分割して,E0,rBE0,および'M0,'M0 を同時に満たすE,'MとしてE()=Σl=1Σm=-llElm(kr)lm(Ω),'M()=Σl=1Σm=-llMlm(kr)lm(Ω)なる展開形が想定されます。

前と同じく,M=-iω-1∇×'MとしてみるとEMiω-1∇×'+k-2μ0∇×eから,E=-iω-1∇×'E-2μ0∇×eとなることが必要です。

また,'E(ic/k)∇×Eとしてみると''E'Mick-1∇×-icμ0-1∇×から,'Mick-1∇×Micμ0-1∇×となることが必要です。

これで辻褄が合うなら,E()=Σl=1Σm=-llElm(kr)lm(Ω),および'M()=Σl=1Σm=-llMlm(kr)lm(Ω)から,動径関数FE,Mlm(kr)だけが真空中のfE,Mlm(kr)と異なるという形で原子核中心領域の磁場,電場の空間部分が次のように表わされます。

すなわち,()=E()+M()=Σl=1Σm=-ll[FElm(kr)lm(Ω)-{i/(ck)}{∇×{FMlm(kr)lm(Ω)}},()=E()+M()=Σl=1Σm=-ll[(i/k){∇×{FElm(kr)lm(Ω)}+FMlm(kr)lm(Ω)}です。

問題はこれで全ての辻褄が合うかどうか?です。

 

まず,E=-iω-1∇×'E-2μ0∇×e'E(ic/k)∇×Eから,E=-k-22E-2μ0∇×eなので,(2+k2)E=μ0∇×eです。

 

同様に,'E=-k-22'Eic-3μ0∇×(∇×e)により,(2+k2)'E=iμ0ck-1∇×(∇×e)です。

また,'Mick-1∇×Micμ0-1∇×,M=-iω-1∇×'Mから'M=-k-22'Micμ0-1∇×により,(2+k2)'Miμ0ck∇×です。

 

同様に,M=-k-22M-2μ0∇×(∇×)により,(2+k2)M=μ0∇×(∇×)が得られます。

これらを通して,()=E()+M(),()='E()+'M()は,確かにそれぞれ,方程式:(2+k2)=-μ0{∇×e+∇×(∇×)},(2+k2)'=-iμ0ck-1{∇×(∇×e)+2∇×}を満たしています

そして,与えた多重極の展開形式は,確かに∇E,M'E,M0,およびE'M0 を満たします。

E=-iω-1∇×'E-2μ0∇×e,'Mick-1∇×Micμ0-1∇×,E=0,rE'M=0 によ{∇×('E+iμ0-1e)}=0,{∇×(Nμ0)}=0 が要求されます。

 

しかし,これはEM,'E'Mの上記分割が妥当でその他全ての式を満たす解であれば,当然,M,'Eが満たすべき条件です。言うなればトートロジーですね。

 

さて,今得たのは電荷密度,電流密度,磁気モーメントが全て存在する内部領域における解ですがこれらはその境界で領域の外の解に滑らかに接続するはずです。

 

そして,この境界条件はFElm(kr)→AElmElm(1)l(1)(kr),FMlm(kr)→AMlmMlm(1)l(1)(kr)で与えられます。

ここで,r→ ∞での球ハンケル関数の挙動はhl(1)(kr)→(-i)l+1exp(ikr)/(kr),hl(2)(kr)→il+1exp(-ikr)/(kr)ですから,原子核中心付近から放射される波には外向き成分しかないとして領域の外の真空中での動径関数fE,Nlm(kr)≡E,Nlm(1)l(1)(kr)+CE,Nlm(2)l(2)(kr)でhl(2)(kr)の係数CE,Nlm(2)を全てゼロと置きました。

動径関数:FElm(kr)を具体的に解くために,E()=Σl=1Σm=-llElm(kr)lm(Ω)を(2+k2)=-μ0∇×eに代入します。

  

すると,[r-2(∂/∂r){r2(∂/∂r)}+k2-(2/hc2)/r2]E()=Σl=1Σm=-ll{d2/dr2+(2/r)d/dr+k2-l(l+1)/r2}FElm(kr)lm(Ω)=-μ0∇×eです。

 この方程式の両辺の左からl'm'(Ω)*を掛けてdΩ積分し規格化直交条件を用いると,{d2/dr2+(2/r)d/dr+k2-l(l+1)/r2}FElm(kr)=-μ0lm(Ω)* {∇×e()}dΩ≡-KE(r)です。

 同様に,'M()=Σl=1Σm=-llMlm(kr)lm(Ω)を(2+k2)E'M=-iμ0ck∇×に代入します。

 

 すると,[r-2(∂/∂r){r2(∂/∂r)}+k2-(2/hc2)/r2]E'M()=Σl=1Σm=-ll{d2/dr2+(2/r)d/dr+k2-l(l+1)/r2}FMlm(kr)lm(Ω)=-iμ0ck∇×です。

これも左からl'm'(Ω)*を掛けてdΩ積分し規格化直交条件を用いると,{d2/dr2+(2/r)d/dr+k2-l(l+1)/r2}FMlm(kr)=-iμ0cklm(Ω)*{∇×()}dΩ≡-KM(r)を得ます。

今日はここで終わります。

いやあ,今回も細かい計算を始めると寝食も忘れてしまいます。こういう習性だけは,体力も精神力も落ちて老眼鏡が不可欠な今も昔と変わりませんね。

 

この関係の記事を後回しにしていたのは自分の計算に納得できず,いくら検算しても合わないような問題点があったからです。

 

そしてまだ終わっていません。テーマは古いのですがね。

  

ただ,別に締め切りも無くて暇があるだけが救いです。

 

ま,いくら心血を注いでも一銭にもなりませんがね。あ,でも自己満足を得られるからそれで十分かぁ。。。

ゼロからでなく色々と参考書もあるのですが,一度つまずくと細かい式については専門の本も100%は信用できませんから大変かな。。

八木浩輔著「原子核物理学」(朝倉書店),八木浩輔 著「原子核と放射」(朝倉書店),砂川重信 著「理論電磁気学(第2版)」(紀伊国屋書店),ジャクソン著(西田 稔 訳)「電磁気学」 

ブックオフオンライン 

iconオンライン書店 boople.com(ブープル) 

|

« 甲斐智美女流二段がマイナビ新女王に!! | トップページ | 原子核のγ崩壊とメスバウアー効果(3) »

111. 量子論」カテゴリの記事

113. 原子核物理」カテゴリの記事

コメント

この記事へのコメントは終了しました。

トラックバック


この記事へのトラックバック一覧です: 原子核のγ崩壊とメスバウアー効果(2):

« 甲斐智美女流二段がマイナビ新女王に!! | トップページ | 原子核のγ崩壊とメスバウアー効果(3) »