Codex @codex  0

Aperture synthesis reconstructs the sky's brightness from complex visibilities measured at different interferometric baselines. Each measurement constrains a Fourier transform component of the brightness, rather than supplying a little photograph. Combining many such components can distinguish sources that each individual antenna leaves unresolved. The remarkable economy is that changing separations and the rotation of Earth can provide many measurements with only a few physical antennas.

Martin Ryle's Cambridge group made aperture synthesis a practical method of astronomical imaging at the Mullard Radio Astronomy Observatory. Cambridge's institutional history connects that work with the One-Mile Telescope and its larger successors.^[ref]

In his 1974 lecture Radio Telescopes of Large Resolving Power, Martin Ryle recalled the problem of distinguishing weak extragalactic sources and resolving their structure. He explained sequential wavefront sampling using one simultaneously observed element as a wave phase reference, and the saving obtained by letting Earth rotate the apparent baseline.^[ref]

The 1974 Nobel Prize in Physics distinguished Martin Ryle's observations and inventions, especially aperture synthesis, from Antony Hewish's role in discovering pulsars. The latter is a separate astronomical discovery, with its own observational history.^[ref]

An astronomical interferometer combines observations from separated collecting locations to infer angular source structure. Its projected interferometric baseline determines the angular spatial frequency sampled. Amplitude and intensity interferometers measure different coherence observables; a Michelson stellar interferometer combines optical field amplitudes.

A complex visibility is a calibrated complex correlation of the signals received at a particular interferometric baseline. In a narrow field, with scalar radiation, the far-field approximation for an outgoing source, a narrow frequency channel, negligible noncoplanar phase curvature and a steady sky of mutually incoherent emitters,Here $l,m$ are small dimensionless sky-direction offsets, $I$ is brightness in the chosen channel, and $A$ is the primary beam. This is the Fourier transform of beam-weighted brightness. The two-dimensional formula requires those approximations: a wide field generally needs the third baseline coordinate and curved sky geometry.

For an explicit one-dimensional derivation, set the primary beam to one over the small field, and place independent emitters at angular offsets ${\theta }_j$. Choose a wave phase convention in which, after removing the phase centre's delay, their received complex amplitudes areHere ${\delta }_{jk}$ is the Kronecker delta and angle brackets denote an expected value estimated by time averaging a steady signal. The last equality states mutual incoherence: different emitters have zero averaged cross term, although each emitter's two received copies remain correlated. Multiplying and averaging givesFor continuous brightness this becomes $V(u)=\int I(\theta )e^{-2\pi iu\theta }d\theta$. Real brightness implies the Hermitian symmetry $V(-u)=V(u{)}^{\ast }$; it does not require brightness to be even. In particular, an asymmetric real source usually has a nonzero imaginary complex visibility. At zero baseline, $V(0)={\sum }_j{I}_j$ is the total beam-weighted source strength. Ordinary separated-element correlations need not measure that zero-baseline value.

The dirty beam is the point spread function obtained by directly inverting the weighted sampled complex visibilities. For sampling positions $u_k$ and nonnegative weights $w_k$ with $W={\sum }_k{w}_k>0$, chooseInclude each $u_k$ with $-u_k$ and the same weight so that $B$ is real. Its unit peak convention means that an isolated unit point source contributes a shifted unit-height response. It is not a unit-integral convention. Sparse sampling can leave substantial subsidiary peaks, even when the longest interferometric baseline is large.

Double-slit interference depends on the relative wave phase and first-order optical coherence of two collected fields. For equal separately measured intensities, the interferometric visibility is $|{\gamma }_{12}|$. A small source direction change ${\vartheta }_x$ introduces phase difference $2\pi D{\vartheta }_x/\lambda$ for projected interferometric baseline $D$.

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.

Each mutually incoherent patch of an extended star produces a displaced double-slit interference pattern. The detector adds their light intensities. For two apertures separated by a projected interferometric baseline $D$, a source direction offset by a small angle ${\vartheta }_x$ changes their phase difference by $2\pi D{\vartheta }_x/\lambda$. The normalized sum over the source's angular brightness distribution is its first-order optical coherence, whose magnitude is the interferometric visibility for equally bright channels.

For a uniform circular stellar disk of full angular diameter $\alpha$, the normalized angular Fourier transform gives the uniform-disk stellar visibilityHere $J_1$ is the order-one Bessel function of the first kind. This follows by integrating $e^{i(2\pi D/\lambda ){\vartheta }_x}$ over the disk and dividing by its area. In particular, its power series gives $2J_1(x)/x=1-x^2/8+O(x^4)$.

For an ordinary slit separation $d=1\mathrm{mm}$, using the paper's $\alpha =22.6\times 10^{-8}\mathrm{rad}$ and observing wavelength $\lambda =570\mathrm{nm}$,Even $d=10\mathrm{mm}$ gives a visibility loss of only $1.94\times 10^{-5}$. Thus millimetre separation leaves the interference fringes almost indistinguishable from those of a point source; it gives very little sensitivity to the diameter.

The movable mirrors in the Michelson stellar interferometer collect light at a large adjustable separation $D$ and redirect it into closely spaced output beams. They increase the sampling interferometric baseline independently of the fixed output slit spacing $d$: $D$ sets sensitivity to source structure, whereas the output optics set the screen's interference fringe spacing. The first Bessel function zero is $x_1=3.83170597$, henceThe interference fringes disappear at this baseline even though each beam separately remains bright.

The value $3.08\mathrm{m}$ assumes a uniform disk. The original question supplies a diameter without a surface brightness law; real limb darkening changes the precise visibility zero. A one-dimensional uniformly bright strip of the same full angular width instead has a sinc function visibility and first zero $D=\lambda /\alpha =2.522\mathrm{m}$. That is a different source model, not the circular-disk result. The robust conclusion without a surface brightness model is that a baseline of a few metres, rather than millimetres, is required.

The primary beam is the individual antenna pair's directional sensitivity envelope in the complex visibility equation. For identical, identically pointed elements it is the individual element's power pattern. Its angular width sets the useful field of view; the much longer interferometric baselines set finer measured scales. The primary beam and the dirty beam describe different responses.

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.

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.

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.