Can we apply Rayleigh's criterion for the resolution of a primary mirror of a reflective telescope?
-
$\begingroup$ Rayleigh's criterion of angle of resolution $\endgroup$– Tapan GuptaCommented Feb 20, 2023 at 4:15
-
1$\begingroup$ You can use edit to amend your question directly or to add clarifications $\endgroup$– planetmakerCommented Feb 20, 2023 at 5:20
3 Answers
Can we apply Rayleigh's criterion for the resolution of a primary mirror of a reflective telescope?
Yes, but...
You've asked about the criterion, not the use of the 1.22 equation.
From hyperphysics (and pretty much everywhere else):
The Rayleigh criterion is the generally accepted criterion for the minimum resolvable detail - the imaging process is said to be diffraction-limited when the first diffraction minimum of the image of one source point coincides with the maximum of another.
So as long as you can identify a "maximum" in one spot, and a "first diffraction minimum" in the other, you can define what "resolved" means.
However, to distinguish this from the other answers, this does not mean that you can always necessarily use the $\theta = 1.22 \lambda/d$ equation blindly. The 1.22 comes from the Airy disk which is the diffraction pattern from an optical system with a single, clean, circular aperture.
Some telescopes have pretty strange apertures, central obstructions, secondary mounts or other pathologies in their optics, and so their diffraction patterns differ from what you would calculate for a clear circular aperture's Airy disk.
Below is a quickie calculation for an aperture the size of the Vera C. Ruben observatory (8.360 meters). I've done it twice - once without the huge 5.116 meter central obstruction and once with, at 550 nm. We can see that the first theoretical minimum is at 8.25 mas for the unobstructed aperture, but decreases to 6.4 mas with the obstruction!
That's a ~30% improvement in resolution based on Rayleigh's criterion, just by putting a giant hole in the mirror!
That's not a complete surprise, you can always engineer an aperture to boost the modulation transfer function at one spatial frequency if you sacrifice others - in this case we can see that the additional outer rings (more widely scattered light) are much higher in intensity for the obstructed aperture. See also
- this answer about aperture interferometry
- Under what situations can an aperture mask improve the resolution of a small/medium amateur telescope? Is this demonstrable mathematically?
import numpy as np
import matplotlib.pyplot as plt
r2, r1 = 5.116, 8.360 # meters
lam = 550E-09
n = 1000
ym, xm = np.mgrid[-8.4:8.4:(2*n+1)*1j, -8.4:8.4:(2*n+1)*1j]
rsq = xm**2 + ym**2
ap_clear, ap_obstructed = np.zeros_like(rsq), np.zeros_like(rsq)
ap_clear[(rsq <= r1**2)] = 1.
ap_obstructed[(rsq >= r2**2) * (rsq <= r1**2)] = 1.
theta_max = 5 * lam / (2 * r1)
thetas = np.linspace(0, theta_max, 501)
thetas_arcsec = np.degrees(thetas) * 3600
I_clear, I_obstructed = [], []
for theta in thetas:
phase = 2 * np.pi * xm * theta / lam # avoid complex by symmetry (xm=0 in center)
E_clear = np.sqrt(ap_clear) * np.cos(phase) # sqrt for amplitude
E_obstructed = np.sqrt(ap_obstructed) * np.cos(phase) # sqrt for amplitude
I_c = np.abs(E_clear.sum())**2
I_o = np.abs(E_obstructed.sum())**2
I_clear.append(I_c)
I_obstructed.append(I_o)
I_clear = np.array(I_clear)
I_obstructed = np.array(I_obstructed)
Imax = max(I_clear.max(), I_obstructed.max())
I_clear /= Imax
I_obstructed /= Imax
if True:
extent = [xm.min(), xm.max(), ym.min(), ym.max()]
fig = plt.figure(figsize=[14, 7.5])
gs = fig.add_gridspec(nrows=2, ncols=2)
ax_clear = fig.add_subplot(gs[0, 0])
ax_obstructed = fig.add_subplot(gs[1, 0])
ax_plot = fig.add_subplot(gs[:, 1:])
ax_clear.imshow(ap_clear, cmap='gray', extent=extent)
ax_clear.set_title('clear')
ax_clear.set_ylabel('meters')
ax_obstructed.imshow(ap_obstructed, cmap='gray', extent=extent)
ax_obstructed.set_title('obstructed')
ax_obstructed.set_ylabel('meters')
ax_plot.plot(thetas_arcsec, np.log10(I_clear), '-', label='clear')
ax_plot.plot(thetas_arcsec, np.log10(I_obstructed), '--', label='obstructed')
ax_plot.set_xlim(0, 0.03)
ax_plot.set_xlabel('arcsec')
ax_plot.set_ylabel('log(I)')
plt.suptitle('Vera C. Ruben (5.116, 8.360 meters, λ=550nm)')
plt.legend()
plt.show()
# first minimum
# clear: .00825 arcsec
# obstructed: .00640 arcsec
Rayleigh's criterion concerns whether two point-sources will appear to merge by diffraction into a single dot. It is independent of the method of forming an image. So applies equally to mirror- and lens-based telescopes.
-
$\begingroup$ So ,in a reflector telescope, should I calculate the angle of resolution of primary mirror or eyepiece $\endgroup$ Commented Feb 20, 2023 at 7:13
-
1$\begingroup$ What would you do for a refraction telescope, the objective lens or the eyepiece? $\endgroup$– James KCommented Feb 20, 2023 at 7:22
-
3
Rayleigh’s criterion is good for both refractors and reflectors.