Editing Normal order

{{Short description|Type of operator ordering in quantum field theory}}
{{Other uses}}

In [[quantum field theory]] a product of quantum fields, or equivalently their [[creation and annihilation operators]], is usually said to be '''normal ordered''' (also called '''Wick order''') when all creation operators are to the left of all annihilation operators in the product. The process of putting a product into normal order is called '''normal ordering''' (also called '''Wick ordering'''). The terms '''antinormal order''' and '''antinormal ordering''' are analogously defined, where the annihilation operators are placed to the left of the creation operators.

Normal ordering of a product of quantum fields or [[creation and annihilation operators]] can also be defined in many [[#Alternative definitions|other ways]].  Which definition is most appropriate depends on the expectation values needed for a given calculation.  Most of this article uses the most common definition of normal ordering as given above, which is appropriate when taking [[expectation value]]s using the vacuum state of the [[creation and annihilation operators]].

The process of normal ordering is particularly important for a [[quantum mechanics|quantum mechanical]] [[Hamiltonian (quantum mechanics)|Hamiltonian]]. When quantizing a [[classical mechanics|classical]] Hamiltonian there is some freedom when choosing the operator order, and these choices lead to differences in the [[Zero-point energy|ground state energy]]. That's why the process can also be used to eliminate the infinite vacuum energy of a quantum field.

==Notation==
If <math>\hat{O}</math> denotes an arbitrary product of creation and/or annihilation operators (or equivalently, quantum fields), then the normal ordered form of <math>\hat{O}</math> is denoted by <math>\mathopen{:} \hat{O} \mathclose{:}</math>.

An alternative notation is <math> \mathcal{N}(\hat{O})</math>.

Note that normal ordering is a concept that only makes sense for products of operators.  Attempting to apply normal ordering to a sum of operators is not useful as normal ordering is not a linear operation.

==Bosons==

[[Bosons]] are particles which satisfy [[Bose–Einstein statistics]]. We will now examine the normal ordering of bosonic creation and annihilation operator products.

===Single bosons===
If we start with only one type of boson there are two operators of interest:

* <math>\hat{b}^\dagger</math>: the boson's creation operator.
* <math>\hat{b}</math>: the boson's annihilation operator.

These satisfy the [[commutator]] relationship

:<math>\left[\hat{b}^\dagger, \hat{b}^\dagger \right]_- = 0</math>
:<math>\left[\hat{b}, \hat{b} \right]_- = 0</math>
:<math>\left[\hat{b}, \hat{b}^\dagger \right]_- = 1</math>

where <math>\left[ A, B \right]_- \equiv AB - BA</math> denotes the [[commutator]]. We may rewrite the last one as: <math>\hat{b}\, \hat{b}^\dagger = \hat{b}^\dagger\, \hat{b} + 1.</math>

====Examples====
1. We'll consider the simplest case first. This is the normal ordering of <math>\hat{b}^\dagger \hat{b}</math>:

:<math> {:\,}\hat{b}^\dagger \, \hat{b}{\,:} = \hat{b}^\dagger \, \hat{b}. </math>

The expression <math>\hat{b}^\dagger \, \hat{b}</math> has not been changed because it is ''already'' in normal order - the creation operator <math>(\hat{b}^\dagger)</math> is already to the left of the annihilation operator <math>(\hat{b})</math>.

2. A more interesting example is the normal ordering of <math>\hat{b} \, \hat{b}^\dagger </math>: 
:<math> {:\,}\hat{b} \, \hat{b}^\dagger{\,:} = \hat{b}^\dagger \, \hat{b}. </math>
Here the normal ordering operation has ''reordered'' the terms by placing <math>\hat{b}^\dagger</math> to the left of <math>\hat{b}</math>.

These two results can be combined with the commutation relation obeyed by <math>\hat{b}</math> and <math>\hat{b}^\dagger</math> to get

:<math> \hat{b} \, \hat{b}^\dagger = \hat{b}^\dagger \, \hat{b} + 1 = {:\,}\hat{b} \, \hat{b}^\dagger{\,:} \; + 1.</math>
or
:<math> \hat{b} \, \hat{b}^\dagger -  {:\,}\hat{b} \, \hat{b}^\dagger{\,:} = 1.</math>

This equation is used in defining the contractions used in [[Wick's theorem]].

3. An example with multiple operators is:

:<math> {:\,}\hat{b}^\dagger \, \hat{b} \, \hat{b} \, \hat{b}^\dagger \, \hat{b} \, \hat{b}^\dagger \, \hat{b}{\,:} = \hat{b}^\dagger \, \hat{b}^\dagger \, \hat{b}^\dagger \, \hat{b} \, \hat{b} \, \hat{b} \, \hat{b} = (\hat{b}^\dagger)^3 \, \hat{b}^4.</math>

4. A simple example shows that normal ordering cannot be extended by linearity from the monomials to all operators in a self-consistent way. Assume that we can apply the commutation relations to obtain:

:<math> {:\,}\hat{b} \hat{b}^\dagger{\,:} = {:\,}1 + \hat{b}^\dagger \hat{b}{\,:}.</math>

Then, by linearity,

: <math>{:\,}1 + \hat{b}^\dagger \hat{b}{\,:} = {:\,}1{\,:} + {:\,}\hat{b}^\dagger \hat{b}{\,:} = 
1 + \hat{b}^\dagger \hat{b} \ne \hat{b}^\dagger \hat{b}={:\,}\hat{b} \hat{b}^\dagger{\,:},</math>

a contradiction.

The implication is that normal ordering is not a linear function on operators, but on the [[free algebra]] generated by the operators, i.e. the operators do not satisfy the [[canonical commutation relations]] while inside the normal ordering (or any other ordering operator like [[time-ordering]], etc).

===Multiple bosons===
If we now consider <math>N</math> different bosons there are <math>2 N</math> operators:
* <math>\hat{b}_i^\dagger</math>: the <math>i^{th}</math> boson's creation operator.
* <math>\hat{b}_i</math>: the <math>i^{th}</math> boson's annihilation operator.
Here <math>i = 1,\ldots,N</math>.

These satisfy the commutation relations:
:<math>\left[\hat{b}_i^\dagger, \hat{b}_j^\dagger \right]_- = 0 </math>
:<math>\left[\hat{b}_i, \hat{b}_j \right]_- = 0 </math>
:<math>\left[\hat{b}_i, \hat{b}_j^\dagger \right]_- = \delta_{ij} </math>
where <math>i,j = 1,\ldots,N</math> and <math>\delta_{ij}</math> denotes the [[Kronecker delta]].

These may be rewritten as: 
:<math>\hat{b}_i^\dagger \, \hat{b}_j^\dagger = \hat{b}_j^\dagger \, \hat{b}_i^\dagger </math>
:<math>\hat{b}_i \, \hat{b}_j = \hat{b}_j \, \hat{b}_i </math>
:<math>\hat{b}_i \,\hat{b}_j^\dagger = \hat{b}_j^\dagger \,\hat{b}_i + \delta_{ij}.</math>

====Examples====
1. For two different bosons (<math>N=2</math>) we have
:<math> : \hat{b}_1^\dagger \,\hat{b}_2 : \,= \hat{b}_1^\dagger \,\hat{b}_2 </math>
:<math> : \hat{b}_2 \, \hat{b}_1^\dagger  : \,= \hat{b}_1^\dagger \,\hat{b}_2 </math>

2. For three different bosons (<math>N=3</math>) we have
:<math> : \hat{b}_1^\dagger \,\hat{b}_2 \,\hat{b}_3 : \,= \hat{b}_1^\dagger \,\hat{b}_2 \,\hat{b}_3</math>
Notice that since (by the commutation relations) <math>\hat{b}_2 \,\hat{b}_3 = \hat{b}_3 \,\hat{b}_2</math> the order in which we write the annihilation operators does not matter.

:<math> : \hat{b}_2 \, \hat{b}_1^\dagger \, \hat{b}_3  : \,= \hat{b}_1^\dagger \,\hat{b}_2 \, \hat{b}_3 </math>
:<math> : \hat{b}_3 \hat{b}_2 \, \hat{b}_1^\dagger  : \,= \hat{b}_1^\dagger \,\hat{b}_2 \, \hat{b}_3 </math>

===Bosonic operator functions===
Normal ordering of bosonic operator functions <math>f(\hat n)</math>, with occupation number operator <math>\hat n=\hat b\vphantom{\hat n}^\dagger \hat b</math>, can be accomplished using [[Factorial power|(falling) factorial powers]] <math>\hat n^{\underline{k}}=\hat n(\hat n-1)\cdots(\hat n-k+1)</math> and [[Newton series]] instead of [[Taylor series]]: 
It is easy to show 
<ref name="Hucht">{{cite journal | last1=König | first1=Jürgen | last2=Hucht | first2=Alfred | title=Newton series expansion of bosonic operator functions | journal=SciPost Physics | publisher=Stichting SciPost | volume=10 | issue=1 | date=2021-01-13 | issn=2542-4653 | doi=10.21468/scipostphys.10.1.007 | page=007| arxiv=2008.11139 | bibcode=2021ScPP...10....7K | s2cid=221293056 | doi-access=free }}</ref> 
that factorial powers <math>\hat n^{\underline{k}}</math> are equal to normal-ordered (raw) [[Exponentiation|powers]] <math>\hat n^{k}</math> and are therefore normal ordered by construction,

: <math>
\hat{n}^{\underline{k}}
 = \hat b\vphantom{\hat n}^{\dagger k} \hat b\vphantom{\hat n}^k
 = {:\,}\hat n^k{\,:},
</math>

such that the Newton series expansion 

: <math>
\tilde f(\hat n) = \sum_{k=0}^\infty \Delta_n^k \tilde f(0) \, \frac{\hat n^{\underline{k}}}{k!}
</math>

of an operator function <math>\tilde f(\hat n)</math>, with <math>k</math>-th [[forward difference]] <math>\Delta_n^k \tilde f(0)</math> at <math>n=0</math>, is always normal ordered. Here, the [[Second_quantization#Action_on_Fock_states|eigenvalue equation]] <math>\hat n |n\rangle = n |n\rangle</math> relates <math>\hat n</math> and <math>n</math>.

As a consequence, the normal-ordered Taylor series of an arbitrary function <math>f(\hat n)</math> is equal to the Newton series of an associated function <math>\tilde f(\hat n)</math>, fulfilling

: <math> 
\tilde f(\hat n) = {:\,} f(\hat n) {\,:},
</math>

if the series coefficients of the Taylor series of <math>f(x)</math>, with continuous <math>x</math>, match the coefficients of the Newton series of <math>\tilde f(n)</math>, with integer <math>n</math>,

: <math>
\begin{align}
       f(x) &= \sum_{k=0}^\infty F_k \, \frac{x^k              }{k!}, \\
\tilde f(n) &= \sum_{k=0}^\infty F_k \, \frac{n^{\underline{k}}}{k!}, \\
F_k &= \partial_x^k f(0) = \Delta_n^k \tilde f(0),
\end{align}
</math>

with <math>k</math>-th [[partial derivative]] <math>\partial_x^k f(0)</math> at <math>x=0</math>.
The functions <math>f</math> and <math>\tilde f</math> are related through the so-called '''[[normal-order transform]]''' <math>\mathcal N[f]</math> according to

: <math>
\begin{align}
\tilde f(n) &= \mathcal N_x[f(x)](n) \\
 &= \frac{1}{\Gamma(-n)} \int_{-\infty}^0 \mathrm d x \, e^x \, f(x) \, (-x)^{-(n+1)} \\
 &= \frac{1}{\Gamma(-n)}\mathcal M_{-x}[e^{x} f(x)](-n),
\end{align}
</math>

which can be expressed in terms of the [[Mellin transform]] <math>\mathcal M</math>, see <ref name="Hucht"/> for details.

==Fermions==

[[Fermions]] are particles which satisfy [[Fermi–Dirac statistics]]. We will now examine the normal ordering of fermionic creation and annihilation operator products.

===Single fermions===
For a single fermion there are two operators of interest:

* <math>\hat{f}^\dagger</math>: the fermion's creation operator.
* <math>\hat{f}</math>: the fermion's annihilation operator.

These satisfy the [[anticommutator]] relationships
:<math>\left[\hat{f}^\dagger, \hat{f}^\dagger \right]_+ = 0</math>
:<math>\left[\hat{f}, \hat{f} \right]_+ = 0</math>
:<math>\left[\hat{f}, \hat{f}^\dagger \right]_+ = 1</math>

where <math>\left[A, B \right]_+ \equiv AB + BA</math> denotes the [[anticommutator]]. These may be rewritten as

:<math>\hat{f}^\dagger\, \hat{f}^\dagger =  0 </math>
:<math>\hat{f} \,\hat{f}  =  0 </math>
:<math>\hat{f} \,\hat{f}^\dagger =  1 - \hat{f}^\dagger \,\hat{f} .</math>

To define the normal ordering of a product of fermionic creation and annihilation operators we must take into account the number of [[Transposition (mathematics)|interchange]]s between neighbouring operators. We get a minus sign for each such interchange.

====Examples====
1. We again start with the simplest cases:
:<math> : \hat{f}^\dagger \, \hat{f} : \,= \hat{f}^\dagger \, \hat{f} </math>
This expression is already in normal order so nothing is changed. In the reverse case, we introduce a minus sign because we have to change the order of two operators:

:<math> : \hat{f} \, \hat{f}^\dagger : \,= -\hat{f}^\dagger \, \hat{f} </math>

These can be combined, along with the anticommutation relations, to show
:<math> \hat{f} \, \hat{f}^\dagger \,= 1 - \hat{f}^\dagger \, \hat{f} = 1 + :\hat{f} \,\hat{f}^\dagger :</math>
or
:<math> \hat{f} \, \hat{f}^\dagger -  : \hat{f} \, \hat{f}^\dagger : = 1.</math>

This equation, which is in the same form as the bosonic case above, is used in defining the contractions used in [[Wick's theorem]].

2. The normal order of any more complicated cases gives zero because there will be at least one creation or annihilation operator appearing twice. For example:
:<math> : \hat{f}\,\hat{f}^\dagger \, \hat{f} \hat{f}^\dagger  : \,= -\hat{f}^\dagger \,\hat{f}^\dagger \,\hat{f}\,\hat{f} = 0 </math>

===Multiple fermions===
For <math>N</math> different fermions there are <math>2 N</math> operators:
* <math>\hat{f}_i^\dagger</math>: the <math>i^{th}</math> fermion's creation operator.
* <math>\hat{f}_i</math>: the <math>i^{th}</math> fermion's annihilation operator.
Here <math>i = 1,\ldots,N</math>.

These satisfy the anti-commutation relations:
:<math>\left[\hat{f}_i^\dagger, \hat{f}_j^\dagger \right]_+ = 0 </math>
:<math>\left[\hat{f}_i, \hat{f}_j \right]_+ = 0 </math>
:<math>\left[\hat{f}_i, \hat{f}_j^\dagger \right]_+ = \delta_{ij} </math>
where <math>i,j = 1,\ldots,N</math> and <math>\delta_{ij}</math> denotes the [[Kronecker delta]].

These may be rewritten as:
:<math>\hat{f}_i^\dagger \, \hat{f}_j^\dagger = -\hat{f}_j^\dagger \, \hat{f}_i^\dagger </math>
:<math>\hat{f}_i \, \hat{f}_j = -\hat{f}_j \, \hat{f}_i </math>
:<math>\hat{f}_i \,\hat{f}_j^\dagger = \delta_{ij} - \hat{f}_j^\dagger \,\hat{f}_i .</math>

When calculating the normal order of products of fermion operators we must take into account the number of [[Transposition (mathematics)|interchange]]s of neighbouring operators required to rearrange the expression. It is as if we pretend the creation and annihilation operators anticommute and then we reorder the expression to ensure the creation operators are on the left and the annihilation operators are on the right - all the time taking account of the anticommutation relations.

====Examples====
1. For two different fermions (<math>N=2</math>) we have
:<math> : \hat{f}_1^\dagger \,\hat{f}_2 : \,= \hat{f}_1^\dagger \,\hat{f}_2 </math>
Here the expression is already normal ordered so nothing changes.

:<math> : \hat{f}_2 \, \hat{f}_1^\dagger  : \,= -\hat{f}_1^\dagger \,\hat{f}_2 </math>
Here we introduce a minus sign because we have interchanged the order of two operators.

:<math> : \hat{f}_2 \, \hat{f}_1^\dagger \, \hat{f}^\dagger_2  : \,= \hat{f}_1^\dagger \, \hat{f}_2^\dagger \,\hat{f}_2 = -\hat{f}_2^\dagger \, \hat{f}_1^\dagger \,\hat{f}_2  </math>
Note that the order in which we write the operators here, unlike in the bosonic case, ''does matter''.

2. For three different fermions (<math>N=3</math>) we have
:<math> : \hat{f}_1^\dagger \, \hat{f}_2 \, \hat{f}_3 : \,= \hat{f}_1^\dagger \,\hat{f}_2 \,\hat{f}_3 = -\hat{f}_1^\dagger \,\hat{f}_3 \,\hat{f}_2</math>
Notice that since (by the anticommutation relations) <math>\hat{f}_2 \,\hat{f}_3 = -\hat{f}_3 \,\hat{f}_2</math> the order in which we write the operators ''does matter'' in this case.

Similarly we have
:<math> : \hat{f}_2 \, \hat{f}_1^\dagger \, \hat{f}_3  : \,= -\hat{f}_1^\dagger \,\hat{f}_2 \, \hat{f}_3 = \hat{f}_1^\dagger \,\hat{f}_3 \, \hat{f}_2</math>
:<math> : \hat{f}_3 \hat{f}_2 \, \hat{f}_1^\dagger  : \,= \hat{f}_1^\dagger \,\hat{f}_3 \, \hat{f}_2 = -\hat{f}_1^\dagger \,\hat{f}_2 \, \hat{f}_3 </math>

== Uses in quantum field theory==

The [[vacuum expectation value]] of a normal ordered product of creation and annihilation operators is zero. This is because, denoting the [[vacuum state]] by <math>|0\rangle</math>, the creation and annihilation operators satisfy
:<math>\langle 0 | \hat{a}^\dagger = 0 \qquad \textrm{and} \qquad \hat{a} |0\rangle = 0</math>
(here <math>\hat{a}^\dagger</math> and <math>\hat{a}</math> are creation and annihilation operators (either bosonic or fermionic)).

Let <math>\hat{O}</math> denote a non-empty product of creation and annihilation operators. Although this may satisfy
:<math>\langle 0 | \hat{O} | 0 \rangle \neq 0,</math>
we have
:<math>\langle 0 | :\hat{O}: | 0 \rangle = 0.</math>

Normal ordered operators are particularly useful when defining a quantum mechanical [[Hamiltonian (quantum mechanics)|Hamiltonian]]. If the Hamiltonian of a theory is in normal order then the ground state energy will be zero:
<math>\langle 0 |\hat{H}|0\rangle = 0</math>.

===Free fields===
With two free fields φ and χ,

:<math>:\phi(x)\chi(y):\,\,=\phi(x)\chi(y)-\langle 0|\phi(x)\chi(y)| 0\rangle</math>

where <math>|0\rangle</math> is again the vacuum state. Each of the two terms on the right hand side typically blows up in the limit as y approaches x but the difference between them has a well-defined limit. This allows us to define :φ(x)χ(x):.

===Wick's theorem===
{{Main|Wick's theorem}}

'''Wick's theorem''' states the relationship between the time ordered product of <math>n</math> fields and a sum of 
normal ordered products. This may be expressed for <math>n</math> even as

:<math>\begin{align}
T\left[\phi(x_1)\cdots \phi(x_n)\right]=&:\phi(x_1)\cdots \phi(x_n):
+\sum_\textrm{perm}\langle 0 |T\left[\phi(x_1)\phi(x_2)\right]|0\rangle :\phi(x_3)\cdots \phi(x_n):\\
&+\sum_\textrm{perm}\langle 0 |T\left[\phi(x_1)\phi(x_2)\right]|0\rangle \langle 0 |T\left[\phi(x_3)\phi(x_4)\right]|0\rangle:\phi(x_5)\cdots \phi(x_n):\\ 
\vdots \\
&+\sum_\textrm{perm}\langle 0 |T\left[\phi(x_1)\phi(x_2)\right]|0\rangle\cdots \langle 0 |T\left[\phi(x_{n-1})\phi(x_n)\right]|0\rangle
\end{align}</math>

where the summation is over all the distinct ways in which one may pair up fields. The result for <math>n</math> odd looks the same
except for the last line which reads

:<math>
\sum_\text{perm}\langle 0 |T\left[\phi(x_1)\phi(x_2)\right]|0\rangle\cdots\langle 0 | T\left[\phi(x_{n-2})\phi(x_{n-1})\right]|0\rangle\phi(x_n).
</math>

This theorem provides a simple method for computing vacuum expectation values of time ordered products of operators and was the motivation behind the introduction of normal ordering.

==Alternative definitions==

The most general definition of normal ordering involves splitting all quantum fields into two parts (for example see Evans and Steer 1996)
<math>\phi_i(x)=\phi^+_i(x)+\phi^-_i(x)</math>.  
In a product of fields, the fields are split into the two parts and the <math>\phi^+(x)</math> parts are moved so as to be always to the left of all the <math>\phi^-(x)</math> parts.  In the usual case considered in the rest of the article, the <math>\phi^+(x)</math> contains only creation operators, while the <math>\phi^-(x)</math> contains only annihilation operators.  As this is a mathematical identity, one can split fields in any way one likes.  However, for this to be a useful procedure one demands that the normal ordered product of ''any'' combination of fields has zero expectation value

:<math>\langle :\phi_1(x_1)\phi_2(x_2)\ldots\phi_n(x_n):\rangle=0</math>

It is also important for practical calculations that all the commutators (anti-commutator for fermionic fields) of all <math>\phi^+_i</math> and <math>\phi^-_j</math> are all c-numbers.  These two properties means that we can apply '''Wick's theorem''' in the usual way, turning expectation values of time-ordered products of fields into products of c-number pairs, the contractions.  In this generalised setting, the contraction is defined to be the difference between the time-ordered product and the normal ordered product of a pair of fields.

The simplest example is found in the context of [[thermal quantum field theory]] (Evans and Steer 1996).  In this case the expectation values of interest are statistical ensembles, traces over all states weighted by  <math>\exp (-\beta \hat{H})</math>.  For instance, for a single bosonic quantum harmonic oscillator we have that the thermal expectation value of the number operator is simply the [[Bose–Einstein statistics|Bose–Einstein distribution]]

:<math>\langle\hat{b}^\dagger \hat{b}\rangle
 = \frac{\mathrm{Tr} (e^{-\beta \omega \hat{b}^\dagger \hat{b}} \hat{b}^\dagger \hat{b} )}{\mathrm{Tr} (e^{-\beta \omega \hat{b}^\dagger \hat{b} })}
 = \frac{1}{e^{\beta \omega}-1}
</math>

So here the number operator <math>\hat{b}^\dagger \hat{b}</math> is normal ordered in the usual sense used in the rest of the article yet its thermal expectation values are non-zero.  Applying Wick's theorem and doing calculation with the usual normal ordering in this thermal context is possible but computationally impractical. The solution is to define a different ordering, such that the <math>\phi^+_i</math> and <math>\phi^-_j</math> are ''linear combinations'' of the original annihilation and creations operators.  The combinations are chosen to ensure that the thermal expectation values of normal ordered products are always zero so the split chosen will depend on the temperature.

==References==
{{Reflist}}

* F. Mandl, G. Shaw, Quantum Field Theory, John Wiley & Sons, 1984.
* S. Weinberg, The Quantum Theory of Fields (Volume I) Cambridge University Press (1995)
* T.S. Evans, D.A. Steer, [https://arxiv.org/abs/hep-ph/9601268 Wick's theorem at finite temperature], Nucl. Phys B 474, 481-496 (1996) [https://arxiv.org/abs/hep-ph/9601268 arXiv:hep-ph/9601268]

{{DEFAULTSORT:Normal Order}}
[[Category:Quantum field theory]]