このページ「
量子力学」は、
まだ書きかけです。加筆・訂正など、協力いただける皆様の
編集を心からお待ちしております。また、ご意見などがありましたら、お気軽に
トークページへどうぞ。
結晶を成す物質の内部エネルギーおよび熱容量を求めよう。議論を簡単にするため、結晶構造の単位である単位胞 1 つをとり、これを 1 つの分子と見なす。このような取り扱いは結晶の具体的構造によらない普遍的な性質を議論する上で重要である。結晶を構成する分子は互いに相互作用するが、最も主要な効果を及ぼすのは最近接格子点上の分子であり、より遠距離にある分子同士の相互作用はそれらの間に存在する分子同士の相互作用として含めることができる。ここまでで扱うべき問題はかなり簡素になったが、結晶分子の運動がそれほど激しいものでない場合には(気体分子運動論の考えを援用すれば、この状況は結晶内部の温度が極めて低いことに相当する)、各分子は固定された平衡点近傍を振動していると見なすことができる。この場合、分子 1 つ 1 つの運動は独立なものとして取り扱うことができ、平衡点近傍で運動する分子 1 個の周りのポテンシャルエネルギーは
は、その平衡点を原点として以下のように表すことができる。
![{\displaystyle U={\frac {1}{2}}k_{x}x^{2}+{\frac {1}{2}}k_{y}y^{2}+{\frac {1}{2}}k_{z}z^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a85459e4fb648d6d6dee6d3bccd1667e77b9b5e1)
分子の周りのポテンシャルは
の 3 成分に対応する 3 つの自由度を持っている。
また分子の運動エネルギー
は
![{\displaystyle K={\frac {1}{2}}mv_{x}^{2}+{\frac {1}{2}}mv_{y}^{2}+{\frac {1}{2}}mv_{z}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/68e6a35d61f1dabb0f89e96af54540fd91ab8e5f)
となって
の 3 つの速度成分に対応する 3 つの自由度を持っている。これらの運動エネルギーとポテンシャルエネルギーの和は今、熱振動をする分子 1 個が持つ全エネルギーに対応し、分子のエネルギーの自由度は合わせて 6 と数えることができる。なぜならこのエネルギーは 3 次元空間上を運動する粒子の位置と速度の 6 つの独立変数
によって決定されるからである。
古典的な統計力学において、平衡状態ではエネルギー等分配の法則が成り立つことから、独立に振動する結晶分子からなる系について、自由度 1 つにつき
のエネルギーが分配され、系全体のエネルギー
との間に
![{\displaystyle E=N\times 6\times {\frac {kT}{2}}=3NkT}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c88e208dbd397fff044f4e2e79bb42bd33db7d03)
という関係が成り立つ。ここで
は結晶内部に含まれる結晶分子の数であり、また
はボルツマン定数、
は熱力学温度である(以下、温度とは熱力学温度のことを指すとする)。ボルツマン定数
とアヴォガドロ定数
の積は気体定数
を与える。
![{\displaystyle k={\frac {R}{N_{\mathrm {A} }}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/173f9082237513ab8e63de1386cc84cefb1b42b2)
結晶分子の個数
をアヴォガドロ定数を用いて物質量
に置き換えれば、上述の関係は気体定数を使って以下のように書き直すことができる。
![{\displaystyle E=3NkT=3nN_{\mathrm {A} }{\frac {R}{N_{\mathrm {A} }}}T=3nRT.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/194aaeb7064a68cb1d280eb9e58dada6383c94b8)
気体定数を用いた形式では分子数が現れず、代わりに物質量という量が定義されることに注意しよう。ボルツマン定数を基本定数とする立場では単なる置き換えに過ぎないが、気体定数を基本定数とする場合、ボルツマン定数を用いた形式を与えるには分子の存在をあからさまに認める必要がある。
結晶の1モル当たりの熱容量
は、温度変化に対するエネルギーの増減の割合を全体の物質量で割ったものに相当するから、
![{\displaystyle C={\frac {1}{n}}{\frac {\partial E}{\partial T}}=3R}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1eaf7d8ce82e3fff12e4327f875095b9777fb156)
となる。これは常温 (
) での結晶の比熱の測定値に一致する。この比熱は温度依存性がなく、常温の固体のモル比熱がほとんど一定であることを示す。固体のモル比熱が常温で一定の値を取るという法則はデュロン=プティの法則 (Dulong-Petit law) と呼ばれる。デュロンとプティはこの法則が多くの物質について良い精度で成り立つことを実験的に発見した人物である。
デュロン=プティの法則が成り立つような系について、常温より遥かに低温の領域においても比熱が一定であることが予想されるが、実験により低温領域では比熱は 0 に収束することを示唆する結果が得られており、低温領域での比熱の温度依存性および比熱の値はデュロン=プティの法則から外れることが知られている。
仮に振動数が
の調和振動子のエネルギーは
の整数倍
しか取れないとする(ただし
は負でないとする)。結晶内部の
個の分子をそれぞれ振動数
の調和振動子と見なせることを仮定し、全部で
の自由度を持つ 1 次元調和振動子の集まりとする。
そうすると、断熱理想気体でも各分子のエネルギーが衝突などにより変動するように(気体全体の全エネルギーは一定)、固体の各振動子のエネルギーも
という飛び飛びの値を移り変わっているとする。
そして
個の振動子のエネルギーの平均値は、仮に下記のように「ボルツマン因子を使って計算できるはず」だと仮定する(※ ボルツマン因子について分からなければ、記事『高等学校化学Ⅱ/化学反応の速さ』の反応速度論での説明(高校~大学初級レベル)、または記事『統計力学I ミクロカノニカル集合』のスターリングの公式を用いた統計力学モデルによる説明(大学中級~)を参照。統計力学的には他にも、ラグランジュの未定乗数法を用いてボルツマン因子の導入を行う方法もある)。
1個の振動子がエネルギー
をとる確率を
とし、この確率がボルツマン因子に比例するとする。
![{\displaystyle \operatorname {Pr} (n)={\frac {1}{Z}}e^{-{\frac {\varepsilon _{n}}{kT}}}={\frac {1}{Z}}e^{-{\frac {nh\nu }{kT}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a949f76d30c5e5f360a259f212f16ef166dc897)
この関数が通常の意味の確率であるためには、すべてのエネルギー状態についての和が 1 に規格化されている必要があるため、比例係数の
は、
![{\displaystyle Z=\sum _{m=0}^{\infty }e^{-{\frac {\varepsilon _{m}}{kT}}}=\sum _{m=0}^{\infty }e^{-{\frac {mh\nu }{kT}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b2f9b239df34297170f9624e66ba5a60fcffee92)
とならなければならない(なお、このZのような量子統計計算の規格化のための関数のことを「分配係数」または「状態和」という)。このとき確率
は
![{\displaystyle \operatorname {Pr} (n)={\frac {\exp \left(-{\frac {nh\nu }{kT}}\right)}{\sum _{m=0}^{\infty }\exp \left(-{\frac {mh\nu }{kT}}\right)}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c41adc129ad64de5afd567e90d649811cd1abcf3)
となる(
は指数関数)。エネルギーの期待値
は、
![{\displaystyle {\begin{aligned}\langle \varepsilon \rangle &=\sum _{n=0}^{\infty }\left\{\varepsilon _{n}\operatorname {Pr} (n)\right\}\\&=\sum _{n=0}^{\infty }\left\{nh\nu \left({\frac {\exp \left(-{\frac {nh\nu }{kT}}\right)}{\sum _{m=0}^{\infty }\exp \left(-{\frac {mh\nu }{kT}}\right)}}\right)\right\}\\&={\frac {1}{\sum _{m=0}^{\infty }\exp \left(-{\frac {mh\nu }{kT}}\right)}}\sum _{n=0}^{\infty }\left\{nh\nu \exp \left(-{\frac {nh\nu }{kT}}\right)\right\}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e44c3032fea78732a49839bf96ed7ef37373b3b1)
と表すことができる。ここでボルツマン定数と温度の積の逆数を
とし(これは逆温度と呼ばれる)、エネルギーの期待値を逆温度
に関する微分を用いて表せば、
![{\displaystyle Z(\beta )=\sum _{m=0}^{\infty }\exp \left(-{\frac {\varepsilon _{m}}{kT}}\right)=\sum _{m=0}^{\infty }\exp \left(-{\frac {mh\nu }{kT}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e11d496ce2bfdf52163e1a7669f80ae29a20466)
より、
![{\displaystyle {\begin{aligned}\langle \varepsilon \rangle &=-{\frac {1}{Z(\beta )}}{\frac {d}{d\beta }}Z(\beta )\\&=-{\frac {d}{d\beta }}\ln Z(\beta )\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2ddea46b6466f3a7c71bcd4066602589916bac6e)
を得る。ここで具体的に右辺の対数を計算すれば、等比級数の和の公式を用いて、
![{\displaystyle {\begin{aligned}Z(\beta )&=\sum _{m=0}^{\infty }\left(e^{-\beta h\nu }\right)^{n}\\&=\left(1-e^{-\beta h\nu }\right)^{-1}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/356c4e9b49c13ce43824d979eda06a5ee912b19b)
と書き直せるから、結局エネルギーの期待値は
![{\displaystyle {\begin{aligned}\langle \varepsilon \rangle &={\frac {d}{d\beta }}\ln \left(1-e^{-\beta h\nu }\right)\\&=h\nu {\frac {e^{-\beta h\nu }}{1-e^{-\beta h\nu }}}\\&={\frac {h\nu }{e^{\beta h\nu }-1}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4385200374d0d31290a9a0b68cb2767e3e11b3d)
と表すことができる。
前節で得た調和振動子のエネルギーの期待値について、調和振動子のエネルギー量子
に掛かる関数
![{\displaystyle {\frac {1}{e^{\beta h\nu }-1}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/38b98c6496077683d7bb283d434632e54dfe71fe)
をプランク分布と呼ぶ。温度がエネルギー量子の大きさに比べて充分小さい場合、
より
という関係が成り立ち、プランク分布は、
![{\displaystyle {\frac {1}{e^{\beta h\nu }-1}}\approx e^{-\beta h\nu }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a6315bc1f45ac834508bcba041866800a5cbcb7)
という形に漸近する。
このプランク分布を利用して、結晶内部の比熱を得ることを考える。結晶を独立な調和振動子の集まりと見なす最も簡単な場合について、結晶全体の内部エネルギーがそれぞれの調和振動子のエネルギー期待値の和にほとんど等しいことから、
![{\displaystyle E=3\langle \varepsilon \rangle =3N{\frac {h\nu }{e^{\beta h\nu }-1}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c8f496edd616a8adf15798d6e9da2770651d85b)
と表すことができる。この場合、結晶分子に対する比熱容量は、
![{\displaystyle c={\frac {1}{N}}{\frac {dE}{dT}}={\frac {1}{N}}{\frac {d\beta }{dT}}{\frac {dE}{d\beta }}=3k(\beta h\nu )^{2}{\frac {e^{\beta h\nu }}{(e^{\beta h\nu }-1)^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5235160133a3388d36f395585a183c28bccfb4c5)
となる。この比熱の低温領域での振る舞いは、
![{\displaystyle c=3k(\beta h\nu )^{2}{\frac {e^{\beta h\nu }}{(e^{\beta h\nu }-1)^{2}}}=3k{\frac {(\beta h\nu )^{2}}{e^{\beta h\nu }}}\to 0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4aee8accfee1c8bba86b7ae0dbba72d1906c33b5)
であり、0 へ収束するという点で低温領域における固体比熱の振る舞いと合致する。高温領域において(ここでいう高温とは調和振動子のエネルギー量子に対してであり、固体の融点温度に比べれば依然低温である)、比熱は
![{\displaystyle c=3ke^{\beta h\nu }\left({\frac {\beta h\nu }{e^{\beta h\nu }-1}}\right)^{2}\to 3k}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3677e846edee5b14a5149f73148cb1636482108f)
となる。高温領域の比熱について、分子比熱
を定積モル比熱
に直すと、
![{\displaystyle C=N_{\mathrm {A} }c\to 3N_{\mathrm {A} }k=3R}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5efbfc2f63298e724571f9137d0957b6f494b858)
となり、これはデュロン=プティの法則に一致する。つまり、エネルギーの量子化という手順を踏むことで低温領域の温度依存性を再現しつつ、常温ではデュロン=プティの法則に漸近するような分布を得られたことになる。
ここからはある物理的な定数を持つことが量子力学的にどのような意味を持つかについて考える。物理的な定数とは例えば、ある物体の持つ位置や運動量のことである。古典力学ではある物体の物理的な状態は位置、運動量などを指定することによって得ることが出来、これらの間に特別な関係は無かった。これらはそれぞれの値を適当に取ってもよい量であったのである。
量子力学的にもある物体の物理的状態を定める量は存在しており、そのような量を定めることで物体がどのような状態にあるかを指定することが出来る。問題なのは、ある場合においてこれらの間に特殊な関係があらわれ、それらの量を任意に選ぶことが出来なくなることである。重要な例として、ある物体の位置と運動量は同時に定めることが出来ない。このことは、古典的には位置と運動量について、運動量が物体の位置を動かす微小変換に関する保存量になっていることによる。このことについては解析力学のネーターの定理が詳しいので、詳細を知りたければ参照せよ。とりあえず量子力学においては、運動量はけっして単なる観念上の量ではなく、物質波の波長に関わる実在的な量である、という事が重要である。
ここで、解析力学の知識を援用して、ある物体の持つエネルギーEが古典的にハミルトニアンHという量で表わされることを用いる。
![{\displaystyle E=H}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8dd06d98f9eebe48afcda925e490cfb39f21edd4)
ここで、ある物理的な状態の全てが数え上げられたとしてこれらの状態全体で張られるベクトルを取る。通常、ある物体が持つ物理的な状態は無数のエネルギーを持ち、このような操作は不可能に思える。実際このことは量子力学の発展の初期に大きな数学的な問題となった。しかし、現在ではベクトルの内積の取り方などを工夫することで、この様な作業が実際可能であることが示されている。詳しくはw:ヒルベルト空間などを参照。
このように全ての物理的状態が数え上げられたとするとき、それらの状態はあるエネルギーを持った状態として存在する。例えば、ある状態
がエネルギー
を持っていたとする。数学的にはこの様な状態はある行列
を用いて
![{\displaystyle {\hat {H}}\psi _{1}=E_{1}\psi _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/172d4e6ac3991d2373c915f86a719d95811db9d5)
と表わせる。ここで、
は、全ての数え上げられた物理的な状態を1つの基底として持つような行列として考えられている。更に
は、それぞれの物理的状態に対して対角化されており、
![{\displaystyle \psi _{1},\psi _{2},\psi _{3},\cdots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce8070ca6791bbe95870e533a27f1e14cb2125e4)
などの全ての物理的状態に対して対応するエネルギー
,
,
などを返すものとする。
このような行列
は、実際にあるエネルギーを持つ状態としては、古典的な考え方と変化することは無い。なぜなら、
は、古典的に考えてある力学系の中に存在する物体が持つと考えられるエネルギー値を全て持っているものと考えることが出来るからである。
このため、仮に全ての量子的状態がエネルギーという量だけで特定されるのならば、ある力学系が取り得るエネルギーを全て定めることが量子的状態を全て求めることになる。ここまでの議論をより数学的な用語を用いてまとめると、出て来た量で
は全ての物理的な状態によって張られた行列であり物理的な状態を表わす
は、
がかかることによってE倍されるようなベクトルであるので、
の固有ベクトルであると考えられる。このときエネルギーEは、固有値方程式
![{\displaystyle {\hat {H}}\psi =E\psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5b9995521be4f5e023d9cec5087f07f99fb4516)
の固有値である。
ここまでのことで、それぞれの量子的状態はエネルギーというただ1つの量で完全に区別されることを仮定して来た。実際にはこのことは必ずしも正しくなく、ある2つの量子的状態が等しいエネルギーを持っていることがある。この時には、各々の量子的状態は互いに区別することが出来ない。しかし、ある状態が持っている量子数は通常エネルギーだけはなく、位置、運動量や角運動量なども含まれている。このような量も用いてそれらの量を区別することが出来る。このとき、エネルギーの場合と同様位置x、運動量pなども量子的状態によって張られる「行列」のような物となることが予想される(「作用素」ともいう)。ここでは、それぞれの行列について
,
などを用いて行列とそうでない量を区別する。
(*注意
行列でかかれる量をq-number、行列でかかれない量をc-numberと呼ぶことがある)
ここまでで位置xと運動量pが行列でかかれることが分かった。
ここで、以前の主張である、量子論である物体の位置と運動量を同時に定めることが出来ないということを用いる。エネルギーの時には、物体の物理量を定めることは対応する物理量を表わす行列を対角化できることに対応していた。このことを用いると、物体の2つの物理量を同時に定められないという主張は対応する2つの物理量を同時に対角化するような基底の張り替え方が存在しないことに対応する。実際に数学的な計算をすることで
![{\displaystyle AB=BA}](https://wikimedia.org/api/rest_v1/media/math/render/svg/992c8ea49fdd26b491180036c5a4d879fec77442)
を満たす2つの行列A,Bについてはこれらを同時に対角化するような基底が取れることが知られる。このことは位置と運動量を表わす行列について
![{\displaystyle {\hat {x}}{\hat {p}}\neq {\hat {p}}{\hat {x}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a2bb32c2ee15c2bde680f20280e3f79adc69289)
が成り立つことを示している。ここで、この結果をまとめるために、ある2つの行列A,Bに対してそれらの交換子
を、
![{\displaystyle [A,B]=AB-BA}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4a3b93b316dd0b6b0ab2c71e486c901ddfe6e79a)
で定義する。性質
![{\displaystyle [A,B]=-[B,A]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37b04a066443e96c0b27d4b14f7f5d51c7b2b0c3)
に注意。このとき、この記号を用いると位置と運動量に関して
![{\displaystyle [{\hat {x}},{\hat {p}}]\neq 0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bcc6398a8c45548dbd85495c83e6e3eeda113403)
が成り立つ。実際には実験的にも
と
の交換子に関して
![{\displaystyle [{\hat {x}},{\hat {p}}]=i\hbar }](https://wikimedia.org/api/rest_v1/media/math/render/svg/42dbbd0db710385288536bcf4f4a1b7cceb75d9a)
が成り立つことが知られている。ここで、
は、w:プランク定数と呼ばれ単位は
で与えられる。ここでプランク定数はきわめて小さい数であり、xやpがある程度大きい量を持つような運動については上の方程式の右辺は0と等しいものとして扱ってよく、このときにはxとpは交換可能であり、2つの量を同時に定めることが出来る。このことは古典力学では2つの量を同時に定めることが可能であることに対応しており、上のようなx,pの交換関係が古典力学とよく適合していることが分かる。
ここまでで量子的な状態
が
![{\displaystyle {\hat {H}}\psi =E\psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5b9995521be4f5e023d9cec5087f07f99fb4516)
の固有方程式で与えられることが分かり、xとpには
![{\displaystyle [{\hat {x}},{\hat {p}}]=i\hbar }](https://wikimedia.org/api/rest_v1/media/math/render/svg/42dbbd0db710385288536bcf4f4a1b7cceb75d9a)
の交換関係が存在することが分かった。実際にはxとpについて上のような交換関係を満たすような行列を用意することはしばしば困難である。この困難に対応する
数学的手法として、ヒルベルト空間の1つとしてw:関数空間が存在することがあげられる。関数空間とはある変数に関する関数をベクトルとして取る手法である。このとき、ベクトルの和は関数の和で表わされ、ベクトルの内積は適当な範囲での関数の積分を用いる。詳しくはw:関数空間を参照。
このような仕方で行列とベクトルを取るとき、固有状態である
はある変数の関数となり、それにかかる
,
,
などの行列はその関数にかかるオペレーターとして表わされる。オペレーターの取り方は元の関数にある関数をかけたり、元の関数を微分したりと様々だが、上の交換関係を満たすような方法は、例えば、xについて
![{\displaystyle {\hat {x}}\rightarrow x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/443fe4bf09bd289b9a26037cf491a5fa5586f508)
のかけ算を対応させ、pに対して、
![{\displaystyle {\hat {p}}\rightarrow -i\hbar {\frac {\partial {}}{\partial {x}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f6d1cd2213ea78b63e8660d62648ff09999b0a76)
を対応させることがあげられる。一方、
pについて
![{\displaystyle {\hat {p}}\rightarrow p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/599cf1c8cf182fb2830d39de6bf2918e3c5743e0)
のかけ算を対応させ、xに対して、
![{\displaystyle {\hat {x}}\rightarrow i\hbar {\frac {\partial {}}{\partial {p}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/344011c22beaf5786847e546c5fd101fa5f4f5f1)
としても同様の結果が得られることが知られており、xとpが互いに入れ換え可能であるという解析力学の結果に適合している。
仮に自由に動く質量mの粒子を考えたとき
![{\displaystyle H={\frac {p^{2}}{2m}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cbb4e14fba5e74749c608d24ee397a97e4025747)
に対して
![{\displaystyle {\hat {H}}=-{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial x^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a20bf2204c159d4bf3465228de3974958e415ac)
が対応することが分かる。このとき上の固有方程式は
![{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial x^{2}}}\psi (x)=E\psi (x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8d296f2841fc64930afc7d636a733da24cd3ec0e)
となり
に関する2階微分方程式になる。ここでいうEのように未知数が含まれている形の微分方程式を特に固有関数方程式と呼ぶことがある。
ここまでで、量子的な状態を計算する手法が1つ得られた。まず、ある古典的なハミルトニアンを選び、それに対して
![{\displaystyle p\rightarrow -i\hbar {\frac {\partial {}}{\partial {x}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce6db4b2eb2ee9deb4f312e84ddd66f06f3c9479)
の置き換えをする。これについて
![{\displaystyle {\hat {H}}\psi =E\psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5b9995521be4f5e023d9cec5087f07f99fb4516)
の微分方程式を解き、対応する微分方程式を解くことで、量子的状態が計算される。
ここで、最後の固有値方程式を時間に依存しないシュレーディンガー方程式 (time-independent Schrödinger equation) と呼ぶ。後に量子力学における時間発展の方程式も扱うが、この名称はそれと比較してのことである。
1次元の系で質量mの物体が、
で、
![{\displaystyle V(x)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b25f708b02fa074441d65ca42c37166dff295d41)
を満たし、
,
で、
![{\displaystyle V(x)=\infty }](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e65635048517912a5369bc53a3982cb73822e36)
を満たすとする。このときシュレーディンガー方程式を解いて、この系に対する波動関数を求めよ。
後に波動関数が0で無い地点では、物体が見つかる可能性があることを解説する。
ここで、
![{\displaystyle V(x)=\infty }](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e65635048517912a5369bc53a3982cb73822e36)
の地点で粒子が見つかってしまうと、その地点で粒子は無限に大きいエネルギーを持っていることになってしまうが、無限のエネルギーというようなことは考えづらいので、ここでは、
,
で波動関数
は、0となることにする。また、量子論でよく用いられる関数空間の要請により、波動関数は連続であることが必要となる。そのため、
![{\displaystyle \psi (0)=\psi (a)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/98ed245a16881faa12095e8ec913687f07a2c052)
が必要となる。さて、このときのシュレーディンガー方程式は
で、
V(x)=0より、
![{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial x^{2}}}\psi (x)=E\psi (x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8d296f2841fc64930afc7d636a733da24cd3ec0e)
![{\displaystyle {\frac {\partial ^{2}}{\partial x^{2}}}\psi (x)=-{\frac {2mE}{\hbar ^{2}}}\psi (x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92a2ebc900413d98abebe4a590b70821aa79d463)
となることが分かるが、これを満たす
は
![{\displaystyle \psi (x)=A\sin(kx)+B\cos(kx)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1b6dcac3da878a6c075863d8856e5d2749878a1a)
で与えられることが分かる。ただし、A,Bは任意の定数であり、
![{\displaystyle k={\sqrt {\frac {2mE}{\hbar ^{2}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ee019d80a9b6ba527719a3619e6620fcf85c13a)
とする。更に、
![{\displaystyle \psi (0)=\psi (a)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/98ed245a16881faa12095e8ec913687f07a2c052)
を用いると、
![{\displaystyle \psi (0)=B=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/73c9a889762fbdc051f565a2d5991eb2f432fe2f)
よりB = 0が分かり、
![{\displaystyle \psi (a)=A\sin(ka)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f730490e44301bebc48d0ae31d710ccccf1074b)
より、nを任意の整数として、
![{\displaystyle ka=n\pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/0150908f6378ed27213d6fa923bef722bae7bac0)
が得られる。よって、それぞれの整数nに対して波動関数
は
![{\displaystyle \psi _{n}(x)=A\sin(k_{n}x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/682994e920329bb9300726cb6cdbd97a90a32dd1)
で与えられる。ただし、
![{\displaystyle k_{n}={\frac {\pi n}{a}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a5371ed7dc26a0ac1a76d20fe4ae636c3a4cbbf)
である。また、上で波動関数の前の任意定数をそのままにしたが、これは関数空間では関数そのものがベクトルであり、それの定数倍はベクトルの性質を変えないことからそのままに残したものである。ただし、後に分かる通り、波動関数の規格化を考えると、この定数は1通りに決まることが示される。
更に、この系において粒子が持つことが出来るエネルギー
も、
![{\displaystyle k_{n}={\sqrt {\frac {2mE_{n}}{\hbar ^{2}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/545adbbb3e0eb973d7d88d3edd51da065eec28e0)
を解くことで得ることが出来る。答えは、
![{\displaystyle E_{n}={\frac {\hbar ^{2}}{2m}}k_{n}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c75f12d80161eb9148e559402241410d729984be)
![{\displaystyle ={\frac {\hbar ^{2}}{2m}}\{{\frac {\pi n}{a}}\}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69db1f386f02278471c9c5a6436d4681fb79e5e0)
![{\displaystyle ={\frac {\hbar ^{2}}{2m}}{\frac {\pi ^{2}n^{2}}{a^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a8e28ef616cf3c45d4d4e529a5e25ba48ca56836)
となり、ある整数nに対して
![{\displaystyle n^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ac9810bbdafe4a6a8061338db0f74e25b7952620)
に比例するエネルギーを持つようになる。このように量子的な系が取り得るエネルギーのことをエネルギー準位と呼ぶことがある。
上で波動関数を計算する方法を得た。ここでは、波動関数の性質について考える。上ではある量子論的な状態をベクトルと見て、ある量子力学の演算子をそれらによって張られた行列として扱った。ここでは一般的にある量子論的な状態を、それらを代表する量子論的な量iを用いて
![{\displaystyle |i>}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ba71f57c1143863337f82383716dfd3e7e16cf0)
と書く。後に述べられる通り、この記法はブラケット記法と呼ばれ、元はw:ディラックによるものである。ここで、量子論的な状態を定めるiのような量をその状態の量子数と呼ぶ。例えば、無限大のポテンシャルによって束縛されている粒子では整数nが量子数
となっている。一般には量子数は整数のような離散的な量である場合も任意の実数で与えられる連続的な量であることもある。
ここで、ある状態
と、それと異なる状態
を取る。ただし、これらの状態はハミルトニアン演算子の、互いに異なった固有値を持つ固有ベクトルであるとする。ここで、ハミルトニアンの固有値は必ず実数でなければならないことが分かる。なぜなら、そうでないときにはエネルギーが虚数になるような量子論的状態が存在することになってしまうからである。一般に、複素数の行列要素を持っており、しかもその固有値が実数になる行列の種類として、エルミート行列があげられる(エルミート行列については物理数学Iを参照)。ここでは、ハミルトニアンはエルミート行列で与えられるものとする。一般に量子論の演算子は通常エルミート演算子である。
更に、あるエルミート行列に対してその行列は必ず対角化され、その固有ベクトルは互いに直交することが知られている。この結果を用いると、エルミート演算子であるハミルトニアンの固有ベクトルである
と
は、互いに直交することが知られる。更に、それぞれの状態の長さを適切に変更することで、任意の状態
,
についてこれらの内積を
とすることが出来る。
については、物理数学Iを参照。ここで、状態の長さを調整することを量子状態の規格化と呼ぶ。ただし、慣習的に状態
,
の内積は
のように書くことが多い。この記法を用いると、任意の
,
に対して、
![{\displaystyle <i|j>=\delta ={ij}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d3ceb27fc7b9d82524170074a5b9f1bce4703be2)
が成り立つ。ここで、ある状態
とそれに対応する波動関数f(x)の関係を、
![{\displaystyle f(x)=<x|i>}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ca3e7908058cba69674543328519ad2e356f5a2)
で取る。ここで、
は対応する粒子がちょうどxで表わされる点にある状態である。この記法は、関数空間の内積の定義と、上で述べた量子論的状態の内積の定義を整合的にすることが分かる。このことを述べるためにまず、関数空間の内積について説明する。ここでは、一般的に波動関数がある複素関数であるとして考える。関数空間の性質によるとある元f(x),g(x)を関数空間の元としたとき、ある積分
が存在して、
![{\displaystyle \int f^{*}(x)g(x)dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2bcbe579c30dce460a5b70ea51b08f7e69dcd9f7)
を元f(x),g(x)の内積と呼ぶ。ここで、xについての積分の範囲は、
とする。ただし、無限大のポテンシャルがある場合のように、波動関数が0となる範囲については積分しなくてもよい。このときには積分範囲はより狭い範囲になるのである。ここで、上の記法を用いると
![{\displaystyle \int f^{*}(x)g(x)dx=\int dx<i|x><x|j>}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d1555705ec8c34b5da18808e452f02fde72d1119)
![{\displaystyle =<i|j>=\delta _{ij}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67536285de405ca63a546d4351b2861bff21a732)
となる。ここで、
![{\displaystyle \int dx<i|x><x|j>}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6ea85a06e8ee82d3f6cdb107f4ab336376751dcc)
についてはまず、
は、任意のxについてもともと
の状態にあった粒子が、xで表わされる点を通過して
の状態に変化することを表わしている。ここで、上では全てのxについてその結果を足し合わせているので、結局、その結果は、
の状態にあった粒子が、
の状態に変化すること方法の全てをつくしていると考えるのである。上で得た
![{\displaystyle \int |x><x|=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/60032ba72309c850e2adb66dd65ba8aea3939f26)
のような表式はベクトルの完全性と呼ばれ、このあと頻繁にでてくる性質である。特に、エルミート演算子に対しては対応する固有ベクトルが完全性の要請を満たすことが知られており、あるエルミート演算子の固有ベクトル
に対して、
![{\displaystyle \Sigma _{i}|i><i|=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7fdb04bb2b560d785681779d7c994ad3b3bc9ec)
が知られている。しかし、特に対応するベクトルが無限個あるときにはこの性質の数学的な証明は難しい場合が多い。
さて、上のことから分かる通り、
![{\displaystyle \int f^{*}(x)g(x)dx=<i|j>=\delta _{ij}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a21ad43b3b0948535162893756703d6b371dec21)
となって、量子論的ベクトルの正規化と対応させるために、波動関数の長さも、1つに定める必要があることが分かる。この条件は全ての波動関数
に対して、
![{\displaystyle \int |\psi (x)|^{2}dx=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4fbb79d036f670d577dd4891442b18355939c072)
とすることで満たされる。このことを波動関数の正規化と呼ぶ。
ここまでで粒子がどの状態にいるのかを指定する方法が分かった。それぞれのエネルギーの固有状態は
などの表示で表わされ、それらの量はどれも対応する波動関数を持つのである。ただし、これらの量はどれも正規化されていなければならない。次に粒子がある状態にいるときに、粒子が実際にどの位置にいるのかを知る方法を考える。ここでいう位置とは古典的な座標の意味であり、
あるエネルギー固有値を持った状態にいる粒子が古典的に見たときにはどの位置で発見されるのかという意味である。仮に対応するエネルギーの固有状態が偶然位置の演算子に対しても固有ベクトルとなっていたとすると、その状態は位置の演算子に対してただ1つの値を持つため、その状態にある粒子が発見される位置は決定している。一方、仮に対応するエネルギーの固有状態が位置の演算子に対して固有ベクトルとなっていなかったとすると、そのときにその粒子は様々な位置で発見されるように思える。実際実験的な結果はそのとおりであり、ある位置の固有状態でない状態にあるときその物体は位置の演算子が値を取り得る位置全体で見つかる確率がある。そして、実際にどの位置にあるかは実際に観測をしてみるまでは、知ることが出来ないのである。このことは全く不思議な結果であるが、例えば量子論的なヤングの実験などにおいてこの結果は確かに確認されているのである。
ここで、あるエネルギーの固有状態
からある位置に発見されてその位置にあることが確定している状態に移行する過程は、対応する位置をxとすると、
![{\displaystyle <x|i>}](https://wikimedia.org/api/rest_v1/media/math/render/svg/693b4636ac0487ea7b75bea6f3f1dd9ff33f5e26)
で与えられることが予想される。しかし、この値はちょうどある固有状態に対応する波動関数f(x)であった。
![{\displaystyle <x|i>=f(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ec49c1865e8e52f312bf5113b566408ad77ca33)
このことから、波動関数f(x)は対応するエネルギーの固有状態にある粒子がある場所xに発見される位置に見つかる過程について関係していることがわかる。実際には更に、この量の絶対値を2乗した量が、ちょうどこの対応する状態にある粒子がその位置に見つかる確率となっているのである。
![{\displaystyle P(x)=|f(x)|^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4628b3c22969bb298893b7ae8a5b9e1e0d641302)
しかし、この量はちょうど
![{\displaystyle \int dx|f(x)|^{2}=P(x)=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d15f69837d0594910aa05b63b6c1e5c243510106)
として、波動関数の正規化を行なった量に対応するが、このことはP(x)を確率を表わす量として扱うための条件とも適合しているのである。
波動関数f(x)が、
![{\displaystyle f(x)={\frac {1}{{}^{4}{\sqrt {\pi }}}}e^{-x^{2}/2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f9327be8434e71f0d3af73f5fcb0408399118a9c)
で与えられるとする。このとき、ある点xで粒子が発見される確率を計算せよ。また、この波動関数が正しく正規化されていることを示せ。
ある点xで粒子が発見される確率P(x)について、
![{\displaystyle P(x)=|f(x)|^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4628b3c22969bb298893b7ae8a5b9e1e0d641302)
が成り立つことを用いればよい。よって、
![{\displaystyle P(x)=|f(x)|^{2}={\frac {1}{\sqrt {\pi }}}e^{-x^{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/36e3c5359d4c1f938f1153a2a7089374c5591cc3)
が得られる。更に、ガウス積分を用いて
![{\displaystyle \int _{-\infty }^{\infty }e^{-x^{2}}={\sqrt {\pi }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/813068f0c910866d8aad50e5a7ba56163eb7ea84)
を用いると、
![{\displaystyle \int dxP(x)=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d85c8473bab8e0fea98547bfdc6e711a8c87d3a8)
が得られ、正しい正規化がなされていることが分かる。ガウス積分については
物理数学Iを参照。
実際にはある状態
からある状態
に移行する確率が
![{\displaystyle |<b|a>|^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56cc1757819c83e9cd0e3e693c2cf9569b861eb4)
で与えられることはあるエネルギーの固有状態がある位置に移行する場合だけにとどまらず、より広い場合にあてはまる。特に上の場合について
![{\displaystyle <b|a>}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c3411092c8ad78eac8b8ca3475ba74794a6e9dc2)
をaからbへの確率振幅と呼ぶ。波動関数は対応するエネルギーの固有状態からある位置で表わされる状態への確率振幅といえる。
互いに直交する状態
,
,
がある。
このとき、
(I)
(II)
![{\displaystyle {\frac {1}{\sqrt {2}}}(|1>+|2>)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b1d3a8b05622f852d123892eef6b07be8f68174)
(III)
![{\displaystyle {\frac {1}{\sqrt {3}}}(|1>+|2>+|3>)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e371f10ce2b7ec0d9348d092a41c3d6c9658ec9f)
(IV)
で与えられる量子状態と状態
との確率振幅を求めよ。それぞれの状態が正しく正規化されていることを示せ。
与えられた状態と
との内積を取ればよい。それぞれの1,2,3で表わされるそれぞれの状態は互いに直交していることに注意せよ。 正規化されていることを調べるにはそれぞれの状態の大きさが1となっていることを調べればよい。
(I)
確率振幅は
![{\displaystyle <1|1>=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a58db9426ab05bfca72ff1b0500a5858d382f9bb)
となり、正規化も
![{\displaystyle <1|1>=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a58db9426ab05bfca72ff1b0500a5858d382f9bb)
となり正しいことが分かる。
(II)
![{\displaystyle <1|{\frac {1}{\sqrt {2}}}(|1>+|2>)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec33506ecb1583f452394214f2edd5675ac9a249)
![{\displaystyle ={\frac {1}{\sqrt {2}}}(1+0)={\frac {1}{\sqrt {2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6641f421c9ad9185190345a7cd02c5bd261956ee)
となる。正規化については
![{\displaystyle {\frac {1}{\sqrt {2}}}(<1|+<2|){\frac {1}{\sqrt {2}}}(|1>+|2>)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/624c806068de233e356fe100ace1640737e65978)
![{\displaystyle ={\frac {1}{2}}(1+0+0+1)=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/21f1c1b137c194945e0be43154911cffec017bcc)
となり正しいことが分かる。
(III)
![{\displaystyle <1|{\frac {1}{\sqrt {3}}}(|1>+|2>+|3>)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e786e4d899737276017bf86f59a7228b8f818340)
![{\displaystyle ={\frac {1}{\sqrt {3}}}(1+0+0)={\frac {1}{\sqrt {3}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/21b8d5355a68aeafa2aa60d594d9c7536f706d46)
となる。正規化については
![{\displaystyle {\frac {1}{\sqrt {3}}}(<1|+<2|+<3|){\frac {1}{\sqrt {3}}}(|1>+|2>+|3>)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6dcd743fed613e4663c6df24c429ed985ea78ac7)
![{\displaystyle ={\frac {1}{3}}(1+0+0+0+1+0+0+0+1)=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb2683038e6f03678c0183178c79279578454016)
となり正しいことが分かる。
(IV)
確率振幅は
![{\displaystyle <1|2>=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fdaf9eb23e90ec1f3b18ce74a3ee8874c9081269)
となる。正規化は
![{\displaystyle <2|2>=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3f6b296b6062e440c8b7ada5cda0e003b689845f)
となって正しいことが分かる。
ここで、あるエネルギーの固有状態
と、対応する波動関数f(x)に対して
![{\displaystyle <i|x|i>=\int dxx|f(x)|^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/28172d4c86273cd08438301369a87a59ec5c87b1)
がどのような意味を持つかを考える。ここで、
が、対応する粒子がxで見つかる確率を表わしていることを考えると、上の式はxの期待値を表わす式そのものである。そのため、
のようなx演算子の対角成分は、対応する状態に粒子が存在するときの粒子が見つかる位置の期待値となることが分かる。一方、位置演算子の非対角成分はそれほど簡単な解釈は持っていない。ただし、これらの量は量子力学的な摂動などでよく使われる。詳しくは量子力学IIを参照。
実際の物理的な系は常に時間に依存して変化する。このため、量子的な状態も何らかの仕方で時間依存性を持つ必要がある。ここで、量子論的な系の時間依存性を考える前に、特殊相対論において、時間と空間を統一的に扱う方法を得たことを思いだす。上の議論で空間方向の成分に対しては
![{\displaystyle {\vec {p}}_{i}\rightarrow -i\hbar {\frac {\partial {}}{\partial {x_{i}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c161117dfefc2ccbaaabb7e5a60fcb3430fc741c)
(ただし、i=1,2,3。)
のような置き換えをしたことを考え、エネルギーと時間方向に同じ様な関係があることを考えると、上の置き換えに対応して
![{\displaystyle E\rightarrow i\hbar {\frac {\partial {}}{\partial {t}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2c7b5f9b95552ea6658f4be8934f850a82810f65)
のような置き換えが出来ることが予想される。一方、量子論的な系では特殊相対論的な考え方が適用できるのかということは疑問が残る。例えば、量子論ではある物体が存在する位置は観測をする前に原理的に知ることができず、観測をした瞬間に物体の位置が決定することが知られている。しかし、ある一瞬で物体の位置が決定されるのなら、その瞬間にその物体はもともと物体があった場所から非常に速い速度で移動しているように思え、その速度は光速を超えてしまうように思える。この様な事情を考えると、特殊相対論と量子力学をお互いに適合させることは非常に困難に思え、上のような置き換えをする理由は定かでないように思える。しかし、実際にはこの様な困難を乗り越えて上の2つを適合させる方法は既に知られており、その結果を用いるなら確かに上の置き換えは正しい結果を与えることが知られるのである。詳しくは場の量子論を参照。
上の置き換えを古典的な方程式
![{\displaystyle E=H}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8dd06d98f9eebe48afcda925e490cfb39f21edd4)
に対して用いるなら、量子論的な方程式は
![{\displaystyle i\hbar {\frac {\partial {}}{\partial {t}}}={\hat {H}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/519d30c9aac434a66a01e8b5bc771eddb1d3e6a9)
のようになる。ここでは、一般に系の量子状態を張るベクトルを
と書いて、
上の方程式を、
![{\displaystyle i\hbar {\frac {\partial {}}{\partial {t}}}\Psi ={\hat {H}}\Psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/59184e55d6b807b07a8407e5d45e45ecce38a501)
と書き換える。仮に
が、エネルギーEを持つハミルトニアンの固有状態だった
とする。このとき、上の方程式は
![{\displaystyle i\hbar {\frac {\partial {}}{\partial {t}}}\Psi =E\Psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/86d499738848243daba23701633d1c5518b24323)
となり、この式は通常の方法で解くことが出来る。仮にt=0で、
![{\displaystyle \Psi (t=0)=\psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/d050a6e7d777192737c7ac028438a3d3e2bb27f1)
が成り立つとすると、上の式の解は
![{\displaystyle \Psi (t)=e^{-iEt/\hbar }\psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/a2dccfba35d6bf506a5703f759eba9c622d6cb88)
となる。このことによって、ある時刻
において、あるハミルトニアンの固有状態で張られる状態にいた物体が時間的にどの状態に変化するかが分かったことになる。一方、ハミルトニアンの固有状態は時間に依存しないシュレーディンガー方程式によって計算されることから、どの状態がどのような時間発展をするかは時間に依存しないシュレーディンガー方程式を解くことによって求められることが分かる。
また、時間に依存するシュレーディンガー方程式と、時間に依存しないシュレーディンガー方程式は互いに関連している。仮に、ある状態
の時間発展がある定数Eを用いて、
![{\displaystyle \Psi (t)=\psi \cdot e^{-iEt/\hbar }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c8f820c31eaea89ff1467b32c43c3d20ff6e1cb6)
で書かれたとする。この時この
を時間に依存するシュレーディンガー方程式に代入すると、結果は、
![{\displaystyle E\psi ={\hat {H}}\psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/8b67d4d6136cfbff6eac163690303ae947214706)
となり、時間に依存しないシュレーディンガー方程式に等しくなる。
水素原子模型でのシュレーディンガー方程式の解法
[編集]
このページ「
量子力学」は、
まだ書きかけです。加筆・訂正など、協力いただける皆様の
編集を心からお待ちしております。また、ご意見などがありましたら、お気軽に
トークページへどうぞ。
境界条件が微分方程式を解く際に必要になってくる。球座標に変換したシュレーディンガー方程式に、角度の周期的境界条件などを入れると(角度は一周すると元に戻るという条件)、水素原子のシュレーディンガー方程式が解けるようになり、なんと大学化学で習う「主量子数」、「方位量子数」などの「量子数」などが導出される。
1次元井戸型ポテンシャル
![{\displaystyle {\begin{cases}V(x)=\infty (-\infty <x<0)\\V(x)=0(0\leq x\leq L)\\V(x)=\infty (L<x<\infty )\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/74ef650419553cfdfb8a9ebfb7d584c19a439e49)
を考える。このときのシュレーディンガー方程式は
![{\displaystyle E\psi (x)=-{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi (x)}{dx^{2}}}+V(x)\psi (x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c7caf93cdfdb546d256f1d4d0a2a70b42ddae67)
となる。このとき
の領域
では粒子侵入不可なので、この領域における波動関数は
となる。波動関数
は
でそれぞれ連続なので、
となる。
における波動関数を考えると
(
)
(
)
![{\displaystyle \psi (x)=A\cos kx+B\sin kx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1f5e12de52a831104b0ed571efe5397d6b7e915c)
より
(
)
となる。また
となるようにBを求めると
![{\displaystyle B={\sqrt {\frac {2}{L}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c862c91c42ce6db6d35c39b3d1d0a9bcb9f0f7c4)
となり
![{\displaystyle \psi (x)={\sqrt {\frac {2}{L}}}\sin {\frac {n\pi x}{L}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc4411131f1710b0bb5720817d9647c5d1d75a6f)
となる。またこのときのエネルギーEは
![{\displaystyle E={\frac {\hbar ^{2}k^{2}}{2m}}={\frac {\hbar ^{2}n^{2}\pi ^{2}}{2mL^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bfb4aa9b950ab3e2b17f51e8691a47ce4810320a)
となり、とびとびの値をとることが分かる。
1次元階段型ポテンシャル
![{\displaystyle {\begin{cases}V(x)=0(x<0)\\V(x)=V_{0}(0\leq x)\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/199978c9f3baa82ad18623254f10dee8fecbed81)
を考える。
領域
における波動関数をそれぞれ
とする.
それぞれのシュレディンガー方程式は,
![{\displaystyle E\psi _{-}(x)=-{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi _{-}(x)}{dx^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80e49e5f91d0cebdbcfe8862bc236af2d3bded5d)
![{\displaystyle E\psi _{+}(x)=-{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi _{+}(x)}{dx^{2}}}+V_{0}\psi _{+}(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/81fc0fb5cf24b1cfbde32e93d22852ca25114b84)
となる.
(1)
の場合
![{\displaystyle k_{-}={\sqrt {\frac {2mE}{\hbar }}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4b1e0da1be790b659583591ec2c1261a3da35957)
![{\displaystyle k_{+}={\sqrt {\frac {2m(V_{0}-E)}{\hbar }}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4cc9bd209a05994de42627435cc48613406d8cda)
とすると,シュレディンガー方程式は,
![{\displaystyle {\frac {d^{2}\psi _{-}(x)}{dx^{2}}}=-k_{-}^{2}\psi _{-}(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d769b12846461a39d4271fa03bea844e9a542add)
![{\displaystyle {\frac {d^{2}\psi _{+}(x)}{dx^{2}}}=k_{+}^{2}\psi _{+}(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/887d0c310437c352f4ce8ab0f1d55ed1f551eb82)
解は
![{\displaystyle \psi _{-}(x)=C_{-}e^{-ik_{-}x}+C_{+}e^{ik_{-}x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad7edbe3be5e254ec2a30b98c4a40162add0328a)
![{\displaystyle \psi _{+}(x)=Ce^{-k_{+}x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1cecbe76e63c030f6c85975e5b7f34dfa51a3948)
(
の
の項は発散し,規格化条件を満たさない除外する.)
波動関数が
で滑らかである条件から定数を定める.
![{\displaystyle \psi _{-}(0)=\psi _{+}(0)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f5e98f989a322ebbf7a4140bdf87cf739b899459)
![{\displaystyle {\frac {d\psi _{-}(0)}{dx}}={\frac {d\psi _{+}(0)}{dx}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95548e2231c12d3c6394aec5a6d3ee9cba66839a)
より,
![{\displaystyle C_{-}+C_{+}=C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fb456cf786d1ec34ec4ea0125de465ccd28254d9)
![{\displaystyle -ik_{-}C_{-}+ik_{-}C_{+}=-k_{+}C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/108fe87edadf2570e6ddf38d51c171ee147e2fcf)
シュテルン=ゲルラッハの実験
4は古典物理的な予想値(じっさいの実験結果ではない)。
じっさいの実験結果は、5のように、原子線は、上下の2つの位置に分かれる。けっして、4のように、そのあいだの中間の位置には、ほぼ原子線は当たらない。
この5のように、原子線は上下2つに分裂する。
通常の棒磁石を、電子に近づけても、電子の磁気モーメントが小さいので、ローレンツ力以外の力は電子はほとんど受けないのだが、磁石が通常でない場合は別である。
図のように、磁石によって、不連続で急峻な磁場が発生するとき、電子の上側と下側とで、磁場の強さが異なる。そのため、電子全体としては、力を受けることになる(なお、前提として、電子には「スピン」という磁極のような性質がある、という事を前提にしている)。
N極の先端のとがったかたちをした棒磁石と、S極の先端のくぼんだ棒磁石を用意して、とがった、N極と、くぼんだS極の軸を一致させ、この2つの磁極の間隔をせまくした不対磁極をつくる。この不対磁極のすきまに、銀を熱して蒸発させて細孔などから飛び出させた銀の原子線を打ち込むと、図のように、上または下のどちらかの力を受け、上下の2箇所に分裂する。けっして、ななめ方向には移動しない。これは、原子線そのものが磁化をもっていることの実験的証明である。
銀の原子は中性である。また、仮に電離していていたとしてローレンツ力を受けたとすると、ローレンツ力の方向は、図中の横向きになる。
この実験をシュテルン・ゲルラッハの実験という。銀以外にも、水素の原子線やナトリウム原子線でも同様の実験が行われ、原子線が上方向または下方向のどちらかの力を受けることが確認された。
このように原子線が上下に分裂する理由は、原子線が磁化をもっている事のあらわれであるが、その原子の磁化の由来は、そもそも電子が物性として磁気をもっているからである。そして、電子そのものの磁気のことをスピンという。
原子線の標的になっている場所を見ると(図の「5」の場所)、原子線が上または下の2通りの位置に分裂して当たっていることから、電子の「スピン」も2通りの値であることが予想され、他の物理理論から、電子の「スピン」が実際に2通りである事が分かっている。
量子論の基礎法則は数学的にはむしろ単純で、線形代数に他ならない。但し、扱う系によりベクトル空間の次元が無限大になったり関数空間になったりするので、そこからくる複雑さが大きい。ここでは有限次元(2次元!)の線形代数で完全に扱うことができる系、「電子のスピン」を例にとりながら、基礎法則を導入する。
電子は点粒子であり位置という属性を持つが、実はそれだけでなく「自転する棒磁石」が持つような属性も持っている。棒磁石がN極とS極を結ぶ線を軸として自転しているとしよう。そこに磁場をかけると(1)棒磁石はその向きに応じたエネルギーを持ち、また(2)磁場を軸としてコマのようにプリセッション(首振運動)を行う。(1)は磁石であることから起き、(2)は((1)に加えて)角運動量を持つことから起きる。
電子も「向き」を持っており、磁場の中に入れると(1)その向きに応じたエネルギーをもち、かつ(2)その向きが棒磁石の向きと同様のプリセッションを起こす。このような事実を称して「電子はスピンをもつ」といい、「向き」を「スピンの向き」と呼ぶ。スピンの起源(Diracによる)や物性との関係は重要な主題となるのだが、ここでは量子力学の法則を説明する例として用いるだけなので、単に
「電子は位置だけでなく、向きという属性を持つ」
という点にのみ着目する。ダイナミクスの説明を行う時に物理的な意味
「電子を棒磁石とみなした時の向きであり、かつ電子のもつ角運動量の向き」
を使う。なお、角運動量の大きさは一定(hbar/2)であり、変わりうるのはその向きだけである。
この向きを単位ベクトル
で表し、その成分を
と書く。単位ベクトルなので
。普通に考えるとこの成分(3つの実数、動く範囲はそれぞれ-1から1)を指定すればスピンの状態を記述したことになるはずである。ところがスピンは量子力学の法則に従うのでそうはならない。
スピンの向きの測定方法で代表的なのはStern-Gerlach型実験と呼ばれるもの。空間的な勾配を持つ磁場をかけることでスピンの向きに依存した力が電子にかかるようにする。大雑把に理解するには、棒磁石に磁場をかけることを考えればよい。磁場が一様だとN極とS極にかかる力が正反対かつ同じ大きさになるので全体としてキャンセルしてしまう。そこでz方向に勾配を持つ、つまり上に行くほど強くなる磁場をかけるとしよう。すると磁石の向きがz軸方向の成分をもつ(
)ならば上側の極により大きな力がかかるので全体にかかる力が残る。適切に磁場を設定することで、近似的にszに比例した力
が磁石にかかるようにできる。そのような磁場の中に小さい磁石を飛ばすと、磁石の向きのz成分に比例した力がかかり、軌道がそれる。従って軌道が上下にどれだけそれたかを調べれば磁石の
が分かるのである。もちろん同じ議論がz軸方向以外でも成り立つ。まとめると、
・
方向(
は単位ベクトル)の勾配を持つ磁場を使うことで、
の
成分
を測定することができる。
z方向に勾配を持つ磁場にランダムな向きの磁石をたくさんいれると、たまたまsz=1だったものの軌道は大きく上にそれ、sz=-1だったものは大きく下にそれる。そしてその中間、特にsz=0に近いものでは軌道はほとんどそれない。磁石の出口に磁石がきたことを感知するスクリーンをおけば、そこには上から下までほぼまんべんなく磁石がきた後が残るであろう。
量子力学の特徴1:観測値の離散化(「量子化」)
[編集]
これで舞台は整ったので、磁石の代りに向きがばらばらの電子を通してみる。すると、驚くべきことに磁石の場合とは全く異なる結果になる。スクリーンの一番上(sz=1に対応)と一番下(sz=-1に対応)にしか電子がこないのである。あたかもsz=1の電子とsz=-1の電子しかないかのような結果になってしまう。それではと、今度は装置を90度回して、sxに応じて軌道が変わるようにしてみる。すると今度は(同じ源からきた電子なのに)sx=1とsx=-1のものしかないかのような結果になる。つまり一番右と左にしか電子がこない。
これは測定法を変えても成り立つ一般的な結果である。即ち、
・電子スピン
の
方向成分
を測定すると、その結果は必ず1か-1のどちらかにしかならない。本来連続的な値をとるはずの量
の測定結果が、1,-1という離散的な値だけになる。
この、「連続的であるはず物理量の測定値が離散的になる」ことが微視的スケールの物理の大きな特徴である。但し、あらゆる量の測定値が離散的になるわけではない。例えば電子の位置の測定値は連続的な値をとりえるし、エネルギーの測定値は系によって連続的な場合も離散的な場合も、両者が混合する場合もある。量子力学の大きな成果(の一つ)は、物理量の測定値がとりえる値のセット(スペクトルという)を求めるための首尾一貫した計算法として与えたことにある。
szの測定値は1,-1しか取り得ないと述べた。では、例えばスピンがx方向を向いている場合にszを測ると測定値はどうなるだろうか。常識的な答である0は取れない。しかも1,-1のどちらを取るにしても妙である。答は「1,-1を全くランダムに取る。つまり等確率で1,-1が現れる」となる(一回ごとでは常識的な値0は取らないが、多数回実験を行った平均値としては常識的な値0になる)。このように、測定結果が確率的になることが微視的世界の法則のもう一つの大きな特徴。
但し上の実験を行うには「スピンがx方向を向いた状態」を作らなければならない。どうするか。実験装置自身を「フィルタ」として使う。例えば上記Stern-Gerlachの実験ではsz=1の電子は軌道が上に曲げられ、-1のものは下に曲げられる。よって上に来た電子だけを集めれば、それらは全てsz=1である、といえる。実際、もう一つ実験装置を用意し、最初の装置で上に曲げられた電子だけを通してszをもう一度測る。すると、確かに全ての電子でsz=1になるのである。もちろん最初の装置と二つの装置の間に磁場などを掛けてしまうと向きが変わることもあるが、そのような「スピンの向きを変える要因」がなければ、一度目と二度目の測定値は必ず同じ値になる。以上のことを踏まえ次の定義を行う。
「sz=1の状態を取っている」電子とは、実験によりsz=1という測定値を取り、その後向きを変えられていない電子のこととする。
任意の向きへの一般化は自明であろう。任意の単位ベクトル
に対し、
・
方向を向いた電子とは、
の測定で測定値1を取り、その後向きを変えられていない電子である。
この定義が意味をなす理由は、上記の通りこの定義に従う電子でもう一度
を測定すれば必ず1を取る、という実験的裏づけがあるためである。但しこの定義に従う電子で
と垂直な成分を測っても測定結果は決して0にはならないことを忘れてはならない。
この定義を使うと、上記のランダム性を実験で確かめるには、まずStern-Gerachの装置を横に倒してsxの値に応じて軌道が曲げられるようにする。そして、右(sx=1に対応)に来た電子だけを第2の装置に通してszを測る。すると、半分はsz=1、残り半分はsz=-1となるのである。
より一般的な場合、つまり
とz軸との角度が任意の角度A(
) の場合を述べよう。z方向を向いた電子スピンの
方向の成分
を測るとする。すると、
測定値が1になる確率は
、-1になる確率は
となる。(特に
の場合には上記の結果に帰着することに注意)。
このような観測値の離散化をどう説明するか。一つの素直な考え方は、本来連続的な値をとるものが観測過程での何かのメカニズムで離散的になると考え、そのメカニズムを追求することだろう。根底では連続なものが何かの原因で離散化すると考えるのである。しかし量子力学ではそのような方針はとらない。観測値が離散化することこそが根底の法則、自然の本性であると捉え、それを適切に記述する数学的な法則を与えるのである。以下でその法則を述べるが、これは別の法則(例えば古典力学)から論理的に導かれるものではない。これ自身がいわば「公理」であり、それが正しいかどうかは実験結果を予言・再現できるかにより判定される。特に古典力学や電磁気学はある範囲内で実験と合うことが確立された理論であるから、それらもこの「公理」から導出されなければならない。
以下、まずスピンという特殊な場合について述べる。一般的な場合への拡張は、少なくとも形式上は単純なこととなる。
スピンの量子力学公理1 物理量の数学的表現は行列、測定値が取りうる値はその固有値
[編集]
スピンがもつ物理量はそのx成分sx、y成分sy, z成分szである(一般の方向の成分はこれらの一次結合となる)。これらの物理量はある2x2行列(Pauli行列と呼ばれる)に対応する。sxに対応する行列をXと書き、sy,szそれぞれに対応する行列をY,Zと書くことにすると、
「対応する」の意味は段階的に説明していく。まず重要なのは観測値が取り得る値が決められることである。物理量szを例にとると
szの測定値になりえるのは対応する行列Zの固有値のみ。つまりszの測定値が1, -1に限られるのは、対応するZの固有値が1, -1だけだからである。
sx, syについても同様。どちらも対応する行列X,Yの固有値が1, -1しかないので測定値としては1, -1しか現れない。
一般的に言うと、スピンに限らず物理量はある行列(より一般的には線形演算子)に対応する。つまりその物理量の観測値は対応する行列の固有値に限られる。逆にいうとある物理量がとりうる測定値を知りたいと思ったら、それに対応する行列を求めその固有値を求めればよい、ということである。
では対応する行列(線形演算子)をどう求めるかという疑問が当然浮かぶが、それは扱う物理系個々の問題になる。基本的なのは正準量子化と呼ばれる手法で、古典力学のポアソン括弧と呼ばれるものから演算子が満たすべき交換関係を推測し、それを満たすような演算子を探す。特にそのような手法から角運動量の一般論を展開でき、上で導入したスピンのPauli行列はその特別な場合として得られる。
但し基本的な物理量の行列が分かれば、それらの関数になっている物理量の行列は行列代数により得られる。例えば
という物理量に対応する行列は
になる。その固有値は(普通の線形代数の方法で計算すると)
なので
がとりえる観測値は
のどちらかになる。従って扱う系で基本的な物理量(例えば位置xと運動量p)に対応する演算子が分かれば、その系の任意の物理量(基本的な物理量の関数になっている量、例えばエネルギー p^2/(2m)+V(x))が対応する行列は演算子の代数(線形代数)で分かってしまうのである。
スピンの量子力学公理2 物理状態の数学的表現は複素列ベクトル=ケット、確率は固有ベクトルとの内積の絶対値自乗
[編集]
次に、離散的な観測値(物理量に対応する行列の固有値)のどれが観測されるかの規則を述べる。まず、スピンの状態はスピンの物理量に対応する行列(X,Y,Z)が作用するベクトル、即ち2行1列の複素列ベクトル
で表される。ここで
はどぢらも複素数だが、次の正規化条件を満たすものとする
。このような複素列ベクトルをDiracによる「ブラケット記号」を使い、「ケット」
で表す。
線形代数によると、複素列ベクトル同士の間には自然に内積が定義される。ケット
と
の間の内積
は次で与えられる:
この内積は、成分が実数の場合には普通の実ベクトル同士の内積になるが、複素数の場合には左側の要素に複素共役を取る。こう定義する理由は、「自分自身との内積
」が必ず0以上の実数になるようにするためである。特に上記の規格化条件は
と書くことができる。
では、このケットの物理的意味について述べよう。上で導入した「z向きの状態」、即ちszを観測すれば必ずsz=1の結果を得られる状態は、行列Zの固有値1の固有ベクトル(で規格化されたもの)で記述される。つまり「szを測定すれば、結果が必ずsz=1になる状態を表すケット」を|sz=1>と書くことにすると次が成立:
(実際にはこの成分で1の代わりに絶対値1の任意の複素数
をおいても同じなのだが、その任意性については後で述べる。とりあえず単位固有ベクトルの中で一番成分が簡単なものを選んだと考えてほしい。)
これは一般的に拡張される。即ち、任意の物理量Aに対して、Aを測定したときに確実に測定値aが得られる状態は、数学的にはAに対応する行列(または演算子)の固有値aの固有ベクトル(で規格化されたもの)で記述される。このような状態を、やや略した言い方で物理量Aの固有値aの固有状態と呼ぶ。
スピンの例に戻ると、sxを測定して必ずsx=1の結果を得られる状態|sx=1>は、行列Xの固有値1の固有状態なので
必ずsx=-1の結果を得られる状態は
では、szを測定する際に、状態がszの固有状態でなかったら結果はどうなるだろうか。ここで始めて実験と比べられる記述が現れる。量子力学が与える予言は次の通り。
スピンの状態
がszの固有状態とは限らない場合にszを測定すると、その結果(測定値が1になるか-1になるか)はランダムな事象となり、確率的にしか予言できない。その確率はもとの状態と固有状態の内積の絶対値自乗で与えられる。例えばsz=1となる確率P(sz=1)は
同様に、sx=1になる確率P(sx=1)は
特に前の例として取り上げた、スピンがx方向を向いている場合にszを測定した結果は、
となり、実験結果とあう確率が与えられる。
さらに、実験例として取り上げた「z方向を向いた電子スピンの
方向の成分
を測る」場合の結果を計算してみる。法則から
問題は
だが、
の成分を球座標での方向
を使い
と表すと
その固有値1の固有ベクトル(で規格化されているもの)は
これを使うと次が得られ、実験結果をきちんと再現する計算結果となる。
ここまで状態の数学的な記述と測定の関係を書いた。次は状態の運動法則である。簡単な例としてスピンに一様な磁場をかけた場合を考える。
このページ「
量子力学」は、
まだ書きかけです。加筆・訂正など、協力いただける皆様の
編集を心からお待ちしております。また、ご意見などがありましたら、お気軽に
トークページへどうぞ。