tecktom.com

出てくるもの1:2次の Hamiltonian では古典解がそのまま演算子の解

2026年10月15日11分で読めますAI生成教科書『虎の巻』第4部 第3章 Heisenberg 描像は Hamilton 方程式の量子版
目次
  1. 線形の方程式は古典と同じ解を持つ
  2. 分散も古典の集団と同じ式で動く
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第4部 第3章 §4.3.3 ★★★ H^\hat H が q^\hat q, p^\hat p の2次式までなら Heisenberg 方程式は線形で、古典の解の式に初期値の演算子を入れたものが解になる。期待値と分散は古典の集団と同じ式で動く

線形の方程式は古典と同じ解を持つ

H^\hat H が q^\hat q, p^\hat p の2次式(1次式と定数を含む)のとき、交換子 [q^,H^][\hat q,\hat H], [p^,H^][\hat p,\hat H] は q^\hat q, p^\hat p の1次式と定数になる。したがって Heisenberg 方程式は q^H\hat q_{\mathrm H}, p^H\hat p_{\mathrm H} についての線形の常微分方程式で、係数は数である。古典の Hamilton 方程式も同じ係数を持つ同じ線形方程式なので、古典の解を

(q(t)p(t))=M(t)(qp)+(cq(t)cp(t))\begin{pmatrix}q(t)\\p(t)\end{pmatrix}=M(t)\begin{pmatrix}q\\p\end{pmatrix}+\begin{pmatrix}c_q(t)\\c_p(t)\end{pmatrix}

と書いたとき(M(t)M(t) は 2×22\times2 の行列、cqc_q, cpc_p は外力などに由来する数の項)、量子の解は初期値を演算子に替えた

(q^H(t)p^H(t))=M(t)(q^p^)+(cq(t)cp(t))(4.3.10)\begin{pmatrix}\hat q_{\mathrm H}(t)\\\hat p_{\mathrm H}(t)\end{pmatrix}=M(t)\begin{pmatrix}\hat q\\\hat p\end{pmatrix}+\begin{pmatrix}c_q(t)\\c_p(t)\end{pmatrix}\tag{4.3.10}

である。右辺を左辺の方程式に代入すると、係数が数なので演算子の積の順番の問題が起きず、古典の場合と同じ計算でそのまま成り立つ。第4部第1章で問題にした「積の順序」は、方程式の右辺が2次以上の項を含むときにだけ現れる。

調和振動子 H^=p^22m+12mω2q^2\hat H=\frac{\hat p^2}{2m}+\frac12m\omega^2\hat q^2 では、Heisenberg 方程式 q^˙H=p^H/m\dot{\hat q}_{\mathrm H}=\hat p_{\mathrm H}/m, p^˙H=−mω2q^H\dot{\hat p}_{\mathrm H}=-m\omega^2\hat q_{\mathrm H} の解が

q^H(t)=q^cos⁡ωt+p^mωsin⁡ωt,p^H(t)=−mωq^sin⁡ωt+p^cos⁡ωt(4.3.11)\hat q_{\mathrm H}(t)=\hat q\cos\omega t+\frac{\hat p}{m\omega}\sin\omega t,\qquad \hat p_{\mathrm H}(t)=-m\omega\hat q\sin\omega t+\hat p\cos\omega t\tag{4.3.11}

になる。M(t)M(t) は相空間の(目盛りを mωm\omega で揃えた)回転で、行列式は cos⁡2ωt+sin⁡2ωt=1\cos^2\omega t+\sin^2\omega t=1 である。式 (4.3.8) の交換子が保たれることは、この行列式が 11 であることに当たる。消滅演算子 a^=mω2ℏ(q^+ip^mω)\hat a=\sqrt{\frac{m\omega}{2\hbar}}\bigl(\hat q+\frac{\mathrm{i}\hat p}{m\omega}\bigr) で書くと、式 (4.3.11) は a^H(t)=a^ e−iωt\hat a_{\mathrm H}(t)=\hat a\,\mathrm{e}^{-\mathrm{i}\omega t} という1行にまとまる。古典の複素振幅 mω/2 (q+ip/(mω))\sqrt{m\omega/2}\,\bigl(q+\mathrm{i} p/(m\omega)\bigr) も同じく e−iωt\mathrm{e}^{-\mathrm{i}\omega t} で回り、その絶対値の2乗が作用変数 J=E/ωJ=E/\omega である(第1部第14章)。この対応から出るエネルギー準位は 第4部第4章で扱う。

分散も古典の集団と同じ式で動く

式 (4.3.10) の両辺の期待値をとると、平均値 ⟨q^H⟩\langle \hat q_{\mathrm H} \rangle, ⟨p^H⟩\langle \hat p_{\mathrm H} \rangle は古典の軌道と同じ式で動く。2次の量も同じで、自由粒子 q^H=q^+p^t/m\hat q_{\mathrm H}=\hat q+\hat pt/m の分散は、積の順番に注意して展開すると

(Δq(t))2=(Δq)2+tmCqp+t2m2(Δp)2,Cqp≔⟨q^p^+p^q^⟩−2⟨q^⟩⟨p^⟩(4.3.12)(\Delta q(t))^2=(\Delta q)^2+\frac{t}{m}C_{qp}+\frac{t^2}{m^2}(\Delta p)^2,\qquad C_{qp}\coloneqq\langle \hat q\hat p+\hat p\hat q \rangle -2\langle \hat q \rangle \langle \hat p \rangle\tag{4.3.12}

である。Δ\Delta は時刻 00 の標準偏差、CqpC_{qp} は位置と運動量の相関を対称に並べたもの(古典の共分散の2倍に当たる)である。古典の点の集団を、同じ平均・分散・共分散で用意して自由運動させても、q(t)=q+pt/mq(t)=q+pt/m の分散はまったく同じ式 (4.3.12) で動く。線形の写像で運ばれる分布の分散は、写像の係数と初期の2次の量だけで決まるからである。

量子力学が付け加えるのは、初期の分散に課される制限、すなわち不確定性関係 Δq Δp≥ℏ/2\Delta q\,\Delta p\ge\hbar/2 だけである。等号の成り立つ Gauss 波束(Δq=σ\Delta q=\sigma, Δp=ℏ/(2σ)\Delta p=\hbar/(2\sigma), Cqp=0C_{qp}=0)を式 (4.3.12) に入れると、第4部第2章の式 (4.2.19) の波束の広がり

Δq(t)=σ1+(tτ)2,τ≔2mσ2ℏ(4.3.13)\Delta q(t)=\sigma\sqrt{1+\Bigl(\frac{t}{\tau}\Bigr)^2},\qquad \tau\coloneqq\frac{2m\sigma^2}{\hbar}\tag{4.3.13}

が出る。第2章ではこの式を期待値のモーメントの方程式から出した。ここでは Heisenberg の解 q^H=q^+p^t/m\hat q_{\mathrm H}=\hat q+\hat pt/m の分散を、古典の集団の分散と同じ1行の計算で出したことになる。時間 τ\tau より後は Δq(t)≈σt/τ=Δp t/m\Delta q(t)\approx\sigma t/\tau=\Delta p\,t/m で、運動量のばらつき Δp\Delta p を持つ粒子の集団が速度の差で散らばっていく広がりそのものである。位置を細かく絞った波束(小さい σ\sigma)ほど Δp\Delta p が大きく、早く広がる。

図

図 4.3.2 最小不確定の Gauss 波束の幅 Δq(t)\Delta q(t)(式 (4.3.13)、点線=厳密)。t≪τt\ll\tau では初めの幅 σ\sigma のまま(灰の実線)、t≫τt\gg\tau では運動量のばらつき Δp=ℏ/(2σ)\Delta p=\hbar/(2\sigma) による広がり σt/τ=Δp t/m\sigma t/\tau=\Delta p\,t/m(赤の実線)に近づく。同じ平均・分散を持つ古典の点の集団も、まったく同じ曲線に従う

例題 4.3.3 逆さの振動子で波束は指数関数的に広がる ★★

オリジナル逆さの振動子 H^=p^22m−12mΩ2q^2\hat H=\frac{\hat p^2}{2m}-\frac12m\Omega^2\hat q^2(Ω>0\Omega>0)を考える。原点は不安定な釣り合いの点である。

  1. (1) Heisenberg 方程式を解いて q^H(t)\hat q_{\mathrm H}(t), p^H(t)\hat p_{\mathrm H}(t) を求め、[q^H(t),p^H(t)]=iℏ[\hat q_{\mathrm H}(t),\hat p_{\mathrm H}(t)]=\mathrm{i}\hbar が保たれることを確かめよ。
  2. (2) 時刻 00 に原点を中心とする最小不確定の Gauss 波束(Δq=σ\Delta q=\sigma, Δp=ℏ/(2σ)\Delta p=\hbar/(2\sigma), Cqp=0C_{qp}=0)があるとき、Δq(t)\Delta q(t) を求めよ。特に σ2=ℏ/(2mΩ)\sigma^2=\hbar/(2m\Omega) のとき Δq(t)=σcosh⁡2Ωt\Delta q(t)=\sigma\sqrt{\cosh2\Omega t} となり、Ωt≫1\Omega t\gg1 で eΩt\mathrm{e}^{\Omega t} に比例して広がることを示せ。
  3. (3) (2) の波束で Ωt=2\Omega t=2 のときの Δq/σ\Delta q/\sigma と、幅が初めの 1010 倍になる Ωt\Omega t を求めよ。m=Ω=1m=\Omega=1 の単位で、第3部第1章の例題 3.1.2(相空間の正方形の引き伸ばし)と比べて、この広がりを相空間の言葉で説明せよ。
解答

(1) Heisenberg 方程式は q^˙H=p^H/m\dot{\hat q}_{\mathrm H}=\hat p_{\mathrm H}/m, p^˙H=mΩ2q^H\dot{\hat p}_{\mathrm H}=m\Omega^2\hat q_{\mathrm H} で、右辺は1次式なので、4.3.3節のとおり古典の解の初期値を演算子に替えればよい。

q^H(t)=q^cosh⁡Ωt+p^mΩsinh⁡Ωt,p^H(t)=mΩq^sinh⁡Ωt+p^cosh⁡Ωt.\hat q_{\mathrm H}(t)=\hat q\cosh\Omega t+\frac{\hat p}{m\Omega}\sinh\Omega t,\qquad \hat p_{\mathrm H}(t)=m\Omega\hat q\sinh\Omega t+\hat p\cosh\Omega t .

[q^,q^]=[p^,p^]=0[\hat q,\hat q]=[\hat p,\hat p]=0 を使うと [q^H,p^H]=cosh⁡2Ωt [q^,p^]+sinh⁡2Ωt [p^,q^]=(cosh⁡2Ωt−sinh⁡2Ωt) iℏ=iℏ[\hat q_{\mathrm H},\hat p_{\mathrm H}]=\cosh^2\Omega t\,[\hat q,\hat p]+\sinh^2\Omega t\,[\hat p,\hat q]=(\cosh^2\Omega t-\sinh^2\Omega t)\,\mathrm{i}\hbar=\mathrm{i}\hbar。係数の行列の行列式が cosh⁡2−sinh⁡2=1\cosh^2-\sinh^2=1 であることが、交換関係の保存に当たる。

(2) 平均は 00 のままなので、Δq(t)2=⟨q^H2⟩\Delta q(t)^2=\langle \hat q_{\mathrm H}^2 \rangle である。⟨q^p^+p^q^⟩=Cqp=0\langle \hat q\hat p+\hat p\hat q \rangle =C_{qp}=0 から交差項は消え、

Δq(t)2=σ2cosh⁡2Ωt+(ℏ2mΩσ)2sinh⁡2Ωt.\Delta q(t)^2=\sigma^2\cosh^2\Omega t+\Bigl(\frac{\hbar}{2m\Omega\sigma}\Bigr)^2\sinh^2\Omega t .

σ2=ℏ/(2mΩ)\sigma^2=\hbar/(2m\Omega) なら第2項の係数も σ2\sigma^2 になり、Δq(t)2=σ2(cosh⁡2Ωt+sinh⁡2Ωt)=σ2cosh⁡2Ωt\Delta q(t)^2=\sigma^2(\cosh^2\Omega t+\sinh^2\Omega t)=\sigma^2\cosh2\Omega t。Ωt≫1\Omega t\gg1 では cosh⁡2Ωt≈e2Ωt/2\cosh2\Omega t\approx\mathrm{e}^{2\Omega t}/2 なので Δq(t)≈σeΩt/2\Delta q(t)\approx\sigma\mathrm{e}^{\Omega t}/\sqrt2 となり、自由粒子の tt に比例する広がり(式 (4.3.13))より速く広がる。

(3) Ωt=2\Omega t=2 では Δq/σ=cosh⁡4=27.31=5.23\Delta q/\sigma=\sqrt{\cosh4}=\sqrt{27.31}=5.23。幅が 1010 倍になるのは cosh⁡2Ωt=100\cosh2\Omega t=100 のときで、Ωt=12arcosh⁡100=2.65\Omega t=\frac12\operatorname{arcosh}100=2.65。m=Ω=1m=\Omega=1 では H^=12p^2−12q^2\hat H=\frac12\hat p^2-\frac12\hat q^2 で、(1) の係数の行列は第3部第1章の例題の M(t)M(t) と同じである。変数 (q,p/(mΩ))(q,p/(m\Omega)) で見ると、σ2=ℏ/(2mΩ)\sigma^2=\hbar/(2m\Omega) の波束の分布は半径 σ\sigma の円の形をしている。M(t)M(t) はこれを (1,1)(1,1) 方向に eΩt\mathrm{e}^{\Omega t} 倍、(1,−1)(1,-1) 方向に e−Ωt\mathrm{e}^{-\Omega t} 倍した楕円に運ぶ。t=2t=2 では第3部第1章の正方形の長い対角線と同じく e2=7.39\mathrm{e}^2=7.39 倍である。この楕円を qq 軸に射影した幅が σ(e2Ωt+e−2Ωt)/2=σcosh⁡2Ωt\sigma\sqrt{(\mathrm{e}^{2\Omega t}+\mathrm{e}^{-2\Omega t})/2}=\sigma\sqrt{\cosh2\Omega t} で、(2) の結果と一致する。楕円の面積は det⁡M=1\det M=1 によって変わらず、量子力学が付け加えるのは、その面積が不確定性関係で ℏ/2\hbar/2 より小さくできないという初期条件だけである。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ