HX-XAFS group of Nagoya Univ. SR-Center and Aichi SR

この文書は読取専用です。文書のソースを閲覧することは可能ですが、変更はできません。もし変更したい場合は管理者に連絡してください。
===== 結晶工学スクールテキスト 16版 の松井先生のテキスト冒頭部の読み下し(の試み) =====

===== - はじめに =====

この文章では、応用物理学会結晶工学分科会が毎年開催している結晶工学スクールの 2022年度、2023年度のテキスト(結晶工学スクールテキスト16版)に
掲載された松井純爾先生のテキスト「歪のあるエピタキシャル層に対するXRCシミュレーション 〜高木-Taupin理論〜」の先頭 3ページほどに展開される式を追い、
2波近似の基本式が現れるあたりまでをしっかりと理解することを目標に、筆者なりの式の整理と注釈を行なったものである。

松井先生の表記と違う表記をとると混乱する可能性もあるが、
できるだけ慣れた表記で式を追わないと間違える可能性が高くなることを避けるためあえて一部表記を変える。
そのことも含めて、ここでの記述のルールを最初にまとめる。

  - 時間方向の振動数 $\nu$ の代わりに、角振動数 $\omega$ を使う(表記に出てくる$\pi$の数を減らすため)。
  - 波数(空間方向の振動数)の定義として $k=1/\lambda$ ではなく、$k=2\pi/\lambda$ を使う(これも表記に出てくる $\pi$ の数を減らすため)。
  - 電磁単位系として MKSA を採用する。
    原文は CGS eus もしくはその派生が使われている様に見えるが、単にこれに不慣れなことと、現在参照できる参考資料の多くが MKSA なので不要な混乱を避けるため。
  - 逆格子ベクトルを表す記号として原文は $\hh$ を使っているが、ここでは $\gg$ を使う。これも慣れている方が間違いにくいというだけ。
  - 式番号は原文のママとする。順序を変える必要があっても、並び順に振り直すようなことはしない。
  - ただし、第一章より深追いするつもりはないので 1-○○ の 1- は省略する。
  - 原文で複数行の式に一つの式番号が振られているようなケースで、より細かく各行に式番号が欲しいときには枝番をつける。\\
    例えば、原文で (1-10) だった式を2つに分けてそれぞれ番号を振る場合 10-1, 10-2 とする("1-"は省略)。
  - 原文に無い式を追加するときは、A1, A2,... 等表記する。(式の文字と紛らわしければ再考するかも)
  - ベクトル微分演算の記号としては、$\rm div$, $\rot$, $\rm grad$ 等の表記と $\nabla$ をあえて混在させる。\\
    松井先生の原文は前者を使用しており、意味合いがクリアになるというメリットもあるが、式展開の途中では煩雑になるため。
  - ここで書いた式と原文の式が対応はしているが完全には一致していないとき式番号に"'"を付けた(例えば(2')の様に)。\\
    一致していない理由はその前後に<color blue>青字</color>で書いた(理由というより単なる推測のときもあるが)
  - 不一致の理由が、表記の違いや流儀の違い等どちらでもあり得るというという場合ではなく、(単純な勘違いだったとしても)
    何らかの間違い/勘違いが含まれていると判断した場合には式を<color red>赤字</color>で示した。

ここで使っている MathJax のマクロ定義は[[:mathjax|こちら]]。

===== - 「1-1. 動力学の基本方程式」より =====

==== - 導入・基本 ====

最初の 4つの式はマクセル方程式そのもの。

\[ \nabla \times \EE = -\mu_0 \Del{\HH}{t}                                                            \tag{1} \]
<color red>\[ \nabla \times \HH = \Del{\DD}{t}                                                                   \tag{2'} \]</color>

<color blue>原文では(2)の右辺に$\mu_0$がかかっている。おそらく単純な勘違い。</color>

\[ \nabla \cdot \DD = 0                                                                               \tag{3} \]
\[ \nabla \cdot \HH = 0                                                                               \tag{4} \]

裸の電荷無しと仮定して(3)の右辺は 0、ただし分極は存在すると考えるので $\epsilon \neq \epsilon_0$(式(2)を$\nabla \times \HH = \epsilon \lDel{\EE}{t}$とは書かない)。

$\DD$や、$\HH$ は時間方向には位相はともかくとして、入射光の$\EE$と同じ振動(角周波数$\omega$)をしていると考える。顕に書くと
\[ \DD(t,\rr) = \exp (i\omega t + \phi) \DD(\rr)                                                      \tag{A1} \]
\[ \HH(t,\rr) = \exp (i\omega t + \phi') \HH(\rr)                                                     \tag{A2} \]
従って、時間微分は係数 $i \omega$ をかける操作に置き換えられる。
\[ \Del{\DD}{t} = \Del{\DD(t, \rr)}{t} = i\omega \exp( i\omega t + \phi ) \DD(\rr) = i\omega \DD(t,\rr) = i\omega \DD      \tag{5'-1} \]
\[ \Del{\HH}{t} = \Del{\HH(t, \rr)}{t} = i\omega \exp( i\omega t + \phi' ) \HH(\rr) = i\omega \HH(t,\rr) = i\omega \HH     \tag{5'-2} \]
もちろん、本家の $\EE$ も同じく
\[ \Del{\EE}{t} = i\omega \EE                                                                         \tag{A3} \]

<color blue>原文では右辺に負号がついているが、この負号があると時間変化の因子を $\exp (-i\omega t)$ と考えていることになり、
後に進行波を表記する際 $\exp i(\omega t - \kk\cdot\rr)$ と書くことと整合しないのでココでは負号を外した
((進行波を $\exp i(\kk\cdot\rr - \omega t)$ と書く流儀もあり得る(というか普通?)ので、ここだけ単独で見ると負号があってもおかしくないが、後の記述と統一が取れていないと判断したという意味))。</color>

==== - 分極・感受率 ====

真空でない媒質中に、電場$\EE$があると一般には分極$\PP$が誘起される。
$\EE$によって$\PP$がどれだけ誘起されるかの係数(一般にはテンソル)を感受率 $\chi$ として $\PP = \epsilon_0 \chi \EE$ と書け、
電束$\DD$は $\epsilon_0 \EE$ と分極$\PP$の合計で
\[ \DD = \epsilon_0 \EE + \PP = \epsilon_0( \EE + \chi \EE ) = \epsilon_0(1+\chi) \EE = \epsilon \EE                         \tag{6'}\]
と書ける(従って $\epsilon = \epsilon_0(1+\chi)$)。
<color blue>原文では式(6)に $\epsilon_0$ がかかっていない。これはおそらく CGS eus 系を採用しているからだと思われる。
原文で $\epsilon = 1 + \chi$ となっているのも同じ理由と推測する。</color>

==== - 式変形開始 ====

式(1) と (2) をメインとし、(3)、(4) を補助的に使いながら電磁波の波動方程式を出すための常套手段を踏む。\\
方針としては、一旦 $\DD$ と $\HH$ だけの式にした後、式(1)は両辺$\rot$をとり、式(2)は両辺$\lDel{}{t}$かける。
その結果両式に現れる$\rot \lDel{\HH}{t}$ を消去することで$\DD$だけの($\PP$は残るが)波動方程式にすることを目指す。

式(6)より、
\[ \epsilon_0 \EE = \DD - \PP                                                                         \tag{A4} \]

式(1) 全体に $\epsilon_0$ をかけて 左辺に式(A4) を代入。
\begin{eqnarray}
     \nabla \times \epsilon_0 \EE & = & -\epsilon_0 \mu_0 \Del{\HH}{t}                                              \\
     \nabla \times (\DD - \PP)    & = & -\epsilon_0 \mu_0 \Del{\HH}{t}                                \tag{1変形-1}
\end{eqnarray}
さらに両辺の $\rot$ をとる
\begin{eqnarray}
     \nabla\times\nabla \times (\DD - \PP)     & = & -\epsilon_0 \mu_0 \nabla \times \Del{\HH}{t}                     \\
     \nabla\times\nabla\times\DD - \rot\rot\PP & = & -\epsilon_0 \mu_0 \nabla \times \Del{\HH}{t}     \tag{1変形-2}
\end{eqnarray}
次に式(2)を時間微分する(時間微分と空間微分の順序入れ替えも行う)
\[ \nabla \times \Del{\HH}{t} = \DDel{\DD}{t}                                                         \tag{2変形-1} \]
(2-1)を(1-2)に代入すると
\[ \nabla\times\nabla\times\DD - \rot\rot\PP = -\epsilon_0 \mu_0 \DDel{\DD}{t}                        \tag{A5} \]
ここで、ベクトル解析の一般公式
\[ \nabla\times\nabla\times\AAA = \nabla\nabla\cdot\AAA - \nabla^2\AAA                                \tag{14} \]
を使って
\[ \nabla \times \nabla \times \DD = \nabla \nabla\cdot \DD - \nabla^2 \DD                            \tag{15} \]
とできる。(3) より $\nabla \cdot \DD = 0$ なので、改めて式(15)は、
\[ \nabla \times \nabla \times \DD = - \nabla^2 \DD                                                   \tag{A6} \]
と書ける。次に、$ \Lap \equiv \nabla^2$ と記号を置き換えつつ、A6 を A5 に使って変形すると
<color red>\[ -\Lap \DD - \rot\rot \PP = -\epsilon_0\mu_0 \DDel{\DD}{t} \Rightarrow \Lap \DD = \epsilon_0\mu_0 \DDel{\DD}{t} -\rot\rot \PP         \tag{9'} \]</color>
となる(この式(9')で、「真空なので$\PP = 0$」としたものが、ごく普通の真空中の電磁波の波動方程式)。\\
<color blue>式(9')は原文の(9)と比べて右辺第一項が違う。
式(2)から $\lDel{\DD}{t} = \nabla \times \HH = \rot\HH$ を使って、式(9') の $\DD$ の2回時間微分を一つ戻すと、
\[ \Lap \DD = \epsilon_0\mu_0 \Del{\rot\HH}{t}-\rot\rot \PP \]
と書けるので、おそらくここの $\HH$ と $\DD$ を取り違えられたものと思われる。また、原文で後に出てくる波動方程式は式(9')と一致している。</color>

==== - $\rot\rot\PP$ をどうにかする ====

=== - $\rot\,\rot \PP$ をどうにかする下準備 1 : $\PP$ を $\DD$ で書く ===
$\rot\rot \PP$ をどうにかする下準備として、$\PP$ を $\DD$ で書く。
「分極・感受率」に書いたように $\PP = \epsilon_0 \chi \EE$ で、$\DD = \epsilon \EE$、さらに $\epsilon = \epsilon_0(1+\chi)$ なので、
\[ \PP = \epsilon_0 \chi \EE = \epsilon_0 \chi \frac{\DD}{\epsilon} = \frac{\epsilon_0 \chi}{\epsilon_0(1+\chi)}\DD \simeq \chi \DD 　　　\tag{10} \]
となる。最後の近似は、$\chi \ll 1$ と考えて $1+\chi \simeq 1$ としている。
この様に、$\PP$が$\DD$で書けたので、$\rot\rot \PP$ の代わりに $\rot\rot \chi\DD$ を考える。

=== - $\rot\rot\PP$ をどうにかする下準備 2 : $\DD$ と $\chi$ を級数展開する ===

次に、$\chi$をフーリエ展開する。$\chi$は物質に付随する量で結晶と同じ周期性を持つはずなので、
結晶の逆格子ベクトル$\gg$で展開できると考えれば、その展開係数を $\chi_{\gg}$として、
\[ \chi = \sum_\gg \chig \exp( -i\gg\cdot\rr )                                                    \tag{11'} \]
と書ける。\\
<color blue>式(11')は原文の(11)と比べると$\pi$の有無が違う。
これは逆格子ベクトル$\gg$の定義が $a^* = (b \times c)/a\cdot(b\times c)$の様になっていて、$\pi$が入っていないからだと思われる
(松井先生の波数の定義が$k=1/\lambda$なのに対して、ここでの波数の定義が $k=2\pi/\lambda$なのと同じ違い)。
これは流儀の違いで一貫していれば問題にならない。
ここで負号が何故出てくるのか(あっても良いが普通はない。両者は$\chig$が変わる($\chi'_\gg = \chi_{-\gg}$のように入れ替わる))疑問だが、後半の式展開で納得できる。</color>

続いて $\DD$ は一次波(波数$\kk_0$)と、"存在し得る"二次波(反射波)を全て足したものだと考える。
この時、反射波の波数を$\kk$ と書くことにすると、存在が許されるのは $\kk - \kk_0 = \gg$が満たされる場合だけだというのを思い出し
<color blue>(回折条件なので当然と思ってこの様に解釈したが、動力学の話をしようとしている時に、前提としてこの関係を使って良いのかどうか、恥ずかしながら分からなかった)</color>、
この関係を満たす$\kk$を$\kg$と書くことにする(すなわち $\kg-\kk_0=\gg$ あるいは $\kg = \gg + \kk_0$)。
$\gg_0$を$\gg_0\equiv(0,0,0)$と定義すると、$\kk_0 = \kk_0 + \gg_0 \equiv\kk_{\gg_0}$ も含めて全ての波を $\kg = \gg + \kk_0$ 使って書けるので、

\begin{eqnarray}
   \DD & = & \sum_\gg \Dg \exp\{ i (\omega t - \kg\cdot\rr) \}                                                  \tag{A7}\\
       & = & \exp( i \omega t ) \sum_\gg \Dg \exp( -\kg\cdot\rr)                                                \tag{12の2行目} \\ 
       & = & \exp( i \omega t ) \sum_\gg \Dg \exp\{ -i ( \gg + \kk_0 )\cdot\rr \}                                              \\
       & = & \exp\{ i (\omega t - \kk_0\cdot\rr) \} \sum_\gg \Dg \exp( -i \gg\cdot\rr )                         \tag{12の1行目} \\
\end{eqnarray}
ただし、総和を取る $\gg$の集合に$\gg_0$も含まれる。ここで、$\Dg$は各成分の大きさを表す定数のベクトル。\\
<color blue> 
私の頭だとどうしても式(A7)従って式(12の2行目)が先にでてきてしまうが、式(12の1行目)が先にでてくる松井先生の感覚は
「(基本波も含めて)全ての波は、基本波が $\gg$ で散乱されたもの」というのが芯から納得できている見通しの良さを感じられる。見習いたい。\\
また、後に明らかになるが、この展開は「許される$\kk_g$を与える$\gg$だけを使った展開」と考えるのではなく、「全ての$\gg$を使った展開(普通のフーリエ級数展開)」と考えたほうが筋が良さそう。
その場合も、許されない$\kk_g$に対応する係数$\Dg$を0にするだけ。
</color>

=== - $\rot\,\rot \PP$ をどうにかする下準備 3 : $\chi\DD$ を2つの級数の積で書く ===

ここまでで、$\chi$ と $\DD$ のそれぞれを級数展開した形が得られた。
そこで両者から $\chi\DD$を計算してみる。

\[ \chi = \sum_\gg \chig \exp( -i\gg\cdot\rr )                                                         \tag{11'再掲} \]
\[ \DD = \exp( i\omega t) \sum_\gg \Dg \exp( -i\kg\cdot\rr )                                           \tag{12-2再掲} \]

両式の $\gg$ は独立なので、$\DD$の式の方の $\gg$を$\ggd$と書き直して $\chi\DD$を計算する。
その際、$\chi$の展開はフーリエ展開なので全ての$\gg$に関して総和が取れていれば良く、
事前に任意に指定した逆格子ベクトル$\ggd$分だけずらして足し上げを実行しても問題ないことに留意する。
すなわち、
\[ \chi = \sum_\gg \chigg \exp\{ -i(\gg-\ggd)\cdot\rr \}                                                \tag{11'-2} \]
でも良い。

$\chi$の表式として式(11'-2)を採用すると、
\begin{eqnarray}
\chi\DD & = & [ \sum_\gg \chigg \exp\{ -i(\gg-\ggd)\cdot\rr \} ]
              \{ \exp( i\omega t) \sum_{\ggd} \Dgd \exp( -i\kgd\cdot\rr ) \}                             \tag{A9-1}\\
        & = & \exp(i\omega t)
              \sum_{\ggd}\sum_{\gg} \chigg\Dgd \exp\{ -i(\gg-\gg')\cdot\rr \}\exp( -i\kgd\cdot\rr )      \tag{A9-2}\\
        & = & \exp(i\omega t)
              \sum_{\ggd}\sum_{\gg} \chigg\Dgd \exp\{ -i(\gg-\gg'+\kgd)\cdot\rr \}                       \tag{A9-3}\\
\end{eqnarray}
ここで$\kg = \gg + \kk_0$、$\kgd = \ggd + \kk_0 \Rightarrow \kk_0 = \kgd - \ggd$ だということを思い出すと、
式(A9-3)の指数関数の肩は $\gg - \ggd + \kgd = \gg + \kk_0 = \kg$ と書き換えられる。結果、
\[
   \chi\DD = \exp(i\omega t)
              \sum_{\gg'}\sum_{\gg} \chigg\Dgd \exp(-i\kg\cdot\rr)                                       \tag{A9-4}
\] 
となる。\\
<color blue> $\chi$の定義を$\ggd$だけずらす考え方すごい !! 
また、$\chi$の展開のとき $\exp(-i\kg\cdot\rr)$ の様にマイナスで展開していなければここで指数の肩が$\kg$にまとまることは無かった。そこまで考慮されていることに驚き。</color>\\
原文では $\chi\DD$ を書き下しておらず、式(A9)に相当するものは直接書かれていない。

=== - $\rot\rot\PP$の評価開始 : まずは $\rot\PP=\rot\chi\DD$ ===

$\AAA = \phi\AAA_0$とするとき、ベクトル解析の公式によると
\[ \nabla \times \AAA = \nabla \times ( \phi \AAA_0 ) = \nabla \phi \times \AAA_0 + \phi \nabla\times\AAA_0          \tag{A10-1} \]
である。特に $\AAA_0$が$\rr$の関数でない場合には、第二項が消えて
\[ \nabla \times ( \AAA ) = \nabla \phi \times \AAA_0                                                                \tag{A10-2} \]
となる。
$\phi = \exp( -i \kk\cdot\rr )$の場合には$\nabla \phi = -i \kk \exp( -i \kk \cdot \rr )$ であることを使って、 
\begin{eqnarray}
   \rot \AAA & = & \rot \{ \exp( -i \kk\cdot\rr ) \AAA_0 \} \\
             & = & -i \kk \exp( -i \kk \cdot \rr ) \times \AAA_0 \\
             & = & -i \kk \times \{ \exp(-i \kk \cdot \rr)\AAA_0 \} \\
             & = & -i \kk \times \AAA                                                                                \tag{16}
\end{eqnarray}
となる。

$\exp(i\omega t)$や$\chigg$ は $\rr$ の関数ではないので$\rot$は
$\Dgd\exp(-i\kg\cdot\rr)$の部分にだけかかることを考慮しながら、
式(16)の結果を 式(A9-2) の $\chi\DD$ に使うと、
<color red>
\[
   \rot\chi\DD = \exp(i\omega t) \sum_{\gg,\gg'} \chigg \rot \{ \Dgd \exp( -i\kg\cdot\rr ) \}
               = \exp(i\omega t) \sum_{\gg,\gg'} \chigg (-i\kg) \times \{ \Dgd \exp( -i\kg\cdot\rr ) \}             \tag{A11, 原文17'}
\]
</color>
となる。\\
<color blue>
式(A11)と原文式(17)が対応するはずだが、少し違っている。
全体にかかる因子が$\exp$になっていないのはまず間違いなく単純ミスだが、肝心の $\kg$の外積がない。
さらに原文ではシグマの中の$\exp$の肩の$\kg$が$\hh$になっている。
原文式(18)では、$\kg$になっているので、この点はおそらく、松井先生の筆が滑ったものと思われる。
原文の式(18)は(A11)から導かれるものになるので(A11)が正しいと判断する。
</color>

=== - ついに $\rot\rot\PP = \rot\rot\chi\DD$ ===

式(A11)に再び$\rot$を作用させる。
その結果は先と同様で
\begin{eqnarray}
   \rot\rot\chi\DD & = & \exp(i\omega t)\sum_{\gg,\gg'} \chigg (-i\kg) \times \rot \{ \Dgd \exp( -i\kg\cdot\rr ) \}                     \\
                   & = & \exp(i\omega t)\sum_{\gg,\gg'} \chigg (-i\kg) \times \{ (-i\kg) \times \Dgd \exp( -i\kg\cdot\rr ) \}           \\
                   & = & - \exp(i\omega t)\sum_{\gg,\gg'} \chigg \{\kg \times (\kg \times \Dgd) \} \exp( -i\kg\cdot\rr ) \}          \tag{18} \\
\end{eqnarray}
となる。
ベクトル3重積の公式により式(18)に含まれている3重積は
\[
   \kg \times\ (\kg \times \Dgd) = (\kg\cdot\Dgd)\kg - \kg^2\Dgd                                                                     \tag{19}
\]
となる。
$\kg$方向の単位ベクトルを$\ii_{\kg}$とすると、
右辺第一項は
\[ |\kg|^2|\Dgd|\cos\theta \ii_{\kg} \]
なので、$\Dgd$の$\kg$方向の成分($|\Dgd|\cos\theta \ii_{\kg}$)を$\kg^2$倍したもの($\theta$は$\kg$と$\Dgd$が成す角)。
第二項は
$\Dgd$を$\kg^2$倍したものなので、その差である式(19)は、$\Dgd$の$\kg$に垂直な方向の成分を$\kg^2$倍したもの、但し符号は負になる。やや強引に書くなら
\[
    \kg \times\ (\kg \times \Dgd) = (\kg\cdot\Dgd)\kg - \kg^2\Dgd = \kg^2 \{ (\ii_{\kg}\cdot\Dgd)\ii_{\kg} - \Dgd \} = \kg^2 (\Dgdpara - \Dgd )
\]
従って、
\[
    \kg \times\ (\kg \times \Dgd) = -\kg^2\Dgdp
\]
但し、$\Dgdp$は$\Dgd$の$\kg$に垂直な成分($\Dgdpara$は平行な成分)。原文ではどのベクトルに垂直なのかを示すために$\DD_{\gg'(\gg)}$の様に表記している。\\
<color blue> 原文では三重積の展開形を睨んで$\Dgd$に垂直な成分と看破しているが、なぜそう見えるのはまだわからない(方向がそうだということだけは、元の三重積の形でわかるが、大きさまで直感的に理解するのは難しい)。\\
また、ここで負号がつく(元の$\Dgd$とは逆向きに垂直)のも大事で、ここで負号が出てこないと最後の式がおかしくなる。</color>

これを式(18)に代入すると、
\[
   \rot\rot\chi\DD = \exp(i\omega t)\sum_{\gg,\gg'}\chigg \kg^2\Dgdp \exp( -i\kg\cdot\rr )                                   \tag{A12}
\]

==== - 電磁波の波動方程式に大集合 ====

ここまでの結果を式(9')
\[ \Lap \DD = \epsilon_0\mu_0 \DDel{\DD}{t} - \rot\rot \PP                                                                     \tag{9'再掲} \]
に代入してまとめていく。

まず、式(9')の右辺第一項は、真空の光速は$c = 1/\sqrt{\epsilon_0\mu_0}=\lambda\nu$で
$\DD$は時間微分すると$i\omega$がでる(式(5'-1))ことから
\[ \epsilon_0\mu_0 \DDel{\DD}{t}  =  -\frac{1}{c^2}\omega^2\DD
                                  =  -\frac{4\pi^2\nu^2}{\lambda^2\nu^2}\DD
                                  =  -\frac{4\pi^2}{\lambda^2}\DD
                                  =  -K^2\DD                                                                                      \tag{A13}
\]
と書ける。
但し $K$は真空中での電磁波の波数の絶対値($=2\pi/\lambda$)。
この時点で式(9')は
\[                                         
   \Lap \DD + K^2\DD + \rot\rot \PP = 0                                                                                           \tag{20} 
\]
となる。

次に$\Lap\DD$を計算する。
\[                          
  \DD = \exp( i\omega t) \sum_\gg \Dg \exp( -i\kg\cdot\rr )                                                                    \tag{12-2再掲}
\]  
なので
\begin{eqnarray}
  \Lap\DD & = &   \exp( i\omega t) \sum_\gg \Dg \Lap                                   \exp( -i\kg\cdot\rr ) \\
          & = &   \exp( i\omega t) \sum_\gg \Dg \{ -({k_g}_x^2+{k_g}_y^2+{k_g}_z^2) \} \exp( -i\kg\cdot\rr ) \\
          & = & - \exp( i\omega t) \sum_\gg \Dg \kg^2                                  \exp( -i\kg\cdot\rr )                   \tag{A14}
\end{eqnarray}
となる。

式(20)に式(A14)と式(12-2)および式(A12)を代入すると
\[
          - \exp( i\omega t) \sum_\gg \Dg \kg^2 \exp( -i\kg\cdot\rr )
          + K^2 \exp( i\omega t) \sum_\gg \Dg \exp( -i\kg\cdot\rr )
          + \exp(i\omega t)\sum_{\gg,\gg'}\chigg \kg^2\Dgdp \exp( -i\kg\cdot\rr ) = 0
\]
$\exp(i\omega t)$の因子を落として、全体を$\gg$のシグマにまとめると
\[
          - \sum_\gg \Dg \kg^2 \exp( -i\kg\cdot\rr )
          + K^2 \sum_\gg \Dg \exp( -i\kg\cdot\rr )
          + \sum_{\gg}\sum_{\gg'} \chigg \kg^2 \Dgdp \exp( -i\kg\cdot\rr ) = 0
\]
\[
    \sum_\gg [
          - \Dg \kg^2 \exp( -i\kg\cdot\rr )
          + K^2 \Dg \exp( -i\kg\cdot\rr )
          + \exp( -i\kg\cdot\rr ) \kg^2 \sum_{\gg'} \chigg \Dgdp ] = 0
\]
\[
     \sum_\gg \exp(-i\kg\cdot\rr) [
         - \Dg \kg^2
         + K^2 \Dg
         + \kg^2 \sum_{\gg'} \chigg \Dgdp ] = 0                                                                              \tag{21'}
\]
これが、任意の$\rr$で成立するためには$[]$内が0でなければならない。\\
すなわち、
\[
    -\Dg \kg^2 + K^2 \Dg + \kg^2 \sum_{\gg'} \chigg \Dgdp = 0
\]
並べ直して整理すると
\begin{eqnarray}
    -\Dg \kg^2 + K^2 \Dg & = & - \kg^2 \sum_{\gg'} \chigg \Dgdp  \\
     \frac{\kg^2 - K^2}{\kg^2}\Dg & = & \sum_{\gg'} \chigg \Dgdp                                                   \tag{22'}
\end{eqnarray}
<color blue>
原文式(21)では、右辺は$\Dgdp$ではなく$\Dgd$に戻っている。この点、左辺の$\Dg$は$\kg$に垂直(物質内でも保証できるか ?)ことを考慮して、
右辺も自動的に$\kg$に垂直な成分だけが残る、とか、なにか理屈付けられそうな気がするが今の時点ではわからない。
</color>\\
原文では、ここで右辺の式について「ここで$\sum_{\gg'}$は$\gg'=\gg$を含む全ての$\gg'$について和を取る。」と言及している。
これは、この式の起源が$\chi$のフーリエ級数展開にあって$\gg_0$も含めて全ての空間周波数成分を含むことを改めて確認した文章である。
また、この式がある意味 $\Dg$ のフーリエ級数展開と読めることも説明されていて、非常に趣深い。\\
<color blue>
本文で $\DD$ を級数展開した式(12)を導いた時の、
私の気持ちは、$\Dg$は、全ての $\gg$に関して存在するわけではなく、
$\kk_g = \gg + \kk_0$ を満たす(結果の $\kk_g$ が $|\kk_g| = |\kk_0|$を満たす)ケースだけピックアップして展開した、という気持ちだったが
これを満たさないケースの $\Dg$は$\Dg=0$ だと考えれば、式(12)を単純に $\DD$ のフーリエ級数展開と呼んで良かったかもしれない。
</color>\\

=== - 屈折率 ===

$\chi_0$ はその定義($\chi$のフーリエ展開の0次成分)から、媒質中の平均の感受率。従ってこれで媒質の屈折率が決まる。
屈折率は媒質中での光速($c'$)と真空中での光速($c$)の比だから
\[
    n = \frac{c}{c'} = \frac{\sqrt{1/\epsilon_0\mu_0}}{\sqrt{1/\epsilon\mu}} = \sqrt{\frac{\epsilon\mu}{\epsilon_0\mu_0}}
\]
ここでは $\mu\simeq\mu_0$ と考えており、$\epsilon = \epsilon_0(1+\chi_0)$ なので、
\[
    n = \sqrt{\frac{\epsilon_0(1+\chi_0)\mu_0}{\epsilon_0\mu_0}} = \sqrt{1+\chi_0}
\]
<color blue> 原文はこの逆数になっているが、おそらく間違いでは ? </color>\\
$\chi_0$が十分小さいと考えてこれをテーラー展開すると、
\[
    n = \sqrt{1+\chi_0} = 1 + \frac{1}{2} \chi_0
\]
<color blue>
(そもそも違う式からの導出だが)原文では、$1-|\chi_0|/2$ となっていて、第二項に負号がつくことと $\chi_0$ の絶対値をとっている点が異なる。\\
本文筆者の知識として、物質中の光速(位相速度)はX線の領域では真空中の光速を越えている($n<1$)ことを知っているので、
そのことからすると、ここでの $\chi_0$ は負のはず($n=1+\chi_0/2<1$となるので)。\\
松井先生の式は$\chi_0$の絶対値をとって引いているので($n=1-|\chi_0|/2$)これは確実に 1 より小さい。
1より小さくするために(?)絶対値を取らないといけないことから間接的に、原文でも(当然ではあるが) $\chi_0$は負だと考えていることがわかる。</color>

==== - 2波近似 ====

=== - スカラー化 ===

原文では「波動振幅は大きさのみを扱うからスカラー量として...」と簡単に書かれているが、この「スカラー量」に直す操作も難しい。
この過程で$\Dgdp$が$\Dgd$に戻りそうな気がするが...

どうするのが正解かわからないが、$\Dg$方向の単位ベクトル$\ig$をかける方針でやってみる。

式(22')
\[
     \frac{\kg^2 - K^2}{\kg^2}\Dg  =  \sum_{\ggd} \chigg \Dgdp                                                   \tag{22'再掲}
\]
の両辺に$\ig$をかける
\begin{eqnarray}
     \frac{\kg^2 - K^2}{\kg^2}\Dg\cdot\ig & = & \sum_{\gg'} \chigg \Dgdp\cdot\ig                                                 \\
     \frac{\kg^2 - K^2}{\kg^2}\sDg & = & \sum_{\gg'} \chigg \Dgdp\cdot\ig    
\end{eqnarray}
$\ig$は必ず$\kg$に垂直なので、$\Dgdp$は$\ig$と平行(<color blue>符号がわからない!! 逆向きの可能性は ??</color>)なので
\[
     \frac{\kg^2 - K^2}{\kg^2}\sDg  =  \sum_{\gg'} \chigg \sDgdp                                              
\]
右辺のシグマから $\ggd=\gg$の項を取り出すと
\[
     \frac{\kg^2 - K^2}{\kg^2}\sDg  =  \chi_0\sDg + \sum_{\gg'\neq\gg} \chigg \sDgdp                             \tag{23'}                                           
\]
と書ける。
<color blue>残念ながら$\sDgd$にはならず$\sDgdp$になった。</color>\\

=== - $K$から$k$へ ===

式(24)、(25)を導出するのは断念。途中経過はまるごと注釈に移した((
取り出した項を左辺に移して整理すると、
\begin{eqnarray}
     \frac{\kg^2 - K^2}{\kg^2}\sDg - \chi_0\sDg  & = & \sum_{\gg'\neq\gg} \chigg \sDgdp                                           \\ 
     \frac{\kg^2 - K^2 - \kg^2\chi_0}{\kg^2}\sDg & = & \sum_{\gg'\neq\gg} \chigg \sDgdp                                           \\ 
\end{eqnarray}
<color blue>困った。分母が$\kg^2$なので、分子に追加されるのは$K^2\chi_0$ではなく$\kg^2\chi_0$になってしまった。</color>
\begin{eqnarray}
     \frac{(1-\chi_0)\kg^2 - K^2}{\kg^2}\sDg     & = & \sum_{\gg'\neq\gg} \chigg \sDgdp                                           
\end{eqnarray}
<color blue>どんどん離れていく...。$(1-\chi_0)$なので屈折率とも解釈できない。解釈して $\kg^2 => K_\gg^2$ にしても仕方ないが... わからないので遊んでみる </color>
\begin{eqnarray}
     \frac{(1+\chi_0)\kg^2 -2\chi_0\kg^2 - K^2}{\kg^2}\sDg     & = & \sum_{\gg'\neq\gg} \chigg \sDgdp                             \\
     \frac{K_g^2 -2\chi_0\kg^2 - K^2}{\kg^2}\sDg     & = & \sum_{\gg'\neq\gg} \chigg \sDgdp      
\end{eqnarray}
<color blue>たぶん、$K_g^2 = K^2$なので</color>
\begin{eqnarray}
     \frac{-2\chi_0\kg^2}{\kg^2}\sDg     & = & \sum_{\gg'\neq\gg} \chigg \sDgdp      
\end{eqnarray}
<color blue>なんだこれ?</color>))。

=== - $K$から$k$へは諦めて : 2波近似本番 ===

原文の式(24)、(25)の導出は一旦諦めるが、原文式(26)には(22')から直接いける。

式(22')
\[
     \frac{\kg^2 - K^2}{\kg^2}\Dg = \sum_{\ggd} \chigg \Dgdp                           \tag{22'再掲}
\]
からスタートする。
この式で $\gg=0$ とし シグマからは$\ggd = 0, \gg$ を選ぶと(他は無いものとすると)式(26)が、
また $\gg=\gg$ とし シグマからはおなじく$\ggd = 0, \gg$ を選ぶと(他は無いものとすると)式(27)が得られるのだが、
この操作では、一般の$\gg$と一つ選んだ特定の$\gg$に同じ記号が使われていたりして、場合によっては一瞬混乱の元になる。
ここでは、なるべく紛れなく理解して納得することを目的としているので、一段階余分に途中段階を踏む。

式(22')に残す2つの波として$\gg_1$と$\gg_2$を選ぶものとする。
シグマに現れる $\ggd$にはこの両方を使って(この2つだけを使って)足し上げを行うが、$\gg$ としては、
このどちらか一つを選ぶことになる。結果$\gg=\gg_1$とした式と$\gg=\gg_2$とした式ができる。\\

まず式(22')の右辺を先に展開する。
\[
     \frac{\kk_{\gg}^2 - K^2}{\kk_{\gg}^2}\DD_{\gg}
          = \chi_{\gg-\gg_1} \DD_{\gg_1}^{\perp \gg} + \chi_{\gg-\gg_2} \DD_{\gg_2}^{\perp \gg}      \tag{A15}
\]
これの $\gg$ に $\gg_1$、もしくは $\gg_2$ を代入するとそれぞれ
\begin{eqnarray}
     \frac{\kk_{\gg_1}^2 - K^2}{\kk_{\gg_1}^2}\DD_{\gg_1}
          & = & \chi_{\gg_1-\gg_1} \DD_{\gg_1}^{\perp \gg_1} + \chi_{\gg_1-\gg_2} \DD_{\gg_2}^{\perp \gg_1}      \tag{A16}  \\
     \frac{\kk_{\gg_2}^2 - K^2}{\kk_{\gg_2}^2}\DD_{\gg_2}
          & = & \chi_{\gg_2-\gg_1} \DD_{\gg_1}^{\perp \gg_2} + \chi_{\gg_2-\gg_2} \DD_{\gg_2}^{\perp \gg_2}      \tag{A17}
\end{eqnarray}
が得られる。

次に $\gg_1 = 0$、$\gg_2 = \gg$ とすると、式(A16)、式(A17)はそれぞれ、
\begin{eqnarray}
     \frac{\kk_{0}^2 - K^2}{\kk_{0}^2}\DD_{0}
          & = & \chi_{0-0} \DD_{0}^{\perp 0} + \chi_{0-\gg} \DD_{\gg}^{\perp 0} \\
          & = & \chi_{0} \DD_{0}^{\perp 0} + \chi_{-\gg} \DD_{\gg}^{\perp 0}                           \tag{A16-2} \\

     \frac{\kk_{\gg}^2 - K^2}{\kk_{\gg}^2}\DD_{\gg}
          & = & \chi_{\gg-0} \DD_{0}^{\perp \gg} + \chi_{\gg-\gg} \DD_{\gg}^{\perp \gg} \\
          & = & \chi_{\gg} \DD_{0}^{\perp \gg} + \chi_{0} \DD_{\gg}^{\perp \gg}                        \tag{A17-2}
\end{eqnarray}
となる。

さらに式(A16-2)に現れる $\DD_\gg^{\perp 0}$ の方向は($\kk_0$に垂直な成分なので?)$\DD_0$と同じと考えてよければ
$\DD_{0}^{\perp 0}$は元々$\DD_0$と同じはずなので、
式(A16-2)に出てくる$\DD$は全てスカラ(大きさ)にできる。\\
式(A17-2)でも、$\Dz^{\perp \gg}$ の方向は($\kk_g$に垂直な成分なので?)$\DD_{\gg}$ と同じと考えて良ければ、
$\DD_{\gg}^{\perp \gg}$は元々$\DD_{\gg}$と同じはずなので、
式(A17-2)に出てくる$\DD$は全てスカラ(大きさ)にできる。\\
などと考えると、式(A16-2)、式(A17-2) はそれぞれ、
\begin{eqnarray}
     \frac{\kk_{0}^2 - K^2}{\kk_{0}^2}D_0 & = & \chi_{0} D_0 + \chi_{-\gg} D_{\gg}^{\perp 0}           \tag{A16-3}  \\
     \frac{\kk_{\gg}^2 - K^2}{\kk_{\gg}^2}D_{\gg} & = & \chi_{\gg} D_0^{\perp \gg} + \chi_{0} D_{\gg}  \tag{A17-3}
\end{eqnarray}
となる。

<color blue>これらの式を見てると、$\DD$ と $\DD^{\perp}$ の違いを吸収しているのが偏光因子$C$にも見えるが ? いずれにしても、差異を$C$に押し込むことにすれば</color>
\[
     \frac{\kk_0^2 - K^2}{\kk_0^2} D_0  = \chi_0 D_0 + C \chi_{-\gg} D_\gg          \tag{26a}                                 
\]
\[
     \frac{\kg^2 - K^2}{\kg^2} D_\gg    = C \chi_{\gg} D_0 + \chi_0 D_\gg           \tag{26b}     
\]
<color blue>式番号のダッシュが消える!</color>

===== - 終了 =====

実は、本文を書き始めたのは式(26)を導出するため(納得して正しさを確認するため)であった。
結果、最大限違っていたとしても$D$と$D^\perp$の違いがあるかもしれないが、それを除けば式(26)を導出できた。

一旦はこれで満足するが、今回の読み下しはかなり勉強になったので後日この先も進めてみたい。