Fock state

One specifies a multiparticle state of N non-interacting identical particles by writing the state as a sum of tensor products of N one-particle states. Additionally, depending on the integrality of the particles' spin, the tensor products must be alternating (anti-symmetric) or symmetric products of the underlying one-particle Hilbert space. Specifically:

• Fermions, having half-integer spin and obeying the Pauli exclusion principle, correspond to antisymmetric tensor products.
• Bosons, possessing integer spin (and not governed by the exclusion principle) correspond to symmetric tensor products.

If the number of particles is variable, one constructs the Fock space as the direct sum of the tensor product Hilbert spaces for each particle number. In the Fock space, it is possible to specify the same state in a new notation, the occupancy number notation, by specifying the number of particles in each possible one-particle state.

Let ${\textstyle \left\{\mathbf {k} _{i}\right\}_{i\in I}}$  be an orthonormal basis of states in the underlying one-particle Hilbert space. This induces a corresponding basis of the Fock space called the "occupancy number basis". A quantum state in the Fock space is called a Fock state if it is an element of the occupancy number basis.

A Fock state satisfies an important criterion: for each i, the state is an eigenstate of the particle number operator ${\displaystyle {\widehat {N_{{\mathbf {k} }_{i}}}}}$  corresponding to the i-th elementary state k_i. The corresponding eigenvalue gives the number of particles in the state. This criterion nearly defines the Fock states (one must in addition select a phase factor).

A given Fock state is denoted by ${\displaystyle |n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},..n_{{\mathbf {k} }_{i}}...\rangle }$ . In this expression, ${\displaystyle n_{{\mathbf {k} }_{i}}}$  denotes the number of particles in the i-th state k_i, and the particle number operator for the i-th state, ${\displaystyle {\widehat {N_{{\mathbf {k} }_{i}}}}}$ , acts on the Fock state in the following way:

${\displaystyle {\widehat {N_{{\mathbf {k} }_{i}}}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},..n_{{\mathbf {k} }_{i}}...\rangle =n_{{\mathbf {k} }_{i}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},..n_{{\mathbf {k} }_{i}}...\rangle }$ 

Hence the Fock state is an eigenstate of the number operator with eigenvalue ${\displaystyle n_{{\mathbf {k} }_{i}}}$ .^[2]: 478

Fock states often form the most convenient basis of a Fock space. Elements of a Fock space that are superpositions of states of differing particle number (and thus not eigenstates of the number operator) are not Fock states. For this reason, not all elements of a Fock space are referred to as "Fock states".

If we define the aggregate particle number operator ${\textstyle {\widehat {N}}}$  as

${\displaystyle {\widehat {N}}=\sum _{i}{\widehat {N_{{\mathbf {k} }_{i}}}},}$ 

the definition of Fock state ensures that the variance of measurement ${\displaystyle \operatorname {Var} \left({\widehat {N}}\right)=0}$ , i.e., measuring the number of particles in a Fock state always returns a definite value with no fluctuation.

For any final state ${\displaystyle |f\rangle }$ , any Fock state of two identical particles given by ${\displaystyle |1_{\mathbf {k} _{1}},1_{\mathbf {k} _{2}}\rangle }$ , and any operator ${\displaystyle {\widehat {\mathbb {O} }}}$ , we have the following condition for indistinguishability:^[3]: 191

${\displaystyle \left|\left\langle f\left|{\widehat {\mathbb {O} }}\right|1_{\mathbf {k} _{1}},1_{\mathbf {k} _{2}}\right\rangle \right|^{2}=\left|\left\langle f\left|{\widehat {\mathbb {O} }}\right|1_{\mathbf {k} _{2}},1_{\mathbf {k} _{1}}\right\rangle \right|^{2}}$ .

So, we must have ${\displaystyle \left\langle f\left|{\widehat {\mathbb {O} }}\right|1_{\mathbf {k} _{1}},1_{\mathbf {k} _{2}}\right\rangle =e^{i\delta }\left\langle f\left|{\widehat {\mathbb {O} }}\right|1_{\mathbf {k} _{2}},1_{\mathbf {k} _{1}}\right\rangle }$ 

where ${\displaystyle e^{i\delta }=+1}$  for bosons and ${\displaystyle -1}$  for fermions. Since ${\displaystyle \langle f|}$  and ${\displaystyle {\widehat {\mathbb {O} }}}$  are arbitrary, we can say,

${\displaystyle \left|1_{\mathbf {k} _{1}},1_{\mathbf {k} _{2}}\right\rangle =+\left|1_{\mathbf {k} _{2}},1_{\mathbf {k} _{1}}\right\rangle }$  for bosons and

${\displaystyle \left|1_{\mathbf {k} _{1}},1_{\mathbf {k} _{2}}\right\rangle =-\left|1_{\mathbf {k} _{2}},1_{\mathbf {k} _{1}}\right\rangle }$  for fermions.^[3]: 191

Note that the number operator does not distinguish bosons from fermions; indeed, it just counts particles without regard to their symmetry type. To perceive any difference between them, we need other operators, namely the creation and annihilation operators.

Bosons, which are particles with integer spin, follow a simple rule: their composite eigenstate is symmetric^[4] under operation by an exchange operator. For example, in a two particle system in the tensor product representation we have ${\displaystyle {\hat {P}}\left|x_{1},x_{2}\right\rangle =\left|x_{2},x_{1}\right\rangle }$  .

We should be able to express the same symmetric property in this new Fock space representation. For this we introduce non-Hermitian bosonic creation and annihilation operators,^[4] denoted by ${\displaystyle b^{\dagger }}$  and ${\displaystyle b}$  respectively. The action of these operators on a Fock state are given by the following two equations:

• Creation operator ${\textstyle b_{{\mathbf {k} }_{l}}^{\dagger }}$ :

  ${\displaystyle b_{{\mathbf {k} }_{l}}^{\dagger }|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle ={\sqrt {n_{{\mathbf {k} }_{l}}+1}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}+1,...\rangle }$ ^[4]

• Annihilation operator ${\textstyle b_{{\mathbf {k} }_{l}}}$ :

  ${\displaystyle b_{{\mathbf {k} }_{l}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle ={\sqrt {n_{{\mathbf {k} }_{l}}}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}-1,...\rangle }$ ^[4]

 

The bosonic Fock state creation and annihilation operators are not Hermitian operators.^[4]

The commutation relations of creation and annihilation operators in a bosonic system are

${\displaystyle \left[b_{i}^{\,},b_{j}^{\dagger }\right]\equiv b_{i}^{\,}b_{j}^{\dagger }-b_{j}^{\dagger }b_{i}^{\,}=\delta _{ij},}$ ^[4]

${\displaystyle \left[b_{i}^{\dagger },b_{j}^{\dagger }\right]=\left[b_{i}^{\,},b_{j}^{\,}\right]=0,}$ ^[4]

where ${\displaystyle [\ \ ,\ \ ]}$  is the commutator and ${\displaystyle \delta _{ij}}$  is the Kronecker delta.

Number of particles (N)	Bosonic basis states[6]: 11
0	${\displaystyle |0,0,0...\rangle }$
1	${\displaystyle |1,0,0...\rangle }$ , ${\displaystyle |0,1,0...\rangle }$ , ${\displaystyle |0,0,1...\rangle }$ ,...
2	${\displaystyle |2,0,0...\rangle }$ , ${\displaystyle |1,1,0...\rangle }$ , ${\displaystyle |0,2,0...\rangle }$ ,...
${\displaystyle n}$	${\displaystyle |n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle }$

The number operators ${\textstyle {\widehat {N_{{\mathbf {k} }_{l}}}}}$  for a bosonic system are given by ${\displaystyle {\widehat {N_{{\mathbf {k} }_{l}}}}=b_{{\mathbf {k} }_{l}}^{\dagger }b_{{\mathbf {k} }_{l}}}$ , where ${\displaystyle {\widehat {N_{{\mathbf {k} }_{l}}}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}...\rangle =n_{{\mathbf {k} }_{l}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}...\rangle }$ ^[4]

Number operators are Hermitian operators.

The commutation relations of the creation and annihilation operators ensure that the bosonic Fock states have the appropriate symmetric behaviour under particle exchange. Here, exchange of particles between two states (say, l and m) is done by annihilating a particle in state l and creating one in state m. If we start with a Fock state ${\displaystyle |\psi \rangle =\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},....n_{\mathbf {k} _{m}}...n_{\mathbf {k} _{l}}...\right\rangle }$ , and want to shift a particle from state ${\displaystyle k_{l}}$  to state ${\displaystyle k_{m}}$ , then we operate the Fock state by ${\displaystyle b_{\mathbf {k} _{m}}^{\dagger }b_{\mathbf {k} _{l}}}$  in the following way:

Using the commutation relation we have, ${\displaystyle b_{\mathbf {k} _{m}}^{\dagger }.b_{\mathbf {k} _{l}}=b_{\mathbf {k} _{l}}.b_{\mathbf {k} _{m}}^{\dagger }}$ 

${\displaystyle {\begin{aligned}b_{\mathbf {k} _{m}}^{\dagger }.b_{\mathbf {k} _{l}}\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},....n_{\mathbf {k} _{m}}...n_{\mathbf {k} _{l}}...\right\rangle &=b_{\mathbf {k} _{l}}.b_{\mathbf {k} _{m}}^{\dagger }\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},....n_{\mathbf {k} _{m}}...n_{\mathbf {k} _{l}}...\right\rangle \\&={\sqrt {n_{\mathbf {k} _{m}}+1}}{\sqrt {n_{\mathbf {k} _{l}}}}\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},....n_{\mathbf {k} _{m}}+1...n_{\mathbf {k} _{l}}-1...\right\rangle \end{aligned}}}$ 

So, the Bosonic Fock state behaves to be symmetric under operation by Exchange operator.

•  
  Wigner function of ${\displaystyle |0\rangle }$ 
•  
  Wigner function of ${\displaystyle |1\rangle }$ 
•  
  Wigner function of ${\displaystyle |2\rangle }$ 
•  
  Wigner function of ${\displaystyle |3\rangle }$ 
•  
  Wigner function of ${\displaystyle |4\rangle }$ 

To be able to retain the antisymmetric behaviour of fermions, for Fermionic Fock states we introduce non-Hermitian fermion creation and annihilation operators,^[4] defined for a Fermionic Fock state ${\displaystyle |\psi \rangle =|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle }$  as:^[4]

• The creation operator ${\displaystyle c_{{\mathbf {k} }_{l}}^{\dagger }}$  acts as:

  ${\displaystyle c_{{\mathbf {k} }_{l}}^{\dagger }|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle ={\sqrt {n_{{\mathbf {k} }_{l}}+1}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}+1,...\rangle }$ ^[4]

• The annihilation operator ${\textstyle c_{{\mathbf {k} }_{l}}}$  acts as:

  ${\displaystyle c_{{\mathbf {k} }_{l}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle ={\sqrt {n_{{\mathbf {k} }_{l}}}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}-1,...\rangle }$ 

These two actions are done antisymmetrically, which we shall discuss later.

The anticommutation relations of creation and annihilation operators in a fermionic system are,

${\displaystyle {\begin{aligned}\left\{c_{i}^{\,},c_{j}^{\dagger }\right\}\equiv c_{i}^{\,}c_{j}^{\dagger }+c_{j}^{\dagger }c_{i}^{\,}&=\delta _{ij},\\\left\{c_{i}^{\dagger },c_{j}^{\dagger }\right\}=\left\{c_{i}^{\,},c_{j}^{\,}\right\}&=0,\end{aligned}}}$  ^[4]

where ${\displaystyle {\{\ ,\ \}}}$  is the anticommutator and ${\displaystyle \delta _{ij}}$  is the Kronecker delta. These anticommutation relations can be used to show antisymmetric behaviour of Fermionic Fock states.

Number operators ${\textstyle {\widehat {N_{{\mathbf {k} }_{l}}}}}$  for Fermions are given by ${\displaystyle {\widehat {N_{{\mathbf {k} }_{l}}}}=c_{{\mathbf {k} }_{l}}^{\dagger }.c_{{\mathbf {k} }_{l}}}$ .

${\displaystyle {\widehat {N_{{\mathbf {k} }_{l}}}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}...\rangle =n_{{\mathbf {k} }_{l}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}...\rangle }$ ^[4]

The action of the number operator as well as the creation and annihilation operators might seem same as the bosonic ones, but the real twist comes from the maximum occupation number of each state in the fermionic Fock state. Extending the 2-particle fermionic example above, we first must convince ourselves that a fermionic Fock state ${\displaystyle |\psi \rangle =\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},n_{\mathbf {k} _{3}}...n_{\mathbf {k} _{l}}...\right\rangle }$  is obtained by applying a certain sum of permutation operators to the tensor product of eigenkets as follows:

${\displaystyle \left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},n_{\mathbf {k} _{3}}...n_{\mathbf {k} _{l}}...\right\rangle =S_{-}\left|i_{1},i_{2},i_{3}...i_{l}...\right\rangle ={\frac {1}{\sqrt {N!}}}{\begin{vmatrix}\left|i_{1}\right\rangle _{1}&\cdots &\left|i_{1}\right\rangle _{N}\\\vdots &\ddots &\vdots \\\left|i_{N}\right\rangle _{1}&\cdots &\left|i_{N}\right\rangle _{N}\end{vmatrix}}}$ ^[7]: 16

This determinant is called the Slater determinant.^{[citation needed]} If any of the single particle states are the same, two rows of the Slater determinant would be the same and hence the determinant would be zero. Hence, two identical fermions must not occupy the same state (a statement of the Pauli exclusion principle). Therefore, the occupation number of any single state is either 0 or 1. The eigenvalue associated to the fermionic Fock state ${\displaystyle {\widehat {N_{{\mathbf {k} }_{l}}}}}$  must be either 0 or 1.

Number of particles (N)	Fermionic basis states[6]: 11
0	${\displaystyle |0,0,0...\rangle }$
1	${\displaystyle |1,0,0...\rangle }$ , ${\displaystyle |0,1,0...\rangle }$ , ${\displaystyle |0,0,1...\rangle }$ ,...
2	${\displaystyle |1,1,0...\rangle }$ , ${\displaystyle |0,1,1...\rangle }$ , ${\displaystyle |0,1,0,1...\rangle }$ , ${\displaystyle |1,0,1,0...\rangle }$ ...
...	...

 

Antisymmetric behaviour of Fermionic states under Exchange operator is taken care of the anticommutation relations. Here, exchange of particles between two states is done by annihilating one particle in one state and creating one in other. If we start with a Fock state ${\displaystyle |\psi \rangle =\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},...n_{\mathbf {k} _{m}}...n_{\mathbf {k} _{l}}...\right\rangle }$  and want to shift a particle from state ${\displaystyle k_{l}}$  to state ${\displaystyle k_{m}}$ , then we operate the Fock state by ${\displaystyle c_{\mathbf {k} _{m}}^{\dagger }.c_{\mathbf {k} _{l}}}$  in the following way:

Using the anticommutation relation we have

${\displaystyle c_{\mathbf {k} _{m}}^{\dagger }.c_{\mathbf {k} _{l}}=-c_{\mathbf {k} _{l}}.c_{\mathbf {k} _{m}}^{\dagger }}$ 

${\displaystyle c_{\mathbf {k} _{m}}^{\dagger }.c_{\mathbf {k} _{l}}\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},....n_{\mathbf {k} _{m}}...n_{\mathbf {k} _{l}}...\right\rangle ={\sqrt {n_{\mathbf {k} _{m}}+1}}{\sqrt {n_{\mathbf {k} _{l}}}}\left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},....n_{\mathbf {k} _{m}}+1...n_{\mathbf {k} _{l}}-1...\right\rangle }$ 

but, $${\displaystyle {\begin{aligned}&c_{{\mathbf {k} }_{l}}.c_{{\mathbf {k} }_{m}}^{\dagger }|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},....n_{{\mathbf {k} }_{m}}...n_{{\mathbf {k} }_{l}}...\rangle \\={}-&c_{{\mathbf {k} }_{m}}^{\dagger }.c_{{\mathbf {k} }_{l}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},....n_{{\mathbf {k} }_{m}}...n_{{\mathbf {k} }_{l}}...\rangle \\={}-&{\sqrt {n_{{\mathbf {k} }_{m}}+1}}{\sqrt {n_{{\mathbf {k} }_{l}}}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},....n_{{\mathbf {k} }_{m}}+1...n_{{\mathbf {k} }_{l}}-1...\rangle \end{aligned}}}$$ 

Thus, fermionic Fock states are antisymmetric under operation by particle exchange operators.

In second quantization theory, the Hamiltonian density function is given by

${\displaystyle {\mathfrak {H}}={\frac {1}{2m}}\nabla _{i}\psi ^{*}(x)\,\nabla _{i}\psi (x)}$ ^[3]: 189

The total Hamiltonian is given by

${\displaystyle {\begin{aligned}{\mathcal {H}}&=\int d^{3}x\,{\mathfrak {H}}=\int d^{3}x\psi ^{*}(x)\left(-{\frac {\nabla ^{2}}{2m}}\right)\psi (x)\\\therefore {\mathfrak {H}}&=-{\frac {\nabla ^{2}}{2m}}\end{aligned}}}$ 

In free Schrödinger theory,^[3]: 189

${\displaystyle {\mathfrak {H}}\psi _{n}^{(+)}(x)=-{\frac {\nabla ^{2}}{2m}}\psi _{n}^{(+)}(x)=E_{n}^{0}\psi _{n}^{(+)}(x)}$ 

and

${\displaystyle \int d^{3}x\,\psi _{n}^{(+)^{*}}(x)\,\psi _{n'}^{(+)}(x)=\delta _{nn'}}$ 

and

${\displaystyle \psi (x)=\sum _{n}a_{n}\psi _{n}^{(+)}(x)}$ ,

where ${\displaystyle a_{n}}$  is the annihilation operator.

${\displaystyle \therefore {\mathcal {H}}=\sum _{n,n'}\int d^{3}x\,a_{n'}^{\dagger }\psi _{n'}^{(+)^{*}}(x)\,{\mathfrak {H}}a_{n}\psi _{n}^{(+)}(x)}$ 

Only for non-interacting particles do ${\displaystyle {\mathfrak {H}}}$  and ${\displaystyle a_{n}}$  commute; in general they do not commute. For non-interacting particles,

${\displaystyle {\mathcal {H}}=\sum _{n,n'}\int d^{3}x\,a_{n'}^{\dagger }\psi _{n'}^{(+)^{*}}(x)\,E_{n}^{0}\psi _{n}^{(+)}(x)a_{n}=\sum _{n,n'}E_{n}^{0}a_{n'}^{\dagger }a_{n}\delta _{nn'}=\sum _{n}E_{n}^{0}a_{n}^{\dagger }a_{n}=\sum _{n}E_{n}^{0}{\widehat {N}}}$ 

If they do not commute, the Hamiltonian will not have the above expression. Therefore, in general, Fock states are not energy eigenstates of a system.

The vacuum state or ${\displaystyle |0\rangle }$  is the state of lowest energy and the expectation values of ${\displaystyle a}$  and ${\displaystyle a^{\dagger }}$  vanish in this state:

${\displaystyle a|0\rangle =0=\langle 0|a^{\dagger }}$ 

The electrical and magnetic fields and the vector potential have the mode expansion of the same general form:

${\displaystyle F\left({\vec {r}},t\right)=\varepsilon ae^{i{\vec {k}}x-\omega t}+h.c.}$ 

Thus it is easy to see that the expectation values of these field operators vanishes in the vacuum state:

${\displaystyle \langle 0|F|0\rangle =0}$ 

However, it can be shown that the expectation values of the square of these field operators is non-zero. Thus there are fluctuations in the field about the zero ensemble average. These vacuum fluctuations are responsible for many interesting phenomenon including the Lamb shift in quantum optics.

In a multi-mode field each creation and annihilation operator operates on its own mode. So ${\displaystyle a_{\mathbf {k} _{l}}}$  and ${\displaystyle a_{\mathbf {k} _{l}}^{\dagger }}$  will operate only on ${\displaystyle \left|n_{\mathbf {k} _{l}}\right\rangle }$ . Since operators corresponding to different modes operate in different sub-spaces of the Hilbert space, the entire field is a direct product of ${\displaystyle |n_{\mathbf {k} _{l}}\rangle }$  over all the modes:

${\displaystyle \left|n_{\mathbf {k} _{1}}\right\rangle \left|n_{\mathbf {k} _{2}}\right\rangle \left|n_{\mathbf {k} _{3}}\right\rangle \ldots \equiv \left|n_{\mathbf {k} _{1}},n_{\mathbf {k} _{2}},n_{\mathbf {k} _{3}}...n_{\mathbf {k} _{l}}...\right\rangle \equiv \left|\{n_{\mathbf {k} }\}\right\rangle }$ 

The creation and annihilation operators operate on the multi-mode state by only raising or lowering the number state of their own mode:

${\displaystyle {\begin{aligned}a_{{\mathbf {k} }_{l}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle &={\sqrt {n_{{\mathbf {k} }_{l}}}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}-1,...\rangle \\a_{{\mathbf {k} }_{l}}^{\dagger }|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}},...\rangle &={\sqrt {n_{{\mathbf {k} }_{l}}+1}}|n_{{\mathbf {k} }_{1}},n_{{\mathbf {k} }_{2}},n_{{\mathbf {k} }_{3}}...n_{{\mathbf {k} }_{l}}+1,...\rangle \end{aligned}}}$ 

We also define the total number operator for the field which is a sum of number operators of each mode:

${\displaystyle {\hat {n}}_{\mathbf {k} }=\sum {\hat {n}}_{\mathbf {k} _{l}}}$ 

The multi-mode Fock state is an eigenvector of the total number operator whose eigenvalue is the total occupation number of all the modes

${\displaystyle {\hat {n}}_{\mathbf {k} }|\{n_{\mathbf {k} }\}\rangle =\left(\sum n_{\mathbf {k} _{l}}\right)|\{n_{\mathbf {k} }\}\rangle }$ 

In case of non-interacting particles, number operator and Hamiltonian commute with each other and hence multi-mode Fock states become eigenstates of the multi-mode Hamiltonian

${\displaystyle {\hat {H}}\left|\{n_{\mathbf {k} }\}\right\rangle =\left(\sum \hbar \omega \left(n_{\mathbf {k} _{l}}+{\frac {1}{2}}\right)\right)\left|\{n_{\mathbf {k} }\}\right\rangle }$ 

Single photons are routinely generated using single emitters (atoms, ions, molecules, Nitrogen-vacancy center,^[8] Quantum dot^[9]). However, these sources are not always very efficient, often presenting a low probability of actually getting a single photon on demand; and often complex and unsuitable out of a laboratory environment.

Other sources are commonly used that overcome these issues at the expense of a nondeterministic behavior. Heralded single photon sources are probabilistic two-photon sources from whom the pair is split and the detection of one photon heralds the presence of the remaining one. These sources usually rely on the optical non-linearity of some materials like periodically poled Lithium niobate (Spontaneous parametric down-conversion), or silicon (spontaneous Four-wave mixing) for example.

The Glauber–Sudarshan P-representation of Fock states shows that these states are purely quantum mechanical and have no classical counterpart. The ${\displaystyle \scriptstyle \varphi (\alpha )\,}$ ^{[clarification needed]} of these states in the representation is a ${\displaystyle 2n}$ 'th derivative of the Dirac delta function and therefore not a classical probability distribution.

• Coherent states
• Heisenberg limit
• Nonclassical light

• 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