たくのろじぃのメモ部屋

プログラミング・数学・マイクラなどの内容を備忘録として残すブログ。

【量子力学メモ】シュレディンガー方程式 2

前回は空間依存のシュレディンガー方程式を導出しました。

2. 時間に依存するシュレディンガー方程式

まずは波動関数の式です。

{\displaystyle 
\Psi(x, t) = Ae^{i(kx - i \omega t)}
}

前回は空間依存のため、 x偏微分しましたが、今回は時間依存を考えたいので  t偏微分します。

{\displaystyle 
\frac{\partial }{\partial t} \Psi(x, t) = \frac{\partial }{\partial t} Ae^{i(kx - i \omega t)}
}

これを計算すると、

{\displaystyle 
\frac{\partial }{\partial t} \Psi(x, t) = -i \omega Ae^{i(kx - i \omega t)}
}

右辺の  -i \omega にかかっているのは  Ae^{i(kx - i \omega t)} で、波動関数  \Psi(x, t) そのものです。よって、

{\displaystyle 
\frac{\partial }{\partial t} \Psi(x, t) = -i \omega \Psi(x, t)
}

となります。つまり、複素数における角速度  -i\omega を作用させたものが波動関数の時間変化になります。

さて、この式をよく見ると角速度がかかっています。前回同様、光子のエネルギーと運動量を考慮します。

{\displaystyle 
光子のエネルギー : E = \hbar \omega  光子の運動量 : p = \hbar k
}

今回は  \omega が含まれているので、両辺に  \hbar をかけると

{\displaystyle 
\hbar \frac{\partial }{\partial t} \Psi(x, t) = -i \hbar \omega \Psi(x, t)
}

となり、 \hbar \omega E が代入できるので、

{\displaystyle 
\hbar \frac{\partial }{\partial t} \Psi(x, t) = -i E \Psi(x, t)
}

となります。右辺の  -i が邪魔なので、両辺に  i をかけると

{\displaystyle 
i \hbar \frac{\partial }{\partial t} \Psi(x, t) = E \Psi(x, t)
}

となり、エネルギー  E波動関数  \Psi に作用する式になりました。左辺を見てみると、 i \hbar \frac{\partial }{\partial t} \Psi に作用しているので、 i \hbar \frac{\partial }{\partial t} はエネルギー演算子となります。(前回のように分数形式にもできますが、分数形式は扱いにくいので省略します。)

なんと、今回はこれだけで時間依存のシュレディンガー方程式を導出できます。エネルギーは運動エネルギーとポテンシャルエネルギーの和  E = T + V ですが、 T は既に前回、 - \frac{\hbar^{2}}{2m} \frac{\partial^{2}}{\partial x^{2}} と求められたので代入すると

{\displaystyle 
i \hbar \frac{\partial }{\partial t} \Psi(x, t) = - \frac{\hbar^{2}}{2m} \frac{\partial^{2}}{\partial x^{2}} \Psi(x, t) + V(x) \Psi(t, x)
}

となり、波動関数  \Psi(x, t) でくくれば

{\displaystyle 
i \hbar \frac{\partial }{\partial t} \Psi(x, t) = \left[ - \frac{\hbar^{2}}{2m} \frac{\partial^{2}}{\partial x^{2}} + V(x) \right] \Psi(x, t)
}

となります。

これで時間に依存するシュレディンガー方程式を求めることができました。

3. ハミルトニアン

さて、シュレディンガー方程式を導出してきましたが、時間に依存するシュレディンガー方程式を見ると、左辺は演算子で右辺はエネルギーが波動関数に作用しているので、両辺を波動関数 \Psi(x, t) で割っても正しい解は得られません。

なぜなら演算子単体では意味を持たないためです。数値あるいは関数に作用させて初めて意味を持つことになります。例えば、 + 単体では何も意味がないですが、 a + b とすれば意味を持つことになり、「 + a b を足し合わせる(の線形和である)」ことを表現します。これと同じように、右辺に対しても何かしらの演算子を定義して、波動関数の値を求められるようにしていきます。

と、行きたいところですが丁寧にやりすぎると大変なので少し省きます。内容としては解析力学です。

https://blog.takunology.jp/entry/2020/05/24/164156blog.takunology.jp

まず、ラグランジアン L は運動エネルギーとポテンシャルエネルギーの差  L = T - V です。(文字は量子力学に合わせます。)ここで、運動している空間を線形  x軸方向だけと仮定すれば、速さ v は距離  x微分なので  \dot{x} となります。つまり、ラグランジアン

{\displaystyle 
L = \frac{1}{2} m\dot{x}^{2} - V(x)
}

となります。ラグランジアン微分(運動エネルギーを微分)すれば、

{\displaystyle 
\frac{\partial}{\partial \dot{x}} L = \frac{\partial}{\partial \dot{x}} \left( \frac{1}{2} m\dot{x}^{2} - V(x) \right)
}

{\displaystyle 
p_{x} = m\dot{x}
}

となり、運動量となりました。ここで共役運動量を求めると、

{\displaystyle 
\dot{x} = \frac{p_{x}}{m}
}

となります。共役運動量は共役複素数と似たようなものです。運動量の式に共役運動量を代入してみれば  m が約分されて  p_{x} となるのと同じです。本質的な説明にはなっていませんが、ここでは省略します。

ラグランジアンは一般化座標と一般化速度であったのに対し、ハミルトニアンは一般化座標と一般化運動量によって力学を表現することになります。ハミルトニアンの定義は

{\displaystyle 
H = \sum_{i} p_{i} \dot{q_{i}} - L
}

です。 i虚数単位ではないのでご注意を。 p は運動量、 \dot{q_{i}} は一般化座標の時間微分です。今回は座標を x軸のみとっているので、距離 x微分、すなわち速度  \dot{x} と同じです。なので、ラグランジアンを代入すると

{\displaystyle 
H = p_{x} \dot{x} - \left[ \frac{1}{2} m\dot{x}^{2} - V(x) \right]
}

ここで、共役運動量を  \dot{x} に代入すると、

{\displaystyle 
H = p_{x} \frac{p_{x}}{m} - \left[ \frac{1}{2} m \left( \frac{p_{x}}{m} \right)^{2} - V(x) \right]
}

右辺の運動量をまとめ、2乗を展開して  m で約分すると

{\displaystyle 
H = \frac{p_{x}^{2}}{m} - \left[ \frac{p_{x}^2}{2m} - V(x) \right]
}

右辺の括弧を展開すると

{\displaystyle 
H = \frac{p_{x}^{2}}{m} - \frac{p_{x}^2}{2m} + V(x)
}

分母が  m 同士で計算できるので、 \frac{2}{2m} - \frac{1}{2m} と考えれば

{\displaystyle 
H = \frac{p_{x}^{2}}{2m} + V(x)
}

となります。これがハミルトニアンです。

この式、見覚えありませんか?
シュレディンガー方程式の右辺と似たような形をしています。

{\displaystyle 
 \left[ - \frac{\hbar^{2}}{2m} \frac{\partial^{2}}{\partial x^{2}} + V(x) \right] \Psi(x, t)
}

そこで、ハミルトニアンにちょっとしたアプローチを行います。運動量が古典力学のままなので、量子力学に変換します。光子の運動量は、 p = \hbar k でしたので、

{\displaystyle 
H = \frac{\hbar^{2} k^{2}}{2m} + V(x)
}

となります。しかし、これでは波数 k が残っており、二階偏微分もないです。ですが、空間依存するシュレディンガー方程式を解くときに、 k^{2} にするために二階偏微分し、そこに運動量と運動エネルギーの関係式を波数  k について解いたものを代入しました。よって、シュレディンガー方程式によって求めた運動量演算子

{\displaystyle 
p = - i \hbar \frac{\partial }{\partial x}
}

を代入して計算すると、

{\displaystyle 
H = \frac{1}{2m} \left( - i \hbar \frac{\partial }{\partial x} \right)^{2} + V(x)
}

となり、計算すると (-i)^{2} でマイナスが残り、2階偏微分になるので

{\displaystyle 
H = - \frac{\hbar^{2}}{2m} \frac{\partial^{2}}{\partial x^{2}} + V(x)
}

となります。よって、シュレディンガー方程式の右辺をハミルトニアンで表現することができ、最終的には

{\displaystyle 
E \Psi(x, t) = H \Psi(x, t)
}

という形にまとめることができます。これがハミルトニアンを用いたシュレディンガー方程式です。

参考文献

高校数学で分かるシュレディンガー方程式

よくわかる解析力学