Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Anticaustics in a Fabry–Perot interferometer

Open Access Open Access

Abstract

We address the response of a Fabry–Perot interferometer to a monochromatic point source. We calculate the anticaustics (that is, the virtual wavefronts of null path difference) resulting from the successive internal reflections occurring in the system. They turn out to be a family of ellipsoids (or hyperboloids) of revolution, which allows us to reinterpret the operation of the Fabry–Perot interferometer from a geometrical point of view that facilitates comparison with other apparently disparate arrangements, such as Young’s double slit.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. INTRODUCTION

Young’s double slit experiment [1] is the epitome of the phenomenon of interference, and it is still found in any modern optics textbook [2]. Initially, the experiment was a crucial step in establishing the wave theory of light. Ironically, it played a quintessential role in elaborating the concept of wave–particle duality, the “central mystery” of quantum theory [3].

The standard interpretation of Young’s experiment is based on the superposition of the wavefronts emanating from two coherent point sources. Similar arrangements were also developed during the 19th century, including Fresnel’s biprism [4], Lloyd’s mirrors [5] and Billet’s split lens [6], among others. In all of these cases, two point sources (real or virtual) are required to explain the observed intensity pattern.

One might rightly wonder how this picture is modified when many mutually coherent waves are superposed. The Fabry–Perot (FP) interferometer is the paradigm of this situation. Apart from its luminosity, the distinctive feature of this setup is its narrow resonances, which is the basis for its extensive use in high-resolution spectroscopy, interferometry, and laser resonators [7,8]. However, the ordinary approach assumes plane waves, that is, the source is at infinity. Therefore, it seems clear that to appreciate the analogies and differences between two- and multiple-source interference, one should look at the response of the FP to a point source.

In this paper, we take this unorthodox viewpoint. To this end, we analyze the wavefronts generated at each of the multiple reflections occurring in the FP. Although methods for obtaining the local properties of the refracted wavefronts are well established, the global properties are more elusive. The reason is that, in general, wavefronts may be very intricate surfaces, particularly in the neighborhood of a focus, and their analytic expressions might become quite involved [9].

However, it is often possible to extract from a family of wavefronts one that is simpler in form, whose global properties are more comprehensible, but that, nonetheless, is representative of the entire family. It corresponds to the privileged wavefront of zero algebraic optical path difference. This concept has been reinvented under a variety of names; an indicative list of dates where it appears with different designations has been compiled by Ouellette [10]: anticaustic [11], secondary caustic [12], orthogonal trajectory [13], simplest orthotomic curve [14], zero phase-front [15], emerging wavefront of null optical path [16], archetypical wavefront [17], zero-distance phase front [18,19], and phase front [20].

In a comprehensive historical review, Chastang and Farouki [21] endorse the use of the anticaustic, returning to the original Bernouilli’s suggestion. This was strongly supported by the late Emil Wolf [22]. We believe that reintroducing the term is appropriate, and constitutes a tribute to a great geometer and to a great physicist.

The anticaustic has been explicitly calculated in a few simple cases, including refraction by a plane (ellipse/hyperbola) and refraction (Cartesian ovals) and reflection (limaçon of Pascal) by a sphere [10,23]. Although the concept has been extensively used in microwaves [24], it has not received the attention it deserves in optics. We hope that this neglect will be repaired with our discussion of the FP and the notion of an anticaustic will experience a merited revival, given its beauty and potential usefulness.

2. ANTICAUSTICS IN A FABRY–PEROT INTERFEROMETER

We will examine a simplified model of a FP interferometer, which consists of a plane parallel transparent plate of refractive index $n$ and thickness $d$ surrounded by a medium of refractive index ${n^\prime}$. As heralded in the Introduction, we consider that the system is illuminated by the spherical wavefronts arising from an ideal point source ${P_0}$ at a distance $h$ of the first interface. We take a coordinate system with origin at point ${P_0}$ and with the $X$ axis normal to the interfaces, as sketched in Fig. 1. Without loss of generality, the problem will be treated as two dimensional. We will look only at the transmitted pattern, although an analogous treatment can be done for the reflected one.

 figure: Fig. 1.

Fig. 1. Scheme of the FP. The system is illuminated by a monochromatic point source ${P_0}$ located at the origin of the coordinate system. A typical ray defining the anticaustic is marked in bold blue, while the faint lines indicate other rays generated in the multiple reflections.

Download Full Size | PDF

The rays proceeding from point source ${P_0}$ form a homocentric pencil, so they constitute a normal congruence since every ray of the pencil is cut orthogonally by spheres centered on the mutual point of intersection of the rays. After being refracted, the resulting pencil (in general, no longer homocentric) will again form a normal congruence. This is the basis of the Malus–Dupin theorem [2], which guarantees the existence of wavefronts at the exit of the interferometer.

Consider a typical ray from ${P_0}$ striking the first interface at an angle of incidence ${\theta ^\prime}$ after traveling a distance $\ell$. Inside the plate, it experiences multiple reflections: after $2k + 1$ (with $k = 0,1, \ldots$) of these reflections, the ray is transmitted, as indicated in the figure. Finally, let the transmitted ray travel a distance ${\ell ^\prime}$ in the second medium (of index ${n^\prime}$).

By direct inspection, one can check that the coordinates of the endpoint $(x,y)$ of this ray can be expressed as

$$\begin{split}x &={ h + d + {\ell ^\prime}\cos {\theta ^\prime} ,}\\y&={ h\tan {\theta ^\prime} + (2k + 1)d\tan \theta + {\ell ^\prime}\sin {\theta ^\prime} ,}\end{split}$$
where $\theta$ is the angle of refraction in the medium $n$, which is related to the angle of incidence according to Snell’s law ${n^\prime}\sin {\theta ^\prime} = n\sin \theta$.

The optical path length of the considered ray, denoted by $\Delta$ and computed from source ${P_0}$, is

$$\Delta (x,y) = {n^\prime}\ell + (2k + 1)\frac{{nd}}{{\cos \theta}} + {n^\prime}{\ell ^\prime}.$$
To ensure that $(x,y)$ define a wavefront, we have to impose that $\Delta (x,y)$ takes a constant value, much in the spirit of Malus–Dupin. This gives a one-parameter family of wavefronts that constitutes the image of the family of object wavefronts centered at ${P_0}$.

Out of this family, we extract a single wavefront, the anticaustic, for which the optical path length $\Delta$ in Eq. (2) is zero: this fixes at once

$${\ell ^\prime} = - \frac{h}{{\cos {\theta ^\prime}}} - (2k + 1)\frac{{\mathfrak{n}d}}{{\cos \theta}} ,$$
where, for simplicity, we have introduced the relative index $\mathfrak{n} = n/{n^\prime}$, and we have taken into account that $\ell = h/\cos {\theta ^\prime}$. The negative sign in the optical path length indicates that we are dealing with a virtual ray, marked with broken lines in Fig. 2.
 figure: Fig. 2.

Fig. 2. First (red) and second (green) anticaustic in a FP interferometer. The broken lines represent virtual rays emanating from the anticaustics.

Download Full Size | PDF

Replacing this value of ${\ell ^\prime}$ in Eq. (1), we get

$$\begin{split}X & ={ d - (2k + 1)\mathfrak{n}d\frac{{\cos {\theta ^\prime}}}{{\cos \theta}} ,}\\Y & ={ (2k + 1)d\tan \theta - (2k + 1)\mathfrak{n}d\frac{{\sin {\theta ^\prime}}}{{\cos \theta}} ,}\end{split}$$
where we denote by $(X,Y)$ the coordinates of the anticaustic point to distinguish them from the coordinates of the endpoint of the ray.

Finally, we use Snell’s law to eliminate ${\theta ^\prime}$; the result reads

$$\begin{split} X &={ d - (2k + 1)\mathfrak{n}d\sqrt {1 + (1 - {\mathfrak{n}^2})\mathop {\tan}\nolimits^2 \theta} ,}\\ Y &={ (2k + 1)(1 - {\mathfrak{n}^2})d\tan \theta .}\end{split}$$
This is the parametric form of the anticaustic we were looking for [25]. To make things crystal clear, the explicit equation can be easily obtained by eliminating $\theta$ between these two equations:
 figure: Fig. 3.

Fig. 3. (Left) First anticaustic of the FP as in Fig. 2, showing the corresponding caustic, which is the envelope of the refracted rays. We have used $\mathfrak{n} = 1.3$ and $d = 2$ (in arbitrary units). (Right) Propagation of the previous anticaustic at four different times. We have included spherical wavefronts (dotted circles), corresponding to a single point source located at the position ${P_0}$.

Download Full Size | PDF

$$\frac{{{{(X - d)}^2}}}{{{{(2k + 1)}^2}{\mathfrak{n}^2}{d^2}}} + \frac{{{Y^2}}}{{{{(2k + 1)}^2}({\mathfrak{n}^2} - 1){d^2}}} = 1 .$$
This represents a series of ellipses for the usual case $\mathfrak{n} \gt 1$ (${n^\prime} \lt n$), whereas for $\mathfrak{n} \lt 1$ (${n^\prime} \gt n$), they are hyperbolas. These wavefronts are unique and provide exact information about every ray emanating from ${P_0}$ and refracted by the FP.

3. DISCUSSION

The action of the FP can be seen as producing the superposition of this set of ellipses (or hyperbolas) labeled by the integer $k$. In what follows, we concentrate on the case of elliptical wavefronts, although the results can be directly translated to the hyperbolic ones. All these wavefronts are concentric and have common axes of directions, the major ones being coincident with the $X$ axis. The center of all of them is at point $(d,0)$ and thus, at a distance $d$ of the source ${P_0}$. The respective semiaxes are given by

$$\begin{split}{{a_k}}&={ (2k + 1)\mathfrak{n}d ,}\\{{b_k}}&={ (2k + 1)\sqrt {{\mathfrak{n}^2} - 1} \;d .}\end{split}$$
The distance from each focus to the center is ${f_k} = \sqrt {a_k^2 - b_k^2} = (2k + 1)d$. The eccentricity of all these ellipses is the same, viz., $\varepsilon = {f_k}/{a_k} = 1/\mathfrak{n} \lt 1$. Note that these elliptical wavefronts are independent of the distance $h$, so they do not change in shape or position when the plate undergoes a translation, showing remarkable invariance.

In Fig. 2, we have plotted the first two anticaustics obtained by a simple transmission ($k = 0$) and a transmission and two internal reflections ($k = 1$). The physical meaning of these surfaces is apparent from the figure. In practice, only a certain cap is effective in each ellipse. In Fig. 2, we have marked these caps for the plotted rays. The larger the plate, the greater the extension of the cap. The limit of the effective region can be found by looking at the ray that emerges forming a greater angle with the normal to the interfaces; when virtually prolonged, it gives, on the corresponding ellipse, the point that serves as the limit.

Consider, for example, the case of $k = 0$. Suppose a spherical wavefront originates from ${P_0}$ at $t = 0$ and, after refractions, takes the form $W$ at time $t$. By propagating $W$ backward in time in the medium ${n^\prime}$ we get an initial wavefront ${W_0}$ at $t = 0$. Mutatis mutandis, the propagation of ${W_0}$ (without any refraction) yields the true wavefront $W$ at time $t$.

As stressed in the Introduction, the wavefront ${W_0}$ was investigated by Bernouilli. The term anticaustic, he concocted, has a mathematical origin: actually, the evolute of ${W_0}$ (i.e., the locus of all its centers of curvature) is precisely the caustic. Conversely, the wavefront ${W_0}$ is an involute of the caustic [26].

In consequence, the evolutes of Eq. (6) are the caustics generated in the successive internal reflections. The calculation is straightforward, with the result [27]

$${\mathfrak{n}^{2/3}}{(X - d)^{2/3}} + {({\mathfrak{n}^2} - 1)^{1/3}}{Y^{2/3}} = (2k + {1)^{2/3}} {d^{2/3}} .$$
This constitutes a series of astroids, all centered at $(d,0)$, in turn the common center of the ellipses. The amazing properties of these curves have been discussed in detail in the literature [28,29]. In Fig. 3, we plot the first anticaustic ($k = 0$) and its corresponding caustic. Although the astroid is a four-cusped curve, we concentrate on the horizontal left cusp; they appear at positions
$${P_k} = d - (2k + 1)\frac{d}{\mathfrak{n}} ,$$
as one can check directly from Eq. (8). They are thus separated by a distance $2d/\mathfrak{n}$, as it shown in Fig. 4.
 figure: Fig. 4.

Fig. 4. First three anticaustics of the same FP as in Fig. 3, with their corresponding caustics. The source ${P_0}$ (not labeled) is located at the origin of the coordinate system. In the paraxial regime, the equivalent virtual sources ${P_1},{P_2},{P_3}$ appear located at the cusps of the caustics, and they are separated by $2d/\mathfrak{n}$.

Download Full Size | PDF

Once the anticaustic ${W_0}$ is known, the propagated wavefront $W$ is equidistant from (or offset from) ${W_0}$. To obtain the equidistant curve, we consider the (oriented) normals of ${W_0}$: the equidistant curve is formed by the points at a distant $vt$ (with $v$ the velocity of light in the medium) in the forward direction from the points of the initial curve along each normal. This is represented in Fig. 3. For small $t$, the equidistant curve is smooth. But from some value of $t$ on (namely, this critical value being the minimal curvature radius of the curve), the equidistant curve acquires singularities. They were already studied by Cayley [30], and they are semiclassical cusped edges and swallowtails. They can be appreciated in the wavefronts appearing inside the caustic. In physical terms, this means that the successive images by the interfaces of point ${P_0}$ are no single points, but the system presents aberrations. This is also apparent in Fig. 4.

In Fig. 3, we have also included the spherical wavefronts associated with source ${P_0}$. We can see significant differences with the propagated anticaustics. For small aperture angles ${\theta ^\prime}$, we are in the paraxial regime, and the effective caps in the ellipses can be considered locally as spheres. In this case, the FP appears as the interference of virtual point sources ${P_1}$, ${P_2}, \ldots$, separated by a distance $2d/\mathfrak{n}$ [31]. In Fig. 4, we can appreciate that these virtual sources are precisely at the cusps of the caustics previously discussed. This makes clear contact with Young’s experiment and bears a close resemblance with a diffraction grating. There is, however, one important difference: in the FP, the amplitudes of each equivalent virtual source decrease due to the reflection and transmission coefficients in each interface.

Alternatively, these secondary sources can be seen as a far-field uniform linear array antenna, with excitation amplitudes arranged according to a geometric sequence. This configuration can be addressed with the standard methods of antenna theory [32].

Since these secondary sources get fainter farther away from the real source, only an effective number of them contribute. This number will depend on the reflectivity of the interfaces.

As a final curiosity, we note that this picture in terms of anticaustics gives an intuition about what happens when the interfaces are not strictly parallel, a conundrum in this field [2]. Now, points ${P_0}$ and the successive images, ${P_1},{P_2}, \ldots$, are spaced around a circle, and the resulting diffraction grating gives a high-order Bessel function [33], which has associated oscillations on one side of the peak.

4. CONCLUDING REMARKS

In summary, we have made extensive use of the notion of an anticaustic to provide an alternative view of the response of a FP interferometer. The intriguing mathematical properties of these surfaces [23] confirm their relevance in optical problems.

We stress that the benefit of this approach lies not in any inherent advantage in terms of efficiency in solving practical matters. Rather, we expect that the beautiful geometrical picture presented here may be helpful in updating modern views on the operation of such a relevant setup.

Funding

Ministerio de Ciencia e Innovación (PGC2018-099183-B-I00).

Acknowledgment

We acknowledge Emil Wolf for bringing our attention to the concept of anticaustic and the historical details about the names under which this notion has been badged in the literature. We thank two anonymous reviewers for their detailed and thoughtful comments that have helped us to improve the manuscript.

Disclosures

The authors declare no conflicts of interest.

Data availability

No data were generated or analyzed in the presented research.

REFERENCES

1. T. Young, “I. The Bakerian Lecture. Experiments and calculations relative to physical optics,” Philos. Trans. R. Soc. London 94, 1–16 (1804). [CrossRef]  

2. M. Born and E. Wolf, Principles of Optics (Cambridge University, 1999).

3. R. P. Feynman, R. B. Leighton, and M. L. Sands, The Feynman Lectures on Physics (Addison Wesley, 2011).

4. A. J. Fresnel, Oeuvres Complètes d’Augustin Fresnel: Théorie de la Lumière (Imprimerie Impériale, 1866).

5. H. Lloyd, “On a new case of interference of the rays of light,” Trans. R. Irish Acad. 17, 171–177 (1831).

6. F. Billet, Traité d’Optique Physique (Mallet-Bachelier Libraire, 1858).

7. G. Hernandez, Fabry-Perot Interferometers (Cambridge University, 1986).

8. J. M. Vaughan, The Fabry-Perot Interferometer: History, Theory, Practice, and Applications (Adam Hilger, 1989).

9. O. N. Stavroudis, The Optics of Rays, Wavefronts and Caustics (Academic, 1972).

10. P. E. Ouellette, “Geometric optics of a refringent sphere illuminated by a point source: caustics, wavefronts, and zero phase-fronts for every rainbow “k” order,” J. Opt. Soc. Am. A 35, 1–11 (2018). [CrossRef]  

11. J. Bernoulli, “Lineae Cycloidales, Evolutae, Anti-Evolutae, Causticae, Anti-Causticae, Peri-Causticae. Earum usus et simplex relatio ad se invicem, Spira mirabilis,” in Acta Eruditorum, Vol. 5 of Die Werke von Jakob Bernoulli (Birkhäuser, 1999), pp. 62–74.

12. A. Quetelet, “Démonstration et développement des principes fondamentaux de la théorie des caustiques secondaires,” Nouv. Mem. Acad. R. Sci. Bruxelles 5, 5–52 (1829).

13. A. Cayley, “XIII. A memoir upon caustics,” Philos. Trans. R. Soc. London 147, 273–312 (1857). [CrossRef]  

14. R. A. Herman, Treatise on Geometrical Optics (Cambridge University, 1900).

15. J. E. Eaton, “The zero phase-front in microwave optics,” IRE Trans. Antennas Propag. 1, 38–41 (1952). [CrossRef]  

16. R. Damien, Théorème sur les Surfaces d’Onde en Optique Géométrique (Gauthier-Villars, 1955).

17. O. N. Stavroudis, “Refraction of wavefronts: a special case,” J. Opt. Soc. Am. 59, 114–115 (1969). [CrossRef]  

18. S. Cornbleet, Microwave and Optical Ray Geometry (Wiley, 1984).

19. A. V. Gitin, “A simple proof of Damien’s theorem and duality in theory of the zero-distance phase front,” Int. J. Antennas Propag. 2015, 478580 (2015). [CrossRef]  

20. M. Avendaño-Alejo, L. Castañeda, and I. Moreno, “Exact wavefronts and caustic surfaces produced by planar ripple lenses,” Opt. Express 23, 21637–21649 (2015). [CrossRef]  

21. J. C. Chastang and R. T. Farouki, “The mathematical evolution of wavefronts,” Opt. Photon. News 3(1), 20–23 (1992). [CrossRef]  

22. E. Wolf (personal communication, 2006).

23. R. T. Farouki, “The Cartesian ovals,” Math. Intell. 42, 1–11 (2022). [CrossRef]  

24. S. Cornbleet, Microwave and Geometrical Optics (Academic, 1994).

25. M. Martínez-Risco, Oeuvres Complètes (Presses Universitaires de France, 1976).

26. R. T. Farouki, H. P. Moon, and B. Ravani, “Minkowski geometric algebra of complex sets,” Geom. Dedicata 85, 283–315 (2001). [CrossRef]  

27. A. Gray, Modern Differential Geometry of Curves and Surfaces (CRC Press, 1993).

28. E. H. Lockwood, A Book of Curves (Cambridge University, 1963).

29. R. Yates, A Handbook on Curves and Their Properties (Literary Licensing, 2012).

30. V. I. Arnold, The Theory of Singularities and Its Applications (Cambridge University, 1991).

31. L. L. Sánchez-Soto, J. J. Monzón, and G. Leuchs, “The many facets of the Fabry–Perot,” Eur. J. Phys. 37, 064001 (2016). [CrossRef]  

32. C. A. Balanis, Antenna Theory: Analysis and Design, 4th ed. (Wiley, 2016).

33. D. Mcgloin and K. Dholakia, “Bessel beams: diffraction in a new light,” Contemp. Phys. 46, 15–28 (2005). [CrossRef]  

Data availability

No data were generated or analyzed in the presented research.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1.
Fig. 1. Scheme of the FP. The system is illuminated by a monochromatic point source ${P_0}$ located at the origin of the coordinate system. A typical ray defining the anticaustic is marked in bold blue, while the faint lines indicate other rays generated in the multiple reflections.
Fig. 2.
Fig. 2. First (red) and second (green) anticaustic in a FP interferometer. The broken lines represent virtual rays emanating from the anticaustics.
Fig. 3.
Fig. 3. (Left) First anticaustic of the FP as in Fig. 2, showing the corresponding caustic, which is the envelope of the refracted rays. We have used $\mathfrak{n} = 1.3$ and $d = 2$ (in arbitrary units). (Right) Propagation of the previous anticaustic at four different times. We have included spherical wavefronts (dotted circles), corresponding to a single point source located at the position ${P_0}$.
Fig. 4.
Fig. 4. First three anticaustics of the same FP as in Fig. 3, with their corresponding caustics. The source ${P_0}$ (not labeled) is located at the origin of the coordinate system. In the paraxial regime, the equivalent virtual sources ${P_1},{P_2},{P_3}$ appear located at the cusps of the caustics, and they are separated by $2d/\mathfrak{n}$.

Equations (9)

Equations on this page are rendered with MathJax. Learn more.

x=h+d+cosθ,y=htanθ+(2k+1)dtanθ+sinθ,
Δ(x,y)=n+(2k+1)ndcosθ+n.
=hcosθ(2k+1)ndcosθ,
X=d(2k+1)ndcosθcosθ,Y=(2k+1)dtanθ(2k+1)ndsinθcosθ,
X=d(2k+1)nd1+(1n2)tan2θ,Y=(2k+1)(1n2)dtanθ.
(Xd)2(2k+1)2n2d2+Y2(2k+1)2(n21)d2=1.
ak=(2k+1)nd,bk=(2k+1)n21d.
n2/3(Xd)2/3+(n21)1/3Y2/3=(2k+1)2/3d2/3.
Pk=d(2k+1)dn,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.