量子力学Ⅰ/波動関数の解釈 のバックアップの現在との差分(No.3)
更新- バックアップ一覧
- 差分 を表示
- ソース を表示
- バックアップ を表示
- 量子力学Ⅰ/波動関数の解釈 へ行く。
- 追加された行はこの色です。
- 削除された行はこの色です。
[[量子力学I]] [[前の単元 <<<>量子力学Ⅰ/電子の波動方程式]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/箱の中の自由粒子]]~ ~ * 目次 [#k467238b] [[量子力学Ⅰ]]&katex(); #contents * 波動関数の物理的な意味は何か? [#e131650a] シュレーディンガー方程式より、定常状態にある水素原子中の電子の取り得る飛び飛びのエネルギー値 (エネルギー準位と呼ぶ)が &math(\hat H); の固有値として求まった。 シュレーディンガー方程式が正しいとして、そこから求まる波動関数の物理的意味は何であろうか? このとき固有値に対応して固有関数も求まるのであるが、この固有関数=波動関数の物理的意味は何であろうか? ** 話の進み方が奇妙? [#fd91a38a] 近代の物理学ではこのように、「正しそうな方程式が先に求まって、後からその意味するところを考える」 近代の物理学ではこのように、「正しそうな方程式が先に求まって、後からその意味を考える」 という手順を踏むことがしばしばある。 これまで学んできた古典的な物理ではありえないスタンスで、始めは戸惑うかもしれない。 これまでであれば直感的に納得のいく方程式が基礎方程式となっていたから、それを学ぶことに疑問はなかった。 量子力学のように極限的な状況でのみその真価が現れるようなつ物理学を進める場合、 そういった状況に置ける物理現象に対して我々の直感の働かないため、 直感を廃したこういうややこしい手順を踏むしかない、と考えて納得して欲しい。 実際、量子力学が役に立つような分野では、我々の直感とはかけ離れたような物理現象がしばしば観測される。 「直感的に理解できる」方程式が基礎となっていた古典論と比べると、 このようなやりかたは始め奇妙に思えるが、 量子力学の真価が現れるような極小の世界に対しては「我々の直感は正しく働かない」のだからしかたがない、 と考えて納得してほしい。 実際これから直感とはかけ離れた物理現象を学ぶことになる。 + 観測された物理現象を説明できそうで、かつ論理的に矛盾のない基礎方程式を適当にでっち上げる 近代科学ではむしろ「直感」をなるべく廃した、以下のような手法が科学的と考えられている。((例えば [[中谷 宇吉郎 著 『科学の方法』 (岩波新書)>https://www.iwanami.co.jp/book/b267447.html]] )) + 観測された物理現象を説明できそうで、かつ%%%論理的に矛盾のない%%%基礎方程式をでっち上げる + その方程式から何が予測されるかを考えて、新たな測定結果と突き合わせる + 矛盾が生じれば理論に修正が必要となる + 矛盾無く測定結果を説明できている限り、それが正しい理論である + 予想と異なる測定結果が得られた場合には理論に修正が必要となる + なぜそのような方程式が成り立つか、など、%%%''実験事実により検証できない内容については議論しない''%%% このような意味で現在受け入れられているのが、 先に %%求めた%% でっち上げた シュレーディンガー方程式である。 (相対論的効果やスピンに関する効果が現れる現象などに対してはさらなる修正が必要になる→ディラック方程式やさらに進んだ量子力学へ) * 測定される物理現象 [#j6d4afe8] 電子の波動関数を理解する上で最も分かりやすい実験結果の1つに、 電子の波動関数の意味を理解する上で最も分かりやすい実験結果の1つに、 二重スリットを通る電子についての測定が挙げられる。 二重スリットの実験は光学でもよく知られており、ホイヘンスの原理によれば、 スリット後の状況は2つのスリットの位置にそれぞれ点光源があるのと同じであり、 二重スリットの実験は光学でもよく知られており、[[ホイヘンスの原理>Wikipedia:ホイヘンスの原理]]によれば、 スリット通過後の状況は2つのスリットの位置にそれぞれ光源があるのと同じであり、 それらの光の干渉によりスクリーン上に干渉縞が現れる。 当然、片方のスリットを閉じれば干渉縞は消失する。 &attachref(double-slit.png,,25%); &attachref(double-slit.gif,,50%); &attachref(double-slit.png,,25%); &attachref(double-slit.gif,,50%,ogp);&qr(http://dora.bk.tsukuba.ac.jp/~takeuchi/?plugin=attach&refer=%E9%87%8F%E5%AD%90%E5%8A%9B%E5%AD%A6%E2%85%A0%2F%E6%B3%A2%E5%8B%95%E9%96%A2%E6%95%B0%E3%81%AE%E8%A7%A3%E9%87%88&openfile=double-slit.gif); 電子に対しても同様の実験を行うことができて、その結果、光と同様に干渉縞が現れる。 #ref(Double-slit_experiment_results_Tonomura_Preview.jpg,right,around); #ref(Double-slit_experiment_results_Tonomura_Preview.jpg,right,around,75%); 電子が古典的な粒子であればスリットを通った電子は単に直進するのみであるから、 もう一方のスリットが空いていようが閉じていようが もう一方のスリットが開いていようが閉じていようが 1つのスリットを通った電子がスクリーン上に到達する位置は変わらないはずで、 二重スリットの結果は単一スリットの結果の足し算になるはずである。 それにもかかわらず干渉縞が現れることは、確かに電子が波の性質を持っていることを表わしてる。 それにもかかわらず干渉縞が現れることは、確かに電子が''波の性質''を持っていることを表わしてる。 それでは電子に粒子としての性質が無いのかといえばそんなことはない。 それでは電子に''粒子としての性質''が無いのかといえばそんなことはない。 高感度のスクリーンを用いることで、 1つ1つの電子がスクリーン上に到達した位置を記録することができる。 干渉縞が現れる場合にも、1つの電子は1つの点に到達する。 すなわち、電子は確かに粒子なのである。 有名な外村による二重スリットの測定結果は右図のようになり、 1つ1つの電子がスクリーン上に到達する位置はランダムに見える物の、 その確率が波動関数の干渉により波打っていることから、 有名な&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 x,t)); で記述される - 波動関数は空間的に広がりを持ち、また、干渉や回折などの波に特有な性質を現す - 電子の位置を実験的に観測した場合には電子はある一点に見出される - 電子の発見される確率は &math(|\psi(\bm x,t)|^2); に比例する - 電子の運動はシュレーディンガー方程式((あるいはそれを拡張したパウリ方程式やディラック方程式などの波動方程式))を満たす波動関数 $\psi(\bm r,t)$ で記述される - 波動関数は一般に複素関数で、空間的に広がりを持ち、また、干渉や回折などの波に特有な性質を現す - 電子の位置を実験的に観測した場合には電子はある一点に見出され、広がりを持たない - 位置 $\bm r$ の周りの微小体積 $d\bm r$ に電子が発見される確率は $|\psi(\bm r,t)|^2d\bm r$ に比例する ← %%%''ここが重要''%%% 二重スリットの実験に当てはめれば、 - 二重スリットを通る電子の波動関数は、2つのスリットのそれぞれを通る経路の間で干渉を起こす - その結果、波動関数の絶対値の二乗 &math(|\psi(\bm x,t)|^2); に干渉縞が現れる - その結果、波動関数の絶対値の二乗 $|\psi(\bm r,t)|^2$ に濃淡=干渉縞が現れる - 電子がスクリーンに当たり、その位置が記録されることが「観測」にあたる - 観測される電子位置は常に1点だけに定まり、広がりを持たない - &math(|\psi(\bm x,t)|^2); の大きな箇所でより多くの電子が発見されるため、 - $|\psi(\bm r,t)|^2$ の大きな箇所でより多くの電子が発見されるため、 多くの電子について観測を繰り返すことによりスクリーン上に干渉縞が現れる 空間的に広がりを持つ電子が観測により1点に見出される様子は「波動関数の収束」と呼ばれる。 現在の量子力学はなぜ観測により波動関数が収束を起こすのか、 観測しないときに電子は度の位置にあるのか、といった問いには答えない。 「観測によって検証できない内容」は物理の範疇ではないというスタンスである。 現在の量子力学は、なぜ観測により波動関数が収束を起こすのか、 とか、観測しないときに電子はどの位置にあるのか、といった問いには答えない。 「観測によって検証できない命題」は物理学の範疇ではないというスタンスである。 そのかわり、「観測により確かめられる内容」については量子力学は完璧な予想を与える。 (対応するシュレーディンガー方程式が数学的に解ける範囲にある限り) (対応するシュレーディンガー方程式が数学的・数値的に解ける限り → 方程式を正確に解く方法がなくて困ることも多い) また、余裕があれば、上で出てくる「確率」が「誰にとっての確率か」を考えてみると良い。 例えば有名な [[Wikipedia:シュレーディンガーの猫]] の問題では猫の生死が観測するまで 確率的にしか決まらないという状況が出てくるが、猫自身にとって自分の生死は確率に頼らず決定している。 実は量子力学で出てくる「確率」は観測者によって異なり、したがって「波動関数」も観測者によって異なるのである。 #clear * 確率密度関数について [#q2585cd6] 確率・統計の需要を取った学生には学習済みの内容であるが、 確率密度関数についてここで復習しておく。 $|\psi(\bm r,t)|^2$ は確率密度関数になる。 そこで確率密度関数についてここで復習しておく。 ** 確率密度関数の定義 [#x546ffd8] #ref(はじめての誤差論/probability_density_function.png,right,around,50%); 「確率変数 &math(x); が確率密度関数 &math(f(x)); に従う」という意味は、 「確率変数 $x$ が確率密度関数 $f(x)$ に従う」という意味は、 測定毎に異なる値を取る変数 &math(x); があり、1回測定したときに &math(x); が &math(x_a<x<x_b); の範囲に入る確率を 測定毎に異なる値を取る変数 $x$ があり、1回測定したときに $x$ が $x_a<x<x_b$ の範囲に入る確率を &math(P\{x_a<x<x_b\}=\int_{x_a}^{x_b}f(x)dx); $$\begin{aligned}P\{x_a<x<x_b\}=\int_{x_a}^{x_b}f(x)dx\end{aligned}$$ として求められる、ということである。 上図で分かるとおり、上式の右辺の積分は &math(f(x)); と &math(x); 軸に挟まれる部分の面積に相当するため、 &math(f(x)); が大きいことは &math(x); がその付近の値を取りやすいことを表わしており、 &math(f(x)); が小さければ &math(x); はその付近の値を滅多に取らない。 右図で分かるとおり、上式の右辺の積分は $f(x)$ と $x$ 軸に挟まれる影付きの部分の面積に相当するため、 $f(x)$ が大きければ $x$ はその付近の値を取りやすく、 $f(x)$ が小さければ $x$ はその付近の値を滅多に取らない。 当然、すべての &math(x); に対して &math(f(x)\ge 0); であり、 また、測定を行えば必ず何らかの値が得られることから、&math(f(x)); を全範囲にわたって積分した値は常に1になる。 当然、すべての $x$ に対して $f(x)\ge 0$ であり、 また、測定を行えば必ず何らかの値が得られることから、$f(x)$ を全範囲にわたって積分した値は常に1でなければならない。 &math(\int_{-\infty}^{\infty}f(x)dx=1); $$\begin{aligned}\int_{-\infty}^{\infty}f(x)dx=1\end{aligned}$$ ** 期待値 [#hd633e60] 確率変数 &math(x); の測定を仮想的に無限回行った場合に得られる平均値を &math(x); の期待値という。 確率変数 $x$ の測定を仮想的に無限回行った場合に得られるであろう「$x$ の平均値」を $x$ の期待値という。 例えば、サイコロを1回振った際に出る目の期待値は、 例えば、サイコロを1つ振った際に出る目の期待値は、それぞれの目が $1/6$ の確率で現れるのだから、 &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); $$\begin{aligned}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\end{aligned}$$ のように、サイコロの目の値に、個々の目の出る確率を掛けて得られる。 すなわち、離散的な値を取る確率変数 &math(x); の期待値 &math(\overline x); は、 すなわち、離散的な値を取る確率変数 $x$ の期待値 $\overline x$ は、 &math(\overline{x}=\sum_{k} x_k p(x_k)); $$\begin{aligned}\overline x=\sum_{k} \underbrace{x_k\vphantom{p(x_k)}}_\text{値} \underbrace{p(x_k)}_\text{その確率}\end{aligned}$$ として与えられる。ここで &math(p(x)); は &math(x=x_k); となる確率である。 として与えられる。ここで $p(x_k)$ は $x=x_k$ となる確率である。 いま考えている &math(x); が連続値を取る場合にも、これと同様に $x$ が連続値を取る場合には和が積分に変わって、 &math(\overline{x}=\int_{-\infty}^\infty xf(x)\,dx); $$\begin{aligned}\overline x=\int_{-\infty}^\infty x\underbrace{f(x)\,dx}_{p(x)に相当}\end{aligned}$$ により期待値が求まる。測定値が範囲 &math([x,x+dx]); に入る確率を &math(f(x)dx); と表せることに注意せよ。 となる。測定値が範囲 $[x,x+dx]$ に入る確率を $f(x)dx$ と表せることに注意せよ。 ** 標準偏差 [#f028b0d1] ** 分散 [#kef9a92d] 標準偏差は確率変数 &math(x); がその期待値からどれほど大きなバラツキを持つかを表わす指標である。 分散は確率変数 $x$ がその期待値を中心にどれほど大きなバラツキを持つかを表わす指標である。 統計学的には、&math(n); 回の測定値 &math(x_1,x_2,\dots,x_n); の標準偏差は、 統計学的には、$n$ 回の測定値 $x_1,x_2,\dots,x_n$ の分散は、 &math(\sigma_x=\sqrt{\frac{1}{n}\sum_{k=1}^n (x_k-\overline x)^2}); $$\begin{aligned}\sigma^2=\frac{1}{n}\sum_{k=1}^n (x_k-\overline x)^2\end{aligned}$$ として定義される。 これに対して確率変数 &math(x); を仮想的に無限回観測した際の標準偏差の期待値は これに対して確率変数 $x$ を仮想的に無限回観測した際の分散は &math(\sigma_x=\sqrt{\int_{-\infty}^\infty (x-\overline x)^2f(x)\,dx}); $$\begin{aligned}\sigma_x^2=\overline{(x-\overline x)^2} =\int_{-\infty}^\infty (x-\overline x)^2f(x)\,dx\end{aligned}$$ として求められる。これを確率密度関数 &math(f(x)); の標準偏差という。 として求められる。これを確率密度関数 $f(x)$ の分散という。 ** 任意の関数の期待値 [#zebdd06f] 分散を $y=(x-\overline x)^2$ という変数の期待値であるととらえらることもできる点に注意しておこう。 ** 標準偏差 [#f028b0d1] 標準偏差の場合と同様に、確率変数 &math(x); の関数として表わされる任意の値 &math(g(x)); の期待値を、 分散の自乗根を標準偏差 $\sigma=\sqrt{\sigma^2}$ という。 &math(\overline{g(x)}=\int_{-\infty}^\infty g(x)\,f(x)\,dx); $$\begin{aligned}\sigma_x=\sqrt{\overline{(x-\overline x)^2}} =\sqrt{\int_{-\infty}^\infty (x-\overline x)^2f(x)\,dx}\end{aligned}$$ として求められる。&math(g(x)); が線型である場合を除き、 例えば $x$ が粒子の位置を表わすとき、 分散が面積の次元を持つのに対して 標準偏差は $x$ と同じ長さの次元を持つため、 $x$ のバラツキの指標として用いやすい。 &math(\overline{g(x)}\ne g(\overline x)); ** 任意の関数の期待値 [#zebdd06f] であることに注意せよ。 分散の場合と同様に、確率変数 $x$ の関数として表わされる任意の値 $g(x)$ の期待値を、 #clear $$\begin{aligned}\overline{g(x)}=\int_{-\infty}^\infty g(x)\,f(x)\,dx\end{aligned}$$ * 波動関数から各種物性値を取り出すには [#z663cace] として求められる。($g(x)=(x-\overline x)^2$ の場合が分散 $\sigma_x^2$ にあたる) ** 波動関数の規格化 [#i56b385a] $g(x)$ が線型である場合を除き、 シュレーディンガー方程式は線型な方程式だから、 ある関数 &math(\psi_1(\bm x,t)); が解であれば、 任意の定数 &math(A); に対して &math(A\psi_1(\bm x,t)); も解である。 $$\begin{aligned}\overline{g(x)} \ne g(\overline x)\end{aligned}$$ 一方、波動関数が確率密度関数となるためには、 であることに注意せよ。 &math(\iiint|\psi(\bm x,t)|^2\,dx\,dy\,dz=\int|\psi(\bm x,t)|^2\,d^3x=1); ** 複合確率 [#wb9dcaf6] でなければならない。 確率変数 $x$ が区間 $[x, x+dx]$ に入る確率が $f(x)dx$ であり、~ 確率変数 $y$ が区間 $[y, y+dy]$ に入る確率が $g(y)dy$ であり、~ 両者が同時に起きる確率が $h(x,y)dx\,dy$ であるとする。 任意のシュレーディンガー方程式の解 &math(\psi_1(\bm x,t)); に対して、 一般には $h(x,y)dx\,dy\ne f(x)g(y)dx\,dy$ である。 &math(\psi(\bm x,t)=\left\{\int|\psi_1(\bm x,t)|^2\,d^3x\right\}^{-1}\psi_1(\bm x,t)); $h(x,y)=f(x)g(y)$ であるとき $x$ と $y$ は独立である、という。 として、確率密度関数を得る手順を波動関数の規格化と呼ぶ。 $h(x,y)\ne f(x)g(y)$ であるとき $x$ と $y$ には相関がある、という。 ** 全確率の保存 [#i7b418d2] すなわち、「確率密度関数が変数分離できる」ことは、個々の変数に対する「確率分布が独立であること」と同値である。→ $\psi(\bm r,t)=\varphi(\bm r)\tau(t)$ と表せるとき、空間分布は時間と独立である→時間によらない。 ある時刻において規格化した波動関数がシュレーディンガー方程式に従って時間発展する時、 確率密度が増えたり減ったりしてはおかしなことになる。この点を確認しよう。 #clear &math( \frac{\PD}{\PD t}\int|\psi(\bm x,t)|^2\,d^3x &=\int\psi^*\frac{\PD\psi}{\PD t}+\frac{\PD\psi^*}{\PD t}\psi\,d^3x\\ &=\frac{i\hbar}{2m}\int\psi^*(\nabla^2\psi)-(\nabla^2\psi^*)\psi\,d^3x\\ &=\frac{i\hbar}{2m}\int\bm \nabla\cdot[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi]\,d^3x\\ &=\frac{i\hbar}{2m}\int_S[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi]\cdot \bm n\,dS\\ ); * 波動関数から各種物性値を取り出すには [#z663cace] 最後の等式はガウスの定理を用いて体積積分を面積積分に直した。 ここでは積分範囲を無限大に取っているためその表面というのは考えづらいが、 &math(\psi); が有限範囲内のみゼロでない値を取る場合を想定し、 それより大きな範囲を取っていると考えればよい。 シュレーディンガー方程式を解いて波動関数 $\psi(\bm r,t)$ が得られたとする。 この波動関数から電子の運動について何をどのように知ることができるのだろうか? そのような場合、表面 &math(S); 上で積分内の関数はゼロになるため、 ** 準備:波動関数の規格化 [#i56b385a] &math(\frac{\PD}{\PD t}\int|\psi(\bm x,t)|^2\,d^3x=0); シュレーディンガー方程式は''線型な方程式''だから、 ある関数 $\psi_1(\bm r,t)$ が解であれば、 任意の定数 $A$ に対して $A\psi_1(\bm r,t)$ も解である。 シュレーディンガー方程式に従った時間発展では確率密度が保存されることが分かった。 一方、波動関数の絶対値の二乗 $|\psi(\bm r,t)|^2$ が確率密度関数となるためには、 ** 確率密度の流れ [#tff807cf] $$\begin{aligned}\iiint|\psi(\bm r,t)|^2\,d\bm r=1\end{aligned}$$ 上記の計算において &math(S); を有限に取れば、 右辺の積分は領域外から領域内へ入ってくる確率密度であると解釈できる。((ピンと来なければ電磁気学で学んだ電荷の保存則の部分を復習せよ)) すなわち、積分内の符号を反転した でなければならない。 ただし $\iiint\ d\bm r\equiv \iiint\ dx\,dy\,dz$ である。 &math( \bm S(\bm x,t)&=-\frac{i\hbar}{2m}[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi]); 任意のシュレーディンガー方程式の解 $\psi_1(\bm r,t)$ に対して、 あるいはこれを変形して、 $$\begin{aligned}\psi(\bm r,t)=\frac{1}{\sqrt{\iiint|\psi_1(\bm r,t)|^2\,d\bm r}}\psi_1(\bm r,t)\end{aligned}$$ &math( \bm S(\bm x,t)&=\frac{\hbar}{2mi}[\psi^*(\bm \nabla\psi)-\{\psi^*(\bm \nabla\psi)\}^*]\\ &=\frac{\hbar}{2mi}\cdot 2\,\mathrm{Im}\Big[\psi^*(\bm \nabla\psi)\Big]\\ &=\mathrm{Re}\Big[\psi^*\frac{\hbar\bm \nabla}{mi}\psi\Big]\\ ); のように波動関数を規格化し、$|\psi(\bm r,t)|^2$ が確率密度関数となるようにできる。 線形代数II を習った学生は、上式がベクトルをそのノルムで割って規格化することに対応しているのを理解せよ。 が局所的な確率密度の流れを表わす(単位面積当たり、単位時間当たりの流量)。 $$\begin{aligned}\iiint|\psi(\bm r,t)|^2\,d\bm r &=\iiint\Big|\frac{1}{\sqrt C}\psi_1(\bm r,t)\Big|^2\,d\bm r\\ &=\frac{1}{C}\underbrace{\iiint|\psi_1(\bm r,t)|^2,d\bm r}_{=\,C}\\[-3mm] &=1\end{aligned}$$ 多数の電子が同じ波動関数に従って移動しているような場合には、 この &math(\bm S); に素電荷 &math(e); を掛けた物が電流密度になるなど、 物理的な意味も大きい。 ** 座標について [#r1e087ba] 時刻 &math(t); において電子が &math(\bm x); に見出される確率が 波動関数 &math(\psi(\bm x,t)); によって &math(|\psi(\bm x,t)|^2); により表わされるとすれば、 この電子の位置の &math(x); 座標の期待値 &math(\overline x); は $|\psi(\bm r,t)|^2$ が確率密度関数であるならば、 電子の位置の $x$ 座標の期待値 $\overline x(t)$ は &math(\overline x(t)=\iiint x|\psi(\bm x,t)|^2\,dx\,dy\,dz=\int x|\psi(\bm x,t)|^2\,d^3x); $$\begin{aligned}\overline x(t)=\iiint x|\psi(\bm r,t)|^2\,d\bm r\end{aligned}$$ として求められる。 二重スリットの実験のように電子の存在確率が広範囲に広がってしまっている場合には 二重スリットの実験のように電子の存在確率が広範囲に広がっている場合には この期待値にあまり意味はない。一方、ポテンシャルにより束縛された電子や、 波束として空間中を飛ぶ電子ではこれがそのまま電子の位置の時間変化を表わすことになる。 波束として空間中を飛ぶ電子では、この値が電子の位置の時間変化を表わすことになる。 もちろんその場合にも波動関数は広がりを持つから、任意の時刻の波動関数の広がりを評価するために &math(x); 座標の標準偏差の期待値を求めるようなことも意味を持つ。 もちろんその場合にも波動関数は広がりを持つから、任意の時刻の''波動関数の広がり''を評価するために $x$ 座標の標準偏差を求めるようなことも意味を持つ。 &math(\sigma_x(t)=\sqrt{\int (x-\overline x)^2|\psi(\bm x,t)|^2\,d^3x}); $$\begin{aligned}\sigma_x(t)=\sqrt{\overline{(x-\overline x)^2}}=\sqrt{\iiint (x-\overline x)^2|\psi(\bm r,t)|^2\,d\bm r}\end{aligned}$$ このように一旦波動関数が求まれば、座標 &math(\bm x); の関数 &math(g(\bm x)); として表わされる任意の関数について、 &math(\overline{g(\bm x)}=\int g(\bm x)|\psi(\bm x,t)|^2\,d^3x); としてその期待値を求めることが可能である。 ** 運動量について [#o650902c] 量子力学によれば電子の位置が確率的にしか決定されないのと同様に、 その運動量も確率的にしか決定されない。 電子の位置が確率的にしか決定されないのと同様に、 量子力学によれば運動量も確率的にしか決定されない。 波動関数から運動量を得るにはどうしたらよいだろうか? 波動関数は位置の確率密度関数を与えるが、 運動量の確率密度関数はどのようなものになるだろうか? シュレーディンガー方程式をでっち上げた際に、平面波に対して &math(\bm p=\hbar\bm k); および &math(\bm \nabla \psi(\bm x,t)=i\bm k\psi(\bm x,t)); $\bm p=\hbar\bm k$ および $\hbar\bm k\psi(\bm r,t)=\frac{\hbar}{i}\bm \nabla \psi(\bm r,t)$ が成り立つことを利用した。 実は平面波でない場合にも、運動量 &math(\bm p); の期待値を、 実は平面波でない場合にも、運動量 $\bm p$ の期待値を、 &math( \overline{\bm p}(t)=\int \psi(\bm x,t)^* (\hbar\bm \nabla/i) \psi(\bm x,t)\,d^3x ); $$\begin{aligned} \overline{\bm p}(t)=\iiint \psi^*(\bm r,t) \frac{\hbar}{i}\bm \nabla \psi(\bm r,t)\,d\bm r \end{aligned}$$ として求めることができる(こうして求めると実験と合う)。 左から &math(\psi(\bm x,t)^*); がかかっているのは、&math(|\psi|=\psi^*\psi); であることに対応している。 として求めることができる。この式は、 同様に、&math(p); の関数で与えられる任意の物性値 &math(g(\bm p)); の期待値は、 &math(g); の中の &math(\bm p); を &math(\hbar\nabla/i); で置き換えた演算子 &math(\hat g); を作り、 $$\begin{aligned}\bm p\sim\hbar\bm k\sim\frac{\hbar}{i}\bm \nabla\end{aligned}$$ &math( \overline{g(\bm p)}=\int \psi(\bm x,t)^* \hat g(\hbar\nabla/i) \psi(\bm x,t)\,d^3x ); の両辺に $\psi(\bm r,t)$ を付けた とすれば求められる。 $$\begin{aligned}\bm p \psi(\bm r,t)\sim\frac{\hbar}{i}\bm \nabla\psi(\bm r,t)\end{aligned}$$ ** 一般の場合 [#tfcb6813] と、 解析力学の授業において、1つの粒子の運動は座標と運動量によって完全に記述されることを学んだ。 $$\begin{aligned} \overline{\bm r}(t)&=\iiint \bm r|\psi(\bm r,t)|^2\,d\bm r\\ &=\iiint \psi^*(\bm r,t) \bm r \psi(\bm r,t)\,d\bm r \end{aligned}$$ すなわち電子の運動を表わす任意の物性値は座標と運動量の関数として &math(g(\bm x,\bm p)); の形に書ける。波動関数からこのような物性値の期待値を求めるには、 &math(\bm p); を &math(-\hbar\nabla); で置き換えた演算子 &math(\hat g); を作り、 とを見比べると、なんとなく納得できないこともないはず。 &math( \overline{g(\bm x,\bm p)}=\int \psi(\bm x,t)^* \hat g(\bm x,-i\hbar\nabla) \psi(\bm x,t)\,d^3x ); $$\begin{aligned} \overline{\bm p}(t) &=\iiint \psi^*(\bm r,t) \frac{\hbar}{i}\bm \nabla \psi(\bm r,t)\,d\bm r\\ &\sim\iiint \psi^*(\bm r,t) \bm p \psi(\bm r,t)\,d\bm r\\ &=\iiint \bm p |\psi(\bm r,t)|^2\,d\bm r\\ \end{aligned}$$ とすればよい。 という雰囲気だ。 このようにして求めた &math(g(\bm x,\bm p)); の期待値が実数になることはそれほど自明ではない。 後に任意の物理量を表わす演算子がエルミートになることを使って期待値が実数になることを理解する。 と言われても、現時点ではそもそも上記の計算結果が実数になるのかどうかからあやふやであるのだが、 後に学ぶように上記の積分は現実的な条件下で必ず実数値となることを数学的に証明でき、 また、計算から予想される値は実験結果を良く再現する。 ** 確率密度の流れ(再) [#oc57b927] ** 一般の場合 [#tfcb6813] &math(\hbar\bm \nabla/i); は運動量を表わす演算子であり、 &math(\hbar\bm \nabla/mi); はそれを質量 &math(m); で割った物であるから、 &math(\bm v=\bm p/m); よりこれは速度を表わす演算子となる。すなわち、確率密度の流れを 解析力学の授業では、1つの質点の運動は座標と運動量によって完全に記述されることを学んだ。 &math( \bm S(\bm x,t)&=\mathrm{Re}\Big[\psi^*\hat{\bm v}\psi\Big]\\ ); すなわちある時刻の座標と運動量が分かればその前後の質点の運動は完全に決定されるし、 その他の任意の物性値は座標と運動量の関数として $g(\bm r,\bm p)$ の形に書ける。波動関数からこのような物性値の期待値を求めるには、 $\bm p$ を $\hbar\bm\nabla/i$ で置き換えた演算子 $\hat g$ を作り、 と書いておくと覚えやすい。 $$\begin{aligned} \overline{g(\bm r,\bm p)}=\iiint \psi^*(\bm r,t) \hat g(\bm r,\hbar\bm\nabla/i) \psi(\bm r,t)\,d\bm r \end{aligned}$$ ** 不確定性原理 [#j06cf6f5] とすればよい。 量子力学の重要な帰結として「不確定性原理」がある。 不確定性原理は例えば、&math(x); と &math(p_x); が同時に正確に定まるような状態は存在しない、という形で言い表せる。以下、この項では &math(p_x=p); と書き、一次元で考える。 このようにして求めた $g(\bm r,\bm p)$ の期待値も現実的な条件の下で実数になることを後に学ぶ。 *** 例:エネルギーの場合 [#df384d23] 上記を数学的に表わせば、&math(\sigma_x\cdot\sigma_p); に最小値があり、ある一定値以下にはならない、ということになる。 ''ハミルトニアン'' は系のエネルギーを一般座標と一般運動量で書き表したものであった。 &math( (\sigma_x\cdot\sigma_p)^2 =\sigma_x^2\sigma_p^2=\big\langle (x-\langle x\rangle)^2\big\rangle \big\langle (p-\langle p\rangle)^2\big\rangle ); 座標を $\bm r$, 運動量を $\bm p$ として、ハミルトニアンを $H(\bm r,\bm p)$ と表わすと、シュレーディンガー方程式に現れる $\hat H(\bm r,\hbar\bm\nabla/i)=-\frac{\hbar^2}{2m}\nabla^2+V(\bm r,t)$ はちょうどハミルトニアンの $\bm p$ を演算子 $\hbar\bm\nabla/i$ で置き換えたものになっている。 演算子 &math(\alpha=x-\langle x\rangle);、 &math(\beta=p-\langle p\rangle=-i\hbar\frac{d}{dx}-\langle p\rangle); を導入すると、 したがって、系のエネルギーの期待値は &math( (\sigma_x\cdot\sigma_p)^2=\int \psi^*\alpha^2\psi\,dx\int \psi^*\beta^2\psi\,dx ); $$\begin{aligned} \overline \varepsilon=\overline{H(\bm r,\bm p,t)}=\iiint \psi^*(\bm r,t)\hat H\psi(\bm r,t)d\bm r \end{aligned}$$ ここで &math(\alpha^*=\alpha); より、 となる。 &math( \sigma_x^2=\int (\alpha\psi)^*(\alpha\psi)\,dx ); ここで、$\psi(\bm r,t)$ が時間に依存しないシュレーディンガー方程式の解で $\psi(\bm r,t)=\varphi(\bm r)\tau(t)$ と書けるとき、 &math(\beta); についても、 $$ \begin{cases} \hat H\varphi(\bm r)=\varepsilon\varphi(\bm r)\\ \tau^*(t)\tau(t)=1&\tau(t)=e^{-i\omega t} \end{cases} $$ &math( \int \psi^*\beta^2\psi\,dx &=\int \psi^*\left[-i\hbar\frac{d}{dx}-\langle p\rangle\right]^2\psi\,dx\\ &=\int -\hbar^2\psi^*\frac{d^2\psi}{dx^2} +2i\hbar\psi^*\langle p\rangle\frac{d\psi}{dx} +\langle p\rangle|\psi|^2\,dx\\ ); より、 であり、 $$\begin{aligned} \overline \varepsilon&=\iiint \varphi^*(\bm r)\tau^*(t)\hat H\varphi(\bm r)\tau(t)d\bm r\\ &=\varepsilon\iiint |\varphi(\bm r)|^2\,d\bm r\\ &=\varepsilon \end{aligned}$$ &math( \int (\beta\psi)^*(\beta\psi)\,dx &=\int \left[+i\hbar\frac{d\psi^*}{dx}-\langle p\rangle\psi^*\right] \left[-i\hbar\frac{d\psi }{dx}-\langle p\rangle\psi \right]\,dx\\ &=\int \hbar^2\frac{d\psi^*}{dx}\frac{d\psi }{dx} -i\hbar\frac{d\psi^*}{dx}\langle p\rangle\psi +i\hbar\langle p\rangle\psi^*\frac{d\psi }{dx} +\langle p\rangle^2|\psi|^2 \,dx\\ ); となり、先に想像したとおり、エネルギー固有値 $\varepsilon$ に対応する固有定常状態のエネルギーの期待値は固有値と一致することを確かめられる。 部分積分により、 ** 振り返ると [#ld9a038c] &math( \int\psi^*\frac{d^2\psi}{dx^2}\,dx = -\int\frac{d\psi^*}{dx}\frac{d\psi}{dx}\,dx +\underbrace{\left[\psi^*\frac{d\psi}{dx}\right]_{-\infty}^{+\infty}}_{=0} ); 上記では、「どうすれば期待値を求められるか」を説明したが、「なぜそれでいいのか」は説明していない。 実は後者は今後も説明されない。あえて理由を付けるならば、そう考えると実験結果と合うから、そしてなによりもいろいろと''つじつまが合うから''、ということになる。実のところ量子力学のような複雑な理論が「つじつまが合う」ように組み立てられているというのはすごいことで、「いくらほじくり返しても矛盾が生じない」という条件を付けるだけで理論の形は相当程度絞られるのである。 &math( \int\psi^*\frac{d\psi}{dx}\,dx = -\int\frac{d\psi^*}{dx}\psi\,dx +\underbrace{\Big[\psi^*\psi\Big]_{-\infty}^{+\infty}}_{=0} ); 以降では、物理量に対する期待値の持つ性質を学ぶことで、 極微の世界で起きている物理現象を学んでいく。 であるから、 1つ、古典論との対比をしておく。古典論では''ある時刻における粒子の位置と運動量(=2つの初期値)''を決めると、その後の粒子の運動はニュートン方程式によりただ一通りに決定されるのであった。量子力学においては''ある時刻の波動関数''を与えるだけでその後の時間発展がシュレーディンガー方程式により決定される。これは波動関数がその絶対値に位置に関する情報を、位相に運動量に関する情報を含んでいるためである。 &math( \int \psi^*\beta^2\psi\,dx=\int (\beta\psi)^*(\beta\psi)\,dx ); $$ \underbrace{\frac{\partial}{\partial t}\psi(\bm r,t)}_\text{時間変化}=\underbrace{\frac{1}{i\hbar}\Big[-\frac{\hbar^2}{2m}\nabla^2+V(\bm r,t)\Big]\psi(\bm r,t)}_\text{現在の波動関数から計算可能} $$ * 確率密度の時間変化 [#ha5eeb3a] と書ける。すなわち、 ある時刻において規格化されていた波動関数がシュレーディンガー方程式に従って時間発展する時、 全確率密度が増えたり減ったりしてはおかしなことになる。 すなわちある点での確率密度が時間と共に変化するとしても、それは確率の移動(流れ)によるものであり、全確率は保存されていなければならない。 この点を確認しよう。 &math( (\sigma_x\cdot\sigma_p)^2=\int |\alpha\psi|^2\,dx\int |\beta\psi|^2\,dx ); ** 全確率の保存 [#i7b418d2] 一般の &math(f,g); について、 全確率の時間微分を計算する。 &math( \int\Bigg|f-g\frac{\int fg^*dx}{\int|g|^2dx}\Bigg|^2dx\ge 0 ); $$ \begin{aligned} \underbrace{\frac{d}{dt}\ \underbrace{\iiint|\psi(\bm r,t)|^2\,d\bm r}_\text{=\,1\ であるべき}}_\text{=\,0\ であるべき}\ \ % &=\iiint\frac{\partial}{\partial t}\big(\psi^*\psi\big)\,d\bm r\\[-9mm] &=\iiint\left\{\psi^*\frac{\partial\psi}{\partial t}+\frac{\partial\psi^*}{\partial t}\psi\right\}d\bm r\\ &=\iiint\left\{\psi^*\frac{\partial\psi}{\partial t}+\Big(\frac{\partial\psi}{\partial t}\Big)^*\psi\right\}d\bm r\\ &=\iiint\left\{\psi^*\bigg(\frac{1}{i\hbar}\hat H\psi\bigg)+\bigg(\frac{1}{i\hbar}\hat H\psi\bigg)^*\psi\right\}d\bm r\\ &=\frac{1}{i\hbar}\iiint\left\{\psi^*\big(\hat H\psi\big)-\big(\hat H\psi^*\big)\psi\right\}d\bm r\hspace{12mm}\hat H^*=\hat H\ \text{は後に学ぶ}\\ &=\frac{i\hbar}{2m}\iiint\big\{\psi^*(\nabla^2\psi)-(\nabla^2\psi^*)\psi\big\}d\bm r\hspace{10.5mm}V(\bm r)\text{の項は打消しあう}\\ &=\frac{i\hbar}{2m}\iiint\bm \nabla\cdot\big[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi\big]\,d\bm r\hspace{7mm}\bm\nabla\psi^*\cdot\bm\nabla\psi \text{の項は打消しあう}\\ &=\frac{i\hbar}{2m}\iint_S\big[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi\big]\cdot \bm n\,dS\hspace{9mm}\text{ガウスの定理}\\ \end{aligned} $$ は積分の中が常に非負であるから必ず成立する。 等号はある定数 &math(\gamma); に対して &math(x); の全範囲において &math(f=\gamma g); となる場合のみ成り立つ。 最後の等式はガウスの定理を用いて体積積分を面積積分に直した。 ここでは積分範囲を無限大に取っているためその表面というのは考えづらいが、 $\psi$ が有限範囲内のみゼロでない値を取る場合を想定し、 それより大きな範囲を取ると考えれば問題ない。 &math( \int|f|^2\,dx -\int f^*g\,dx\frac{\int fg^*dx}{\int|g|^2dx} -\int fg^*\,dx\frac{\int f^*gdx}{\int|g|^2dx} +\int|g|^2\,dx\left(\frac{\int fg^*dx}{\int|g|^2dx}\right)^2 \ge 0 ); そのような場合、表面 $S$ 上で積分内の関数はゼロになるため、 両辺に &math(\int|g|^2dx\ge 0); を掛けて、 $$\begin{aligned}\frac{d}{dt}\iiint|\psi(\bm r,t)|^2\,d\bm r=0\end{aligned}$$ &math( \int|f|^2dx\int|g|^2dx -\int f^*g\,dx\int fg^*dx \underbrace{-\int fg^*dx\int f^*g\,dx +\left(\int fg^*dx\right)^2}_{=0} \ge 0 ); シュレーディンガー方程式に従った時間発展では全確率が保存されることが分かった((全確率が保存されることは、線形代数学的にはハミルトニアンがエルミートであり、時間発展演算子がユニタリであることにに対応する。発展的内容に興味があればこちらを参照 → [[線形代数を用いた証明>量子力学Ⅰ/波動関数の解釈/メモ#f6db1a8f]]))。これは上で述べた「つじつまが合う」の一例になっている。 したがって、 以下に見るように、閉曲面 $S$ 上で $\psi$ がゼロでないとき、最後の表面積分はこの面を通って外から中へ流れ込む単位時間当たりの確率を表す。 &math(\int|f|^2dx\int|g|^2dx\ge\left|\int f^*g\,dx\right|^2); ** 確率の流れの大きさ(確率流密度) [#tff807cf] を得る。&math(f=\alpha\psi);、&math(g=\beta\psi); とすれば、 実は上記の計算は積分記号なしのまま進められて、 &math( (\sigma_x\cdot\sigma_p)^2 &\ge\left|\int (\alpha\psi)^*\beta\psi\,dx\right|^2\\ &=\left|\int \psi^*\alpha\beta\psi\,dx\right|^2 ); $$\begin{aligned} \frac{\partial}{\partial t}\big|\psi(\bm r,t)\big|^2 &=-\bm \nabla\cdot\bigg\{\frac{\hbar}{2mi}\Big[\psi^*(\bm \nabla\psi)-(\bm \nabla\psi^*)\psi\Big]\bigg\}\\ &=-\bm\nabla\cdot\bm S(\bm r,t) \end{aligned}$$ となる。ここで、&math(\alpha^*=\alpha); を用いた。 右辺は、 が得られる。これはそのまま確率の ''連続の式'' となっており、 $\bm S(\bm r,t)$ が局所的な確率の流れを表わす(単位面積当たり、単位時間当たりの流量) これがピンと来なければ電磁気学で学んだ電荷の保存則を復習せよ。 * エーレンフェストの定理 [#wbe3d8ad] $$ \underbrace{\frac{\partial\rho(\bm r,t)}{\partial t}}_\text{電荷密度の増加速度} =\underbrace{-\bm\nabla\cdot\bm i(\bm r,t)\vphantom{\frac{\partial\rho(\bm r,t)}{\partial t}}}_\text{流れ込む電流量} $$ 初期状態において電子の存在確率があまり広がっておらず、 その広がりに対してポテンシャル &math(V(\bm x,t)); の変化が十分に緩やかであれば、 電子の運動は古典論から予想されるものと等しくなるはずである。 このことを確かめてみよう。 これらは、「確率や電荷は保存しているため、ある領域から流れ出したなら、その分だけその領域での密度が減る」という「確率や電荷の保存則」を表している。 まず、電子の位置座標の時間変化を求める。 電流密度 $\bm i$ に面積をかけることで、その面を通る電流を求められたのと同様に、 確率流密度 $\bm S$ に面積をかけることで、その面を通る単位時間当たりの確率を求められる。 &math( \frac{d}{dt}\langle x\rangle &=\frac{d}{dt}\int \psi^*x\psi\,d^3x\\ &=\int\psi^* x\frac{\PD\psi}{\PD t}+\frac{\PD\psi^*}{\PD t}x\psi\,d^3x\\ ); 上の形はさらに変形できて、 シュレーディンガー方程式により &math(\PD/\PD t); を書き直して、 $$ \begin{aligned} \bm S(\bm r,t)&=\frac{\hbar}{2mi}\big[\psi^*(\bm \nabla\psi)-\{\psi^*(\bm \nabla\psi)\}^*\big] \hspace{10mm}z-z^*=2\,\text{Im}\,z\ \text{の形}\\ &=\frac{\hbar}{2mi}\cdot 2\,\text{Im}\Big[\psi^*(\bm \nabla\psi)\Big]\\ &=\text{Re}\Big[\psi^*\frac{\hbar\bm \nabla}{mi}\psi\Big]\\ &=\text{Re}\Big[\psi^*\frac{\hat{\bm p}}{m}\psi\Big]\\ &=\text{Re}\Big[\psi^*\hat{\bm v}\psi\Big]\\ \end{aligned} $$ &math( \frac{d}{dt}\langle x\rangle &=\int \psi^*x\frac{1}{i\hbar}\Big(-\frac{\hbar^2}{2m}\nabla^2\psi+V\psi\Big)\,d^3x\\ &+\int \frac{1}{-i\hbar}\Big(-\frac{\hbar^2}{2m}\nabla^2\psi^*+V\psi^*\Big)x\psi\,d^3x\\ ); 運動量演算子 $\hat{\bm p}$ を質量 $m$ で割った物は 速度の演算子と考えられるため $\hat{\bm v}=\hat{\bm p}/m$ と置いた。 &math(V); を含む項は互いに打ち消し合って、 確率の流れが速度演算子で書けるのは納得がいくところである。 &math( \frac{d}{dt}\langle x\rangle &=\frac{i\hbar}{2m}\int \psi^*x(\nabla^2\psi)-(\nabla^2\psi^*)x\psi\,d^3x\\ ); 多数の電子が同じ波動関数に従って移動しているような場合には、 この $\bm S$ に素電荷 $e$ を掛けた物が電流密度になるなど、 物理的な意味も大きい。 ここで部分積分により、 覚えておくべき特徴: - 平面波に対しては $\bm S(\bm r,t)=\frac{\hbar\bm k}{m}|\psi|^2$ である - 波動関数が実関数として表わされるとき、括弧内は純虚数となって $S(\bm r,t)=0$ つまり定常状態となる。 &math( \int(\nabla^2\psi^*)x\psi\,d^3x &=-\int (\bm \nabla\psi^*)\cdot\bm \nabla(x\psi)\,d^3x +\underbrace{\int_S(\bm \nabla\psi^*\cdot x\psi)\cdot\bm n\,dS}_{=0}\\ &=\int \psi^*\nabla^2(x\psi)\,d^3x +\underbrace{\int_S(\psi^*\bm\nabla(x\psi))\cdot\bm n\,dS}_{=0}\\ ); ~ [[前の単元 <<<>量子力学Ⅰ/電子の波動方程式]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/箱の中の自由粒子]]~ ~ したがって、 * 質問・コメント [#c0695bd0] &math( \frac{d}{dt}\langle x\rangle &=\frac{i\hbar}{2m}\int \psi^*[x(\nabla^2\psi)-\nabla^2(x\psi)]\,d^3x\\ &=-\frac{i\hbar}{2m}\int \psi^*\frac{\PD}{\PD x}\psi\,d^3x\\ &=\frac{1}{m}\langle p_x\rangle ); #article_kcaptcha **無題 [#g4d8537a] > (&timetag(2023-08-29T03:59:18+09:00, 2023-08-29 (火) 12:59:18);)~ ~ 僕らがいる時点で~ 観測者が現れる確率は0には決してならない~ 途方もない時間笑をかけて~ この世界は物質化するのかもしれない~ を得る。ここで、 // &math( \int \psi^*\nabla^2(x\psi)\,d^3x &=\int \psi^*\Big[(\PD/\PD x)^2+(\PD/\PD y)^2+(\PD/\PD z)^2\Big](x\psi)\,d^3x\\ &=\int \psi^*\Big[(\PD/\PD x)\psi+x(\PD/\PD x)^2\psi+x(\PD/\PD y)^2\psi+x(\PD/\PD z)^2\psi\Big]\,d^3x\\ &=\int \psi^*\Big[(\PD/\PD x)\psi+x\nabla^2\psi\Big]\,d^3x ); #comment_kcaptcha を用いた。一方、運動量の時間変化は、 **この世界は [#g9e02629] > (&timetag(2023-08-29T03:56:24+09:00, 2023-08-29 (火) 12:56:24);)~ ~ この世界は観測した時に物質化し~ 観測者が現れる可能性が0の時~ それは目などの観測デバイスを持たないものからすれば~ 全てがただの波になるということだ~ つまり地球はただのエネルギーの集まりであり~ 物質化する前の観測者がいない状態の話を~ 定義することはナンセンスだ~ &math( \frac{d}{dt}\langle p_x\rangle &=\frac{d}{dt}\int\psi^*\frac{\hbar}{i}\frac{\PD}{\PD x}\psi\,d^3x\\ &=\frac{\hbar}{i}\int\psi^*\frac{\PD}{\PD x}\frac{\PD}{\PD t}\psi+\left(\frac{\PD}{\PD t}\psi^*\right)\frac{\PD}{\PD x}\psi\,d^3x\\ &=\frac{\hbar}{i}\int\psi^*\frac{\PD}{\PD x}\frac{1}{i\hbar}\left(-\frac{\hbar^2}{2m}\nabla^2\psi+V\psi\right)\,d^3x +\frac{\hbar}{i}\int\frac{1}{-i\hbar}\left(-\frac{\hbar^2}{2m}\nabla^2\psi^*+V\psi^*\right)\frac{\PD}{\PD x}\psi\,d^3x\\ ); // ここで、 #comment_kcaptcha &math( \int\psi^*\frac{\PD}{\PD x}\nabla^2\psi\,d^3x &=\int\psi^*\nabla^2\frac{\PD\psi}{\PD x}\,d^3x\\ &=-\int\bm\nabla\psi^*\cdot\bm\nabla\frac{\PD\psi}{\PD x}\,d^3x +\underbrace{\int_S\left(\psi^*\bm\nabla\frac{\PD\psi}{\PD x}\right)\cdot\bm n\,dS}_{=0}\\ &=\int(\nabla^2\psi^*)\frac{\PD\psi}{\PD x}\,d^3x -\underbrace{\int_S\left(\bm\nabla\psi^*\cdot\frac{\PD\psi}{\PD x}\right)\cdot\bm n\,dS}_{=0}\\ &=\int(\nabla^2\psi^*)\frac{\PD\psi}{\PD x}\,d^3x ); より &math(\nabla^2); を含む項は打ち消し合って、 &math( \frac{d}{dt}\langle p_x\rangle &=-\int\psi^*\left(\frac{\PD}{\PD x}(V\psi)-V\frac{\PD}{\PD x}\psi\right)\,d^3x\\ &=-\int\psi^*\frac{\PD V}{\PD x}\psi\,d^3x\\ &=-\left\langle\frac{\PD V}{\PD x}\right\rangle ); となり、&math(\bm x,\bm p); の期待値が古典論の運動方程式 &math(\frac{d\bm x}{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: 174220 (from 2010/06/03),
today: 27,
yesterday: 0