header image
  • プログラム・ソフト
  • 物理学・工学
  • ガジェット
  • 当サイトについて
  • お問い合わせ
  • プログラム・ソフト
  • 物理学・工学
  • ガジェット
  • 当サイトについて
  • お問い合わせ
  • HOME
  • 物理学・工学
  • 光学
  • 偏光をザックリと数式的に扱ってみる

偏光をザックリと数式的に扱ってみる

投稿日:2020年3月8日 更新日:2022年6月15日

  • Tags: [数学],

nocap

内容

  • 直線偏光、円偏光とは
  • 偏光の状態と操作の数学的記述方法(ジョーンズベクトル、行列)
  • 内容
  • 動機
  • 直線偏光・円偏光とは
    • 直線偏光のイメージ
    • 直線偏光の角度を変える
    • 円偏光を作り出す
  • 数学的な取り扱い(ジョーンズ計算法)
    • 波の数式表現
  • 横成分・縦成分の表現
    • 偏光の操作
  • まとめ

動機

2枚の偏光板を90度回転させて重ねると光が通らない、いわゆるクロスニコルは知っており、直線偏光のイメージは頭の中にありました。しかし円偏光と聞いても、少し前までなんだそれは状態でした。最近なんとなく掴めてきて、偏光状態や偏光の操作を数学的に表せることもこともわかってきましたので、ここで書き出してみようと思います。偏波にも共通することも述べると思いますが、ここでは光に限定して話を進めていきます。

直線偏光・円偏光とは

直線偏光のイメージ

光は電磁波の一種で、波として扱うことができます。波は光の進行方向に対して垂直に振動するわけですが、”垂直に振動”といってもいろいろな振動方向が考えられます(下図の偏光板を通る前)。このいろいろな振動方向のうち、ある一方向の波のみを切り取ったものが、直線偏光となります(下図の偏光板を通った後)。直線偏光は一方向にのみ振動しているため、光が進む正面から振動をみると直線にみえます。

偏光板による直線偏光の切り出し
Fig.1. 偏光板による直線偏光の切り出し

直線偏光の角度を変える

直進する一本の光の中に、いろいろな振動方向の光が含まれるとき、この光を通す偏光板(直線偏光に切り出すフィルタ)の角度を変えれば、角度の異なる直線偏光が得られます。

では既に偏光板で切り出した直線偏光の角度を変えるにはどうしたらいいでしょうか。それには、縦方向の直線偏光と横方向の直線偏光の成分の比率を変えることで実現できます。

ベクトルのように考えて縦成分と横成分の大きさが決まれば、縦方向と横方向ベクトルの和である斜め方向のベクトルが決まり、その斜め角度も決定できます。例えば、水平線から45度傾いた直線偏光であれば、縦成分と横成分の強度比を1:1にすることで表現できます(下図)。このように縦成分と横成分の強度比から直線偏光の角度が決定し、この比率を変えれば角度を変えることができるということになります。

縦・横方向の直線偏光への分解
Fig.2. 縦・横方向の直線偏光への分解

円偏光を作り出す

それでは円偏光の話に入っていきます。円偏光は縦成分の波と横成分の波の位相がπ/2(つまり波の山と山の間隔λ÷4)ずれているときにみられます。つまり横成分と縦成分のそれぞれの波の山(最大値)になるタイミングがずれると、縦成分と横成分の和となるベクトルは回転するようになります。縦と横成分の波の山が揃っていると直線を描き、λ/4ずれていると円を描き、それ以外だと楕円を描きます。円・楕円の場合、縦と横のどちらの山が先に来るかによって、右回転・左回転が決まります。

縦・横方向の波の位相のずれによって生じる円偏光
Fig.3. 縦・横方向の波の位相のずれによって生じる円偏光

数学的な取り扱い(ジョーンズ計算法)

波の数式表現

波といったらまず思い浮かぶのはサイン波(コサイン波)でしょう。 しかし工学的によく利用されるのは$e^{i\theta}$のような複素数です($i$ :虚数)。これは有名なオイラーの公式

$$\cos\theta+i\sin\theta=e^{i\theta} $$

の通り、コサインに $-i\sin\theta$ という虚数部分がくっついたものであると解釈でき、単純なコサインよりも計算が便利になることが多々あるため利用されます。これを複素正弦波と呼びます。z方向に進む光があるとき、波は $e^{i(\omega t-kz)}$( $\omega$ : 角周波数, $k$: 波数)で示せます。横方向・縦方向それぞれの振幅を $E_x$, $E_y$とすると、それぞれ $$\left(\begin{array}{c}E_x\\ E_y\end{array}\right)=\left(\begin{array}{c}A_xe^{i(\omega t-kz)}\\ A_ye^{i(\omega t-kz+\phi)}\end{array}\right)$$

($A_x$, $A_y$: $E_x$, $E_y$ の振幅の大きさ, $\phi$: $E_x$ と $E_y$ の位相差)になります。

横成分・縦成分の表現

横方向の振幅と縦方向の振幅を数学的に示したものとして、ジョーンズベクトルというものがあります。これは縦成分・横成分それぞれの振幅を$x$, $y$として表したものになります。加えて縦・横成分の位相差の情報も虚部に含まれています。縦・横それぞれの成分を

$$\displaystyle \left(\begin{array}{c}E_x\\ E_y\end{array}\right)=\left(\begin{array}{c}A_xe^{i(\omega t-kz)}\\ A_ye^{i(\omega t-kz+\phi)}\end{array}\right) =\left(\begin{array}{c}A_xe^{i(\omega t-kz)}\\ A_ye^{i(\omega t-kz)}e^{i\phi}\end{array}\right)$$

としたとき、共通項をひとくくりにして $$\displaystyle \left(\begin{array}{c}E_x\\ E_y\end{array}\right) =e^{i(\omega t-kz)}\left(\begin{array}{c}A_x\\ A_ye^{i\phi}\end{array}\right)$$

になり、この$\displaystyle \left(\begin{array}{c}A_x\\ A_ye^{i\phi}\end{array}\right)$の部分がジョーンズベクトルになります。$E_x$ と $E_y$ の位相差 $\phi$ が$\pi/2$ の場合、$e^{i\pi/2}=\cos(\pi/2)+i\sin(\pi/2)=0+i=i$となり、ジョーンズベクトルは $$\displaystyle \left(\begin{array}{c}1\\ i\end{array}\right)$$

となります。このような形で位相情報を積で示せるのが複素数の利点の1つだといえます。なおベクトルの大きさを1に規格化するのが一般的であるため $$\displaystyle \left(\begin{array}{c}E_x\\ E_y\end{array}\right) =\frac{1}{\sqrt{A_x^2+A_y^2}}e^{i(kz-\omega t)}\left(\begin{array}{c}A_x\\ A_ye^{i\phi}\end{array}\right)$$

とした$\displaystyle \frac{1}{\sqrt{A_x^2+A_y^2}}\left(\begin{array}{c}A_x\\ A_ye^{i\phi}\end{array}\right)$部分がジョーンズベクトルになります。

偏光の操作

横と縦のそれぞれの波の位相差を変えることによって、直線偏光の角度を変えたり直線偏光から円偏光に変えたりする働きをするのが波長板です。この働きを数学的に記述するものに、ジョーンズ行列(ジョーンズマトリックス)があります。これはある偏光状態を示すジョーンズベクトルに作用することによって、別の偏光状態を意味するジョーンズベクトルに変換する働きをします。例えば、x軸から45度傾いた方向に振動している直線偏光のジョーンズベクトルは

$$\displaystyle \frac{1}{\sqrt{2}}\left(\begin{array}{c}1\\ 1\end{array}\right)$$

であり、これにジョーンズ行列$\displaystyle \left(\begin{array}{cc}1&0\\ 0&i\end{array}\right)$を作用させると $$\displaystyle \left(\begin{array}{cc}1&0\\ 0&i\end{array}\right)\frac{1}{\sqrt{2}}\left(\begin{array}{c}1\\ 1\end{array}\right)=\frac{1}{\sqrt{2}}\left(\begin{array}{c}1\\ i\end{array}\right)$$

となり、円偏光に変換されたことになります。 $\displaystyle \left(\begin{array}{c}1\\ i\end{array}\right)$では、$x$が1に対して$y$が$i$、すなわちx軸方向の波に対してy軸方向の波の位相差が $\pi/2$ であることになり、円偏光を意味しています。このように直線偏光と円偏光とを変換することのできる光学素子は $\lambda/4$波長板といいます。また直線偏光の傾きや、円偏光・楕円変更の回転方向を変えることができる $\lambda/2$ 波長板というのも存在します。

まとめ

縦方向と横方向の波の強度比、位相差を変えることで直線偏光や円偏光を表現することができる 偏光状態は、縦・横方向のそれぞれの波の強度と位相情報を示したジョーンズベクトルによって記述できる 偏光状態は波長板などの光学素子によって変換でき、この作用はジョーンズ行列によって記述できる こうやって数学的に波を扱ってみると、複素数が便利であることが実感できますね。

“偏光をザックリと数式的に扱ってみる” への1件のコメント

  1. 偏光(偏波)をインタラクティブなグラフで描画 | Y Lab Desk より:
    2020年3月11日 11:21 PM

    […] Pythonにはいろいろなグラフ作成ライブラリがありますが、中でもインタラクティブな(グリグリ動かせる)グラフが作れるplotlyというライブラリがあります。興味が湧き、前回の記事で扱ったような偏光の様子がみられれぱ、より理解が深まるのでは、ということでplotlyを使って偏光を表すグラフを作ってみました。 […]

    返信

コメントを残す コメントをキャンセル

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です

CAPTCHA


フーリエ変換のイメージと数式の意味をわかりやすく説明したい(複素数verも)

フーリエ級数展開からフーリエ変換~非周期関数への拡張へのイメージ

【plotly】偏光板、波長板を回転したときの偏波をグリグリしてみる

連続波形から離散波形へのフーリエ変換の拡張とサンプリング条件

【動くグラフ】物理的な波の扱いのパラメータと複素数表示

偏光(偏波)をインタラクティブなグラフで描画

【動くグラフ】群速度・位相速度と群速度分散によるパルス広がり

  • プロフィール

  • Yoshiharu

  • 現在、ひよっこ光学系エンジニアをやっている一般人です。我流でプログラムをいじったりしています。備忘録、知識整理を兼ねてアウトプットをしていく予定です。

    • 【ラズパイ】格安な抵抗式・静電容量式の土壌湿度センサーを実際に使って比較

      2022.11.23

      【ラズパイ】格安な抵抗式・静電容量式の土壌湿度センサーを実際に使って比較


    • ELEGOO Neptune 2Sをレビュー!2万円台3Dプリンターの実力を実際に使って評価!

      2022.11.03

      ELEGOO Neptune 2Sをレビュー!2万円台3Dプリンターの実力を実際に使って評価!


    • 自作ラズパイゼロケースを3Dプリンター用にモデリング【FreeCAD】

      2022.10.13

      自作ラズパイゼロケースを3Dプリンター用にモデリング【FreeCAD】


    • 【初心者向け】FreeCADの基本的な使い方と手順の流れ

      2022.10.11

      【初心者向け】FreeCADの基本的な使い方と手順の流れ


    • 【Python】ラズパイでLチカ:LEDの取り扱い~抵抗値の計算

      2022.09.07

      【Python】ラズパイでLチカ:LEDの取り扱い~抵抗値の計算


  • アーカイブ

    • 2022年11月
    • 2022年10月
    • 2022年9月
    • 2022年8月
    • 2022年5月
    • 2022年4月
    • 2022年3月
    • 2022年1月
    • 2021年11月
    • 2021年9月
    • 2021年8月
    • 2021年7月
    • 2021年2月
    • 2021年1月
    • 2020年12月
    • 2020年10月
    • 2020年9月
    • 2020年8月
    • 2020年7月
    • 2020年3月
    • 2020年2月
    • 2020年1月
  • カテゴリー

    • Vim
    • ガジェット
    • プログラム・ソフト
    • 光学
    • 情報処理
    • 物理学・工学
    • プログラム・ソフト
    • 物理学・工学
    • ガジェット
    • 当サイトについて
    • お問い合わせ
    Y Lab Desk

    Copyright © Y Lab Desk, 2022 All Rights Reserved.