量子力学の様々な構成方法 第2部 第1章 §2.1.3 ★★★ 古典の運動方程式を表の方程式として要請する。エネルギーの表は対角で、時間微分はエネルギーの表との交換子になる
運動方程式を表の上で要請する
Heisenberg が古典から持ち込んだ3つ目の要請は、運動方程式の形である。1次元でポテンシャル V の中を動く質量 m の粒子について、
x˙=mp,p˙=−V′(x)(2.1.13)
を表の等式として要請する。V′(x) は、V′ が多項式なら積の規則 (2.1.12) で計算し、一般の関数はその極限で定める。古典では x(t) の時間発展を解くことがこの方程式の意味だったが、表の時間依存は式 (2.1.8) で既に決めてある。したがって式 (2.1.13) は、振幅 xnm と振動数 ωnm が満たすべき代数方程式になる。
エネルギーの表は対角になる
エネルギーの表を H=p2/(2m)+V(x) と置く。定常状態はエネルギーが一定の状態だから、H は時間によらないはずである。ところが H も遷移の表なので、その要素は Hnmeiωnmt と振動する。時間によらないためには、ωnm=0 の要素がすべて0でなければならない。準位に縮退がなければ、H は対角になる。対角要素 Hnn を状態 n のエネルギー En とみなすと、振動数条件 (2.1.3) は
Hnm=Enδnm,ℏωnm=Hnn−Hmm(2.1.14)
と書ける。ここで「H が時間によらない」は要請として置いた。運動方程式 (2.1.13) から H˙=0 を導くには、次節の交換関係が要る(Born と Jordan は交換関係を使ってエネルギーの保存を証明した)。
式 (2.1.14) を使うと、時間微分を表の積で書ける。対角な H との交換子の要素は(xH−Hx)nm=xnm(Em−En)=−ℏωnmxnm なので、
x˙=iℏ1(xH−Hx)(2.1.15)
が成り立つ。どの表についても同じである。これが Heisenberg の運動方程式で、古典の x˙={x,H} と同じ形をしている。一般のハミルトニアンでのこの式と、古典の正準方程式との対応は第2部第3章で扱う。
対角要素の時間微分は0
式 (2.1.8) から、どんな表 F についても
(dtdF)nn=iωnnFnn=0(2.1.16)
である。古典では、周期運動で有界な量 F の時間微分は、1周期で平均すると0になる(T1∫0TF˙dt=TF(T)−F(0)=0)。対角要素は古典の時間平均の量子版だったので(例題 2.1.2)、式 (2.1.16) は「時間微分の時間平均は0」の量子版である。この性質から、古典の時間平均についての関係式が表の対角要素の関係式として移ってくる。その代表がビリアル定理で、次の例題で導く。
例題 2.1.3 ビリアル定理を表の対角要素で出す ★★
オリジナルF=21(xp+px) と置く。
- (1) 運動方程式 (2.1.13) を使って F˙ を求め、式 (2.1.16) から 2Tnn=(xV′(x))nn(T=p2/(2m))を導け。
- (2) 調和振動子 V=21mω2x2 では Tnn と Vnn の関係はどうなるか。
- (3) ℏ=m=ω=1 となる単位で V=21x2+x4 とする。基底状態のエネルギーは E0=0.80377、(x2)00=0.25715 である(数値対角化の値)。T00 と (x4)00 を求めよ。
解答
(1) Leibniz 則と式 (2.1.13) から
F˙=21(x˙p+xp˙+p˙x+px˙)=mp2−21(xV′(x)+V′(x)x).
V′(x) は x の多項式なので x と可換で、xV′(x)=V′(x)x。対角要素をとると、式 (2.1.16) から0=2Tnn−(xV′(x))nn。古典の時間平均のビリアル定理 2T=xV′ と同じ形である。
(2) xV′=mω2x2=2V なので Tnn=Vnn。エネルギー En=Tnn+Vnn は運動と位置の部分に半分ずつ分かれる。
(3) xV′=x2+4x4 だから 2T00=(x2)00+4(x4)00。一方 E0=T00+21(x2)00+(x4)00。第2式を4倍して第1式を引くと 4E0−2T00=4T00+(x2)00、すなわち
T00=64E0−(x2)00=64×0.80377−0.25715=0.49298,
(x4)00=E0−T00−21(x2)00=0.80377−0.49298−0.12858=0.18221。どちらも数値対角化で直接求めた値と5桁で一致する。4次の項があると T00>V00=0.31079 となり、(2) の等分は成り立たない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com