摂動法(その1)
量子閉じ込めシュタルク効果は量子井戸の量子準位が電界印加によって変化することによって起こることを前回説明しました。準位の変化は量子力学的に検証しなければなりません。準位が電界の影響によってどう変化するかを計算するにはシュレディンガー方程式を解かなければなりませんが、電界によるポテンシャルエネルギーの変化を加味すると解析的には解けなくなってしまいます。
そこで近似計算が行われるのですが、その近似法として使われるのが摂動法です。これまでこの方法については説明していませんので、この機会に紹介したいと思います。
摂動法というのは量子力学のために作られた方法ではなく、古くは天文学の星の軌道計算などに使われてきた方法です。例えば太陽の周りを回る地球の軌道は太陽と地球の間の万有引力によって大体決まりますが、厳密には他の惑星、金星や火星や木星などによる引力の影響をわずかながら受けて少し軌道が変わります。このわずかな変化はまともに複数の天体からの引力を考慮して厳密に計算するのは困難です。
そこで用いられたのが摂動法です。摂動は英語のpurtervationの訳語ですが、purterveとは乱すという意味です。つまり摂動法とはわずかな乱れを計算する近似方法という意味です。基本的には乱れがないとき、例えば地球と太陽しかない場合の解がわかっているとし、乱れが生じたとき、例えば火星の引力を考慮する場合、その乱れ、変化は火星の影響を考えない場合に対してわずかなものと考えます。そのときの解は乱れのない場合の分かっている解をもとに、それに微小な項を級数として加えた関数で表せるであろうと仮定します。そしてこの級数の係数を合理的に決定すれば、近似解が得られるという考え方です。
ここでは古典力学から解き明かすような一般的な方法にまで立ち返らず、量子力学のための近似方法を説明します。
時間に依存しない定常状態で、固有エネルギーは1つだけ、つまり縮退のない1次元の系を考えます。この場合のシュレディンガー方程式は波動関数をψ、固有エネルギーをE、ポテンシャルエネルギーをVとして
\psi&space;=0)
と書けますが、これを書き直して
\psi&space;=E\psi)
とすると、左辺のカッコ内がハミルトン演算子(ハミルトニアン)Hですから
(1)
となります。シュレディンガー方程式の解析解が得られ、Ψが既知であるような基準となる状態(例えばV=0としてもよいです)の(1)式を
(2)
と書きます。右肩の添え字(0)は基準状態(摂動のない状態)を表し、波動関数の添え字nは量子数(つまり(2)式の解は複数あるのでそれに番号をつける)で、面倒ですが後で登場する他の状態と区別するために付けてあります。
ここで外力が加わるなどしてポテンシャルエネルギーが変化したとします。(1)式よりハミルトニアンが変化しますから、これを
}+\lambda&space;\boldsymbol{H}{}')
と書きます。λは定数で0~1の値をとるとします。シュレディンガー方程式は
(3)
となります。ここでψnが求めたいけれども解が得られない波動関数です。λが小さいとすればこれを使ってψnとEnを展開した関数で求めたい解は近似できると考えます。その級数を
(4)
(5)
と書きます。これを(3)式に代入すると、
}\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}+\lambda&space;\left&space;(&space;\boldsymbol{H^{\left&space;(&space;0&space;\right&space;)}}\psi&space;_{n}^{\left&space;(&space;1&space;\right&space;)}+\boldsymbol{H{}'}\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}&space;\right&space;)+\lambda&space;^{2}\left&space;(&space;\boldsymbol{H^{\left&space;(&space;0&space;\right&space;)}}\psi&space;_{n}^{\left&space;(&space;2&space;\right&space;)}+\boldsymbol{H{}'}\psi&space;_{n}^{\left&space;(&space;1&space;\right&space;)}&space;\right&space;)+\cdot&space;\cdot&space;\cdot)
}&space;\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}&space;+\lambda&space;{\left&space;(&space;E_{n}^{\left&space;(&space;0&space;\right&space;)}}\psi&space;_{n}^{\left&space;(&space;1&space;\right&space;)}+E_{n}^{\left&space;(&space;1&space;\right&space;)}\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}&space;\right&space;))
(6)
となりますが、λの各次数について(6)式の両辺の係数は等しいはずですから、次数ごとに整理すると
<0次>
}}-E_{n}^{\left&space;(&space;0&space;\right&space;)}&space;\right&space;)\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}=0)
<1次>
(7)
<2次>
}-E_{n}^{\left&space;(&space;0&space;\right&space;)}\right&space;)\psi&space;_{n}^{\left&space;(&space;2&space;\right&space;)})
(8)
となります。0次は(2)式と同じ式になります。
つぎに固有エネルギーの第1近似を求めます。(7)式を用います。両辺に左から複素共役関数ψn(0)*をかけて積分し、スカラー積(内積)をつくります。すると(7)式左辺は
左辺
*}\left&space;(&space;\boldsymbol{H}&space;-E_{n}^{\left&space;(&space;0&space;\right&space;)}\right&space;)\psi&space;_{n}^{\left&space;(&space;1&space;\right&space;)}dx)
となります。エルミート演算子の性質を使って
*}\left&space;(&space;\boldsymbol{H}&space;-E_{n}^{\left&space;(&space;0&space;\right&space;)}\right&space;)\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}dx&space;\right&space;\}^{*})
と変形し、(2)式の関係を使えば
*}\left&space;(&space;E_{n}^{\left&space;(&space;0&space;\right&space;)}&space;-E_{n}^{\left&space;(&space;0&space;\right&space;)}\right&space;)\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}dx&space;\right&space;\}^{*}=0)
となり、左辺は0であることがわかります。
つぎに(7)式の右辺です。
右辺
*}\left&space;(&space;-\boldsymbol{H}{}'+E_{n}^{\left&space;(&space;1&space;\right&space;)}\right&space;)\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}dx)
*}&space;\boldsymbol{H}{}'\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}dx+E_{n}^{\left&space;(&space;1&space;\right&space;)}\int&space;\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)*}\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}dx)
*}&space;\boldsymbol{H}{}'\psi&space;_{n}^{\left&space;(&space;0&space;\right&space;)}dx+E_{n}^{\left&space;(&space;1&space;\right&space;)})
となるので、左辺が0であるので
(9)
が得られます。
長くなるので今回はここまでとします。
そこで近似計算が行われるのですが、その近似法として使われるのが摂動法です。これまでこの方法については説明していませんので、この機会に紹介したいと思います。
摂動法というのは量子力学のために作られた方法ではなく、古くは天文学の星の軌道計算などに使われてきた方法です。例えば太陽の周りを回る地球の軌道は太陽と地球の間の万有引力によって大体決まりますが、厳密には他の惑星、金星や火星や木星などによる引力の影響をわずかながら受けて少し軌道が変わります。このわずかな変化はまともに複数の天体からの引力を考慮して厳密に計算するのは困難です。
そこで用いられたのが摂動法です。摂動は英語のpurtervationの訳語ですが、purterveとは乱すという意味です。つまり摂動法とはわずかな乱れを計算する近似方法という意味です。基本的には乱れがないとき、例えば地球と太陽しかない場合の解がわかっているとし、乱れが生じたとき、例えば火星の引力を考慮する場合、その乱れ、変化は火星の影響を考えない場合に対してわずかなものと考えます。そのときの解は乱れのない場合の分かっている解をもとに、それに微小な項を級数として加えた関数で表せるであろうと仮定します。そしてこの級数の係数を合理的に決定すれば、近似解が得られるという考え方です。
ここでは古典力学から解き明かすような一般的な方法にまで立ち返らず、量子力学のための近似方法を説明します。
時間に依存しない定常状態で、固有エネルギーは1つだけ、つまり縮退のない1次元の系を考えます。この場合のシュレディンガー方程式は波動関数をψ、固有エネルギーをE、ポテンシャルエネルギーをVとして
と書けますが、これを書き直して
とすると、左辺のカッコ内がハミルトン演算子(ハミルトニアン)Hですから
となります。シュレディンガー方程式の解析解が得られ、Ψが既知であるような基準となる状態(例えばV=0としてもよいです)の(1)式を
と書きます。右肩の添え字(0)は基準状態(摂動のない状態)を表し、波動関数の添え字nは量子数(つまり(2)式の解は複数あるのでそれに番号をつける)で、面倒ですが後で登場する他の状態と区別するために付けてあります。
ここで外力が加わるなどしてポテンシャルエネルギーが変化したとします。(1)式よりハミルトニアンが変化しますから、これを
と書きます。λは定数で0~1の値をとるとします。シュレディンガー方程式は
となります。ここでψnが求めたいけれども解が得られない波動関数です。λが小さいとすればこれを使ってψnとEnを展開した関数で求めたい解は近似できると考えます。その級数を
と書きます。これを(3)式に代入すると、
となりますが、λの各次数について(6)式の両辺の係数は等しいはずですから、次数ごとに整理すると
<0次>
<1次>
<2次>
となります。0次は(2)式と同じ式になります。
つぎに固有エネルギーの第1近似を求めます。(7)式を用います。両辺に左から複素共役関数ψn(0)*をかけて積分し、スカラー積(内積)をつくります。すると(7)式左辺は
左辺
となります。エルミート演算子の性質を使って
と変形し、(2)式の関係を使えば
となり、左辺は0であることがわかります。
つぎに(7)式の右辺です。
右辺
となるので、左辺が0であるので
が得られます。
長くなるので今回はここまでとします。
この記事へのコメント