量子力学Ⅰ/波動関数の解釈 のバックアップ差分(No.7)

更新


  • 追加された行はこの色です。
  • 削除された行はこの色です。
[[量子力学I]]

#contents

* 波動関数の物理的な意味は何か? [#e131650a]

シュレーディンガー方程式が正しいとして、そこから求まる波動関数の物理的意味は何であろうか?

** 話の進み方が奇妙? [#fd91a38a]

近代の物理学ではこのように、「正しそうな方程式が先に求まって、後からその意味を考える」
という手順を踏むことがしばしばある。

「直感的に理解できる」方程式が基礎となっていた古典論と比べると、
このようなやりかたは始め奇妙に思えるが、
量子力学の真価が現れるような極小の世界に対しては我々の「直感の働かない」のだからしかたがない、
量子力学の真価が現れるような極小の世界に対しては「我々の直感の働かない」のだからしかたがない、
と考えて納得してほしい。
実際これから直感とはかけ離れた物理現象を学ぶことになる。

+ 観測された物理現象を説明できそうで、かつ論理的に矛盾のない基礎方程式を適当にでっち上げる
近代科学ではむしろ「直感」をなるべく廃した、以下のような手法が科学的手法と考えられている。

+ 観測された物理現象を説明できそうで、かつ%%%論理的に矛盾のない%%%基礎方程式をでっち上げる
+ その方程式から何が予測されるかを考えて、新たな測定結果と突き合わせる
+ 矛盾が生じれば理論に修正が必要となる
+ 矛盾無く測定結果を説明できている限り、それが正しい理論である

このような意味で現在受け入れられているのが、
先に %%求めた%% でっち上げた シュレーディンガー方程式である。

* 測定される物理現象 [#j6d4afe8]

電子の波動関数の意味を理解する上で最も分かりやすい実験結果の1つに、
二重スリットを通る電子についての測定が挙げられる。

二重スリットの実験は光学でもよく知られており、ホイヘンスの原理によれば、
スリット後の状況は2つのスリットの位置にそれぞれ点光源があるのと同じであり、
スリット後の状況は2つのスリットの位置にそれぞれ光源があるのと同じであり、
それらの光の干渉によりスクリーン上に干渉縞が現れる。
当然、片方のスリットを閉じれば干渉縞は消失する。

&attachref(double-slit.png,,25%);      
&attachref(double-slit.gif,,50%);

電子に対しても同様の実験を行うことができて、その結果、光と同様に干渉縞が現れる。

#ref(Double-slit_experiment_results_Tonomura_Preview.jpg,right,around,75%);

電子が古典的な粒子であればスリットを通った電子は単に直進するのみであるから、
もう一方のスリットが開いていようが閉じていようが
1つのスリットを通った電子がスクリーン上に到達する位置は変わらないはずで、
二重スリットの結果は単一スリットの結果の足し算になるはずである。

それにもかかわらず干渉縞が現れることは、確かに電子が波の性質を持っていることを表わしてる。

それでは電子に粒子としての性質が無いのかといえばそんなことはない。
高感度のスクリーンを用いることで、
1つ1つの電子がスクリーン上に到達した位置を記録することができる。
干渉縞が現れる場合にも、1つの電子は1つの点に到達する。
すなわち、電子は確かに粒子なのである。

有名な&ruby(とのむら){外村};による二重スリットの測定結果は右図のようになり、
有名な&ruby(とのむら){外村};による二重スリットの測定結果は右図のようになる。
1つ1つの電子がスクリーン上に到達する位置はランダムに見える物の、
その確率が波動関数の干渉により波打つために、
多数の電子に対して観測を繰り返すことにより干渉縞が現れてくる。(from [[Electron double-slit experiment (dr. Tonomura, Hitachi Research) - ATLAS@CERN>http://www.learningwithatlas-portal.eu/en/node/93624]])

* コペンハーゲン解釈 [#of5d8123]

現在受け入れられている波動関数の解釈は以下のような物である。コペンハーゲンにあるボーア研究所で主に構築されたため、コペンハーゲン解釈と呼ばれる。((歴史的にはこの解釈がなされた時点で上記のような顕著な実験結果は得られておらず、この解釈が受け入れられるまでにはかなりの紆余曲折があった))

- 電子の時間発展はシュレーディンガー方程式を満たす波動関数 &math(\psi(\bm r,t)); で記述される
- 波動関数は空間的に広がりを持ち、また、干渉や回折などの波に特有な性質を現す
- 電子の位置を実験的に観測した場合には電子はある一点に見出される
- 位置 &math(\bm r); に電子が発見される確率は &math(|\psi(\bm r,t)|^2); に比例する

二重スリットの実験に当てはめれば、

- 二重スリットを通る電子の波動関数は、2つのスリットのそれぞれを通る経路の間で干渉を起こす
- その結果、波動関数の絶対値の二乗 &math(|\psi(\bm r,t)|^2); に干渉縞が現れる
- 電子がスクリーンに当たり、その位置が記録されることが「観測」にあたる
- 観測される電子位置は常に1点だけに定まり、広がりを持たない
- &math(|\psi(\bm r,t)|^2); の大きな箇所でより多くの電子が発見されるため、
多くの電子について観測を繰り返すことによりスクリーン上に干渉縞が現れる

空間的に広がりを持つ電子が観測により1点に見出される様子は「波動関数の収束」と呼ばれる。

現在の量子力学は、なぜ観測により波動関数が収束を起こすのか、
とか、観測しないときに電子はどの位置にあるのか、といった問いには答えない。
「観測によって検証できない内容」は物理の範疇ではないというスタンスである。

そのかわり、「観測により確かめられる内容」については量子力学は完璧な予想を与える。
(対応するシュレーディンガー方程式が数学的に解ける範囲にある限り)

#clear

* 確率密度関数について [#q2585cd6]

波動関数の絶対値の二乗 &math(|\psi(\bm r,t)|^2); は時刻 &math(t); において位置 &math(\bm r); 
に電子を見出す確率密度関数となる。

確率・統計の需要を取った学生には学習済みの内容であるが、
確率密度関数についてここで復習しておく。

** 確率密度関数の定義 [#x546ffd8]

#ref(はじめての誤差論/probability_density_function.png,right,around,50%);

「確率変数 &math(x); が確率密度関数 &math(f(x)); に従う」という意味は、

測定毎に異なる値を取る変数 &math(x); があり、1回測定したときに 
&math(x); が &math(x_a<x<x_b); の範囲に入る確率を

 &math(P\{x_a<x<x_b\}=\int_{x_a}^{x_b}f(x)dx);

として求められる、ということである。

右図で分かるとおり、上式の右辺の積分は &math(f(x)); と 
&math(x); 軸に挟まれる部分の面積に相当するため、
&math(f(x)); が大きいことは &math(x); がその付近の値を取りやすいことを表わしており、
&math(f(x)); が小さければ &math(x); はその付近の値を滅多に取らない。

当然、すべての &math(x); に対して &math(f(x)\ge 0); であり、
また、測定を行えば必ず何らかの値が得られることから、&math(f(x)); を全範囲にわたって積分した値は常に1になる。

 &math(\int_{-\infty}^{\infty}f(x)dx=1);

** 期待値 [#hd633e60]

確率変数 &math(x); の測定を仮想的に無限回行った場合に得られる平均値を
&math(x); の期待値という。

例えば、サイコロを1回振った際に出る目の期待値は、

 &math(1\cdot\frac{1}{6}+2\cdot\frac{1}{6}+3\cdot\frac{1}{6}+4\cdot\frac{1}{6}+5\cdot\frac{1}{6}+6\cdot\frac{1}{6}=3.5);

のように、サイコロの目の値に、個々の目の出る確率を掛けて得られる。

すなわち、離散的な値を取る確率変数 &math(x); の期待値 &math(\langle x\rangle); は、

 &math(\langle x\rangle=\sum_{k} x_k p(x_k));

として与えられる。ここで &math(p(x)); は &math(x=x_k); となる確率である。

いま考えている &math(x); が連続値を取る場合にも、これと同様に

 &math(\langle x\rangle=\int_{-\infty}^\infty xf(x)\,dx);

により期待値が求まる。測定値が範囲 &math([x,x+dx]); に入る確率を &math(f(x)dx); と表せることに注意せよ。

** 標準偏差 [#f028b0d1]

標準偏差は確率変数 &math(x); がその期待値からどれほど大きなバラツキを持つかを表わす指標である。

統計学的には、&math(n); 回の測定値 &math(x_1,x_2,\dots,x_n); の標準偏差は、

 &math(\sigma_x=\sqrt{\frac{1}{n}\sum_{k=1}^n (x_k-\langle x\rangle)^2});

として定義される。

これに対して確率変数 &math(x); を仮想的に無限回観測した際の標準偏差の期待値は

 &math(\sigma_x=\sqrt{\int_{-\infty}^\infty (x-\langle x\rangle)^2f(x)\,dx});

として求められる。これを確率密度関数 &math(f(x)); の標準偏差という。

** 任意の関数の期待値 [#zebdd06f]

標準偏差の場合と同様に、確率変数 &math(x); の関数として表わされる任意の値 &math(g(x)); の期待値を、

 &math(\langle g(x)\rangle=\int_{-\infty}^\infty g(x)\,f(x)\,dx);

として求められる。&math(g(x)); が線型である場合を除き、

 &math(\langle g(x)\rangle\ne g(\langle x\rangle));

であることに注意せよ。

** 複合確率 [#wb9dcaf6]

確率変数 &math(x); が区間 &math([x, x+dx]); に入る確率が &math(f(x)dx); であり、~
確率変数 &math(y); が区間 &math([y, y+dy]); に入る確率が &math(g(y)dy); であり、~
両者が同時に起きる確率が &math(h(x,y)dxdy); であるとする。

「&math(x); と &math(y); が独立ならば」&math(h(x,y)=f(x)g(y)); である。

「&math(x); と &math(y); に相関があれば」&math(h(x,y)\ne f(x)g(y)); である。

#clear

* 波動関数から各種物性値を取り出すには [#z663cace]

シュレーディンガー方程式を解いて波動関数 &math(\psi(\bm r,t)); が得られたとする。
この波動関数から電子の運動について何をどのように知ることができるのだろうか?

** 波動関数の規格化 [#i56b385a]

シュレーディンガー方程式は線型な方程式だから、
ある関数 &math(\psi_1(\bm r,t)); が解であれば、
任意の定数 &math(A); に対して &math(A\psi_1(\bm r,t)); も解である。

一方、波動関数の絶対値の二乗が確率密度関数となるためには、

 &math(\iiint|\psi(\bm r,t)|^2\,dx\,dy\,dz=\iiint|\psi(\bm r,t)|^2\,d\bm r=1);

でなければならない。

任意のシュレーディンガー方程式の解 &math(\psi_1(\bm r,t)); に対して、

 &math(A=\sqrt{\iiint|\psi_1(\bm r,t)|^2\,d\bm r});

を掛けることで

 &math(\psi(\bm r,t)=A\psi_1(\bm r,t));

のように波動関数の規格化し、&math(|\psi(\bm r,t)|^2); が確率密度関数になるようにできる。

** 全確率の保存 [#i7b418d2]
後に見るように、
[[時間に依存するシュレーディンガー方程式は全確率を保存する>量子力学Ⅰ/確率密度の保存]]ため、
ある時刻で規格化された波動関数は任意の時刻で規格化された物になる。

ある時刻において規格化した波動関数がシュレーディンガー方程式に従って時間発展する時、
確率密度が増えたり減ったりしてはおかしなことになる。この点を確認しよう。

 &math(
\frac{d}{dt}\iiint|\psi(\bm r,t)|^2\,d\bm r
&=\iiint\psi^*\frac{\PD\psi}{\PD t}+\frac{\PD\psi^*}{\PD t}\psi\,d\bm r\\
&=\frac{i\hbar}{2m}\iiint\psi^*(\nabla^2\psi)-(\nabla^2\psi^*)\psi\,d\bm r\\
&=\frac{i\hbar}{2m}\iiint\bm \nabla\cdot\big[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi\big]\,d\bm r\\
&=\frac{i\hbar}{2m}\iint_S\big[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi\big]\cdot \bm n\,dS\\
);

最後の等式はガウスの定理を用いて体積積分を面積積分に直した。
ここでは積分範囲を無限大に取っているためその表面というのは考えづらいが、
&math(\psi); が有限範囲内のみゼロでない値を取る場合を想定し、
それより大きな範囲を取ると考えれば問題ない。

そのような場合、表面 &math(S); 上で積分内の関数はゼロになるため、

 &math(\frac{d}{dt}\iiint|\psi(\bm r,t)|^2\,d\bm r=0);

シュレーディンガー方程式に従った時間発展では全確率密度が保存されることが分かった。

** 確率密度の流れ [#tff807cf]

上記の計算において積分領域を有限に取れば、
右辺の積分は領域外から領域内へ入ってくる確率密度であると解釈できる。((ピンと来なければ電磁気学で学んだ電荷の保存則の部分を復習せよ))
すなわち、積分内の符号を反転した

 &math(
\bm S(\bm r,t)&=-\frac{i\hbar}{2m}\big[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi\big]);

あるいはこれを変形した、

 &math(
\bm S(\bm r,t)&=\frac{\hbar}{2mi}\big[\psi^*(\bm \nabla\psi)-\{\psi^*(\bm \nabla\psi)\}^*\big]\\
&=\frac{\hbar}{2mi}\cdot 2\,\mathrm{Im}\Big[\psi^*(\bm \nabla\psi)\Big]\\
&=\mathrm{Re}\Big[\psi^*\frac{\hbar\bm \nabla}{mi}\psi\Big]\\
);

が局所的な確率密度の流れを表わす(単位面積当たり、単位時間当たりの流量)。

多数の電子が同じ波動関数に従って移動しているような場合には、
この &math(\bm S); に素電荷 &math(e); を掛けた物が電流密度になるなど、
物理的な意味も大きい。

** 座標について [#r1e087ba]

&math(|\psi(\bm r,t)|^2); が確率密度関数であるならば、
電子の位置の &math(x); 座標の期待値 &math(\langle x(t)\rangle); は

 &math(\langle x(t)\rangle=\iiint x|\psi(\bm r,t)|^2\,d\bm r);

として求められる。

二重スリットの実験のように電子の存在確率が広範囲に広がっている場合には
この期待値にあまり意味はない。一方、ポテンシャルにより束縛された電子や、
波束として空間中を飛ぶ電子ではこれがそのまま電子の位置の時間変化を表わすことになる。

もちろんその場合にも波動関数は広がりを持つから、任意の時刻の波動関数の広がりを評価するために 
&math(x); 座標の標準偏差の期待値を求めるようなことも意味を持つ。

 &math(\sigma_x(t)=\sqrt{\iiint (x-\langle x\rangle)^2|\psi(\bm r,t)|^2\,d\bm r});

このように一旦波動関数が求まれば、座標 &math(\bm r); の関数 &math(g(\bm r)); として表わされる任意の関数について、

 &math(\langle g(\bm r)\rangle=\iiint g(\bm r)|\psi(\bm r,t)|^2\,d\bm r);

としてその期待値を求めることが可能である。

もう1つ例を挙げよう。波動関数 &math(\psi(\bm r,t)); に従った空間分布を持つ電子が
あるポテンシャル &math(U(\bm r,t)); の中に存在すると、
ポテンシャルエネルギーは場所により異なるから、
この電子のポテンシャルエネルギーを何らかの方法で測った場合、
一般にその値は確率的にしか決まらない。すなわち、測定毎に異なる値を取りうる。

しかしその期待値を波動関数から求めることができて、

 &math(\langle U(t)\rangle=\iiint U(\bm r,t)|\psi(\bm r,t)|^2\,d\bm r);

となる。

** 運動量について [#o650902c]

量子力学によれば、
電子の位置が確率的にしか決定されないのと同様に、
その運動量も確率的にしか決定されない。

位置の確率密度関数は波動関数により与えるとして、
運動量の確率分布はどのようなものになるだろうか?

シュレーディンガー方程式をでっち上げた際に、平面波に対して
&math(\bm p=\hbar\bm k); および &math(\bm \nabla \psi(\bm r,t)=i\bm k\psi(\bm r,t)); 
が成り立つことを利用した。

実は平面波でない場合にも、運動量 &math(\bm p); の期待値を、

 &math(
\langle \bm p(t)\rangle=\iiint \psi^*(\bm r,t) (\hbar\bm \nabla/i) \psi(\bm r,t)\,d\bm r
);

として求めることができる(こうして求めると実験と合う)。
左から &math(\psi^*(\bm r,t)=\psi(\bm r,t)^*); がかかっているのは、&math(|\psi|^2=\psi^*\psi); 
であることに対応している。

同様に、&math(\bm p); の関数で与えられる任意の物性値 &math(g(\bm p)); の期待値は、
&math(g); の中の &math(\bm p); を &math(\hbar\bm\nabla/i); で置き換えた演算子 &math(\hat g); を作り、

 &math(
\langle g(\bm p)\rangle=\iiint \psi^*(\bm r,t) \hat g(\hbar\bm\nabla/i) \psi(\bm r,t)\,d\bm r
);

とすれば求められる。

** 一般の場合 [#tfcb6813]

解析力学の授業において、1つの粒子の運動は座標と運動量によって完全に記述されることを学んだ。

すなわち電子の運動を表わす任意の物性値は座標と運動量の関数として 
&math(g(\bm r,\bm p)); の形に書ける。波動関数からこのような物性値の期待値を求めるには、
&math(\bm p); を &math(\hbar\bm\nabla/i); で置き換えた演算子 &math(\hat g); を作り、

 &math(
\langle g(\bm r,\bm p)\angle=\iiint \psi^*(\bm r,t) \hat g(\bm r,\hbar\bm\nabla/i) \psi(\bm r,t)\,d\bm r
\langle g(\bm r,\bm p)\rangle=\iiint \psi^*(\bm r,t) \hat g(\bm r,\hbar\bm\nabla/i) \psi(\bm r,t)\,d\bm r
);

とすればよい。

このようにして求めた &math(g(\bm r,\bm p)); の期待値が実数になることはそれほど自明ではない。
後に任意の物理量を表わす演算子がエルミートになることを使って期待値が実数になることを理解する。

** 確率密度の流れ(再) [#oc57b927]
** 演習: エネルギーの場合 [#df384d23]

&math(\hbar\bm \nabla/i); は運動量を表わす演算子であり、
&math(\hbar\bm \nabla/mi); はそれを質量 &math(m); で割った物であるから、
&math(\bm v=\bm p/m); よりこれは速度を表わす演算子となる。すなわち、確率密度の流れを
系のエネルギーを一般座標と、それに対応する一般運動量で書き表したものを、
その系のハミルトニアンと呼ぶ。

 &math(
\bm S(\bm r,t)&=\mathrm{Re}\Big[\psi^*\hat{\bm v}\psi\Big]\\
);
座標を &math(\bm r);, 運動量を &math(\bm p); として、ハミルトニアンを 
&math(H(\bm r,\bm p); と表わすと、シュレーディンガー方程式に現れる 
&math(\hat H(\bm r,\hbar\bm\nabla/i)=\frac{\hbar}{2m}\nabla^2+U(\bm r,t)); 
はちょうどハミルトニアンの &math(\bm p); を演算子 &math(\hbar\bm\nabla/i); 
で置き換えた物になっている。

と書いておくと覚えやすい。
したがって、系のエネルギーの期待値は

** 不確定性原理 [#j06cf6f5]

量子力学の重要な帰結として「不確定性原理」がある。

不確定性原理は例えば、「&math(x); と &math(p_x); が同時に正確に定まるような状態は存在しない」という形で言い表せる。

[[計算の詳細は別掲するとして>量子力学I/不確定性原理]]、シュレーディンガー方程式より導かれる結果は次のようになる。

 &math(
\sigma_x\cdot\sigma_{p_x} = 
\sqrt{\big\langle (x-\langle x\rangle)^2\big\rangle}
\sqrt{\big\langle (p_x-\langle p_x\rangle)^2\big\rangle} \ge \frac{\hbar}{2}
\overline E=\overline{H(\bm r,\bm p,t)}=\iiint \psi^*(\bm r,t)\hat H(\bm r,\hbar\bm\nabla/i,t)\psi(\bm r,t)d\bm r
);

どんな波動関数に対してもこの不等式が成り立つことは、
そもそも電子自身が「上記の不等式以上に定まった位置や運動量を持つことがない」
ことを示している。
となる。

このことと、位置や運動量を測定する際にどうしても有限の誤差を生じるため、
位置や運動量を「測定によって正確に決定できない」ということとは、
しばしば混同して議論されるが、本来異なる話であるため分けて考える必要がある。
さて、時間に依らないシュレーディンガー方程式の固有値 &math(E_k); に対する固有関数を
&math(\psi_k(\bm r,t)); であるとする。

実際、ハイゼンベルグが当初議論したのは位置 &math(x); を測定する際の測定誤差 &math(\epsilon_x); と、
その位置の測定により運動量 &math(p); が受ける撹乱 &math(\eta_p); の積 &math(\epsilon_x\eta_p); 
に関するものであり、ハイゼンベルグの思考実験においては &math(\epsilon_x\eta_p\ge\hbar/2); 
が得られた。
この &math(\psi_k); に対して、

しばしば量子ゆらぎに関する不確定性  &math(\sigma_x\cdot\sigma_{p_x}\ge\hbar/2); と、~
測定の及ぼす撹乱に関する不確定性 &math(\epsilon_x\eta_p\ge\hbar/2); は混同されてきたが、~
特に後者については測定のしかたによっては正しくなく、
上記の不等式を下回る実験が可能であることが近年報告されている。
((http://dx.doi.org/10.1103/PhysRevA.67.042105))
(1) エネルギーの期待値が &math(E_k); となることを示せ。

*** 波動関数から理解する不確定性の意味 [#d6a6e62e]
(2) エネルギーの分散 (&math(\sigma_E^2); : 標準偏差の自乗) がゼロになることを示せ。

運動量の標準偏差がゼロであるような波動関数はどのようなものであろうか?
このように、&math(\psi_k); に対しては測定ごとにエネルギーの測定値が揺らぐことはなく、
常に一定値 &math(E_k); が得られることになる。

&math(\hat p-\langle p\rangle); がエルミートであることを用いて式変形すれば、
** 解答 [#y53b9710]

 &math(
\sigma_p^2
&=\iiint\psi^*\left(\hat p-\langle p\rangle\right)^2\psi\,d\bm r\\
&=\iiint\left(\hat p\psi-\langle p\rangle\psi\right)^*
        \left(\hat p\psi-\langle p\rangle\psi\right)\,d\bm r\\
&=\iiint\left|\hat p\psi-\langle p\rangle\psi\right|^2\,d\bm r\\
&=0
);
(1)

となり、すなわち至る所で

 &math(\hat p\psi=\frac{\hbar}{i}\frac{\PD\psi}{\PD x}=\langle p\rangle\psi);

が要求される。これは、

 &math(\psi(\bm x,t)=\psi(\bm 0,t)e^{i\langle p\rangle x/\hbar});

を表わしており、このとき &math(|\psi(\bm x,t)|^2=|\psi(\bm 0,t)|^2); となり
波動関数の絶対値は &math(x); によらず一定値を取る。すなわち波動関数は &math(x); の全範囲に広がってしまう。

* エーレンフェストの定理 [#wbe3d8ad]

初期状態において電子の存在確率があまり広がっておらず、
その広がりに対してポテンシャル &math(V(\bm r,t)); の変化が十分に緩やかであれば、
電子の運動は古典論から予想されるものと等しくなるはずである。
このことを確かめてみよう。

まず、電子の位置座標の時間変化を求める。

 &math(
\frac{d}{dt}\langle x\rangle
&=\frac{d}{dt}\iiint \psi^*x\psi\,d\bm r\\
&=\iiint\psi^* x\frac{\PD\psi}{\PD t}+\frac{\PD\psi^*}{\PD t}x\psi\,d\bm r\\
\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
);

シュレーディンガー方程式により &math(\PD/\PD t); を書き直して、
(2)

 &math(
\frac{d}{dt}\langle x\rangle
&=\iiint \psi^*x\frac{1}{i\hbar}\Big(-\frac{\hbar^2}{2m}\nabla^2\psi+V\psi\Big)\,d\bm r\\
&+\iiint \frac{1}{-i\hbar}\Big(-\frac{\hbar^2}{2m}\nabla^2\psi^*+V\psi^*\Big)x\psi\,d\bm r\\
\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(V); を含む項は互いに打ち消し合って、
上で見たように &math(\overline E=E_k); であるから、&math(\sigma_E=0); となる。

 &math(
\frac{d}{dt}\langle x\rangle
&=\frac{i\hbar}{2m}\iiint \psi^*x(\nabla^2\psi)-(\nabla^2\psi^*)x\psi\,d\bm r\\
);

ここで部分積分により、

 &math(
\iiint(\nabla^2\psi^*)x\psi\,d\bm r 
&=-\iiint (\bm \nabla\psi^*)\cdot\bm \nabla(x\psi)\,d\bm r
  +\underbrace{\int_S(\bm \nabla\psi^*\cdot x\psi)\cdot\bm n\,dS}_{=0}\\
&=\iiint \psi^*\nabla^2(x\psi)\,d\bm r
  +\underbrace{\int_S(\psi^*\bm\nabla(x\psi))\cdot\bm n\,dS}_{=0}\\
);

したがって、

 &math(
\frac{d}{dt}\langle x\rangle
&=\frac{i\hbar}{2m}\iiint \psi^*[x(\nabla^2\psi)-\nabla^2(x\psi)]\,d\bm r\\
&=\frac{i\hbar}{2m}\iiint \psi^*\left(-2\frac{\PD}{\PD x}\right)\psi\,d\bm r\\
&=\frac{1}{m}\langle p_x\rangle
);

を得る。ここで、

 &math(
\nabla^2(x\psi)
&=\left(\frac{\PD^2}{\PD x^2}+\frac{\PD^2}{\PD y^2}+\frac{\PD^2}{\PD z^2}\right)(x\psi)\\
&=\frac{\PD\psi}{\PD x}+\frac{\PD\psi}{\PD x}+x\frac{\PD^2\psi}{\PD x^2}+x\frac{\PD^2\psi}{\PD y^2}+x\frac{\PD^2\psi}{\PD z^2}\\
&=\left(2\frac{\PD}{\PD x}+x\nabla^2\right)\psi
);

を用いた。一方、運動量の時間変化は、

 &math(
\frac{d}{dt}\langle p_x\rangle
&=\frac{d}{dt}\iiint\psi^*\frac{\hbar}{i}\frac{\PD}{\PD x}\psi\,d\bm r\\
&=\frac{\hbar}{i}\iiint\psi^*\frac{\PD}{\PD x}\frac{\PD}{\PD t}\psi+\left(\frac{\PD}{\PD t}\psi^*\right)\frac{\PD}{\PD x}\psi\,d\bm r\\
&=\frac{\hbar}{i}\iiint\psi^*\frac{\PD}{\PD x}\frac{1}{i\hbar}\left(-\frac{\hbar^2}{2m}\nabla^2\psi+V\psi\right)\,d\bm r
+\frac{\hbar}{i}\iiint\frac{1}{-i\hbar}\left(-\frac{\hbar^2}{2m}\nabla^2\psi^*+V\psi^*\right)\frac{\PD}{\PD x}\psi\,d\bm r\\
);

ここで、

 &math(
\iiint\psi^*\frac{\PD}{\PD x}\nabla^2\psi\,d\bm r
&=\iiint\psi^*\nabla^2\frac{\PD\psi}{\PD x}\,d\bm r\\
&=-\iiint\bm\nabla\psi^*\cdot\bm\nabla\frac{\PD\psi}{\PD x}\,d\bm r
+\underbrace{\int_S\left(\psi^*\bm\nabla\frac{\PD\psi}{\PD x}\right)\cdot\bm n\,dS}_{=0}\\
&=\iiint(\nabla^2\psi^*)\frac{\PD\psi}{\PD x}\,d\bm r
-\underbrace{\int_S\left(\bm\nabla\psi^*\cdot\frac{\PD\psi}{\PD x}\right)\cdot\bm n\,dS}_{=0}\\
&=\iiint(\nabla^2\psi^*)\frac{\PD\psi}{\PD x}\,d\bm r
);

より &math(\nabla^2); を含む項は打ち消し合って、

 &math(
\frac{d}{dt}\langle p_x\rangle
&=-\iiint\psi^*\left(\frac{\PD}{\PD x}(V\psi)-V\frac{\PD}{\PD x}\psi\right)\,d\bm r\\
&=-\iiint\psi^*\frac{\PD V}{\PD x}\psi\,d\bm r\\
&=-\left\langle\frac{\PD V}{\PD x}\right\rangle
);

となり、&math(\bm r,\bm p); の期待値が古典論の運動方程式

 &math(\frac{d\bm r}{dt}=\frac{\bm p}{m});

 &math(\frac{d\bm p}{dt}=-\bm\nabla V);

を満たすことが示された。

巨視的極限に於いてシュレーディンガー方程式が古典論の運動方程式を与えるという
この定理をエーレンフェストの定理と呼ぶ。

* メモ [#if68c58e]

上記干渉のアニメーションを求めるための Mathematica コード:

 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分単位)がかかる。

* 質問・コメント [#j43557b0]

#article_kcaptcha


Counter: 174241 (from 2010/06/03), today: 6, yesterday: 0