量子力学Ⅰ/波動関数の解釈/メモ のバックアップの現在との差分(No.2)

更新


  • 追加された行はこの色です。
  • 削除された行はこの色です。
* 目次 [#obe40a19]

#contents

* 概要 [#d7d6b3b3]

[[量子力学I/波動関数の解釈]] に関するメモです。
[[量子力学Ⅰ/波動関数の解釈]] に関するメモです。

* 測定される物理現象 [#q0464429]
** Mathematica ソース [#if68c58e]

上記干渉のアニメーションを表示する。
干渉のアニメーションを表示する。

 LANG:mathematica
 anim = Table[
    Plot3D[
      Re[ Exp[I(Sqrt[(x - d)^2 + y^2] - t)]/((x - d)^2 + y^2) + 
          Exp[I(Sqrt[(x + d)^2 + y^2] - t)]/((x + d)^2 + y^2) ] /. d -> 20.25,
     {x, -30, 30}, {y, 0, 60}, 
     Mesh -> None, Mesh -> Automatic, ViewPoint -> 1.4 {1.2, 1, 1}, 
     PlotPoints -> 401, ImageSize -> Large, 
     PlotRange -> {{-30, 30}, {0, 60}, {-0.01, 0.01}}
   ], {t, 0, 2 Pi, 2 Pi/40}
 ];
 Export["double-slit.gif", anim, "GIF"]

2015-02-23 時点において、この計算にはかなりの時間(10分単位)がかかる。

* 波動関数から各種物性値を取り出すには [#z2a7962b]
** 解答:エネルギーの場合  [#y53b9710]
* 確率密度の時間変化 [#i58843b4]

(1)
** 全確率の保存 [#n19bd59e]

 &math(
\overline E &=\iiint \psi_k^*(\bm r,t)\hat H\psi_k(\bm r,t)d\bm r\\
&=\iiint \psi_k^*(\bm r,t)E_k\psi_k(\bm r,t)d\bm r\\
&=E_k\iiint |\psi_k(\bm r,t)|^2d\bm r\\
&=E_k
);
*** 線形代数を用いた証明 [#f6db1a8f]

(2)
シュレーディンガー方程式を変形すれば、

 &math(
\sigma_E^2&=\iiint \psi_k^*(\bm r,t)(\hat H-\overline E)^2\psi_k(\bm r,t)d\bm r\\
&=\iiint \psi_k^*\Big(\bm r,t)(\hat H^2\psi_k(\bm r,t)-2\hat H\overline E\psi_k(\bm r,t)+\overline E^2\psi_k(\bm r,t)\Big)d\bm r\\
&=\iiint \psi_k^*\Big(\bm r,t)(E_k^2\psi_k(\bm r,t)-2E_k\overline E\psi_k(\bm r,t)+\overline E^2\psi_k(\bm r,t)\Big)d\bm r\\
&=(E_k-\overline E)^2\iiint |\psi_k(\bm r,t)|^2d\bm r\\
&=(E_k-\overline E)^2
);
 &math(\frac{\PD}{\PD t}\psi(\bm r,t)=(-i\hat H/\hbar)\psi(\bm r,t));

上で見たように &math(\overline E=E_k); であるから、&math(\sigma_E^2=0); すなわち 
&math(\sigma_E=0); となる。
と書ける。したがって、&math(\Delta t); が小さく &math(t); から &math(t+\Delta t); の間に 
&math(\hat H(\hat{\bm r},\hat{\bm p},t)); が変化しないと見なせるとき、

 &math(\psi(\bm r,t+\Delta t)=e^{-i\Delta t\hat H(\hat{\bm r},\hat{\bm p},t)/\hbar}\psi(\bm r,t));

と表せることになる。

&math(\hat H); がエルミートであれば &math(-\Delta t\hat H/\hbar); もエルミートであるから、

 &math(\hat U=e^{-i\Delta t\hat H(\hat{\bm r},\hat{\bm p},t)/\hbar});

の部分はユニタリーになる。→ [[[&math(H); がエルミート行列の時、&math(e^{iH}); がユニタリーになることの証明]>量子力学Ⅰ/線形代数の復習#a5808b18]]

ユニタリー変換の性質より、

 &math(\|\psi(\bm r,t+\Delta t)\|^2=\|\hat U\psi(\bm r,t)\|^2=\|\psi(\bm r,t)\|^2); 

である。

この等式が任意の時刻に成り立つことから、シュレーディンガー方程式に従った波動関数の時間発展では、
確率密度が保存される((時間発展演算子はユニタリーである、とも言える))。


Counter: 3252 (from 2010/06/03), today: 1, yesterday: 0