水素原子の「電子軌道」の波動関数,シュレディンガー方程式の解

E.関数R(ラゲール陪多項式)


関数R(r)の式(B-3)を解きます。
関数R(r)
再掲

まず並べ替え、整理します。

関数と変数を次のようにおきます。

これより次のようなことが言えます。
αとは
αは定数です。後ほどα^2のところで詳しく述べます。
ρとは
「A.2体問題」で出てきた円柱座標の変数ρとは全く関係ありません。ここではρ=2αrで、2点間の距離rの単純に2α倍の変数です。
式(E-1)の第1項はρとS(ρ)で次のように書き替えられます。

式(E-1)に代入します。

α^2を次のようにおきます。
エネルギーE=E[n, ℓ, m]はマイナス値です。E[n, ℓ, m]は、主量子数n、方位量子数ℓ、磁気量子数mで決まるエネルギー準位の値であり、定数です。E[n, ℓ, m]は、関数Rの変数rとは全く関係なく、他の2個の変数θ、φとも全く関係なく、定数です。換算質量μ、プランク定数hも定数です。つまりα^2とは、主量子数n、方位量子数ℓ、磁気量子数mで決まる値であり、定数です。
式(E-4)に、式(E-2)・ 式(E-3) ・式(E-5)を代入すると次のようになります。

さらに書き替えます。

次のようにおきます。

これを代入します。

両辺を4α^2で割ります。

第1項を書き替えます。

関数S(ρ)を、さらに次のような初項がゼロでないベキ級数の関数に置き換えます。

この関数S(ρ)をρで微分します。

もう一度ρで微分します。

式(E-6)にこれらを代入すると次のようになります。

式(E-9)は、ラゲールの陪微分方程式 (Associated Laguerre’s Differential Equation ) と呼ばれています。この式(E-9)を解いていきます。


上式はラゲールの微分方程式(Laguerre’s Differential Equation )です。このラゲールの微分方程式は、ルジャンドルの微分方程式ほど知られていませんが、掲載されている微分方程式の本もあります。このラゲールの微分方程式を経由して式(E-9)を解く方法もあります。しかし今回は、式(E-9)ラゲールの陪微分方程式を、級数を使った数学的帰納法で解いていきます。最初に漸化式を求めます。そして次に数学的帰納法により解いていきます。

漸化式を求めます。
L(ρ)は次のような無限級数とおきました。
再掲

ρで微分します。

さらにρでもう一度微分します。

これらを式(E-9)に代入すると次のようになります。

整理します。


ここでρ^νの係数だけをピックアップしてみる。
第1項はν=ν+1とおけばρ^νとなる。ρ^νの係数は次のようになります。

第2項もν=ν+1とおけばρ^νとなる。ρ^νの係数は次のようになります。

第3項と第4項はν=νとおけばρ^νとなる。ρ^νの係数は次のようになります。


ρ^ν乗で無限級数を書き替えます。

左辺の級数を展開し、ρの累乗別に整理すると次のようになります。

ρにかかわらず左辺が0になるためには、ρのすべての次数の係数が0でなければならない。一般項で表すと次のようになります。

ゆえに次のような漸化式が得られます。


有限性の確認
ここまでで、L(ρ)を無限級数とおき、漸化式は得られました。
最終的に得る波動方程式の解R(r)=S(ρ)=e(^-ρ/2)ρ(^ℓ) L(ρ)はρ→∞の時、有限つまり収束する必要があるわけですが、ここで確認してみます。
L(ρ) は無限級数ですがν→∞の時e(^ρ)に漸近し収束します。

R(r) はν→∞の時e(^ρ/2)ρ(^ℓ)となります。

しかしρ→∞の時、R(r) =e(^ρ/2)ρ(^ℓ) は、無限大になる。つまり発散し、このR(r)は波動方程式の解としては受け入れられない。


このままではR(r)は受け入れられないので少し変更します。
R(r)はe(^-ρ/2)ρ(^ℓ) 部分だけならρ→∞で収束します。R(r)の残りの級数部分L(ρ)もρ→∞で収束するようにします。無限項数の無限級数L(ρ)から、有限項数の有限級数L(ρ)に変えます。項数が有限であれば級数L(ρ)はρ→∞で収束します。

L(ρ)を有限級数に変えます。
漸化式において、λ-ℓ-1=ν’とすれば、有限な項数の級数となる。 ここで、書き換えをおこないν’=n’と書き換える。さらに次のようにλ=nと書き換えると、n-ℓ-1=n’となり、級数L(ρ)は次のように表されます。


最高次数の係数はa(ν’)= a(n’)= a(n-ℓ-1)、最高次数はρ(^ν’)= ρ(^n’)= ρ(^n-ℓ-1)となる。すべての項を展開すると次のように表されます。


このL(ρ)を、R(r)の式に代入し、再確認すると次のようになります。

この式の右辺の各項はρ→∞で個々に収束します。有限の項数なのでR(r)も収束する。 つまりL(ρ)を有限の項数としたときR(r)も有限であり波動方程式の解として受け入れられることがわかります。
最初に
「変数ρ[半径]の関数L(ρ)を、人為的に無限級数」とおきました。関数L(ρ)は、値の変化するものが、変数ρ[半径, 0~∞]と項数ν[1~∞; つまり無限の項数]の2個あります。
*ここでは私的に、項数を仮に変数として [ν]と表し関数の中に入れ、関数L(ρ)を、L(ρ, [ν])と表示します。

ここでの変更点
変数ρ[半径, 0~∞]部分は自然界の現象であり、変数ρは必ず0~∞を考える必要があります。
一方、項数νは最初に人為的に無限級数と仮定しました。ゆえに項数νは人為的に有限級数と変更可能です。
ここでは、「変数ρ[半径, 0~∞]と変数ν[項数, 1~ν’ ; つまり有限の項数]の関数L(ρ, [ν])を有限級数」と変更します.

現在の状態の再確認
有限な項数νの有限級数の関数L(ρ, [ν])とe(^-ρ/2)ρ(^ℓ)を掛けたもの、つまり波動関数R(r)は変数ρ[半径]が無限大∞まで大きくなっても発散しない。ゆえにR(r)=e(^-ρ/2)ρ(^ℓ)×L(ρ, [ν])は波動関数として受け入れられる。
数学的帰納法により、有限級数L(ρ)を求めます。
(i) まず最初に、最高次数 の係数 を次のように規定する。

(ii) 漸化式で、ν→ν-1と置き換えると次のように書き換えられる。


ここで、λ-ℓ-1=n’(=ν’)とすれば、a(ν’+1)= a(ν’+2)= a(ν’+3)・・・=0となり、 λ=nと書き換えれば、n-ℓ-1=n’となる。

最高次数ρ(^ν’)の係数aν’は、ν’=n’と書き換えるとaν’= an’となる。
つまり上の式でν’をn’と書き換えると次のようになる。

さらに、n’=n-ℓ-1を代入すると次のようになる。

ここで最高次数ρ(^ν’)= ρ(^n’)= ρ(^n-ℓ-1)の係数:

を代入すると次のようになる。


(iii) 漸化式で、ν→ν-2と漸化式で置き換えると…
(iv) 漸化式で、ν→ν-3と漸化式で置き換えると…
これらを整理して書き並べていくと次にようになる。

数学的帰納法により一般項は次にようになる。

これより級数は次のように表される。


k’=(n-ℓ-1)-kと置き換えて、初項と末項を入れ替えます。

このL (ρ)を、R(r)の式に代入すると次のようになる。

以上でrの波動関数R(r)が得られました。しかし規格化されていません。次に規格化します。


規格化された波動関数R(r)
関数R(r)の規格化の式は次のものです。 
再掲

R(r)は実数の関数でありR*(r)= R(r)である。ゆえに次のように簡略化できます。

ここで関数R(r)に規格化定数Nを追加し、規格化された関数R(r)を次のようにおき直します。

規格化の式に代入すると次のようになります。


定積分部分は次の値であることが知られています。


これを代入すると次のようになります。

ここでは規格化定数として次を採用します。


ゆえに規格化された関数R(r)は次のように得られます。

表現が異なる3種類を書いておきます。

●ρで表わした規格化された関数R
関数 を解いた直後の形。変数ρは主量子数nとボーア半径a0に対して比で表わされています。


●σで表わした規格化された関数R
変数σはボーア半径a0に対して比で表わされています。


●rで表わした規格化された関数R
α=Z/(na0)およびr=2αρより、変数をρ→rと戻す。
球面座標の変数r(電子と陽子の2点間距離)で、これが求めた関数 。




製本版をamazon.co.jpのPODで出版しました。
量子力学入門 水素原子の電子軌道 : シュレディンガー波動方程式 (改訂)
ここでは掲載できなかったラゲール(Laguerre)の規格化定数(normalization constant)を、
付録Ⅱ. 関数Rの規格化」 (pp159-166) に掲載しています。
 ISBN
 9784802095792
 

戻る<< 目次 >>次へ