tecktom.com

当てる道具(1):双極子近似と、縮約した模型のゲージ

2026年10月15日9分で読めますAI生成教科書『虎の巻』第18章 2準位原子と双極子近似・回転波近似
目次
  1. 長さのゲージの Hamiltonian
  2. 捨てた項の大きさ
  3. 2準位に縮約するとゲージで結果が変わる
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第18章 §18.2 ★★★ 電場を原子の位置の値で置き換える。捨てた磁気双極子と電気四重極は αfs/2\alpha_{\mathrm{fs}}/2 と ka0/2ka_0/2 の桁。2準位に縮約すると A⋅p^\bm{A}\cdot\hat{\bm{p}} と −d^⋅E-\hat{\bm{d}}\cdot\bm{E} は共鳴でだけ一致する

長さのゲージの Hamiltonian

電子(電荷 −e-e)と光の場の最小結合は H^=(p^+eA(r^,t))2/(2m)+V(r^)\hat H=(\hat{\bm{p}}+e\bm{A}(\hat{\bm{r}},t))^2/(2m)+V(\hat{\bm{r}}) である(Coulomb ゲージで、光のスカラーポテンシャルは 00)。原子の大きさ a0a_0 が波長 λ\lambda よりはるかに小さいとき、A(r,t)\bm{A}(\bm{r},t) を原子の中心 r0\bm{r}_0 の値 A(r0,t)\bm{A}(\bm{r}_0,t) で置き換える。これが双極子近似で、小さいパラメータは ka0=2πa0/λka_0=2\pi a_0/\lambda である。置き換えたあとの A\bm{A} は位置に依らないので、波動関数の位相を ψ→e−ier⋅A(r0,t)/ℏψ\psi\to\mathrm{e}^{-\mathrm{i} e\bm{r}\cdot\bm{A}(\bm{r}_0,t)/\hbar}\psi と時刻ごとに変えるゲージ変換(Göppert-Mayer の変換)で A\bm{A} を消すことができ、

H^′=p^22m+V(r^)−d^⋅E(r0,t),E=−∂A∂t,d^=−er^(18.5)\hat H'=\frac{\hat{\bm{p}}^2}{2m}+V(\hat{\bm{r}})-\hat{\bm{d}}\cdot\bm{E}(\bm{r}_0,t),\qquad\bm{E}=-\frac{\partial \bm{A}}{\partial t},\quad\hat{\bm{d}}=-e\hat{\bm{r}}\tag{18.5}

になる。導出は 解析力学から様々な理論分野を組み立てる 第4部第12章 にあり、ここでは結果だけを使う。式 (18.5) を長さのゲージ、もとの A⋅p^\bm{A}\cdot\hat{\bm{p}} の形を速度のゲージと呼ぶ。2つの Hamiltonian はユニタリ変換で移り合うので、すべての準位を残して解く限り、物理量の値は同じになる。

捨てた項の大きさ

E(r,t)\bm{E}(\bm{r},t) を r0\bm{r}_0 のまわりで展開した次の項は、電場の空間変化と原子の電荷分布の結合(電気四重極)と、光の磁場 B\bm{B} と電子の軌道・スピンの磁気モーメントの結合(磁気双極子)である。大きさを、電気双極子の結合 ea0E0ea_0E_0 との比で見積もる。電気四重極は電場が原子の大きさの範囲で変わる割合 ka0ka_0 の分だけ小さく、展開の係数 1/21/2 を含めて比は ka0/2=πa0/λka_0/2=\pi a_0/\lambda の桁である。磁気双極子は、Bohr 磁子 μB=eℏ/(2m)\mu_{\mathrm B}=e\hbar/(2m) と光の磁場 B0=E0/cB_0=E_0/c の積なので、比は

μBB0ea0E0=ℏ2mca0=αfs2≈3.65×10−3(18.6)\frac{\mu_{\mathrm B}B_0}{ea_0E_0}=\frac{\hbar}{2mca_0}=\frac{\alpha_{\mathrm{fs}}}{2}\approx3.65\times10^{-3}\tag{18.6}

になる。ここで αfs=ℏ/(mca0)≈1/137\alpha_{\mathrm{fs}}=\hbar/(mca_0)\approx1/137 は微細構造定数で、分極率 α\alpha と区別するために添字 fs を付けた。遷移の速さは行列要素の2乗に比例するので、許された電気双極子遷移に比べて磁気双極子遷移は 10−510^{-5}、可視光の電気四重極遷移は 10−710^{-7} の桁で遅い。deg=0\bm{d}_{eg}=0 になる遷移(同じ偶奇性の準位の間など)では、これらの項が主になる。光格子時計の時計の遷移や、水素の 1s→2s1s\to2s の2光子遷移は、電気双極子で禁じられた遷移を使うことで線幅を極めて狭くしたものである(第24章・第28章、第24章・第28章)。

2準位に縮約するとゲージで結果が変わる

2つのゲージがユニタリ変換で移り合うことは、すべての準位を残したときにだけ成り立つ。変換 e−ier⋅A/ℏ\mathrm{e}^{-\mathrm{i} e\bm{r}\cdot\bm{A}/\hbar} は位置の演算子の関数で、2準位の部分空間を部分空間のまま保たないからである。2準位に縮約した模型で比べてみる。速度のゲージで A\bm{A} の1次の項は (e/m)A(r0,t)⋅p^(e/m)\bm{A}(\bm{r}_0,t)\cdot\hat{\bm{p}} で(Coulomb ゲージでは p^⋅A=A⋅p^\hat{\bm{p}}\cdot\bm{A}=\bm{A}\cdot\hat{\bm{p}})、運動量の行列要素は交換関係 [H^A, r^]=−iℏp^/m\left[\hat H_{\mathrm A},\,\hat{\bm{r}}\right]=-\mathrm{i}\hbar\hat{\bm{p}}/m から

⟨e∣p^∣g⟩=imω0⟨e∣r^∣g⟩(18.7)\langle e | \hat{\bm{p}} | g \rangle =\mathrm{i} m\omega_0\langle e | \hat{\bm{r}} | g \rangle\tag{18.7}

になる(解析力学から様々な理論分野を組み立てる 第4部第12章 と同じ式)。電場 E=Re⁡[E0e−iωt]\bm{E}=\operatorname{Re}[\bm{E}_0\mathrm{e}^{-\mathrm{i}\omega t}] に対するベクトルポテンシャルは A=Re⁡[A0e−iωt]\bm{A}=\operatorname{Re}[\bm{A}_0\mathrm{e}^{-\mathrm{i}\omega t}]、A0=E0/(iω)\bm{A}_0=\bm{E}_0/(\mathrm{i}\omega) なので、e−iωt\mathrm{e}^{-\mathrm{i}\omega t} の項の行列要素は

em A0⋅⟨e∣p^∣g⟩=em⋅E0iω⋅imω0⟨e∣r^∣g⟩=−ω0ω deg⋅E0(18.8)\frac em\,\bm{A}_0\cdot\langle e | \hat{\bm{p}} | g \rangle =\frac em\cdot\frac{\bm{E}_0}{\mathrm{i}\omega}\cdot\mathrm{i} m\omega_0\langle e | \hat{\bm{r}} | g \rangle =-\frac{\omega_0}{\omega}\,\bm{d}_{eg}\cdot\bm{E}_0\tag{18.8}

になる。長さのゲージの −deg⋅E0-\bm{d}_{eg}\cdot\bm{E}_0 と比べて ω0/ω\omega_0/\omega 倍で、2つは共鳴 ω=ω0\omega=\omega_0 でだけ一致する。共鳴から離れた光では、縮約した2準位の模型の Rabi 周波数が、どちらのゲージで縮約したかで違う。

どちらが正しいかは、捨てた準位の寄与まで含めた計算と比べて決まる。速度のゲージには A\bm{A} の2次の項 e2A2/(2m)e^2\bm{A}^2/(2m) もあり、これは双極子近似のもとですべての準位を同じだけずらす。この項と、ほかの準位を通る2次の寄与を足すと、Thomas–Reiche–Kuhn の総和則によって長さのゲージの結果に戻ることを示すことができる。長さのゲージの模型では、捨てた準位の寄与は遷移の振動子強度で決まるので、共鳴に近い遷移が1つの振動子強度の大部分を持っていれば小さい。このため、量子光学では2準位の模型を長さのゲージで作るのが標準になっている。この巻もそれに従う。共鳴の近くでは ω0/ωL=1−Δ/ωL\omega_0/\omega_{\mathrm L}=1-\Delta/\omega_{\mathrm L} なので、違いは離調の割合 ∣Δ∣/ω0\lvert \Delta\rvert /\omega_0 の桁にとどまる。

例題 18.2 運動量の行列要素と、ゲージの違いの大きさ ★★

オリジナル例題 18.1 の水素の波動関数を使う(原子単位、ℏ=m=a0=1\hbar=m=a_0=1、ω0=3/8\omega_0=3/8)。

  1. (1) ⟨2p,0∣∂z∣1s⟩=∫ψ2p,0 ∂zψ1s d3r\langle 2p,0 | \partial_z | 1s \rangle =\int\psi_{2p,0}\,\partial_z\psi_{1s}\,\mathrm{d}^3r を計算し、式 (18.7) が成り立つことを確かめよ。
  2. (2) 遷移の波長が 780 nm780\,\mathrm{nm} の原子に波長 1064 nm1064\,\mathrm{nm} の光を当てる(第1章の総合問題の光トラップ)。2準位の模型を速度のゲージで作ると、長さのゲージで作った模型に比べて Rabi 周波数は何倍になるか。光シフト(Ω2\Omega^2 に比例する)は何倍になるか。
  3. (3) (2) と、式 (18.6) の磁気双極子の比を比べ、遠い離調の光で縮約の誤差がどれほど大きいかを述べよ。
解答

(1) ∂ze−r=−(z/r)e−r=−cos⁡θ e−r\partial_z\mathrm{e}^{-r}=-(z/r)\mathrm{e}^{-r}=-\cos\theta\,\mathrm{e}^{-r} なので

⟨2p,0∣∂z∣1s⟩=−1π⋅142π∫0∞r3e−3r/2dr∫cos⁡2θ dΩ=−142 π⋅6(3/2)4⋅4π3=−16281.\langle 2p,0 | \partial_z | 1s \rangle =-\frac{1}{\sqrt\pi}\cdot\frac{1}{4\sqrt{2\pi}}\int_0^\infty r^3\mathrm{e}^{-3r/2}\mathrm{d} r\int\cos^2\theta\,\mathrm{d}\Omega=-\frac{1}{4\sqrt2\,\pi}\cdot\frac{6}{(3/2)^4}\cdot\frac{4\pi}{3}=-\frac{16\sqrt2}{81}.

一方 C=1282/243C=128\sqrt2/243 に ω0=3/8\omega_0=3/8 を掛けると ω0C=162/81\omega_0C=16\sqrt2/81 で、⟨2p,0∣∂z∣1s⟩=−ω0C\langle 2p,0 | \partial_z | 1s \rangle =-\omega_0C。p^z=−i∂z\hat p_z=-\mathrm{i}\partial_z なので ⟨2p,0∣p^z∣1s⟩=iω0C\langle 2p,0 | \hat p_z | 1s \rangle =\mathrm{i}\omega_0C となり、式 (18.7)(原子単位で m=1m=1)と一致する。

(2) 式 (18.8) から、Rabi 周波数は ω0/ωL=λL/λ0=1064/780=1.364\omega_0/\omega_{\mathrm L}=\lambda_{\mathrm L}/\lambda_0=1064/780=1.364 倍、光シフトはその2乗の 1.8611.861 倍になる。

(3) 磁気双極子のように原子の構造に由来する補正は 10−310^{-3} の桁だが、ゲージの選び方による違いは 86%86\% で、比べものにならないほど大きい。共鳴から 27%27\% も離れた光では、「2つの準位だけを残す」という縮約そのものが不確かになり、どのゲージで縮約したかが結果を左右する。この場合は、ほかの準位も足した分極率(第11章の振動子強度の和、第11章)で計算するのが確かである。

近似の裏付け

型 B(極限) | 小さいもの ka0=2πa0/λ→0ka_0=2\pi a_0/\lambda\to0(原子が波長より十分小さい) | 証明 要確認(多重極展開の第1項。次の項の大きさは ka0/2ka_0/2 と αfs/2\alpha_{\mathrm{fs}}/2 の桁と見積もったが、打ち切りの誤差を評価する定理の形は原論文で確かめていない) | 行き先 近似の数理 第7章

ほかの巻では

木下の勉強ステータス

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

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

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