1
$\begingroup$

Can we apply Rayleigh's criterion for the resolution of a primary mirror of a reflective telescope?

$\endgroup$
2
  • $\begingroup$ Rayleigh's criterion of angle of resolution $\endgroup$ Commented Feb 20, 2023 at 4:15
  • 1
    $\begingroup$ You can use edit to amend your question directly or to add clarifications $\endgroup$ Commented Feb 20, 2023 at 5:20

3 Answers 3

3
$\begingroup$

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


simulated diffraction pattern of Vera C. Ruben (5.116, 8.360 meters, λ=550nm)

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
$\endgroup$
3
$\begingroup$

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.

$\endgroup$
3
  • $\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 K
    Commented Feb 20, 2023 at 7:22
  • 3
    $\begingroup$ Objective lens...? $\endgroup$ Commented Feb 21, 2023 at 2:01
2
$\begingroup$

Rayleigh’s criterion is good for both refractors and reflectors.

$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .