量子力学Ⅰ/固有値と期待値 のバックアップ差分(No.22)
更新- バックアップ一覧
- 現在との差分 を表示
- ソース を表示
- バックアップ を表示
- 量子力学Ⅰ/固有値と期待値 へ行く。
- 追加された行はこの色です。
- 削除された行はこの色です。
[[前の単元 <<<>量子力学Ⅰ/調和振動子]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/物理量の固有関数]]~ * 目次 [#q790942d] [[量子力学Ⅰ]] &mathjax(); &katex(); #contents * 概要 [#ca866d3c] シュレーディンガー方程式を解いて物理量の期待値を求める問題は 線形代数の知識を使って理解すると見通しが良い。 以下、線形代数IIで学んだ関数空間の考え方が量子力学でどのように生かされるかを学ぶ。 * 関数ベクトル・線形演算子 [#a7c68936] #multicolumns ''[数ベクトル]'' &math(\bm a); あるいは &math(\bm a(t)); #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); 「線形演算子の積」は正確に言えば「線形演算子の合成」のことだ。&math(\hat A\hat B); は &math(\hat A\circ\hat B); のことで、「&math(\hat B); を行ってから &math(\hat A); を行う」という合成演算を表す。 #multicolumns(end) * シュレーディンガー方程式・線形演算子の固有値 [#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(\psi); に対して演算子 &math(\hat O); で表される物理量を観測した際の期待値は、 &math(\overline O(t)=(\psi,\hat O\psi)=\iiint\psi^*(\bm r,t)\hat O\psi(\bm r,t)d\bm r); と求められるが、&math(\psi); が &math(\hat O); の固有関数で、&math(\hat O\psi=\lambda\psi); のとき、 期待値は固有値になり、 &math(\overline O=(\psi,\hat O\psi)=(\psi,\lambda\psi)=\lambda\|\psi\|=\lambda); 分散、標準偏差はゼロになる。 &math(\overline{\sigma_O^2} &=\bigl(\psi,(\hat O-\overline O)^2\psi\bigr) =\bigl(\psi,(\hat O-\lambda)^2\psi\bigr) =\bigl(\psi,\ \hat O^2\psi-2\lambda\hat O\psi+\lambda^2\psi\bigr) =\bigl(\psi,\ \lambda^2\psi-2\lambda^2\psi+\lambda^2\psi\bigr) =0 ); すなわち、演算子の固有関数に対して物理量を観測すれば、 常に固有値に等しい値が得られることになる。 - 一般に固有値は複素数値をとるが「物理量」は実数となるべきである - 固有関数ではない波動関数に対しては観測結果は確率的にしか決まらない~ → 確率分布はどのようなものになる? こういったことを考える際に、物理量を表わす演算子がエルミートであることが重要になる。 * エルミート演算子 [#l1d2ced7] #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); その空間に含まれる任意の関数 &math(\psi_1,\psi_2); に対して成り立たなければならない。 #multicolumns(end) #multicolumns ''[エルミート行列]'' &math(H); がエルミート行列とは &math(H^\dagger=H); つまり 任意の &math(\bm a,\bm b); に対して &math((\bm a,H\bm b)=(H\bm a,\bm b)); を満たす行列のこと。 #multicolumns ''[エルミート演算子]'' &math(\hat H); がエルミート演算子とは &math(\hat H^\dagger=\hat H); つまり その関数空間の任意の &math(\psi_1,\psi_2); に対して &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) * エルミート演算子の固有値問題 [#u191050a] エルミート行列の固有値問題は以下の特徴を持っていた。 - 固有値はすべて実数 - 異なる固有値に属する固有ベクトルは互いに直交 - 固有ベクトルを列ベクトルとするユニタリー行列を用いて対角化が可能 この、エルミート行列を対角化するユニタリー行列の各列ベクトルは正規直交系を為し、 なおかつ &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); #multicolumns ''[固有関数による展開]'' 任意の関数を固有関数 &math(\{\varphi_i(\bm r)\}); で展開できる(完全性)。 &math(\varphi(\bm r)=\sum_{i=1}^\infty c_i\varphi_i(\bm r)); #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] ** 2つの固有関数の重ね合わせ [#xe11ca94] ある物理量を表す演算子 &math(\hat O); に対して、 &math(\hat O\varphi_1=\lambda_1\varphi_1); &math(\hat O\varphi_2=\lambda_2\varphi_2); であり、&math(\varphi_1,\varphi_2); は &math((\varphi_i,\varphi_j)=\delta_{ij}); (正規直交) を満たすとする。 この &math(\varphi_1,\varphi_2); を混ぜ合わせて新しい波動関数を作ろう。 &math(\varphi=a\varphi_1+b\varphi_2); この関数のノルムは、 &math( \|\varphi\|^2 &=(a\varphi_1+b\varphi_2, a\varphi_1+b\varphi_2)\\ &=a^*a(\varphi_1,\varphi_1)+a^*b(\varphi_1,\varphi_2)+b^*a(\varphi_2,\varphi_1)+b^*b(\varphi_2,\varphi_2)\\ &=|a|^2+|b|^2\\ ); したがって、&math(|a|^2+|b|^2=1); を満たす実数 &math(a,b); に対して これは正規化された関数となる。(ピタゴラスの定理で理解できることに気づける?) この &math(\varphi); で表される状態に対して &math(\hat O); を観測した際の期待値は、 &math( (\varphi,\hat O\varphi) &=(a\varphi_1+b\varphi_2,\ a\lambda_1\varphi_1+b\lambda_2\varphi_2)\\ &=\lambda_1|a|^2+\lambda_2|b|^2 ); &math(|a|^2+|b|^2=1); の制約の下、&math(|a|^2); を大きくする、すなわち &math(\varphi_1); を混ぜる量を増やすと、期待値は &math(\varphi_1); の固有値 &math(\lambda_1); に近づき、 逆に減らすともう一方の期待値 &math(\lambda_2); に近づく。 ここでは導くことはしないが、このとき実際に観測される値は &math(\lambda_1); か &math(\lambda_2); のどちらか一方で、 - &math(O=\lambda_1); となる確率は &math(|a|^2); - &math(O=\lambda_2); となる確率は &math(|b|^2); となっている。 これを一般化すれば、 ** 一般の場合 [#y88cc546] 任意の状態 $\varphi$ はエルミート演算子 $\hat O$ の固有関数が作る正規直交完全系 $\{\varphi_n\}$ で展開可能であり、その係数を $c_n$ とすれば、 $$ \varphi=\sum_{n=0}^\infty c_n\varphi_n $$ と書ける。この $\varphi$ に対して物理量 $O$ を観測したならば、$O=\lambda_n$ が得られる可能性は $$ P\{O=\lambda_n\}=|c_n|^2 $$ と表せる。ただし、$\{\varphi_n\}$ の正規直交性を用いれば $c_n$ は $$ c_n=\iiint \varphi_n^*(\bm r)\varphi(\bm r)d\bm r $$ として求められ、また、$\varphi$ が規格化されている限り、 $$ \begin{aligned} \iiint \varphi^*(\bm r)\varphi(\bm r)d\bm r &=\sum_{m=0}^\infty\sum_{n=0}^\infty c_m^*c_n\underbrace{\iiint\varphi_m^*(\bm r)\varphi_n(\bm r)d\bm r}_{\delta_{mn}}\\ &=\sum_{n=0}^\infty |c_n|^2\\ &=\sum_{n=0}^\infty P\{O=\lambda_n\}\\ &=1 \end{aligned} $$ いずれかの $\lambda_n$ が観測される確率を全て足すと1になることが保証される。 * 演習:物理量を表わす演算子のエルミート性 [#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)); がエルミート演算子であることを示せ。 (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); はエルミート演算子ではない。 [[>>> 解答例はこちら>@量子力学Ⅰ/固有値と期待値/メモ#f4fd82da]] ** 解説:物理量を表わす演算子のエルミート性 [#c961a723] 上でも述べた通り、量子力学で「観測可能」とされる物理量に対応する演算子は必ずエルミートになる。 ならないとどう困るかはこの後で見る。 逆に、物理量を演算子で表わした際にエルミートにならないような物理量は直接観測可能ではない。 例えば、位置 &math(\hat x); と運動量 &math(\hat p); の積 &math(xp); は直接観測可能ではない。 これは、位置 &math(\hat x); と運動量 &math(\hat p); を同時に、 正確に測定する手段がないという意味で、そのため積を直接は観測できないことになる。 観測可能性については後に ''不確定性原理'' のところで詳しく学ぶ。 一方、&math(\frac{1}{i}\frac{\PD}{\PD x}); がエルミートになるのは波動関数空間に &math(\phi(-\infty)=\phi(\infty)=0); (固定端条件) の境界条件を付いているためであった。 波動関数が有限領域 &math([a,b]); に閉じ込められている場合も含めれば、領域の端 &math(a,b); で、 &math(\phi(a)=\phi(b)=0); (固定端条件) が条件となる。この条件を少し緩めて、 &math(\phi(a)=\phi(b)); (周期的境界条件) としても &math(\frac{1}{i}\frac{\PD}{\PD x}); はエルミートになるが、 さもなければこの演算子がエルミートになることはない。 すなわち、&math(\frac{1}{i}\frac{\PD}{\PD x}); は任意の関数を含む「全関数空間」 においてはエルミートでないが、 「周期的境界条件を満たす」という条件を付けた部分空間においてはエルミートとなる。 このため量子力学の問題では現実では成り立たないにもかかわらず「周期的境界条件」を仮定して問題が設定されることがあるが、戸惑わないように。 * 「観測可能な物理量」の実数性 [#d81900f4] &math(\hat O); がエルミート演算子であれば、固有値 &math(\lambda_i); はすべて実数であるから、 その期待値 &math(\overline O); も必ず実数になることが分かる。 ~ ~ [[前の単元 <<<>量子力学Ⅰ/調和振動子]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/物理量の固有関数]]~ * 質問・コメント [#w546cd07] #article_kcaptcha
Counter: 63206 (from 2010/06/03),
today: 1,
yesterday: 0