tecktom.com

結合と近似:双極子近似と回転波近似

2026年10月15日10分で読めますAI生成教科書『虎の巻』第1章 序章:モードと2準位系で光を読む
目次
  1. 双極子近似
  2. 回転波近似
  3. 何を近似したか
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第1章 §1.3 ★★★ 結合 −d^⋅E-\hat{\bm{d}}\cdot\bm{E} に2つの B 型の近似(ka0→0ka_0\to0 と Ω/ω0→0\Omega/\omega_0\to0)を当てると、2準位系と光の問題は 2×22\times2 の時間に依らない行列になる

双極子近似

原子の大きさは Bohr 半径 a0=5.29×10−11 ma_0=5.29\times10^{-11}\,\mathrm m の桁で、可視光の波長 λ≈500 nm\lambda\approx500\,\mathrm{nm} より4桁小さい。光の電場 E(r,t)=Re⁡[E0ei(k⋅r−ωt)]\bm{E}(\bm{r},t)=\operatorname{Re}[\bm{E}_0\mathrm{e}^{\mathrm{i}(\bm{k}\cdot\bm{r}-\omega t)}] は、原子の大きさの範囲ではほぼ一定とみなせる。原子の中心 r0\bm{r}_0 のまわりで eik⋅r=eik⋅r0[1+ik⋅(r−r0)+⋯ ]\mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}=\mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}_0}[1+\mathrm{i}\bm{k}\cdot(\bm{r}-\bm{r}_0)+\cdots] と展開して第1項だけを残すと、原子と光の相互作用は

H^int=−d^⋅E(r0,t)(1.11)\hat H_{\mathrm{int}}=-\hat{\bm{d}}\cdot\bm{E}(\bm{r}_0,t)\tag{1.11}

になる。これが双極子近似で、小さいパラメータは ka0=2πa0/λka_0=2\pi a_0/\lambda である。展開の第2項は磁気双極子と電気四重極の相互作用にあたり、deg=0\bm{d}_{eg}=0 になる遷移(選択則で禁制の遷移)では第2項が主になる。最小結合 (p^+eA)2/(2m)(\hat{\bm{p}}+e\bm{A})^2/(2m) から式 (1.11) への書き換え(Göppert-Mayer のゲージ変換)は解析力学から様々な理論分野を組み立てる 第4部第12章 にある。

回転波近似

2準位系に式 (1.11) を入れると、Hamiltonian は H^=ℏω02σz−ℏΩcos⁡(ωLt)(σ++σ−)\hat H=\frac{\hbar\omega_0}{2}\sigma_z-\hbar\Omega\cos(\omega_{\mathrm L}t)(\sigma_++\sigma_-) になる。ここで ωL\omega_{\mathrm L} はレーザーの角周波数、Ω=deg⋅E0/ℏ\Omega=\bm{d}_{eg}\cdot\bm{E}_0/\hbar は Rabi 周波数である。cos⁡(ωLt)=12(e−iωLt+eiωLt)\cos(\omega_{\mathrm L}t)=\frac12(\mathrm{e}^{-\mathrm{i}\omega_{\mathrm L}t}+\mathrm{e}^{\mathrm{i}\omega_{\mathrm L}t}) と分けると、結合は4つの項 σ+e−iωLt\sigma_+\mathrm{e}^{-\mathrm{i}\omega_{\mathrm L}t}・σ−eiωLt\sigma_-\mathrm{e}^{\mathrm{i}\omega_{\mathrm L}t}・σ+eiωLt\sigma_+\mathrm{e}^{\mathrm{i}\omega_{\mathrm L}t}・σ−e−iωLt\sigma_-\mathrm{e}^{-\mathrm{i}\omega_{\mathrm L}t} を含む。原子の自由な時間発展では σ±\sigma_\pm が e±iω0t\mathrm{e}^{\pm\mathrm{i}\omega_0t} で回るので、前の2つは e±i(ω0−ωL)t\mathrm{e}^{\pm\mathrm{i}(\omega_0-\omega_{\mathrm L})t} でゆっくり回り、後の2つは e±i(ω0+ωL)t\mathrm{e}^{\pm\mathrm{i}(\omega_0+\omega_{\mathrm L})t} で速く回る。速く回る項は、Rabi 周波数の時間尺度 1/Ω1/\Omega で平均すると O(Ω/ω0)\mathcal{O}(\Omega/\omega_0) の寄与しか残さない。これを捨てるのが回転波近似(RWA)である。

回転波近似のあとで、角周波数 ωL\omega_{\mathrm L} で回る座標(eiωLtσz/2\mathrm{e}^{\mathrm{i}\omega_{\mathrm L}t\sigma_z/2} で状態を変換した座標)に移ると、時間依存が消えて

H^=−ℏΔ2σz−ℏΩ2σx,Δ=ωL−ω0(1.12)\hat H=-\frac{\hbar\Delta}{2}\sigma_z-\frac{\hbar\Omega}{2}\sigma_x,\qquad\Delta=\omega_{\mathrm L}-\omega_0\tag{1.12}

になる(この本の流儀では離調 Δ\Delta はレーザー −- 原子で、Δ<0\Delta<0 が赤い離調)。2×22\times2 の時間に依らない行列なので、時間発展は対角化で求まる。t=0t=0 に ∣g⟩| g \rangle にあった原子が ∣e⟩| e \rangle に見つかる確率は

Pe(t)=Ω2Ω′2sin⁡2Ω′t2,Ω′=Ω2+Δ2(1.13)P_e(t)=\frac{\Omega^2}{\Omega'^2}\sin^2\frac{\Omega't}{2},\qquad\Omega'=\sqrt{\Omega^2+\Delta^2}\tag{1.13}

で、共鳴(Δ=0\Delta=0)では sin⁡2(Ωt/2)\sin^2(\Omega t/2) で上下の準位を行き来する。これが Rabi 振動で、第18・19章の主題である(第18章・第19章)。同じ 2×22\times2 の行列は磁場中のスピンの歳差と同じ形をしている(解析力学から様々な理論分野を組み立てる 第4部第5章)。

図 1.1 は、共鳴で Ω/ω0=0.1\Omega/\omega_0=0.1 のとき、回転波近似の式 (1.13) と、速く回る項を残したまま Schrödinger 方程式を数値で解いた結果を比べたものである。厳密な解は回転波近似の曲線のまわりで周期 π/ω0\pi/\omega_0 の小さな振動をし、山の位置も少しずれる。Ω/ω0=0.1\Omega/\omega_0=0.1 でずれは最大 0.0240.024 で、実際の原子と可視光ではΩ/ω0\Omega/\omega_0 は 10−810^{-8} の桁(例題 1.3)なので、ずれは図に描けないほど小さい。

図

図 1.1 共鳴(Δ=0\Delta=0)で Ω/ω0=0.1\Omega/\omega_0=0.1 の2準位系の Rabi 振動。実線は回転波近似の式 (1.13)、点線は速く回る項を残した H^=ℏω02σz−ℏΩcos⁡(ω0t)σx\hat H=\frac{\hbar\omega_0}{2}\sigma_z-\hbar\Omega\cos(\omega_0t)\sigma_x の Schrödinger 方程式を数値で解いた結果である。点線は実線のまわりで角周波数 2ω02\omega_0 の小さな振動をする。

何を近似したか

この節までに置いた近似は3つある。1つ目は原子の多くの準位から2つだけを残したこと、2つ目は電場を原子の位置の値で置き換えたこと(双極子近似)、3つ目は速く回る項を捨てたこと(回転波近似)である。さらに、この章では光を古典の電場 E(t)\bm{E}(t) として扱った(半古典の扱い)。光を量子化したモードとして扱う必要が出るのは第6部で、そこでは式 (1.12) の Ω\Omega が光子の演算子に置き換わる(第34章)。

近似の裏付け

型 A(射影) | 小さいもの ∣g⟩| g \rangle と ∣e⟩| e \rangle が張る2次元の部分空間(ほかの準位への離調が Ω\Omega と Δ\Delta より十分大きい) | 証明 要確認(ほかの準位への漏れの大きさを評価する定理の形は確かめていない) | 行き先 近似の数理 第4章

近似の裏付け

型 B(極限) | 小さいもの ka0=2πa0/λ→0ka_0=2\pi a_0/\lambda\to0(原子が波長より十分小さい) | 証明 要確認(多重極展開の第1項。展開の打ち切りの誤差の評価は原論文で確かめていない) | 行き先 近似の数理 第7章

近似の裏付け

型 B(極限) | 小さいもの Ω/ω0→0\Omega/\omega_0\to0 かつ ∣Δ∣/ω0→0\lvert \Delta\rvert /\omega_0\to0(速く回る項を時間平均で捨てる) | 証明 条件付き(時間 1/Ω1/\Omega の桁の範囲で誤差が O(Ω/ω0)\mathcal{O}(\Omega/\omega_0)。周期的な係数の常微分方程式に対する平均化の定理による。強い結合では Bloch–Siegert のずれが残る) | 行き先 近似の数理 第7章

例題 1.3 弱いレーザー光の Rabi 周波数と2つの小さいパラメータ ★★

オリジナル遷移の波長 λ0=780 nm\lambda_0=780\,\mathrm{nm}、双極子の行列要素 d=2.5 ea0d=2.5\,ea_0 の2準位原子(与えられた値として使う)に、強度 I=10 mW cm−2I=10\,\mathrm{mW\,cm^{-2}} の共鳴の光を当てる。強度と電場の振幅の関係は I=12cε0E02I=\frac12c\varepsilon_0E_0^2(真空)とする。

  1. (1) E0E_0 と Rabi 周波数 Ω\Omega を求めよ。Ω/2π\Omega/2\pi を MHz で答えよ。
  2. (2) 回転波近似の小さいパラメータ Ω/ω0\Omega/\omega_0 と、双極子近似の小さいパラメータ ka0ka_0 を求めよ。
  3. (3) 原子を ∣g⟩| g \rangle から ∣e⟩| e \rangle へ移しきる光のパルス(π\pi パルス)の長さを求めよ。
解答

(1) I=100 W m−2I=100\,\mathrm{W\,m^{-2}} なので E0=2I/(cε0)=200/(2.998×108×8.854×10−12)≈274 V m−1E_0=\sqrt{2I/(c\varepsilon_0)}=\sqrt{200/(2.998\times10^8\times8.854\times10^{-12})}\approx274\,\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}\approx2.12\times10^{-29}\,\mathrm{C\,m} で、Ω=dE0/ℏ=2.12×10−29×274/1.0546×10−34≈5.52×107 s−1\Omega=dE_0/\hbar=2.12\times10^{-29}\times274/1.0546\times10^{-34}\approx5.52\times10^7\,\mathrm{s^{-1}}、Ω/2π≈8.78 MHz\Omega/2\pi\approx8.78\,\mathrm{MHz} である。

(2) ω0=2πc/λ0≈2.41×1015 s−1\omega_0=2\pi c/\lambda_0\approx2.41\times10^{15}\,\mathrm{s^{-1}} なので Ω/ω0≈2.3×10−8\Omega/\omega_0\approx2.3\times10^{-8}。ka0=2π×5.29×10−11/7.8×10−7≈4.3×10−4ka_0=2\pi\times5.29\times10^{-11}/7.8\times10^{-7}\approx4.3\times10^{-4}。どちらも小さく、2つの近似はこの条件で十分に成り立つ。回転波近似の小さいパラメータは双極子近似のものより4桁小さい。

(3) 共鳴で Pe=sin⁡2(Ωt/2)P_e=\sin^2(\Omega t/2) が 11 になるのは Ωt=π\Omega t=\pi のときで、t=π/Ω≈57 nst=\pi/\Omega\approx57\,\mathrm{ns} である。光の周期 2.6 fs2.6\,\mathrm{fs} の約 2×1072\times10^7 倍の時間をかけて、原子はゆっくりと上の準位へ移る。

木下の勉強ステータス

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

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

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