Photon statistics

Last updated July 23, 2022

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.^[1] 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,^[1]

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 ^[1]

|\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 ${\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.^[2] Using the intensity distribution together with Mandel's formula ^[3] 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.^[3] 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^[4]

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
<
n
>
=
6
{\displaystyle \langle n\rangle =6}
. PoissonAndBEDistribution.jpg — 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$ .^[1] 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.^[5] A technique used to measure the sub-Poissonian structure of light is a homodyne intensity correlation scheme.^[6] 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.^[5]

Related Research Articles

In quantum mechanics, the uncertainty principle is any of a variety of mathematical inequalities asserting a fundamental limit to the accuracy with which the values for certain pairs of physical quantities of a particle, such as position, x, and momentum, p, can be predicted from initial conditions.

In quantum mechanics, a density matrix is a matrix that describes the quantum state of a physical system. It allows for the calculation of the probabilities of the outcomes of any measurement performed upon this system, using the Born rule. It is a generalization of the more usual state vectors or wavefunctions: while those can only represent pure states, density matrices can also represent mixed states. Mixed states arise in quantum mechanics in two different situations: first when the preparation of the system is not fully known, and thus one must deal with a statistical ensemble of possible preparations, and second when one wants to describe a physical system which is entangled with another, as its state can not be described by a pure state.

In thermodynamics, the Helmholtz free energy is a thermodynamic potential that measures the useful work obtainable from a closed thermodynamic system at a constant temperature (isothermal). The change in the Helmholtz energy during a process is equal to the maximum amount of work that the system can perform in a thermodynamic process in which temperature is held constant. At constant temperature, the Helmholtz free energy is minimized at equilibrium.

In physics, specifically in quantum mechanics, a coherent state is the specific quantum state of the quantum harmonic oscillator, often described as a state which has dynamics most closely resembling the oscillatory behavior of a classical harmonic oscillator. It was the first example of quantum dynamics when Erwin Schrödinger derived it in 1926, while searching for solutions of the Schrödinger equation that satisfy the correspondence principle. The quantum harmonic oscillator arise in the quantum theory of a wide range of physical systems. For instance, a coherent state describes the oscillating motion of a particle confined in a quadratic potential well. The coherent state describes a state in a system for which the ground-state wavepacket is displaced from the origin of the system. This state can be related to classical solutions by a particle oscillating with an amplitude equivalent to the displacement.

In physics, the S-matrix or scattering matrix relates the initial state and the final state of a physical system undergoing a scattering process. It is used in quantum mechanics, scattering theory and quantum field theory (QFT).

A granular material is a conglomeration of discrete solid, macroscopic particles characterized by a loss of energy whenever the particles interact. The constituents that compose granular material are large enough such that they are not subject to thermal motion fluctuations. Thus, the lower size limit for grains in granular material is about 1 μm. On the upper size limit, the physics of granular materials may be applied to ice floes where the individual grains are icebergs and to asteroid belts of the Solar System with individual grains being asteroids.

In quantum field theory, the LSZ reduction formula is a method to calculate S-matrix elements from the time-ordered correlation functions of a quantum field theory. It is a step of the path that starts from the Lagrangian of some quantum field theory and leads to prediction of measurable quantities. It is named after the three German physicists Harry Lehmann, Kurt Symanzik and Wolfhart Zimmermann.

In quantum mechanics, a two-state system is a quantum system that can exist in any quantum superposition of two independent quantum states. The Hilbert space describing such a system is two-dimensional. Therefore, a complete basis spanning the space will consist of two independent states. Any two-state system can also be seen as a qubit.

In mathematics, the Fubini–Study metric is a Kähler metric on projective Hilbert space, that is, on a complex projective space CPⁿ endowed with a Hermitian form. This metric was originally described in 1904 and 1905 by Guido Fubini and Eduard Study.

A quasiprobability distribution is a mathematical object similar to a probability distribution but which relaxes some of Kolmogorov's axioms of probability theory. Quasiprobabilities share several of general features with ordinary probabilities, such as, crucially, the ability to yield expectation values with respect to the weights of the distribution. They can however violate the σ-additivity axiom: integrating them over does not necessarily yield probabilities of mutually exclusive states. Indeed, quasiprobability distributions also counterintuitively have regions of negative probability density, contradicting the first axiom. Quasiprobability distributions arise naturally in the study of quantum mechanics when treated in phase space formulation, commonly used in quantum optics, time-frequency analysis, and elsewhere.

In statistical mechanics, the two-dimensional square lattice Ising model is a simple lattice model of interacting magnetic spins. The model is notable for having nontrivial interactions, yet having an analytical solution. The model was solved by Lars Onsager for the special case that the external magnetic field H = 0. An analytical solution for the general case for $has yet to be found.$

The Sudarshan-Glauber P representation is a suggested way of writing down the phase space distribution of a quantum system in the phase space formulation of quantum mechanics. The P representation is the quasiprobability distribution in which observables are expressed in normal order. In quantum optics, this representation, formally equivalent to several other representations, is sometimes championed over alternative representations to describe light in optical phase space, because typical optical observables, such as the particle number operator, are naturally expressed in normal order. It is named after George Sudarshan and Roy J. Glauber, who worked on the topic in 1963. Despite many useful applications in laser theory and coherence theory, the Glauber–Sudarshan P representation has the drawback that it is not always positive, and is not a true probability function.

Monte Carlo in statistical physics refers to the application of the Monte Carlo method to problems in statistical physics, or statistical mechanics.

Resonance fluorescence is the process in which a two-level atom system interacts with the quantum electromagnetic field if the field is driven at a frequency near to the natural frequency of the atom.

Photon antibunching generally refers to a light field with photons more equally spaced than a coherent laser field, a signature being signals at appropriate detectors which are anticorrelated. More specifically, it can refer to sub-Poissonian photon statistics, that is a photon number distribution for which the variance is less than the mean. A coherent state, as output by a laser far above threshold, has Poissonian statistics yielding random photon spacing; while a thermal light field has super-Poissonian statistics and yields bunched photon spacing. In the thermal (bunched) case, the number of fluctuations is larger than a coherent state; for an antibunched source they are smaller.

In probability theory and directional statistics, a wrapped normal distribution is a wrapped probability distribution that results from the "wrapping" of the normal distribution around the unit circle. It finds application in the theory of Brownian motion and is a solution to the heat equation for periodic boundary conditions. It is closely approximated by the von Mises distribution, which, due to its mathematical simplicity and tractability, is the most commonly used distribution in directional statistics.

Thermal fluctuations Random temperature-influenced deviations of particles from their average state

In statistical mechanics, thermal fluctuations are random deviations of a system from its average state, that occur in a system at equilibrium. All thermal fluctuations become larger and more frequent as the temperature increases, and likewise they decrease as temperature approaches absolute zero.

For many paramagnetic materials, the magnetization of the material is directly proportional to an applied magnetic field, for sufficiently high temperatures and small fields. However, if the material is heated, this proportionality is reduced. For a fixed value of the field, the magnetic susceptibility is inversely proportional to temperature, that is

The Mandel Q parameter measures the departure of the occupation number distribution from Poissonian statistics. It was introduced in quantum optics by Leonard Mandel. It is a convenient way to characterize non-classical states with negative values indicating a sub-Poissonian statistics, which have no classical analog. It is defined as the normalized variance of the boson distribution:

The Dolinar Receiver is a device based upon the Kennedy receiver that may be used to discriminate between two or more low-amplitude coherent states of light using displacements and adaptive measurements. The ability to discriminate signals encoded in coherent light has applications in communications where losses are unavoidable, such as transmission along fiber optics cable, through the atmosphere, or across deep space.

References

1 2 3 4 M. Fox, Quantum Optics: An Introduction, Oxford University Press, New York, 2006
↑ I. Deutsch, Quantum Optics Course Fall 2015, http://info.phys.unm.edu/~ideutsch/Classes/Phys566F15/Lectures/Phys566_Lect02.pdf. Retrieved 9 December 2015
1 2 Mandel, L (1959-09-01). "Fluctuations of Photon Beams: The Distribution of the Photo-Electrons". Proceedings of the Physical Society. IOP Publishing. 74 (3): 233–243. doi:10.1088/0370-1328/74/3/301. ISSN 0370-1328.
↑ J. W. Goodman, Statistical Optics, Wiley, New York, (1985) 238-256, 466-468
1 2 Schulte, Carsten H. H.; Hansom, Jack; Jones, Alex E.; Matthiesen, Clemens; Le Gall, Claire; Atatüre, Mete (2015-08-31). "Quadrature squeezed photons from a two-level system". Nature. Springer Science and Business Media LLC. 525 (7568): 222–225. arXiv: 1506.06827 . doi:10.1038/nature14868. ISSN 0028-0836.
↑ Vogel, Werner (1995-05-01). "Homodyne correlation measurements with weak local oscillators". Physical Review A. American Physical Society (APS). 51 (5): 4160–4171. doi:10.1103/physreva.51.4160. ISSN 1050-2947.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[Fox2006-1] 1 2 3 4 M. Fox, Quantum Optics: An Introduction, Oxford University Press, New York, 2006

[Deutsch-2] I. Deutsch, Quantum Optics Course Fall 2015, http://info.phys.unm.edu/~ideutsch/Classes/Phys566F15/Lectures/Phys566_Lect02.pdf. Retrieved 9 December 2015

[Mandel1959-3] 1 2 Mandel, L (1959-09-01). "Fluctuations of Photon Beams: The Distribution of the Photo-Electrons". Proceedings of the Physical Society. IOP Publishing. 74 (3): 233–243. doi:10.1088/0370-1328/74/3/301. ISSN 0370-1328.

[Good1985-4] J. W. Goodman, Statistical Optics, Wiley, New York, (1985) 238-256, 466-468

[Schulte2015-5] 1 2 Schulte, Carsten H. H.; Hansom, Jack; Jones, Alex E.; Matthiesen, Clemens; Le Gall, Claire; Atatüre, Mete (2015-08-31). "Quadrature squeezed photons from a two-level system". Nature. Springer Science and Business Media LLC. 525 (7568): 222–225. arXiv: 1506.06827 . doi:10.1038/nature14868. ISSN 0028-0836.

[Vogel1995-6] Vogel, Werner (1995-05-01). "Homodyne correlation measurements with weak local oscillators". Physical Review A. American Physical Society (APS). 51 (5): 4160–4171. doi:10.1103/physreva.51.4160. ISSN 1050-2947.

[1]

[2]

[3]

[4]

[5]

[6]