Abstract
When illuminated by a white light source, the discrete resonances of a Fabry-Pérot interferometer (FP) provide a broad bandwidth, comb-like spectrum useful for frequency calibration. We report on the design, construction, and laboratory characterization of two planar, passively stabilized, low finesse (≈40) FPs spanning 380 nm to 930 nm and 780 nm to 1300 nm, with nominal free spectral ranges of 20 GHz and 30 GHz respectively. These instruments are intended to calibrate astronomical spectrographs in radial velocity searches for extrasolar planets. By tracking the frequency drift of three widely-separated resonances in each FP, we measure fractional frequency drift rates as low as 1 × 10−10 day−1. However, we find that the fractional drift rate varies across the three sample wavelengths, such that the drift of two given resonance modes disagrees with the ratio of their mode numbers. We explore possible causes of this behavior, as well as quantify the temperature and optical power sensitivity of the FPs. Our results demonstrate the advancement of Fabry-Pérot interferometers as robust and frequency-stable calibrators for astronomical and other broad bandwidth spectroscopy applications, but also highlight the need for chromatic characterization of these systems.
© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
The Fabry-Pérot interferometer [1] is one of the most basic and important optical instruments, with widespread use in both technical and fundamental applications ranging from telecommunications to gravitational wave detection [2]. Many of its earliest uses were in astronomical spectroscopy; analyzing the interferometer’s fringe pattern while observing the Orion nebula allowed Fabry and Buisson to measure gas compositions and radial velocities [3]. Extending this rich history, the Fabry-Pérot interferometer continues to play a critical role in astronomical spectroscopy today.
In this work we investigate vacuum gap Fabry-Pérot interferometers (here referred to as FPs; etalons; cavities) constructed with plane-parallel mirrors to be used as wavelength references for the calibration of high precision astronomical spectrographs [4–8]. When illuminated by a broadband light source the FP functions as a spectral filter and provides a dense, broad bandwidth, comb-like spectrum of lines spaced in frequency by the cavity free spectral range. For the astronomical spectroscopy we envision, the output of the FP is collected in an optical fiber that feeds the input of a spectrograph. The spectrograph also receives stellar light in a second fiber, and in its focal plane, the FP transmission spectrum then provides a calibration grid against which wavelength shifts in the stellar spectrum are measured.
The FPs studied here are intended to serve as calibrators in exoplanet detection surveys that employ the Doppler radial velocity (RV) technique [9], wherein an unseen orbiting planet induces a periodic Doppler shift of the spectrum of the host star. The most important quality of the FP in this context is its dimensional stability, which is equivalent to the stability of its frequency modes, and is tied to its optical and mechanical properties. For the most demanding RV spectroscopy aimed at detecting a terrestrial mass extrasolar planet (exoplanet) in the habitable zone of the host star, the required RV precision is in the range of $\approx$0.1 m s$^{-1}$ to 0.3 m s$^{-1}$, corresponding to a fractional calibration stability $\lesssim 3\times 10^{-9}$ over timescales of years. In spite of recent progress with even the best cryogenic optical cavities [10,11] this level of stability remains uncertain over months to years, though it is now achievable with laser frequency combs (LFCs) [12–17]. However passive FPs do remain a valuable tool for spectrograph calibration because of their broad spectral coverage ($\geq 1$ octave) and high stability ($\sim 10^{-10}$ fractional) over a period of days, as well as their simplicity and cost relative to a LFC. In this role the passive FP offers a secondary, independent calibration system provided its slow (ideally predictable) drift can be characterized with other frequency or wavelength references on timescales of weeks to months.
For robust use as a wavelength reference, the frequency stability of the FP transmission spectrum must be quantified and its sensitivities to environmental and experimental parameters understood. Here we describe laboratory measurements to characterize the frequency stability and drift of two FP calibrators of similar design and construction (low coefficient of thermal expansion materials, vacuum gap interferometric cavities, plane-parallel mirrors), but spanning different wavelength ranges as determined by the reflectivity response of their mirror coatings. The first is a near-infrared (780 nm to 1300 nm) FP calibrator designed for the Habitable Zone Planet Finder (HPF) [18] spectrograph at the 10 m Hobby-Eberly Telescope at McDonald Observatory. The second is a visible (380 nm to 930 nm) FP calibrator for use with the NEID [19] spectrograph at the 3.5 m WIYN Telescope at Kitt Peak National Observatory. Using an auxiliary LFC and a suite of scanning continuous wave lasers, we precisely and accurately monitor the frequency of multiple FP resonances over periods as long as 1 month. This allows us to characterize the stability of the FP over timescales relevant for astronomical applications. In both FP systems we observe fractional drift rates that are predominantly linear at the level of $\approx 5\times 10^{-10}$ day$^{-1}$.
Our most significant observation is a wavelength-dependent drift of the FP modes that deviates from the prediction that $\Delta \omega /\omega$ be constant. This suggests that knowledge of the frequency and drift rate of a single resonance of a planar FP may not be sufficient to characterize the relative drifts of all other modes, as has been proposed for FPs serving as astronomical calibrators [20,21]. We find that drift rates vary from 50 kHz day$^{-1}$ to 800 kHz day$^{-1}$ ($10^{-10} - 10^{-9}$ fractional; $0.1$ m s$^{-1}$ to 1.0 m s$^{-1}$ RV equivalent) across wavelengths spanning hundreds of nm. We suspect this behavior can be traced to sensitivities of the angular illumination, parallelism and flatness of the FPs; however the exact cause is not yet identified.
In the following we present the instrument designs (Section 3.1) and measurement techniques (Section 3.2) we employ, results of the FP stability measurements (Section 4), and discussion of the design and environmental factors that influence this stability. Our results demonstrate the high stability performance of simple, plane-parallel FP etalons, while emphasizing the importance of additional characterization and mode structure modeling of these instruments to further progress their use as robust spectrograph calibrators.
2. Background
The resonance condition for a FP requires that the round trip optical phase shift be an integer multiple of $2\pi$. This includes both the propagation phase $k2L$ and the reflective mirror phase shift $2\phi _r$ such that
where $m$ is an integer, $L$ the cavity length, $c$ the speed of light, and $k=\omega n/c$ the frequency-dependent propagation constant. We set the index of refraction $n=1$ in vacuum. This relationship implies that the discrete resonant frequencies $\omega _m$ (the cavity longitudinal modes) are where $\phi _r$ is generally frequency dependent.For an ideal, dispersionless FP cavity the frequency response of the resonance modes to changes in the cavity length $\Delta L$ can be illustrated by a simple ‘rubber band’ model in which higher frequency modes undergo larger frequency shifts, but for which the fractional frequency shift $\Delta \omega _m / {\omega _m} = -\Delta L / L$ across the spectrum is constant. The number of modes between two resonance frequencies predicts their relative frequency response to a changing cavity length, analogous to equidistant marks on a rubber band whose separations disperse nonuniformly as the band is stretched. Within this simplified model, the ratio of the frequencies of two modes is independent of $\Delta L$ and given simply by the ratio of their mode numbers $\omega _m/\omega _l=m/l$. Similarly, in response to a cavity length change the fractional shift in frequency $\Delta \omega _m$ is equal to the ratio of mode numbers $\Delta \omega _m / \Delta \omega _l=m/l$.
A more realistic model of the FP cavity must also account for dispersion of the cavity medium and mirror coatings, as well as wavelength-dependent phase shifts arising from imperfections in the surface figure of the mirrors and their parallelism. In addition, cavity length changes occur on short (minute to hour) timescales in response to thermal fluctuations and mechanical stresses, as well as long (day to year) timescales due to a gradual aging of the etalon spacer material and slow thermal and mechanical changes in the environment. In FP cavities built from ULE glass and Zerodur that are intended for laser stabilization, the latter behavior has been seen to result in an exponential change of the cavity resonance frequencies that asymptotically approaches a constant. This ’drift’ rate is typically 1 kHz day$^{-1}$ to 100 kHz day$^{-1}$ over years [22–24] and on timescales of days – weeks can be well approximated as linear. While resonance mode frequency response to various perturbations has been studied previously, to our knowledge a comprehensive model to account for the sum of their behaviors does not yet exist. We expand on previous studies by experimentally identifying, over broad bandwidths and in two independent cavities, the deviation of cavity mode frequencies from the ideal case of $\Delta \omega _m / \Delta \omega _l=m/l$.
3. Experimental setup
3.1 Etalon design
The principle design requirements for etalons as RV spectrograph calibrators are broad bandwidth ($\approx 1$ octave), free spectral range (FSR=$c/(2L)$) of order 10 GHz, moderate finesse ($\sim 10$ to 100), and excellent passive short-term frequency stability. The broad bandwidth is needed to provide calibration across the full spectral range of modern echelle spectrographs, while the FSR is chosen to provide the maximum spectral calibration information. This is typically 3 to 5 times the spectrograph’s resolution, e.g., for a spectrograph with resolution $R = \lambda / \Delta \lambda \approx 50\ 000$ centered at $\lambda = 1 \mu$m the preferred FSR is $\approx 20$ GHz. With the FSR fixed, the etalon’s finesse then determines the linewidth $\Delta \nu$, with a preference for $\Delta \nu < R$ to allow measurement of the spectrograph’s point spread function. This is balanced by the requirement that $\Delta \nu$ be sufficiently large for the etalon to efficiently pass the largest photon flux from a white light illuminator, and together these conditions set an etalon finesse $\mathscr {F} \approx 5 - 100$. This is much lower than typical values in etalons used for atomic spectroscopy and frequency metrology, where $10^5 \lesssim \mathscr {F} \lesssim 5 \times 10^5$ [23,24]. However the desired short-term stability and drift is comparable to the state of the art, representing a relatively unexplored region of design parameter space. While recent technological demonstrations and in situ (at the telescope, ’on-sky’) measurements with low finesse FPs have made evident the improved precision capabilities they can offer for RV spectrographs, there have been few published studies (though see [8,25–27]) quantifying the performance of these FP calibrators in detail and verifying theoretically predicted behavior.
Based on these criteria for an RV spectrograph calibrator, we have designed, built and characterized two FP etalons whose mechanical and optical parameters are given in Table 1. Fig. 1 shows the etalons’ schematic and housing, as well as our characterization approach. Both etalons are fabricated by Light Machinery from ultra-low expansion (ULE) glass with specified zero crossing of the coefficient of thermal expansion near 30 $^{\circ}$C. There is a center bore along the optical axis and a small side vent hole. The etalon mirrors, also made from ULE glass, are wedged and have broad bandwidth anti-reflective (AR) coatings on the outer surfaces to minimize parasitic residual etalons. The mirrors are optically contacted to the spacer with the wedges aligned such that the outer surfaces are also parallel to reduce additional prismatic effects over the large bandwidth of both cavities. The larger dimensions of the visible etalon (relative to the NIR etalon) were chosen to improve alignment tolerances.
The housings for both etalons are similar in construction, and so we discuss the general features, while noting important differences here and in Table 1. The etalons are held in temperature stabilized and evacuated aluminum chambers of dimensions $\approx 40$ cm $\times 40$ cm $\times 20$ cm. Fiber feedthroughs are used to transmit light into the vacuum, and the broad bandwidth light is coupled into and out of the etalons with commercially available off-axis parabolic mirrors (OAPs). The NIR etalon uses 780HP FP/APC single mode fibers, while the visible etalon employs endlessly single mode fibers (ESMFs; NKT Photonics LMA-5) feeding the OAPs.
This inner stage is surrounded by thermal shielding in the evacuated chamber with an ion pump and commercial bench-top pump controller holding vacuum pressure $\leq 2 \times 10^{-7}$ Torr. A commercial bench-top temperature controller is used in conjunction with a thermoelectric heater (TEC) to hold the base of the thermal shielding at a nominal temperature of $32.5\ ^{\circ}\textrm {C}$ for the NIR etalon and $34.3 \ ^{\circ}\textrm {C}$ for the visible etalon. Independent thermistors positioned around the thermal shield are used to verify that the temperature is stable to better than $5 \times 10^{-3}\ ^{\circ}\textrm {C}$. Lab tests indicate the setpoint is within 0.5 $^{\circ}$C of the coefficient of thermal expansion (CTE) zero-crossing for the visible etalon but $\lesssim 5$ $^{\circ}$C from the zero-crossing for the NIR etalon. Fiber-to-fiber throughput at sampled wavelengths of 780 nm, 1064 nm and 1319 nm in the NIR etalon (780 nm, 532 nm and 660 nm in the visible etalon) is $>50\%$. A cutaway mechanical drawing indicative of both instruments is shown in Fig. 1(a).
3.2 Laboratory characterization
To characterize the frequency stability of the FP cavities we monitor the frequencies of multiple cavity resonances relative to a self-referenced optical frequency comb. Full details of this technique are presented in [28]; here we give a brief summary while noting differences from our previous work. A near infrared, self-referenced, octave spanning LFC [29] calibrates the scanning of continuous wave (CW) lasers at 780 nm, 1064 nm and 1319 nm with the NIR etalon (780 nm, 532 nm and 660 nm with the visible etalon, the latter two being the 1064 and 1319 nm laser light frequency doubled in nonlinear crystals) that are transmitted through the FP to simultaneously measure and track the frequencies of the resonance modes in real time. The CW laser power incident on both etalons at all wavelengths is $< 1$ mW. Figure 1(c)–(d) give a schematic overview for both setups. The CW lasers are periodically scanned across their respective FP resonances (of width $\approx 500$ MHz to 900 MHz) by applying a sawtooth frequency scan with 40 s period for the NIR etalon, 100 s for the visible etalon. After transmission the beams are de-multiplexed in fiber and sent to separate detectors where the transmitted power from each laser is photodetected and digitized. In addition, our NIR etalon setup includes a portion of the 780 nm laser light being sent through a vapor cell containing rubidium (Rb), and saturated absorption of the Rb ${\rm D_2}$ lines are observed. These spectroscopic signals provide a secondary frequency reference for the 780 nm light to verify our primary frequency axis calibration technique (see below).
In separate detector channels the CW lasers are heterodyned with the spectrally broadened output of a self-referenced Er:fiber LFC with repetition rate $f_{\textrm {rep}}=250$ MHz. In both setups a portion of the LFC near 1560 nm is frequency doubled in a nonlinear crystal to perform this heterodyne with the 780 nm laser. The LFC is stabilized to an in-house hydrogen maser, which is part of the NIST timescale. This ensures absolute optical uncertainty $<100$ Hz for all times $> 100$ s.
As described in greater detail in [28], our electronic mixing scheme provides a calibration ‘tick’ each time the CW laser moves 62.5 MHz, or at precisely $f_{\textrm {rep}} / 4$. The known frequency spacing of the calibration ticks gives discrete points that map acquisition time to relative frequency. We fit a tenth order polynomial to this data to interpolate the frequency axis between calibration ticks and remove nonlinearities in the laser frequency scan [30].
This calibrates the frequency axis for a single scan of the 780 nm laser to $<1.5$ MHz and the 1064 and 1319 nm lasers to $<100$ kHz in the NIR etalon setup. The larger uncertainty in the 780 nm signal is a consequence of the increased sensitivity of the distributed feedback laser source to noise on the sweep voltage controlling the laser current (contrasted with the 1064 nm and 1319 nm nonplanar ring oscillator lasers swept by thermal tuning). In the visible etalon setup the calibration uncertainty at 780 nm is the same, at 532 nm is $<300$ kHz and at 660 nm is $<1$ MHz. The higher uncertainty in the latter two relative to the NIR setup is due to scanning the lasers over a larger range, doubling the frequency axis and thus doubling its error to calibrate the frequency doubled 1064 nm and 1319 nm light sent into the visible etalon, and lower calibration tick SNR in the visible etalon setup.
Using the respective frequency axis for each wavelength, the FP transmission resonance data are fit with a Lorentzian profile whose centroid is assigned as the peak resonance frequency. Example acquisitions of both the calibration tick signal and FP resonance scan at each of the three wavelengths in the NIR etalon are shown in Fig. 2 (Fig. 3 in the visible etalon). The random (white noise) component of the residuals to the Lorentzian fits suggest an achievable noise-limited centroid measurement precision at or better than $\approx 2$ MHz in a single scan. However Fig. 9 also demonstrates stable structure of unknown origin in the Lorentzian fit residuals, discussed in Appendix A.
Validation of our method for calibrating the frequency axes of etalon resonance scans is achieved by sending a portion of the 780 nm CW light used to track the NIR etalon resonance through a rubidium (Rb) vapor cell to simultaneously track the frequencies of several Doppler-free hyperfine transitions. Rubidium was chosen for its convenience and well known spectroscopy. Fig. 4(a) shows a single trace of the $^{87}\textrm {Rb}\ \textrm {D}_2$ and $^{85}\textrm {Rb}\ \textrm {D}_2$ lines over which the laser is scanned, and Fig. 4(b) shows an example fit of two of those transitions. Figure 4(c) shows the frequency of one of these, the $^{87}\textrm {Rb}\ \textrm {D}_2\ \textrm {F}' = 3,2$ crossover transition, over the full period for which we track the 780 nm etalon resonance. In Fig. 4(d) the Allan deviation of this measurement time series is decreasing out to the longest integration times, in agreement with the expectation that the vapor cell should provide a Rb spectrum stable at the level of 100 kHz. Consistent with the Allan deviation, a linear fit in Fig. 4(c) shows no clear drift of the $^{87}\textrm {Rb}\ \textrm {D}_2\ \textrm {F}' = 3,2$ crossover transition greater than 10 kHz day$^{-1}$. The fluctuations seen in Fig. 4(c) and (d) are consistent with the achieved signal-to-noise ratio and known systematic frequency shifts of the rubidium lines driven primarily by laser power and room temperature variations [31–34]. This secondary frequency reference for the 780 nm light therefore independently confirms the precision of our frequency axis calibration at an order of magnitude below the etalon drift rates we measure in the NIR etalon and at least a factor of 4 below drift rates in the visible etalon (Section 4).
4. Results
While the etalons presented here are, by design, inherently simple opto-mechanical systems, the output spectrum is still sensitive to cavity temperature to some degree. We thus perform experiments under three sets of conditions. In the first, we monitor the NIR etalon resonance frequency at 1064 nm while additionally sending a broadband ($\approx 1$ nm) 976 nm laser into the cavity at optical powers of 2 mW to 20 mW to test resonance frequency response to optical heating $df/dP$ (Section 4.1). In the second experiment, we step the temperature in the NIR etalon vacuum enclosure in $2\ ^{\circ}\textrm {C}$ increments to examine resonance mode thermal response $df/dT$ at 1064 nm and 1319 nm (Section 4.1). In the third experiment, we monitor FP resonance frequencies under low ($<500\ \mu \textrm {W}$) total input power while the cavity is under temperature control to track the drift rate at 780 nm, 1064 nm and 1319 nm in the NIR etalon (780 nm, 532 nm and 660 nm in the visible etalon). This allows us to measure the resonance frequency drift $df/dt$ as the etalon materials age over weeks (Section 4.2). In all experiments we continuously record the etalon transmission peaks and heterodyne frequency marker data series.
4.1 Optical power and temperature sensitivity
In the first experiment we scan the 1064 nm CW across the etalon resonance at $<200\ \mu \textrm {W}$ power into the cavity while additionally sending in 2 mW to 20 mW of broadband 976 nm light spanning $\approx 1$ nm. In each of four successive steps of the broadband source power, we measure the resonance frequency at fixed power to obtain a reference value, step the power, allow the cavity to equilibrate for 12 hr (we empirically find this time ensures the resonance frequency has asymptotically approached a new equilibrium value), and measure the induced shift in the resonance frequency. Fig. 5(a) shows the resulting frequency response. A simple linear fit yields a relation for the frequency shift induced by optical heating that reaches 10 % of the ambient frequency drift presented in Section 4.2 ($df/dt = 260$ kHz day$^{-1}$) at $P_{\textrm {incident}} \approx 4.23$ mW. Input optical powers $<1$ mW thus should not induce heating that affects measurement of the slow background drift. In trials measuring the ambient drift (Section 4.2) we consequently maintain $<500\ \mu \textrm {W}$ total input power so as not to appreciably affect calculated drift rates.
In the second experiment we step the temperature setpoint in the NIR cavity from 32.5 $^{\circ}$C to 38.5 $^{\circ}$C in 2 $^{\circ}$C increments. Figure 5(b) shows the resonance response at 1064 nm and 1319 nm, with trials run for 25 hr (again an empirically determined appropriate time) to allow the resonance frequencies to equilibrate at the new temperature setpoint. These measurements yield simple estimates of the cavity CTE: $2.8 \times 10^{-8}\ ^{\circ}\textrm {C}^{-1}$ and $3.4 \times 10^{-8}\ ^{\circ}\textrm {C}^{-1}$ using the 1064 nm and 1319 nm data respectively. Based on the specifications of the ULE glass used for the spacer, we would expect a lower CTE in this temperature range. However the CTE of the full cavity assembly also depends on the mirror substrates as well as the exact details of the mechanical holding structure. These temperature step trials also indicate we are not at the CTE zero-crossing (where the cavity length response to temperature is minimized) in the NIR etalon, as the frequency response to temperature in Fig. 5(b) is roughly constant over a 6 $^{\circ}$C span. By contrast the same experiment performed with the visible etalon shows negligible effect on the resonance frequencies, indicating this cavity is near the CTE zero-crossing.
To assess resonance mode thermal response to slower and lower amplitude ambient temperature fluctuations, Fig. 5(c) – (f) shows the typical relationship between frequency of the 1064 nm mode and the ambient temperature when the etalon temperature is held constant. This 25 hr segment shows a $\approx 5\ \textrm {mK}\ ^{\circ} \textrm{C}^{-1}$ response of the witness thermistor to lab temperature fluctuations. Yet Fig. 5(f) demonstrates that this weak tracking of the witness temperature with lab temperature causes no discernible systematic effect on the resonance frequency. This is consistent with what one would expect from the theoretical CTE, where a 5 mK change would correspond to a $\approx 50$ kHz resonance frequency shift. We thus do not find thermal sensitivity to variations in lab temperature to substantially affect our measured drift rates.
4.2 Frequency drift
The frequency (or wavelength) stability of the FP etalon modes over time is the most important characteristic that our measurements allow us to address. During these measurements the temperature of the etalon thermal shield is being actively controlled at its fixed point, but otherwise the outer vacuum housing experiences the ambient laboratory environment, with the lab temperature fluctuating over a maximum range of $\approx 1$ $^{\circ}$C. Using the techniques described in Section 3.2 we measure the drift of the NIR resonance frequencies at 780 nm, 1064 nm and 1319 nm over periods $>1$ week. Results are shown in Fig. 6, with a linear fit to the measured mode frequencies over time yielding observed drift rates of $df/dt =$783(7) kHz day$^{-1}$, 260(1) kHz day$^{-1}$ and 171(1) kHz day$^{-1}$ (see Table 2 for a summary of these and the following quantities). These correspond to a fractional frequency stability of the respective etalon modes of $2 \times 10^{-9}$, $9 \times 10^{-10}$ and $8 \times 10^{-10}$ day$^{-1}$. A simple model would predict that the ratio of the drift rates should scale as the ratio of their mode numbers ($df_{780} / df_{1064} = 1.37$, $df_{1064} / df_{1319} = 1.24$ and $df_{780} / df_{1319} = 1.70$; Section 2). However the observed ratios differ markedly from this at 3.02(0.03), 1.52(0.01) and 4.57(0.05) respectively. This discrepancy cannot be explained by uncertainties in our measurement approach, and possible physical mechanisms are discussed below. After removing a linear fit to the drift we calculate the Allan deviation, a measure of the frequency stability as a function of averaging time. As shown in Fig. 6, this indicates that under the assumption of linear drift we are able to determine the line center at all sample resonances to better than 300 kHz for averaging periods of 1 hr up to the maximum averaging time of $\approx 100$ hr. Importantly the different scatter in the resonance measurements and the computed Allan deviations are indicative of the measurement SNR rather than an intrinsic property of a specific resonance.
Analogously Fig. 7 shows the ambient drift in the visible etalon at 780 nm, 532 nm and 660 nm. Here the respective observed drift rates, $df/dt =$36(1) kHz day$^{-1}$, 129(2) kHz day$^{-1}$ and 48(2) kHz day$^{-1}$, correspond to a frequency stability of the respective etalon modes of $9 \times 10^{-11}$ day$^{-1}$, $2 \times 10^{-10}$ day$^{-1}$ and $1 \times 10^{-10}$ day$^{-1}$ (see Table 2 for a summary). Again a simple model of the FP would predict drift rate ratios of $df_{532} / df_{780} = 1.47$, $df_{532} / df_{660} = 1.24$ and $df_{660} / df_{780} = 1.18$, but instead we observe 3.60(0.11), 2.66(0.12) and 1.36(0.07) respectively. The larger Allan deviations at a given averaging time relative to those in the NIR etalon data are due to lower SNR in the visible etalon measurement setup. In particular the measurements of the 660 nm resonance data show increased fluctuations when compared to the measurements at other wavelengths, likely due to increased sensitivity to laser amplitude noise and laboratory temperature perturbations.
The ability to accurately predict the drift rate of the FP cavity resonances is critical for their use as standalone frequency calibrators. Fig. 8(a)–(b) show how the measured drift rates in both etalons evolve as we fit a straight line for the linear drift to successively larger portions of the cumulative datastream. As would be expected, the ability to predict the drift rates improves with measurement time. And while the data for the NIR cavity is limited in duration, the visible etalon data shows evidence of a nonlinear change in the drift rate, at least over the measurement time of $\approx 1$ month. We note that measurements with the visible etalon were begun within a few weeks of the mirrors being optically contacted. This slowing of the drift rate could be indicative of the cavity "relaxing" in its mechanical structure and/or a slow change in the cavity length itself. Particularly encouraging for the visible cavity is that the drift rates at all three colors are evolving to values $\lesssim 100$ kHz day$^{-1}$, equivalent to fractional drift rates $\lesssim 2\times 10^{-10}$ day$^{-1}$. For astronomical RV calibrations this corresponds to 5 cm s$^{-1}$ day$^{-1}$.
5. Discussion
The experiments and data reported above demonstrate the excellent frequency stability that can be obtained from broad bandwidth FP etalons designed for astronomical spectrograph calibration. Slow, predictable and largely linear drifts should enable such FP etalons to serve as efficient standalone calibrators on timescales of days – weeks, and cross calibration with other sources (lamps, gas cells or frequency combs) should provide enhanced long-term stability. However the most striking observation of our work is the discrepancy between the theoretical and observed drift ratios, i.e., an observed chromatic variation in the linear drift that does not fit the simple model of a FP etalon. Our measurements indicate a wavelength dependence of the FP modes beyond that predicted by the ratio of mode numbers, which cannot be accounted for by the uncertainties in the present laboratory measurements. The significance of our current result is that knowledge (and control) of the frequency of a single etalon mode does not reliably predict the frequencies of other modes, especially over broad bandwidths. This observed behavior could arise from one or a combination of multiple aspects of the etalon’s construction and illumination. While we do not have conclusive results on the exact mechanism, our experiments do allow us to suggest some possibilities.
One of the most striking properties of the plane-parallel FP cavity is its sensitivity to input alignment – that is, the degree of parallelism between the input beam and the cavity’s optical axis. In experiments conducted before the visible etalon was placed in vacuum, the CW lasers were coupled into a single mode fiber, the output of which was collimated and aligned to the cavity by maximizing the power reflected back into the same optical fiber. We then scanned the incidence angle onto the FP in one dimension by adjusting a kinematic mirror mount before the etalon and sampled the resulting resonance frequency shifts at our three colors. In this experiment we measured $\sim 1$ MHz $\mu$rad$^{-1}$ sensitivity of the mode frequencies at 780 nm, 660 nm and 532 nm to input alignment. This suggests that wavelength-dependent alignment errors $\leq 1\ \mu$rad varying differently in time could result in frequency shifts with magnitude comparable to that observed in our ambient drift measurements. For example, although care was taken in collimation of the broad bandwidth light incident on the cavity, small differences in angular divergence could lead to different colors sampling different spatial regions of the cavity. This could then couple to spatial variations in the surface figure of the mirrors or to a relaxing of the cavity that is not radially or azimuthally symmetric relative to the cavity axis. Such a hypothesis is difficult to model or verify experimentally, but it would be informative if revealed in other FP cavities of similar construction that are in use or intended for similar application. We note that in a previous experiment [28] we did recover the expected dispersion stability $\Delta \omega _{1064} / \Delta \omega _{1319}$ at the $10^{-3}$ level in a fiber Fabry-Pérot using the measurement technique and data analysis pipeline we employ here. Though note for the fiber Fabry-Pérot, the fiber is single mode at all wavelengths, thus ensuring ideal spatial overlap of all wavelengths.
Another possible cause of the observed anomalous wavelength dependence of the etalon drift is the frequency-dependent phase upon reflection. For example if the reflective phase shift is strongly varying near one or more of the three test wavelengths we measure, then as the cavity length changes by $\Delta L$ the different resonant frequencies would be given by
6. Conclusions
We have described the construction and characterization of two evacuated, broadband, planar Fabry-Pérot etalons that are to be used as spectral calibrators for precision astronomical RV spectroscopy. Using an auxiliary laser frequency comb, we track three resonances in each etalon over a time frame of 1 week to 1 month. We observe that the drift rates of the cavity resonances are slow and essentially linear over this same time frame, with magnitudes $< 2\times 10^{-9}$ day$^{-1}$. However we also find that the ratios of drift rates among the monitored resonances differ significantly from the theoretical expectation. The cause is likely a combination of factors including alignment sensitivity and cavity imperfections over extremely broad bandwidths, as well as unknown factors related to the details of the dispersion and stability of the cavity mirror coatings. An important outcome of our work is that tracking one resonance mode of such a planar etalon is likely insufficient to provide full information about the frequency drift of other modes via a simple ratio of the frequencies. This may be of importance to the range of users in the astronomical community who have adopted etalon-based calibrators.
Appendix A. Etalon resonance residual structure
Residuals to a Lorentzian fit of the etalon resonances in both the NIR and visible cavities show stable structure of unknown origin. In Fig. 9(a) – (c), 2D histograms show this systematic structure at the $<1\%$ level over 12 hr ($\approx 500$ stacked 40 s resonance scans) in the NIR etalon resonance fit residuals at 1064 nm, 1319 nm and 780 nm respectively. We verified that the structure is not due to approximating the resonance Airy profile by a Lorentzian; for the data in Fig. 9, the difference of an Airy fit and our Lorentzian fit at each of the 3 test wavelengths does not exceed $10^{-5}$ fractional of the resonance profile’s amplitude (i.e., is a factor of 1000 below the systematic residual structure). The stable residual structure is instead likely a consequence of imperfections that distort the transmission profile, such as mirror surface defects that would cause a position-dependent mirror separation [36–38] or slight mirror non-parallelism that would produce an asymmetric line profile [36,39]. It is unclear whether this structure is a consequence of the same mechanism(s) causing the anomalous wavelength dependence of the etalon drift. Figure 9(d) shows the similarity of the residual structure to the third derivative of the Lorentzian fit to a single resonance scan at 1064 nm. Parasitic etalons are also seen distorting the stable residual structure in Fig. 9(e) – (f). In many cases these could be removed by adjusting the upstream optics. For those that could not be removed, we test their efficacy to bias our resonance fits and explain the anomalous wavelength dependence of the etalon drift. For the data in Fig. 9(a) for example, we take the residual amplitude with the largest count at each bin along the frequency axis and subtract this from the raw resonance profile in each of the $\approx 500$ scans used to produce the plot. For each scan we then fit a Lorentzian to this data from which the stable residual structure is subtracted and compare the fitted linecenter with that obtained by fitting to the raw data. The difference between the fitted linecenters is $<1$ kHz, demonstrating that this stable residual structure does not systematically affect our measured drift rates.
Funding
National Aeronautics and Space Administration (1547612); National Science Foundation (AST-1006676, AST-1126413, AST-1310875, AST-1310885, AST-1517592).
Acknowledgments
The mention of trade names and specific products is for scientific information only and does not imply an endorsement of such products by NIST. This work was performed by SPH [in part] under contract with the Jet Propulsion Laboratory (JPL) funded by NASA through the Sagan Fellowship Program executed by the NASA Exoplanet Science Institute. We thank Rich Fox and Daniel Nicolodi for their thoughtful comments on the manuscript. JJ thanks A. Bagga for his discussions on the work.
Disclosures
MG, MN: Stable Laser Systems (E).
References
1. C. Fabry and A. Perot, “On a New Form of Interferometer,” Astrophys. J. 13, 265 (1901). [CrossRef]
2. J. M. Vaughan, The Fabry-Perot interferometer. History, theory, practice and applications (1989).
3. H. Buisson, C. Fabry, and H. Bourget, “An application of interference to the the study of the Orion nebula,” Astrophys. J. 40, 241–258 (1914). [CrossRef]
4. F. Wildi, B. Chazelas, and F. Pepe, “A passive cost-effective solution for the high accuracy wavelength calibration of radial velocity spectrographs,” in Ground-based and Airborne Instrumentation for Astronomy IV, vol. 8446 of Proc. Soc. Photo-Opt. Instrum. Eng. (2012), p. 84468E.
5. T. M. McCracken, C. A. Jurgenson, D. A. Fischer, R. A. Stoll, A. E. Szymkowiak, J. Bradford, and W. Rutter, “Single-lock: a stable Fabry-Perot based wavelength calibrator,” in Ground-based and Airborne Instrumentation for Astronomy V, vol. 9147 of Proc. Soc. Photo-Opt. Instrum. Eng. (2014), p. 91473L.
6. S. Halverson, S. Mahadevan, L. Ramsey, F. Hearty, J. Wilson, J. Holtzman, S. Redman, G. Nave, D. Nidever, M. Nelson, N. Venditti, D. Bizyaev, and S. Fleming, “Development of fiber fabry-perot interferometers as stable near-infrared calibration sources for high resolution spectrographs,” Publ. Astron. Soc. Pac. 126(939), 445–458 (2014). [CrossRef]
7. A. Reiners, R. K. Banyal, and R. G. Ulbrich, “A laser-lock concept to reach cm-precision in doppler experiments with fabry-pérot wavelength calibrators,” Astron. Astrophys. 569, A77 (2014). [CrossRef]
8. F. F. Bauer, M. Zechmeister, and A. Reiners, “Calibrating echelle spectrographs with Fabry-Pérot etalons,” Astron. Astrophys. 581, A117 (2015). [CrossRef]
9. C. Lovis and D. Fischer, Radial Velocity Techniques for Exoplanets (2010), pp. 27–53.
10. C. Hagemann, C. Grebing, C. Lisdat, S. Falke, T. Legero, U. Sterr, F. Riehle, M. J. Martin, and J. Ye, “Ultrastable laser with average fractional frequency drift rate below 5 × 10−19/s,” Opt. Lett. 39(17), 5102 (2014). [CrossRef]
11. J. M. Robinson, E. Oelker, W. R. Milner, W. Zhang, T. Legero, D. G. Matei, F. Riehle, U. Sterr, and J. Ye, “Crystalline optical cavity at 4 k with thermal-noise-limited instability and ultralow drift,” Optica 6(2), 240–243 (2019). [CrossRef]
12. M. T. Murphy, T. Udem, R. Holzwarth, A. Sizmann, L. Pasquini, C. Araujo-Hauck, H. Dekker, S. D’Odorico, M. Fischer, T. W. Hänsch, and A. Manescau, “High-precision wavelength calibration of astronomical spectrographs with laser frequency combs,” Mon. Not. R. Astron. Soc. 380(2), 839–847 (2007). [CrossRef]
13. D. A. Braje, M. S. Kirchner, S. Osterman, T. Fortier, and S. A. Diddams, “Astronomical spectrograph calibration with broad-spectrum frequency combs,” Eur. Phys. J. D 48(1), 57–66 (2008). [CrossRef]
14. C.-H. Li, A. J. Benedick, P. Fendel, A. G. Glenday, F. X. Kärtner, D. F. Phillips, D. Sasselov, A. Szentgyorgyi, and R. L. Walsworth, “A laser frequency comb that enables radial velocity measurements with a precision of 1cms−1,” Nature 452(7187), 610–612 (2008). [CrossRef]
15. T. Steinmetz, T. Wilken, C. Araujo-Hauck, R. Holzwarth, T. W. Hänsch, L. Pasquini, A. Manescau, S. D’Odorico, M. T. Murphy, T. Kentischer, W. Schmidt, and T. Udem, “Laser Frequency Combs for Astronomical Observations,” Science 321(5894), 1335–1337 (2008). [CrossRef]
16. G. G. Ycas, F. Quinlan, S. A. Diddams, S. Osterman, S. Mahadevan, S. Redman, R. Terrien, L. Ramsey, C. F. Bender, B. Botzer, and S. Sigurdsson, “Demonstration of on-sky calibration of astronomical spectra using a 25 GHz near-IR laser frequency comb,” Opt. Express 20(6), 6631 (2012). [CrossRef]
17. T. Wilken, G. L. Curto, R. A. Probst, T. Steinmetz, A. Manescau, L. Pasquini, J. I. González Hernández, R. Rebolo, T. W. Hänsch, T. Udem, and R. Holzwarth, “A spectrograph for exoplanet observations calibrated at the centimetre-per-second level,” Nature 485(7400), 611–614 (2012). [CrossRef]
18. S. Halverson, S. Mahadevan, L. Ramsey, R. Terrien, A. Roy, C. Schwab, C. Bender, F. Hearty, E. Levi, S. Osterman, G. Ycas, and S. Diddams, “The habitable-zone planet finder calibration system,” in Ground-based and Airborne Instrumentation for Astronomy V, vol. 9147 of Proc. Soc. Photo-Opt. Instrum. Eng. (2014), p. 91477Z.
19. C. Schwab, A. Rakich, Q. Gong, S. Mahadevan, S. P. Halverson, A. Roy, R. C. Terrien, P. M. Robertson, F. R. Hearty, E. I. Levi, A. J. Monson, J. T. Wright, M. W. McElwain, C. F. Bender, C. H. Blake, J. Stürmer, Y. V. Gurevich, A. Chakraborty, and L. W. Ramsey, “Design of NEID, an extreme precision Doppler spectrograph for WIYN,” in Ground-based and Airborne Instrumentation for Astronomy VI, vol. 9908 of Proc. Soc. Photo-Opt. Instrum. Eng. (2016), p. 99087H.
20. C. Schwab, J. Stürmer, Y. V. Gurevich, T. Führer, S. K. Lamoreaux, T. Walther, and A. Quirrenbach, “Stabilizing a fabry-perot etalon peak to 3 cm s−1 for spectrograph calibration,” Publ. Astron. Soc. Pac. 127(955), 880–889 (2015). [CrossRef]
21. J. Stürmer, A. Seifahrt, C. Schwab, and J. L. Bean, “Rubidium-traced white-light etalon calibrator for radial velocity measurements at the cm s−1 level,” J. Astron. Telesc. Instrum. Syst. 3(2), 025003 (2017). [CrossRef]
22. M. Zhu and J. L. Hall, “Short and long term stability of optical oscillators,” in Proceedings of the 1992 IEEE Frequency Control Symposium, (1992), pp. 44–55.
23. C. Salomon, D. Hils, and J. L. Hall, “Laser stabilization at the millihertz level,” J. Opt. Soc. Am. B 5(8), 1576–1587 (1988). [CrossRef]
24. D. G. Matei, T. Legero, S. Häfner, C. Grebing, R. Weyrich, W. Zhang, L. Sonderhouse, J. M. Robinson, J. Ye, F. Riehle, and U. Sterr, “1.5 µ m Lasers with Sub-10 mHz Linewidth,” Phys. Rev. Lett. 118(26), 263202 (2017). [CrossRef]
25. F. Cersullo, F. Wildi, B. Chazelas, and F. Pepe, “A new infrared Fabry-Pérot-based radial-velocity-reference module for the SPIRou radial-velocity spectrograph,” Astron. Astrophys. 601, A102 (2017). [CrossRef]
26. S. Schäfer, E. W. Guenther, A. Reiners, J. Winkler, M. Pluto, and J. Schiller, “Two Fabry-Pérots and two calibration units for CARMENES,” in Ground-based and Airborne Instrumentation for Astronomy VII, vol. 10702 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series (2018), p. 1070276.
27. F. Cersullo, A. Coffinet, B. Chazelas, C. Lovis, and F. Pepe, “New wavelength calibration for echelle spectrographs using Fabry-Pérot etalons,” Astron. Astrophys. 624, A122 (2019). [CrossRef]
28. J. Jennings, S. Halverson, R. Terrien, S. Mahadevan, G. Ycas, and S. A. Diddams, “Frequency stability characterization of a broadband fiber Fabry-Pérot interferometer,” Opt. Express 25(14), 15599 (2017). [CrossRef]
29. G. Ycas, S. Osterman, and S. A. Diddams, “Generation of a 660-2100 nm laser frequency comb based on an erbium fiber laser,” Opt. Lett. 37(12), 2199 (2012). [CrossRef]
30. One may expect the nonlinearity in the laser scans to induce a signal noise (‘ringing’) in the calibration tick signal as the heterodyne beat between the CW laser and comb passes through the RF bandpass windows used to generate the ticks. We find no such discernible ringing and thus are confident that our Gaussian fits to individual ticks are not systematically biasing the separation between ticks over the nonlinear scan, i.e., that we are justified in assuming the separation between the fitted Gaussians to neighboring ticks is 62.5 MHz.
31. H. Furuta and M. Ohtsu, “Evaluations of frequency shift and stability in rubidium vapor stabilized semiconductor lasers,” Appl. Opt. 28(17), 3737–3743 (1989). [CrossRef]
32. C. Affolderbach and G. Mileti, “Tuneable, stabilised diode lasers for compact atomic frequency standards and precision wavelength references,” Opt. Lasers Eng. 43(3-5), 291–302 (2005). [CrossRef]
33. W. Loh, M. T. Hummon, H. F. Leopardi, T. M. Fortier, F. Quinlan, J. Kitching, S. B. Papp, and S. A. Diddams, “Microresonator brillouin laser stabilization using a microfabricated rubidium cell,” Opt. Express 24(13), 14513–14524 (2016). [CrossRef]
34. T. P. Dinneen, C. D. Wallace, and P. L. Gould, “Narrow linewidth, highly stable, tunable diode laser system,” Opt. Commun. 92(4-6), 277–282 (1992). [CrossRef]
35. F. Riehle, “Use of optical frequency standards for measurements of dimensional stability,” Meas. Sci. Technol. 9(7), 1042–1048 (1998). [CrossRef]
36. R. M. Hill, “Some Fringe-broadening Defects in a Fabry-Perot Étalon,” Opt. Acta 10(2), 141–152 (1963). [CrossRef]
37. P. A. Wilksch, “Instrument function of the Fabry-Perot spectrometer,” Appl. Opt. 24(10), 1502–1511 (1985). [CrossRef]
38. K. P. Reardon and F. Cavallini, “Characterization of Fabry-Perot interferometers and multi-etalon transmission profiles. The IBIS instrumental profile,” Astron. Astrophys. 481(3), 897–912 (2008). [CrossRef]
39. D. P. Mahapatra and S. K. Mattoo, “Exact evaluation of the transmitted amplitude for a Fabry-Perot interferometer with surface defects,” Appl. Opt. 25(10), 1646–1649 (1986). [CrossRef]