量子力学Ⅰ/固有値と期待値 の履歴(No.34)

更新


前の単元 <<<              量子力学Ⅰ              >>> 次の単元

目次

量子力学Ⅰ

概要

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

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

関数ベクトル・線形演算子

[数ベクトル]

a\bm aa(t)\bm a(t)

[波動関数]

φ(r)\varphi(\bm r)ψ(r,t)\psi(\bm r,t)

ここで、r\bm r はベクトルの成分を指定する指標であり、 tt はベクトルの時間変化を表す。

[ベクトルの成分]

ana_nan(t)a_n(t)  (整数 nn に対応する成分)

[波動関数の成分]

φ(r)\varphi(\bm r)ψ(r,t)\psi(\bm r,t)  (実数ベクトル r\bm r に対応する成分)

[内積]

(a,b)=nanbn\displaystyle(\bm a,\bm b)=\sum_n a_n^*b_n

[積分]

(φ1,φ2)=φ1(r)φ2(r)dr\displaystyle(\varphi_1,\varphi_2)=\iiint \varphi_1^*(\bm r)\varphi_2(\bm r)d\bm r

[ノルムの二乗]

a2=(a,a)\|\bm a\|^2=(\bm a,\bm a)

[全確率密度]

φ(r)φ(r)dr=φ(r)2dr0\displaystyle\iiint \varphi^*(\bm r)\varphi(\bm r)d\bm r=\iiint |\varphi(\bm r)|^2d\bm r\ge 0

内積の定義にある複素共役のおかげで必ず 0\ge 0 になる。 規格化された波動関数のノルムは1である。

[規格化]

ea=1aa\displaystyle\bm e_a=\frac{1}{\|\bm a\|}\bm a

[規格化]

φ2(r)=φ1(r)φ1(r)2dr\displaystyle\varphi_2(\bm r)=\frac{\varphi_1(\bm r)}{\sqrt{\iiint |\varphi_1(\bm r)|^2d\bm r}}

[行列]

b=Aa\bm b=A\bm a

A(αa+βb)=αAa+βAbA(\alpha\bm a+\beta\bm b)=\alpha A\bm a+\beta A\bm b

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

[線型演算子]

g(r)=A^f(r)g(\bm r)=\hat Af(\bm r)

A^(aφ1(r)+bφ2(r))=aA^φ1(r)+bA^φ2(r)\hat A\big(a\varphi_1(\bm r)+b\varphi_2(\bm r)\big) =a\hat A\varphi_1(\bm r)+b\hat A\varphi_2(\bm r)

関数を別の関数に変換する線型変換を、演算子 を左から掛ける形で書く。 すると線形演算子には行列のように分配法則を適用できる。

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

[行列の積の非可換性]

一般に、ABx/=BAxAB\bm x\ne BA\bm x である。

[演算子の非可換性]

一般に、A^B^ψ/=B^A^ψ\hat A\hat B\psi\ne \hat B\hat A\psi である。

例:x^:φ(r)xφ(r), p^x:φ(r)hixφ(r)\hat x:\varphi(\bm r)\to x\varphi(\bm r),\ \hat p_x:\varphi(\bm r)\to \frac{\hbar}{ i}\frac{\partial}{\partial x}\varphi(\bm r) なら

x^p^xφ(r)p^xx^φ(r)=ihφ(r)\hat x\hat p_x\varphi(\bm r)-\hat p_x\hat x\varphi(\bm r)=i\hbar\varphi(\bm r)

あるいは

x^p^xp^xx^=ih\hat x\hat p_x-\hat p_x\hat x=i\hbar

p^xx^\hat p_x\hat x のような「線形演算子の積」は「線形演算子の合成」のことだ。A^B^\hat A\hat BA^B^\hat A\circ\hat B のことで、「B^\hat B を行ってから A^\hat A を行う」という合成演算を表す。

p^xx^=p^xx^\hat p_x\hat x=\hat p_x\circ\hat xは、xxをかけてからxxで微分する演算。

シュレーディンガー方程式・線形演算子の固有値

[ベクトル方程式]

ihta(t)=AH(t)a(t)\displaystyle i\hbar\frac{\partial}{\partial t}\bm a(t)=A_H(t)\bm a(t)

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

ihtψ(r,t)=H^(r^,p^,t)ψ(r,t)\displaystyle i\hbar\frac{\partial}{\partial t}\psi(\bm r,t)=\hat H(\hat{\bm r},\hat{\bm p},t)\psi(\bm r,t)

[固有値方程式]

AHa=εaA_H\bm a=\varepsilon\bm a

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

H^(r^,p^)φ(r)=εφ(r)\hat H(\hat{\bm r},\hat{\bm p})\varphi(\bm r)=\varepsilon\varphi(\bm r)

[AHA_Hの固有値・固有ベクトル]

ε=ε1,ε2,\varepsilon=\varepsilon_1,\varepsilon_2,\dots

a=a1,a2,\bm a=\bm a_1,\bm a_2,\dots

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

ε=ε1,ε2,\varepsilon=\varepsilon_1,\varepsilon_2,\dots

φ(r)=φ1(r),φ2(r),\varphi(\bm r)=\varphi_1(\bm r),\varphi_2(\bm r),\dots

固有関数の物理量は固有値そのものである

波動関数 φ(r)\varphi(\bm r) に対して演算子 O^\hat O で表される物理量を観測した際の期待値・分散は、

Oφ=(φ,O^φ)=φ(r)O^φ(r)dr\begin{aligned}\overline O_\varphi=(\varphi,\hat O\varphi)=\iiint\varphi^*(\bm r)\hat O\varphi(\bm r)d\bm r\end{aligned}

σOφ2=(φ,(O^O)2φ)\begin{aligned}\sigma_{O\varphi}^2 &=\bigl(\varphi,(\hat O-\overline O)^2\varphi\bigr) \end{aligned}

と求められるが、O^φ=λφ\hat O\varphi=\lambda\varphi のとき、

Oφ=(φ,O^φ)=(φ,λφ)=λφ2=λ\begin{aligned}\overline O_\varphi=(\varphi,\hat O\varphi)=(\varphi,\lambda\varphi)=\lambda\|\varphi\|^2=\lambda\end{aligned}

σOφ2=0\begin{aligned}\sigma_{O\varphi}^2&=0 \end{aligned}

期待値は λ\lambda で、分散はゼロになる。

すなわち、演算子の固有関数に対して物理量を観測すれば、 常に固有値に等しい値が得られる。

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

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

エルミート演算子

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

任意の a,bV\bm a,\bm b\in V に対して

(a,Ab)=(Aa,b)(\bm a,A\bm b)=(A^\dagger\bm a,\bm b)

となるとき AA^\daggerAA のエルミート共役と言い、

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

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

任意の φ1,φ2U\varphi_1,\varphi_2\in U に対して

φ1(r)A^φ2(r)dr=(A^φ1(r))φ2(r)dr\displaystyle\iiint \varphi_1^*(\bm r)\hat A\varphi_2(\bm r)d\bm r=\iiint \big(\hat A^\dagger\varphi_1(\bm r)\big)^*\varphi_2(\bm r)d\bm r

となるとき A^\hat A^\daggerA^\hat A のエルミート共役と言う。

[エルミート行列]

H=HH^\dagger=H つまり 任意の a,bV\bm a,\bm b\in V に対して

(a,Hb)=(Ha,b)(\bm a,H\bm b)=(H\bm a,\bm b)

を満たす行列のこと。

[エルミート演算子]

H^=H^\hat H^\dagger=\hat H つまり 任意の φ1,φ2U\varphi_1,\varphi_2\in U に対して

φ1(r)H^φ2(r)dr=(H^φ1(r))φ2(r)dr\displaystyle\iiint \varphi_1^*(\bm r)\hat H\varphi_2(\bm r)d\bm r=\iiint \big(\hat H\varphi_1(\bm r)\big)^*\varphi_2(\bm r)d\bm r

を満たす演算子のこと。

エルミート行列・演算子の固有値問題

エルミート行列・演算子では、

  • 固有値はすべて実数
  • 異なる固有値に属する固有ベクトル・関数は互いに直交
  • 固有ベクトル・関数により正規直交基底・完全系を作れる

という特徴を持つ。→ 線形代数II/固有値問題・固有空間・スペクトル分解#b253f18d

以下では HHH^\hat H をエルミートとする。

[固有値は実数]

Hx=λx,  x/=0H\bm x=\lambda\bm x,\ \ \bm x\ne \bm 0 なら λR\lambda\in\mathbb R

[固有値は実数]

H^φ(r)=λφ(r),  φ/=0\hat H\varphi(\bm r)=\lambda\varphi(\bm r),\ \ \varphi\ne 0 なら λR\lambda\in\mathbb R

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

nn 個のベクトル {ei}\{\bm e_i\} を、

Hei=λieiH\bm e_i=\lambda_i\bm e_i

かつ

(ei,ej)=δij(正規直交)(\bm e_i,\bm e_j)=\delta_{ij}\hspace{1cm} \text{(正規直交)}

となるように選べる。

それらを並べたユニタリ行列 U=(e1 e2  en)U=\Big(\bm e_1\ \bm e_2\ \dots\ \bm e_n\Big) により、UHU=ΛU^\dagger HU=\Lambda は対角行列になる。

これは eiHej=(ei,Hej)=(ei,λjej)=λiδij\bm e_i^\dagger H\bm e_j=(\bm e_i,H\bm e_j)=(\bm e_i,\lambda_j\bm e_j)=\lambda_i\delta_{ij} に対応する。

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

完全系 {φi(r)}\{\varphi_i(\bm r)\} を、

H^φi(r)=λiφi(r)\hat H\varphi_i(\bm r)=\lambda_i\varphi_i(\bm r)

かつ

φi(r)φj(r)dr=δij(正規直交)\displaystyle\iiint\varphi_i^*(\bm r)\varphi_j(\bm r)\,d\bm r=\delta_{ij}\hspace{1cm} \text{(正規直交)}

となるように選べる。

それらに対して、

φi(r)H^φj(r)dr=λjδij\displaystyle\iiint\varphi_i^*(\bm r)\hat H\varphi_j(\bm r)\,d\bm r=\lambda_j\delta_{ij}

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

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

a=i=1ciei\displaystyle\bm a=\sum_{i=1}^\infty c_i\bm e_i

ただし ci=(ei,a)c_i=(\bm e_i,\bm a)

[固有関数による展開]

任意の関数を固有関数 {φi(r)}\{\varphi_i(\bm r)\} で展開可(完全性)

φ(r)=i=1ciφi(r)\displaystyle\varphi(\bm r)=\sum_{i=1}^\infty c_i\varphi_i(\bm r)

ただし ci=φi(r)φ(r)dr\displaystyle c_i=\iiint\varphi_i^*(\bm r)\varphi(\bm r)\,d\bm r

固有関数でない一般の関数に対する測定値

2つの固有関数の重ね合わせ

正規直交する2つの固有関数、

O^φ1=λ1φ1\begin{aligned}\hat O\varphi_1=\lambda_1\varphi_1\end{aligned}

O^φ2=λ2φ2\begin{aligned}\hat O\varphi_2=\lambda_2\varphi_2\end{aligned}

を混ぜ合わせて新しい波動関数を作る。

φ=aφ1+bφ2\begin{aligned}\varphi=a\varphi_1+b\varphi_2\end{aligned}

この関数のノルムは、

φ2=(aφ1+bφ2,aφ1+bφ2)=aa(φ1,φ1)+ab(φ1,φ2)+ba(φ2,φ1)+bb(φ2,φ2)=a2+b2\begin{aligned} \|\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\\ \end{aligned}

したがって、a2+b2=1|a|^2+|b|^2=1 を満たす実数 a,ba,b に対して これは正規化された関数となる。(ピタゴラスの定理で理解できることに気づける?)

この φ\varphi で表される状態に対して O^\hat O を観測した際の期待値は、

(φ,O^φ)=(aφ1+bφ2, aλ1φ1+bλ2φ2)=λ1a2+λ2b2\begin{aligned} (\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 \end{aligned}

a2+b2=1|a|^2+|b|^2=1 の制約の下、a2|a|^2 を大きくする、すなわち φ1\varphi_1 を混ぜる量を増やすと、期待値は φ1\varphi_1 の固有値 λ1\lambda_1 に近づき、 逆に減らすともう一方の期待値 λ2\lambda_2 に近づく。

このとき実際に観測される値は λ1\lambda_1λ2\lambda_2 のどちらか一方で、

  • O=λ1O=\lambda_1 となる確率は a2|a|^2
  • O=λ2O=\lambda_2 となる確率は b2|b|^2
  • 両方を合わせると確率は a2+b2=1|a|^2+|b|^2=1

となっている。

発展: O=λO=\lambda' となる確率は、

P{O=λ}(λ)={1λ=λ0λ/=λ P_{\{O=\lambda'\}}(\lambda)=\begin{cases} 1&\lambda=\lambda'\\ 0&\lambda\ne\lambda'\\ \end{cases}

の期待値により与えられるので、上記の内容を各自確かめてみよ。

一般の場合

任意の状態 φ\varphi はエルミート演算子 O^\hat O の固有関数が作る正規直交完全系 {φn}\{\varphi_n\} で展開可能であり、その係数を cnc_n とすれば、

φ=n=0cnφn \varphi=\sum_{n=0}^\infty c_n\varphi_n

と書ける。この φ\varphi に対して物理量 OO を観測したならば、O=λnO=\lambda_n が得られる可能性は

P{O=λn}=cn2 P\{O=\lambda_n\}=|c_n|^2

と表せる。ただし、{φn}\{\varphi_n\} の正規直交性を用いれば cnc_n

cn=φn(r)φ(r)dr c_n=\iiint \varphi_n^*(\bm r)\varphi(\bm r)d\bm r

として求められ、また、φ\varphi が規格化されている限り、

φ(r)φ(r)dr=m=0n=0cmcnφm(r)φn(r)drδmn=n=0cn2=n=0P{O=λn}=1 \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}

いずれかの λn\lambda_n が観測される確率を全て足すと1になることが保証される。

エルミート演算子の持つ線形代数的な性質により、量子力学的な期待値の理論が非常にうまく構成されることに注目せよ。

演習:物理量を表わす演算子のエルミート性

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

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

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

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

(3) 演算子 p^:φ(x)hiddxφ(x)\displaystyle\hat p:\varphi(x)\mapsto \frac{\hbar}{i}\frac{d}{dx}\varphi(x) がエルミート演算子であることを示せ。

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

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

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

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

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

>>> 解答例はこちらScan me!

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

上でも述べた通り、量子力学で「観測可能」とされる物理量に対応する演算子は必ずエルミートになる。 エルミートにならないと、上で見た以下の性質を利用できない。

  • 固有値は実数になる
  • 任意の関数を固有状態の一次結合として表せる
  • その係数から観測値の確率分布が分かる

古典論との比較から言うと、 物理量を演算子で表わした際にエルミートにならないような物理量は直接観測可能ではない。 例えば、位置 x^\hat x と運動量 p^\hat p の積 xpxp は直接観測可能ではない。

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

観測できないことと、そもそもその値が定まるような状態(固有状態)が存在しないこととは 本来異なるのであるがこのあたりについては後に 不確定性原理 のところで詳しく学ぶ。

一方、p^x=1ix\displaystyle\hat p_x=\frac{1}{i}\frac{\partial}{\partial x} がエルミートになるのは波動関数空間に

ϕ()=ϕ()=0  (固定端条件)\begin{aligned}\phi(-\infty)=\phi(\infty)=0\text{  (固定端条件)}\end{aligned}

の境界条件を付けたことによるものであった。 波動関数が有限領域 [a,b][a,b] に閉じ込められている場合も含めれば、領域の端 a,ba,b で、

ϕ(a)=ϕ(b)=0  (固定端条件)\begin{aligned}\phi(a)=\phi(b)=0\text{  (固定端条件)}\end{aligned}

が条件となる。この条件を少し緩めて、

ϕ(a)=ϕ(b)  (周期的境界条件)\begin{aligned}\phi(a)=\phi(b)\text{  (周期的境界条件)}\end{aligned}

としても p^x\hat p_x はエルミートになるが、 さもなければこの演算子がエルミートになることはない。

すなわち、p^x\hat p_x は任意の関数を含む「全関数空間」 においてはエルミートでないが、 「周期的境界条件を満たす」という条件を付けた部分空間においてはエルミートとなる。

このため量子力学の問題では現実では成り立たないにもかかわらず「周期的境界条件」を仮定して問題が設定されることがあるが、これは数学的に扱いやすく、なおかつ p^x\hat p_x がエルミートとなるような関数空間を定義することを目的としたものであるから戸惑わないように。



前の単元 <<<              量子力学Ⅰ              >>> 次の単元

質問・コメント




ハミルトニアン演算子について

()

ハミルトニアンはどのような変換を表す演算子と考えられますか。
また、ハミルトニアンを波動関数に作用させることは、何を意味しますか。

  • ハミルトニアンを使ってエネルギーの期待値を求める方法について見直してみるのがよいと思います。 -- 武内(管理人)

波動関数が固有関数である場合とその重ね合わせである場合

()

本文に書かれている内容の解釈は以下のようで正しいですか。

1. 波動関数がある物理量に対応する演算子の固有関数となっている場合
・固有値は存在する。
さらに、
・期待値は固有値と等しい。
・期待値の分散はゼロである。
したがって、物理量の測定値は常に固有値と等しくなる。

2. 波動関数がある物理量に対応する演算子の固有関数の重ね合わせとなっている場合
・固有値は存在しない。
さらに、
・期待値は固有値の加重平均となる。
・期待値の分散は非ゼロである。
このとき、物理量の測定値は常にいずれかの固有値となり、その中間の値とはならない。

ここで、二つの疑問があります。

一つは、1.の場合より、正準共役な物理量に対応する演算子が可換でないことは自明といえますか。なぜなら、演算子が可換であるとき、同時固有関数が存在し、二つの物理量に固有値が存在することは不確定性原理に反するように思うからです。

もう一つは、2.の場合に測定値が常にいずれかの固有値となり、その中間の値とならないことはどのようにして示されますか。

  • 解釈の部分、固有値は演算子に固有のもので、波動関数に属するものではないため「固有値は存在(する/しない)」の部分に違和感を覚えました。その他は問題ないと思います。1つ目の疑問について、自明ではないかもしれませんが、非可換なのは事実です。正準共役な物理量の交換関係について調べてみるのが良いと思います。2つ目の疑問について、O=λ1O=\lambda_1 のときに 11 を取り、O/=λ1O\ne\lambda_1 のときにゼロとなるような関数の期待値を求めることで示せるように思います。 -- 武内(管理人)
  • 2つ目の疑問を抱く人は他にもいると考えて本文を改訂しました。質問ありがとうございました。#xe11ca94 の節の最後の部分です。 -- 武内(管理人)

固有値の意味

()

量子力学では、波動関数に演算子を作用させることで物理量の固有値を得ることを学びました。そのようにして得られる固有値はどのような意味を持っているのでしょうか?
抽象的な質問ですが、よろしくお願いします。

固有関数の物理量は固有値そのものである の項について

()

第1式の被積分関数に含まれる関数はψでなく、φではないでしょうか?
また、演算子Oの平均に関する第3式で内積を求める関数もψではなく、φではないでしょうか?

  • ありがとうございます、ご指摘の通りでした。本文を修正いたしました。 -- 武内(管理人)

Counter: 71622 (from 2010/06/03), today: 3, yesterday: 12