光の量子論13

※「光の量子論12」からの続きです。

(※余談):昨日(12/7(土))は,1学年年下ですが,全員69歳の障害者

友達２人に誘われて昼から都営新宿線住吉駅から健康な足で

徒歩４分(私は徒歩10分くらい)の「ティアラ江東」で開かれた江東区

の「第39回障害者福祉大会」というイベントに行ってきました。

(私,1985年から1994年まで塩浜や豊洲付近の運河のそばの

分譲マンション10階に９年間住んでいて,巣鴨に住む前は江東区民

だったので江東区に全く無縁というワケではありませんが。。)

最近は,近くの買い物と月１の病院以外は,ほぼ寝た切りに近い

ので現地まで行って３時間程度のイベントを見るだけなのですが,

身体がたえられるか少し心配でしたが,まあ,身体障害要介護２でも

なぜか要支援認定らしい「林家こん平」師匠よりは,ましな程度に

動けると思うので鑑賞だけなら大丈夫だろうということで,一応,

アフリカの太鼓とか歌や踊りとかを見て楽しめました。

ただ,39回も福祉をやっている割には会場も,まわりも階段が多く,

来年はパラリンピックもある,といいながら1995年のサリン事件

以来,警備,セキュリティに重きを置きすぎてるのか,来年に備えて

工事中段階とはいえ都内の駅には,ゴミ箱もベンチも少なく,足が

不自由な身にはバリアが多いのが気になりました。

「向こう３軒両隣り」とか「渡る世間に鬼はない」とかよりも

「ヒトを見たらドロボーと思え」「敷居またげば適ばかり」の

都会的思考の方が上にきて,99匹の健康な子羊を守るため,

迷える１匹の子羊などは,ほったらかしの,弱者にやさしくない

街つくりに向かっているようです。

オリンピックというと,多くの外国客への見栄で,ホームレス隠し

や,一時的に臭いモノにはフタをする,という,北京五輪や上海

万博時の中国や前の東京五輪時のわが国でもあった恥ずかしい

行為が,今回の東京オリ・パラでもあるカモしれない,と心配して

います。

最近の私は,自分は,スグにいなくなる世界なので関係ないクセに,

環境回復して,地球というか生命体の滅亡を遅らせるに役立つだろう

「人工光合成」や「プラスチックの分解,無害化」のバクテリアに興味

を持って調べています。(余談終わり※)

さて,本題です。

前回は第2章 原子・放射相互作用の量子力学の§2.10

(放射減衰を伴なうRabi振動)の項の残りの部分を記述

して終わりました。

今回はその続き次の節からです。

• 2.11衝突広がり

さて,§2.8,§2.9で論じた放射広がりと飽和広がりは,入射光

ビームで励起された原子遷移にとって避けられない性質でした。

それ以外の線幅を広げる機構の重要性は,原子気体の物理的

条件に左右されます。

このうち,原子の運動に起因する主要な効果は「衝突広がり」と,

「ドプラー(Doppler)広がり」の機構です。これらを本節と次節で

論じます。

まず,ここでは「衝突広がり」の機構を説明しますが.完全という

わけではなく,どうにか間に合う程度の詳しさで,そのプロセスを

考察します。原子気体における原子間衝突の発生は.一種の

確率過程(stochastic process)です。

「気体運動論」によると,原子が時刻τから,時刻τ＋ｄτまでの

間に自由飛行を継続する確率はｄτに比例し,それは,ｐ(τ)ｄτ

＝(1/τ₀)exp(－τ/τ₀)ｄτ.(2.131)で与えられると考えられます。

ただし,τ₀は衝突なしで飛行する平均自由飛行時間で,これは

1/τ₀＝(4ｄ²Ｎ/Ｖ){π/(βＭ)}^1/2.(2,132)なる式で与えられる

ことがわかっています。(1/τ₀)は衝突頻度(単位時間当りの衝突回数)

ここで,ｄは衝突する２原子の中心間距離(＝衝突径数)であり.Ｍ

は気体の１原子の質量です。βはＴを絶対温度ｋ_ＢをBoltzmann定数

として,β＝1/ｋ_ＢＴを意味します。

※(注13-1):これらの根拠を説明します。

まず,時刻:に自由飛行していた１個の原子が時刻τ＋ｄτに衝突

せずに自由飛行で残っている確率は,(1－dτ/τ₀)倍となります。

そこで,時刻τに自由運動している原子の個数をＮ^(τ)とすると,

その変化は,ｄＮ^＝Ｎ^(τ＋ｄτ)－Ｎ^(τ)＝－(ｄτ/τ₀)Ｎ＾(τ)を

満足します。そこで,これを積分するとＮ^(τ)＝Ｎ^(0)exp(－τ/τ₀)

です。

気体原子総数をＮとして,時刻τに衝突せず残っている確率ｐ(τ)は,

ｐ(τ)＝Ｎ^(τ)/Ｎで与えられ,∫₀^∞Ｎ^(τ)ｄτ＝Ｎ,or ∫₀^∞ｐ(τ)ｄτ

＝1となる規格化条件により,指数分布:ｐ(τ)＝(1/τ₀)exp(－τ/τ₀)

が得られ,飛行時間τの平均は,確かに,＜τ＞＝­∫τｐ(τ)ｄτ＝τ₀

となります。そして,次に.この平均自由飛行時間τ₀を評価します。

対象の気体を構成する気体原子(分子)の数密度をｎ＝Ｎ/Ｖとする

とき,ある１原子(分子)が自由直線運動で長さＬを通過する間に,

それと衝突することが可能な距離が衝突径数ｄであり,これを半径

として長さがＬの体積:πｄ²Ｌの円筒内にあるの原子(分子)数は

ｎπｄ²Ｌなので,この１原子の,この間の衝突回数がｎπｄ²Ｌに等しい

ことになります。

この原子(分子)の平均速度を＜ｖ>とすると,その飛行時間がτなら,

平均飛行距離は,Ｌ＝＜ｖ＞τです。

そこで,平均衝突回数:ｎπｄ²Ｌ＝ｎπｄ²＜ｖ＞τが丁度１となる

飛行時間をτ＝τ₀と定義すれば,τ₀＝1/(ｎπｄ²＜ｖ＞),または,

1/τ₀＝ｎπｄ²＜ｖ＞が成立します。

このとき,Ｌ₀＝＜ｖ＞τ₀とすると,これは１原子(分子)が全く衝突

せずに飛行できる平均距離を示しているので平均自由行程と

呼ばれます。ところが.希薄な気体分子(原子)は.Maxwellの速度

分布:f(ｖ)ｄ³ｖ＝{Ｍ/(2πｋ_ＢＴ)}^3/2exp(―Ｍｖ²/2ｋ_ＢＴ)ｄ³ｖに従う

ことがわかっているため,平均速度の大きさ＜ｖ＞は,ｆ(ｖ)が速度ｖ

空間で球対称故,＜ｖ＞＝∫₀^∞ｖｆ(ｖ)ｖ²ｄｖ/∫₀^∞ｖｆ(ｖ)ｖ²ｄｖなる計算式

から得ることができます、このｆ(ｖ)に,上のMaxwell分布式を代入

すれば,結局,＜ｖ＞＝(8ｋ_ＢＴ/πＭ)^1/2＝{(8/(πβＭ)}^1/2

が得られます。

ただし,厳密には同じ質量Ｍの原子の２体衝突なので,1個の質量Ｍ

ではなく,これを相対運動の換算質量Ｍ^＊＝Ｍ/2に置換する必要

があり,＜ｖ＞＝{8/(πβＭ^＊)}^1/2＝4/(πβＭ)^1/2が正しい評価式です。

これによって1/τ₀＝ｎπｄ²＜ｖ>＝4ｄ²ｎ/{π/(βＭ)}^1/2(ｎ＝Ｎ/Ｖ)

という求める表式が得られました。

以下では,上の計算で実際に＜ｖ＞＝{8/(πβＭ^＊)}^1/2

＝4/(πβＭ)^1/2が得られる,数学的根拠を示しておきます。

まず,速度ベクトルｖの極座標表示を(ｖ,θ,φ)(ｖ＝｜ｖ｜)として

ｖに対する確率密度が,ｐ(ｖ)＝Ｃｆ(ｖ,θ,φ)で与えられる場合を

考えます。確率密度の満たすべき条件:∫ｐ(ｖ)ｄ³ｖ＝1により,

規格化定数:Ｃは,Ｃ＝1/[∫ｆ(ｖ,θ,φ)ｖ²ｄｖ(sinθ)ｄθｄφ]で

与えられます。特にｆが球対称でｆ(ｖ,θ,φ)＝ｆ(ｖ)と表わせる場合,

規格化定数は,Ｃ^-1＝∫ｆ(ｖ,θ,φ)sinθｖ²ｄｖｄθｄφ

＝4π∫₀^∞ｆ(ｖ)ｖ²ｄｖで決まります。

この場合,速度の大きさ:ｖ＝|ｖ|の平均値は,

＜ｖ＞＝＜|ｖ|＞＝∫ｖｐ(ｖ)ｄ³ｖで与えられますが,

これは,＜ｖ＞＝Ｃ∫ｖｆ(ｒ,θ,φ)ｖ²ｄｖ(sinθ)ｄθｄφ

＝4π∫₀^∞{ｖｆ(ｖ)}ｖ²ｄｖ/{4π∫₀^∞ｆ(ｖ)ｖ²ｄｖ}

＝∫₀^∞ｖ³ｆ(ｖ)ｄｒ/{∫₀^∞ｖ²ｆ(ｒ)ｄｒ}

で得られることになります。

そこで,ｆ(ｖ)＝Ｃexp(－αｖ²)と表わされる場合には,

＜ｖ＞＝∫₀^∞ｖ³exp(－αｖ²)ｄｖ/∫₀^∞ｖ²exp(－αｖ²)ｄｖと

なります。後は分子,分母の指数関数を含む積分の評価だけです。

まず.Ｉ(α)＝∫₀^∞exp(－αｘ²)ｄｘと置くと,Gauss分布の積分

公式から,Ｉ(α)＝(1/2)(π/α)^1/2です。

これから,－ｄＩ/ｄα＝∫₀^∞ｘ²exp(－αｘ²)ｄｘ＝(1/4)π^1/2α^-3/2

を得ます。

一方,Ｊ(α)＝∫₀^∞ｘexp(－αｘ²)ｄｘとして,これを求めます。

ｕ＝αｘ²と変数置換すると,ｄｕ＝αｘｄｘより,

Ｊ(α)＝∫₀^∞ｘexp(－αｘ²)ｄｘ

＝(1/2)α^－1∫₀^∞exp(－ｕ)ｄｕ＝(1/2)α^－1となります。

そこで,I(α)の場合と同じく両辺をαで微分すると,

－ｄＪ/ｄα＝∫₀^∞ｘ³exp(－αｘ²)ｄｘ＝(1/2)α^－2を得ます。

以上から,＜ｖ＞＝∫₀^∞ｖ³exp(－αｖ²)ｄｖ

/{∫₀^∞ｖ²exp(－αｖ²)ｄｖ＝(1/2)α^－2/{(1/4)π^1/2α^-3/2}

＝2(πα)^-1/2となることがわかりました。

これに,Maxwell速度分布:f(ｖ)ｄ³ｖ＝{Ｍ^＊/(2πｋ_ＢＴ)}^3/2

exp(―Ｍｖ²/2ｋ_ＢＴ)ｄ³ｖを適用するために,

ｆ(ｖ)＝Ｃexp(－αｖ²),α＝Ｍ^＊/(2ｋ_ＢＴ),(Ｍ^＊＝Ｍ/2)を

代入することから,＜ｖ＞＝2(πα)^-1/2＝(8ｋ_ＢＴ/πＭ^＊)^1/2

＝4/(πβＭ)^1/2が得られました。

したがって,結局,自由飛行時間τ₀の求める表式,

1/τ₀＝ｎπｄ²＜ｖ>＝4ｄ²ｎ/{π/(βＭ)}^1/2が得られました。

(注13-1終わり※)

さて,原子のエネルギー準位と波動関数に及ぼす衝突の

効果は甚だ複雑です。

衝突する2原子の相互作用によって,エネルギー準位が衝突期間

中にシフトし,波動関数は非摂動原子の波動関数の,ある１次結合

になります。衝突している期間が十分短かいなら,その期間内の

光の吸収,放出は無視できます。そこで,以下では,衝突期間は

平均自由飛行時間τ₀に比較して,極く短かいと仮定します。

衝突には,「非弾性衝突」i.e.原子の状態があるエネルギー

準位から,他のエネル.ギー準位に変化するものがあります。

この効果は,原子占位数の減衰の速さを,もう１つ補うこと

で表わされ.光学Bloch方程式には,放射による減衰速度γの

他に,衝突による減衰速度を適当に加えた形で現われるもの

があります。

一方,「弾性衝突」では,原子は前と同じエネルギー準位に

留まり,その効果は原子波動関数の位相(phase)を変える

だけです。これはフェイズシフト(phase-shift)と呼ばれます。

波動関数:Ψ(ｒ,ｔ)＝Ｃ₁(ｔ)Ψ₁(ｒ,ｔ)＋Ｃ₂(ｔ)Ψ₂(ｒ,ｔ)

の係数:Ｃ₁とＣ₂の位相が変われば,原子密度行列:ρ_ij＝Ｃ_iＣ_j^＊

のρ₁₂＝Ｃ₁Ｃ₂^＊とρ₂₁＝Ｃ₂Ｃ₁^＊は変化しますが,対角要素である

ρ₁₁＝Ｃ₁Ｃ₁^＊＝|Ｃ₁|²と,ρ₂₂＝Ｃ₂Ｃ₂^＊＝|Ｃ₂|²は変化しません。

かくして,弾性的位相遮断衝突で(2.115)の,ｄρ₁₂/ｄｔ

＝(iΩ/2)exp{－i(ω₀－ω)ｔ}(ρ₁₁－ρ₂₂)－γρ₁₂には,減衰

速度が余分に入り込みますが,対角要素の方程式(2.114)の,

ｄρ₂₂/ｄｔ＝(－iΩ^＊/2)exp{i(ω₀－ω)ｔ}ρ₁₂

＋(iΩ/2)exp{－i(ω₀－ω)ｔ}ρ₂₁－2γρ₂₂,の方は不変です。

これらは,広範囲の物理的状況に対して線幅を広げる際だった

効果を及ぼし,また,第３章の時間に依存する光の性質の議論では

非常に便利なモデルを提供してくれます。

ここで,弾性衝突による非対角要素の減衰効果をγ_collと記すこと

にします。この量は,明らかに衝突頻度:1/τ₀に比例し,§3.3で,

実際にγ_coll＝1/τ₀.(2.133)となることを示しますが,これと類似の

関係:2γ＝Ａ₂₁＝1/τ₀.(2.102)との間には因子２の違いがあること

には注意が必要です。

なお,以下では,衝突広がりの話は弾性衝突に限ることにします。

すると,光学Bloch方程式:(2.115)は,次のように修正されます。

つまり,ｄρ₁₂/ｄｔ＝(iΩ/2)exp{－i(ω₀－ω)ｔ}(ρ₁₁－ρ₂₂)

－γ^ρ₁₂.(2.134)です。γ^＝γ＋γ_coll..(2.135)です。

この結果,以前「光の量子論11」で求めた定常状態の式

(2.119)のρ₂₂＝(|Ω|²/4)/{(ω₀－ω)²＋γ²＋|Ω|²/2},

および,(2.120)のρ₁₂＝－exp{－i(ω₀－ω)ｔ}

×[(Ω/2)ω₀－ω－iγ)/{(ω₀－ω)²＋γ^²＋|Ω|²/2}

は,ρ₂₂＝{(γ^/γ)|Ω|²/4}

/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}.(2.136),および,

ρ₁₂＝－exp{－i(ω₀－ω)ｔ}[(Ω/2)ω₀－ω－iγ＾)

/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}(2.137)に一般化

されます。

※(注13-2);上記を証明します。

[証明]:前の(注11-1)での(2.119)(2.120)の証明と同じく,

ρ₁₂,ρ₂₁を,ρ~₁₂＝exp{i(ω₀－ω)ｔ}ρ₁₂,

ρ~₂₁＝]exp{－i(ω₀－ω)ｔ}ρ₂₁と置換して振動型因子を消去

した３方程式系から出発します。前は.

ｄρ₂₂/ｄｔ＝(－iΩ^＊/2)ρ~₁₂＋(iΩ/2)ρ~₂₁－2γρ₂₂.

ｄρ~₁₂/ｄｔ＝(iΩ/2)(1－2ρ₂₂)－γρ~₁₂－i(ω₀－ω)ρ~₁₂,

ｄρ~₂₁/ｄｔ＝(－iΩ^＊/2)(1－2ρ₂₂)－γρ~₂₁＋i(ω₀－ω)ρ~₂₁

でしたが,これらが下の２つでは,γをγ＾に変更した方程式系:

ρ₂₂/ｄｔ＝(－iΩ^＊/2)ρ~₁₂＋(iΩ/2)ρ~₂₁－2γρ₂₂,

ｄρ~₁₂/ｄｔ＝(iΩ/2)(1－2ρ₂₂)－γ^ρ~₁₂－i(ω₀－ω)ρ~₁₂,

ｄρ~₂₁/ｄｔ＝(－iΩ^＊/2)(1－2ρ₂₂)－γ^ρ~₂₁＋i(ω₀－ω)ρ~₂₁

に修正されます。

定常状態では.これら３方程式の左辺の変化速度を全てゼロと置き,

前の(注11-1)の証明と同じ道をたどって,修正しようと思います。

まず,ｄρ₂₂/ｄｔ＝0から,ρ₂₂＝{－i/(4γ)}{Ω^＊ρ~₁₂－Ωρ~₂₁}

(実数)を得ます。これは前と同じです。

一方,ｄρ~₁₂/ｄｔ＝0,または,ρ~₂₁/ｄｔ＝0より

(－Ω^＊/2)(1－2ρ₂₂)＝(ω₀－ω－iγ^)ρ~₂₁

(Ω/2)(1－2ρ₂₂)＝－(ω₀－ω＋iγ^)ρ~₁₂,

または,(|Ω|²/2)(1－2ρ₂₂)＝－(ω₀－ω＋iγ^)Ω^＊ρ~₁₂

です。

ここで,Ω^＊ρ~₁₂＝ａ＋iｂ(ａ,ｂは実数)と置くと,

(|Ω|²/2)(1－2ρ₂₂)＝－ａ(ω₀－ω)＋ｂγ^,かつ,

0＝－ｂ(ω₀－ω)―ａγ^を得るため,

ａ＝(－|Ω|²/2)(1₁－2ρ₂₂)[(ω₀－ω)/{(ω₀－ω)²＋γ^²}],

かつ,ｂ＝(－|Ω|²/2)(1₁－2ρ₂₂)[－γ/{(ω₀－ω)²＋γ^²}],

です。これから,前には,いくつかの計算手続きの結果として,.

Ω^＊ρ~₁₂－Ωρ~₁₂^＊＝２iｂ＝i|Ω|²(1－2ρ₂₂)

×[γ/{(ω₀－ω)²＋γ²}],を得ましたが,今回の修正では,

これが,Ω^＊ρ~₁₂－Ωρ~₁₂^＊＝２iｂ

＝i|Ω|²(1－2ρ₂₂)[γ^/{(ω₀－ω)²＋γ^²}]に変わります。

これを最初の式:ρ₂₂＝{－i/(4γ)}{Ω^＊ρ~₁₂－Ωρ~₁₂^＊}

に代入すれば,4iγρ₂₂＝iγ^|Ω|²/{(ω₀－ω)²＋γ^²}

－2iγ^ρ₂₂|Ω|²/{(ω₀－ω)²＋γ^²}となるため,

4iγρ₂₂{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}/{(ω₀－ω)²＋γ^²}

＝iγ^|Ω|²/{(ω₀－ω)²＋γ^²}を得ます。

それ故,ρ₂₂＝{(γ^/γ)|Ω|²/4}

/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}が得られます。

よって,1－2ρ₂₂

＝{(ω₀－ω)²＋γ^²}/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}

ですから,Ω^＊ρ~₁₂＝(－|Ω|²/2)(1－2ρ₂₂)

×[(ω₀－ω－iγ^)/{(ω₀－ω)²＋γ^²}]

＝(－|Ω|²/2)(ω₀－ω－iγ^)

/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}です。

故に,ρ~₁₂＝(－Ω/2)(ω₀－ω－iγ^)

/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}ですから,

ρ₁₂＝－exp{－i(ω₀－ω)ｔ}[(Ω/2)(ω₀－ω－iγ^)

/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}も得られます。

[証明終わり］　(注13-2終わり※)

ここまでくれば,感受率χを導出し直すのは簡単であり,

(2.122)のχ(ω)＝{Ｎｅ²|Ｄ₁₂|²/(3ε₀ｈ_cＶ)}

×[(ω₀－ω＋iγ)/{(ω₀－ω)²＋γ²＋|Ω|²/2}]

の一般化は,χ(ω)＝{Ｎｅ²|Ｄ₁₂|²/(3ε₀ｈ_cＶ)}

×[(ω₀－ω＋iγ^)/{(ω₀－ω)²＋γ^²＋(γ^/γ)|Ω|²/2}]

(2.138)　になることがわかります。

こうして,放射,飽和,衝突による広がりを含む原子遷移の

線幅は,(2.123)の2(γ²＋|Ω|²/2)^1/2を一般化した表式:,

2{γ^²＋(γ^/γ)|Ω|²/2}^1/2.(2.139)となります。

今後のほとんどでは.飽和広がりが無視できて,上の(2.139)

が.2γ＾＝2γ＋2γ_coll.(2.140)に帰着するような,かなり弱い

ビームについてのみ議論します。

　この場合,吸収線の形はLorentz型であり,幅は放射による寄与

と衝突による寄与を単に加えることで得られます。

衝突時間τ₀の典型的な値は室温で10 ⁵ Pascalの圧力に相当する

気体密において,τ₀ ～3×10 ^-11ｓ.(2.141) です。

この値に対する衝突による線幅:2γ_coll＝τ₀^-1×2 ～ 6.7×10 ¹⁰ s^-1

は,(2.57)の角周波数の線幅:2γ＝Ａ₂₁ ～ 6×10 ⁸ｓ^-1のＡの大きさ

から導かれる周波数の自然線幅:Δｆ＝Δω/(2π)～10 ⁸ Hzという,

放射による線幅の,約100倍で与えられます。

• 2.12 Doppler(ドプラー)広がり

気体内の原子の速度には,ある広がりがあるので,原子が吸収,

放出できる光の周波数にはドプラー効果を通じて,それに対応

したある分布が生じます。当面は,線幅を広げる他の原因を全て

無視し,ｚ軸に平行に伝播する角周波数ωの光の吸収を考察します。

　低い準位Ｅ₁にある原子の速度がｖ₁であるとき,光:ｈ_cωを吸収

して高い準位:Ｅ₂に上がります。ｋを,大きさがｋ＝ω/ｃでｚに

向かうベクトルとすると.その光の運動量はｐ＝ｈ_cｋであり,Ｅ₂に

おける原子の速度はｖ₂となります。

原子質量をＭとすると,Ｍｖ₁＋ｈ_cｋ＝Ｍｖ₂.(2.142).かつ,

Ｅ₁＋(1/2)Ｍｖ₁²＋ｈ_cω＝Ｅ₂＋(1/2)Ｍｖ₂².(2.143)です。

特に,ｖ₁＝0,ｖ₂＝0のケースの光子周波数をω₀とすると,

ｈ_cω₀＝Ｅ₂－Ｅ₁.(2.144)です。.

故に,ｈ_cω₀＝ｈ_cω＋(Ｍ/2)(ｖ₁²－ｖ₂²)です。

一方,(2.142)からｖ₂＝ｖ₁＋(ｈ_cｋ/Ｍ)なので,

ｖ₂²={ｖ₁＋(ｈ_cｋ/Ｍ)}²＝ｖ₁²＋ｈ_c²ｋ²/Ｍ²＋2ｈ_cｖ₁ｋ/Ｍ

なので,すぐ前の式に代入するとｖ₂が消去されて,

ｈ_cω₀＝ｈ_cω－ｈ_cｖ₁ｋ－ｈ_c²ｋ²/(2Ｍ).(2.145)となります。

ところが,ｋはｚ軸に平行でｋ＝|ｋ|＝ω/ｃなので,

ｖ₁ｋ＝ｖ_1ｚ(ω/ｃ)であり,ｈ_c²ｋ²＝ｈ_c²ω²/ｃ²です。

したがって,ω₀＝ω－(ｖ_1ｚ/ｃ)ω－ｈ_cω²/(2Ｍｃ²)

(2.146)を得ます。

この右辺の量の代表的な大きさのオーダーは,ｖ_1z/ｃ

～10 ^-5より,ｈ_cω/(2Ｍｃ²)～10 ^-9.(2.147)ですから,

ωとω₀の差は僅かで,右辺最後の項は無視できるため,

ω＝ω₀/{1－(ｖ_1ｚ/ｃ)}～ ω₀{1＋(ｖ_1ｚ/ｃ)}.(2.148)

という形になります。

(※ つまり,ｖ_1ｚ ～ｃ(ω－ω₀)/ω₀です。)

すなわち,この光は周波数ωが原子の初速に応じた

ドプラー・シフトだけ,ω⁰からずれている場合に限り,吸収可能

であることを示しています。

したがって,気体内の全原子の吸収する光の周波数ωの分布は,

まず,原子速度のｚ成分がｖ_ｚとｖ_ｚ＋ｄｖ_ｚの間にある相対確率

が,絶対温度がＴのとき,Maxwellの速度分布則により,

exp{－Ｍｖ_ｚ²/(2ｋ_ＢＴ)}ｄｖ_ｚとなるため,これにｖ_ｚ＝ｃ(ω－ω₀)

を代入することによってexp{－Ｍｃ²(ω－ω₀)²/(2ω₀²ｋ_ＢＴ)}

×(ｃ/ω₀)ｄω.(2.149)で与えられることがわかります。

この吸収光の周波数分布はGauss型曲線として知られています。

Gauss曲線のピークはω＝ω₀にありますが,この曲線では,

exp{－Ｍｃ²(ω－ω₀)²/(2ω₀²ｋ_ＢＴ)}＝1/2.(2.150)を満たす

周波数のところ,つまり,ω＝ω₀±ω₀{2ｋ_ＢＴln2/(Ｍｃ²)}^1/2

のとき,強度が最大値の半分になります。。

　よって,この半値幅:2Δがドプラー広がりの全幅ですが,これ

は,2Δ＝2ω₀{2ｋ_ＢＴln2/(Ｍｃ²)}^1/2.(2.151)と表わされます。

この幅は,一般のGauss分布の平均２乗幅:δを用いて表わす

こともできます。すなわち,δとΔの関係が,同じオーダーで

δ＝ω₀{2ｋ_ＢＴ/(Ｍｃ²)}^1/2＝Δ/(2ln2)^1/2～Δ/1.18.(2.152)と

表わせます。

すぐ前の衝突広がりの節では,幅が2γ_collで,これを衝突頻度

で,γ_coll＝τ₀^-1と表わすとき,2.141)のτ₀ ～3×10 ^-11ｓを得た

のと同じく,ドプラー広がりの幅のパラメータもΔ＝τ₀^-1として

評価した場合,ドプラー幅は,ほぼ,衝突幅に等しいと評価されます。

Gauss曲線は,Ｆ_Ｇ(ω)＝(2πδ²)^-1/2exp{－(ω－ω₀)²/(2δ²)}

(2.153)で定義されますが,このとき,∫_－∞^∞Ｆ_Ｇ(ω)ｄω＝1(2.154)

が満足されます。

これに対して,前の放射広がりや衝突広がりで現われた周波数

分布の形状のLorentz型曲線は,Ｆ_Ｌ(ω)＝(γ/π)/{(ω－ω₀)²＋γ²}

で定義されますが,これも∫_－∞^∞Ｆ_Ｌ(ω)ｄω＝1の条件を満たすよう

に係数が決められています。

今回も§2.12がここで終わるので,ここまでにします。(つづく)

(参考文献)：Rodney Loudon 著

(小島忠宣・小島和子 共訳)

「光の量子論第２版」(内田老鶴舗)