光線行列解析のソースを表示

'''光線行列解析'''（こうせんぎょうれつかいせき、英 : ray transfer matrix analysis）は、[[レーザー]]等の光学系の設計に用いられる[[光線]]追跡手法の一種である。'''ABCD行列解析'''とも呼ばれる。これは、光学系を光線[[行列]]により記述し、光線を表す[[ベクトル空間|ベクトル]]にこの行列を乗算することで、系内を進む光線を追跡する手法である。加速器物理学においても、[[粒子加速器]]中の磁界を進む粒子を追跡するため、同様の解析手法が用いられる。

以下で述べる手法は、[[近軸近似]]を用いている。すなわち、すべての光線は光軸に対して小さな角度で進み、また光軸の近くにあることを仮定している<ref>タンジェンシャル光線に対する厳密な光線追跡 : [http://spie.org/Documents/ETOP/1991/389_1.pdf Modernizing the Teaching of Advanced Geometric Optics]</ref>。

== 光線行列の定義 ==
[[File:RayTransferMatrixDefinitions.svg|thumb|300px|光線追跡(ABCD)行列解析では、光学要素（ここでは厚レンズ）を、入射面での<math>(x_1, \theta_1)</math> から 出射面での <math>(x_2, \theta_2)</math> への変化として扱う。]]

光線追跡は、入射面と出射面と呼ばれる、光軸に垂直な2つの面に基づいて行われる。一般性を失うことなく、光軸が固定座標系の ''z'' 軸に一致しているとすることができる。系に入射する光線が入射面と交わる点において、光線は光軸から ''x''<sub>1</sub> だけ離れており、光線と光軸のなす角は ''θ''<sub>1</sub> であるとする。光線が系内を進み、出射面と交わる点において、光線は光軸から ''x''<sub>2</sub> だけ離れており、光線と光軸のなす角は ''θ''<sub>2</sub> であるとする。また、入射・出射面中の媒質の[[屈折率]]を ''n''<sub>1</sub>,''n''<sub>2</sub> とする。

このとき、これらの量は、以下の式により関係づけられる。

:<math> \begin{bmatrix} x_2 \\ \theta_2 \end{bmatrix} = \begin{bmatrix} A & B \\ C & D \end{bmatrix} \ \begin{bmatrix} x_1 \\ \theta_1 \end{bmatrix} </math>

ここで、

:<math>A = {x_2 \over x_1 } \bigg|_{\theta_1 = 0} \qquad B = {x_2 \over \theta_1 } \bigg|_{x_1 = 0}</math>

また、

:<math>C = {\theta_2 \over x_1 } \bigg|_{\theta_1 = 0} \qquad D = {\theta_2 \over \theta_1 } \bigg|_{x_1 = 0}</math>

これは、入射面と出射面における光線ベクトルが、2つの面の間の光学系を表す光線行列 '''M''' によって関係づけられることを表している。[[黒体輻射]]に基づく[[熱力学]]的な議論より、光線行列 '''M''' の[[行列式]]は、屈折率の比により表されることが示される。

:<math>\det(\mathbf{M}) = A D - B C = {   n_1  \over  n_2 }  </math>

したがって、入射面と出射面が同一の媒質中にあるか、あるいは異なる媒質であるが同一の屈折率を有している場合には、'''M''' の行列式は1となる。

文献<ref>{{cite book
 |title= Introduction to matrix methods in optics
 |last1= Gerrard |first1= Anthony
 |last2= Burch |first2 = James M.
 |year=1994
 |publisher=Courier Dover
 |location=
 |isbn= 9780486680446
 |pages=
 |url=https://books.google.com/books?id=naUSNojPwOgC
}} </ref>によっては、光線行列について異なった定義を用いている。そこでは光学方向余弦 ''n'' sin''θ'' を ''θ'' の代わりに用いている。これにより、特に屈折に関わる光線行列に変更が生じる。

同様の手法は電気回路の解析にも用いられる。[[二端子対回路]]を参照のこと。

== 例 ==
* 例えば、2つの面の間が自由空間であれば、光線行列は以下で与えられる。

:<math> \mathbf{S} = \begin{bmatrix} 1 & d \\ 0 & 1 \end{bmatrix} </math>

''d'' は光軸に沿って測った面間の距離である。光線追跡の方程式は以下のようになる。

:<math> \begin{bmatrix} x_2 \\ \theta_2 \end{bmatrix} = \mathbf{S} \begin{bmatrix} x_1 \\ \theta_1 \end{bmatrix} </math>

よって、入射光線と出射光線の間の関係は、

:<math> \begin{bmatrix} x_2 & = & x_1 + d\theta_1 \\
\theta_2 & = & \theta_1 \end{bmatrix} </math>

* 他の簡単な例として、[[薄レンズ]]がある。光線行列は以下で与えられる。

:<math> \mathbf{L} = \begin{bmatrix} 1 & 0 \\ \frac{-1}{f} & 1 \end{bmatrix} </math>

''f'' はレンズの[[焦点距離]]である。光学要素の組み合わせを記述するには、各要素を表す光線行列を掛け合わせて、複合光学系を表す光線行列を得れば良い。例えば、長さ ''d'' の自由空間の後に焦点距離 ''f'' のレンズがある場合には、

:<math>\mathbf{L} \mathbf{S} = \begin{bmatrix} 1 & 0 \\ \frac{-1}{f} & 1\end{bmatrix}
\begin{bmatrix} 1 & d \\ 0 & 1 \end{bmatrix}
= \begin{bmatrix} 1 & d \\ \frac{-1}{f} & 1-\frac{d}{f} \end{bmatrix} </math>

行列の積は交換可能でないので、この行列はレンズの後に自由空間がある場合の光線行列と一致しない。

:<math> \mathbf{S L} =
\begin{bmatrix} 1 & d \\ 0 & 1 \end{bmatrix}
\begin{bmatrix} 1 & 0 \\ \frac{-1}{f} & 1\end{bmatrix}
= \begin{bmatrix} 1-\frac{d}{f} & d \\ \frac{-1}{f} & 1 \end{bmatrix} </math>

== 光線行列の表 ==

{| border="1" cellspacing="0" cellpadding="2"
|- style="background-color: #AAFFCC"
! 光学要素
! 光線行列
! 備考
|- 
| 自由空間／一定の屈折率の媒質中の伝搬
| align="center" |<math>\begin{bmatrix} 1 & d\\ 0 & 1 \end{bmatrix} </math>
| ''d'' : 伝搬距離
|- 
| 平面境界での屈折
| align="center" | <math>\begin{bmatrix} 1 & 0 \\ 0 & \frac{n_1}{n_2} \end{bmatrix} </math>
| ''n''<sub>1</sub> : 入射面の屈折率<br/>
''n''<sub>2</sub> : 出射面の屈折率
|- 
| 曲面境界での屈折
| align="center" | <math>\begin{bmatrix} 1 & 0 \\ \frac{n_1-n_2}{R \cdot n_2} & \frac{n_1}{n_2} \end{bmatrix} </math>
| ''R'' = 曲率半径, ''R'' > 0 は凸面(曲率中心が境界の後ろ)を表す<br/>
''n''<sub>1</sub> : 入射面の屈折率<br/>
''n''<sub>2</sub> : 出射面
|- 
| 平面鏡での反射
| align="center" | <math> \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix} </math>
| 光軸に垂直な鏡でのみ成立
|- 
| 曲面鏡での反射
| align="center" | <math> \begin{bmatrix} 1 & 0 \\ -\frac{2}{R_e} & 1 \end{bmatrix} </math>
|''R''<sub>e</sub> = ''R'' cos''θ''  : タンジェンシャル面(水平方向)における実効曲率半径<br/>
''R''<sub>e</sub> = ''R''/cos''θ'' : サジタル面(垂直方向)における実効曲率半径<br/>
''R'' : [[曲率半径]], R > 0 は凸面<br/>
''θ'' : 水平面に対する鏡の入射角
|- 
| 薄レンズ
| align="center" | <math> \begin{bmatrix} 1 & 0 \\ -\frac{1}{f} & 1 \end{bmatrix} </math>
| ''f'' : 焦点距離,  ''f'' > 0 は凸／正(収束)レンズを表す<br/>
焦点距離がレンズの厚みよりも十分大きい場合のみ成立
|-
| 厚レンズ
| align="center" | <math>\begin{bmatrix} 1 & 0 \\ \frac{n_2-n_1}{R_2n_1} & \frac{n_2}{n_1} \end{bmatrix} \ \begin{bmatrix} 1 & t \\ 0 & 1 \end{bmatrix} \ \begin{bmatrix} 1 & 0 \\ \frac{n_1-n_2}{R_1n_2} & \frac{n_1}{n_2} \end{bmatrix}</math>
| ''n''<sub>1</sub> : レンズ外の屈折率,　''n''<sub>2</sub> : レンズ内の屈折率 <br/>
''R''<sub>1</sub> : 第1面の曲率半径,　''R''<sub>2</sub> : 第2面の曲率半径<br/>
''t'' : レンズの厚さ
|- 
| 単一直角プリズム
| align="center" | <math> \begin{bmatrix} k & \frac{d}{nk} \\ 0 & \frac{1}{k} \end{bmatrix} </math>
| ''k'' = cos''ψ''/cos''φ'' : ビーム拡大率<br/>
''φ'' : 入射角, ''ψ'' : 屈折角, ''d'' : プリズム長<br/>
''n'' : プリズム材質の屈折率<ref name=TLO>{{cite book | title = Tunable Laser Optics | author = [[F. J. Duarte]] | publisher = Elsevier-Academic | location = New York |  year = 2003 }} Chapter 6.</ref>

|- 
| ''r'' 個の直角プリズムからなるビームエキスパンダ
| align="center" | <math> \begin{bmatrix} M & B \\ 0 & \frac{1}{M} \end{bmatrix} </math>
| ''M'' はビーム拡大倍率を表し、 <math>M = k_1 k_2 k_3...k_r</math> で与えられる
''B'' は光の総伝搬距離を表す<ref name=TLO />
|}

== 共振器の安定性 ==
光線行列解析は、レーザーで使われるような[[光共振器]]内での光の状態を表現するのに有効である。最も簡単な例として、[[反射率]] 100%、曲率半径 ''R'' の2枚の鏡が、距離 ''d'' だけ離れて対向している共振器を考える。これは、焦点距離 ''f'' = ''R''/2 の同一の薄レンズが距離 ''d'' だけ離れて一列に並んでいるのと等価である。これは共振器に等価なレンズ導波路として知られている。導波路の各繰り返し要素の光線行列は、上で述べたとおり、

:<math>\mathbf{M} =\mathbf{L}\ \mathbf{S} = \begin{bmatrix} 1 & d \\ \frac{-1}{f} & 1-\frac{d}{f} \end{bmatrix} </math>

光線行列解析を用いればは、導波路(と、それに等価な共振器)の安定性を判定することができる。つまり、光が導波路内で周期的に集光され、導波路に沿って進むためには、いかなる条件が必要かを調べることができる。このためには、系の"固有光線"を求めれば良い。すなわち、導波路の繰り返し要素について、入射光線ベクトルに係数 ''λ'' をかけたものが、出射光線ベクトルに一致すれば良い。よって、

:<math> \mathbf{M} \begin{bmatrix} x_1 \\ \theta_1 \end{bmatrix} = \begin{bmatrix} x_2 \\ \theta_2 \end{bmatrix} = \lambda \begin{bmatrix} x_1 \\ \theta_1 \end{bmatrix} </math>

これは[[固有値]]方程式を与える。

:<math> \left[ \mathbf{M} - \lambda\mathbf{I} \right] \begin{bmatrix} x_1 \\ \theta_1 \end{bmatrix} = 0 </math>

'''I''' は2×2の[[単位行列]]である。

光線行列の固有値を求める。式

:<math>\operatorname{det}  \left[  \mathbf{M} - \lambda\mathbf{I} \right] = 0  </math>

は、[[固有多項式|固有方程式]]

:<math> \lambda^2 - \operatorname{tr}(\mathbf{M}) \lambda + \operatorname{det}( \mathbf{M}) = 0 </math>

に書き換えられる。ここで、

:<math>  \operatorname{tr} ( \mathbf{M} )   =    A + D   =   2 - { d \over f }  </math>

は光線行列の[[跡 (線型代数学)|トレース]]であり、また

:<math>\operatorname{det}(\mathbf{M}) = AD - BC  = 1 </math>

は光線行列の行列式である。ここで、安定性パラメータ

:<math> g \ \stackrel{\mathrm{def}}{=}\   { \operatorname{tr}(\mathbf{M}) \over 2 } = 1 - { d \over 2 f } </math>

を導入すると、固有方程式は、

:<math> \lambda^2 - 2g \lambda + 1 = 0  </math>

[[二次方程式の解の公式]]より、固有値は、

:<math> \lambda_{\pm}  =  g \pm \sqrt{g^2 - 1} \, </math>

ここで、''N'' 回繰り返し要素を通過した後の光線を考える。

:<math> \begin{bmatrix} x_N \\ \theta_N \end{bmatrix} = \lambda^N \begin{bmatrix} x_0 \\ \theta_0 \end{bmatrix} </math>

導波路が安定であるならば、いかなる光線も光軸より遠方へと離れては行かないので、 <math>\lambda^N</math> は発散しない。ここで、 <math> g^2 > 1</math> としよう。このとき、2つの固有値はいずれも実数である。 <math> \lambda_+ \lambda_- = 1</math> であるので、一方の固有値は絶対値で1よりも大きくならなくてはならない。これは、この固有ベクトルに対応する光線は収束しないことを表している。ゆえに、安定な導波路においては、<math> g^2 \le 1</math>であり、また固有値は複素数にて表される。

:<math> \lambda_{\pm} = g \pm i \sqrt{1 - g^2} = \cos(\phi) \pm i \sin(\phi) =  e^{\pm i \phi} </math>

ただし<math> g = \cos{\phi}</math>と置換した。

<math> g^2 < 1</math>において、固有値 <math>\lambda_+</math>, <math>\lambda_-</math> に属する固有ベクトルを <math>r_+</math>, <math>r_-</math> とする。 <math>\lambda_+</math> ≠ <math>\lambda_- </math> より固有ベクトルは直交しているので、全ベクトル空間を張る。よって、入射光の光線ベクトルは、 <math>c_+</math>, <math>c_-</math> を用いて

:<math> c_+ r_+ + c_- r_-  </math>

と表される。 ''N'' 回繰り返し要素を通過した後の出力光は、

:<math> \mathbf{M}^N (c_+ r_+ + c_- r_-) = \lambda_+^N c_+ r_+ + \lambda_-^N c_- r_- = e^{i N \phi} c_+ r_+ + e^{- i N \phi} c_- r_- </math>

これは周期関数を表している。

== 光線行列によるガウシアンビームの追跡 ==
光線行列による定式化は、[[ガウシアンビーム]]を表現するのにも便利である<ref>{{cite journal|last1=Rashidian vaziri|first1=M R|title=New ducting model for analyzing the Gaussian beam propagation in nonlinear Kerr media and its application to spatial self-phase modulations|journal=Journal of Optics|volume=15|issue=3|doi=10.1088/2040-8978/15/3/035202|url=https://iopscience.iop.org/article/10.1088/2040-8978/15/3/035202|bibcode=2013JOpt...15c5202R}}</ref>。波長 ''λ''<sub>0</sub> 、曲率半径 ''R'' (発散方向を正、収束方向を負とする)、ビーム径を ''w''、屈折率を ''n'' とすると、複素ビームパラメータ ''q''を
:<math> \frac{1}{q} = \frac{1}{R} - \frac{i\lambda_0}{\pi n w^2} </math>

により定義することができる。ただし ''R'', ''w'', ''q'' は位置の関数である。ビームの軸を ''z'' 方向とし、ビームウエストの位置を ''z''<sub>0</sub> 、[[レイリー長]]を ''z''<sub>R</sub> とすると、この式は以下と等価である。

:<math> q = (z-z_0) + i z_R</math>

''q'' で表されるビームの伝搬は、光線行列を用いて以下のように表される。

:<math> \begin{bmatrix} q_2 \\ 1 \end{bmatrix} = k \begin{bmatrix} A & B \\ C & D \end{bmatrix} \ \begin{bmatrix} q_1 \\ 1 \end{bmatrix} </math>

''k'' は、光線ベクトルの第2要素を1に維持するための規格化定数である。式を展開すると、

:<math> q_2 = k(Aq_1 + B) ,\qquad 1 =   k(Cq_1 + D) \,  </math>

第1式を第2式で割り、規格化定数を消去すると、

:<math> q_2 =\frac{Aq_1+B}{Cq_1+D}</math>

あるいは、逆数を用いて書き直すと、

:<math>  { 1 \over q_2 }    = {  C + D/q_1  \over  A + B/q_1 }   </math>

=== 例 : 自由空間 ===
自由空間を距離 ''d'' だけ進むビームを考える。光線行列は、

: <math>\begin{bmatrix}A&B\\C&D\end{bmatrix}=\begin{bmatrix}1&d\\0&1\end{bmatrix}</math>

したがって、

: <math>q_2 =\frac{Aq_1+B}{Cq_1+D} = \frac{q_1+d}{1} = q_1+d</math>

これは、ガウシアンビームが <math> q = (z-z_0) + i z_R</math> と表されることと合致している。ビームの伝搬に伴い、曲率半径とビーム径の両方が変化する。

=== 例 : 薄レンズ ===
焦点距離 ''f'' の薄レンズを通過するビームを考える。光線行列は、

: <math>\begin{bmatrix}A&B\\C&D\end{bmatrix}=\begin{bmatrix}1&0\\-1/f&1\end{bmatrix}</math>

したがって、

: <math>q_2 =\frac{Aq_1+B}{Cq_1+D} = \frac{q_1}{-\frac{q_1}{f}+1} </math>
: <math>\frac{1}{q_2} = \frac{-\frac{q_1}{f}+1}{q_1} =\frac{1}{q_1}-\frac{1}{f}</math>

1/''q'' の実部のみが変化している。すなわち、曲率半径の逆数は 1/''f'' だけ減少しているが、ビーム径は薄レンズの両側で変化しない。

== 高次の行列 ==
3×3、4×4、6×6といった高次の行列も光学解析に用いられている<ref>W. Brouwer, ''Matrix Methods in Optical Design'' (Benjamin, New York, 1964).</ref><ref>[[Anthony E. Siegman|A. E. Siegman]], ''Lasers'' (University Science Books, Mill Valley, 1986).</ref><ref>H. Wollnik, ''Optics of Charged Particles'' (Academic, New York, 1987).</ref>。特に、4×4行列は、[[超短パルス|フェムト秒レーザー]]の[[プリズムコンプレッサー|パルス圧縮]]のためのプリズム列の設計や解析に用いられている<ref name=TLO />。

== 出典 ==
{{reflist}}

== 参考文献 ==
* Bahaa E. A. Saleh and Malvin Carl Teich (1991). Fundamentals of Photonics. New York: John Wiley & Sons. Section 1.4, pp. 26--36.

== 関連項目 ==
* {{仮リンク|転送行列法 (光学)|en|Transfer-matrix method (optics)}}
* {{仮リンク|線形正準変形|en|Linear canonical transformation}}

== 外部リンク ==
* [https://web.archive.org/web/20010507191110/http://physics.tamuk.edu/~suson/html/4323/thick.html#Matrix Thick lenses (Matrix methods)]
* [http://www.photonics.byu.edu/ABCD_Matrix_tut.phtml ABCD Matrices Tutorial]  複雑な系の光線行列の計算例
* [http://www.photonics.byu.edu/ABCD_Calc.phtml ABCD Calculator]  光線行列計算のための対話型の計算サイト
* {{Wayback|url=https://play.google.com/store/apps/details?id=com.dmt195.simple.abcd.optical.designer |title=Simple Optical Designer (Android App) |date=20150212133414}} 光線行列計算を扱うAndroidアプリケーション

{{デフォルトソート:こうせんきようれつかいせき}}
[[Category:幾何光学]]