We systematically characterize the Fabry-Pérot resonator. We derive the generic Airy distribution of a Fabry-Pérot resonator, which equals the internal resonance enhancement factor, and show that all related Airy distributions are obtained by simple scaling factors. We analyze the textbook approaches to the Fabry-Pérot resonator and point out various misconceptions. We verify that the sum of the mode profiles of all longitudinal modes is the fundamental physical function that characterizes the Fabry-Pérot resonator and generates the Airy distribution. Consequently, the resonator losses are quantified by the linewidths of the underlying Lorentzian lines and not by the measured Airy linewidth. Therefore, we introduce the Lorentzian finesse which provides the spectral resolution of the Lorentzian lines, whereas the usually considered Airy finesse only quantifies the performance of the Fabry-Pérot resonator as a scanning spectrometer. We also point out that the concepts of linewidth and finesse of the Airy distribution of a Fabry-Pérot resonator break down at low reflectivity. Furthermore, we show that a Fabry-Pérot resonator has no cut-off resonance wavelength. Finally, we investigate the influence of frequency-dependent mirror reflectivities, allowing for the direct calculation of its deformed mode profiles.
© 2016 Optical Society of America
The Fabry-Pérot resonator which was invented in 1899  has proven very useful as a high-finesse interferometer in uncountable spectroscopic applications. Since 1960, it has also formed the fundamental basis for a large class of open resonators that have enabled laser oscillation. Lasers with resonator configurations that provide strong light confinement and accordingly high optical gain, such as semiconductor and dielectric channel waveguide lasers, as well as fiber lasers, are particularly successful examples. Due to their high gain, these configurations do not only accommodate high propagation losses and allow for high transmission losses through one of their resonator mirrors, but are in fact often preferably operated with high transmission, because the outcoupling efficiency, i.e., the fraction of useful transmission losses among the overall resonator losses, is one of the key parameters determining the laser slope efficiency. Driven by the quest for miniaturization and mass fabrication, in recent years also various integrated-optical resonator structures with low quality factor Q, be it because of low feedback or high intrinsic losses or both, have emerged. Therefore, low-reflectivity and high-loss Fabry-Pérot resonators have become very common tools nowadays.
The Fabry-Pérot resonator has been extensively investigated, experimentally as well as theoretically. Its theoretical descriptions and proposed approximations and simplifications stem from times when high-finesse spectroscopic applications were dominating and lasers either did not exist or still provided rather low gain, hence necessitating low-loss resonators. However, when the losses become high, discrepancies between the theoretical approaches surface, common approximations turn invalid, and even the definitions of typical parameters break down or prove inappropriate. All these deficiencies go unnoticed when the resonator losses are low. In fact, the class of low-reflectivity and high-loss Fabry-Pérot resonators is not well described by the approximations applied in the common literature.
In this paper, we analyze the performance and relevant parameters of the Fabry-Pérot resonator, with reference to the book on the Fabry-Pérot resonator by Vaughan [2,3] and the standard text books by Siegman [4–7], Svelto [8,9], and Saleh and Teich [10–12]. We demonstrate that the sum of the mode profiles of all longitudinal modes, which are Lorentzian-shaped in case of frequency-independent losses, is the fundamental spectral function that characterizes the Fabry-Pérot resonator. It physically generates the Airy distribution of the Fabry-Pérot resonator, even at low reflectivity and for different or frequency-dependent reflectivities of the two mirrors. Based on this fact, we emphasize that the resonator losses are quantified by the linewidths of the underlying Lorentzian lines and not by the measured Airy linewidth, and we define the Lorentzian finesse which quantifies the resolution of the Lorentzian lines, whereas the usually considered Airy finesse only quantifies the performance of the Fabry-Pérot resonator as a scanning spectrometer. We also point out that the concepts of linewidth and finesse of the Airy distribution of a Fabry-Pérot resonator break down at a mirror-reflectivity product equaling (0.172)2. Moreover, we investigate the influence of frequency-dependent mirror reflectivity, resulting in deformed mode profiles that can be analytically determined.
2. Resonator losses and outcoupled light
2.1 Definition of resonator losses
Throughout this paper, we assume a two-mirror Fabry-Pérot resonator of geometrical length ℓ, homogeneously filled with a medium of refractive index nr. Both, ℓ and nr are assumed to vary insignificantly over the frequency range of interest. The round-trip time tRT of light travelling in the resonator with speed c = c0/nr, where c0 is the speed of light in vacuum, is given by4] and others named the ratios δout,i(ν) the “mirror coupling coefficients”, referred to as the “delta notation”.
Other discrete losses are, e.g., diffraction losses at the mirrors of finite lateral dimensions, absorption and scattering losses in the mirror coatings, diffraction losses at the edges of an active laser medium that is shorter than the resonator length, or losses due to other optical elements inside the resonator. Furthermore, continuous intensity losses per unit length may be induced by the medium, e.g. scattering losses due to material imperfections, waveguide propagation losses due to interface roughness, or absorption losses. All these losses shall be neglected in this paper, because they do not change the fundamental findings presented here. The photon-decay time τc(ν) of the resonator is then given by
In the two subsequent Sections 3 and 4, the outcoupling losses are assumed to be independent of frequency. The case in which τc varies with frequency due to a frequency dependence of the mirror reflectivities Ri is treated in Section 5.
2.2 Photon decay
The resonator can be investigated by inserting light as a short pulse or continuously. In the transient regime, the number φ(t) of photons at frequency ν, present inside the resonator at time t, is described via the photon-decay time τc(ν) by the differential rate equation
In contrast, a steady-state or continuous-wave (cw) regime where the decaying photon number is continuously replenished by the incoupling of photons from an external source, such that φ(t) = φs at all times, is described by the equation
2.3 Resonance frequencies, free spectral range, and spectral line shapes
With ϕ(ν) quantifying the single-pass phase shift between the mirrors, the round-trip phase shift at frequency ν accumulates to3]
According to Eq. (5), light at frequency ν oscillating inside the resonator decays out of the resonator with a time constant of τc(ν). If the resonator losses are independent of frequency, the photon-decay time τc(ν) is the same at all frequencies. The decaying electric field at frequency νq is represented by a damped harmonic oscillation with an initial amplitude of Eq,s and a decay-time constant of 2τc. In phasor notation, it can be expressed as13]Eq. (15) compensates for the broadening of the Lorentzian line as the photon-decay time becomes shorter, such that the integral of the Lorentzian line is normalized to unity. Introducing the full-width-at-half-maximum (FWHM) linewidth Δνc of the Lorentzian spectral line shape, we obtain
When repeating the above Fourier transformation for all the modes with mode index q in the resonator, one obtains the full mode spectrum of the resonator.
3. Airy distributions of the Fabry-Pérot resonator
In this Section, we demystify the Airy distribution. We derive the generic Airy distribution, show that it represents the spectral dependence of the internal resonant enhancement factor, and demonstrate that all other Airy distributions describing the circulating, back-circulating, transmitted, and back-transmitted intensity (Fig. 1) only include simple scaling factors, depending on whether the light incident upon mirror 1 or its fraction launched into the resonator is considered as a reference. Most importantly, we verify that the physical origin of the Airy distribution is the sum of mode profiles of the longitudinal resonator modes.
3.1 Generic Airy distribution: The internal resonance enhancement factor
The response of the Fabry-Pérot resonator is most easily derived by use of the circulating-field approach , as displayed in Fig. 1. This approach assumes a steady state and derives the Airy distributions via the electric field Ecirc circulating inside the resonator. In fact, Ecirc is the field propagating in the forward direction from mirror 1 to mirror 2 after interference between the field ERT that is circulating after one round trip, i.e., after having suffered outcoupling losses at both mirrors, and the field Elaun launched through the first mirror.Eqs. (56)-(57) of Appendix A yieldsFig. 2. At the resonance frequencies νq, where sin(ϕ) equals zero, the internal resonance enhancement factor is
The internal resonance enhancement factor Acirc integrated over one free spectral range, i.e., 2ϕ in the interval [−π, … π] or ϕ in the interval [−π/2, … π/2], is (see Appendix B)Fig. 3(a).Fig. 3(b). The spectral range of enhancement decreases from ΔνFSR/2 for R1R2 = 0 to ΔνFSR/3 for R1R2 → 1.
3.2 Other Airy distributions
Once the internal resonance enhancement, the generic Airy distribution of Eq. (19), is established, all other Airy distributions, i.e., the observed light intensities relative to the launched or initial intensity from the light source (Fig. 1) can then be straight-forwardly deduced by simple scaling factors. Since the intensity launched into the resonator equals the transmitted fraction of the intensity incident upon mirror 1,
The back-transmitted intensity Iback cannot be measured, because also the initially back-reflected light adds to the backward-propagating signal. The measurable case of the intensity resulting from the interference of both backward-propagating electric fields is treated in Appendix A. At the resonance frequencies, the back-emitted electric field Eback destructively interferes with the electric field Erefl,1 initially back-reflected.
Like Acirc represents the internal resonance enhancement factor with respect to Ilaun, of Eq. (29) represents the external resonance enhancement factor with respect to Iinc. The external resonance enhancement factor at the resonance frequencies, , and integrated over one free spectral range, , respectively, areEq. (23).
In experimental situations, often light is transmitted through a Fabry-Pérot resonator in order to characterize the resonator or to use it as a scanning interferometer (see later in Section 4). Therefore, an often applied Airy distribution is of Eq. (31). It describes the fraction Itrans of the intensity Iinc of a light source incident upon mirror 1 that is transmitted through mirror 2, see Fig. 1. is displayed in Fig. 4 (solid lines) for different values of the reflectivities R1 = R2. Its peak value at the resonance frequencies νq is5] by considering an additional phase shift of during each transmission through a mirror (alternatively, but equivalently, often a phase shift of during reflection from the low-refractive-index side of a mirror is assumed)  and by use of Eqs. (56)-(58) of Appendix A:Eq. (31). Alternatively, it can be obtained via the round-trip-decay approach  by tracing the infinite number of round trips that the incident electric field exhibits after entering the resonator and accumulating the electric field transmitted in all round trips. The field transmitted after the first propagation and the smaller and smaller fields transmitted after each consecutive propagation through the resonator areEq. (37), therefore one obtains the same Airy distribution of Eq. (38). Svelto  uses the alternative convention for the phase shift occurring at the mirrors, hence obtains no minus sign in the numerator of Eq. (40), which becomes irrelevant when squaring the electric field to calculate the intensity. has also been derived in the impulse-response approach . This approach is identical to the round-trip-decay approach , except that it assumes a specific light source that has a delta function in time, which Fourier transforms to an infinite frequency spectrum, such that is probed over the entire spectral range.
3.3 Airy distribution as a sum of mode profiles
In Fig. 4 one observes that, at high reflectivity, there is almost perfect agreement between the spectral shape of the Airy distribution (solid purple line) and its underlying Lorentzian lines (dashed purple line), i.e., the former is rather well represented by the latter. This fact has prompted Saleh and Teich  to propose that in this case the Airy linewidth ΔνAiry of a Fabry-Pérot resonator is similar to the linewidth Δνc = 1/(2πτc) of its underlying Lorentzian lines, both defined as FWHM (black line). However, as is generally well known, with decreasing reflectivity the linewidth of the Airy distribution (solid lines) broadens faster than that of the underlying Lorentzian lines (dashed lines).
Svelto  attributes this discrepancy to Eq. (5) being only an approximation, thereby implicating that also Eq. (15) is only an approximation, such that the Airy linewidth ΔνAiry of a Fabry-Pérot resonator can only at high reflectivity be approximated by the linewidth Δνc of its underlying Lorentzian lines. We will demonstrate that the discrepancy has a different reason.
According to Koppelmann , Bayer-Helms  “showed that the Airy distribution can be represented exactly” by the sum of Lorentzian spectral line shapes times a calibration factor. Firstly, while being literally correct, this statement is physically misleading, secondly, the calibration factor used in  remains unexplained, thirdly, the equivalence is shown only for equal reflectivities, R1 = R2, and finally, the equivalence is not investigated for non-Lorentzian spectral line shapes. Kumer and Uplinger  compared the sum of Lorentzian spectral line shapes to an Airy distribution that is different from Eq. (19) and Eq. (31) and, not surprisingly, these authors found a difference.
Here we verify that the Airy distribution is nothing else but the sum of the mode profiles of the longitudinal resonator modes, thereby revealing the physical origin of the Airy distribution. Our approach starts from the electric field Ecirc circulating inside the resonator, considers the exponential decay in time of this field through both mirrors of the resonator, see Fig. 1, Fourier transforms it to frequency space according to Eq. (14) to obtain the normalized spectral line shapes of Eq. (15), divides it by the round-trip time tRT to account for how the total circulating electric-field intensity is longitudinally distributed in the resonator and coupled out per unit time, resulting in the emitted mode profiles,Eq. (28). Exploiting the derivation given in Eqs. (64)-(72) of Appendix C, the sum of emitted mode profiles of Eq. (41) yieldsEq. (28), with Eq. (19) inserted. Each spectral line shape of Eq. (15) and mode profile of Eq. (41) is extended over the infinite frequency range, consequently light at a specific frequency ν inside the resonator excites all longitudinal modes of the resonator. However, the contributions from different longitudinal modes to the light at frequency ν do not interfere with each other, because all optical modes are orthogonal with each other. For this reason, the sum in Eq. (42) is over the intensity mode profiles rather than over the electric fields. Because of resonant enhancement of the launched light the peak value of Aemit at the resonance frequencies νq,Eqs. (23) and (59) the energy of the system is conserved.
The observation that with decreasing R1 and R2 the linewidth of the resulting Airy distribution in Fig. 4 is increasingly broader than the linewidth of the underlying Lorentzian lines simply arises from the fact that one sums up mode profiles (with the same linewidth as the Lorentzian lines) that resonate at different frequencies. It does not constitute a discrepancy, as has often been proposed.
The derivation shown here demonstrates that―from a physical point of view―the spectral line shapes and mode profiles are the fundamental spectral functions that characterize the Fabry-Pérot resonator and their sum quantifies its spectral response. The circulating-field and round-trip-decay approaches only implicitly exploit this fact, thereby obscuring the physical nature of the Airy distribution. As we will see in Section 4, this fundamental understanding has direct consequences for the definitions of linewidth and finesse.
The same simple scaling factors of Eqs. (23) and (24) that provide the relations between the individual Airy distributions, see Eqs. (19) and (25)-(33), also provide the relations among and the other mode profiles:
3.4 0th-order mode: no cut-off wavelength!
The resonance frequencies of the Fabry-Pérot resonator are given by Eq. (12). Particularly, q = 0, resulting in νq = 0, is a physical solution of Eq. (12), and the mode profiles summed in Eq. (42) include this 0th-order mode . Consequently, in wavelength space there exists no cut-off wavelength at the resonance line with q = 1, but an additional resonance occurs at infinite wavelength. As long as the resonator losses are low (purple line in Fig. 5), in most cases this resonance will go unnoticed. However, when the losses become significant and the mode profiles broaden, this resonance manifests itself in a steady increase of the Airy distribution in the long-wavelength range (other lines in Fig. 5). While this effect looks quite dramatic in Fig. 5, one should keep in mind that the natural energy or frequency units are, when converted to wavelength units, expanded by λ2 at the long-wavelength side.
4. Lorentzian linewidth and finesse versus Airy linewidth and finesse
The most commonly accepted definitions of spectral resolution are the Rayleigh criterion and the Taylor criterion. The Rayleigh criterion  applies to diffraction spots and proposes that two intensity maxima can be resolved if the maximum of the first spot is spectrally superimposed on the first minimum of the second spot. For spectral lines that do not exhibit minima, such as the Lorentzian lines underlying the Airy distribution of a Fabry-Pérot resonator, as well as the Airy distribution itself within one free spectral range from νq − ΔνFSR/2 to νq + ΔνFSR/2, the Taylor criterion  is more appropriate. It proposes that two spectral lines can be resolved if the individual lines cross at half intensity. In the case of two identical, symmetric spectral lines their peaks would then be separated by their FWHM. The Taylor criterion is utilized in the following.
4.1 Characterizing the Fabry-Pérot resonator: Lorentzian linewidth and finesse
When launching light into the Fabry-Pérot resonator in a non-scanning experiment, i.e., at fixed resonator length (and fixed angle of incidence), the Fabry-Pérot resonator becomes the object of investigation. The spectral line shapes, Lorentzian lines, and mode profiles are the fundamental functions. By measuring the sum of mode profiles, the Airy distribution, one can derive the total loss of the Fabry-Pérot resonator via recalculating the Lorentzian linewidth Δνc of Eq. (16), displayed (blue line) relative to the free spectral range in Figs. 6(a) and 6(c).
The underlying Lorentzian lines can be resolved as long as the Taylor criterion is obeyed (Fig. 7). Consequently, one can define a parameter which we call the Lorentzian finesse of a Fabry-Pérot resonator:Figs. 6(b) and 6(d). This parameter has received little attention in the scientific literature, other than that at high reflectivity, where ΔνAiry ≈Δνc, see the purple line in Fig. 4, it can serve as a reasonable approximation of the Airy finesse FAiry (see Section 4.2). Yet, the Lorentzian finesse Fc has a fundamental physical meaning: it describes how well the Lorentzian lines underlying the Airy distribution can be resolved when measuring the Airy distribution. A Fabry-Pérot resonator generating single-longitudinal-mode laser light is characterized by its Lorentzian linewidth and finesse.
Since Δνc exists for any mirror reflectivity, the definition of the Lorentzian finesse does not break down at a critical value. However, at the point whereEq. (45) plays an essential role in the characterization of low-reflectivity or otherwise high-loss Fabry-Pérot resonators.
4.2 Scanning the Fabry-Pérot resonator: Airy linewidth and finesse
A different situation occurs when the Fabry-Pérot resonator is used as a scanning interferometer, i.e., at varying resonator length (or angle of incidence), to spectroscopically distinguish spectral lines at different frequencies within one free spectral range. In this case several Airy distributions, each one created by an individual spectral line, must be resolved. Therefore, now the Airy distribution becomes the underlying fundamental function and the measurement delivers a sum of Airy distributions. The parameters that properly quantify this situation are the Airy linewidth ΔνAiry and the Airy finesse FAiry.
On either side of the peak of Eq. (35) located at νq, the transmitted intensity decreases to when the phase shift ϕ changes by the amount Δϕ and, accordingly, sin2(ϕ) changes from 0, such that in the denominator of Acirc of Eq. (19)Eqs. (10) and (11) to calculate ϕ = πν/ΔνFSR, resulting in Δϕ = π(ΔνAiry/2)/ΔνFSR, then provides the FWHM linewidth ΔνAiry of the Airy distribution ,Figs. 6(a) and 6(c) in direct comparison with the Lorentzian linewidth Δνc. The ratio between ΔνAiry of Eq. (48) and Δνc of Eq. (16) is displayed in Fig. 6(e).
The concept of defining the linewidth of the Airy peaks as FWHM breaks down at ΔνAiry = ΔνFSR (solid red line in Fig. 4), because at this point the Airy linewidth instantaneously jumps to an infinite value. For lower reflectivity values R1R2, the FWHM linewidth of the Airy peaks is undefined and other definitions or concepts would have to be utilized to describe the situation. The limiting case occurs atFig. 4). A number of fiber and channel waveguide lasers that produce very high gain due to their strong pump confinement are operated in resonators that rely solely on the Fresnel reflection on both ends of the resonator, which is typically below this value. I.e., we are operating numerous laser resonators whose passive Airy linewidth ΔνAiry is undefined. Nevertheless, no problem arises from it, because the appropriate parameters in the situation of the lasing Fabry-Pérot resonator are its Lorentzian linewidth Δνc and finesse Fc.
Definitions can never be right or wrong. They can only be more or less meaningful and, in the worst case, contradict with another more or less meaningful definition. The finesse of the Airy distribution of a Fabry-Pérot resonator, which is displayed as the green curve in Figs. 6(b) and 6(d) in direct comparison with the Lorentzian finesse Fc, is properly defined asFig. 8). Consequently, this definition of the Airy finesse is consistent with the Taylor criterion of the resolution of a spectrometer and is, therefore, denoted by an exclamation mark in Eq. (50). Since the concept of the FWHM linewidth breaks down at ΔνAiry = ΔνFSR, consequently the Airy finesse is defined only until FAiry = 1, see Fig. 6(d), because the arcsin function in Eq. (50) cannot produce values above π/2.
Generally, one can argue that the Taylor criterion leaves some ambiguity to the definition of the limit of spectral resolution, because it does not state whether it requires the absence or allows for the presence of additional resolvable spectral lines. In the example of Fig. 8, the sum of two adjacent “resolvable” Airy distributions (dashed gray line) is better resolvable than the sum of all “resolvable” Airy distributions (dashed black line), because in the latter case additional Airy distributions contribute, see the dashed lines within the red circle in Fig. 8. If either of the two options defines the limit of spectral resolution, the other cannot.
4.3 Questionable approximations and definitionsFigs. 6(a) and 6(c), deviates from the correct curve at low reflectivities and incorrectly does not break down when ΔνAiry > ΔνFSR. This approximation is then typically also inserted into Eq. (50) to calculate the Airy finesse , resulting in2] and Siegman  even defined the Airy finesse by its approximation of Eq. (52),Eq. (53) does not comply with the Taylor criterion, it is denoted by a question mark. Saleh and Teich [12,11] also proposed Eq. (53) for the Airy finesse, but from their derivation it remains unclear whether they consider it as a definition or as an approximation of Eq. (50).
Whereas the definitions of Airy linewidth and finesse break down at ΔνAiry = ΔνFSR and FAiry = 1, the Lorentzian linewidth Δνc and finesse Fc do not break down at a critical value. Moreover, the latter reach their Taylor limit at higher losses than the former. Therefore, useful and well defined information can be extracted from non-scanning Fabry-Pérot measurements beyond the limit of the Airy linewidth and finesse (which are irrelevant quantities for non-scanning measurements).
5. Response to frequency-dependent reflectivity
In the previous Sections we described a Fabry-Pérot resonator whose mirror reflectivities are independent of frequency. We showed that the Airy distribution is nothing else but the sum of its underlying mode profiles. Now we consider the mirror reflectivities as general functions of frequency, Ri(ν), such that the photon-decay time τc(ν) of Eq. (3) becomes a function of frequency. As a result, the spectral line shapes of Eq. (15), the mode profiles of Eq. (41) and all other mode profiles, as well as the Airy distribution Aemit(ν) of Eq. (42) and all other Airy distributions, are spectrally modified. In Eqs. (16) and (17), Δνc turns into a local function of frequency, thereby losing its meaning as the FWHM linewidth. Nevertheless, all Airy distributions fundamentally remain the sums of their corresponding mode profiles. Consequently, even for frequency-dependent reflectivities Ri(ν) we can calculate the spectral line shape and mode profile of each mode directly from Eqs. (15) and (41), respectively, and obtain all other mode profiles via the simple scaling factors of Eq. (44). These spectral line shapes and mode profiles are the most general expressions for the longitudinal modes of a Fabry-Pérot resonator with frequency-dependent reflectivities.
In order to illustrate the deformation of the mode profiles and their resulting Airy distributions in a resonator with frequency-dependent reflectivities, we consider two different resonators. Both have air (nr = 1) as their internal medium and identical multilayer mirrors with the reflectivity curve shown in both, Figs. 9(a) and 9(b), which is symmetric around its maximum value of 98% at 300 THz (vertical black dashed lines in Fig. 9). These resonators are distinguished only by their lengths of ℓ1 = 999.808 µm (left-hand side in Fig. 9) and ℓ2 = 999.924 µm (right-hand side in Fig. 9). Their resonance frequencies (vertical dashed lines) with mode indexes q = 2000 (red lines), q = 2001 (green lines), and q = 2002 (blue lines) occur at νq = 299.85 THz, 300.00 THz, and 300.15 THz (left-hand side of Fig. 9) and νq = 299.815 THz, 299.965 THz, and 300.115 THz (right-hand side of Fig. 9). Consequently, the nearest resonance frequency with mode index q = 2001 coincides with the reflectivity maximum at 300 THz, leading to a symmetric situation (left-hand side), or is spectrally offset by 35 GHz, leading to an asymmetric situation (right-hand side). In order to account only for the effects of reflectivity, we neglect the penetration depth of the light inside the mirrors, which can, in principle, also be frequency dependent.
In Figs. 9(c) and 9(d), the mode profiles γq,emit of the three nearest modes with mode indexes q = 2000, q = 2001, and q = 2002 (red, green, and blue solid lines), the numerical summation of 6 × 106 mode profiles, which were summed at various frequencies (magenta dots), and the Airy distribution Aemit (black solid line) are displayed. In both situations, the sum of mode profiles and the Airy distribution show excellent agreement (and would show perfect agreement if it were possible to extend the summation over the infinite sum of mode profiles). In the symmetric situation of Fig. 9(c), the main peak of the central mode profile (green solid line) as well as the Airy distribution (black solid line) reach a peak value of ∼99, but are only shown until the value of 4 to better resolve the details of the mode profiles. In the asymmetric situation of Fig. 9(d), the mode profile of the nearest mode with q = 2001 (green solid line) and the Airy distribution (black solid line) become strongly asymmetric. The main peak of the nearest mode profile as well as the Airy distribution lies between the resonance frequency (green dashed line) and the reflectivity maximum (black dashed line). The main peak of the mode profile has a FWHM linewidth of 7.2 GHz that is significantly narrower than the FWHM linewidth of 118.7 GHz that would occur if the photon-decay time τc were constant with frequency and determined by the reflectivity value of 8.3% at the mode’s resonance frequency of ν2001 = 299.965 THz. The spectral dependence of the mode profiles γq,emit is determined by the internal resonance enhancement factor Acirc and the simple scaling factor 1−R1(ν)R2(ν) of Eqs. (28) or (44), which both depend on frequency. In particular, the internal resonant enhancement factor Acirc depends on frequency through the phase ϕ(ν) and through the reflectivity curves R1(ν) and R2(ν) of the two resonator mirrors. Acirc (black solid line) is shown in Figs. 9(e) and 9(f). The peak value (not displayed in the figures) is 2500 in the symmetric situation and 7.56 in the asymmetric situation. Also shown is the Airy distribution Atrans (magenta solid line), whose peak value reaches 50 and 2.2 in the symmetric and asymmetric situation, respectively. The spectral contents of Aemit, Acirc, and Atrans and the frequencies of their main peaks differ from each other, because the scaling factors between these Airy distributions are frequency dependent. Besides Aemit, the distributions Acirc and Atrans are only two of the other Airy distributions related to the launched intensity Ilaun. We leave it to the reader to calculate Ab-circ and Aback and the mode profiles of all these Airy distributions and verify that, in all cases, their sums generate the corresponding Airy distributions.
Among the Airy distributions related to the incident intensity Iinc, the most utilized one, (black solid line), is shown in Figs. 9(g) and 9(h), together with the three nearest mode profiles (red, green, and blue solid lines) and the sum of 6 × 106 mode profiles (magenta dots). At the resonance frequencies becomes unity. In addition, whenever the reflectivity curve in Figs. 9(a) and 9(b) approaches zero, approaches unity, but obviously these peaks do not represent resonances. In Fig. 9(g), the resonance with mode index q = 2001 is located in the center of the reflection band and identified as a resonance line, whereas in Fig. 9(h) it almost coincides with a frequency of zero reflection and, thus, does not appear as a separate peak. Like for , those non-resonant peaks occur for all other Airy distributions. The Airy distributions (black solid line) and (magenta solid line) are displayed in Figs. 9(i) and 9(j). In the symmetric situation, extends to 50 and to 1.98, whereas in the asymmetric situation extends to 2.2 and to 1.42.
In all cases, the mode profiles are spectrally modified by the reflectivity curve. They obtain large values around their resonance frequencies νq, but their spectral center of mass usually does not coincide with their resonance frequency.
One can numerically sum a large number M = 2m + 1, for q ranging from –m to m, of mode profiles and reproduce the corresponding Airy distribution. Examples are shown in Figs. 9(c) and 9(d) for the mode profiles γq,emit and the Airy distribution Aemit and in Figs. 9(g) and 9(h) for the mode profiles and the Airy distribution . The error ϵ decreases towards zero as M increases to infinity. Empirically we find that for sufficiently large values of M, such that the product m × ΔνFSR is larger than the frequency at which the error is evaluated, the error ϵ follows a monomial function ϵ(M) = εM −1. When the contributions of those modes resonating at frequencies that are very distant from the frequency under investigation become too small for the numerical addition to work properly, eventually one runs into numerical errors. These can be reduced by using variable precision arithmetic or by special summation algorithms such the Kahan summation .
Now we assume, firstly, that the mirror reflectivities change slowly with frequency, such that they are almost constant over each free spectral range, Ri(νq − ΔνFSR/2) ≈Ri(νq) ≈Ri(νq + ΔνFSR/2). Secondly, we assume that the spectral line shape of each mode is mostly confined to a free spectral range, i.e., its linewidth is Δνq < ΔνFSR. These two assumptions would provide a justification to approximate the mode profiles by Lorentzian-shaped mode profiles. However, these Lorentzian-shaped mode profiles would be different from each other, Δνq1 ≠ Δνq2, because Ri(νq1) ≠ Ri(νq2). When summing these Lorentzian-shaped mode profiles, the error we would make in the evaluation of the Fabry-Pérot transmission would never decrease to zero with increasing number of summation terms, but would converge to a constant value. Furthermore, for specific frequencies the error may also assume negative values which may, in specific situations, lead to the unrealistic outcome of the relative transmission exceeding unity. We verified numerically that the proposed approximation would lead to errors below 5 × 10−6 only for slowly varying reflectivities (with respect to the free spectral range) in the range of 0.02% variation per THz, and for large reflectivity values near 80% reflectivity, such that the spectral mode shape is mostly contained within a free spectral range. For reflectivity variations of ~0.085% (or 0.02%) per THz and minimum reflectivities reaching 17.2%, the error stays below 1 × 10−4 (or 1 × 10−5).
As we have demonstrated in this paper, the understanding that the Airy distribution describing the spectral transmission of a Fabry-Pérot resonator physically originates in the sum of mode profiles of the longitudinal resonator modes has fundamental consequences. The resonator losses are related to the linewidth of the Lorentzian lines rather than the linewidth of the Airy distribution. Hence, a new parameter, the Lorentzian finesse, becomes important. Whereas the concept of Airy linewidth and finesse breaks down at low reflectivity, the Lorentzian linewidth and finesse are still well defined, thus allowing for a correct characterization of the Fabry-Pérot resonator. Once the internal resonance enhancement, equaling the generic Airy distribution that characterizes the light intensity that is forward circulating inside the resonator, is known, all other Airy distributions of back-circulating, transmitted, back-transmitted, and total emitted light intensity can be derived by simple scaling factors. Furthermore, in the case of frequency-dependent mirror reflectivities, the deformed spectral line shapes and mode profiles of the underlying modes can be derived from the same simple equations. Also in this generalized situation, each sum of mode profiles generates the corresponding Airy distribution.
Especially the characterization of Fabry-Pérot resonators with low reflectivity, high loss, or strongly frequency-dependent mirror reflectivities, such as distributed-Bragg-reflector and distributed-feedback resonators as well as any type of micro-resonator, which nowadays are routinely applied in high-gain lasers and nano-optics, respectively, will profit from the fundamental insight provided in this paper.
A. Derivation of Airy distribution of reflected intensity
In the situation of Fig. 1, the back-transmitted electric field is
B. Calculation of the integral resonant enhancement
The integral of the internal resonant enhancement factor Acirc of Eq. (19) can be written as
C. Derivation of the sum of mode profiles
With the definitionsEq. (2) and assumed a resonator that exhibits only outcoupling losses, we can write the spectral line shapes (normalized to unit integral) of Eq. (15) in terms of δout, which can be a general function of frequency, as follows22]Eq. (64) into Eq. (70) yieldsEq. (42).
All authors acknowledge financial support by the ERC Advanced Grant “Optical Ultra-Sensor” No. 341206 from the European Research Council.
References and links
1. C. Fabry and A. Pérot, “Théorie et applications d’une nouvelle méthode de spectroscopie interférentielle,” Ann. de Chim. et de Phys. 16(7), 115–144 (1899).
2. J. M. Vaughan, The Fabry-Pérot Interferometer (Bristol and Philadelphia, 1989), Ch. 3.2.3., pp. 97−102.
3. J. M. Vaughan, The Fabry-Pérot Interferometer (Bristol and Philadelphia, 1989), Ch. 3.3.2., p. 105.
4. A. E. Siegman, Lasers (University Science Books, 1986), Ch. 11.4, pp. 428−430.
5. A. E. Siegman, Lasers (University Science Books, 1986), Ch. 11.3, pp. 413–428.
6. A. E. Siegman, Lasers (University Science Books, 1986), Ch. 11.5, pp. 432–440.
7. A. E. Siegman, Lasers (University Science Books, 1986), Ch. 11.1, pp. 398–408.
8. O. Svelto, Principles of Lasers, 5th ed. (Springer, 2010), Ch. 4.5.1, pp. 142−146.
9. O. Svelto, Principles of Lasers, 5th ed. (Springer, 2010), Ch. 5.3, pp. 169−171.
10. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, 2nd ed. (Wiley-Interscience, 2007), pp. 571−572.
11. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, 2nd ed. (Wiley-Interscience, 2007), Ch. 7.1B, pp. 254−257.
12. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, 2nd ed. (Wiley-Interscience, 2007), Ch. 2.5B, pp. 62−66.
13. M. Eichhorn and M. Pollnau, “Spectroscopic foundations of lasers: Spontaneous emission into a resonator mode,” IEEE J. Sel. Top. Quantum Electron. 21(1), 486 (2015). [CrossRef]
14. J. O. Stoner, “PEPSIOS purely interferometric high-resolution scanning spectrometer. III. Calculation of interferometer characteristics by a method of optical transients,” J. Opt. Soc. Am. 56(3), 370–376 (1966). [CrossRef]
15. G. Koppelmann, “Multiple-beam interference and natural modes in open resonators,” in Progress in Optics, Vol. 7, E. Wolf ed. (1969), Ch. 1, pp. 1−66.
16. F. Bayer-Helms, “Analyse von Linienprofilen. I Grundlagen und Messeinrichtungen,” Z. Angew. Phys. 15, 330–338 (1963).
18. J. M. Llorens, J. Buencuerpo, and P. A. Postigo, “Absorption features of the zero frequency mode in an ultra-thin slab,” Appl. Phys. Lett. 105(23), 231115 (2014). [CrossRef]
19. J. W. Strutt, “Investigations in optics, with special reference to the spectroscope,” Philos. Mag. 8(49), 261–274, 403–411, 477–486 (1879). [CrossRef]
20. I. Juvells, A. Carnicer, J. Ferré-Borrull, E. Martín-Badosa, and M. Montes-Usategui, “Understanding the concept of resolving power in the Fabry-Perot interferometer using a digital simulation,” Eur. J. Phys. 27(5), 1111–1119 (2006). [CrossRef]
21. W. Kahan, “Further remarks on reducing truncation errors,” Appl. Opt. 8(1), 40 (1965).
22. R. Remmert, “Convergent series of meromorphic functions,” in Theory of Complex Functions, 4th ed., S. Axler and F. W. Gehring, eds. (Springer, 1998), Ch. 11, pp. 321−341.