量子力学Ⅰ/固有値と期待値 のバックアップ差分(No.7)
更新- バックアップ一覧
- 現在との差分 を表示
- ソース を表示
- バックアップ を表示
- 量子力学Ⅰ/固有値と期待値 へ行く。
- 追加された行はこの色です。
- 削除された行はこの色です。
[[量子力学Ⅰ]] &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(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(\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); において、「観測可能な物理量」を表わす演算子はすべてエルミート演算子になる。 現実的な問題では常にこの境界条件は満たされるため、このことは非常に重要である。 (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(\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(E\bm a=A_H\bm a); #multicolumns ''[時間に依存しないシュレーディンガー方程式]'' &math(E\varphi(\bm r)=\hat H(\hat{\bm r},\hat{\bm p})\varphi(\bm r)); #multicolumns(end) #multicolumns ''[$A_H$の固有値・固有関数]'' &math(E=E_1,E_2,\dots); &math(\bm a=\bm a_1,\bm a_2,\dots); #multicolumns ''[エネルギー固有値・固有関数]'' &math(E=E_1,E_2,\dots); &math(\varphi(\bm r)=\varphi_1(\bm r),\varphi_2(\bm r),\dots); #multicolumns(end) * エルミート演算子の固有値問題 [#u191050a] 時間に依存しないシュレーディンガー方程式はそのままハミルトニアン演算子に対する固有値問題であったが、 これと同様に任意の「観測可能な物理量」に対応するエルミート演算子に対して固有値問題を考えることができる。 これと同様に任意の「観測可能な物理量」に対応するエルミート演算子に対して固有値問題を考えることができる。(位置演算子、運動量演算子、角運動量演算子・・・) エルミート行列では、 - 固有値はすべて実数 - 異なる固有値に属する固有ベクトルは互いに直交 - ユニタリー行列により対角化が可能 エルミート行列を対角化するユニタリー行列の各列ベクトルは 正規直交系を為し、なおかつ &math(n); 次元空間を張るため、 正規直交基底となるのであった。 であった。 さらに、 エルミート行列を対角化するユニタリー行列の各列ベクトルは正規直交系を為し、 なおかつ &math(n); 次元空間を張るため、''正規直交基底'' をなすのであった。 同様に、エルミート演算子では、 - 固有値はすべて実数 - 異なる固有値に属する固有関数は互いに直交 - 固有関数により正規直交完全系を作れる ** 固有関数による正規直交完全系の例 [#x9fabaa9] 以下では &math(A); や &math(\hat A); をエルミートであるものとする。 無限大の井戸ポテンシャルに対する &math(\varphi_n(x)=\sqrt{\frac{\,2\,}{a}}\sin\Big(\frac{n\pi}{a} x\Big)); は正規直交完全系を為す 調和振動子に対する &math(\varphi_n(x)=\sqrt{\frac{1}{2^nn!\sqrt\pi}\sqrt{\frac{\hbar}{m\omega}}}H_n\Big(\sqrt{\frac{m\omega}{\hbar}}x\Big)\exp\Big(-\frac{1}{2}\frac{m\omega}{\hbar}x^2\Big)); は正規直交完全系を為す * エルミート演算子の期待値 [#w21362f1] ここでは &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\bm \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)=\lambda_j\delta_{ij}=\lambda_i\delta_{ij}); これは &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(\{\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); で与えられる。 #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) * 「観測可能な物理量」の実数性 [#d81900f4] &math(\hat A); がエルミート演算子であれば、固有値 &math(\lambda_i); はすべて実数であるから、 その期待値 &math(\overline A); も必ず実数になることが分かる。 * 固有関数の物理量は固有値そのものである [#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); が得られることになる。 固有関数以外では物理量が確率的にしか決まらないことと対比させて理解せよ。
Counter: 63811 (from 2010/06/03),
today: 38,
yesterday: 0