Uniform-disk stellar visibility (source code)

= Uniform-disk stellar visibility
{title2=$\gamma(B)=2J_1(x)/x,\quad x=\pi B\theta/\lambda$}

Stipulate a distant circular source of uniform angular brightness, full <angular diameter> $\theta\ll1$ 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 $\boldsymbol\alpha$ from the disk centre, the baseline-dependent <phase difference> is $2\pi\mathbf b_\perp\cdot\boldsymbol\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:
$$
\gamma(B)=\frac{1}{\pi a^2}\int_{|\boldsymbol\alpha|\leq a}
 e^{i2\pi\mathbf b_\perp\cdot\boldsymbol\alpha/\lambda}\,d^2\boldsymbol\alpha,
\qquad a=\frac\theta2.
$$
Choose the horizontal coordinate along $\mathbf b_\perp$, put $k=2\pi B/\lambda$, and integrate along vertical chords of the disk:
$$
\gamma(B)=\frac{2}{\pi a^2}\int_{-a}^{a}\sqrt{a^2-u^2}\,e^{iku}\,du
=\frac{2J_1(ka)}{ka}=\frac{2J_1(x)}x,
\qquad x=\frac{\pi B\theta}{\lambda}.
$$
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 series
$$
J_1(x)=\sum_{m=0}^{\infty}\frac{(-1)^m(x/2)^{2m+1}}{m!(m+1)!},
\qquad \gamma(B)=1-\frac{x^2}{8}+\frac{x^4}{192}+O(x^6).
$$
At $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\approx3.831705970$, so
$$
B_0=\frac{x_1\lambda}{\pi\theta}\approx\frac{1.21967\lambda}{\theta},
\qquad \theta\approx\frac{1.21967\lambda}{B_0}.
$$
This 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\times10^{-8}\,\mathrm{rad}$ and $B_0\approx25.1575\,\mathrm m$. At $B=10\,\mathrm m$, the model gives $x\approx1.52309$ and $V^2\approx0.542783$.

\Image[/stellar-uniform-disk-correlation-baseline.png]
{title=Stipulated uniform-disk model: signed field coherence, squared visibility and finite-gate intensity correlation versus <interferometric baseline>}
{description=Original calculated curves for a 5-mas disk at 500 nm, with no measured points. The lower panel uses one polarization, independent signal fractions 0.8 in each channel, and the explicitly specified exponential temporal correlation and 50-ps gate below. It is not a fit to Sirius.}

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.