# Photon statistics

Photon statistics is the theoretical and experimental study of the statistical distributions produced in photon counting experiments, which use photodetectors to analyze the intrinsic statistical nature of photons in a light source. In these experiments, light incident on the photodetector generates photoelectrons and a counter registers electrical pulses generating a statistical distribution of photon counts. Low intensity disparate light sources can be differentiated by the corresponding statistical distributions produced in the detection process.
Three regimes of statistical distributions can be obtained depending on the properties of the light source: Poissonian, super-Poissonian, and sub-Poissonian. The regimes are defined by the relationship between the variance and average number of photon counts for the corresponding distribution. Both Poissonian and super-Poissonian light can be described by a semi-classical theory in which the light source is modeled as an electromagnetic wave and the atom is modeled according to quantum mechanics. In contrast, sub-Poissonian light requires the quantization of the electromagnetic field for a proper description and thus is a direct measure of the particle nature of light.

## Poissonian light

In classical electromagnetic theory, an ideal source of light with constant intensity can be modeled by a spatially and temporally coherent electromagnetic wave of a single frequency. Such a light source can be modeled by,
    $E(x,t)=E_{0}\sin(kx-\omega t+\phi )$
where $\omega$ is the frequency of the field and $\phi$ is a time independent phase shift.
The analogue in quantum mechanics is the coherent state
    $|\alpha \rangle =\sum _{n=0}^{\infty }{\frac {{\alpha }^{n}}{\sqrt {n!}}}e^{\frac {{-\left\vert {\alpha }\right\vert }^{2}}{2}}|n\rangle$
By projecting the coherent state onto the Fock state $|n\rangle$, we can find the probability $P_{n}$ of finding $n$ photons using the Born rule, which gives
    $P_{n}={\frac {{\left\vert \alpha \right\vert }^{2n}}{n!}}e^{{-\left\vert \alpha \right\vert }^{2}}={\frac {{\langle n\rangle }^{n}}{n!}}e^{-\langle n\rangle }$
The above result is a Poissonian distribution with variance ${\Delta n}^{2}=\langle n\rangle$ which is a distinct feature of the coherent state.

## Super-Poissonian light

Light that is governed by super-Poissonian statistics exhibits a statistical distribution with variance ${\Delta n}^{2}>\langle n\rangle$. An example of light that exhibits super-Poissonian statistics is thermal light. The intensity of thermal light fluctuates randomly and the fluctuations give rise to super-Poissonian statistics, as shown below by calculating the distribution of the intensity fluctuations. Using the intensity distribution together with Mandel's formula which describes the probability of the number of photon counts registered by a photodetector, the statistical distribution of photons in thermal light can be obtained.
Thermal light can be modeled as a collection of $N$ harmonic oscillators. Suppose the $j$-th oscillator emits an electromagnetic field $E_{j}(t)=E_{0}e^{-i\omega t}e^{i\phi _{j}(t)}$ with phase $\phi _{j}(t)$. Using the theory of superposition of fields, the total field produced by the $N$ oscillators is
    $E(t)=\sum _{j=1}^{N}E_{j}(t)=\sum _{j=1}^{N}E_{0}e^{-i\omega t}e^{i\phi _{j}(t)}$
After pulling out all the variables that are independent of the summation index $j$, a random complex amplitude can be defined by
    $\beta (t)=\sum _{j=1}^{N}e^{i\phi _{j}(t)}=b(t)e^{i\Phi (t)}$
where $\beta (t)$ was rewritten in terms of its magnitude $b(t)=\left\vert \beta \right\vert$ and its phase $e^{i\Phi (t)}$. Because the oscillators are uncorrelated, the phase of the superposed field will be random. Therefore, the complex amplitude $\beta (t)$ is a stochastic variable. It represents the sum of the uncorrelated phases of the oscillators which models the intensity fluctuations in thermal light. On the complex plane, it represents a two dimensional random walker with $N$ representing the steps taken. For large $N$ a random walker has a Gaussian probability distribution. Thus, the joint probability distribution for the real and imaginary parts of the complex random variable $\beta (t)$ can be represented as,
    $p(\beta (t))=\left[{\frac {1}{\sqrt {2\pi \sigma ^{2}}}}e^{\frac {-{Re(\beta (t))}^{2}}{2\sigma ^{2}}}\right]\left[{\frac {1}{\sqrt {2\pi \sigma ^{2}}}}e^{\frac {-{Im(\beta (t))}^{2}}{2\sigma ^{2}}}\right]$
    $={\frac {1}{2\pi \sigma ^{2}}}e^{\frac {-\left({Re(\beta (t))}^{2}+{Im(\beta (t))}^{2}\right)}{2\sigma ^{2}}}={\frac {1}{2\pi \sigma ^{2}}}e^{\frac {-{\left\vert \beta (t)\right\vert }^{2}}{2\sigma ^{2}}}$
After $N$ steps, the expectation value of the radius squared is $\langle {\left\vert \beta (t)\right\vert }^{2}\rangle =2\sigma ^{2}=N$. The expectation value $\langle \left\vert \beta (t)\right\vert \rangle =0$ which can be thought of as all directions being equally likely. Rewriting the probability distribution in terms of $N$ results in
    $p(\beta (t))={\frac {1}{\pi N}}e^{\frac {-{\left\vert \beta (t)\right\vert }^{2}}{N}}$
With the probability distribution above, we can now find the average intensity of the field (where several constants have been omitted for clarity)
    $\langle I(t)\rangle =\langle E^{*}(t)E(t)\rangle$
    $={E_{0}}^{2}{\langle \left\vert \beta (t)\right\vert }^{2}\rangle =N{E_{0}}^{2}$
The instantaneous intensity of the field $I(t)$ is given by
    $I(t)=E^{*}(t)E(t)=E_{0}^{2}{\left\vert \beta (t)\right\vert }^{2}$
Because the electric field and thus the intensity are dependent on the stochastic complex variable $\beta (t)$. The probability of obtaining an intensity in between $I$ and $I+dI$ is
    $p(I(t))dI=p(I(\beta (t)))d^{2}\beta$
where $d^{2}\beta$ is the infinitesimal element on the complex plane. This infinitesimal element can be rewritten as
    $d^{2}\beta =2\pi \left\vert \beta (t)\right\vert d\left\vert \beta (t)\right\vert$
The above intensity distribution can now be written as
    $p(I(t))=p(\left\vert \beta (t)\right\vert )2\pi \left\vert \beta (t)\right\vert {\left({\frac {dI}{d\left\vert \beta (t)\right\vert }}\right)}^{-1}$
    $={\frac {1}{\pi N}}e^{\frac {-{\left\vert \beta (t)\right\vert }^{2}}{N}}2\pi \left\vert \beta (t)\right\vert {\left(2E_{0}^{2}\left\vert \beta (t)\right\vert \right)}^{-1}$
    $={\frac {1}{NE_{0}^{2}}}e^{\frac {-E_{0}^{2}{\left\vert \beta (t)\right\vert }^{2}}{NE_{0}^{2}}}={\frac {1}{\langle I(t)\rangle }}e^{\frac {-I(t)}{\langle I(t)\rangle }}$
This last expression represents the intensity distribution for thermal light. The last step in showing thermal light satisfies the variance condition for super-Poisson statistics is to use Mandel's formula. The formula describes the probability of observing n photon counts and is given by
    $P(n)=\int _{0}^{\infty }{\frac {{\left(\epsilon I\right)}^{n}}{n!}}e^{-\epsilon I}P(I)dI$
The factor $\epsilon ={\frac {\eta }{h\nu }}$ where $\eta$ is the quantum efficiency describes the efficiency of the photon counter. A perfect detector would have $\eta =1$. $I$ is the intensity incident on an area A of the photodetector and is given by
    Comparison of the Poisson and Bose-Einstein distributions. The Poisson distribution is characteristic of coherent light while the Bose-Einstein distribution is characteristic of thermal light. Both distribution have the same expectation value $\langle n\rangle =6$.$I=\iint \limits _{A}\int _{t}^{t+\tau }\,I(x,y,t')dt'\,dx\,dy$
On substituting the intensity probability distribution of thermal light for P(I), Mandel's formula becomes
    $P(n)=\int _{0}^{\infty }{\frac {{\left(\epsilon I\right)}^{n}}{n!}}e^{-\epsilon I}{\frac {e^{\frac {-I}{\langle I\rangle }}}{\langle I\rangle }}dI={\frac {{\epsilon }^{n}}{n!\langle I\rangle }}\int _{0}^{\infty }{I}^{n}e^{-\left(\epsilon +{\frac {1}{\langle I\rangle }}\right)I}dI$
Using the following formula to evaluate the integral
$\int _{0}^{\infty }x^{n}e^{-ax}dx={\frac {n!}{a^{n+1}}}\quad \left(n=0,1,2,..a>0\right)$
The probability distribution for n photon counts from a thermal light source is
    $P(n)={\frac {1}{\left(\langle n\rangle +1\right)}}{\left({\frac {\langle n\rangle }{\langle n\rangle +1}}\right)}^{n}$
where $\langle n\rangle =\epsilon I$ is the average number of counts. This last distribution is known as the Bose-Einstein distribution. The variance of the distribution can be shown to be
    ${\Delta n}^{2}=\langle n\rangle +{\langle n\rangle }^{2}$
In contrast with the Poisson distribution for a coherent light source, the Bose-Einstein distribution has ${\Delta n}^{2}>\langle n\rangle$ characteristic of thermal light.

## Sub-Poissonian light

Light that is governed by sub-Poisson statistics cannot be described by classical electromagnetic theory and is defined by ${\Delta n}^{2}<\langle n\rangle$. The advent of ultrafast photodetectors has made it possible to measure the sub-Poissonian nature of light. An example of light exhibiting sub-Poissonian statistics is squeezed light. Recently researchers have shown that sub-Poissonian light can be induced in a quantum dot exhibiting resonance fluorescence. A technique used to measure the sub-Poissonian structure of light is a homodyne intensity correlation scheme. In this scheme a local oscillator and signal field are superimposed via a beam splitter. The superimposed light is then split by another beam splitter and each signal is recorded by individual photodetectors connected to correlator from which the intensity correlation can be measured. Evidence of the sub-Poissonian nature of light is shown by obtaining a negative intensity correlation as was shown in.