# Optical phase space

Optical phase diagram of a coherent state's distribution across phase space.

In quantum optics, an optical phase space is a phase space in which all quantum states of an optical system are described. Each point in the optical phase space corresponds to a unique state of an optical system. For any such system, a plot of the quadratures against each other, possibly as functions of time, is called a phase diagram. If the quadratures are functions of time then the optical phase diagram can show the evolution of a quantum optical system with time.

An optical phase diagram can give insight into the properties and behaviors of the system that might otherwise not be obvious. This can allude to qualities of the system that can be of interest to an individual studying an optical system that would be very hard to deduce otherwise. Another use for an optical phase diagram is that it shows the evolution of the state of an optical system. This can be used to determine the state of the optical system at any point in time.

## Background information

When discussing the quantum theory of light, it is very common to use an electromagnetic oscillator as a model.[1] An electromagnetic oscillator describes an oscillation of the electric field. Since the magnetic field is proportional to the rate of change of the electric field, this too oscillates. Such oscillations describe light. Systems composed of such oscillators can be described by an optical phase space.

Let u(x,t) be a vector function describing a single mode of an electromagnetic oscillator. For simplicitity, it is assumed that this electromagnetic oscillator is in vacuum. An example is the plane wave given by

${\displaystyle \mathbf {u} (\mathbf {x} ,t)=\mathbf {u_{0}} e^{i(\mathbf {k} \cdot \mathbf {x} -wt)}}$

where u0 is the polarization vector, k is the wave vector, w the frequency, and A${\displaystyle \cdot }$B denotes the dot product between the vectors A and B. This is the equation for a plane wave and is a simple example of such an electromagnetic oscillator. The oscillators being examined could either be free waves in space or some normal mode contained in some cavity.

A single mode of the electromagnetic oscillator is isolated from the rest of the system and examined. Such an oscillator, when quantized, is described by the mathematics of a quantum harmonic oscillator.[1] Quantum oscillators are described using creation and annihilation operators ${\displaystyle {\widehat {a}}^{\dagger }}$ and ${\displaystyle {\widehat {a}}}$. Physical quantities, such as the electric field strength, then become quantum operators.

In order to distinguish a physical quantity from the quantum mechanical operator used to describe it, a "hat" is used over the operator symbols. Thus, for example, where ${\displaystyle E_{i}}$ might represent (one component of) the electric field, the symbol ${\displaystyle {\widehat {E}}_{i}}$ denotes the quantum-mechanical operator that describes ${\displaystyle E_{i}}$. This convention is used throughout this article, but is not in common use in more advanced texts, which avoid the hat, as it simply clutters the text.

In the quantum oscillator mode, most operators representing physical quantities are typically expressed in terms of the creation and annihilation operators. In this example, the electric field strength is given by:

${\displaystyle {\widehat {E}}_{i}=u_{i}^{*}(\mathbf {x} ,t){\widehat {a}}^{\dagger }+u_{i}(\mathbf {x} ,t){\widehat {a}}}$[2]

(where xi is a single component of x, position). The Hamiltonian for an electromagnetic oscillator is found by quantizing the electromagnetic field for this oscillator and the formula is given by:

${\displaystyle {\widehat {H}}=\hbar \omega ({\widehat {a}}^{\dagger }{\widehat {a}}+1/2)}$[2]

where ${\displaystyle \omega }$ is the frequency of the (spatio-temportal) mode. The annihilation operator is the bosonic annihilation operator and so it obeys the canonical commutation relation given by:

${\displaystyle [{\widehat {a}},{\widehat {a}}^{\dagger }]=1}$

The eigenstates of the annihilation operator are called coherent states:

${\displaystyle {\widehat {a}}|\alpha \rangle =\alpha |\alpha \rangle }$

It is important to note that the annihilation operator is not Hermitian; therefore its eigenvalues ${\displaystyle \alpha }$ can be complex. This has important consequences.

Finally, the photon number is given by the operator ${\displaystyle {\widehat {N}}={\widehat {a}}^{\dagger }{\widehat {a}},}$ which gives the number of photons in the given (spatial-temporal) mode u.

Operators given by

$\displaystyle \widehat q = \tfrac 1 { 2}(\widehat a^\dagger + \widehat a)$

and

${\displaystyle {\widehat {p}}={\tfrac {i}{2}}({\widehat {a}}^{\dagger }-{\widehat {a}})}$

are called the quadratures and they represent the real and imaginary parts of the complex amplitude represented by ${\displaystyle {\widehat {a}}}$.[1] The commutation relation between the two quadratures can easily be calculated:

{\displaystyle {\begin{aligned}\left[{\widehat {q}},{\widehat {p}}\right]&={\tfrac {i}{4}}[{\widehat {a}}^{\dagger }+{\widehat {a}},{\widehat {a}}^{\dagger }-{\widehat {a}}]\\&={\tfrac {i}{4}}([{\widehat {a}}^{\dagger },{\widehat {a}}^{\dagger }]-[{\widehat {a}}^{\dagger },{\widehat {a}}]+[{\widehat {a}},{\widehat {a}}^{\dagger }]-[{\widehat {a}},{\widehat {a}}])\\&={\tfrac {i}{4}}(-(-1)+1)\\&={\tfrac {i}{2}}\end{aligned}}}

This looks very similar to the commutation relation of the position and momentum operator. Thus, it can be useful to think of and treat the quadratures as the position and momentum of the oscillator although in fact they are the "in-phase and out-of-phase components of the electric field amplitude of the spatial-temporal mode", or u, and have nothing really to do with the position or momentum of the electromagnetic oscillator (as it is hard to define what is meant by position and momentum for an electromagnetic oscillator).[1]

The eigenstates of the quadrature operators ${\displaystyle {\widehat {q}}}$ and ${\displaystyle {\widehat {p}}}$ are called the quadrature states. They satisfy the relations:

• ${\displaystyle {\widehat {q}}|q\rangle =q|q\rangle }$ and ${\displaystyle {\widehat {p}}|p\rangle =p|p\rangle }$
• ${\displaystyle \langle q|q'\rangle =\delta (q-q')}$ and ${\displaystyle \langle p|p'\rangle =\delta (p-p')}$
• ${\displaystyle \int _{-\infty }^{\infty }|q\rangle \langle q|dq=1}$ and ${\displaystyle \int _{-\infty }^{\infty }|p\rangle \langle p|dp=1}$

as these form complete basis sets.

### Important result

The following is an important relation that can be derived from the above which justifies our interpretation that the quadratures are the real and imaginary parts of a complex ${\displaystyle \alpha }$ (i.e. the in-phase and out-of-phase components of the electromagnetic oscillator)

${\displaystyle \langle \alpha |{\widehat {q}}|\alpha \rangle =2^{-1/2}(\langle \alpha |{\widehat {a}}^{\dagger }|\alpha \rangle +\langle \alpha |{\widehat {a}}|\alpha \rangle )=2^{-1/2}(\alpha ^{*}\langle \alpha |\alpha \rangle +\alpha \langle \alpha |\alpha \rangle )}$

The following is a relationship that can be used to help evaluate the above and is given by:

${\displaystyle \langle \alpha '|\alpha \rangle =e^{(-1/2)(|\alpha '|^{2}+|\alpha |^{2})+\alpha '^{*}\alpha }}$[1]

This gives us that:

${\displaystyle \langle \alpha |{\widehat {q}}|\alpha \rangle =2^{-1/2}(\alpha ^{*}+\alpha )=q_{\alpha }}$
${\displaystyle \langle \alpha |{\widehat {p}}|\alpha \rangle =i2^{-1/2}(\alpha ^{*}-\alpha )=p_{\alpha }}$ by a similar method as above.
${\displaystyle \alpha =2^{-1/2}(\langle \alpha |{\widehat {q}}|\alpha \rangle +i\langle \alpha |{\widehat {p}}|\alpha \rangle )=2^{-1/2}(q_{\alpha }+ip_{\alpha })}$

Thus, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "/mathoid/local/v1/":): \alpha is just a composition of the quadratures.

Another very important property of the coherent states becomes very apparent in this formalism. A coherent state is not a point in the optical phase space but rather a distribution on it. This can be seen via

${\displaystyle q_{\alpha }=\langle \alpha |{\widehat {q}}|\alpha \rangle }$

and

${\displaystyle p_{\alpha }=\langle \alpha |{\widehat {p}}|\alpha \rangle }$.

These are only the expectation values of ${\displaystyle {\widehat {q}}}$ and ${\displaystyle {\widehat {p}}}$ for the state ${\displaystyle |\alpha \rangle }$.

It can be shown that the quadratures obey Heisenberg's Uncertainty Principle given by:

${\displaystyle \Delta q\Delta p\geq 1/2}$[1] (where ${\displaystyle \Delta q}$ and ${\displaystyle \Delta p}$ are the variances of the distributions of q and p, respectively)

This inequality does not necessarily have to be saturated and a common example of such states are squeezed coherent states. The coherent states are Gaussian probability distributions over the phase space localized around ${\displaystyle \alpha }$.

## Operators on phase space

It is possible to define operators to move the coherent states around the phase space. These can produce new coherent states and allow us to move around phase space.

### Phase-shifting operator

Phase shifting operator acting on a coherent state rotating it by an angle ${\displaystyle \theta }$ in phase space.

The phase-shifting operator rotates the coherent state by an angle ${\displaystyle \theta }$ in the optical phase space. This operator is given by:

${\displaystyle {\widehat {U}}(\theta )=e^{-i\theta {\widehat {N}}}}$ [1]

The important relationship

${\displaystyle {\widehat {U}}(\theta )^{\dagger }{\widehat {a}}{\widehat {U}}(\theta )={\widehat {a}}e^{-i\theta }}$

is derived as follows:

${\displaystyle d/d\theta ({\widehat {U}}^{\dagger }{\widehat {a}}{\widehat {U}})=i{\widehat {N}}{\widehat {U}}^{\dagger }{\widehat {a}}{\widehat {U}}-i{\widehat {U}}^{\dagger }{\widehat {a}}{\widehat {U}}{\widehat {N}}={\widehat {U}}^{\dagger }i[{\widehat {N}},{\widehat {a}}]{\widehat {U}}}$
${\displaystyle ={\widehat {U}}^{\dagger }i({\widehat {a}}^{\dagger }{\widehat {a}}{\widehat {a}}-{\widehat {a}}{\widehat {a}}^{\dagger }{\widehat {a}}){\widehat {U}}={\widehat {U}}^{\dagger }i[{\widehat {a}}^{\dagger },{\widehat {a}}]{\widehat {a}}{\widehat {U}}=-i{\widehat {U}}^{\dagger }{\widehat {a}}{\widehat {U}}}$

and solving this differential equation yields the desired result.

Thus, using the above it becomes clear that

${\displaystyle {\widehat {U}}(\theta )|\alpha \rangle =|\alpha e^{-i\theta }\rangle }$,

or a rotation by an angle theta on the coherent state in phase space. The following illustrates this more clearly:

${\displaystyle {\widehat {a}}({\widehat {U}}|\alpha \rangle )={\widehat {U}}{\widehat {a}}e^{-i\theta }|\alpha \rangle }$

(which is obtained using the fact that the phase-shifting operator is unitary

${\displaystyle {\widehat {a}}({\widehat {U}}|\alpha \rangle )={\widehat {U}}\alpha e^{-i\theta }|\alpha \rangle =\alpha e^{-i\theta }({\widehat {U}}|\alpha \rangle )}$

Thus,

${\displaystyle (\alpha e^{-i\theta },{\widehat {U}}|\alpha \rangle )}$

is the eigenpair of

${\displaystyle {\widehat {a}}{\widehat {U}}|\alpha \rangle }$.

From this it is possible to see that

${\displaystyle (\alpha e^{-i\theta }=2^{-1/2}[q_{\alpha }cos(\theta )+p_{\alpha }sin(\theta )]+i2^{-1/2}[-q_{\alpha }csin(\theta )+p_{\alpha }cos(\theta )],{\widehat {U}}|\alpha \rangle =|\alpha e^{-i\theta }\rangle )}$

which is another way of expressing the eigenpair which more clearly illustrates the effects of the phase-shifting operator on coherent states.

### Displacement operator

Displacement operator acting on a coherent state displacing it by some value Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "/mathoid/local/v1/":): \alpha in phase space.

The displacement operator takes a coherent state and moves it (by some value) to another coherent state somewhere in the phase space. The displacement operator is given by:

${\displaystyle {\widehat {D}}(\alpha )=e^{\alpha {\widehat {a}}^{\dagger }-\alpha ^{*}{\widehat {a}}}}$

The relationship

${\displaystyle {\widehat {D}}^{\dagger }{\widehat {a}}{\widehat {D}}={\widehat {a}}+\alpha }$.[1]

can be derived quite easily.

To do this, take an infinitesimal displacement ${\displaystyle \delta \alpha }$.

The operators ${\displaystyle {\widehat {D}}}$ and ${\displaystyle {\widehat {D}}^{\dagger }}$ can be expanded using

${\displaystyle e^{X}=\sum _{k=0}^{\infty }{\frac {X^{k}}{k!}}}$

and look at the first order terms and ignore all the higher order terms (all higher order terms are very close to zero as ${\displaystyle \delta \alpha }$ is very small).

Thus:

${\displaystyle {\widehat {D}}^{\dagger }(\delta \alpha ){\widehat {a}}{\widehat {D}}(\delta \alpha )=\sum _{i,j}(\delta \alpha ^{*}{\widehat {a}}-\delta \alpha {\widehat {a}}^{\dagger })^{i}{\widehat {a}}(\delta \alpha {\widehat {a}}^{\dagger }-\delta \alpha ^{*}{\widehat {a}})^{j}/i!j!}$
${\displaystyle ={\widehat {a}}+(\delta \alpha ^{*}{\widehat {a}}-\delta \alpha {\widehat {a}}^{\dagger }){\widehat {a}}+{\widehat {a}}(\delta \alpha {\widehat {a}}^{\dagger }-\delta \alpha ^{*}{\widehat {a}})+O(\delta \alpha ^{2},(\delta \alpha ^{*})^{2})}$ (but as given above, the higher order terms are very close to zero and therefore neglected)
${\displaystyle ={\widehat {a}}+{\widehat {a}}(\delta \alpha {\widehat {a}}^{\dagger }-\delta \alpha ^{*}{\widehat {a}})-(\delta \alpha {\widehat {a}}^{\dagger }-\delta \alpha ^{*}{\widehat {a}}){\widehat {a}}}$
${\displaystyle ={\widehat {a}}+[{\widehat {a}},\delta \alpha {\widehat {a}}^{\dagger }-\delta \alpha ^{*}{\widehat {a}}]={\widehat {a}}+\delta \alpha [{\widehat {a}},{\widehat {a}}^{\dagger }]-\delta \alpha ^{*}[{\widehat {a}},{\widehat {a}}]}$ (use the identity from above)
${\displaystyle {\widehat {D}}^{\dagger }(\delta \alpha ){\widehat {a}}{\widehat {D}}(\delta \alpha )={\widehat {a}}+\delta \alpha }$

The above identity can be applied repeatedly in the following fashion to derive the following:

${\displaystyle ({\widehat {D}}^{\dagger }(\delta \alpha ))^{k}{\widehat {a}}({\widehat {D}}(\delta \alpha ))^{k}={\widehat {a}}+k\delta \alpha }$

Thus, the above identity suggests that repeated use of the displacement operator generates translations in phase space.

#### Important consequence

The following is an important consequence of the displacement vector.

First note that the displacement operator is a unitary operator. Use

${\displaystyle {\widehat {D}}^{\dagger }(\delta \alpha ){\widehat {a}}{\widehat {D}}(\delta \alpha )={\widehat {a}}+\delta \alpha }$

to get:

${\displaystyle {\widehat {a}}{\widehat {D}}(-\alpha )|\alpha \rangle ={\widehat {D}}(-\alpha )({\widehat {a}}-\alpha )|\alpha \rangle }$
${\displaystyle ={\widehat {D}}(-\alpha )({\widehat {a}}|\alpha \rangle -\alpha |\alpha \rangle )}$
${\displaystyle ={\widehat {D}}(-\alpha )(\alpha |\alpha \rangle -\alpha |\alpha \rangle )}$

Thus,

${\displaystyle {\widehat {a}}({\widehat {D}}(-\alpha )|\alpha \rangle )=0}$

or it follows that

${\displaystyle {\widehat {D}}(-\alpha )|\alpha \rangle =|0\rangle }$

${\displaystyle |\alpha \rangle ={\widehat {D}}(\alpha )|0\rangle }$.