フォック状態

量子力学において**フォック状態**（フォックじょうたい、）または**数状態**（すうじょうたい、）、または**粒子数状態**（りゅうしすうじょうたい）とは、粒子（または量子）の数が明確に定義されたフォック空間のベクトルである量子状態のこと。ソビエトの物理学者ウラジミール・フォックにちなんで名づけられた。
また多体系や量子場をフォック状態で表すことを**フォック表示**、**占有数表示**などと呼ぶ。量子光学では**光子数状態**あるいは**光子数確定状態**とも呼ばれる。

フォック状態は量子力学の第二量子化形式において重要な役割を果たす。

粒子表現は、ポール・ディラックがボース粒子について、パスクアル・ヨルダンとユージン・ウィグナーがフェルミ粒子について詳細に扱ったのが最初である。

## 1つのモードの場合
### ボース粒子の場合

最も単純な１つのモードのボース粒子フォック状態を考える（ここで考えている状態を**1粒子状態**と呼ぶこともあるが、この状態を多粒子状態と考えることもできるため（後述）、ここでは単にモードと呼ぶことにする）。これはエネルギー的には１つの調和振動子と等価である。

以下の交換関係を満たす生成消滅演算子を定義する。
: [\hat{b}, \hat{b}^\dagger ] = 1 
数演算子とよばれるエルミート演算子を以下で定義する。
:\hat{N} \equiv \hat{b}^\dagger \hat{b}
この固有値方程式は、
:\hat{N}|n\rangle=n|n\rangle
この固有値Nは非負の整数である。また、この固有ベクトル|n\rangleに生成消滅演算子が作用すると、
:\begin{align}
\hat{b}^{\dagger}|n\rangle &= \sqrt{n +1} |n+1\rangle, \\
\hat{b}|n\rangle &= \sqrt{n} |n-1\rangle
\end{align}
また固有値0の場合の固有状態|0\rangleを真空状態と呼び、これに消滅演算子が作用すると次のようになる。
:\hat{b}|0\rangle=0
数演算子の固有状態は、真空状態から生成演算子をくり返し作用することで作ることができる。これをボース粒子における**フォック状態**と呼ぶ。
:|n\rangle=\frac{1}{\sqrt{n!}}(\hat{b}^\dagger)^n|0\rangle
調和振動子では、数演算子とハミルトニアンは互いに交換する。
:\hat{H}=\hbar\omega\left(\hat{N}+\frac{1}{2}\right)
:[\hat{N},\hat{H}]=0
よってフォック状態|n\rangleは調和振動子のハミルトニアンの固有状態:|\psi_n\rangleでもある（同時固有状態）。
:\hat{H} | n \rangle = \hbar\omega\left(n+\frac{1}{2}\right)| n \rangle
つまり状態|n\rangleにおける粒子数とエネルギーの測定値には量子的なバラつき（ゆらぎ）は無い。

#### 粒子的な解釈
以上のことから、次のような再解釈を行うことができる。
*\hat{b}^\dagger, \hat{b}は、あるエネルギー（調和振動子では\hbar\omega）を持つ粒子を生成・消滅させる演算子である。
*ハミルトニアンの固有値・固有状態を表すnはその粒子の**粒子数**（あるいは**占有数**）であり、演算子\hat{N}は粒子数を表すオブザーバブルである。
このように考えることで元々は1粒子状態であった調和振動子の状態も、その励起数だけ粒子がある多粒子状態となる。

### フェルミ粒子の場合
以下の反交換関係を満たす生成消滅演算子を定義する。
:\{\hat{c},\hat{c}^\dagger\}=1
:\{\hat{c},\hat{c}\}=\{\hat{c}^\dagger,\hat{c}^\dagger\}=0
また数演算子とよばれるエルミート演算子を以下で定義する。
:\hat{N} \equiv \hat{c}^\dagger \hat{c}
この固有値方程式は、
:\hat{N}|n\rangle=n|n\rangle
真空状態を以下で定義する。
:\hat{c}|0\rangle=0
真空状態 |0\rangleに生成演算子をくり返し作用させると、\{\hat{c}^{\dagger},\hat{c}^{\dagger}\}=2(\hat{c}^{\dagger})^2=0より、
:\hat{c}^{\dagger}| 0 \rangle = | 1 \rangle
:(\hat{c}^{\dagger})^2| 0 \rangle =\hat{c}^{\dagger} | 1 \rangle= 0
つまり最大占有数は1で、1個以上のフェルミ粒子は同じ状態を占有できない。

逆に |1\rangleにくり返し消滅演算子を作用させると、\{\hat{c},\hat{c}\}=2(\hat{c})^2=0より
:\hat{c}| 1 \rangle = | 0 \rangle
:(\hat{c})^2|1\rangle=\hat{c}| 0 \rangle = 0
つまり粒子数は0以下になれない。この|0\rangle,|1\rangleをフェルミ粒子におけるフォック状態と呼ぶ。

ボース粒子の場合と同様に、この場合も粒子的な解釈を行うことができる。ただし粒子の生成は1つまでしかできない。

## 同種粒子の不可弁別性
1つのモードでのフォック状態を合成し、2つ以上のモード（多重モード）でのフォック状態を作る。このフォック状態は、同種粒子の不可弁別性を満たすような形でなければならない。

### 2つのモードの場合
合成した固有状態に交換演算子\hat{P}が作用すると、ボース粒子の場合は対称的、フェルミ粒子の場合は反対称的になければならない。

たとえばテンソル積表現での2粒子系においては次のようになる。
:\hat{P}\left|x_1, x_2\right\rangle  = \left|x_2, x_1\right\rangle

2つの同種粒子のフォック状態|1_, 1_\rangleについて、次の不可弁別性が成り立つ。よって
* ボース粒子では|1_, 1_\rangle=+|1_, 1_\rangle
* フェルミ粒子では|1_, 1_\rangle=-|1_, 1_\rangle
となる。

:

ここで数演算子は、ボース粒子とフェルミ粒子を区別せずに（つまり対称性を考慮せずに）粒子を数えるだけの演算子であることに注意。
この２種類の粒子の差を見るには、生成消滅演算子が必要となる。

このように複数の同種モードがある場合、ボース粒子のフォック状態は対称性を、フェルミ粒子のフォック状態は反対称性を持たなければならない。
これらの性質を満たすため、以下で述べるように、多重モードのフォック状態をテンソル積を用いて構成する。
テンソル積は、粒子がフェルミ粒子かボース粒子かによって、基となる1粒子ヒルベルト空間の交代積または対称積でなければならない。

## 多重モードのボース粒子フォック状態

### 生成消滅演算子の導入

この新しいフォック空間表現では、同じ対称性の性質を表現しなければならない。

そのため各モードiのボース粒子の生成消滅演算子を、以下の交換関係を満たす演算子\hat{b}_i^\dagger、\hat{b}_iとして定義する。

:[\hat{b}_i, \hat{b}^\dagger_j] \equiv \hat{b}_i \hat{b}^\dagger_j - \hat{b}^\dagger_j \hat{b}_i = \delta_{i j}
:[\hat{b}^\dagger_i, \hat{b}^\dagger_j] = [\hat{b}_i, \hat{b}_j] = 0 
ここで\delta_{i j}はクロネッカーのデルタである。生成消滅演算子はエルミート演算子ではない。

:

### 粒子数演算子
各モードの粒子数演算子を次のように定義する。
:\hat{N}_i\equiv \hat{b}_i^\dagger\hat{b}_i
これはエルミート演算子である。この固有値・固有ベクトルは以下を満たす。
:\hat{N}_i|n_i\rangle=n_i|n_i\rangle
また全数演算子を、各モードiの数演算子の和として定義する。
:\hat{n} = \sum_i \hat{n}_{i} 
### 状態のテンソル積
多重モードにおける**フォック状態**| n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangleを、全数演算子の固有ベクトルと定義する。固有値は全てのモードの粒子数（占有数）の和である。
:\hat{n} | n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle = \left( \sum_i n_{i} \right) | n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle

多重モードフォック状態は、各モードでの|n_i\rangle  の直積（テンソル積）で表される。
:| n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle
\equiv | n_1 \rangle | n_{2} \rangle \cdots| n_{i} \rangle \cdots


ここで各モードにおける粒子的な解釈を行うことで、多重モード全体の状態を各モードの**粒子数**（あるいは**占有数**）を用いて表すことができる。

### 生成消滅演算子の作用
多重モードのフォック状態におけるそれぞれの生成消滅演算子は、それ自身のモードにのみ作用する。たとえば\hat{b}_{i}と\hat{b}^\dagger_{i}は、|n_i\rangleにだけ作用する。

:\begin{align}
\hat{b}_{i} |n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle
&= \sqrt{n_{i}} | n_{1}, n_{2}, \cdots, n_{i} - 1, \cdots \rangle
\\
\hat{b}^\dagger_{i} | n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle
&= \sqrt{n_{i} + 1} | n_{1}, n_{2}, \cdots, n_{i} + 1, \cdots \rangle
\end{align}  

つまり生成消滅演算子の多重モード状態への作用は、それら自身のモードの粒子数（あるいは占有数）を１だけ増加または減少させるだけである。異なるモードに対応する演算子はヒルベルト空間の異なる部分空間に作用する。

たとえば、
* 真空状態(どの状態にも粒子が無い状態） |0_{1}, 0_{2}, \cdots, 0_{i}, \cdots \rangleにモードiの生成消滅演算子が作用すると、
::\begin{align}
\hat{b}^{\dagger}_{i} | 0_{1}, 0_{2}, \cdots, 0_{i}, \cdots \rangle
&= | 0_{1}, 0_{2}, \cdots, 1_{i}, \cdots \rangle \\
\hat{b}_{l} | 0_{1}, 0_{2}, \cdots, 0_{i}, \cdots \rangle
&=0
\end{align} 

* 生成演算子を真空状態に作用させることで、どんなフォック状態も作ることができる。
:: |n_1, n_2, \cdots\rangle=\frac{(b^{\dagger}_1)^{n_1}}{\sqrt{n_1!}}\frac{(b^{\dagger}_2)^{n_2}}{\sqrt{n_2!}}\cdots|0_1, 0_2 ,\cdots\rangle

### 数演算子の作用
*i*番目のモードの粒子数演算子\hat{N}_{i}は、ボース粒子フォック状態に以下のように作用する。
:\hat{N}_{i} | n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle
= n_{i} | n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle
つまりフォック状態も粒子数演算子の固有ベクトルになっている。

フォック状態は全粒子数の固有状態であるため、全粒子数の測定値の分散は\operatorname{Var}(\hat{N})=0となる。すなわちフォック状態における全粒子数の測定は、常にゆらぎが無く確定値を与える。
よって、
:\hat{N}=\sum_i \hat{N}_{i}

### 数演算子と交換するハミルトニアン
各モード間の相互作用が無い系を考える。この系の全ハミルトニアン\hat{H}は、各モードのハミルトニアン\hat{H}_iの和で表される。
:\hat{H}=\sum_i\hat{H}_i
全ハミルトニアンと全数演算子は交換する。よって多重モードフォック状態は多重モードハミルトニアンの固有状態でもある（同時固有状態）。
:\hat{H} | n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle = \left( \sum_i \hbar \omega \left(n_{i} + \frac{1}{2} \right) \right) | n_{1}, n_{2}, \cdots, n_{i}, \cdots \rangle

各モードで考えると、モードiのハミルトニアン\hat{H}_iの固有状態は、粒子数演算子\hat{N}_{i}の固有状態にもなっている（厳密にはさらに位相因子を選択しなければならない）。

ここで粒子的な解釈を行うことで、\hat{N}_{i}の固有状態に対応する固有値n_iは、モードiのハミルトニアン\hat{H}_iのn番目の固有状態|\psi_n\rangleにおける粒子数を与える。

### フォック空間
フォック状態はエルミート演算子である粒子数演算子の固有ベクトルであるため、正規直交基底をなす。
このフォック状態（とそれらの線形結合）から成る空間を**[[フォック空間]]**という。
フォック空間での基底であるフォック状態は、「占有数基底」とも呼ばれる。
フォック空間は、それぞれの粒子数におけるテンソル積ヒルベルト空間の直和となる。

フォック空間のベクトルの中で、粒子数が異なる状態の重ね合わせであるもの（たとえばコヒーレント状態など）は、数演算子の固有状態ではないためフォック状態ではない。よってフォック空間の全てのベクトルが「フォック状態」と呼ばれる訳ではない。

### ボース粒子フォック状態の対称性
ボース粒子フォック状態が粒子の交換によって対称性を示すことを確認する。

ここで２つの状態間の粒子の交換は、ある状態でのある粒子を消滅させ、別の状態での粒子を生成させることで行われる。
フォック状態|\cdots, n_m ,\cdots, n_l ,\cdots \rangle
から出発して状態lから状態mへの粒子をシフトさせたい場合は、交換関係より\hat{b}_{m}^\dagger \hat{b}_{l}
= \hat{b}_{l} \hat{b}_{m}^\daggerであるため、

:\begin{align}
\hat{b}_{m}^\dagger \hat{b}_{l} |  \cdots, n_{m}, \cdots, n_{l}, \cdots \rangle
&= \hat{b}_{l} \hat{b}_{m}^\dagger | \cdots, n_{m}, \cdots, n_{l}, \cdots \rangle \\
&=\sqrt{n_{m} + 1} \sqrt{n_{l}} |  \cdots, n_{m} + 1, \cdots, n_{l} - 1, \cdots \rangle
\end{align}

よってボース粒子フォック状態は粒子の交換において対称的である。

## 多重モードのフェルミ粒子フォック状態
### フェルミ粒子の生成消滅演算子
フェルミ粒子の反対称性を保持するために、フェルミ粒子の生成消滅演算子を導入する。
フェルミ粒子フォック状態  |n_{1}, n_{2}, \cdots, n_{l}, \cdots \rangleに**生成演算子**を次のように作用させる。
:\hat{c}^\dagger_{l} | n_{1}, n_{2}, \cdots, n_{l}, \cdots \rangle
= \sqrt{n_{l} +1} | n_{1}, n_{2}, \cdots, n_{l} + 1, \cdots \rangle  
**消滅演算子**は次のように作用する。
:\hat{c}_{l} | n_{1}, n_{2}, \cdots, n_{l}, \cdots \rangle
= \sqrt{n_{l}} | n_{1}, n_{2}, \cdots, n_{l} - 1, \cdots \rangle  
これら２つの演算子の作用は反対称的に行われる（後述）。

#### 演算子の反交換関係
フェルミ粒子系における生成消滅演算子の反交換関係は、
:\{ \hat{c}_i, \hat{c}^\dagger_j \} \equiv \hat{c}_i \hat{c}^\dagger_j + \hat{c}^\dagger_j \hat{c}_i = \delta_{i j}
:\{ \hat{c}^\dagger_i, \hat{c}^\dagger_j \} = \{ \hat{c}_i, \hat{c}_j \} = 0 
ここで{\{ \ , \ \} } は反交換子、\delta_{i j}はクロネッカーのデルタである。これらの反交換関係は、フェルミ粒子フォック状態の反対称性を表すために用いられる。

### 数演算子の作用
フェルミ粒子の数演算子は次のように与えられる。
:\hat{N}_l = \hat{c}^\dagger_{l} \hat{c}_{l}
これがフェルミ粒子フォック状態に作用すると、
:\hat{N}_l |n_{1}, n_{2}, \cdots, n_{l}, \cdots \rangle
= n_{l} | n_{1}, n_{2}, \cdots, n_{l}, \cdots \rangle  

#### 最大占有数
生成消滅演算子や数演算子の作用はボース粒子の場合と同じであるように見えるが、フェルミ粒子フォック状態の最大占有数から違いが生じる。
上述の２つのフェルミ粒子の例を拡張し、
次のように置換演算子のある特定の和を固有ケットのテンソル積に適用することで、
フェルミ粒子フォック状態| n_{1}, n_{2}, \cdots, n_{l}, \cdots \rangleが得られることをまず確かめておかなければならない。

:|n_{1}, n_{2}, \cdots, n_{l}, \cdots \rangle
= S_- |i_1, i_2, \cdots, i_l, \cdots \rangle
= \frac{1}{\sqrt{N!}} \begin{vmatrix}
|i_1\rangle_1 & \cdots & |i_1\rangle_{N} \\
\vdots & \ddots & \vdots \\
|i_{N}\rangle_1 & \cdots & |i_{N}\rangle_{N}
\end{vmatrix} 
この行列式はスレーター行列式と呼ばれる。
もし1粒子状態でも同じものがあればスレーター行列式の2つの行は同じであり、行列式は0になる。
これは2つの同種フェルミ粒子が同じ状態を占めないことをあらわしている。
よっていかなる単一状態の占有数も0または1のどちらかである。
フェルミ粒子フォック状態に関連する固有値\hat{N}_{l}は0または1である。

### フェルミ粒子フォック空間の基底

### 具体的なフォック状態への作用
| n_{1}, n_{2}, \cdots, n_{\beta}, n_{\alpha}, \cdots \rangleで表される多重モードフェルミ粒子フォック状態では、
::\hat{c}_ | n_{1}, n_{2}, \cdots, n_, n_, \cdots \rangle
= (-1)^{\sum_{\beta 
ここで(-1)^{\sum_{\beta は、ヨルダン-ウィグナーストリングと呼ばれ、含まれる1粒子状態の順序に依存し、全ての前に来る状態のフェルミ粒子占有数を足し合わせる。

### フェルミ粒子フォック状態の反対称性
交換演算子のもとでのフェルミ粒子状態の反対称性は、反交換関係によるものである。
ここでは2つの状態間の粒子の交換は、ある状態のある粒子を消滅させ、別の粒子を生成させることで行われる。
フォック状態| n_{1}, n_{2}, \cdots, n_{m} \cdots, n_{l}, \cdots \rangleから出発して状態lから状態 mへ粒子をシフトさせたい場合、反交換関係より\hat{c}_{m}^\dagger \hat{c}_{l} = - \hat{c}_{l} \hat{c}_{m}^\daggerであるため、
:\hat{c}_{m}^\dagger \hat{c}_{l} | n_{1}, n_{2}, \cdots, n_{m}, \cdots, n_{l}, \cdots \rangle
= \sqrt{n_{m} +1} \sqrt{n_{l}} | n_{1}, n_{2}, \cdots, n_{m} + 1, \cdots, n_{l} - 1, \cdots \rangle
しかし、 
:\begin{align}
\hat{c}_{l} \hat{c}_{m}^\dagger | n_{1}, n_{2}, \cdots, n_{m}, \cdots, n_{l}, \cdots \rangle
&= - \hat{c}_{m}^\dagger \hat{c}_{l} | n_{1}, n_{2}, \cdots, n_{m}, \cdots, n_{l}, \cdots \rangle \\
&= - \sqrt{n_{m} + 1} \sqrt{n_{l}} | n_{1}, n_{2}, \cdots, n_{m} + 1, \cdots, n_{l} - 1, \cdots \rangle
\end{align}
よってフェルミ粒子フォック状態は、粒子の交換において反対称性である。

## フォック状態は一般的にエネルギー固有状態では無い
第二量子化理論では、ハミルトニアン密度関数は次のように与えられる。
:\mathfrak{H} = \frac{1}{2m} \nabla_{i}\psi^{*}(x)\nabla_{i}\psi(x) 
全ハミルトニアンは、次のように与えられる。
:\mathcal{H} = \int d^{3}x\,\mathfrak{H}=\int d^{3}x \psi^{*}(x)\left(-\tfrac{\nabla^2}{2m}\right)\psi(x) \quad (\therefore \mathfrak{H}=-\tfrac{\nabla^2}{2m})

自由粒子のシュレーディンガー方程式 は、
:\mathfrak{H}\psi_{n}^{(+)}(x)=-\tfrac{\nabla^2}{2m}\psi_{n}^{(+)}(x)=E_{n}^0\psi_{n}^{(+)}(x)
この解は直交性を満たす。
:\int d^{3}x\, \psi_{n}^{(+)^{*}}(x)\psi_{n'}^{(+)}(x)=\delta_{nn'}
また消滅演算子を\hat{a}_nとして、次の関係がある。
:\psi(x)=\sum_n \hat{a}_n \psi_{n}^{(+)}(x)
よって、
:\therefore \mathcal{H} = \sum_{n,n'}\int d^{3}x \, \hat{a}^{\dagger}_{n'}\psi_{n'}^{(+)^{*}}(x)\mathfrak{H} \hat{a}_n \psi_{n}^{(+)}(x)

相互作用しない粒子においてのみ、\mathfrak{H}と\hat{a}_nが交換する。しかし一般の場合にはこれらは交換しない。
相互作用しない粒子では、
:\mathcal{H}=\sum_{n,n'}\int d^{3}x \, \hat{a}^{\dagger}_{n'}\psi_{n'}^{(+)^{*}}(x)E^0_{n}\psi_{n}^{(+)}(x) \hat{a}_n 
= \sum_{n,n'}E^0_{n} \hat{a}^{\dagger}_{n'} \hat{a}_n\delta_{nn'}
= \sum_{n}E^0_{n} \hat{a}^{\dagger}_{n} \hat{a}_n
= \sum_{n}E^0_{n} \hat{N}
これらが交換しない場合、ハミルトニアンは上述の表現を持たない。
よって一般的にフォック状態は系のエネルギー固有状態ではない。

## 真空ゆらぎ
真空状態|0\rangleは最低エネルギーの状態で、\hat{a}と\hat{a}^{\dagger}の期待値はこの状態では0になる。
:\hat{a} | 0 \rangle = \langle 0 | \hat{a}^{\dagger} = 0
電磁場とベクトルポテンシャルは同じ一般形のモード展開を持つ。
:F(\boldsymbol{r}, t) = \varepsilon a e^{i \boldsymbol{\psi} \cdot \boldsymbol{r} - \omega t} + hc 
これらの場の演算子の期待値が真空状態では0になることを見るのは簡単である。
:\langle0|F|0\rangle = 0 
しかし、これらの場の演算子の二乗の期待値は0ではないことを示すことができる。
このように0アンサンブル平均についての場におけるゆらぎが存在する。
これらの真空ゆらぎは、量子光学におけるラムシフトなど多くの興味深い現象の原因となる。

## 単一フォトン状態の源
単一フォトンは、シングルエミッター(原子、窒素-空孔中心、 量子ドット)を用いてごく普通に生成される。
しかしこれらの源はいつも効率的であるとは限らず(要求に応じた単一フォトンを実際に得る確率は低い)、しばしば複雑で、実験環境からは適さない。
非決定的ふるまいを犠牲にして、これらの問題を克服する別の源が一般的に用いられる。
歓迎された単一フォトン源は確率的２フォトン源で、そこからペアが分かれ、1フォトンの検出は残りのフォトンの存在を歓迎する。
これらの源はたとえば周期的に分極したニオブ酸リチウム(自発的パラメトリック下方変換)やシリコン(自発的四波混合)など、通常いくつかの材料の非線形光学作用に依存している。

## 非古典的ふるまい
フォック状態のは、これらの状態が純粋に量子力学的な状態であり、古典的対応物は存在しないことを示している。
この表現におけるこれらの状態の\varphi(\alpha)はディラックのデルタ関数の2n次導関数であり、よって古典的な確率分布ではない。

## 脚注

## 関連項目
*コヒーレント状態
*ハイゼンベルク極限
*非古典光

## 外部リンク
*Vladan Vuletic of MIT has used an ensemble of atoms to produce a Fock state (a.k.a. single photon) source (PDF)
* Produce and measure a single photon state (Fock state) with an interactive experiment QuantumLab