場の量子論への準備１（自由電子）

『量子力学の初歩』より続き統一理論へつながる場の量子論を理解するために必要となる概念を理解することを目標にしよう。自由電子の波束の記述を通して位置と運動量の関係を理解し、相対性理論的量子論につなげる。

尚、各章節間で波動関数のグラフにおける振幅や、数式中のいくつかの係数に関して、あまり一貫性を考慮していないので、そのつもりで読んでいただきたい。

1. 自由電子と波束

1.1 計算機代数的アプローチ

量子力学の初歩ではポテンシャル井戸の中に束縛された電子の波動関数の表し方を理解し、原子の中も同様の仕組みで電子が束縛されていることがイメージ出来るようになった。

一方、霧箱の中を飛行機雲の様な、軌跡を残しながら電子が飛んでいく様子を観察することができる。ニュートン力学では『質量ｍの電子が外力の働かない状態で速度ｖ（つまり運動量p=mv）で運動する場合、その軌道はｘ＝ｖｔで表される』と簡単に言い表すことが出来る。しかしこれを以下の波動関数の形式で表すことは、一見困難に思える。

\[\small \Psi(x,t)=A e^{i (k x-\omega t)}\\\small \tt \tag{Eq1.1-1}\]

原子やポテンシャル井戸の狭い空間に閉じ込めらていた場合の解であった。従って関数の分布範囲が定まり正規化により振幅が求められる。なのでこの関数ψは空間に対する条件を与えて解かないかぎりAは求まらず”絵”に描くことが出来ない。しかし仮にAを適当に与えて描画し、動的な変化を考察することはこの関数ψの性質を理解する上で有効だろう。空間時刻t＝0における関数ψの空間分布をプロットすると

Fig1.1-1

このように左右に際限なく広がってしまい、電子の位置を示すことができない。当然、波数kを指定してプロットしたのだから、この図から波数を観測することにより運動量の方は知ることはできる。

この関数ψに光子の運動をあてはめてみよう。光子は質量が無いので、相対性理論によりエネルギーと運動量の関係はE=cp となる。そして光子の運動速度は絶えず光速度定数ｃとなる。つまり、

\[\small \begin{cases} i\hbar\partial_t\psi(x,t) =E=c p\\ \frac{\large \hbar }{ \large i}\small\partial_x\psi(x,t)=p \end{cases} \]

運動量pで等式をむすぶと、

\[\small \partial_t\psi(x,t)=-c \space \partial_x\psi(x,t) \tt \tag{Eq1.1-2} \]

の微分方程式を得る。つまり『光子の場合、波動関数の時間的変位はその点の空間微分（隣接点間の波動成分の傾き）により決まる』という事だ。

光子では光速度一定の原理により、 ω=2π, ν=c/λ, k=2π/λ の関係から、その波長、つまり光の色が決まると波数ｋと角周波数は ω=ck/λ と定まる。適当にAとｋを決めれば、波動関数Ae^i(kx-ωt)をプロットできる。自然単位系としてC=1 ,ℏ=1とすると、

\[\small A e^{i(k x-\omega t)}=e^{ik( x-t )} \]

以下の時計が並んだ様な図は、上段が Aeⁱ^k^(x-t) を時刻ｔで空間サンプリングしたもので、下段がt+Δtで空間サンプリングしたものを示している。

円の中の青いポインターが空間ｘの各点での波動の成分を表している。赤いポインタは以下の\[\small \partial_t\psi\ =-\partial_x\psi\]の計算結果を表している。以降この１サンプル分の時計の様な図を”複素円”と呼ぼう。

この例では、x軸上、左から右にかけて振幅は変わらず、各点でポインターの位相だけが+Δθ（左周りに増加）している。

↓Δｔ後

Fig1.1-2

Δt後、下段の各点の位相は上段での青ポインタは少しだけ右に回転している。さて、Eq1.1-2はサンプリング間隔⊿x、⊿ｔと指標 \(n_x, n_t=1,2,\dots\) を用いて以下の様に\(\psi(\space)\)から\(\psi[\space]\)へと離散化できる。\[\small \left\{\psi((n_t + 1)\Delta t,n_x \Delta x)-\psi(n_t \Delta t,n_x \Delta x)\right\}/\Delta t\\\small \quad= -\left\{\psi(n_t \Delta t,x+n_x\Delta x)-\psi(n_t\Delta t,n_x \Delta x)\right\}/\Delta x\\\small \downarrow\\\small \left\{\psi[n_t+1,n_x]-\psi[n_t,n_x]\right\}/\Delta t\\\small \quad=-\left\{\psi[n_t,n_x+1]-\psi[n_t,n_x]\right\}/\Delta x \]

と書ける。これを少し変形して、

\[\small \psi[n_t+1,n_x]-\psi[n_t,n_x]\\\small \quad\quad\quad\quad\quad=-\left\{\psi[n_t,n_x+1]-\psi[n_t,n_x]\right\}\Delta t/\Delta x \]

見づらいので以下の様にｘとｔの差分を示す関数を定義すると、

\[\small \psi_x'[n_t,n_x]\equiv \psi[n_t,n_x+1]-\psi[n_t,n_x],\\\small \psi_t'[n_t,n_x]\equiv \psi[n_t+1,n_x]-\psi[n_t,n_x]. \]

以下のように時間-空間の差分関係を表すことができる。

\[\small \psi_t'[n_t,n_x]=-\space\psi_x'[n_t,n_x]\Delta t/\Delta x \tt \tag{Eq1.1-3} \]

この式をFig1.1-2にあてはめてみよう。以下の図はFig1.1-2の拡大表示したものだ。左上端の複素円はFig1.1-2の1行6列の複素円から4個分を抜き出してある。

Fig1.1-3

図中の（1）から（3）までを説明すると、

(1)の上向緑ベクトルは右隣りの複素円との差分 ψ’_x を表している。各点の位相は少しずつ右に回っている。つまりψ’_x を表す緑ベクトルは空間方向に振った場合の位相の回転速度ベクトルと言える。
(2)は(1)のψ’_x(上向緑ベクトル)に -Δt/Δx を掛けて算出された時間差分を表している。この量は下段の下向緑ベクトルに対応し、図の通りΔt後の位相を右に回転させる。
(3)の赤ポインターは常に青ポインタ―（波動成分）と共に表示される時間差分の成分であり、下段に記載した下向緑ベクトルと同じものになる。

このように振幅固定で位相をすこしずつ左回転させた空間分布を持つ波動関数は、時間発展でも振幅が変わらず、位相だけが少しずつ右回転していくことになる。

つまり時刻t₀における空間成分さえ与えられれば、その先波動関数がどのように変わっていくかが決まっていることになる。

このFig1.1-3での時間発展の仕組みを数値計算によりシミュレーションしてみよう。ψ(x,t)の成分を時間 t=0で、ｘ上120点サンプリングする。以降Δtごとに、全120点にEq1.1-3の計算を繰り返し適用した。以下をがその結果となる。アニメーションでは初期サンプルを表示後、Δt毎の計算結果を8回分繰り返し表示ている。

Fig1.1-4

今ある波の形状そのものが、波を進める原動力になっていると言えよう。この例で光子が波として表現された場合に、どのような仕組みで進行波となるのかが理解できたと思う。

一方、電子の場合は

\[\small \begin{cases} i\hbar\partial_t\psi(x,t) =E\\ \partial^2_x\psi(x,t)=-k^2=-\frac{\large 2m}{\large \hbar^2}E \end{cases} \]

Eで等式をむすぶと、

\[\small \hbar\partial_t\psi(x,t)=i\frac{\hbar^2}{2m}\partial^2_x\psi(x,t) \tt \tag{Eq1.1-4} \]

となる。つまり『波動関数の時間的変位は、その点の2階空間微分（曲がり具合）により決まる』と言うことだ。尚、光子の場合と違って虚数 i が掛かっていることに注意したい。時間発展をもたらす位相成分（赤いポインター）は、空間二階微分により現在の位相（青いポインター）に対し180度位相が反転し、さらに”i” により90度右回転させられて、光子の場合と同様に元のポインターに対し左90度の位置に直交することになる。

1/(2m)=1、 ℏ=1 とし、これを差分方程式に置き換えると

\[\small\begin{eqnarray} \partial_t\psi(x,t)&=&i\partial^2_x\psi(x,t)\\\ \downarrow\\\small (\psi[n_x,n_t+1]-\psi&[&n_x,n_t])= \\\small i\space \Delta t \space(\psi[n_x+1,n_t]&+&\psi[n_x-1,n_t] -2\psi[n_x,n_t])/\Delta^2 x \\\small \downarrow \\\small \psi’_t [n_x,n_t]&=&i\space \psi_x” [n_x,n_t]\Delta t/\Delta^2 x \end{eqnarray} \]

\[\tt \tag{Eq1.1-5}\]

以下のように光子の場合と同じ図を使ってこのEq1.1-7がどのように時間的発展で機能するのかを説明しよう。

Fig1.1-5

図中の（1）から（3）までの詳細は、

(1)の左向緑ベクトルは左右隣りの複素円における成分との差分から計算されるψ’’_x を表している。円運動なので加速度は円の中心に向かい、当複素円の青ポインターの成分と逆向きになっている。

Fig1.1-6

(2)はEq1.1-7にっ基づき、(1)のψ”_x(左向緑ベクトル)に i Δt/Δx を掛けて算出された時間差分を表している。この量は下段の下向緑ベクトルに対応し、図の通りΔt後の位相を右に回転させる。
(3)の赤ポインターは常に青ポインタ―（波動成分）と共に表示される時間差分の成分であり、下段に記載した下向緑ベクトルと同じものになる。

光子の場合と同様に Eq1.1-7 を使ってt=0からの時間発展をシミュレーションしてみよう。以下の通り、光子とは異なる仕組みで位相は回転させられるが、電子の進行波を再現できている。

Fig1.1-7

以上の様に光子は位相回転速度を原動力とし時間発展し、電子は回転による加速度（遠心力の基）を原動力として時間発展をするというイメージがつかめたと思う。いずれの場合も、その時間発展は同じ振幅を保って位相を回して、進行波となっている。この様な波を、ここでは平面進行波とよぼう。

それなら、『平面進行波の一部を切り取って持ってきても、そいつはかってに動いてくれるのでは？』と期待したくなる。

切り取り方にも色々考えられるが、まずは最も安易な方法で切り取ってみよう。電子の位置x₀=0とすると、ｘ＝0を中心にもつ高さ１の方形波U(x)を波動関数に掛けること、つまりU(x)ψ(x)とすることだ。

Fig1.1-8a

この例では中央 x=0 で幅±2としたが、以下の図のように波動関数を切り取ることができた。

Fig1.1-8b

この波動関数をψ_uとしよう。これを周波数領にフーリエ変換すると、運動量(の1/ℏ倍）を示す波動関数となり、波形は以下の様になる。

Fig1.1-8c

この様にｋは切り取る前には幅の無いδ関数から、大きく広がりを持ってしまい、定数で表せな無くなってしまった。運動量ｋが変わると、エネルギーＥが変わり、それに応じて周波数νが変わってしまう。従って以下の様にωは運動量ｋの関数とる。

従って時間発展を表現しようとすると異なるωに対してそれぞれの波動関数を割り振って合成する必要があることが分かる。しかし、このユニット関数のだと運動量の広がりがあまりにも大きく解析が複雑で自由電子の物理モデルとしてはいまいちだ。方形波をフーリエ変換した波形で分かるように、方形波自身が周波数領域で大きな広がりを持っている。

Fig1.1-9a

位置も運動量もうまく局在化させられる切り取り方を見つけたいところだ。方形波の代わりに以下のガウス関数を使ってみよう。

\[\small g(x)=e^{-(\frac{x-x_0}{2\sigma})^2} \]

Fig1.1-9b

ガウス関数はフーリエ変換後もガウス関数を保ち、方形波と異なり分散が小さい。

Fig1.1-9c

方形波 u(x) を g(x) を置き換えると以下の通りとなる。

Fig1.1-9d

Fig1.1-9e

この波動モデルの数式を完成させる前に、このモデルがどのように動いてくれそうかを差分方程式によるシミュレーションにより明らかにしてみよう。初期値として以下のt=0の波動関数が示す空間分布を与え、Δt後の時間的発展を計算する。

\[\small \Psi(x,0)=e^{-(\frac{x-x_0}{2\sigma})^2} e^{i k_0 x} e^{t=0} \tt \tag{Eq1.1-6} \]

光子波束の時間発展シミュレーション（ -ψ’_x[n_x,n_t] のたしこみループ)の結果は以下のとおり。

Fig1.1-10

期待通り、切り取った波動関数でも光子の波束が形を変えずに移動していくことが分かる。尚、誤差の累積に関する配慮を全くしていないのでこの程度の短い時間範囲で打ち切っている。さて波数一定・振一定の単純な波動関数が進行波となる仕組みは理解できたのが、ガウス関数での切り取りにより波形がかなり複雑になった。フーリエ変換で確認した通り色々な波数が混ざり合ってこの波束を構成いている。この様な状態でどのような力が働き波束を進めているのだろうか？

以下、Fig1.1-10での波の盛り上がった近辺を、適当な間隔に間引いて複素円にプロットした。波束のピークは中央付近の複素円となる。

Fig1.1-11a　波束の中央部

ここで注目すべきは、ピークの左側では時間差分ポインタ（赤）と波動成分ポインタ（青）の交差角が90度より少し開いて、逆に右側では90度より少し閉じているという事だ。以下は右側のさらに先を拡大表示したものだが、しっかり狭まって行くことが確認できる。

Fig1.1-11b　波束の右端（進行前縁）

こちらは左側で90度より開いている。

Fig1.1-11c　波束の左端（進行後縁）

Fig1.1-11c　波束の右端（進行後方向）は振幅を減少させる方向に働き、右側では逆にを振幅を増大させるように働くという事だ。この仕組みが波束のピークを右にに進めているということだ。

同様に電子波束の場合もシミュレーション（ i ψ”_x[n_x,n_t] のたしこみループ)をしてみよう。

Fig1.1-12

光子と同様に波束は移動している。しかし光子とは異なり、移動しながら波束の幅が広がっていく様子がうかがえる。

Fig1.1-13a　波束の中央部

右端よりさらに先を拡大表示してみる。右に進めば進むほど光子以上に時間差分ポインタと波動成分ポインタとの交差角は狭くなっていく傾向が確認できる。

Fig1.1-13b　波束の右端

左も同様に光子以上に開き方が大きくなっている。

Fig1.1-13c　波束の左端

以上により光子も電子も波束を移動させるときに何が起こっているかある程度理解できたと思う。

続いて、なぜ光子波束が形を変えずに進行出来て、電子波束は広がりながら進行するのか、そのメカニズムを考察してみよう。

波動関数ψをガウス関数を使って切り取った。切り取る前の波動関数ψに対して時間差分の操作 -ψ’_x[n_x,n_t] を行うと、確認した通り波は右に進む。ところで切り取りに使ったガウス関数に対しても時間差分の操作 -ψ’_x[n_x,n_t] を行うと、何がおきるか？

以下の図の通り、緑の破線は波動関数を微分しマイナス符号を掛けた、足しこまれるべき時間変位微小量を表示している。

Fig1.1-14a

これを元のガウス関数に足しこむと、右側面が盛り上がり左側面が押し下げられる（右半分は右傾斜なので足しこまれる量は実数で正符号だ。左半分は負符である）ので、右にずれる。従って波は（位相回転しないで）実数のまま以下の通り少しずつ右に移動する。

Fig1.1-14b

切り取った波束は、このガウス関数と、元の平面波である波動関数との合成関数と言えるので、同様に形を変えずに右に進むと考えれば納得できるだろう。

一方電子波束の場合、ガウス関数に時間差分の操作 +i ψ”_x[n_x,n_t] （波動関数を2階微分しi符号で位相を90度回した微小量を足しこむ）を行うと以下の様にガウス関数の形状が持つ曲率の正負の変化（凹凸）に応じて微小加算量の虚数成分（紫の線）がでこぼこと正負の値を持つ。

Fig1.1-15a

例えば右半分だけでも、外寄りの右側は[+]に働き、内寄りの左側は[-]に働く。この点が右半分は単純に[+]だけだった光子の場合と異なる。

この微小加算量が以下のアニメーションの通り、実数内に固定されていたガウス関数を、虚数方向にふらふらと揺さぶる。

Fig1.1-15b

上記以上の時間経過は誤差が大きくなるので、この後導出する時間発展の関数を先出して、プロットしよう。

以下Fig1.1-15cは上記シミュレーションの直後の状態であり、さらにその下のFig1.1-15dは5倍ほど時間が経過したものだ。

Fig1.1-15c

山の両斜面は自身の凹凸形状により端に行くほど位相の回転力が与えられ、波動が広がって行く様子が分かる。まるで山と積み上げた石がごろごろと崩れ落ちるように崩壊していくようだ。

Fig1.1-15d

尚、グレーの線は振幅（包絡線と同じ）をプロットしており、ひしゃげてしまっているが、ガウス関数の形は保っている。

このガウス関数と平面進行波の合成関数である電子波束は、上記の様に内在するガウス関数の崩壊によって、広がりながら移動して行くと考えると、納得できるだろう。

1.2 解析的アプローチ

シミュレーションでは限界があるので波束の時間的発展を含めた波動関数を求めてみよう。結論から言うと、（少々複雑にデコレーションされた感があるが）以下の式が答えと言う事の様だ。＊参照URL（田内修＠筑波大）

\[\small \psi (x,t)=\sqrt{\frac{1}{\sqrt{2\pi}\sigma_{x0}(1+i\xi t)}}\exp\left[\frac{-x^2/4\sigma_{x0}^2+i(k_0x-\omega_{k0} t)}{1+i\xi t}\right], \]

\[\small \xi=\frac{\hbar}{2m\sigma_{x0}^2 },\space\omega_{k0}=\frac{\hbar k_0^2}{2m} \tt \tag{Eq1.2-1} \]

シミュレーションではEq1.1-6を t=0 として利用したのだが、この式では時間的発展の表現を組み込んであるので任意の時間を指定して空間分布を知ることが可能となる。

　自然単位系（ｍ＝1、ℏ＝1）とし、適当な運動量(p=ℏkより、k₀に値を代入)を与えて、以下の様に、時間ｔを０から増やしながら順次プロットすると、確かにx₀付近にいた波動関数が時間経緯と共に右に移動していく様子が分かる。シミュレーションの限界だったFig1.1-12の続きを見ることになる。

時刻１（t＝0で、位置ｘ＝0 に局在）

時刻 2（右に移動しながら、幅が広がり始める）

時刻 3（さらに右に移動しながら、広がりが大きくなり高さは低くなる）

Fig1.2-1

包絡線のピークが電子のニュートン力学的な位置: x を表し x=vt に従う。もちろん v=p/m となる。尚、波束ピークの移動速度は群速度と呼ばれる。

波動関数の広がりから分かる通り、何と電子は時間が経つと、どこにいるのか分からなくなってしまうという事だ。電子波を無理やり局在化させて粒に見立てても、時間がたつとどんどん広がって波に帰ってしまうと言うことなのか？

続いてEq1.2-1を導出してみよう。流れとしてはシミュレーションと同様に、以下の様に時刻ｔ=0を始点とする。

\[\small \psi_x(x,0)=g_x(x) e^{i k x-0 t} =e^{-\large(\frac{ x}{2\sigma})^2} e^{i k x} \]

\[\small \tt \tag{Eq1.2-2} \]

この式は運動量p=ℏkを持つ平面進行波の波導関数にガウス関数を掛けて無理やり局在化させたものであった。簡略のため、σ＝１とし、任意の運動量k₀を与えると

\[\small \psi_x(x,0) =e^{-\large(\frac{ x}{2})^2} e^{i k_0 x} \]

以下は k₀=3の運動量を持つ場合の例となる。包絡線はガウス関数だが、複素数成分がねじれて運動量をため込んだイメージだろうか。

Fig1.1-3 g_x(x)e^ik₀x

因みに、静止している状態はk₀=0でありFig1.1-14aと同じ形状を持つガウス関数となり、その成分は実数のままで、複素円は回転しない。

この『ねじをまかれた(？)ガウス関数波束』を波数ｋでフーリエ変換してみると、

\[\small \begin{eqnarray} \psi(x,t)&=&\frac{1}{\sqrt{2\pi}}\int^\infty_{-\infty} [g_x(x) e^{ik_0x}]e^{-ikx}dx\\\small &=&\frac{1}{\sqrt{2\pi}}\int^\infty_{-\infty}[e^{-\large(\frac{ x}{2})^2}e^{ik_0x}]e^{-ikx}dx\\\small &= &\sqrt{2}e^{-(k-k_0)^2}\equiv g_k(k) \end{eqnarray} \]

以下の様に k=3にピークを持つガウス関数g_k(k)となる。

Fig1.1-4　g_k(k)

k₀の単一周波数しか持たなかった平面波は、切取りという操作により、様々な波数を持つ波の足し合わせで表される様になってしまった。このk₀＝３を中心に広がったガウス関数は、それぞれの波に振幅を表している。

具体例としてピーク周辺の5点ｋ=｛2.0, 2.5, 3.0, 3.5, 4.0｝を選んで{g_k(2.0), g_k(2.5), g_k(3.0), g_k(3.5), g_k(4.0)｝の様に振幅係数のサンプルをとり、それぞれのg_k(k)e^ikx をプロットすると、以下の様に振幅の異なる平面波が現れる。もちろんｋ＝３の平面波の振幅が最大となる。

Fig1.1-5

この5つのサンプルを足し合わせたものが一番下のグラフだ。たった5つのサンプルをとっただけだが、結構きれいにFig1.1-3 の波束を再現できていると思う。

この様に運動量空間のガウス関数は、運動量(空間周波数)の異なる波が、どれくらい含まれているかを示しており、このサンプルを細かく採って足し合わせる（つまり積分をとる）とFig1.1-3が再現できることになる。

以上の様に周波数領域にフーリエ変換された波動関数はp=ℏkの関係により、運動量空間で定義される波動関数と言える。波動関数を示す数式も以下のとおり、どちらの定義でも使えるようになっている。

\[\small \psi_x(x,t)\equiv A e^{i k x-\omega t}\equiv \psi_k(k,t) \]

両者は等価の情報を持ち、同じ量子状態を空間領域で表示するか、運動量領域で表示するかという『表示の仕方』が違うだけだと考えよう。

ここまで来ると、あとはFig1.1-4に時間の項を追加するだけだ。以下の様に運動量空間の波動関数においてt を復活させて、エネルギーと運動量の関係から（ℏを含めすべての係数を１としているので）ω=k²と置けば良い。

\[\small \psi_k(k, t)=g_k(k) e^{i k x-i\omega t} =g_k(k) e^{i k x-i k^2 t} \]

\[\small \tt \tag{Eq1.2-3} \]

結果は以下の通りで、運動量の低い（例えばk=2）はゆっくり進行し、運動量の高い（例えばk=4）では早く進行している。

Fig1.1-6

K=2 と k=4 では ω=k² の関係より4倍速いと思うかもしれないが、ωが2倍によりλが1/2となっているのでv=ωλにより、結局2倍の速さとなっている。

6番目のグラフはこれ等5個の波を以下のように足し合わせた結果の波束となり、波束の方も移動することを確認できる。

\[\small \begin{eqnarray} \psi_{c5}(x,t)&=& \frac{1}{5}\sum_{\scriptsize k=(2, \space 2.5,\space 3,\space 3.5,\space 4)}g_k(k) e^{ikx} \\\small &=& \frac{1}{5}\sum_{\scriptsize k=(2, \space 2.5,\space 3,\space 3.5,\space 4)}\sqrt{2}e^{-(k-k_0)^2-i k^2 t} e^{ikx} \end{eqnarray} \]

\[\small \tt \tag{Eq1.2-4} \]

前述の通り、ピークの進行波（k=3）の速度に対して、k=2は遅く、k=4は早い。そのため波束の後縁では、波長が長く進行速度の遅いの波となり、そして前縁では、波長が短く進行速度の速い波となる。

この様に波束は移動しながらも両端が広がって行くことが理解できたのでは思う。

さて、この6番目の移動する波束のグラフを、『数式』で記述することが目的であった。そのためには以下の通り、Eq1.2-4の足し合わせの代わりに以下のように積分すればよい。

\[\small \begin{eqnarray} \psi(x,t)&=&\frac{1}{\sqrt{2\pi}}\int^\infty_{-\infty} g_k(k) e^{ikx}dk \\\small &=&\frac{1}{\sqrt{2\pi}}\int^\infty_{-\infty}\sqrt{2}e^{-(k-k_0)^2-i k^2 t} e^{ikx}dk \end{eqnarray} \]

\[\small \tt \tag{Eq1.2-5} \]

これはフーリエ逆変換の式と同じものになっていることが分かると思う。積分の頭にある係数はフーリエ変換での慣例的なものを使った。

Fig1.1-6の正確な時間発展の様子を見たいだけなら、この積分を数式処理に任せて運動量、つまりkを指定し、時間ｔを刻みながら、xを振って空間分布のグラフを複数回プロットすればよい。しかし数式処理（Mathematica ）の吐き出した積分結果はかなりごちゃごちゃして理解に苦しむので、もう少し整理を進めよう。

まず、積分式

\[\small \frac{1}{\sqrt{\pi}}\int^\infty_{-\infty}e^{-(k-k_0)^2-i k^2 t} e^{ikx}dk \]

の指数部を抜き出して以下の様に[ak+b]²の項にまとまるように整理する。

\[\small \begin{eqnarray} & &-(k-k_0)^2-i k^2 t+ikx\\\small &=&-k^2+2k k_0-k_0^2-ik^2 t+ikx\\\small &=&-(1+it)k^2+(ix+2k_0)k-k_0^2 \\\small&=&-\left[(\sqrt{1+it}\space k)^2-\frac{\sqrt{1+it}(ix+2k_0)k}{\sqrt{1+it}}+\frac{(ix+2k_0)^2}{(2\sqrt{1+it})^2}\right]\\\small&&+\frac{(ix+2k_0)^2}{(2\sqrt{1+it})^2}-k_0^2\\\small&=&-\left[\sqrt{1+it}\space k-\frac{ix+2k_0}{2\sqrt{1+it}}\right]^2+\frac{-x^2/4+i( k_0 x – k_0^2) t}{1+it} \end{eqnarray}\\\small \]

ガウス積分の形態として以下がある。

\[\small \int^\infty_{-\infty}e^{-(a k+b)^2}\space dk= \sqrt{\frac{\pi }{a^2}},\space Re(a^2)\geq 0 \]

これにより\(\small a\rightarrow\sqrt{1+i t} \) と置き換えると

\[\small \begin{eqnarray} &&\frac{1}{\sqrt{\pi}}exp\left[\frac{-x^2/4+i( k_0 x – k_0^2) t}{1+it}\right]\\\small &\times&\int^\infty_{-\infty}exp\left[-\left(\sqrt{1+it}\space k-\frac{ix+2k_0}{2\sqrt{1+it}}\right)^2\right]dk \end{eqnarray}\\\ \downarrow\\\small \frac{1}{\sqrt{\pi}}exp\left[\frac{-x^2/4+i( k_0 x – k_0^2) t}{1+it}\right]\sqrt{\frac{\pi}{1+it}} \]

となり以下の時間発展解を得る。

\[\small \frac{1}{\sqrt{\pi}}\int^\infty_{-\infty}e^{-(k-k_0)^2-i k^2 t} e^{i k x}dk\\\small = \frac{1}{\sqrt{\pi(1+i\space t)}}exp\left[\frac{-x^2/4+i \left(k_0 x-k_0^2 t\right)}{1+i \space t}\right] \]

ｔに掛かったk²がEq1.2-1におけるωに対応する。なお、Eq1.2-1と異なり下記の様に簡略化してあるので振幅等異なってなっている。

\[ \scriptsize(Eq1.2-1):\small\frac{1}{\sqrt{\sqrt{2\pi}\sigma_{x0}(1+i\xi t)}}\exp\left[\frac{-x^2/4\sigma_{x0}^2+i(k_0x-\omega_{k0} t)}{1+i\xi t}\right], \]

\[\small \xi=\frac{\hbar}{2m\sigma_{x0}^2 }\rightarrow1,\space\omega_{k0}=\frac{\hbar k_0^2}{2m}\rightarrow k_0^2,\space\sigma_{x0}\rightarrow 1,\space \sqrt{\sqrt{2\pi }}\rightarrow \sqrt{\pi} \]

分散等を指定したい場合は、オリジナルの数式Eq1.2-2を使用のこと。

1.3 まとめ

自由空間における電子と光子の運動が量子力学ではどのように記述されるか、主に波動関数の数式の変形を追いかけた。直観的な理解を得られるように数値計算によるシミュレーションで補足した。結局、（正直、わけが分からなくなったのだが）得られたイメージととしては

光子波束は位相の回転速度を原動力としてどこまでも転がって行く。
電子波束は転がるために加速度を利用しているが、悲しいかな質量を持つゆえに広がって消えてしまう。（『おーい、電子や。おまえどこにいる？』,　電子：『どこだか分からなくなりました。ただ運動量 k₀ 前後で走っている事だけは分かってます』)
質量の有無で波動の伝搬する仕組みが変わっていが、波動関数は空間と時間を覆う連続体の様な実在だ。そうすると、それは物理的な基盤であり時空の構成要素なのかも？
電子や光子を『ここに在る物理実態』と認識するためには局在化させる操作が必要となる。つまり全時空に広がった波動関数という物理基盤を時間空間を定めて『摘まむ(？)とか、つっつく(？)』ようなイメージだが、これが量子物理学で言う所の『測定』という行為になるのだろう。

以上は私の感想であり、あなたがどう解釈するかはあなた次第だと思う。量子を記述する数式は合っていると思うが、数式化されたからと言って、量子がなんであるかが理解できたかは別のようだ。

相対性理論に比べて量子力学や波動関数などは、原子など限られたミクロの現象を扱うときに利用される、なんとも辛気臭い学問というイメージがあった。しかし以上の様に『なんと波動関数は宇宙の時空を構成する壮大な物理基盤なのか？』という風に認識に変わった。

アインシュタインは相対性理論と同様に時空の幾何学をもって統一理論を構築しようとしたがうまくいかなかったとのこと。しかし多くの量子物理学者が量子論（波動力学や、ふれなかったが等価なマトリックス力学）を元に統一理論の構築に進んでいった。おぼろげながらも、その動機が理解できた様な気がする。

続いては、統一理論への第一歩となる量子力学の相対論的拡張となる。

つづく→　場の量子論への準備２（相対論的量子論）

Post Views: 25