~ Quantum mechanics of elastic colliding particles ~
簡単なところから始めるために,壁に衝突して完全反射する1粒子の量子力学を考えることにします.
まずはじめに,時刻 の初期状態の波動関数が
であったと仮定します.つまり,確率密度 は
というように, の範囲にだけ,ある一定の値を持ち,それ以外は であるという状態にあると考えます.このは規格化定数で,
でなければなりません.そして,この についてもう一つ仮定を付け加えます.
要するに, はを波数とした時の波長である,ということです.まとめると,初期状態の波動関数は
であるということです.
粒子の位置 の期待値は
で与えられます.この状態での値は
また, の期待値は
したがって,の標準偏差は
となります.
粒子の運動量の期待値は
で与えられるのですが,今考えている状態の波動関数は
だったので,波束の境目のところで微分ができません.そこで,「運動量表示の波動関数」というものを考えます.それは,通常の「座標表示の波動関数」から
で求められます.この関数は元の関数と
の関係にあります.要するに,フーリエ変換と,その逆変換の関係です.実際に計算すると,として
この運動量表示の波動関数から,運動量の確率分布が得られます.
fig. EC-1 ()
運動量の期待値は
で与えられますが,上の図より明らかに確率分布がをピークにして左右対称なので,運動量の期待値は
であることがわかります.
次に,運動量の2乗の期待値を考えますが,これは
という積分になります.この形のままではよくわからないので,変数変換してみます.とすれば,でなので
となります.このはとおけば
fig.EC-2 (青:第1項,橙:第3項,紫:1項と3項の和)
すると,二番目の項の被積分関数は奇関数なので積分した値は,三番目の項は積分値は有限で ,ところが一番目の項は積分が無限大に発散しています.ということは,運動量の二乗の期待値は
ということです.したがって,運動量の標準偏差も無限大
ということになります.
それでは,初期状態が指定された後の時刻で状態がどのようになるか考えていきます.手始めに最も簡単な,ポテンシャルのない自由粒子の場合を考えます.シュレーディンガー方程式は
ですので,たとえば初期状態が
であれば,
と仮定して
すなわち
となります.また,たとえば初期状態が
であれば
となります.したがって
であれば
なので
そこで,とおいてこれを計算していくと,一般に
なので,上の積分では
とおけて
となるので
となります.もし,時間の変数を とおけば
とシンプルに書けます.
任意の関数はδ関数を使って
と表せますから,初期状態の波動関数も
と書けます.これは で表される状態の重ね合わせと見ることができます.ある一つの で表される状態は,時間 の後に
という状態に変わるのですから,任意の初期状態の波動関数の時間 の後の状態もこれの重ね合わせの状態になります.つまり
となります.ここで
と定義すれば,自由粒子の波動関数は,初期状態が与えられれば
で求められる,ということになります.
ここで求めた積分方程式
は,時間を2分割してみると,
ですから,
一方,当然
なので関数について
という,重要な性質があることがわかります.これはもともとシュレーディンガー方程式に波動関数の初期状態を関数とした時の解がであることから導かれるので,自由粒子でなくても,一般にポテンシャルがある場合のについても成り立つ話です.
さて,今は時間を2分割しましたが,少し考えると,これは何分割しても同じように成り立つものであることがわかります.例えば分割した場合の式は
のようになります.この式のイメージは,時間を で等分した場合で考えると分かりやすいです.
ここで
とおくと,
すると, より指数関数の積の部分は
と一つにまとまります.それで,から までの多重積分を.
とみなし.一つの記号
で表すことにします.
積分される関数は,自由粒子のラグランジアン
より,からに至る,ある経路に対する作用
の形になっているので,次のような式で表すことができます.
これが量子力学の経路積分の方程式です.今回の自由粒子の場合は,具体的には次のようになります.
ここで少し寄り道して,シュレーディンガー方程式と経路積分の方程式が同等なものであることを確認してみます.
以下,ブラ・ケット記法を使って考えます.
量子状態 を表すケット・ベクトル は,位置の演算子 の固有状態
を使って
と表せます.ここで,位置の固有状態は
ですので,波動関数は
です.あらためて
と書けば,状態 に対するシュレーディンガー方程式も
のようになります.この方程式を解くと
となります.このは初期状態のことですから,初期状態が位置の固有状態の場合
この状態が時刻の時に位置にある確率振幅は,つまり波動関数は
これは時刻に位置の固有状態にあったものが,時刻に位置の固有状態に遷移する確率振幅を意味しているので,
とかけます.これが前節までにもとめた積分核
で
となります.そこで,前節のを考えてみましょう.
この式の各々のが
となります.そこで,運動量の固有状態を考えると,これは
というような完全系をとるので,ハミルトニアンが位置と運動量の演算子の関数だったものが
のように,ハミルトニアンを位置と運動量の固有値の関数に置き換えることができます.すると,なので
したがって,
が得られます.ここでハミルトニアンは
の形をしているものと仮定すると
すると,一般に
なので
以上より,
が得られます.ハミルトニアンに対応するラグランジアンは
なので,定数を
として書き直すと
となります.この結果より,
そして,から までの多重積分を.
とみなし.一つの記号
で表すことにすれば,からに至る,あるひとつの経路に対する作用は
でしたから,
は次のように表します.
波動関数の時間発展を表す方程式
の積分核をこのように求めることを経路積分の方法と呼びます.
すでに求めた自由粒子の積分核は
でしたが,自由粒子の古典力学で実現する経路の作用,いわゆる古典的作用が
だったことを考えると,あらゆる経路について積分したにも関わらず,古典的経路の作用の成分しか効いていないようにみえます.
これはどういうことかというと,もともと経路積分は「可能な限りすべての経路」について
だったのですが,古典的経路で作用が停留値を持つ時(つまり,古典的に方程式が解けて軌道が求まる時)は,
の因子が停留値とその周辺の軌道の分だけが,重ね合わせの干渉によって残る成分となり,停留力大きく離れたところの軌道とその周辺の軌道との干渉は,位相の変化が速いため干渉によって実質相殺されてしまうのです.
作用を古典的作用とその変分に分けて
と書くと,
となります.この式のは途中の経路に関係なく得られます.ポテンシャルの存在によっての項が重ね合わされて干渉するのですが,こと自由粒子に限って言うと,
この中のが古典経路の場合常にゼロなので
となります.したがって自由粒子の場合は積分核が
に等しくなるのです.
そこで,いよいよで のポテンチャルがあった場合を考えます.この場合古典経路が下の図のように2本あります.
一つは自由粒子と同じもの,もう一つは壁で速度が反転して に到達する経路です.
すると,積分核全体は二つの和に分解できて
このは自由粒子の積分核そのもので
もう一つのは壁で反射した時間 で分割した二つの自由粒子の積分核の積に位相の変化の因子がかかったもので表されます.
この時,反射した前後での速度は
となっています.
ところが,
なので,
ですから,壁の反対側の仮想的な場所に向かう自由粒子の積分核が出てきます.つまり,は
であることがわかります.ここに出てきた定数 を定めるには つまり,ちょうど壁にぶつかるところへの積分核を考えれば良いです.このとき波動関数の絶対値は必ずゼロになるので,もゼロになります.ということは,
でなければなりませんから,この定数は
だったことがわかります.
以上をまとめると,で絶壁ポテンシャルがある場合の経路積分の積分核は自由粒子の積分核
を用いて
で表される,ということです.
前回求めた結果を,波数による積分の形で求めてみましょう.
これは物理的には,例えば与えられたハミルトニアンに対するシュレーディンガー方程式で,エネルギー固有値と固有関数 が求められた場合で考えますと,
この時,任意の波動関数は
であらわされます.ここで,係数のは
で求められます.これをブラ・ケット記法で書くと
また, より任意の波動関数は
で表わされます.
さて,経路積分の積分核は
でした.上の式でとすると,
つまり,積分核は
というように,エネルギー固有値と固有関数を使った和で表されます.
このように固有値が離散的な場合は和で表されますが,固有値が連続的な値を取る場合,つまりシュレーディンガー方程式が連続的なパラメータを使って
となる場合は,積分核は
のように,パラメータによる積分で表されます.
すると,自由粒子の場合は として
ですから,自由粒子の積分核は,
となって,前回求めた積分核が
として得られます.この式をもう少し計算すると
となりますから,として,エネルギーの固有値はで,そのエネルギーの固有状態はであるといえます.この式を見ると, で確かにはゼロになっていることがよくわかります.
それでは§1で考えた,初期状態の波動関数が
である場合の時間発展について,経路積分の式を使って計算してみたいと思います.ただ,この式のままではの場合,壁から離れる方向に動いていくので,初期状態の波動関数は
に変更して考えていきます.それで,波動関数の時間発展を表す方程式は
で,このは
として
で表されます.また,
なので,一般に
が成り立ちます.
それでは
を計算していきます.1項目を
ここで,
したがって,
また,
なので,2項目を
ここで,
なので
さて,これらの積分はどう計算したら良いものだろうか?
フレネル積分とは次のように定義された関数である.
この図の青線が ,赤線が である.
ここで,
のフーリエ逆変換,すなわち
を考えると,
なので
となる.
より一般的に,
のフーリエ逆変換
を考えると, と置き直して
この積分は,関数
に等しいので,まとめると
である.
これで,§10の積分が実行できそうだ.
§10の結果を改めて書くと
ここで,
とおけば,求めるべき任意の時刻の波動関数は
で表わされて,以下のようになる.
これらの項はそれぞれ,§11で求めた
で計算できる.ここでとはフレネル積分で
と定義された関数としている.
一つ一つ計算すると,
次の項は
その次の項は
最後の項は
あとは,これら4つの式をつなげれば,波動関数が求まるということになる.
§12の結果をまとめて一気に書くと
となる.
この関数の時間発展の様子を具体的に見てみよう.波動関数の初期状態のパラメータを
と定めて,その上で時刻を表す変数を
に書き換えよう.すると,この状態での任意の時刻の波動関数は
となる.
この関数のグラフをMathStudio(http://mathstud.io/) を使ってプロットしてみた.
う〜〜ん..波束がくずれ過ぎちゃって,壁に反射するイメージが全然出てないですね..
初期状態の選定が失敗だったようです.
次回は初期状態の波動関数を変えて,どうなるか見てみたいと思います.
前回の結果は予想以上に,パッとしないものでした.
そもそも初期状態を§1で定めた形にしたのは,「こうすれば計算が楽になるかな〜」という期待があってのことでした.計算は全然楽ではなかったので,自分としては不満が残っています.
そこで,初期状態の波動関数を別の形にして,もとめた積分核を使って波動関数の時間発展がどうなるかいろいろ試したいと思います.
初期状態の関数形はいくつか候補があるのですが,とりあえず試したい関数形が二つほどあるので,順番に調べたいと思います.
..さて,平面波を表す波動関数は,規格定数を無視すると
で,このとき運動量はの値を持つ固有状態です.
そこで,運動量がを中心にの成分が一様に分布している状態を考えます.
つまり,(以下,規格定数は無視し続けます)
この波動関数の絶対値の2乗,つまり確率密度は
となります.たとえばのときのグラフは
となります.確率密度はを中心におおよその範囲に分布しているのが見えます.
運動量の分布の標準偏差が,位置の分布の標準偏差がなので,この状態は
という特徴を持ちます.いま求めた波動関数は
で,これは,位置の平均がの点にある状態なので,任意の点を中心に分布している場合は,この関数を方向にずらした
(ここでは,ただの定数なので無くしました)
になります.
さて,この関数を初期状態として計算してもいいのですが,この関数はの領域で,ゼロではない値を持っています.それはよくないです.もっといけないのは,ちょうど壁のあるところで波動関数の絶対値は常にゼロになっていなくてはいけないのが,そうでないところです.
ここで,思い出してしてみると,そもそも任意の時刻の波動関数は積分核を使って,
で,計算されるのでした.
そこで,初期状態の波動関数が,別のある複素関数を使った奇関数(の半分のようなもの)であったと仮定します.
(ここではの領域も定義された関数とします)
すると,単純に上の式に代入してみると,こんな感じになります.
ところが,
なので,
ということになります.
まとめると,
初期状態の波動関数が
で与えられたとき,任意の時刻の波動関数は
として,
で求められます.そしてもし,関数が別の関数をつかって
で与えられていたのなら,任意の時刻の波動関数は関数を使って
で求められます.
ということなのですが,この方法が正しいなら,いろいろと応用ができそうです.
それでは,今回想定している初期状態での具体的な計算を,次回以降で行っていきます
次のように定義された関数
から作られる奇関数
を使って,初期状態の波動関数が
であるとき,
として,任意の時間の波動関数は
となりますが,これは
としても同じ結果が得られます.
注:ここでは波動関数の規格化は省略して考えています.
ここで,
とおけば,
ですので,まず,を求めましょう.
簡単のため
とおくと,
なので,
である.
しかしこれは,またしても「フレネル積分」になってしまっている.
せっかく,簡単に計算できる状態を探していたはずなのに,これではぜんぜん簡単ではないではないか.
もうこうなったら,初期状態は,最もシンプルな「ガウス分布」で考えるのがいいだろう.
数学が知りたいのではない.ザックリとでいいから,物理が知りたいのである.
気分をリセットして,計算が簡単なガウス分布で考えてみましょう.計算の基本となる 関数を
とます.もし波動関数がなら,それは
であり,位置の期待値が,運動量の期待値がのガウス分布で表わされる状態になります.
このとき確率分布は,を中心にだいたいの範囲にかたまっているので,
のように,付近の確率分布がほぼゼロとみなせるように,とを値が定められていれば,「近似的には初期状態はである」ような状態を設定して波動関数の時間発展が計算ができます.
さて,次のように定義された関数
から作られる奇関数
を使って,初期状態の波動関数が
であるとき,
として,任意の時間の波動関数は,の範囲で
となりますが,これは
としても同じ結果が得られます.
注:今回以降も波動関数の規格化は省略して考えています.
この波動関数は二つの部分の和
で表わされています.つまり
と
という状態の重ね合わせと「 の範囲では」同じものになります.このの式は,の式に対して
と符号を置き換えたもですので,実際の計算の手間は,
だけで済むでしょう.計算のために,とを
と置き換えます.すると,は
となります.これは,
となり,公式
を使えば,
として,
となります.それでは,この式のと を計算していきましょう.
なので
次に,
を計算します.
を計算しますが,その前に
としておきます.
また,
と書いておけば,は
すると,
であり,
なので,は
ここで,
と置き換えてみると
すると,
であり,
なので,
だったので,戻して書けば
となります.この式が求められれば,あとは簡単に進められます.
の式は,の式に対して
と符号を置き換えたものでしたから,
となります.したがって,求めるべき任意の時刻の波動関数はとして,
の範囲で
および,
の範囲で
となりました.(ただし,規格化定数は考えていません.)
最後に,この波動関数の時間発展の様子を具体的に観察すつるために,数値計算した様子を動画にしてみました.
パラメータの値はとしました.
画面中央の再生ボタンを押すと,私が手で時刻の数値を動かした時の波動関数の絶対値の2乗,つまり確率分布の挙動を見ることができます.
なお,今回も使用したツールはhttp://mathstud.io/です.
実際に使用したスクリプトはこちらにあります.
http://mathstud.io/qB32qM
時刻がスライドバーになっていますので,試しに動かしてみて下さい.
それではごきげんよう.