tecktom.com

出てくるもの(2):位相速度・群速度と遅い光

2026年10月15日15分で読めますAI生成教科書『虎の巻』第11章 Lorentz 振動子と局所場
目次
  1. パルスの包絡線の速さ
  2. 吸収の線の中の群速度
  3. 結合した2つの振動子と遅い光
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第11章 §11.4 ★★★ パルスの包絡線は群屈折率 ng=n+ω dn/dωn_{\mathrm g}=n+\omega\,\mathrm{d} n/\mathrm{d}\omega で決まる速さで進む。2つの振動子を結合させると、吸収の谷の中に急な正常分散が開き、ngn_{\mathrm g} は ωp2/Ωc2\omega_{\mathrm p}^2/\Omega_{\mathrm c}^2 まで大きくなる

パルスの包絡線の速さ

平面波の位相は ωt−kz\omega t-k z が一定の面で、その速さは位相速度 vp=ω/k=c/nv_{\mathrm p}=\omega/k=c/n である。ところが情報やエネルギーを運ぶのは、有限の長さを持つパルスである。中心の角周波数 ωc\omega_{\mathrm c} のまわりに広がった波を

E(z,t)=∫dω2π E~(ω) ei(k(ω)z−ωt),k(ω)=n(ω)ωcE(z,t)=\int\frac{\mathrm{d}\omega}{2\pi}\,\tilde E(\omega)\,\mathrm{e}^{\mathrm{i}(k(\omega)z-\omega t)},\qquad k(\omega)=\frac{n(\omega)\omega}{c}

と書き(この巻の Fourier の流儀)、吸収を無視して k(ω)k(\omega) を ωc\omega_{\mathrm c} のまわりで k≈kc+k′(ω−ωc)+12k′′(ω−ωc)2k\approx k_{\mathrm c}+k'(\omega-\omega_{\mathrm c})+\frac12k''(\omega-\omega_{\mathrm c})^2 と展開する。2次の項を落とすと

E(z,t)≈ei(kcz−ωct) E(t−k′z),E(τ)=∫dω2π E~(ω) e−i(ω−ωc)τE(z,t)\approx\mathrm{e}^{\mathrm{i}(k_{\mathrm c}z-\omega_{\mathrm c}t)}\,\mathcal E(t-k'z),\qquad\mathcal E(\tau)=\int\frac{\mathrm{d}\omega}{2\pi}\,\tilde E(\omega)\,\mathrm{e}^{-\mathrm{i}(\omega-\omega_{\mathrm c})\tau}

で、包絡線 E\mathcal E は形を変えずに速さ 1/k′1/k' で進む。この速さを群速度と呼ぶ。k′=d(nω/c)/dωk'=\mathrm{d}(n\omega/c)/\mathrm{d}\omega を計算して

vg=(dkdω)−1=cng,ng=n+ωdndω=n−λdndλ(11.13)v_{\mathrm g}=\Bigl(\frac{\mathrm{d} k}{\mathrm{d}\omega}\Bigr)^{-1}=\frac{c}{n_{\mathrm g}},\qquad n_{\mathrm g}=n+\omega\frac{\mathrm{d} n}{\mathrm{d}\omega}=n-\lambda\frac{\mathrm{d} n}{\mathrm{d}\lambda}\tag{11.13}

を得る。ngn_{\mathrm g} を群屈折率と呼ぶ(λ\lambda は真空中の波長で、ω d/dω=−λ d/dλ\omega\,\mathrm{d}/\mathrm{d}\omega=-\lambda\,\mathrm{d}/\mathrm{d}\lambda)。正常分散の範囲では dn/dω>0\mathrm{d} n/\mathrm{d}\omega>0 なので ng>nn_{\mathrm g}>n で、パルスは位相より遅れて進む。

落とした2次の項 k′′=d2k/dω2k''=\mathrm{d}^2k/\mathrm{d}\omega^2 は群速度分散(GVD)と呼ばれ、周波数ごとに群速度が違うためにパルスが広がる度合いを表す。波長で書くと

k′′=λ32πc2d2ndλ2(11.14)k''=\frac{\lambda^3}{2\pi c^2}\frac{\mathrm{d}^2n}{\mathrm{d}\lambda^2}\tag{11.14}

である。包絡線が exp⁡(−t2/(2τ02))\exp(-t^2/(2\tau_0^2)) の Gauss 型のパルスは、距離 zz 進むと幅が τ01+(k′′z/τ02)2\tau_0\sqrt{1+(k''z/\tau_0^2)^2} に広がる。

図

図 11.3 光学ガラス N-BK7 の屈折率 nn(実線)と群屈折率 ngn_{\mathrm g}(破線)。SCHOTT 社のデータシートの Sellmeier の係数(式 (11.18))から計算した。可視域では紫外の共鳴に近づくほど nn も ng−nn_{\mathrm g}-n も大きくなる。

例題 11.5 ガラスの中のフェムト秒パルス ★★

オリジナル光学ガラス N-BK7 の屈折率は、波長 800 nm800\,\mathrm{nm} で n=1.5108n=1.5108、dn/dλ=−0.01984 μm−1\mathrm{d} n/\mathrm{d}\lambda=-0.01984\,\mu\mathrm m^{-1}、d2n/dλ2=0.04925 μm−2\mathrm{d}^2n/\mathrm{d}\lambda^2=0.04925\,\mu\mathrm m^{-2} である(式 (11.18) の Sellmeier の式から計算した値)。

  1. (1) 群屈折率 ngn_{\mathrm g} を求め、厚さ 10 mm10\,\mathrm{mm} の板を通るときの、位相の遅れの時間と包絡線の遅れの時間を比べよ。
  2. (2) 群速度分散 k′′k'' を fs2/mm\mathrm{fs^2/mm} の単位で求めよ。
  3. (3) τ0=10 fs\tau_0=10\,\mathrm{fs} の Gauss 型のパルスがこの板を通ると、幅は何倍になるか。
解答

(1) 式 (11.13) から ng=1.5108+0.8×0.01984=1.5267n_{\mathrm g}=1.5108+0.8\times0.01984=1.5267。板の通過時間は位相が Ln/c=50.39 psLn/c=50.39\,\mathrm{ps}、包絡線が Lng/c=50.92 psLn_{\mathrm g}/c=50.92\,\mathrm{ps} で、包絡線は位相より 0.53 ps0.53\,\mathrm{ps} 遅れて出てくる。この差は、波長 800 nm800\,\mathrm{nm} の光の周期 2.67 fs2.67\,\mathrm{fs} の約 200200 倍にあたる。

(2) 式 (11.14) で λ3 d2n/dλ2=0.512 μm3×0.04925 μm−2=2.522×10−8 m\lambda^3\,\mathrm{d}^2n/\mathrm{d}\lambda^2=0.512\,\mu\mathrm m^3\times0.04925\,\mu\mathrm m^{-2}=2.522\times10^{-8}\,\mathrm m なので

k′′=2.522×10−8 m2π×(2.998×108 m/s)2=4.47×10−26 s2/m=44.7 fs2/mm.k''=\frac{2.522\times10^{-8}\,\mathrm m}{2\pi\times(2.998\times10^8\,\mathrm{m/s})^2}=4.47\times10^{-26}\,\mathrm{s^2/m}=44.7\,\mathrm{fs^2/mm}.

k′′>0k''>0 で、赤い成分が青い成分より速く進む(正常分散の範囲の群速度分散)。

(3) k′′z=44.7×10=447 fs2k''z=44.7\times10=447\,\mathrm{fs^2}、k′′z/τ02=4.47k''z/\tau_0^2=4.47 なので、幅は 1+4.472=4.58\sqrt{1+4.47^2}=4.58 倍の 45.8 fs45.8\,\mathrm{fs} になる。(1) の遅れ 0.53 ps0.53\,\mathrm{ps} はパルスの位置をずらすだけだが、k′′k'' は形を変える。10 fs10\,\mathrm{fs} のパルスは帯域が広いので、1 cm1\,\mathrm{cm} のガラスで幅が4倍を超える。短いパルスのレーザーで、プリズムの対や分散を持つ鏡で負の k′′k'' を足して打ち消すのはこのためで、第23章(第23章)で扱う。

吸収の線の中の群速度

正常分散の範囲から共鳴に近づくと ngn_{\mathrm g} は大きくなる。ところが共鳴の中心の異常分散の範囲では dn/dω<0\mathrm{d} n/\mathrm{d}\omega<0 で、例題 11.3 の蒸気では中心の傾きから ng≈ω0 dn/dω=−ωp2/γ2=−8.1×104n_{\mathrm g}\approx\omega_0\,\mathrm{d} n/\mathrm{d}\omega=-\omega_{\mathrm p}^2/\gamma^2=-8.1\times10^4 になる。群速度が真空の光速より大きいか負になるが、これは信号が光速を超えることを意味しない。この範囲では光が 59 μm59\,\mu\mathrm m ごとに 1/e1/\mathrm{e} に吸収され、パルスの前の部分より後ろの部分が強く削られて、包絡線のピークが前へずれて見えるだけである。式 (11.13) を導くときに吸収と k′′k'' を落としたので、この範囲ではそもそも包絡線が形を保たない。信号の先端が光速 cc を超えないことは、χ(ω)\chi(\omega) が上半面で解析的であること(因果律)から言え、第12章で扱う。

結合した2つの振動子と遅い光

群屈折率を大きくするには、吸収を小さく保ったまま dn/dω\mathrm{d} n/\mathrm{d}\omega を大きくすればよい。1つの Lorentz 振動子ではこの2つは両立しない。急な分散は共鳴の中心にしかなく、そこは吸収も最大だからである。2つ目の振動子を結合させると、この2つが両立する。光が直接駆動する振動子1(減衰 γ1\gamma_1)に、光と結合しない振動子2(減衰 γ2≪γ1\gamma_2\ll\gamma_1)を、同じ固有角周波数 ω0\omega_0 で結合させる:

x¨1+γ1x˙1+ω02x1=−emEloc+ω0Ωcx2,x¨2+γ2x˙2+ω02x2=ω0Ωcx1.(11.15)\ddot x_1+\gamma_1\dot x_1+\omega_0^2x_1=-\frac{e}{m}E_{\mathrm{loc}}+\omega_0\Omega_{\mathrm c}x_2,\qquad\ddot x_2+\gamma_2\dot x_2+\omega_0^2x_2=\omega_0\Omega_{\mathrm c}x_1 .\tag{11.15}

Ωc\Omega_{\mathrm c} は結合の強さを角周波数で表したものである。定常解を求めて、共鳴の近くで ω02−ω2−iγkω≈2ω0(ω0−ω−iγk/2)\omega_0^2-\omega^2-\mathrm{i}\gamma_k\omega\approx2\omega_0(\omega_0-\omega-\mathrm{i}\gamma_k/2) とおき、離調 δ=ω−ω0\delta=\omega-\omega_0 を使うと、薄い気体の感受率は

χ(δ)=ωp22ω0 1−δ−iγ1/2+Ωc2/4δ+iγ2/2(11.16)\chi(\delta)=\frac{\omega_{\mathrm p}^2}{2\omega_0}\,\frac{1}{-\delta-\mathrm{i}\gamma_1/2+\dfrac{\Omega_{\mathrm c}^2/4}{\delta+\mathrm{i}\gamma_2/2}}\tag{11.16}

になる(Ωc=0\Omega_{\mathrm c}=0 で式 (11.11) に戻る)。

γ2→0\gamma_2\to0 とすると、分母の第3項は δ=0\delta=0 で発散するので、χ(0)=0\chi(0)=0 になる。共鳴の中心で吸収も屈折も消え、媒質は透明になる。振動子1を駆動する力と、振動子2から結合を通して戻ってくる力が、δ=0\delta=0 でちょうど打ち消すからである。δ\delta が小さいとき分母は Ωc2/(4δ)\Omega_{\mathrm c}^2/(4\delta) に近いので χ≈2ωp2δ/(ω0Ωc2)\chi\approx2\omega_{\mathrm p}^2\delta/(\omega_0\Omega_{\mathrm c}^2) で、n−1=χ/2n-1=\chi/2 の傾きと式 (11.13) から

ng≈ω0dndω=ωp2Ωc2(11.17)n_{\mathrm g}\approx\omega_0\frac{\mathrm{d} n}{\mathrm{d}\omega}=\frac{\omega_{\mathrm p}^2}{\Omega_{\mathrm c}^2}\tag{11.17}

を得る。透明な窓の中の分散は正常分散で、その傾きは結合を弱くするほど急になる。窓の幅は、分母の Ωc2/(4∣δ∣)\Omega_{\mathrm c}^2/(4\lvert \delta\rvert ) が γ1/2\gamma_1/2 と同じになる ∣δ∣≈Ωc2/(2γ1)\lvert \delta\rvert \approx\Omega_{\mathrm c}^2/(2\gamma_1) で決まり、結合を弱くするほど狭くなる。γ2≠0\gamma_2\ne0 では中心の吸収は完全には消えず、結合が無い場合の 1/(1+Ωc2/(γ1γ2))1/(1+\Omega_{\mathrm c}^2/(\gamma_1\gamma_2)) 倍に下がる。図 11.4 が式 (11.16) の形である。

図

図 11.4 結合した2つの振動子の感受率 (11.16)(Ωc=0.8 γ1\Omega_{\mathrm c}=0.8\,\gamma_1、γ2=0.01 γ1\gamma_2=0.01\,\gamma_1)。破線が Im⁡χ\operatorname{Im}\chi(吸収)、実線が Re⁡χ\operatorname{Re}\chi(屈折)、点線は結合が無いときの Im⁡χ\operatorname{Im}\chi。中心に吸収の谷が開き、その中で Re⁡χ\operatorname{Re}\chi が急な正の傾きを持つ。

これは量子力学の3準位原子で起こる電磁誘起透明化(EIT)の古典版で、結合した2つの振動子の模型は C. L. Garrido Alzar, M. A. G. Martinez, P. Nussenzveig (2002)、Am. J. Phys. 70, 37 が詳しく調べている。原子では、振動子2の役を、光で直接つながらない別の準位が担い、Ωc\Omega_{\mathrm c} は別のレーザー(結合光)の Rabi 周波数である。量子論の取り扱いは第19章(第19章)の光 Bloch 方程式の上で行う。L. V. Hau, S. E. Harris, Z. Dutton, C. H. Behroozi (1999)、Nature 397, 594 は、ナトリウム原子の極低温の気体でこの窓を使い、光パルスの群速度を 17 m/s17\,\mathrm{m/s}(ng≈1.8×107n_{\mathrm g}\approx1.8\times10^7)まで下げた。

例題 11.6 遅い光の窓と遅延時間 ★★★

オリジナル共鳴の波長の近くで ωp=3.99×1010 s−1\omega_{\mathrm p}=3.99\times10^{10}\,\mathrm{s^{-1}}(N=1.0×1018 m−3N=1.0\times10^{18}\,\mathrm{m^{-3}}、f=0.5f=0.5 に相当)、γ1=2π×10 MHz\gamma_1=2\pi\times10\,\mathrm{MHz}、γ2→0\gamma_2\to0 の媒質を考える。

  1. (1) 群速度を 17 m/s17\,\mathrm{m/s} にする Ωc\Omega_{\mathrm c} を求めよ。
  2. (2) そのときの窓の半幅 Ωc2/(2γ1)\Omega_{\mathrm c}^2/(2\gamma_1) を周波数(Hz\mathrm{Hz})で求めよ。
  3. (3) 長さ L=200 μmL=200\,\mu\mathrm m の媒質での遅延時間 τ=Lng/c\tau=Ln_{\mathrm g}/c と窓の半幅の積が、結合の無い媒質の中心での αabsL/2\alpha_{\mathrm{abs}}L/2 に等しいことを示し、数値を出せ。
解答

(1) ng=c/vg=2.998×108/17=1.76×107n_{\mathrm g}=c/v_{\mathrm g}=2.998\times10^8/17=1.76\times10^7。式 (11.17) から Ωc=ωp/ng=3.99×1010/4199=9.50×106 s−1\Omega_{\mathrm c}=\omega_{\mathrm p}/\sqrt{n_{\mathrm g}}=3.99\times10^{10}/4199=9.50\times10^6\,\mathrm{s^{-1}}(2π×1.51 MHz2\pi\times1.51\,\mathrm{MHz})。結合は線幅 γ1\gamma_1 より弱い。

(2) Ωc2/(2γ1)=(9.50×106)2/(2×6.28×107)=7.18×105 s−1\Omega_{\mathrm c}^2/(2\gamma_1)=(9.50\times10^6)^2/(2\times6.28\times10^7)=7.18\times10^5\,\mathrm{s^{-1}}、周波数では 2π2\pi で割って 114 kHz114\,\mathrm{kHz}。線幅 10 MHz10\,\mathrm{MHz} の吸収の線の中に、その約 1/901/90 の幅の透明な窓が開く。パルスの帯域がこの窓に収まるには、パルスの長さは数 μs\mu\mathrm s 以上でなければならない。

(3) τ⋅Ωc2/(2γ1)=Lcωp2Ωc2⋅Ωc22γ1=Lωp22cγ1\tau\cdot\Omega_{\mathrm c}^2/(2\gamma_1)=\dfrac{L}{c}\dfrac{\omega_{\mathrm p}^2}{\Omega_{\mathrm c}^2}\cdot\dfrac{\Omega_{\mathrm c}^2}{2\gamma_1}=\dfrac{L\omega_{\mathrm p}^2}{2c\gamma_1} で、Ωc\Omega_{\mathrm c} が消える。例題 11.3 (2) の αabs=ωp2/(cγ1)\alpha_{\mathrm{abs}}=\omega_{\mathrm p}^2/(c\gamma_1) と比べて、これは αabsL/2\alpha_{\mathrm{abs}}L/2 に等しい。数値は τ=200×10−6×1.76×107/(2.998×108)=11.8 μs\tau=200\times10^{-6}\times1.76\times10^7/(2.998\times10^8)=11.8\,\mu\mathrm s、積は 8.458.45。結合を弱めて光を遅くすると窓が同じ割合で狭くなり、遅延と帯域の積は媒質の光学的な厚さだけで決まる。同じ帯域のパルスをより長く遅らせるには、原子の数を増やして光学的な厚さを大きくする必要がある。

ほかの巻では

木下の勉強ステータス

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

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

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