エーレンフェストの定理 のバックアップの現在との差分(No.1)

更新


  • 追加された行はこの色です。
  • 削除された行はこの色です。
[[量子力学I/波動関数の解釈]]
[[前の単元 <<<>量子力学Ⅰ/物理量の固有関数]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/不確定性原理]]~

* エーレンフェストの定理 [#wbe3d8ad]
&katex();

* 物理量期待値の時間変化 [#j882b192]

時間を含むシュレーディンガー方程式は、波動関数の時間発展をそのまま記述したものになっていた。

$$i\hbar\frac{\partial}{\partial t}\psi=\hat H\psi$$

つまり、

$$\psi(\bm r,t+\Delta t)=\underbrace{\psi(\bm r,t)+\frac{1}{i\hbar}\hat H\psi(\bm r,t)\cdot\Delta t}_{\psi(\bm r,t)\text{のみから求まる}}$$

同様の考えで、ある物理量を表わす演算子を &math(\hat O); として、
この期待値 &math(\langle \hat O\rangle); の時間変化を考える。

$$\begin{aligned}
\frac{d}{dt}\langle\hat O\rangle
&=\frac{d}{dt}\iiint \psi^*\hat O\psi\,d\bm r\\
&=\iiint \frac{\partial}{\partial t}\Big(\psi^*\hat O\psi\Big)\,d\bm r\\
&=\iiint\bigg[\psi^* \hat O\frac{\partial\psi}{\partial t}+\psi^*\underbrace{\frac{\partial\hat O}{\partial t}}_{=\,0}\psi+\frac{\partial\psi^*}{\partial t}\hat O\psi\bigg]d\bm r\\
\end{aligned}$$

ここで $\frac{\PD\hat O}{\PD t}=0$ として真ん中の項は落とす(&math(O); が顕わに時間に依存している場合にはこの項も残る)。

シュレーディンガー方程式:&math(\frac{\PD\psi}{\PD t}=\frac{1}{i\hbar}\hat H\psi); により、

$$\begin{aligned}
\frac{d}{dt}\langle\hat O\rangle
&=\iiint\bigg[\psi^* \hat O\Big(\frac{1}{i\hbar}\hat H\Big)\psi+\Big(\frac{1}{i\hbar}\hat H\psi\Big)^*\hat O\psi\bigg]d\bm r\\
&=\iiint\bigg[\frac{1}{i\hbar}\psi^* \hat O\hat H\psi+\frac{1}{-i\hbar}\Big(\hat H\psi\Big)^*\hat O\psi\bigg]d\bm r\\
&=\frac{1}{i\hbar}\iiint\bigg[\psi^* \hat O\hat H\psi-\Big(\hat H\psi\Big)^*\hat O\psi\bigg]d\bm r\\
&=\frac{1}{i\hbar}\iiint\bigg[\psi^* \hat O\hat H\psi-\psi^*\hat H\hat O\psi\bigg]d\bm r\hspace{1.5cm}\because\ \hat H\ はエルミート\\
&=\frac{1}{i\hbar}\langle \hat O\hat H-\hat H\hat O \rangle
\end{aligned}$$

したがって、期待値の時間変化は演算子とハミルトニアンとの交換関係の期待値で表わされる。

 &math(
\frac{d}{dt}\langle\hat O\rangle
&=\frac{1}{i\hbar}\langle \hat O\hat H-\hat H\hat O \rangle
);

ここから、演算子とハミルトニアンが交換するならば、その物理量は時間に依存しない保存量となることが分かる。


ここに現れた &math(\hat H\hat O-\hat O\hat H); に代表される、
&math(\hat A\hat B-\hat B\hat A); の形を &math(\hat A); と &math(\hat B); との交換子と呼び、
しばしば &math(\hat A\hat B-\hat B\hat A=[\hat A,\hat B]); と書く。
また、交換子で表される &math(\hat A,\hat B); の関係を交換関係と呼ぶ。

演算子の期待値を求める計算にはこれからも交換子が頻出する。

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

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

まず、電子の位置座標の時間変化を求める。
まず、電子の位置座標((求めているのは期待値の変化であるが、ここでは電子はあまり広がっていないと考えているため、期待値はそのまま電子の位置と見なせる。))の時間変化を求める((&math(x,t); は独立なパラメータであるため、&math(\frac{\PD x}{\PD t}=0); である。))。

 &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\\
(1) &math(\hat O=x); と置き、

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

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

 &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\\
);
(2) 交換関係 &math(x\hat p_x-\hat p_xx=i\hbar); を用いて、

&math(V); を含む項は互いに打ち消し合って、
 &math(\hat p_x^2x&=x\hat p_x^2-2i\hbar\hat p_x);

を導け。

(3) (2) および &math(x\hat p_y-\hat p_yx=0, x\hat p_z-\hat p_zx=0); を用いて、

 &math(x\hat p^2-\hat p^2x&=2i\hbar\hat p_x);

を導け。

(4) (1)、(3) を用いて、

 &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\\
&=\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
&=\frac{\left\langle p_x\right\rangle}{m}
);

を得る。2つ目の等号は運動量演算子がエルミートであることを用いた。3つ目の等号は、
を導け。

 &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
);
次に運動量の時間変化について考える。

を用いた。一方、運動量の時間変化は、
(5) &math(\hat p_x); と &math(\hat p^2); が交換することを導け。

 &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\\
);
(6) &math(\hat p_xV&=\frac{\hbar}{i}\frac{\PD V}{\PD x}+V\hat p_x); を導け。

ここで、
(7) &math(\frac{d}{dt}\langle p_x\rangle&=-\left\langle\frac{\PD V}{\PD x}\right\rangle); を導け。

 &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); を含む項は打ち消し合って、
#hr

 &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
);
$y,z$ に対しても同様の結果が得られるため、
$\bm r,\bm p$ の期待値が古典論の運動方程式

となり、&math(\bm r,\bm p); の期待値が古典論の運動方程式
$$\frac{d\bm r}{dt}=\frac{\bm p}{m}$$

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

 &math(\frac{d\bm p}{dt}=-\bm\nabla V);
に対応して、

$$\frac{d\langle\hat{\bm r}\rangle}{dt}=\frac{\langle\hat{\bm p}\rangle}{m}$$

$$\frac{d\langle\hat{\bm p}\rangle}{dt}=-\langle\bm\nabla V\rangle$$

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

巨視的極限に於いてシュレーディンガー方程式が古典論の運動方程式を与えるという
この定理をエーレンフェストの定理と呼ぶ。
シュレーディンガー方程式が巨視的極限(波動関数の広がりのスケールに対してポテンシャルの変化が十分にゆっくりである条件)において古典論の運動方程式に一致する、という この定理をエーレンフェストの定理と呼ぶ。
一方、波動関数の広がりが大きく、その中でポテンシャルが大きく変化してしまう場合には量子力学的効果が強く表れることになる。

上で求めた、

$$
\frac{d}{dt}\langle\hat O\rangle
=\frac{1}{i\hbar}\langle [\hat O,\hat H] \rangle
$$

という式は、演算子自体が時間を含む場合には

$$
\frac{d}{dt}\langle\hat O\rangle
=\frac{1}{i\hbar}\langle [\hat O,\hat H]\rangle+\langle\tfrac{\partial \hat O}{\partial t} \rangle
$$

となるが、この式は [[ハミルトン力学における力学変数 $F$ の時間発展>解析力学/ポアソン括弧式#m6d75ca9]] を表す、

$$
\frac{d}{dt} F=\big\{F,H\big\}+\frac{\partial F}{\partial t}
$$

のような式に現れるポアソン括弧 

$$
\big\{F,H\big\}=\sum_{i=1}^n \Big[\frac{\partial F}{\partial q_i}\frac{\partial H}{\partial p_i}-\frac{\partial F}{\partial p_i}\frac{\partial H}{\partial q_i}\Big]
$$

を交換子 $\tfrac{1}{i\hbar}[\cdot,\cdot]$ で置き換えた形になっており、ハイゼンベルクの運動方程式と呼ばれる。

この対応関係は非常に重要なものであるが、解析力学(ハミルトン力学)に関する最低限の知識がないと
何を言っているのかわからないと思うので、もし足りない人は [[解析力学]] などをテキストとして
解析力学の基本的な部分を見直してほしい。

~
[[前の単元 <<<>量子力学Ⅰ/物理量の固有関数]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/不確定性原理]]~

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

#article_kcaptcha


Counter: 18132 (from 2010/06/03), today: 23, yesterday: 0