量子力学Ⅰ/群速度と波束の崩壊 のバックアップ差分(No.11)

更新


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

#contents

* ある位置に局在する、有限の運動量を持つ波束 [#b0c6b8f8]

-&math(x=x_0); に存在する運動量 &math(p=p_0); の古典粒子

に対応するのは、

-&math(x=x_0); ___付近___に局在する波数 &math(k=k_0=p_0/\hbar); の波動関数

である(波動関数が幅を持たなければ波長や波数を定義できないことに注意せよ)。

このような局在した波動関数を「波束」と呼ぶ。

以下、&math(x_0=0); として波束を作り、その運動を考えよう。

波数 &math(k=k_0); の波動関数は、

#ref(wave-packet-shape.png,right,around,25%);

#ref(wave-packet.png,right,around,50%);

 &math(
\varphi_{k_0}(x)=\frac{1}{\sqrt{2\pi}}e^{ik_0x}
);

ガウス関数を掛けて「波束」にする。

 &math(
\varphi(x)
&=\frac{1}{\sqrt{\sqrt{2\pi}\sigma_{x0}}}e^{-x^2/4\sigma_{x0}^2}\,e^{ik_0x}\\
);

ここでは空間分布が &math(x=0); を中心とする標準偏差 &math(\sigma_{x0}); 
の___正規化された___ガウス関数になるよう係数を設定した。

 &math(
|\varphi(x)|^2=\frac{1}{\sqrt{2\pi}\sigma_{x0}}e^{-x^2/2\sigma_{x0}^2}
);

この関数の波数は &math(k_0); と言って良いだろうか?
(波数が決まるのは平面波のみのはずだが・・・)

波動関数を波数毎の成分に分解するにはフーリエ変換すればいい:

 &math(
\varphi(k)&=\int_{-\infty}^\infty \frac{1}{\sqrt{2\pi}}e^{-ikx}\varphi(x)dx\\
&=\frac{1}{\sqrt{2\pi\sqrt{2\pi}\sigma_{x0}}}\int_{-\infty}^\infty e^{-x^2/4\sigma_{x0}^2}e^{-i(k-k_0)x}dx\\
&=\sqrt{\frac{2\sigma_{x0}}{\pi\sqrt{2\pi}}}\underbrace{\int_{-\infty}^\infty e^{-\{x/2\sigma_{x0}+i\sigma_{x0}(k-k_0)\}^2}\frac{dx}{2\sigma_{x0}}}_{\sqrt \pi}e^{-\sigma_{x0}^2(k-k_0)^2}\\
&=\sqrt{\frac{2\sigma_{x0}}{\sqrt{2\pi}}}e^{-4\sigma_{x0}^2(k-k_0)^2/4}
);

逆フーリエ変換は、

 &math(
\varphi(x)=\int_{-\infty}^\infty \underbrace{\varphi(k)}_{係数}\,
\underbrace{\left(\frac{1}{\sqrt{2\pi}}e^{ikx}\right)}_{波数kの成分}dk
);

であるから、&math(\varphi(x)); に含まれる波数 &math(k); の成分の係数が &math(\varphi(k)); である。

したがって、&math(\varphi); の波数が &math(k); となる確率は、

 &math(
|\varphi(k)|^2&=\frac{2\sigma_{x0}}{\sqrt{2\pi}}e^{-4\sigma_{x0}^2(k-k_0)^2/2}
);

&math(k); の確率分布は &math(k_0); を中心に、&math(\sigma_{k0}=1/2\sigma_{x0}); 
のガウシアンとなる。

「波束」にしたことにより、&math(k=k_0); の成分だけでなく、
おおよそ &math(k=k_0\pm\sigma_{x0}); の成分を含んでしまったことになる。
おおよそ &math(k=k_0\pm\sigma_{k0}); の成分を含んでしまったことになる。
空間的に狭い範囲に局在した波束ほど、幅広い範囲の波数の重ね合わせになることにも注意せよ。

この関数に対して、

 &math(
\sigma_{x}\cdot\sigma_{p}=\sigma_{x}\cdot\hbar\sigma_{k}=\hbar/2
);

は不確定性原理で与えられる最小値をとる。このような波束は「最小波束」と呼ばれる。

ここでは波束 &math(\varphi(x)); の空間分布をガウス関数としたために最小波束が得られた。
- 空間分布がガウス関数でなければ、&math(\sigma_{x}\cdot\sigma_{p}>\hbar/2); となる。
- 空間分布がガウス関数であっても、位相の分布が上記と異なる場合には、やはり最小波束にはならない。

* 自由な波束の運動 [#nbb01d97]

自由な粒子では、

 &math(\varepsilon=\frac{p^2}{2m}); に対応して &math(\hbar\omega_k=\frac{\hbar^2 k^2}{2m});

がその分散関係を与えるのであった。
がその分散関係((&math(\omega); と &math(k); の関係))を与えるのであった。

したがって、波数 &math(k); の成分 &math(\varphi_k(x)=\frac{1}{\sqrt{2\pi}}e^{ikx}); の位相は 
&math(e^{-i\omega_kt}=e^{-i\hbar^2 k^2t/2m}); で回転する。

すなわち、上記最小波束の時間発展は、

 &math(
\psi(x,t)&=\int_{-\infty}^{\infty}\varphi(k)\underbrace{\frac{1}{\sqrt{2\pi}}e^{ikx}e^{-i\omega_kt}}_{\psi_k(x,t)}dk\\
\psi(x,t)&=\int_{-\infty}^{\infty}\varphi(k)\underbrace{\frac{1}{\sqrt{2\pi}}e^{ikx}e^{-i\omega_kt}}_{\psi_k(x,t)=\varphi_k(x)e^{-i\omega_kt}}dk\\
&=\dots\\
&=\sqrt{\frac{1}{\sqrt{2\pi}\sigma_{x0}(1+i\xi t)}}\exp\left[\frac{-x^2/4\sigma_{x0}^2+i(k_0x-\omega_{k0} t)}{1+i\xi t}\right]\\
);

となる(→ [[計算の詳細>量子力学Ⅰ/群速度と波束の崩壊/メモ]])。
となる(→ [[計算の詳細>@量子力学Ⅰ/群速度と波束の崩壊/メモ]])。
ただし、&math(\xi=\frac{\hbar}{2m\sigma_{x0}^2});、&math(\omega_0=\frac{\hbar k_0^2}{2m});

このとき、

 &math(
|\psi(x,t)|^2
&=\frac{1}{\sqrt{2\pi}\sigma_{x0}\sqrt{1+\xi^2 t^2}}
\exp\left[\frac{-\{x-(\hbar k_0/m) t\}^2}{2\sigma_{x0}^2(1+\xi^2t^2)}\right]\\
);

であり、これは 

 &math(x=x_0-(\hbar k_0/m) t);

を中心とする、標準偏差

 &math(\sigma_x=\sigma_{x0}\sqrt{1+\xi^2 t^2}=\sigma_{x0}\sqrt{1+\left(\frac{\hbar}{2m\sigma_{x0}^2}\right)t^2}); 

のガウス関数になっている。

すなわちこの関数は古典論で期待されるのと等しい速度 
&math(v_G=\hbar k_0/m=p_0/m); で進みながら、
徐々に幅が広がり、高さがつぶれていく。

幅の広がる速さは &math(\sigma_{x0}); が小さいほど速いことも注目に値する。

* 分散と波束の崩壊 [#b99dfdf6]

上記の結果より、

- 波束の中心は古典的に期待される一定速度で移動する
- 波束の幅は徐々に広がる

なぜ時間と共に波束が広がるのかを理解するために、
&math(\sigma_{x0}=0.7);, &math(\hbar/2m=1);, &math(k_0=5); とした場合の 
&math(\psi(x,t)); の「実部」を &math(0\le t\le 2); の範囲でプロットした。

 &attachref(wave-packet.gif);

注目すべきは広がった後の波束において、
後縁部に比べて前縁部の波数が大きい(波長が短い)ことである。

なぜか?

平面波一般に対して

 &math(\omega_k=kv_k);

であることと(&math(e^{ik(x-vt)}=e^{i(kx-\omega t)}); を思い出せ)、
自由な電子の分散関係((&math(\omega_k); と &math(k); との関係を分散関係と呼ぶ))が、

 &math(\hbar \omega_k=\frac{\hbar^2k^2}{2m});

であることにより、

 &math(v_k=\frac{\hbar k}{2m});

を得る。

すなわち、自由空間の平面波は「波数の大きな物ほど早く進む」。
これが波束の形状が変化する理由である。

上で見たように波束は「異なる波数を持つ平面波の重ね合わせ」で構成される。
&math(t=0); の最小波束では、
すべての波数成分の位相が位置 &math(x=0); で揃っていた。
その後、波数の小さなものほど遅く、波数の大きなものほど早く移動した結果、
前縁部がより大きな波数を持つ、幅の広がった波束に変化したのである。

&math(\sigma_{k0}=1/2\sigma_{x0}); から分かるとおり、
&math(\sigma_{x0}); が小さい波束ほど、広い範囲の波数成分を含んでいる。
このため空間的に鋭いパルスほど、速やかに幅が広がる。

このように、異なる波数成分が異なる速度を持つ状況を「分散がある」という。
逆に異なる波数成分がすべて同じ速度を持つ場合を「分散がない」という。

分散がない場合の分散関係は

 &math(\omega\propto k);

具体的には上記の通り、

 &math(\omega=vk);  (ただし &math(v); は定数)

である。

* 位相速度と群速度 [#je14e75b]

上の波束の平均的な波数は &math(k=k_0); で、その波数成分の伝播速度は

 &math(v_{k_0}=\frac{\hbar k_0}{2m});

である。このことと対応して、波束の中心部において &math(\mathrm{Re}(\psi)); がゼロになる点(より一般には「等しい位相を持つ点」)が移動する速度は、

 &math(v_{\phi}=\frac{\hbar k_0}{2m}); (「位相速度 (phase velocity)」と呼ぶ)

となる。一方、「波束の中心」が移動する速度は

 &math(v_G=\frac{\hbar k_0}{m}); (「群速度 (group velocity)」 と呼ぶ)((波束は異なる波長を持つ波の一群と見なせることを上で見た。その波の一群が全体として持つ速度が群速度である。))

であるから、自由な粒子に対しては群速度が位相速度の2倍となる。

この違いは群速度が

 &math(v_G(k_0)=\left .\frac{\PD \omega_k}{\PD k}\right|_{k_0}=\frac{\hbar^2k_0}{m}=2v_{k_0});

で与えられるためである(→ [[詳しい導出>量子力学Ⅰ/群速度と波束の崩壊/メモ#b368cfc6]] 結果が上記の計算と一致することを確認せよ)。

したがって、一般に分散がある場合、つまり &math(\omega\not\propto k); の場合には、

 &math(v_{k_0}=\frac{\omega_{k_0}}{v_{k_0}}\ne\left .\frac{\PD\omega_k}{\PD k}\right|_{k_0}=v_G);

となる。

分散がない場合には位相速度と群速度とは一致するから、
波束の概形が崩れないだけでなく、波束中の搬送波の位相も時間に対して変化しない。

 &attachref(wave-packet2.gif);

同様に、真空中や大気中を進む光の分散はほぼゼロと見なせるが、
物質中では分散がゼロでなくなるため、超短光パルスのパルス幅は物質中で徐々に広がる。

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

#article_kcaptcha


Counter: 66745 (from 2010/06/03), today: 4, yesterday: 0