Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Problem 7.1 Consider a radiating time-harmonic point source in r1=(x1,0,0)\mathbf{r}_1=(x_1,0,0). The complex field at the observation point r=(x,y,z)\mathbf{r}=(x,y,z), where z>0z>0, is given by:

U1(x,y,z)=Q1eik(x−x1)2+y2+z2(x−x1)2+y2+z2.\begin{align*} U_1(x,y,z) = Q_1 \frac{e^{ik \sqrt{(x-x_1)^2 + y^2 + z^2}}}{\sqrt{(x-x_1)^2 + y^2 + z^2}}. \end{align*}

where Q1Q_1 is a given complex number whose modulus is proportional to the source strength.

(a) Derive that for fixed x1x_1 and sufficiently large z>0z>0 such that the source-dependent phase kx12/(2z)kx_1^2/(2z) is negligible, the field can be approximated by

U1,far(x,y,z)=Q1eikzzeikx2+y22ze−ikx1xz.\begin{align*} U_{1,far}(x,y,z) = Q_1 \frac{e^{ik z}}{z} e^{i k\frac{x^2+ y^2}{2z}} e^{-i \frac{k x_1 x}{z}}. \end{align*}

(b) Let there be a second point source at r2=(x2,0,0)\mathbf{r}_2=(x_2,0,0), with complex source strength Q2Q_2 with ∣Q2∣=∣Q1∣|Q_2|=|Q_1|. We assume that both point sources are coherent. This means that there is ϕ\phi such that Q2=Q1eiϕQ_2=Q_1 e^{i \phi}. Show that the field at (x,y,z)(x,y,z) for z>0z>0 large due to the two point sources can be written as

U(x,y,z)≈U1,far(x,y,z)(1+eiϕeikΔxxz),\begin{align*} U(x,y,z) \approx U_{1,far}(x,y,z)\left( 1 + e^{i\phi} e^{i \frac{k \Delta x x}{z}} \right), \end{align*}

where Δx=x1−x2\Delta x = x_1-x_2.

(c) For which small angles θ≃x/z\theta\simeq x/z does the far-field intensity vanish? Show that the distance between the point sources Δx\Delta x can be determined from the angular separation of the zeros of the intensity. Does the angular separation depend on the phase difference ϕ\phi?

(d) What should be the phase difference between the point sources such that for x/z=0x/z=0 on the screen at large distance zz the intensity vanishes?

Problem 7.2 Consider two slits of equal width aa in a non-transparent screen of thickness dd in the plane z=0z=0 as shown in Figure 1. The screen is illuminated by a plane wave with unit amplitude and propagating in the positive zz-direction. In the second slit there is a piece of glass with refractive index nn and thickness dd. In the first slit there is vacuum. Neglect interface reflections, absorption, and the small lateral displacement inside the glass; treat it as an ideal phase-only insert.

Two equal slits separated by b in an opaque screen, one containing glass, with diffracted light measured to the right.

Figure 1:Two slits with centers separated along xx and long sides parallel to yy in a dark screen of thickness dd. The lower slit is filled with glass, the upper is in vacuum.

(a) If the field immediately behind slit 1 has complex amplitude equal to 1, explain that the field immediately behind slit 2 is given by

eiϕ\begin{align*} e^{ i \phi} \end{align*}

with ϕ=k0(n−1)d\phi=k_0(n-1)d, where k0=2π/λ0k_0=2\pi/\lambda_0 is the vacuum wave number.

(b) Derive (using Eq. (2) or in another way) that the Fraunhofer intensity pattern on a screen along the line y=0y=0 at large distance zz is given by (up to factors that do not depend on x/zx/z).

Ifar(x,0,z)=2a2z2[sin⁡(kax2z)kax2z]2[1+cos⁡(kbxz+ϕ)].\begin{align*} I_{far}(x,0,z) = 2 \frac{a^2}{z^2} \left[ \frac{\sin \left(\frac{kax}{2z}\right)}{\frac{kax}{2z}}\right]^2 \left[ 1 + \cos\left( \frac{k b x}{z} +\phi\right)\right]. \end{align*}

In deriving this result you may omit all factors that are independent of x/zx/z and y/zy/z. If you use Eq. (3) you may take Q1=1Q_1=1, Q2=eiϕQ_2=e^{i \phi}.

(c) Make a sketch of this intensity pattern, showing the zeros and maxima as a function of θ=x/z\theta=x/z when a=2λa=2 \lambda, b=4λb=4 \lambda and ϕ=−π/2\phi=-\pi/2. Explain where the envelope and the other factor that depends on x/zx/z are caused by.

Problem 7.3 We consider the optical set-up shown at the left of Figure 2 where a point source at point rs=(a,0,0)\mathbf{r}_s=(a, 0,0) is above a plane mirror in the plane x=0x=0. The time-harmonic field emitted by the point source without the mirror being present is in complex notation given by:

Us(x,y,z,t)=eik(x−a)2+y2+z2−iωt(x−a)2+y2+z2.\begin{align*} U_s(x,y,z,t)= \frac{e^{i k \sqrt{ (x-a)^2 + y^2 + z^2}-i \omega t}}{\sqrt{(x-a)^2 + y^2 + z^2}}. \end{align*}
A source and an aperture above a horizontal mirror each send direct and reflected light to a distant screen.

Figure 2:Lloyd mirror configuration with a point source (left) and a rectangular aperture in a dark screen (right), above a mirror and with a screen at distance zz where the field is observed.

(a) Let UrU_r be the field reflected by the mirror. Assume that the mirror is perfect so that the total field Us(x,y,z)+Ur(x,y,z)U_s(x,y,z)+U_r(x,y,z) is zero on the surface of the mirror, i.e. when x=0x=0. Show that the reflected field UrU_r can be considered to be emitted by a point source in (−a,0,0)(-a,0,0), which is the image of the original point source by the mirror, and which is out of phase with the original point source and has the same strength.

(b) Consider the field on a screen at z>0z>0. According to (a) the field in the presence of the mirror can be considered to be radiated by two point sources, namely at (a,0,0)(a,0,0) and (−a,0,0)(-a,0,0), that are of equal strength and out of phase with each other. Assume that zz is so large that the spherical waves emitted by these point sources and arriving at the screen can both be considered to be plane waves. Derive that the point on the screen (x,0,z)(x,0,z) with smallest x>0x>0 where the field is zero is given by

x=λ2az,\begin{align*} x= \frac{\lambda}{2a} z, \end{align*}

where λ\lambda is the wavelength. In your derivation use path length differences of interfering rays and make a drawing.

(c) What happens with this zero and with the fringe pattern on the screen when the perfect mirror is replaced by a dielectric such as a piece of glass?

(d) Suppose now that there is a second point source at (2a,0,0)(2a,0,0) above the mirror and suppose that it has the same strength and is in phase with the point source in (a,0,0)(a,0,0). Assume again that the mirror is perfectly reflecting and derive (again by considering path length differences and using a drawing) that the smallest x>0x>0 for which a zero occurs at point (x,0,z)(x,0,z) on the screen is given by

x=λ3az\begin{align*} x= \frac{\lambda}{3a}z \end{align*}

(e) Derive the smallest x>0x>0 for which the field is zero at (x,0,z)(x,0,z) when the two point sources at (a,0,0)(a,0,0) and (2a,0,0)(2a,0,0) are mutually incoherent. Use again path length differences and make a drawing.

(f) Next suppose that there is a square aperture with center at (a,0,0)(a,0,0) and sides of length b<ab<a parallel to the xx- and yy-directions in an opaque (i.e. dark) screen above the mirror as shown at the right of Figure 2. The aperture is illuminated by a plane wave that propagates parallel to the zz-axis, hence the field in the aperture has constant phase and amplitude. Compute the smallest positive xx for which a zero occurs at (x,0,z)(x,0,z) on the screen at large distance z>0z>0. Use again path length differences and a drawing in your derivation.

Problem 7.4 Note: to answer the following questions it is NOT necessary to compute diffraction integrals.

(a) Consider two equally strong point sources which with respect to a coordinate system (x,y,z)(x,y,z) are at (−a/2,0,0)(-a/2,0,0) and (a/2,0,0)(a/2,0,0), where the zz-axis is the optical axis. Suppose the point sources are mutually coherent, emit in phase, and have separation a≫λa\gg\lambda so that the first zero is within the paraxial range. The maximum intensity on a screen at large (i.e. Fraunhofer) distance is then on the optical axis. Show that the smallest angle with the optical axis at which there is a zero on this screen is given by λ/(2a)\lambda/(2a).

(b) What is the smallest angle with the optical axis at which there is a zero on the screen when the two point sources emit with phase difference π/2\pi/2?

(c) Are there any zeros on the screen when the two point sources are mutually incoherent?

(d) Consider now two identical apertures in an opaque screen at z=0z=0. One aperture has its center at (−a/2,0,0)(-a/2,0,0) and the other has its center at (a/2,0,0)(a/2,0,0). The apertures are illuminated by a time-harmonic plane wave at perpendicular incidence to the screen (i.e. propagating parallel to the zz-axis). Show that the two-aperture interference factor vanishes at θ≃λ/(2a)\theta\simeq\lambda/(2a), whatever their identical shape. Could zeros of the single-aperture envelope occur at smaller angles?

(e) Suppose that the plane wave is incident at some angle different from 90o90^o. Let its complex field be given by

U(x,z)=ei(kxx+kzz)\begin{align*} U(x,z)= e^{i (k_x x + k_z z)} \end{align*}

where kx2+kz2=k\sqrt{k_x^2+k_z^2}=k. Suppose that kxa=π/2k_x a =\pi/2 (modulo 2π2\pi). Find the signed angle of the interference zero closest to the optical axis, assuming the single-aperture envelope does not vanish first. Explain how the oblique illumination shifts the pattern.

(f) Now imagine that both (identical) apertures are filled with glass plates with thickness that varies with position. The two glass plates are identical and they are identically positioned in each of the two apertures. Imagine that we illuminate the apertures with a plane wave at perpendicular incidence. The field transmitted by each aperture is now a rather complicated function of position, however the transmitted fields behind both apertures are identical. Does the two-aperture interference zero remain at θ≃λ/(2a)\theta\simeq\lambda/(2a), or do the identical glass plates shift it? What can change about the envelope? Explain your answer.

Problem 7.5 Bessel beams.

Suppose there is a mask in the entrance pupil of radius aa of a positive thin lens with image focal length fif_i with a thin ring-shaped aperture at r=br=b with width Δr\Delta r. If a plane wave with amplitude AA is at perpendicular incidence on the mask, the field immediately behind the mask is given by

UBessel(x,y)={A,       if b−Δr<x2+y2<b,0,       otherwise\begin{align*} U_{Bessel}(x,y)= \left\{ \begin{array}{l}A, \;\;\; \text{ if } b-\Delta r < \sqrt{x^2+y^2}<b, \\0, \;\;\; \text{ otherwise} \end{array}\right. \end{align*}

(a) Use the integral

∫02πeiζcos⁡ψdψ=2πJ0(ζ).\begin{align*} \int_0^{2\pi} e^{i \zeta \cos \psi} \mathrm{d}\psi = 2\pi J_0(\zeta). \end{align*}

to derive that for sufficiently small Δr\Delta r, the field in the focal plane, apart from the common lens-propagation prefactor and quadratic phase, is approximately

UBessel(x,y,fi)=2πAbΔrJ0(kbrfi).\begin{align*} U_{Bessel}(x,y,f_i)= 2\pi A b \Delta r J_0\left(k\frac{br}{f_i}\right). \end{align*}

(b) The beam obtained this way is called a Bessel beam. Explain why this beam has a very long focal depth.

(c) In the thin-ring approximation Δr≪b\Delta r\ll b, suppose that a full circular pupil illuminated with unit amplitude and the ring pupil carry the same incident power. Show that then the amplitude of the Bessel beam is given by

A=a2bΔr\begin{align*} A = \frac{a}{\sqrt{2 b \Delta r}} \end{align*}

(d) Derive that the ratio of the field amplitudes in the focal point of the Bessel beam: UBessel(0,0,fi)U_{Bessel}(0,0,f_i) and the Airy spot: UAiry(0,0,fi)U_{Airy}(0,0,f_i) is given by

UBesselUAiry=2bΔra.\begin{align*} \frac{U_{Bessel}}{U_{Airy}} = \frac{\sqrt{2 b \Delta r}}{ a}. \end{align*}

If b=ab=a and Δr=0.1a\Delta r=0.1 a this becomes

UBesselUAiry=2Δra≈0.45,\begin{align*} \frac{U_{Bessel }}{U_{Airy}} = \sqrt{\frac{2 \Delta r}{ a}} \approx 0.45, \end{align*}

The finite ring in the plotted example gives about 0.44 when its exact area is used; the displayed formula uses the thin-ring approximation. This is the case shown in Figure 3.

Radial plot comparing an Airy central peak with a lower Bessel peak and stronger Bessel side lobes.

Figure 3:Amplitude in the focal plane of a Bessel beam and of an Airy spot with the same total energy. The lens pupil has radius a=10000λa=10000\lambda, the ring aperture of the Bessel beam case is at the outer edge of the pupil (b=ab=a) and has width Δr=0.1a\Delta r = 0.1 a and the NA=0.1\text{NA}=0.1.

(e) The Bessel beam has stronger side lobes than the Airy spot. Explain the reason.

Problem 7.6 Stellar interferometry.

Observe a star through a narrow spectral band centered at vacuum wavelength λ\lambda. Model it as a distant spatially incoherent source with angular brightness B(θx,θy)B(\theta_x,\theta_y), where (θx,θy)(\theta_x,\theta_y) specifies the transverse propagation direction of light arriving at Earth. Assume paraxial propagation, negligible temporal decorrelation across the instrument paths, and equal collection efficiencies. Let two collectors have transverse positions r1\mathbf r_1 and r2\mathbf r_2 and baseline b=r1−r2\mathbf b=\mathbf r_1-\mathbf r_2.

(a) Show that one source element at direction θ=(θx,θy)\boldsymbol\theta=(\theta_x,\theta_y) contributes a plane-wave factor exp⁡[i(2π/λ)θ⋅r]\exp[i(2\pi/\lambda)\boldsymbol\theta\cdot\mathbf r] at collector position r\mathbf r, apart from a common phase.

(b) Use the independence of different source elements to derive the mutual coherence

Γ(b)∝∬B(θ)ei(2π/λ)b⋅θ d2θ.\Gamma(\mathbf b)\propto \iint B(\boldsymbol\theta) e^{i(2\pi/\lambda)\mathbf b\cdot\boldsymbol\theta} \,\mathrm d^2\boldsymbol\theta.

State the normalized complex visibility V(b)=Γ(b)/Γ(0)V(\mathbf b)=\Gamma(\mathbf b)/\Gamma(\mathbf 0).

(c) If the two collected beams have equal mean intensity, how is the measured fringe contrast related to ∣V(b)∣|V(\mathbf b)|? What additional phase information is needed to reconstruct an image from measurements at many baselines?

(d) Estimate the angular scale probed by a largest baseline Bmax⁡B_{\max}. Explain why limited baseline coverage and measurement noise also affect the reconstructed image.