tecktom.com

当てる道具:相互作用描像と1次の摂動

2026年10月15日14分で読めますAI生成教科書『虎の巻』第4部 第12章 光と物質
目次
  1. 相互作用描像
  2. 1次の摂動
  3. 正しさの範囲:2準位の厳密解と比べる
  4. 古典の対応物:定数変化法
  5. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第4部 第12章 §4.12.2 ★★★ 解ける H^0\hat H_0 の時間発展に乗った描像で、残りの H^1\hat H_1 だけで状態を動かす。1次の遷移確率は古典の振動子が外力から吸うエネルギーと同じ関数になる

相互作用描像

第4部第3章(第4部第3章)で、Hamiltonian を H^0+H^1(t)\hat H_0+\hat H_1(t) と分けたときの相互作用描像を導入した。状態と摂動を

∣ψI(t)⟩=eiH^0t/ℏ∣ψ(t)⟩,H^1,I(t)=eiH^0t/ℏH^1(t)e−iH^0t/ℏ(4.12.8)| \psi_{\mathrm I}(t) \rangle =\mathrm{e}^{\mathrm{i}\hat H_0t/\hbar}| \psi(t) \rangle ,\qquad \hat H_{1,\mathrm I}(t)=\mathrm{e}^{\mathrm{i}\hat H_0t/\hbar}\hat H_1(t)\mathrm{e}^{-\mathrm{i}\hat H_0t/\hbar}\tag{4.12.8}

と定めると、

iℏddt∣ψI(t)⟩=H^1,I(t)∣ψI(t)⟩(4.12.9)\mathrm{i}\hbar\frac{\mathrm{d} }{\mathrm{d} t}| \psi_{\mathrm I}(t) \rangle =\hat H_{1,\mathrm I}(t)| \psi_{\mathrm I}(t) \rangle\tag{4.12.9}

となる。H^1=0\hat H_1=0 なら ∣ψI⟩| \psi_{\mathrm I} \rangle は動かない。H^0\hat H_0 による時間発展を状態から取り除き、H^1\hat H_1 による変化だけを残したのが相互作用描像である。

状態を H^0\hat H_0 の固有状態で ∣ψI(t)⟩=∑ncn(t)∣n⟩| \psi_{\mathrm I}(t) \rangle =\sum_nc_n(t)| n \rangle と展開すると、cnc_n は「H^0\hat H_0 だけで動いたときの位相 e−iEnt/ℏ\mathrm{e}^{-\mathrm{i} E_nt/\hbar} を除いた振幅」で、式 (4.12.9) は

iℏ c˙f(t)=∑n⟨f∣H^1(t)∣n⟩ eiωfntcn(t),ωfn=Ef−Enℏ(4.12.10)\mathrm{i}\hbar\,\dot c_f(t)=\sum_n\langle f | \hat H_1(t) | n \rangle \,\mathrm{e}^{\mathrm{i}\omega_{fn}t}c_n(t),\qquad \omega_{fn}=\frac{E_f-E_n}{\hbar}\tag{4.12.10}

になる。eiH^0t/ℏ\mathrm{e}^{\mathrm{i}\hat H_0t/\hbar} を ⟨f∣\langle f | と ∣n⟩| n \rangle に当てると位相 ei(Ef−En)t/ℏ\mathrm{e}^{\mathrm{i}(E_f-E_n)t/\hbar} が出ることを使った。

1次の摂動

時刻 00 に状態 ∣i⟩| i \rangle にいたとする(cn(0)=δnic_n(0)=\delta_{ni})。H^1\hat H_1 が小さいなら、右辺の cn(t)c_n(t) を初めの値で置き換えてよく、

cf(t)≈1iℏ∫0t⟨f∣H^1(t′)∣i⟩ eiωfit′dt′(f≠i)(4.12.11)c_f(t)\approx\frac{1}{\mathrm{i}\hbar}\int_0^t\langle f | \hat H_1(t') | i \rangle \,\mathrm{e}^{\mathrm{i}\omega_{fi}t'}\mathrm{d} t'\qquad(f\neq i)\tag{4.12.11}

を得る。式 (4.12.9) を積分方程式 ∣ψI(t)⟩=∣i⟩+1iℏ∫0tH^1,I(t′)∣ψI(t′)⟩dt′| \psi_{\mathrm I}(t) \rangle =| i \rangle +\frac{1}{\mathrm{i}\hbar}\int_0^t\hat H_{1,\mathrm I}(t')| \psi_{\mathrm I}(t') \rangle \mathrm{d} t' に書き直して、右辺に繰り返し代入した級数(Dyson 級数)の1次の項にあたる。2次以降の項には時刻を順に並べた積 H^1,I(t1)H^1,I(t2)⋯\hat H_{1,\mathrm I}(t_1)\hat H_{1,\mathrm I}(t_2)\cdots(t1>t2>⋯t_1>t_2>\cdots)が現れる。第4部第2章(第4部第2章)で述べた、時刻を含む H^(t)\hat H(t) の時間発展を時間順序積で書く式を、H^1\hat H_1 について展開したものである。

摂動が角振動数 ω\omega で振動するときは、エルミートな H^1(t)\hat H_1(t) を

H^1(t)=H^ω e−iωt+H^ω† eiωt(4.12.12)\hat H_1(t)=\hat H_\omega\,\mathrm{e}^{-\mathrm{i}\omega t}+\hat H_\omega^\dagger\,\mathrm{e}^{\mathrm{i}\omega t}\tag{4.12.12}

と書く。平面波 (4.12.6) なら H^ω=−q2mA0eik⋅r^e⋅p^\hat H_\omega=-\frac{q}{2m}A_0\mathrm{e}^{\mathrm{i}\bm{k}\cdot\hat{\bm{r}}}\bm{e}\cdot\hat{\bm{p}} である(e⋅k=0\bm{e}\cdot\bm{k}=0 なので eik⋅r^\mathrm{e}^{\mathrm{i}\bm{k}\cdot\hat{\bm{r}}} と e⋅p^\bm{e}\cdot\hat{\bm{p}} は交換する)。式 (4.12.11) に入れると、e−iωt\mathrm{e}^{-\mathrm{i}\omega t} の項は ei(ωfi−ω)t′\mathrm{e}^{\mathrm{i}(\omega_{fi}-\omega)t'} を、eiωt\mathrm{e}^{\mathrm{i}\omega t} の項は ei(ωfi+ω)t′\mathrm{e}^{\mathrm{i}(\omega_{fi}+\omega)t'} を積分することになる。Ef>EiE_f>E_i(吸収)で ω\omega が ωfi\omega_{fi} に近いとき、前者の積分は tt に比例して大きくなり、後者は振動して O(1/ω)\mathcal{O}(1/\omega) にとどまる。前者だけを残し、離調を Δ=ω−ωfi\Delta=\omega-\omega_{fi}(光の角振動数 −- 遷移の角振動数。記法の正本の流儀)とおくと、∣∫0te−iΔt′dt′∣2=sin⁡2(Δt/2)/(Δ/2)2\lvert \int_0^t\mathrm{e}^{-\mathrm{i}\Delta t'}\mathrm{d} t'\rvert ^2=\sin^2(\Delta t/2)/(\Delta/2)^2 から

Pf←i(t)=∣cf(t)∣2=∣⟨f∣H^ω∣i⟩∣2ℏ2 sin⁡2(Δt/2)(Δ/2)2(4.12.13)P_{f\leftarrow i}(t)=\lvert c_f(t)\rvert ^2=\frac{\lvert \langle f | \hat H_\omega | i \rangle \rvert ^2}{\hbar^2}\,\frac{\sin^2(\Delta t/2)}{(\Delta/2)^2}\tag{4.12.13}

を得る。共鳴(Δ=0\Delta=0)では P=∣⟨f∣H^ω∣i⟩∣2t2/ℏ2P=\lvert \langle f | \hat H_\omega | i \rangle \rvert ^2t^2/\hbar^2 で、時間の2乗で増える。Ef<EiE_f<E_i の場合は、逆に eiωt\mathrm{e}^{\mathrm{i}\omega t} の項が共鳴し、H^ω\hat H_\omega を H^ω†\hat H_\omega^\dagger に替えた同じ式になる。これが誘導放出で、4.12.5節で使う。

正しさの範囲:2準位の厳密解と比べる

1次の式 (4.12.13) は確率なので、11 を超えたら明らかに使えない。どこまで使えるかを、厳密に解ける2準位の場合で確かめる。双極子の形(4.12.4節)で H^ω=−d^⋅eE0/2\hat H_\omega=-\hat{\bm{d}}\cdot\bm{e} E_0/2 とし、Rabi 周波数を Ω=deg⋅eE0/ℏ\Omega=\bm{d}_{eg}\cdot\bm{e}E_0/\hbar(位相を選んで実数)とおくと∣⟨e∣H^ω∣g⟩∣=ℏΩ/2\lvert \langle e | \hat H_\omega | g \rangle \rvert =\hbar\Omega/2 で、共鳴の1次の答えは Pe=(Ωt/2)2P_e=(\Omega t/2)^2 である。一方、共鳴しない項を落とした2準位の Schrödinger 方程式はic˙e=Ω2cg\mathrm{i}\dot c_e=\frac\Omega2c_g, ic˙g=Ω2ce\mathrm{i}\dot c_g=\frac\Omega2c_e(共鳴のとき)となり、厳密解は Pe=sin⁡2(Ωt/2)P_e=\sin^2(\Omega t/2) である(Rabi 振動。光学と量子エレクトロニクス 第19章で詳しく扱う。光学と量子エレクトロニクス 第19章)。sin⁡2x≈x2\sin^2x\approx x^2 だから、両者は Ωt≪1\Omega t\ll1 で一致する。図 4.12.1 のとおり、1次の答えが厳密解の 1.11.1 倍になるのは Ωt=1.06\Omega t=1.06、Pe=0.26P_e=0.26 のときで(厳密解が1次の答えの 0.90.9 倍になる所を基準にとると Ωt=1.118\Omega t=1.118、Pe=0.281P_e=0.281 で、光学と量子エレクトロニクス 第18章はこちらの基準で書いている。光学と量子エレクトロニクス 第18章)、遷移確率が 0.10.1 程度までなら1次の式は数 %\% の精度で正しい。

図

図 4.12.1 共鳴の2準位系で、上の準位にいる確率 PeP_e。実線が1次の摂動(近似)、点線が厳密解(共鳴しない項を落とした2準位の Schrödinger 方程式の解)。Ωt≲1\Omega t\lesssim1 で両者は重なり、1次の答えは Ωt=π\Omega t=\pi で 2.472.47 と確率の上限 11 を超える

近似の裏付け

型 B(極限) | 小さいもの x=∥H^1∥t/ℏ→0x=\lVert \hat H_1\rVert t/\hbar\to0(摂動の大きさと時間の積が ℏ\hbar に比べて小さい) | 証明 証明済み(H^1\hat H_1 が有界な演算子のとき、Dyson 級数は絶対収束し、2次以降の和の大きさは ex−1−x≤12x2ex\mathrm{e}^x-1-x\le\frac12x^2\mathrm{e}^x 以下。p^\hat{\bm{p}} を含む H^1\hat H_1 は非有界なので、遷移にかかわる有限個の準位に限ったときにこの評価を使う) | 行き先 近似の数理 第2章

古典の対応物:定数変化法

古典力学では、解ける H0H_0 の解の積分定数(作用角変数や複素振幅)を新しい変数に選び、残りの H1H_1 による積分定数の変化を追う方法がこれに当たる。解ける部分の解に乗った座標に移る操作は、第1部第12章(第1部第12章)で述べた「時間発展そのものが正準変換である」ことを使った正準変換で、相互作用描像の eiH^0t/ℏ\mathrm{e}^{\mathrm{i}\hat H_0t/\hbar} はその量子版である。次の例題で、調和振動子についてこの対応を具体的に確かめる。

例題 4.12.2 外力を受ける振動子:古典の吸収エネルギーと量子の遷移確率 ★★★

オリジナル質量 mm、角振動数 ω0\omega_0 の調和振動子に外力 F(t)F(t) がはたらき、H=p22m+12mω02q2−F(t)qH=\frac{p^2}{2m}+\frac12m\omega_0^2q^2-F(t)q とする。F(t)F(t) は時刻 00 から加わる。

  1. (1) 第4部第4章(第4部第4章)の複素振幅 z=(mω0q+ip)/2mω0z=(m\omega_0q+\mathrm{i} p)/\sqrt{2m\omega_0} について、z˙={z,H}\dot z=\{z,H\} を計算せよ。時刻 00 に静止していた振動子(z(0)=0z(0)=0)が時刻 tt までに得るエネルギー ω0∣z(t)∣2\omega_0\lvert z(t)\rvert ^2 を FF で表せ。
  2. (2) 量子の振動子が時刻 00 に基底状態 ∣0⟩| 0 \rangle にいたとする。H^1=−F(t)q^\hat H_1=-F(t)\hat q の1次の摂動で第1励起状態 ∣1⟩| 1 \rangle に移る確率 P1←0P_{1\leftarrow0} を求め、ℏω0P1←0\hbar\omega_0P_{1\leftarrow0} が (1) の古典のエネルギーと等しいことを示せ。
  3. (3) F(t)=F0cos⁡ωtF(t)=F_0\cos\omega t で、ω\omega が ω0\omega_0 に近いとき、(1) のエネルギーが F028msin⁡2(Δt/2)(Δ/2)2\frac{F_0^2}{8m}\frac{\sin^2(\Delta t/2)}{(\Delta/2)^2}(Δ=ω−ω0\Delta=\omega-\omega_0)で近似できることを示せ。
解答

(1) {z,q}=12mω0{ip,q}=−i2mω0\{z,q\}=\frac{1}{\sqrt{2m\omega_0}}\{\mathrm{i} p,q\}=-\frac{\mathrm{i}}{\sqrt{2m\omega_0}} と、第4部第4章の {z,ω0zz∗}=−iω0z\{z,\omega_0zz^*\}=-\mathrm{i}\omega_0z から

z˙={z,H}=−iω0z−F(t){z,q}=−iω0z+iF(t)2mω0\dot z=\{z,H\}=-\mathrm{i}\omega_0z-F(t)\{z,q\}=-\mathrm{i}\omega_0z+\frac{\mathrm{i} F(t)}{\sqrt{2m\omega_0}}

である。F=0F=0 の解 e−iω0t\mathrm{e}^{-\mathrm{i}\omega_0t} の係数を変化させて z=w(t)e−iω0tz=w(t)\mathrm{e}^{-\mathrm{i}\omega_0t} とおくと、w˙=iF(t)2mω0eiω0t\dot w=\frac{\mathrm{i} F(t)}{\sqrt{2m\omega_0}}\mathrm{e}^{\mathrm{i}\omega_0t} で、ww の変化は外力だけで決まる(これが定数変化法で、ww が相互作用描像の cnc_n にあたる)。w(0)=0w(0)=0 から

ω0∣z(t)∣2=ω0∣w(t)∣2=12m∣∫0tF(t′) eiω0t′dt′∣2.\omega_0\lvert z(t)\rvert ^2=\omega_0\lvert w(t)\rvert ^2=\frac{1}{2m}\Bigl|\int_0^tF(t')\,\mathrm{e}^{\mathrm{i}\omega_0t'}\mathrm{d} t'\Bigr|^2 .

外力が無くなった後の振動子のエネルギーは ω0∣z∣2\omega_0\lvert z\rvert ^2 そのものなので、これが外力から吸い取ったエネルギーである。式は厳密で、近似を使っていない。これは第2部第1章の式 (2.1.8)(ΔE=∣∫Feiω0tdt∣2/2m\Delta E=\lvert \int F\mathrm{e}^{\mathrm{i}\omega_0t}\mathrm{d} t\rvert ^2/2m)と同じものである。量子の ⟨H^0⟩\langle \hat H_0 \rangle は第4部第3章の4.3.7節の総合問題 (4) で出した。

(2) 式 (4.12.11) で ⟨1∣H^1(t)∣0⟩=−F(t)⟨1∣q^∣0⟩=−F(t)ℏ/(2mω0)\langle 1 | \hat H_1(t) | 0 \rangle =-F(t)\langle 1 | \hat q | 0 \rangle =-F(t)\sqrt{\hbar/(2m\omega_0)}(第4部第4章の行列要素)、ω10=ω0\omega_{10}=\omega_0 だから

P1←0=1ℏ2ℏ2mω0∣∫0tF(t′)eiω0t′dt′∣2,ℏω0P1←0=12m∣∫0tF(t′)eiω0t′dt′∣2.P_{1\leftarrow0}=\frac{1}{\hbar^2}\frac{\hbar}{2m\omega_0}\Bigl|\int_0^tF(t')\mathrm{e}^{\mathrm{i}\omega_0t'}\mathrm{d} t'\Bigr|^2,\qquad \hbar\omega_0P_{1\leftarrow0}=\frac{1}{2m}\Bigl|\int_0^tF(t')\mathrm{e}^{\mathrm{i}\omega_0t'}\mathrm{d} t'\Bigr|^2 .

(1) と同じである。ℏ\hbar は消え、量子1個のエネルギー ℏω0\hbar\omega_0 と遷移確率の積が、古典の吸収エネルギーに一致する。q^\hat q は準位を1つしか動かさないので、1次では ∣1⟩| 1 \rangle 以外への遷移は起こらない。

(3) F0cos⁡ωt′ eiω0t′=F02(e−iΔt′+ei(ω0+ω)t′)F_0\cos\omega t'\,\mathrm{e}^{\mathrm{i}\omega_0t'}=\frac{F_0}2\bigl(\mathrm{e}^{-\mathrm{i}\Delta t'}+\mathrm{e}^{\mathrm{i}(\omega_0+\omega)t'}\bigr) の第2項は積分しても O(F0/ω0)\mathcal{O}(F_0/\omega_0) にとどまり、第1項の積分は∣∫0te−iΔt′dt′∣2=sin⁡2(Δt/2)/(Δ/2)2\lvert \int_0^t\mathrm{e}^{-\mathrm{i}\Delta t'}\mathrm{d} t'\rvert ^2=\sin^2(\Delta t/2)/(\Delta/2)^2 で tt とともに大きくなる。第1項だけを残すと、エネルギーは F028msin⁡2(Δt/2)(Δ/2)2\frac{F_0^2}{8m}\frac{\sin^2(\Delta t/2)}{(\Delta/2)^2} になる。m=ω0=1m=\omega_0=1, F0=10−3F_0=10^{-3}, ω=0.99\omega=0.99, t=200t=200 で運動方程式を数値で解くと、エネルギーは 3.578×10−33.578\times10^{-3} で、この近似式の 3.54×10−33.54\times10^{-3} と約 1%1\% で一致する。量子の式 (4.12.13) と同じ関数である。

木下の勉強ステータス

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

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

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