解析力学から様々な理論分野を組み立てる 第4部 第3章 §4.3.3 ★★★ H^ が q^, p^ の2次式までなら Heisenberg 方程式は線形で、古典の解の式に初期値の演算子を入れたものが解になる。期待値と分散は古典の集団と同じ式で動く
線形の方程式は古典と同じ解を持つ
H^ が q^, p^ の2次式(1次式と定数を含む)のとき、交換子 [q^,H^], [p^,H^] は q^, p^ の1次式と定数になる。したがって Heisenberg 方程式は q^H, p^H についての線形の常微分方程式で、係数は数である。古典の Hamilton 方程式も同じ係数を持つ同じ線形方程式なので、古典の解を
(q(t)p(t))=M(t)(qp)+(cq(t)cp(t))
と書いたとき(M(t) は 2×2 の行列、cq, cp は外力などに由来する数の項)、量子の解は初期値を演算子に替えた
(q^H(t)p^H(t))=M(t)(q^p^)+(cq(t)cp(t))(4.3.10)
である。右辺を左辺の方程式に代入すると、係数が数なので演算子の積の順番の問題が起きず、古典の場合と同じ計算でそのまま成り立つ。第4部第1章で問題にした「積の順序」は、方程式の右辺が2次以上の項を含むときにだけ現れる。
調和振動子 H^=2mp^2+21mω2q^2 では、Heisenberg 方程式 q^˙H=p^H/m, p^˙H=−mω2q^H の解が
q^H(t)=q^cosωt+mωp^sinωt,p^H(t)=−mωq^sinωt+p^cosωt(4.3.11)
になる。M(t) は相空間の(目盛りを mω で揃えた)回転で、行列式は cos2ωt+sin2ωt=1 である。式 (4.3.8) の交換子が保たれることは、この行列式が 1 であることに当たる。消滅演算子 a^=2ℏmω(q^+mωip^) で書くと、式 (4.3.11) は a^H(t)=a^e−iωt という1行にまとまる。古典の複素振幅 mω/2(q+ip/(mω)) も同じく e−iωt で回り、その絶対値の2乗が作用変数 J=E/ω である(第1部第14章)。この対応から出るエネルギー準位は 第4部第4章で扱う。
分散も古典の集団と同じ式で動く
式 (4.3.10) の両辺の期待値をとると、平均値 ⟨q^H⟩, ⟨p^H⟩ は古典の軌道と同じ式で動く。2次の量も同じで、自由粒子 q^H=q^+p^t/m の分散は、積の順番に注意して展開すると
(Δq(t))2=(Δq)2+mtCqp+m2t2(Δp)2,Cqp:=⟨q^p^+p^q^⟩−2⟨q^⟩⟨p^⟩(4.3.12)
である。Δ は時刻 0 の標準偏差、Cqp は位置と運動量の相関を対称に並べたもの(古典の共分散の2倍に当たる)である。古典の点の集団を、同じ平均・分散・共分散で用意して自由運動させても、q(t)=q+pt/m の分散はまったく同じ式 (4.3.12) で動く。線形の写像で運ばれる分布の分散は、写像の係数と初期の2次の量だけで決まるからである。
量子力学が付け加えるのは、初期の分散に課される制限、すなわち不確定性関係 ΔqΔp≥ℏ/2 だけである。等号の成り立つ Gauss 波束(Δq=σ, Δp=ℏ/(2σ), Cqp=0)を式 (4.3.12) に入れると、第4部第2章の式 (4.2.19) の波束の広がり
Δq(t)=σ1+(τt)2,τ:=ℏ2mσ2(4.3.13)
が出る。第2章ではこの式を期待値のモーメントの方程式から出した。ここでは Heisenberg の解 q^H=q^+p^t/m の分散を、古典の集団の分散と同じ1行の計算で出したことになる。時間 τ より後は Δq(t)≈σt/τ=Δpt/m で、運動量のばらつき Δp を持つ粒子の集団が速度の差で散らばっていく広がりそのものである。位置を細かく絞った波束(小さい σ)ほど Δp が大きく、早く広がる。

図 4.3.2 最小不確定の Gauss 波束の幅 Δq(t)(式 (4.3.13)、点線=厳密)。t≪τ では初めの幅 σ のまま(灰の実線)、t≫τ では運動量のばらつき Δp=ℏ/(2σ) による広がり σt/τ=Δpt/m(赤の実線)に近づく。同じ平均・分散を持つ古典の点の集団も、まったく同じ曲線に従う
例題 4.3.3 逆さの振動子で波束は指数関数的に広がる ★★
オリジナル逆さの振動子 H^=2mp^2−21mΩ2q^2(Ω>0)を考える。原点は不安定な釣り合いの点である。
- (1) Heisenberg 方程式を解いて q^H(t), p^H(t) を求め、[q^H(t),p^H(t)]=iℏ が保たれることを確かめよ。
- (2) 時刻 0 に原点を中心とする最小不確定の Gauss 波束(Δq=σ, Δp=ℏ/(2σ), Cqp=0)があるとき、Δq(t) を求めよ。特に σ2=ℏ/(2mΩ) のとき Δq(t)=σcosh2Ωt となり、Ωt≫1 で eΩt に比例して広がることを示せ。
- (3) (2) の波束で Ωt=2 のときの Δq/σ と、幅が初めの 10 倍になる Ωt を求めよ。m=Ω=1 の単位で、第3部第1章の例題 3.1.2(相空間の正方形の引き伸ばし)と比べて、この広がりを相空間の言葉で説明せよ。
解答
(1) Heisenberg 方程式は q^˙H=p^H/m, p^˙H=mΩ2q^H で、右辺は1次式なので、4.3.3節のとおり古典の解の初期値を演算子に替えればよい。
q^H(t)=q^coshΩt+mΩp^sinhΩt,p^H(t)=mΩq^sinhΩt+p^coshΩt.
[q^,q^]=[p^,p^]=0 を使うと [q^H,p^H]=cosh2Ωt[q^,p^]+sinh2Ωt[p^,q^]=(cosh2Ωt−sinh2Ωt)iℏ=iℏ。係数の行列の行列式が cosh2−sinh2=1 であることが、交換関係の保存に当たる。
(2) 平均は 0 のままなので、Δq(t)2=⟨q^H2⟩ である。⟨q^p^+p^q^⟩=Cqp=0 から交差項は消え、
Δq(t)2=σ2cosh2Ωt+(2mΩσℏ)2sinh2Ωt.
σ2=ℏ/(2mΩ) なら第2項の係数も σ2 になり、Δq(t)2=σ2(cosh2Ωt+sinh2Ωt)=σ2cosh2Ωt。Ωt≫1 では cosh2Ωt≈e2Ωt/2 なので Δq(t)≈σeΩt/2 となり、自由粒子の t に比例する広がり(式 (4.3.13))より速く広がる。
(3) Ωt=2 では Δq/σ=cosh4=27.31=5.23。幅が 10 倍になるのは cosh2Ωt=100 のときで、Ωt=21arcosh100=2.65。m=Ω=1 では H^=21p^2−21q^2 で、(1) の係数の行列は第3部第1章の例題の M(t) と同じである。変数 (q,p/(mΩ)) で見ると、σ2=ℏ/(2mΩ) の波束の分布は半径 σ の円の形をしている。M(t) はこれを (1,1) 方向に eΩt 倍、(1,−1) 方向に e−Ωt 倍した楕円に運ぶ。t=2 では第3部第1章の正方形の長い対角線と同じく e2=7.39 倍である。この楕円を q 軸に射影した幅が σ(e2Ωt+e−2Ωt)/2=σcosh2Ωt で、(2) の結果と一致する。楕円の面積は detM=1 によって変わらず、量子力学が付け加えるのは、その面積が不確定性関係で ℏ/2 より小さくできないという初期条件だけである。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com