tecktom.com

出てくるもの(1):第二高調波・位相整合・パラメトリック増幅

2026年10月15日15分で読めますAI生成教科書『虎の巻』第17章 非線形光学
目次
  1. 基本波の減らない第二高調波
  2. 複屈折による位相整合
  3. 擬似位相整合
  4. 基本波が減るとき
  5. パラメトリック増幅
  6. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第17章 §17.3 ★★★ 第二高調波の強度は L2 sinc2(ΔkL/2)L^2\,\mathrm{sinc}^2(\Delta kL/2) に比例する。正常分散が作る不整合を、複屈折か χ(2)\chi^{(2)} の周期的な反転で消す。高い周波数の光を当てると、低い周波数の2つの波は指数関数的に増える

基本波の減らない第二高調波

第二高調波への変換が小さい間は、基本波の振幅 AωA_\omega を一定とおいてよい。式 (17.8) の第1式を 00 から LL まで積分すると

A2ω(L)=iωχ(2)2n2ωcAω2 eiΔkL−1iΔkA_{2\omega}(L)=\frac{\mathrm{i}\omega\chi^{(2)}}{2n_{2\omega}c}A_\omega^2\,\frac{\mathrm{e}^{\mathrm{i}\Delta kL}-1}{\mathrm{i}\Delta k}

で、∣eiΔkL−1∣=2∣sin⁡(ΔkL/2)∣\lvert \mathrm{e}^{\mathrm{i}\Delta kL}-1\rvert =2\lvert \sin(\Delta kL/2)\rvert を使い、強度に直すと

I2ω=ω2(χ(2))2L22nω2n2ωc3ε0 Iω2 sinc2(ΔkL2),sinc u=sin⁡uu(17.10)I_{2\omega}=\frac{\omega^2(\chi^{(2)})^2L^2}{2n_\omega^2n_{2\omega}c^3\varepsilon_0}\,I_\omega^2\,\mathrm{sinc}^2\Bigl(\frac{\Delta kL}{2}\Bigr),\qquad\mathrm{sinc}\,u=\frac{\sin u}{u}\tag{17.10}

を得る。変換効率 η=I2ω/Iω\eta=I_{2\omega}/I_\omega は基本波の強度 IωI_\omega に比例するので、パルスで強度を上げるほど効率が上がる。位相が整合していれば(Δk=0\Delta k=0)I2ωI_{2\omega} は結晶の長さの2乗に比例する。これは、各点の分極が作る波が同じ位相で足し合わさり、振幅が LL に比例して増えるからである。

Δk≠0\Delta k\ne0 では、I2ωI_{2\omega} は ∣Δk∣L=2π\lvert \Delta k\rvert L=2\pi ごとに 00 に戻る。最初の極大までの距離

Lc=π∣Δk∣=λ4∣n2ω−nω∣(17.11)L_{\mathrm c}=\frac{\pi}{\lvert \Delta k\rvert }=\frac{\lambda}{4\lvert n_{2\omega}-n_\omega\rvert }\tag{17.11}

をコヒーレンス長と呼ぶ(λ\lambda は基本波の真空の波長。Δk=2kω−k2ω=(2ω/c)(nω−n2ω)\Delta k=2k_\omega-k_{2\omega}=(2\omega/c)(n_\omega-n_{2\omega}) を使った)。添字の c はこの節だけの記号である。長さ LcL_{\mathrm c} を越えて結晶を長くしても、例題 17.3 の逆変換が始まって第二高調波は減る(図 17.1 の点線)。

第11章で見たとおり、透明な結晶では共鳴が紫外にあり、可視から近赤外は正常分散の範囲にある。したがって n2ω>nωn_{2\omega}>n_\omega が普通に成り立ち、Δk≠0\Delta k\ne0 である。屈折率の差が 0.080.08 なら、波長 1 μm1\,\mu\mathrm m の基本波で LcL_{\mathrm c} は約 3 μm3\,\mu\mathrm m にすぎない。P. A. Franken らの最初の第二高調波の実験(17.5節)は、位相を整合させる方法が知られる前のもので、第二高調波はごく弱かった。

複屈折による位相整合

第15章(第15章)で扱う一軸結晶では、光学軸に垂直に偏った常光線の屈折率 non_{\mathrm o} は進む向きに依らず、光学軸を含む面内に偏った異常光線の屈折率は、光学軸と進む向きのなす角 θ\theta によって

1ne(θ)2=cos⁡2θno2+sin⁡2θne2\frac{1}{n_{\mathrm e}(\theta)^2}=\frac{\cos^2\theta}{n_{\mathrm o}^2}+\frac{\sin^2\theta}{n_{\mathrm e}^2}

と変わる(ここでは第15章の結果を引く)。ne<non_{\mathrm e}<n_{\mathrm o} の結晶(負の一軸結晶)で、基本波を常光線に、第二高調波を異常光線に取ると、条件 ne2ω(θ)=noωn_{\mathrm e}^{2\omega}(\theta)=n_{\mathrm o}^{\omega} を満たす角 θpm\theta_{\mathrm{pm}} が存在しうる。2つの屈折率の差を、分散でなく複屈折で埋めるのである。これを複屈折による位相整合と呼び、2つの基本波の光子が同じ偏光で入る形をタイプ I と呼ぶ。J. A. Giordmaine(1962、Phys. Rev. Lett. 8, 19)と P. D. Maker ら(1962、同誌 8, 21)が独立に示した。

擬似位相整合

もう1つの方法は、χ(2)\chi^{(2)} の符号を距離 LcL_{\mathrm c} ごとに反転させることである。強誘電体の結晶では、強い電場を掛けて分極の向きを領域ごとに反転させ、χ(2)\chi^{(2)} の符号を周期 Λ=2Lc\Lambda=2L_{\mathrm c} で入れ替えることができる。逆変換が始まる所で χ(2)\chi^{(2)} の符号が変わると、例題 17.3 の相対位相 θ\theta が π\pi だけ飛び、第二高調波は再び増える区間に入る。

χ(2)\chi^{(2)} が ±χ(2)\pm\chi^{(2)} を周期 Λ\Lambda で繰り返す矩形波は、Fourier 級数で χ(2)(z)=χ(2)∑m2imπei2πmz/Λ\chi^{(2)}(z)=\chi^{(2)}\sum_{m}\frac{2}{\mathrm{i} m\pi}\mathrm{e}^{\mathrm{i}2\pi mz/\Lambda}(mm は正負の奇数)と書くことができる。Λ=2π/∣Δk∣\Lambda=2\pi/\lvert \Delta k\rvert に選ぶと、m=±1m=\pm1 のうち一方の項の ei2πmz/Λ\mathrm{e}^{\mathrm{i}2\pi mz/\Lambda} が eiΔkz\mathrm{e}^{\mathrm{i}\Delta kz} を打ち消し、それ以外の項は不整合のまま残るので平均して増えない。したがって

Λ=2π∣Δk∣=2Lc,χeff(2)=2πχ(2)(17.12)\Lambda=\frac{2\pi}{\lvert \Delta k\rvert }=2L_{\mathrm c},\qquad\chi^{(2)}_{\mathrm{eff}}=\frac{2}{\pi}\chi^{(2)}\tag{17.12}

で、効率は完全な位相整合の (2/π)2=0.405(2/\pi)^2=0.405 倍になる。これを擬似位相整合と呼ぶ。J. A. Armstrong, N. Bloembergen, J. Ducuing, P. S. Pershan(1962、Phys. Rev. 127, 1918)が結合波方程式の論文の中で提案した。複屈折による位相整合と違い、偏光と進む向きを自由に選べるので、χ(2)\chi^{(2)} のテンソルの最も大きな成分を使える。

図

図 17.1 基本波の減らない第二高調波の成長(式 (17.8) を数値で積分)。位相整合(点線)では z2z^2 で増える。不整合(細かい点線)では LcL_{\mathrm c} ごとに増減を繰り返し、最大でも (2/π)2=0.405(2/\pi)^2=0.405 にとどまる。擬似位相整合(実線)は χ(2)\chi^{(2)} の符号を LcL_{\mathrm c} ごとに反転させたもので、平均して ((2/π)z/Lc)2((2/\pi)z/L_{\mathrm c})^2 で増え、z=4Lcz=4L_{\mathrm c} で 6.486.48 になる。

基本波が減るとき

変換が進むと基本波が減る。Δk=0\Delta k=0、nω=n2ω=nn_\omega=n_{2\omega}=n として式 (17.8) を厳密に解く。入口で Aω=A0A_\omega=A_0(実数)、A2ω=0A_{2\omega}=0 とすると、例題 17.3 のとおり A2ω=iBA_{2\omega}=\mathrm{i} B(BB は実数)の形が保たれ、dB/dz=κAω2\mathrm{d} B/\mathrm{d} z=\kappa A_\omega^2、dAω/dz=−κBAω\mathrm{d} A_\omega/\mathrm{d} z=-\kappa BA_\omega(κ=ωχ(2)/(2nc)\kappa=\omega\chi^{(2)}/(2nc))になる。Manley–Rowe の関係から Aω2+B2=A02A_\omega^2+B^2=A_0^2 なので、解は

Aω(z)=A0 sech(Γz),B(z)=A0tanh⁡(Γz),Γ=κA0=ωχ(2)A02nc(17.13)A_\omega(z)=A_0\,\mathrm{sech}(\Gamma z),\qquad B(z)=A_0\tanh(\Gamma z),\qquad\Gamma=\kappa A_0=\frac{\omega\chi^{(2)}A_0}{2nc}\tag{17.13}

である。変換効率は η=tanh⁡2(ΓL)\eta=\tanh^2(\Gamma L) で、ΓL→∞\Gamma L\to\infty で 11 に近づく。基本波の減らない近似の効率は (ΓL)2(\Gamma L)^2 で、式 (17.10) の Δk=0\Delta k=0 の場合と一致し、tanh⁡2u=u2−23u4+⋯\tanh^2u=u^2-\frac23u^4+\cdots なので、近似の相対誤差は 23η\frac23\eta の程度である。

パラメトリック増幅

式 (17.7) で、高い周波数 ω3\omega_3 の強い波(ポンプ光)を当て、ω1\omega_1 の弱い波(シグナル光)を入れる。差周波の分極から ω2=ω3−ω1\omega_2=\omega_3-\omega_1 の波(アイドラー光)が生まれる。ポンプ光が減らないとして A3A_3 を一定とおき、Δk=0\Delta k=0 とすると、第1式を zz で微分して第2式の複素共役を入れ、

d2A1dz2=iκ1A3dA2∗dz=κ1κ2∣A3∣2A1\frac{\mathrm{d}^2A_1}{\mathrm{d} z^2}=\mathrm{i}\kappa_1A_3\frac{\mathrm{d} A_2^*}{\mathrm{d} z}=\kappa_1\kappa_2\lvert A_3\rvert ^2A_1

を得る。解は指数関数的に増える。入口でアイドラー光が無い(A2(0)=0A_2(0)=0)とき、

A1(z)=A1(0)cosh⁡(gz),A2(z)=iκ2A3gA1(0)∗sinh⁡(gz),g2=κ1κ2∣A3∣2=ω1ω2(χ(2))2I32n1n2n3c3ε0(17.14)A_1(z)=A_1(0)\cosh(gz),\qquad A_2(z)=\mathrm{i}\frac{\kappa_2A_3}{g}A_1(0)^*\sinh(gz),\qquad g^2=\kappa_1\kappa_2\lvert A_3\rvert ^2=\frac{\omega_1\omega_2(\chi^{(2)})^2I_3}{2n_1n_2n_3c^3\varepsilon_0}\tag{17.14}

である。シグナル光の強度は G=cosh⁡2(gL)G=\cosh^2(gL) 倍になり、アイドラー光の光子の流れは Manley–Rowe の関係どおり、増えたシグナル光の光子の流れ sinh⁡2(gL)\sinh^2(gL) 倍と同じだけ生まれる。これをパラメトリック増幅と呼ぶ。位相が整合していないときは gg が g2−(Δk/2)2\sqrt{g^2-(\Delta k/2)^2} に置き換わり、∣Δk∣>2g\lvert \Delta k\rvert >2g では増幅が止まって振動になる。

第二高調波との違いは、高い周波数から低い周波数へエネルギーが流れることと、増え方が z2z^2 でなく指数関数になることである。シグナル光が増えるとアイドラー光の源の分極が大きくなり、増えたアイドラー光がまたシグナル光の源の分極を大きくするという正のフィードバックがあるからである。パラメトリック増幅の媒質を光共振器(第5章、第5章)に入れ、利得が共振器の損失を上回ると発振する。これが光パラメトリック発振器で、ω1\omega_1 と ω2\omega_2 の配分を位相整合の条件(結晶の角度や温度)で連続的に変えられる波長可変の光源になる。シグナル光を入れなくても、真空の揺らぎからシグナル光とアイドラー光が対で生まれる。この自発パラメトリック下方変換は場を量子化して初めて書ける過程で、第33章(第33章)で光子対の源として、第31章(第31章)でスクイーズド状態の源として扱う。

例題 17.4 コヒーレンス長と擬似位相整合の周期 ★★★

オリジナルある結晶の屈折率を、波長 1064 nm1064\,\mathrm{nm} で nω=2.156n_\omega=2.156、532 nm532\,\mathrm{nm} で n2ω=2.234n_{2\omega}=2.234 とする(偏光は同じ)。

  1. (1) 第二高調波のコヒーレンス長 LcL_{\mathrm c} を求めよ。
  2. (2) 擬似位相整合に必要な反転の周期 Λ\Lambda を求めよ。
  3. (3) 長さ 10 mm10\,\mathrm{mm} の結晶を (a) 位相整合なし (b) 擬似位相整合 で使うとき、第二高調波の強度は完全な位相整合の何倍になるか。
解答

(1) 式 (17.11) から Lc=1.064 μm/(4×0.078)=3.41 μmL_{\mathrm c}=1.064\,\mu\mathrm m/(4\times0.078)=3.41\,\mu\mathrm m。

(2) Λ=2Lc=6.82 μm\Lambda=2L_{\mathrm c}=6.82\,\mu\mathrm m。

(3) (a) 式 (17.10) の sinc2(ΔkL/2)\mathrm{sinc}^2(\Delta kL/2) は LL を LcL_{\mathrm c} の何倍に取るかで 00 と (2Lc/(πL))2(2L_{\mathrm c}/(\pi L))^2 の間を振動する。最大でも (2×3.41×10−3/(π×10))2=4.7×10−8(2\times3.41\times10^{-3}/(\pi\times10))^2=4.7\times10^{-8} 倍である。(b) 式 (17.12) から (2/π)2=0.405(2/\pi)^2=0.405 倍。10 mm10\,\mathrm{mm} は LcL_{\mathrm c} の約 29002900 倍なので、反転の周期を μm\mu\mathrm m の精度で 14661466 周期そろえる加工が要る。

例題 17.5 複屈折で位相を整合させる角度 ★★

オリジナル負の一軸結晶の屈折率を noω=1.655n_{\mathrm o}^\omega=1.655、no2ω=1.675n_{\mathrm o}^{2\omega}=1.675、ne2ω=1.555n_{\mathrm e}^{2\omega}=1.555 とする。

  1. (1) タイプ I の位相整合 ne2ω(θ)=noωn_{\mathrm e}^{2\omega}(\theta)=n_{\mathrm o}^\omega を満たす角 θpm\theta_{\mathrm{pm}} を求めよ。
  2. (2) 同じ結晶で、第二高調波も常光線に取ると位相整合できないのはなぜか。
解答

(1) 本文の ne(θ)n_{\mathrm e}(\theta) の式に条件を入れて sin⁡2θ\sin^2\theta について解くと

sin⁡2θpm=(noω)−2−(no2ω)−2(ne2ω)−2−(no2ω)−2=0.36510−0.356430.41356−0.35643=0.1517,\sin^2\theta_{\mathrm{pm}}=\frac{(n_{\mathrm o}^\omega)^{-2}-(n_{\mathrm o}^{2\omega})^{-2}}{(n_{\mathrm e}^{2\omega})^{-2}-(n_{\mathrm o}^{2\omega})^{-2}}=\frac{0.36510-0.35643}{0.41356-0.35643}=0.1517,

θpm=22.9∘\theta_{\mathrm{pm}}=22.9^\circ。

(2) 常光線の屈折率は向きに依らず、正常分散で no2ω>noωn_{\mathrm o}^{2\omega}>n_{\mathrm o}^\omega が成り立つので、どの向きでも Δk≠0\Delta k\ne0 である。異常光線の屈折率は θ=0\theta=0 の no2ω=1.675n_{\mathrm o}^{2\omega}=1.675 から θ=90∘\theta=90^\circ の 1.5551.555 まで動き、その途中で 1.6551.655 を横切るので、整合する角が存在する。複屈折の大きさ no2ω−ne2ω=0.12n_{\mathrm o}^{2\omega}-n_{\mathrm e}^{2\omega}=0.12 が分散 no2ω−noω=0.020n_{\mathrm o}^{2\omega}-n_{\mathrm o}^\omega=0.020 より大きいことが条件である。

例題 17.6 パラメトリック増幅の利得 ★★

オリジナル波長 532 nm532\,\mathrm{nm} のポンプ光(強度 I3=1 GW cm−2I_3=1\,\mathrm{GW\,cm^{-2}})で、波長 1550 nm1550\,\mathrm{nm} のシグナル光を増幅する。結晶は χ(2)=4 pm V−1\chi^{(2)}=4\,\mathrm{pm\,V^{-1}}、3つの波で屈折率 n=1.66n=1.66、長さ L=5 mmL=5\,\mathrm{mm}、位相は整合しているとする。

  1. (1) アイドラー光の波長を求めよ。
  2. (2) 利得の係数 gg と、シグナル光の強度の増幅率 GG を求めよ。
  3. (3) シグナル光の光子1個あたり、出口で何個のアイドラー光の光子があるか。
解答

(1) 1/λ2=1/λ3−1/λ11/\lambda_2=1/\lambda_3-1/\lambda_1 から λ2=810 nm\lambda_2=810\,\mathrm{nm}。

(2) ω1=2πc/λ1=1.215×1015 s−1\omega_1=2\pi c/\lambda_1=1.215\times10^{15}\,\mathrm{s^{-1}}、ω2=2.326×1015 s−1\omega_2=2.326\times10^{15}\,\mathrm{s^{-1}}、I3=1013 W m−2I_3=10^{13}\,\mathrm{W\,m^{-2}} を式 (17.14) に入れて g=455 m−1=0.455 mm−1g=455\,\mathrm{m^{-1}}=0.455\,\mathrm{mm^{-1}}、gL=2.28gL=2.28。G=cosh⁡2(2.28)=24.2G=\cosh^2(2.28)=24.2。

(3) 入口のシグナル光の光子の流れを 11 とすると、出口のシグナル光は cosh⁡2(gL)=24.2\cosh^2(gL)=24.2、アイドラー光は sinh⁡2(gL)=23.2\sinh^2(gL)=23.2 である。増えたシグナル光の光子 23.223.2 個と同じ数のアイドラー光の光子が生まれ、ポンプ光の光子が同じ数だけ消えた(Manley–Rowe の関係)。

木下の勉強ステータス

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

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

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