tecktom.com

出てくるもの(2):光子との結合と、放出の n+1n+1

2026年10月15日9分で読めますAI生成教科書『虎の巻』第30章 つなぎ:場の量子化と光子
目次
  1. 結合の Hamiltonian
  2. 吸収は 、放出は
  3. 放出の が自然放出
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第30章 §30.4 ★★★ −d^⋅E^-\hat{\bm{d}}\cdot\hat{\bm{E}} に回転波近似を当てると ℏg0(σ+a^+σ−a^†)\hbar g_0(\sigma_+\hat a+\sigma_-\hat a^\dagger)。吸収の行列要素は n\sqrt n、放出は n+1\sqrt{n+1} で、+1+1 が自然放出。光子が多いと古典の Rabi 周波数に戻る

結合の Hamiltonian

物質のモデルは第18章の2準位原子 H^A=ℏω02σz\hat H_{\mathrm A}=\frac{\hbar\omega_0}{2}\sigma_z で、結合は双極子近似の −d^⋅E^(r0)-\hat{\bm{d}}\cdot\hat{\bm{E}}(\bm{r}_0)(r0\bm{r}_0 は原子の位置)である。第18章との違いは、E\bm{E} が古典の関数から式 (30.8) の演算子に替わることだけである。原子と共鳴に近い1つのモードだけを残し、deg\bm{d}_{eg} を実数でモードの偏光に平行、u(r0)\bm{u}(\bm{r}_0) を実数に選ぶと、Schrödinger 描像で

−d^⋅E^(r0)=−i dEph(r0) (σ++σ−)(a^−a^†),Eph(r0)=ℏω2ε0 ∣u(r0)∣-\hat{\bm{d}}\cdot\hat{\bm{E}}(\bm{r}_0)=-\mathrm{i}\,dE_{\mathrm{ph}}(\bm{r}_0)\,(\sigma_++\sigma_-)(\hat a-\hat a^\dagger),\qquad E_{\mathrm{ph}}(\bm{r}_0)=\sqrt{\frac{\hbar\omega}{2\varepsilon_0}}\,\lvert \bm{u}(\bm{r}_0)\rvert

になる(d=∣deg∣d=\lvert \bm{d}_{eg}\rvert)。右辺を展開した4つの項のうち、σ+a^\sigma_+\hat a(原子が上がり光子が1個減る)と σ−a^†\sigma_-\hat a^\dagger(原子が下がり光子が1個増える)は相互作用描像で e±i(ω0−ω)t\mathrm{e}^{\pm\mathrm{i}(\omega_0-\omega)t} とゆっくり回り、σ+a^†\sigma_+\hat a^\dagger と σ−a^\sigma_-\hat a は e±i(ω0+ω)t\mathrm{e}^{\pm\mathrm{i}(\omega_0+\omega)t} と速く回る。回転波近似で後の2つを捨て、モードの位相の原点を選び直して a^→ia^\hat a\to\mathrm{i}\hat a とすると

H^=ℏω02σz+ℏω a^†a^+ℏg0(σ+a^+σ−a^†),g0=dEph(r0)ℏ(30.18)\hat H=\frac{\hbar\omega_0}{2}\sigma_z+\hbar\omega\,\hat a^\dagger\hat a+\hbar g_0\bigl(\sigma_+\hat a+\sigma_-\hat a^\dagger\bigr),\qquad g_0=\frac{dE_{\mathrm{ph}}(\bm{r}_0)}{\hbar}\tag{30.18}

を得る(場の零点のエネルギーは定数なので省いた)。g0g_0 を結合定数と呼ぶ。これが第34章で扱う Jaynes–Cummings の模型の Hamiltonian である。

近似の裏付け

型 B(極限) | 小さいもの g0n+1/ω0→0g_0\sqrt{n+1}/\omega_0\to0 かつ ∣ω−ω0∣/ω0→0\lvert \omega-\omega_0\rvert /\omega_0\to0(速く回る σ+a^†\sigma_+\hat a^\dagger と σ−a^\sigma_-\hat a を捨てる) | 証明 要確認(古典の場では周期的な係数の常微分方程式に対する平均化の定理で条件付きに示せる。量子化した場の Rabi 模型で、光子数の上限を置かずに誤差を押さえる定理の形は原論文で確かめていない) | 行き先 近似の数理 第7章

吸収は n\sqrt n、放出は n+1\sqrt{n+1}

a^∣n⟩=n∣n−1⟩\hat a| n \rangle =\sqrt n| n-1 \rangle、a^†∣n⟩=n+1∣n+1⟩\hat a^\dagger| n \rangle =\sqrt{n+1}| n+1 \rangle から、式 (30.18) の結合の行列要素は

⟨e,n−1∣H^∣g,n⟩=ℏg0n(吸収),⟨g,n+1∣H^∣e,n⟩=ℏg0n+1(放出)(30.19)\langle e,n-1 | \hat H | g,n \rangle =\hbar g_0\sqrt n\quad(\text{吸収}),\qquad\langle g,n+1 | \hat H | e,n \rangle =\hbar g_0\sqrt{n+1}\quad(\text{放出})\tag{30.19}

である。結合は ∣e,n⟩| e,n \rangle と ∣g,n+1⟩| g,n+1 \rangle の2つの状態だけをつなぎ、その組の中で問題は第18章の 2×22\times2 の行列と同じ形になる。共鳴で ∣e,n⟩| e,n \rangle から出発すると

Pe(t)=cos⁡2(g0n+1 t),Ωn=2g0n+1(30.20)P_e(t)=\cos^2\bigl(g_0\sqrt{n+1}\,t\bigr),\qquad\Omega_n=2g_0\sqrt{n+1}\tag{30.20}

で、Rabi 振動の角周波数 Ωn\Omega_n は光子数で決まる。n=0n=0(真空)でも Ω0=2g0\Omega_0=2g_0 で振動し、これを真空 Rabi 振動と呼ぶ(第34章)。

古典の場の振幅 E0E_0 のモードのエネルギーは、時間平均で ε0E02Vm/2\varepsilon_0E_0^2V_{\mathrm m}/2(電場と磁場の分を合わせる)である。これを nℏωn\hbar\omega に等しいとおくと E0=2EphnE_0=2E_{\mathrm{ph}}\sqrt n で、古典の Rabi 周波数は Ω=dE0/ℏ=2g0n\Omega=dE_0/\hbar=2g_0\sqrt n になる。式 (30.20) の n+1\sqrt{n+1} と古典の n\sqrt n の比は 1+1/n\sqrt{1+1/n} で、光子が多いモードでは2つを区別することができない。第4・5部の半古典の扱いは、この n→∞n\to\infty の極限を使っていた(光子数が揺らいだレーザー光での正確な戻り方は第31章で扱う)。

放出の +1+1 が自然放出

光子数 nn のモードへの放出の確率は、行列要素の2乗 ∣ℏg0∣2(n+1)\lvert \hbar g_0\rvert ^2(n+1) に比例し、吸収は ∣ℏg0∣2n\lvert \hbar g_0\rvert ^2n に比例する。放出のうち nn に比例する分は、すでにある光子が原子を下へ押す誘導放出で、吸収と同じ係数を持つ。残りの 11 は光子が無くても起こる放出で、これが自然放出である。第20章で Einstein が熱平衡の条件から導いた3つの過程(B12B_{12}・B21B_{21}・A21A_{21})は、行列要素の n\sqrt n と n+1\sqrt{n+1} から1つの規則として出る。

原子が多くのモードに結合しているときは、原子の遷移の近くにある各モードへの速さを足す。モードの光子数の平均が nBn_{\mathrm B} の熱放射なら、

Wem=A21(nB(ω0)+1),Wabs=A21 nB(ω0)(30.21)W_{\mathrm{em}}=A_{21}\bigl(n_{\mathrm B}(\omega_0)+1\bigr),\qquad W_{\mathrm{abs}}=A_{21}\,n_{\mathrm B}(\omega_0)\tag{30.21}

になる(縮退の無い2準位)。ここで A21A_{21} は真空のモードに放出する速さで、すべてのモードについての和を第34章で実行する。u(ω)=ℏω3π2c3nBu(\omega)=\frac{\hbar\omega^3}{\pi^2c^3}n_{\mathrm B} を使うとB21u(ω0)=A21nBB_{21}u(\omega_0)=A_{21}n_{\mathrm B} は A21=ℏω03π2c3B21A_{21}=\frac{\hbar\omega_0^3}{\pi^2c^3}B_{21} と同じ式で、記法の正本と第20章の関係に一致する。第20章では Planck の分布を前提にこの関係を出したが、この章では逆に、行列要素の n+1n+1 と nn から熱平衡の占有の比が出る(例題 30.4)。

自然放出の光は、どのモードに出るかも位相も決まっていない。レーザーの発振のモードに入る自然放出の光子が、発振の線幅を決める(第22章の Schawlow–Townes の線幅)。場を古典の関数にすると n+1n+1 が nn になるので、半古典のレーザー理論ではこの線幅が出ない。

例題 30.4 n+1n+1 と nn から熱平衡を出す ★★

オリジナル

  1. (1) モードの光子数が幾何分布 p(n)=(1−q)qnp(n)=(1-q)q^n(q=e−ℏω/kBTq=\mathrm{e}^{-\hbar\omega/ k_{\mathrm B} T}、解析力学から様々な理論分野を組み立てる 第5部第3章)のとき、⟨n⟩=nB\langle n \rangle =n_{\mathrm B} と ⟨n+1⟩=nB+1\langle n+1 \rangle =n_{\mathrm B}+1 を示し、吸収と放出の速さがつり合う原子の占有の比 Pe/PgP_e/P_g を求めよ。共鳴(ω=ω0\omega=\omega_0)とする。
  2. (2) 例題 30.1 の光の微小共振器の腹に、d=2.5 ea0d=2.5\,ea_0 の2準位原子(遷移の波長 780 nm780\,\mathrm{nm})を置く。g0/2πg_0/2\pi を求めよ。
  3. (3) 光子数 nn の状態の Rabi 角周波数 2g0n+12g_0\sqrt{n+1} と、同じエネルギーの古典の場の 2g0n2g_0\sqrt n の違いが 1%1\% 未満になる最小の nn を求めよ。
解答

(1) ∑nn(1−q)qn=q/(1−q)=1/(q−1−1)=nB\sum_n n(1-q)q^n=q/(1-q)=1/(q^{-1}-1)=n_{\mathrm B}、したがって ⟨n+1⟩=nB+1\langle n+1 \rangle =n_{\mathrm B}+1。吸収の速さは Pg⟨n⟩P_g\langle n \rangle、放出の速さは Pe⟨n+1⟩P_e\langle n+1 \rangle に比例し、比例の係数は同じ ∣g0∣2\lvert g_0\rvert ^2 なので、つり合いの条件は

PePg=nBnB+1=1/(ex−1)ex/(ex−1)=e−x,x=ℏω0kBT\frac{P_e}{P_g}=\frac{n_{\mathrm B}}{n_{\mathrm B}+1}=\frac{1/(\mathrm{e}^{x}-1)}{\mathrm{e}^{x}/(\mathrm{e}^{x}-1)}=\mathrm{e}^{-x},\qquad x=\frac{\hbar\omega_0}{ k_{\mathrm B} T}

で、Boltzmann 因子になる。放出に +1+1 が無ければ比は 11 になり、原子は場と同じ温度にならない。

(2) Eph=3.83×103 V m−1E_{\mathrm{ph}}=3.83\times10^3\,\mathrm{V\,m^{-1}}、d=2.5×1.602×10−19×5.292×10−11=2.12×10−29 C md=2.5\times1.602\times10^{-19}\times5.292\times10^{-11}=2.12\times10^{-29}\,\mathrm{C\,m} で、g0=dEph/ℏ=7.69×108 s−1g_0=dE_{\mathrm{ph}}/\hbar=7.69\times10^{8}\,\mathrm{s^{-1}}、g0/2π=122 MHzg_0/2\pi=122\,\mathrm{MHz}。

(3) (n+1)/n<1.01\sqrt{(n+1)/n}<1.01 は 1+1/n<1.02011+1/n<1.0201、n>49.75n>49.75 なので、n=50n=50 から。光子が50個入れば、この模型の Rabi 振動は古典の場の計算と 1%1\% 以内で一致する。

木下の勉強ステータス

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

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

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