平面波

**平面波**（へいめんは、）とは、等位相面が波数ベクトルを法線ベクトルとする等値平面から成る周期関数のことである。

## 平面波の定義
平面波と呼ばれる関数には、「時間変数を持たない平面波」と、「時間変数を持つ平面波」がある。「時間変数を持たない平面波」は、周期関数のフーリエ級数展開や、フーリエ変換、時間発展のないシュレーディンガー方程式の計算に用いられる。「時間変数を持つ平面波」は、波動方程式の解として現れる。

通常、「時間変数を持たない平面波」と、「時間変数を持つ平面波」は、区別されずに混同されて用いられるが、異なるものなので、曖昧さを回避する観点から区別が必要な場合には、用語を使い分けることにする。それぞれの用語の定義は以下に行う。

また、本稿では、「時間変数を持たない平面波」と、「時間変数を持つ平面波」の総称として「平面波」という用語を用いることにする。

### 時間変数を持たない平面波
実数または複素数に値を取る実  変数関数  が時間変数を持たない平面波であるとは、周期  の実1変数の周期関数  と、波数ベクトルと言われる **'  次元実定数ベクトル （但し ）を用いて、

:<math>\Psi(\boldsymbol{x})=f(2\pi \boldsymbol{k} \cdot \boldsymbol{x})</math>

と表されることを意味する。

=== 時間変数を持つ平面波 ===
時間変数を持つ平面波は、[[波動方程式]]の固有解に現れる。

実数または複素数に値を取る関数  が時間変数を持つ平面波であるとは、空間変数  （ 次元実数ベクトル）と時間変数''t'' （実数）と、周期  の実1変数の周期関数  と、波数ベクトル **'（*d* 次元実定数ベクトル、但し ）と、角振動数  を用いて、

:\Phi(\boldsymbol{x},t)=f(2\pi (\boldsymbol{k} \cdot \boldsymbol{x} - \omega t))

であることを意味する。

尚、本稿では、時間変数と空間変数を のように分ける。つまり、変数の最後の成分を時間変数と考える。

### 時間変数を持つ平面波と、時間変数を持たない平面波
物理的には、空間変数  と時間変数  は異なるものであるが、数学ではどちらも単なる変数である。この意味において、 次元の時間変数を持つ平面波は、 変数の時間変数を持たない平面波と見做すことができる。

時間変数を持つ平面波

:\Phi(\boldsymbol{x},t)=f(2\pi (\boldsymbol{k} \cdot \boldsymbol{x} - \omega t))

に対して、新たに  を、空間成分 **' と、時間成分  を並べた  次元の実数ベクトルとする。即ち、

:<math>\boldsymbol{K} = \left( \begin{array}{c}
k_1     \\
\vdots  \\
k_d \\
-\omega
\end{array} \right)</math>

とする。但し、 は、波数ベクトル  の第  成分を意味する。

又、とする。このとき、

:<math>\Phi(\boldsymbol x,t)= f(2\pi (\boldsymbol k\cdot \boldsymbol x-\omega t)) = f(2\pi \boldsymbol K\cdot \boldsymbol X)</math>

のように書くことが出来る。この意味において、 次元の時間変数を持つ平面波は、 変数の時間変数を持たない平面波と見做すことができた。

== 正弦平面波 ==
正弦平面波は、[[正弦波]]の多次元への拡張の1つで、代表的な平面波である。正弦平面波には、実正弦平面波と複素正弦平面波がある。正弦平面波のことを単に平面波ということもあるが、正弦平面波ではない平面波もある。
=== 実正弦平面波の一般式 ===
実正弦平面波は、数学的には振幅 、波数ベクトル 、位相項  の3つの定数/定数ベクトルで特徴付けられる。一般に  次元の実正弦平面波は、時間変数を持たない形で書くと

:<math>A\cos(2\pi(\boldsymbol{K}\cdot\boldsymbol{X}+\delta))</math>

時間変数を持つ形で書くと

:<math>A\cos (2\pi(\boldsymbol{k} \cdot \boldsymbol{x} -\omega  t+ \delta) )</math>

で表される。

ここで、波数ベクトルや時間・空間変数は、それぞれ

:<math>\boldsymbol{k} = \begin{pmatrix}
k_1    \\
\vdots \\
k_d    
\end{pmatrix},
\quad
\boldsymbol{x} = \begin{pmatrix}
x_1    \\
\vdots \\
x_d    
\end{pmatrix},
\qquad
\boldsymbol{K}= \begin{pmatrix}
k_1     \\
\vdots  \\
k_d \\
-\omega
\end{pmatrix},
\quad
\boldsymbol{X}=\begin{pmatrix}
x_1     \\
\vdots  \\
x_d \\
t
\end{pmatrix}</math>

である。

===複素正弦平面波の一般式===
実正弦平面波は重ね合わせの計算などが面倒であることから、計算上のテクニックとして、実正弦平面波の値域を[[オイラーの公式]]を用いて複素数域に拡張した複素正弦波が発案された。古典物理では、複素平面正弦波は実正弦平面波の重ね合わせを計算するための便宜にすぎないが、量子力学では複素平面正弦波を用いなければ説明がつかない現象があるため、計算上の便宜のためだけのものではない。

複素正弦平面波は数学的には、振幅 （複素定数）、波数ベクトル （実定数ベクトル）、位相項 （実定数）の3つの定数/定数ベクトルで特徴付けられる。一般に、  次元の複素正弦平面波は、

:<math>A\exp i(2\pi(\boldsymbol{K} \cdot \boldsymbol{X} +\delta) )</math>

の形で表される。

=== 複素正弦平面波を用いた実正弦平面波の重ね合わせ ===
を実定数（ただし ）としたときに、重ね合わせ

:<math>\sum_{j=1}^m a_j \cos i \theta_j</math>

を計算する問題を考える。

オイラーの公式より、複素数をベクトルのように表記して
:<math>a_j \exp i \theta_j = a_j \begin{pmatrix}
\cos \theta_j \\
\sin \theta_j 
\end{pmatrix} </math>　(2-1)

と見なすことができる。

式(2-1)の右辺に、ベクトルの[[平行四辺形則]]を適用すると

:<math>\begin{cases}
a= \left| \sum\limits_{j=1}^m a_j \exp i \theta_j \right| \\
\theta = \arg  \sum\limits_{j=1}^m a_j \exp i \theta_j  \\
\end{cases}</math>　(2-2)

としたときに、

:<math>a\cos \theta = \sum_{j=1}^m a_j \cos i \theta_j</math>

が成り立つ。従って、重ね合わせ(1)を計算する問題は、式(2-2)の2つの式を求める問題に帰着される。ここで、  は実定数なので、

:<math>( \exp ik \theta )^* = \exp (-ik \theta)</math>

が成り立つ。 は[[複素共役]]を意味する。このことに注意して、 の展開を行うと

:<math>a^2 = \sum_{j=1}^m a_j \exp i \theta_j \sum_{l=1}^m a_l \exp (-i \theta_j )= \sum_{j,l=1}^m a_j a_l \exp i( \theta_j - \theta_l )</math>　(2-3)

が成立する。式(2-3)と、条件''a<sub>j</sub>'' &ge; 0 を考え併せると、式(2-2)は、

:<math>\begin{cases}
a^2 = \sum\limits_{j,l=1}^m a_j a_l \exp i( \theta_j - \theta_l ) \\
\theta = \arg \sum\limits_{j=1}^m a_j \exp i \theta_j 
\end{cases}</math> (2-2’)

と変形できる。従って、重ね合わせを計算する問題は、式(2-2’)を求める問題に帰着される。計算上の便宜としての複素正弦波を持ち出す最大の理由は、式(2-2)から(2-2’)（特に振幅の関係式）が導き出せることにある。

一般には、これ以上簡単な形に変形することは難しいが、いくつかの特殊な場合には振幅の項あるいは位相項の片方あるいは両方がより簡単な形になる。例えば

:<math>\begin{align}
& \theta_1 = \boldsymbol{k} \cdot \boldsymbol{x} - \omega t \\
& \theta_2 = \boldsymbol{k} \cdot \boldsymbol{x} -\omega t+ \delta \\
& m=2 \end{align}</math>

のときには、

:<math>a= \sqrt

以下、<math>2\pi\textbf{E}_{1},\cdots ,2\pi\textbf{E}_{d}</math>
が周期である場合に関して、d次元のフーリエ級数展開について、帰納的に説明していく。ここで、<math>\textbf{E}_{j}</math>は単位行列の第j列ベクトルを意味する。

=== 2次元のフーリエ級数展開 ===
簡単のためdを2とした場合について考えてみる。即ち、2変数関数<math>f(x_1,x_2)</math> が、周期<math>2\pi\textbf{E}_{1},\cdots ,2\pi\textbf{E}_{2}</math>を持つ場合に、1変数のフーリエ級数に帰着することを考える。尚、周期性の定義等、用語の定義を知らずとも、計算の流れのみから本ケースの証明は理解が可能であると思われるため、定義などは後回しにする（一般の場合を考える際に、再定義する）。以下の定理が成り立つ：

\left(\sum_{z_2\in\mathbb{Z}} C_{2,(z_1,z_2)} \exp(2\pi\textbf{i} z_2x_2) \exp(2\pi \textbf{i}z_1x_1) \right)\\
&= \sum_{z_1\in\mathbb{Z}}\left(\sum_{z_2\in\mathbb{Z}} C_{2,(z_1,z_2)} \exp(2\pi\textbf{i} (z_1x_1+z_2x_2))  \right)\\
&= \sum_{z_1\in\mathbb{Z}}\sum_{z_2\in\mathbb{Z}} C_{2,(z_1,z_2)} \exp(2\pi\textbf{i} (z_1x_1+z_2x_2))
\end{align}</math>

を得る。ここで、
:<math>C_{2,(z_1,z_2)} := \int_{{x_2}=0}^{{x_2}=2\pi} C_{1,z_1}(x_2)\exp(2\pi \textbf{i}z_2x_2)\, \mathrm dx_2 </math>
:<math>=\int_{{x_2}=0}^{{x_2}=2\pi} \left(\int_{x_1=0}^{x_1=2\pi} f_{[x_2]}(x_1) \exp(2\pi \textbf{i}{z_1}x_1)\, \mathrm dx_1 \right)\exp(2\pi \textbf{i}z_2x_2)\, \mathrm dx_2 </math>
:<math>=\int_{{x_2}=0}^{{x_2}=2\pi} \left(\int_{x_1=0}^{x_1=2\pi} f(x_1,x_2) \exp(2\pi \textbf{i}{z_1}x_1)\, \mathrm dx_1 \right)\exp(2\pi \textbf{i}z_2x_2)\, \mathrm dx_2 </math>
:<math>=\int_{{x_2}=0}^ C_{2,(z_1,z_2)} \exp(2\pi\textbf{i} (z_1x_1+z_2x_2))
\end{align}</math>

を得る。ここで、
:<math>C_{2,(z_1,z_2)} := \int_{{x_2}=-\pi}^{{x_2}=\pi} C_{1,z_1}(x_2)\exp(2\pi \textbf{i}z_2x_2)\, \mathrm dx_2 </math>
:<math>=\int_{{x_2}=-\pi}^{{x_2}=\pi} \left(\int_{x_1=-\pi}^{x_2=\pi} f_{[x_2]}(x_1) \exp(2\pi \textbf{i}{z_1}x_1)\, \mathrm dx_1 \right)\exp(2\pi \textbf{i}z_2x_2)\, \mathrm dx_2 </math>
:<math>=\int_{{x_2}=-\pi}^{{x_2}=\pi} \left(\int_{x_1=-\pi}^{x_2=\pi} f(x_1,x_2) \exp(2\pi \textbf{i}{z_1}x_1)\, \mathrm dx_1 \right)\exp(2\pi \textbf{i}z_2x_2)\, \mathrm dx_2 </math>
:<math>=\int_{{x_2}=-\pi}^{{x_2}=\pi} \int_{x_1=-\pi}^{x_2=\pi} f(x_1,x_2) \exp(2\pi \textbf{i}{z_1}x_1)\exp(2\pi \textbf{i}z_2x_2)\, \mathrm dx_1  \mathrm dx_2 </math>
:<math>=\int_{{x_2}=-\pi}^^{d}</math> を、 の周期格子という。

簡単な計算から以下の定理が判る。



この定理により、周期行列が存在するような''d'' 重周期関数の問題は、すべて、標準正方格子を周期格子として持つような周期関数の問題に帰着されることが判る。

=== 平面波の周期性 ===
平面波の周期性について、以下の命題が成り立つ。



即ち、命題1は、 の直交補空間の点は皆、波数 **' の平面波  の周期であることを主張している。

以下の定理より、 重周期関数  と同じ  重周期を持つ平面波を沢山作る方法が与えられる。

内積の双線形性より、
:\left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}+ \textbf{T}_{j}\right\rangle
:\left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}\right\rangle + \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|}  \textbf{T}_{j}\right\rangle

で、前節の命題から、
: \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|}  \textbf{T}_{j}\right\rangle = 2\pi

従って、
:\left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}+ \textbf{T}_{j}\right\rangle= \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}\right\rangle + 2\pi

である。従って、
:\Phi_{\textbf{G}({z}_{1},{z}_{2},...,{z}_{d})}(\textbf{x}+ \textbf{T}_{j})=\exp\left(\textbf{i} \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}+ \textbf{T}_{j}\right\rangle \right)
:=\exp\left(\textbf{i} \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}\right\rangle +\textbf{i} \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|}  \textbf{T}_{j}\right\rangle\right)
:=\exp\left(\textbf{i} \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}\right\rangle +2\pi\textbf{i} \right)
:=\exp\left(\textbf{i} \left\langle \textbf{G}({z}_{1},{z}_{2},...,{z}_{d})  \mathrel{|} \textbf{x}\right\rangle \right)

となる。
### 平面波展開

:h(\textbf{x}):=f(\frac{1}{2\pi}T \textbf{x})
と定めると、2\pi\textbf{E}_1,\cdots 2\pi\textbf{E}_dはhの周期である。従って、
:{C}_{d}(z_1,\cdots,z_d):=\int_{{x_d}=0}^{{x_d}=2\pi} \cdots \int_{{x_1}=0}^

## 量子論における平面波
非相対論的な量子論では、自由粒子のエネルギー固有状態は平面波となる。また自由粒子のハミルトニアンと運動量が可換であるため、運動量の固有状態も平面波である。つまりエネルギーと運動量についての同時固有関数となっている。量子論においても平面波は、基底関数として様々な場面で用いられるが、本来1に規格化されるべき2乗積分が有限の値を持たないこと、時間的・空間的に無限の彼方まで広がっており非現実的であること等の問題も抱えている。

## 第一原理バンド計算における平面波

波動関数は、基底関数で展開した形で記述することができる。この時に用いられる基底の1つに平面波基底（）がある。バンド計算における表式化が比較的簡単で（それ故、プログラムも構築し易い）力やストレスの計算も他の基底（局在基底など）を使った場合より容易に実現が可能である。また、平面波基底では、Pulay補正項の問題が回避できることも利点のひとつである。

欠点として、例えば波動関数や電荷密度への寄与の s, p, d 軌道毎への分割や、ユニットセル内の特定の原子の電荷を求めることが困難になることが挙げられる。

## 脚注

### 注釈

### 出典

## 関連項目
*球面波
*直交化された平面波
*第一原理バンド計算
*局在基底
*混合基底