# Hanbury Brown and Twiss effect

(Redirected from Hanbury-Brown and Twiss effect)

In physics, the Hanbury Brown and Twiss (HBT) effect is any of a variety of correlation and anti-correlation effects in the intensities received by two detectors from a beam of particles. HBT effects can generally be attributed to the wave–particle duality of the beam, and the results of a given experiment depend on whether the beam is composed of fermions or bosons. Devices which use the effect are commonly called intensity interferometers and were originally used in astronomy, although they are also heavily used in the field of quantum optics.

## History

In 1956, Robert Hanbury Brown and Richard Q. Twiss published A test of a new type of stellar interferometer on Sirius, in which two photomultiplier tubes (PMTs), separated by about 6 meters, were aimed at the star Sirius. Light collected into the PMTs using mirrors from searchlights. An interference effect was observed between the two intensities, revealing a positive correlation between the two signals, despite the fact that no phase information was collected. Hanbury Brown and Twiss used the interference signal to determine the apparent angular size of Sirius, claiming excellent resolution.

Also, in the field of particle physics, Goldhaber et al. performed an experiment in 1959 in Berkeley and found an unexpected angular correlation among identical pions, discovering the ρ0 resonance, by means of ${\displaystyle \rho ^{0}\to \pi ^{-}\pi ^{+}}$ decay.[1] From then on, the HBT technique started to be used by the heavy-ion community to determine the space–time dimensions of the particle emission source for heavy-ion collisions. For recent developments in this field, see for example the review article by Lisa.[2]

An example of an intensity interferometer that would observe no correlation if the light source is a coherent laser beam, and positive correlation if the light source is a filtered one-mode thermal radiation. The theoretical explanation of the difference between the correlations of photon pairs in thermal and in laser beams was first given by Roy J. Glauber, who was awarded the 2005 Nobel Prize in Physics "for his contribution to the quantum theory of optical coherence".

The original HBT result met with much skepticism in the physics community. Although intensity interferometry had been widely used in radio astronomy, where Maxwell's equations are valid, at optical wavelengths the light would be quantised into a relatively small number of photons. Many physicists worried that the correlation was inconsistent with the laws of thermodynamics. Some even claimed that the effect violated the uncertainty principle. Hanbury Brown and Twiss resolved the dispute in a neat series of articles (see References below) that demonstrated, first, that wave transmission in quantum optics had exactly the same mathematical form as Maxwell's equations, albeit with an additional noise term due to quantisation at the detector, and second, that according to Maxwell's equations, intensity interferometry should work. Others, such as Edward Mills Purcell immediately supported the technique, pointing out that the clumping of bosons was simply a manifestation of an effect already known in statistical mechanics. After a number of experiments, the whole physics community agreed that the observed effect was real.

The original experiment used the fact that two bosons tend to arrive at two separate detectors at the same time. Morgan and Mandel used a thermal photon source to create a dim beam of photons and observed the tendency of the photons to arrive at the same time on a single detector. Both of these effects used the wave nature of light to create a correlation in arrival time – if a single photon beam is split into two beams, then the particle nature of light requires that each photon is only observed at a single detector, and so an anti-correlation was observed in 1986[where?]. Finally, bosons have a tendency to clump together, giving rise to Bose–Einstein correlations, while fermions due to the Pauli exclusion principle, tend to spread apart, leading to Fermi–Dirac (anti)correlations. Bose–Einstein correlations have been observed between pions, kaons and photons, and Fermi–Dirac (anti)correlations between protons, neutrons and electrons. For a general introduction in this field, see the textbook on Bose–Einstein correlations by Richard M. Weiner[3] A difference in repulsion of Bose–Einstein condensate in the "trap-and-free fall" analogy of the HBT effect[4] affects comparison.

## Wave mechanics

The HBT effect can, in fact, be predicted solely by treating the incident electromagnetic radiation as a classical wave. Suppose we have a monochromatic wave with frequency ${\displaystyle \omega }$ on two detectors, with an amplitude ${\displaystyle E(t)}$ that varies on timescales slower than the wave period ${\displaystyle 2\pi /\omega }$. (Such a wave might be produced from a very distant point source with a fluctuating intensity.)

Since the detectors are separated, say the second detector gets the signal delayed by a time ${\displaystyle \tau }$, or equivalently, a phase ${\displaystyle \phi =\omega \tau }$; that is,

${\displaystyle E_{1}(t)=E(t)\sin(\omega t),}$
${\displaystyle E_{2}(t)=E(t-\tau )\sin(\omega t-\phi ).}$

The intensity recorded by each detector is the square of the wave amplitude, averaged over a timescale that is long compared to the wave period ${\displaystyle 2\pi /\omega }$ but short compared to the fluctuations in ${\displaystyle E(t)}$:

{\displaystyle {\begin{aligned}i_{1}(t)&={\overline {E_{1}(t)^{2}}}={\overline {E(t)^{2}\sin ^{2}(\omega t)}}={\tfrac {1}{2}}E(t)^{2},\\i_{2}(t)&={\overline {E_{2}(t)^{2}}}={\overline {E(t-\tau )^{2}\sin ^{2}(\omega t-\phi )}}={\tfrac {1}{2}}E(t-\tau )^{2},\end{aligned}}}

where the overline indicates this time averaging. For wave frequencies above a few terahertz (wave periods less than a picosecond), such a time averaging is unavoidable, since detectors such as photodiodes and photomultiplier tubes cannot produce photocurrents that vary on such short timescales.

The correlation function ${\displaystyle \langle i_{1}i_{2}\rangle (\tau )}$ of these time-averaged intensities can then be computed:

{\displaystyle {\begin{aligned}\langle i_{1}i_{2}\rangle (\tau )&=\lim _{T\to \infty }{\frac {1}{T}}\int \limits _{0}^{T}i_{1}(t)i_{2}(t)\,\mathrm {d} t\\&=\lim _{T\to \infty }{\frac {1}{T}}\int \limits _{0}^{T}{\tfrac {1}{4}}E(t)^{2}E(t-\tau )^{2}\,\mathrm {d} t.\end{aligned}}}

Most modern schemes actually measure the correlation in intensity fluctuations at the two detectors, but it is not too difficult to see that if the intensities are correlated, then the fluctuations ${\displaystyle \Delta i=i-\langle i\rangle }$, where ${\displaystyle \langle i\rangle }$ is the average intensity, ought to be correlated, since

{\displaystyle {\begin{aligned}\langle \Delta i_{1}\Delta i_{2}\rangle &={\big \langle }(i_{1}-\langle i_{1}\rangle )(i_{2}-\langle i_{2}\rangle ){\big \rangle }=\langle i_{1}i_{2}\rangle -{\big \langle }i_{1}\langle i_{2}\rangle {\big \rangle }-{\big \langle }i_{2}\langle i_{1}\rangle {\big \rangle }+\langle i_{1}\rangle \langle i_{2}\rangle \\&=\langle i_{1}i_{2}\rangle -\langle i_{1}\rangle \langle i_{2}\rangle .\end{aligned}}}

In the particular case that ${\displaystyle E(t)}$ consists mainly of a steady field ${\displaystyle E_{0}}$ with a small sinusoidally varying component ${\displaystyle \delta E\sin(\Omega t)}$, the time-averaged intensities are

{\displaystyle {\begin{aligned}i_{1}(t)&={\tfrac {1}{2}}E_{0}^{2}+E_{0}\,\delta E\sin(\Omega t)+{\mathcal {O}}(\delta E^{2}),\\i_{2}(t)&={\tfrac {1}{2}}E_{0}^{2}+E_{0}\,\delta E\sin(\Omega t-\Phi )+{\mathcal {O}}(\delta E^{2}),\end{aligned}}}

with ${\displaystyle \Phi =\Omega \tau }$, and ${\displaystyle {\mathcal {O}}(\delta E^{2})}$ indicates terms proportional to ${\displaystyle (\delta E)^{2}}$, which are small and may be ignored.

The correlation function of these two intensities is then

{\displaystyle {\begin{aligned}\langle \Delta i_{1}\Delta i_{2}\rangle (\tau )&=\lim _{T\to \infty }{\frac {(E_{0}\delta E)^{2}}{T}}\int \limits _{0}^{T}\sin(\Omega t)\sin(\Omega t-\Phi )\,\mathrm {d} t\\&={\tfrac {1}{2}}(E_{0}\delta E)^{2}\cos(\Omega \tau ),\end{aligned}}}

showing a sinusoidal dependence on the delay ${\displaystyle \tau }$ between the two detectors.

## Quantum interpretation

Photon detections as a function of time for a) antibunching (e.g. light emitted from a single atom), b) random (e.g. a coherent state, laser beam), and c) bunching (chaotic light). τc is the coherence time (the time scale of photon or intensity fluctuations).

The above discussion makes it clear that the Hanbury Brown and Twiss (or photon bunching) effect can be entirely described by classical optics. The quantum description of the effect is less intuitive: if one supposes that a thermal or chaotic light source such as a star randomly emits photons, then it is not obvious how the photons "know" that they should arrive at a detector in a correlated (bunched) way. A simple argument suggested by Ugo Fano [Fano, 1961] captures the essence of the quantum explanation. Consider two points ${\displaystyle a}$ and ${\displaystyle b}$ in a source that emit photons detected by two detectors ${\displaystyle A}$ and ${\displaystyle B}$ as in the diagram. A joint detection takes place when the photon emitted by ${\displaystyle a}$ is detected by ${\displaystyle A}$ and the photon emitted by ${\displaystyle b}$ is detected by ${\displaystyle B}$ (red arrows) or when ${\displaystyle a}$'s photon is detected by ${\displaystyle B}$ and ${\displaystyle b}$'s by ${\displaystyle A}$ (green arrows). The quantum mechanical probability amplitudes for these two possibilities are denoted by ${\displaystyle \langle A|a\rangle \langle B|b\rangle }$ and ${\displaystyle \langle B|a\rangle \langle A|b\rangle }$ respectively. If the photons are indistinguishable, the two amplitudes interfere constructively to give a joint detection probability greater than that for two independent events. The sum over all possible pairs ${\displaystyle a,b}$ in the source washes out the interference unless the distance ${\displaystyle AB}$ is sufficiently small.

Two source points a and b emit photons detected by detectors A and B. The two colors represent two different ways to detect two photons.

Fano's explanation nicely illustrates the necessity of considering two-particle amplitudes, which are not as intuitive as the more familiar single-particle amplitudes used to interpret most interference effects. This may help to explain why some physicists in the 1950s had difficulty accepting the Hanbury Brown and Twiss result. But the quantum approach is more than just a fancy way to reproduce the classical result: if the photons are replaced by identical fermions such as electrons, the antisymmetry of wave functions under exchange of particles renders the interference destructive, leading to zero joint detection probability for small detector separations. This effect is referred to as antibunching of fermions [Henny, 1999]. The above treatment also explains photon antibunching [Kimble, 1977]: if the source consists of a single atom, which can only emit one photon at a time, simultaneous detection in two closely spaced detectors is clearly impossible. Antibunching, whether of bosons or of fermions, has no classical wave analog.

From the point of view of the field of quantum optics, the HBT effect was important to lead physicists (among them Roy J. Glauber and Leonard Mandel) to apply quantum electrodynamics to new situations, many of which had never been experimentally studied, and in which classical and quantum predictions differ.