Stellar intensity interferometry

Stellar intensity interferometry measures a star's angular diameter by comparing fluctuations of detected light intensity at separated collectors. Each collector has its own photodetector; an analog multiplier or a digital correlator compares the electrical outputs. Changing the interferometric baseline changes the sampled spatial optical coherence. An unresolved star produces correlated fluctuations; a resolved extended source supplies different mixtures of emission directions, reducing their correlation.

This makes a useful telescope out of a surprising observable: the fluctuations discarded by an ordinary average-brightness measurement. Optical beams need not travel to a common combiner. A relative wave phase perturbation $E_2\mapsto e^{i\varphi }{E}_2$ leaves $|E_2{|}^2$ unchanged, so pure phase errors do not directly change the local light intensity. Electronic delay alignment remains necessary; intensity changes from atmospheric scintillation, collection losses, and instrumental noise still require controls.

Let $\mathbf{b}$ join the collectors and let ${\mathbf{s}}_0$ be the unit line of sight toward the star. The relevant interferometric baseline is its projection onto the plane of the sky:The reference arrival-time difference has magnitude $|\mathbf{b}\cdot {\mathbf{s}}_0|/c$, with speed of light $c$. Correct that delay, with a sign fixed by channel order, before averaging the residual correlation at $\tau =0$. The physical collector separation equals $B$ only when it is perpendicular to ${\mathbf{s}}_0$; tracking a star can change the projection even for fixed collectors.

In their 1956 Sirius report, Robert Hanbury Brown and Richard Q. Twiss used two 1.56-m searchlight reflectors feeding photomultipliers. The photocurrent fluctuations passed through 5–45-MHz amplifiers; an adjustable delay aligned arrival times. A linear mixer multiplied the outputs; an integrating motor accumulated their product. A second motor monitored its root-mean-square fluctuation to control gain changes. Separate lamps illuminated capped photomultipliers before and after observations to test electronic coupling and drift; blank-sky runs tested correlated sky background. They fitted a uniformly bright disk and an unknown zero-baseline scale to relative correlations at different interferometric baselines. Broadband response and possible baseline-dependent systematic errors limited the diameter inference.^[ref]

This astronomical test extends the preceding laboratory Hanbury Brown and Twiss effect to two independently collected samples of a star. The historical motor reading was an instrument-specific correlation statistic. It should not be identified numerically with the normalized $g^{(2)}$ below without the corresponding measurement calibration.

For one selected polarization of waves, write the detected complex fields as $E_1(t),E_2(t+\tau )$ and their light intensities as $I_j=|E_j{|}^2$, absorbing a common conversion constant. With stationary ensemble averages and nonzero ${\bar{I}}_j=⟨I_j⟩$, defineHere second-order optical coherence concerns distinct detector channels, with the corresponding normal ordering in the quantum formulation. Independent detection has $g^{(2)}=1$. Subtracting the means gives the normalized covarianceThus multiplying raw photocurrents contains a large mean-product contribution; multiplying their fluctuations isolates the small excess. Independent detector shot noise does not supply a mean cross-channel excess, although it increases measurement uncertainty.

The displayed Siegert relation requires a zero-mean, jointly circular complex Gaussian random field, the usual chaotic thermal light model. Circular means invariance under a common arbitrary complex phase rotation; in particular $⟨E_i{E}_j⟩=0$. Apply the Gaussian moment pairing theorem to the real and imaginary parts:The remaining pairing vanishes by circularity. Division proves the Siegert relation. Gaussianity is an assumption about fourth moments, not a consequence of merely measuring first-order optical coherence; a coherent state or photon antibunching cannot be interpreted by this thermal formula.

For equal-intensity beams combined in a common optical mode with matched wave polarizations, interferometric visibility is $V=|g_{12}^{(1)}|$. Consequently ideal thermal stellar intensity interferometry measures $g^{(2)}-1=V^2$, while the full $g^{(2)}$ includes the independent-count reference one. Unpolarized thermal light with two independent equally populated wave polarizations has half this excess: adding their covariances doubles the numerator while multiplying total means quadruples the denominator. Finite detector response further suppresses the peak and must be calibrated.^[ref]

Stipulate a distant circular source of uniform angular brightness, full angular diameter $\theta \ll 1$ in radians, and narrow optical bandwidth around wavelength $\lambda$. Distinct source patches emit independently. Assume each collector is small enough not to resolve the disk individually. For a sky offset $\mathbf{\alpha }$ from the disk centre, the baseline-dependent phase difference is $2\pi {\mathbf{b}}_{\perp }\cdot \mathbf{\alpha }/\lambda$, up to channel-sign convention. Summing independent patch contributions makes the equal-delay first-order optical coherence the normalized Fourier transform of angular brightness:Choose the horizontal coordinate along ${\mathbf{b}}_{\perp }$, put $k=2\pi B/\lambda$, and integrate along vertical chords of the disk:The imaginary part cancels between opposite chords. Here $J_1$ is the order-one Bessel function of the first kind. One way to obtain the last equality is to integrate the exponential's power series term by term on the finite disk. Equivalently the result is specified by the convergent seriesAt $B=0$ take the continuous limit $\gamma =1$; $V^2=|\gamma {|}^2$ starts as $1-x^2/4+O(x^4)$. This is why very short interferometric baselines give little leverage on angular diameter. Larger disks lose spatial optical coherence at shorter baselines.

The first positive zero of $J_1$ is $x_1\approx 3.831705970$, soThis is a full diameter convention. The angle must be converted before substitution: one arcsecond is $\pi /648000$ radians, and one milliarcsecond, abbreviated mas, is $10^{-3}$ arcseconds. For the stipulated $\theta =5\mathrm{mas}$ and $\lambda =500\mathrm{nm}$, $\theta =2.4240684\times 10^{-8}\mathrm{rad}$ and $B_0\approx 25.1575\mathrm{m}$. At $B=10\mathrm{m}$, the model gives $x\approx 1.52309$ and $V^2\approx 0.542783$.

Beyond $B_0$, $\gamma$ changes sign but $V^2$ stays nonnegative and develops weak further lobes. A single measured squared visibility can therefore match several angular diameters; sampling several interferometric baselines constrains the fit. Limb darkening, spectral averaging, or resolved collector apertures change the uniform-disk prediction. A uniform-disk diameter is a parameter of that brightness model; converting it to a physical stellar radius also requires distance and an appropriate brightness profile.

An exact teaching model makes the peak suppression explicit. Assume the narrowband field correlation separates as $g_{12}^{(1)}(B,u)=\gamma (B)g_t^{(1)}(u)$ after delay compensation. Average each channel's light intensity over the same rectangular gate of duration $T$. The two time integrals have triangular relative-delay weightIntegrating the Siegert relation against this weight givesThe factor $p$ is one for a selected polarization of waves and $1/2$ for two independent equally populated polarizations. The signal fraction ${\rho }_j=s_j/(s_j+b_j)$ uses mean stellar signal $s_j$ and independent additive background $b_j$ in each detector. The photon correlation loss and background model then supplies ${\rho }_1{\rho }_2$. This assumes no correlated sky background or electronic coupling.

For the original model $|g_t^{(1)}(u){|}^2=e^{-|u|/{\tau }_0}$, elementary integration givesThe area-based optical coherence time of this intensity envelope is ${\tau }_c=\int |g_t^{(1)}(u){|}^{2}du=2{\tau }_0$. For $T\gg {\tau }_0$, $\beta \approx {\tau }_c/T$; for $T\ll {\tau }_0$, $\beta \to 1$. With the chosen ${\tau }_0=1\mathrm{ps}$, $T=50\mathrm{ps}$, $p=1$, and ${\rho }_1={\rho }_2=0.8$, $\beta \approx 0.0392$ and $g_{\mathrm{obs}}^{(2)}(0,0)=1.025088$. At $B=10\mathrm{m}$ the stipulated disk gives approximately $1.013617$. A nearly unit measured $g^{(2)}$ can therefore conceal a substantial $V^2$.

This rectangular-gate calculation specifies a counting model. An analog instrument's amplifier bandwidth and impulse response instead weight the fluctuations it actually passes; the Sirius amplifiers were bandpass devices. Both require measuring the response and calibrating the zero-baseline excess. When $p$, ${\rho }_j$, and the temporal response are unchanged, the ratio of excesses recovers $V^2$:A sufficiently unresolved calibration source can supply the denominator; dividing the full $g^{(2)}$ values would give the wrong result. Narrower optical bandwidth lengthens optical coherence time, whereas wider usable electronic bandwidth can resolve faster fluctuations; these are distinct bandwidths. Background changes, unequal spectral responses, drift, or averaging over changing interferometric baselines can invalidate a common calibration. Resolving a disk means establishing a reproducible decline of calibrated correlation, supported by those controls.