量子力学Ⅰ/固有値と期待値 のバックアップソース(No.15)

更新

[[前の単元 <<<>量子力学Ⅰ/調和振動子]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/物理量の固有関数]]~

* 目次 [#q790942d]

[[量子力学Ⅰ]]
&mathjax();

#contents

* 概要 [#ca866d3c]

シュレーディンガー方程式を解いて物理量の期待値を求める問題は
線形代数の知識を使って理解すると見通しが良い。

以下、線形代数IIで学んだ関数空間の考え方が量子力学でどのように生かされるかを学ぶ。

* 確率密度・期待値 [#a7c68936]

#multicolumns
''[数ベクトル]''

&math(\bm a);

#multicolumns
''[波動関数]''

&math(\psi(\bm r,t));

#multicolumns(end)

#multicolumns
''[内積]''

&math((\bm a,\bm b));

#multicolumns
''[積分]''

&math(\iiint \psi_1^*(\bm r,t)\psi_2(\bm r,t)d\bm r);

#multicolumns(end)

#multicolumns
''[ノルムの二乗]''

&math(\|\bm a\|^2=(\bm a,\bm a));

#multicolumns
''[全確率密度]''

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

規格化された波動関数のノルムは1である。
#multicolumns(end)

#multicolumns
''[規格化]''

&math(\bm e_a=\frac{1}{\|\bm a\|}\bm a);

#multicolumns
''[規格化]''

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

#multicolumns(end)

#multicolumns
''[行列]''

&math(\bm b=A\bm a);

&math(A(\alpha\bm a+\beta\bm b)=\alpha A\bm a+\beta A\bm b);

数ベクトルを数ベクトルに変換する任意の線型変換は、行列のかけ算で表現できる。

#multicolumns
''[線型演算子]''

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

&math(\hat A\big(a\psi_1(\bm r,t)+b\psi_2(\bm r,t)\big)
=a\hat A\psi_1(\bm r,t)+b\hat A\psi_2(\bm r,t));

関数を別の関数に変換する線型変換を、''演算子'' を左から掛ける形で書く。

例:ハミルトニアン &math(\hat H=-\frac{\hbar^2}{2m}\nabla^2+V(\bm r));

#multicolumns(end)

#multicolumns
''[行列の積の非可換性]''

一般に、&math(AB\bm x\ne BA\bm x); である。

#multicolumns
''[演算子の非可換性]''

一般に、&math(\hat A\hat B\psi\ne \hat B\hat A\psi); である。

例:&math(\hat A:\psi(x)\to \frac{d}{dx}\psi(x),\ \hat B:\psi(x)\to x\psi(x)); なら

 &math(\frac{d}{dx}(x\psi)=\psi+x\frac{d}{dx}\psi=\Big(1+x\frac{d}{dx}\Big)\psi\ne x\frac{d}{dx}\psi);

#multicolumns(end)

#multicolumns
''[行列のエルミート共役]''

&math((\bm a,A\bm b)=(A^\dagger\bm a,\bm b));

&math(A^\dagger=(A^T)^*);   転置の複素共役

#multicolumns
''[演算子のエルミート共役]''

&math(\iiint \psi_1^*(\bm r,t)\hat A\psi_2(\bm r,t)d\bm r=\iiint \big(\hat A^\dagger\psi_1(\bm r,t)\big)^*\psi_2(\bm r,t)d\bm r);

#multicolumns(end)

#multicolumns
''[エルミート行列]''

&math(H); がエルミート行列とは &math(H^\dagger=H); つまり

&math((\bm a,H\bm b)=(H\bm a,\bm b));

を満たす行列のこと。

#multicolumns
''[エルミート演算子]''

&math(\hat H); がエルミート演算子とは &math(\hat H^\dagger=\hat H); つまり

&math(\iiint \psi_1^*(\bm r,t)\hat H\psi_2(\bm r,t)d\bm r=\iiint \big(\hat H\psi_1(\bm r,t)\big)^*\psi_2(\bm r,t)d\bm r);

を満たす演算子のこと。

#multicolumns(end)

* 演習:物理量を表わす演算子のエルミート性 [#qf83b503]

境界条件 &math(\varphi(-\infty)=\varphi(\infty)=0); を満たす一次元波動関数 &math(\varphi(x)); 
の集合 &math(U); は線形空間と見なせる。~
→ 和と複素数倍に対して境界条件は保存されるから、これらの演算に対して &math(U); は閉じている。

&math(U); において、「観測可能な物理量」を表わす演算子はすべてエルミート演算子になる。
現実的な問題では常にこの境界条件は満たされるため、このことは非常に重要である。

(1) 演算子 &math(\hat f:\varphi(x)\mapsto f(x)\varphi(x)); 
(ただし &math(f(x)); は実数関数つまり &math(f(x)\in\mathbb R);) 
のエルミート共役が &math(\hat f); 自身になること、すなわち &math(\hat f); 
がエルミート演算子であることを示せ。(&math(f(x)=x); と置けば、演算子 
&math(\hat x:\varphi(x)\mapsto x\varphi(x)); もエルミートであることがわかる)

(2) 演算子 &math(\frac{d}{dx}); のエルミート共役が &math(-\frac{d}{dx});
となることを示せ。部分積分を使い、上記の境界条件を用いるとよい。

(3) 演算子 &math(\hat p:\varphi(x)\mapsto \frac{\hbar}{i}\frac{d}{dx}\varphi(x)); 
のエルミート共役が &math(\hat p); 自身になること、すなわち &math(\hat p); 
がエルミート演算子であることを示せ。

(4) エルミート演算子 &math(\hat X); の&ruby(べき){冪}; 
&math(\hat X^n:\varphi(x)\mapsto \underbrace{\hat X\hat X\dots\hat X}_{n}\varphi(x)\rule{0pt}{1.8em}); 
がエルミート演算子となることを示せ。

(5) エルミート演算子 &math(\hat X,\hat Y); の和 
&math(\hat X+\hat Y:\varphi(x)\mapsto \hat X\varphi(x)+\hat Y\varphi(x)); 
がエルミート演算子となることを示せ。

(6) エルミート演算子 &math(\hat X,\hat Y); の積 
&math(\hat X\hat Y:\varphi(x)\mapsto \hat X\big(\hat Y\varphi(x)\big)); は 
&math(\hat X\hat Y-\hat Y\hat X=0); すなわち &math(\hat X); と &math(\hat Y); 
が可換でない限り ''エルミート演算子ではない'' ことを示せ。

(7) ハミルトニアン演算子 
&math(\hat H=-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+V(x)=\frac{\hat p^2}{2m}+V(x)); 
がエルミート演算子であることを示せ。

(8) 上記演算子 &math(\hat x,\hat p); について、
&math(\hat x\hat p-\hat p\hat x=i\hbar); 
であることを示せ。
すなわち &math(\hat p\hat x); や &math(\hat x\hat p); はエルミート演算子ではない。

** 解説:物理量を表わす演算子のエルミート性 [#c961a723]

上記の通り、量子力学で「観測可能」とされる物理量は必ずエルミートになる。

逆に、物理量を演算子で表わした際にエルミートにならないような物理量は観測可能ではない。
例えば、位置 &math(\hat x); と運動量 &math(\hat p); の積 &math(xp); は観測可能ではない。

これは、位置 &math(\hat x); と運動量 &math(\hat p); を同時に、
正確に測定する手段がないという意味で、そのため積を観測できないことになる。

観測可能性については後に ''不確定性原理'' のところで詳しく学ぶ。


* シュレーディンガー方程式・固有値 [#h29e3be1]

#multicolumns

''[ベクトル方程式]''

&math(i\hbar\frac{\PD}{\PD t}\bm a(t)=A_H(t)\bm a(t));

#multicolumns

''[シュレーディンガー方程式]''

&math(i\hbar\frac{\PD}{\PD t}\psi(\bm r,t)=\hat H(\hat{\bm r},\hat{\bm p},t)\psi(\bm r,t));

#multicolumns(end)
#multicolumns

''[固有値方程式]''

&math(\varepsilon\bm a=A_H\bm a);

#multicolumns

''[時間に依存しないシュレーディンガー方程式]''

&math(\varepsilon\varphi(\bm r)=\hat H(\hat{\bm r},\hat{\bm p})\varphi(\bm r));

#multicolumns(end)
#multicolumns

''[$A_H$の固有値・固有ベクトル]''

&math(\varepsilon=\varepsilon_1,\varepsilon_2,\dots);

&math(\bm a=\bm a_1,\bm a_2,\dots);

#multicolumns

''[エネルギー固有値・固有関数]''

&math(\varepsilon=\varepsilon_1,\varepsilon_2,\dots);

&math(\varphi(\bm r)=\varphi_1(\bm r),\varphi_2(\bm r),\dots);

#multicolumns(end)

* 固有関数の物理量は固有値そのものである [#h3f6434b]

期待値は固有値になり、

&math(\overline A
&=\iiint\varphi^*(\bm r)\hat A\varphi(\bm r)\,d\bm r\\
&=\iiint\varphi^*(\bm r)\lambda\varphi(\bm r)\,d\bm r\\
&=\lambda\iiint\varphi^*(\bm r)\varphi(\bm r)\,d\bm r\\
&=\lambda
);

分散、標準偏差はゼロになる。

&math(\overline{\sigma_A^2}
&=\iiint\varphi^*(\bm r)(\hat A-\overline A)^2\varphi(\bm r)\,d\bm r\\
&=\iiint\varphi^*(\bm r)(\hat A-\lambda)^2\varphi(\bm r)\,d\bm r\\
&=\iiint\varphi^*(\bm r)\big(\hat A^2\varphi(\bm r)-2\lambda\hat A\varphi(\bm r)+\lambda^2\varphi(\bm r)\big)\,d\bm r\\
&=\iiint\varphi^*(\bm r)\big(\lambda^2\varphi(\bm r)-2\lambda^2\varphi(\bm r)+\lambda^2\varphi(\bm r)\big)\,d\bm r\\
&=0
);

すなわち、ある物理量を表わす演算子 &math(A); の固有関数 &math(\varphi(\bm r));
に対して物理量 &math(A); を観測すれば、常に固有値 &math(\lambda); が得られることになる。

- 一般に波動関数は複素数値をとるが「物理量」は実数となるべきである
- 固有関数ではない波動関数に対しては観測結果は確率的にしか決まらない~
→ 確率分布はどのようなものになる?

こういったことを考える際に、物理量を表わす演算子がエルミートであることが重要になる。

* エルミート演算子の固有値問題 [#u191050a]

時間に依存しないシュレーディンガー方程式はそのままハミルトニアン演算子に対する固有値問題であったが、
これと同様に任意の「観測可能な物理量」に対応するエルミート演算子に対して固有値問題を考えることができる。(位置演算子、運動量演算子、角運動量演算子・・・)

エルミート行列では、
- 固有値はすべて実数
- 異なる固有値に属する固有ベクトルは互いに直交
- 固有ベクトルを列ベクトルとするユニタリー行列を用いて対角化が可能

であった。

この、エルミート行列を対角化するユニタリー行列の各列ベクトルは正規直交系を為し、
なおかつ &math(n); 次元空間を張るから、&math(n); 次元数ベクトル空間の ''正規直交基底'' 
を作ることになる。

同様に、エルミート演算子では、
- 固有値はすべて実数
- 異なる固有値に属する固有関数は互いに直交
- 固有関数により正規直交完全系を作れる

以下では &math(A); や &math(\hat A); をエルミートとする。

#multicolumns
''[固有値は実数]''

&math(A\bm x=\lambda\bm x); ならば &math(\lambda\in\mathbb R);

#multicolumns
''[固有値は実数]''

&math(\hat A\varphi(\bm r)=\lambda\varphi(\bm r)); ならば &math(\lambda\in\mathbb R);

#multicolumns(end)

#multicolumns
''[固有ベクトルで正規直交系基底を作れる]''

&math(n); 個のベクトル &math(\{\bm e_i\}); を、

&math(A\bm e_i=\lambda_i\bm e_i);

かつ

&math((\bm e_i,\bm e_j)=\delta_{ij});

となるように選ぶことができる。

それらを並べた &math(U=\Big(\bm e_1\ \bm e_2\ \dots\ \bm e_n\Big)); 
により、&math(U^\dagger AU=\Lambda); を対角行列にできる。

これは &math(\bm e_i^\dagger A\bm e_j=(\bm e_i,A\bm e_j)=(\bm e_i,\lambda_j\bm e_j)=\lambda_j\delta_{ij}); 
となることに対応する。

#multicolumns
''[固有関数で正規直交完全系を作れる]''

完全系 &math(\{\varphi_i(\bm r)\}); を、

&math(\hat A\varphi_i(\bm r)=\lambda_i\varphi_i(\bm r));

かつ

&math(\iiint\varphi_i^*(\bm r)\varphi_j(\bm r)\,d\bm r=\delta_{ij});

となるように選ぶことができる。

それらに対して、

&math(\iiint\varphi_i^*(\bm r)\hat A\varphi_j(\bm r)\,d\bm r=\lambda_j\delta_{ij}=\lambda_i\delta_{ij});

#multicolumns(end)

#multicolumns

''[固有ベクトルによる展開]''

任意のベクトルを固有ベクトル &math(\{\bm e_i\}); で展開できる(&math(n); 本の一次独立なベクトルは基底をなすから)。

&math(\bm a=\sum_{i=1}^\infty c_i\bm e_i); 

このとき、

&math((\bm a,A\bm a)&=\sum_{i=1}^\infty\sum_{j=1}^\infty c_i^*c_j(\bm e_i,A\bm e_j)\\
&=\sum_{i=1}^\infty\sum_{j=1}^\infty c_i^*c_j\,\lambda_i\delta_{ij}\\
&=\sum_{i=1}^\infty \lambda_i|c_i|^2
);

#multicolumns
''[固有関数による展開]''

任意の関数を固有関数 &math(\{\varphi_i(\bm r)\}); で展開できる(完全性)。

&math(\varphi(\bm r)=\sum_{i=1}^\infty c_i\varphi_i(\bm r));

このとき、

&math(\overline A&=\iiint \varphi^*(\bm r)\hat A\varphi(\bm r)d\bm r\\
&=\sum_{i=1}^\infty \lambda_i|c_i|^2
);

展開係数 &math(c_i); から期待値を求められる。

より具体的には、測定値が &math(\lambda_i); となる確率は
&math(|c_i|^2); で与えられる。

測定値が &math(\lambda_i); となる確率を、

&math(
\mathrm{Eq}(x,y)=\begin{cases}
1&x=y\\
0&x\ne y
\end{cases});

を用いて次のように表せると考えれば理解しやすい。

&math(
\overline{\mathrm{Eq}(A,\lambda_i)}
&=\sum_{j=1}^\infty \mathrm{Eq}(\lambda_j,\lambda_i)|c_j|^2\\
&=|c_i|^2\\
);

#multicolumns(end)

#multicolumns
''[正規直交系による展開係数]''

上記の &math(c_i); は &math(\bm a); に左から &math(\bm e_i); を掛けて、

&math(c_i
&=(\bm e_i,\bm a)\\
&=\sum_{j=1}^\infty c_j(\bm e_i,\bm e_j) \\
&=c_i\\
);

として求められる。

#multicolumns
''[正規直交系による展開係数]''

上記の &math(c_i); は &math(\varphi(\bm r)); に左から &math(\varphi_i(\bm r)); を掛けて、

&math(c_i=\iiint\varphi_i^*(\bm r)\varphi(\bm r)\,d\bm r);

として求められる。

#multicolumns(end)

* 簡単な例 [#bc173cd9]

 &math(\hat O\varphi_1=\lambda_1\varphi_1);

 &math(\hat O\varphi_2=\lambda_2\varphi_2);

であるとする。

&math(\varphi_1,\varphi_2); が正規直交であれば、
&math(a^2+b^2=1); を満たす実数 &math(a,b); に対して、

 &math(\varphi=a\varphi_1+b\varphi_2);

は正規化された関数となる。なぜなら、

 &math(
\|\varphi\|^2
&=(\varphi,\varphi)\\
&=a^2\|\varphi_1\|^2+ab\cancel{(\varphi_1,\varphi_2)}+ba\cancel{(\varphi_2,\varphi_1)}+b^2\|\varphi_2\|^2\\
&=a^2+b^2\\
&=1
);

この &math(\varphi); から &math(a,b); を取り出すには、

- &math((\varphi_1,\varphi)=a\|\varphi_1\|^2+b\cancel{(\varphi_1,\varphi_2)}=a);
- &math((\varphi_2,\varphi)=a\cancel{(\varphi_2,\varphi_1)}+b\|\varphi_1\|^2=b);

また、

- &math(O=\lambda_1); となる確率は &math(a^2);
- &math(O=\lambda_2); となる確率は &math(b^2);

であるから、

 &math(\langle\hat O\rangle=a^2\lambda_1+b^2\lambda_2);

* 「観測可能な物理量」の実数性 [#d81900f4]

&math(\hat A); がエルミート演算子であれば、固有値 &math(\lambda_i); はすべて実数であるから、
その期待値 &math(\overline A); も必ず実数になることが分かる。

~
~
[[前の単元 <<<>量子力学Ⅰ/調和振動子]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/物理量の固有関数]]~

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

#article_kcaptcha

Counter: 63198 (from 2010/06/03), today: 9, yesterday: 0