tecktom.com

出てくるもの(1):能動モード同期と Gauss 型のパルス

2026年10月15日10分で読めますAI生成教科書『虎の巻』第23章 モード同期と光周波数コム
目次
  1. 1周の主方程式
  2. 定常なパルスは調和振動子の基底状態
  3. 周波数の上での見方
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第23章 §23.3 ★★ 1周ごとの変化を主方程式にまとめると、利得の帯域の有限さは ∂t2\partial_t^2、変調は t2t^2 の項になり、定常なパルスは調和振動子の基底状態と同じ Gauss 型になる。幅は変調の周波数と利得の帯域の幾何平均の逆数に比例する

1周の主方程式

共振器の中を回るパルスを、出口の鏡の位置で1周ごとに観測する。TT を1周ごとに trtt_{\mathrm{rt}} ずつ進む遅い時間、tt をパルスとともに動く速い時間とし、包絡線を a(T,t)a(T,t) と書く(∣a∣2\lvert a\rvert ^2 を出力の単位で測る)。1周の間に包絡線が受ける変化が小さいとき、a(T+trt,t)−a(T,t)≈trt ∂Taa(T+t_{\mathrm{rt}},t)-a(T,t)\approx t_{\mathrm{rt}}\,\partial_Ta と置いて、1周に起こることを足し合わせる。

利得媒質は、第19章(第19章)の均一に広がった線を持つとする。線の中心 ωc\omega_{\mathrm c} から ω\omega だけずれた成分の振幅の利得は g/(1+ω2/Ωg2)g/(1+\omega^2/\Omega_{\mathrm g}^2) で、Ωg\Omega_{\mathrm g} は利得の線の半値半幅(均一な広がりでは 1/T21/T_2)である。パルスのスペクトルが線の幅より狭いとき、g(1−ω2/Ωg2)g(1-\omega^2/\Omega_{\mathrm g}^2) と展開できる。包絡線の成分 e−iωt\mathrm{e}^{-\mathrm{i}\omega t} に対して ∂t2\partial_t^2 は −ω2-\omega^2 をかけるので、−ω2-\omega^2 を ∂t2\partial_t^2 に置き換えると、利得は g(1+Ωg−2∂t2)g(1+\Omega_{\mathrm g}^{-2}\partial_t^2) という作用になる。損失 ll(鏡の透過と散乱)、変調器の振幅の損失 M(1−cos⁡ωmt)M(1-\cos\omega_{\mathrm m}t) を合わせ、どれも1周あたり 11 よりずっと小さいとすると、

trt ∂a∂T=[g−l+gΩg2∂2∂t2−M(1−cos⁡ωmt)]a(23.9)t_{\mathrm{rt}}\,\frac{\partial a}{\partial T}=\Bigl[g-l+\frac{g}{\Omega_{\mathrm g}^2}\frac{\partial ^2}{\partial t^2}-M(1-\cos\omega_{\mathrm m}t)\Bigr]a\tag{23.9}

を得る。これを能動モード同期の主方程式と呼ぶ(Haus (2000), IEEE J. Sel. Top. Quantum Electron. 6, 1173 の形)。変調の角周波数は ωm=ωrep\omega_{\mathrm m}=\omega_{\mathrm{rep}} に合わせる。そうしないと、パルスは1周ごとに変調の違う位相に来て、同じ場所に留まらない。

定常なパルスは調和振動子の基底状態

パルスが変調の損失の最も小さい時刻 t=0t=0 の近くにあり、幅が変調の周期よりずっと短いとすると、M(1−cos⁡ωmt)≈Mωm2t2/2M(1-\cos\omega_{\mathrm m}t)\approx M\omega_{\mathrm m}^2t^2/2 と展開できる。定常な解を a(T,t)=a(t) eλT/trta(T,t)=a(t)\,\mathrm{e}^{\lambda T/t_{\mathrm{rt}}} と置くと、

gΩg2 a′′−Mωm22 t2a=−(g−l−λ) a(23.10)\frac{g}{\Omega_{\mathrm g}^2}\,a''-\frac{M\omega_{\mathrm m}^2}{2}\,t^2a=-(g-l-\lambda)\,a\tag{23.10}

になる。これは調和振動子の Schrödinger 方程式と同じ形で、第4章(第4章)で高次のガウスビームを出したときと同じく、解は Hermite–Gauss 関数である。基底の解 a=exp⁡(−t2/(2τa2))a=\exp(-t^2/(2\tau_a^2)) を入れると a′′=(t2/τa4−1/τa2)aa''=(t^2/\tau_a^4-1/\tau_a^2)a なので、t2t^2 の係数と定数の項がそれぞれつり合って

τa4=2gMωm2Ωg2,g−l−λ=gΩg2τa2=gM2 ωmΩg(23.11)\tau_a^4=\frac{2g}{M\omega_{\mathrm m}^2\Omega_{\mathrm g}^2},\qquad g-l-\lambda=\frac{g}{\Omega_{\mathrm g}^2\tau_a^2}=\sqrt{\frac{gM}{2}}\,\frac{\omega_{\mathrm m}}{\Omega_{\mathrm g}}\tag{23.11}

を得る。強度の半値全幅は τp=2ln⁡2 τa=1.665 τa\tau_{\mathrm p}=2\sqrt{\ln2}\,\tau_a=1.665\,\tau_a である。nn 番目の Hermite–Gauss の解は g−l−λn=(2n+1)g/(Ωg2τa2)g-l-\lambda_n=(2n+1)g/(\Omega_{\mathrm g}^2\tau_a^2) を満たすので、λn\lambda_n は基底の解が最も大きい。利得が飽和して λ0=0\lambda_0=0 まで下がると、ほかの解は1周ごとに減衰し、基底の Gauss 型のパルスだけが残る。Gauss 型の定常なパルスと式 (23.11) の幅の形は、Kuizenga と Siegman が利得と変調の自己無撞着な条件から求めた(Kuizenga, Siegman (1970), IEEE J. Quantum Electron. QE-6, 694)。

式 (23.10) の2つの項の役割は次のとおりである。変調の損失はパルスの中心から外れた時刻ほど大きいので、パルスは1周ごとに時間の上で両端が削られて短くなる。利得は線の中心から外れた周波数の成分ほど小さいので、スペクトルは1周ごとに両端が削られ、パルスは時間の上で広がる。両者がつり合う幅が τa\tau_a である。τa∝(ωmΩg)−1/2\tau_a\propto(\omega_{\mathrm m}\Omega_{\mathrm g})^{-1/2} なので、パルスの幅は利得の帯域の逆数 1/Ωg1/\Omega_{\mathrm g} と変調の周期の幾何平均の程度にとどまり、帯域の逆数までは縮まない。変調の深さ MM への依存は M−1/4M^{-1/4} と弱く、幅を半分にするには MM を16倍にする必要がある。

周波数の上での見方

同じ方程式を縦モードの振幅で書くこともできる。a(t)=∑ℓAℓe−iℓωmta(t)=\sum_\ell A_\ell\mathrm{e}^{-\mathrm{i}\ell\omega_{\mathrm m}t} を式 (23.9) に入れ、cos⁡ωmt e−iℓωmt\cos\omega_{\mathrm m}t\,\mathrm{e}^{-\mathrm{i}\ell\omega_{\mathrm m}t} が ℓ±1\ell\pm1 の項を作ることを使うと、展開する前の式から

trt dAℓdT=[g(1−ℓ2ωm2Ωg2)−l−M]Aℓ+M2(Aℓ+1+Aℓ−1)(23.12)t_{\mathrm{rt}}\,\frac{\mathrm{d} A_\ell}{\mathrm{d} T}=\Bigl[g\Bigl(1-\frac{\ell^2\omega_{\mathrm m}^2}{\Omega_{\mathrm g}^2}\Bigr)-l-M\Bigr]A_\ell+\frac M2\bigl(A_{\ell+1}+A_{\ell-1}\bigr)\tag{23.12}

が得られる。各モードは自分の利得で増え、変調によってとなりのモードと結合する。これは ℓ\ell を格子点とする1次元の強束縛のモデルに、中心から離れるほど深くなる放物線のポテンシャルを足したものと同じ形である。23.1節で述べた「側帯がとなりのモードを駆動する」ことを、式にしたものである。AℓA_\ell を ℓ\ell の滑らかな関数とみなして Aℓ+1+Aℓ−1−2Aℓ≈∂ℓ2AA_{\ell+1}+A_{\ell-1}-2A_\ell\approx\partial_\ell^2A と置くと、式 (23.10) を Fourier 変換した形に戻る。

図

図 23.2 1周期の変調の損失(点線=厳密な 1−cos⁡ωmt1-\cos\omega_{\mathrm m}t、実線=放物線の近似)と、例題 23.3 の定常なパルスの強度(点線、ピークを 11 にした)。パルスは損失の底の幅 0.07 rad0.07\,\mathrm{rad} ほどの所にあり、そこでは放物線の近似と厳密な損失が重なる

例題 23.3 能動モード同期のパルスの幅 ★★

オリジナル例題 23.1 のレーザー(νrep=99.9 MHz\nu_{\mathrm{rep}}=99.9\,\mathrm{MHz})の利得の線を、半値全幅 120 GHz120\,\mathrm{GHz} の Lorentz 型(Ωg=2π×60 GHz\Omega_{\mathrm g}=2\pi\times60\,\mathrm{GHz})とする。1周の振幅の利得を g=0.1g=0.1、変調の深さを M=0.2M=0.2、ωm=ωrep\omega_{\mathrm m}=\omega_{\mathrm{rep}} とする。

  1. (1) 定常なパルスの τa\tau_a と強度の半値全幅を求め、位相がそろった平らなスペクトルの幅 7.38 ps7.38\,\mathrm{ps}(例題 23.1)と比べよ。
  2. (2) ωmτa\omega_{\mathrm m}\tau_a を求め、変調の損失を放物線で近似したことの妥当性を述べよ。
  3. (3) 定常な状態で、利得 gg は損失 ll をどれだけ上回っているか。
解答

(1) 式 (23.11) に ωm=2π×99.9 MHz\omega_{\mathrm m}=2\pi\times99.9\,\mathrm{MHz}, Ωg=2π×60 GHz\Omega_{\mathrm g}=2\pi\times60\,\mathrm{GHz} を入れると

τa=(2×0.10.2×(6.28×108 s−1)2×(3.77×1011 s−1)2)1/4=65.0 ps,\tau_a=\Bigl(\frac{2\times0.1}{0.2\times(6.28\times10^8\,\mathrm{s^{-1}})^2\times(3.77\times10^{11}\,\mathrm{s^{-1}})^2}\Bigr)^{1/4}=65.0\,\mathrm{ps},

強度の半値全幅は 1.665×65.0=108 ps1.665\times65.0=108\,\mathrm{ps} である。平らなスペクトルの 7.38 ps7.38\,\mathrm{ps} より15倍ほど長い。定常なパルスのスペクトルの幅は、時間帯域幅積 0.4410.441 から 4.1 GHz4.1\,\mathrm{GHz} で、利得の帯域 120 GHz120\,\mathrm{GHz} のごく中心の部分しか使っていない。その中の縦モードは 4040 本ほどである。

(2) ωmτa=6.28×108×65.0×10−12=0.041\omega_{\mathrm m}\tau_a=6.28\times10^8\times65.0\times10^{-12}=0.041。1−cos⁡x=x2/2−x4/24+⋯1-\cos x=x^2/2-x^4/24+\cdots の次の項は、パルスの幅の範囲では放物線の項の x2/12≈1.4×10−4x^2/12\approx1.4\times10^{-4} 倍にすぎない。厳密な cos⁡\cos のまま式 (23.9) の定常な解を数値で求めると、強度の2乗平均の幅は Gauss 型の値と 10−310^{-3} 以下の違いで一致する(verify_ch23.py)。

(3) 式 (23.11) で λ=0\lambda=0 とすると g−l=gM/2 ωm/Ωg=0.1×(99.9 MHz/60 GHz)=1.67×10−4g-l=\sqrt{gM/2}\,\omega_{\mathrm m}/\Omega_{\mathrm g}=0.1\times(99.9\,\mathrm{MHz}/60\,\mathrm{GHz})=1.67\times10^{-4} である。この小さな余りが、パルスが変調と利得の帯域から受ける1周あたりの損失を補っている。

ほかの巻では

木下の勉強ステータス

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

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

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