## Abstract

Improved wavelength calibrators for high-resolution astrophysical spectrographs will be essential for precision radial velocity (RV) detection of Earth-like exoplanets and direct observation of cosmological deceleration. The astro-comb is a combination of an octave-spanning femtosecond laser frequency comb and a Fabry-Pérot cavity used to achieve calibrator line spacings that can be resolved by an astrophysical spectrograph. Systematic spectral shifts associated with the cavity can be 0.1-1 MHz, corresponding to RV errors of 10-100 cm/s, due to the dispersive properties of the cavity mirrors over broad spectral widths. Although these systematic shifts are very stable, their correction is crucial to high accuracy astrophysical spectroscopy. Here, we demonstrate an *in-situ* technique to determine the systematic shifts of astro-comb lines due to finite Fabry-Pérot cavity dispersion. The technique is practical for implementation at a telescope-based spectrograph to enable wavelength calibration accuracy better than 10 cm/s.

©2010 Optical Society of America

## 1. Introduction

High precision wavelength calibrators for astrophysical spectrographs will be key components of new precision radial velocity (RV) observations, including the search for Earth-like extra-solar planets (exoplanets) [1] and direct observation of cosmological acceleration [2, 3]. Recent work has demonstrated the potential of octave-spanning femtosecond laser frequency combs [4] (“astro-combs”) to serve as wavelength calibrators for astrophysical spectrographs providing RV sensitivity down to 1 cm/s [5–11]. Exoplanet searches place stringent demands upon such calibrators. For example, the RV amplitude of the reflex motion of a solar-mass star induced by an Earth-mass planet in an orbit within the habitable zone is about 10 cm/s. The current state of the art astrophysical spectrograph, HARPS, has demonstrated stellar RV sensitivity ≈ 60 cm/s [1], largely limited by its thorium argon lamp calibrator [12, 13]. These calibrators are limited by their unevenness in line spacing and intensity as well as the slow variation of their line wavelengths with time. An astro-comb provides emission lines with uniform intensity and controllable spacing, which can be referenced to atomic frequency standards and the global positioning system (GPS), yielding excellent long-term stability and reproducibility.

To date, astro-combs consist of a combination of an octave-spanning femtosecond laser frequency comb (source comb) and a Fabry-Pérot cavity (FPC), see Fig. 1. Spectral lines generated by the source comb are spaced by the pulse repetition rate (*f _{r}*), currently ≈ 1 GHz, which results in a line spacing too dense to be resolved by broadband astrophysical spectrographs [10]. The FPC serves as a mode filter with a free spectral range (FSR) set to an integer multiple of the repetition rate, FSR=

*Mf*, with

_{r}*M*≈ 10−100, depending on the spectrograph resolution. Ideally, the FPC passes only every

*M*source comb spectral line, providing thousands of calibration lines well matched to a practical spectrograph’s resolution, with fractional frequency uncertainty limited only by the stability of the RF reference used to stabilize the source comb and FPC. This frequency uncertainty can be < 10

^{th}^{−11}using commonly available atomic clock technology, which corresponds to ~ 3 kHz uncertainty in the optical frequency or 0.3 cm/s in RV precision. However, because the spectrograph fails to resolve neighboring source comb lines, finite suppression of these neighboring lines by the FPC affects the lineshape and potentially the centroid of measured astro-comb lines. For example, in the results presented here, source comb modes neighboring the astro-comb line, with intensities after passing through the FPC that differ by 0.1% of the main astro-comb peak, shift the measured line centroid by 1 MHz, which corresponds to an RV systematic error of 1 m/s. In practice, such systematic RV shifts are inevitable over spectral bandwidths of 1000 Å due to the dispersive properties of the mirrors of the FPC. Although these systematic shifts can be very stable over timescales of years, the correction of such shifts is crucial to high accuracy astrophysical spectroscopy.

In this paper, we demonstrate an *in-situ* technique to determine the systematic shifts of astro-comb lines due to FPC dispersion, which can be applied at a telescope-based spectrograph to enable wavelength calibration accuracy better than 10 cm/s. By measuring the intensity of astro-comb lines as the FPC length is adjusted, we determine (i) the offset of each FPC resonance from the nearest source comb line; (ii) FPC finesse as a function of wavelength; and (iii) the intensity of the astro-comb lines and their neighboring (normally suppressed) source comb lines. These parameters can be determined quickly and reliably over the full 1000 Å bandwidth of the astro-comb with only ≈ 50 measurements at slightly different FPC lengths, and can be performed quickly (< 1 hour) and reliably. The measurement has also been performed with a lower resolution commercial optical spectrum analyzer with consistent results. The astro-comb line characterization technique presented here builds on past work in which femtosecond lasers coupled to swept cavities were used to study both the medium in which the cavity was immersed and the cavity mirrors [14, 15].

## 2. Model

Imperfect suppression of unwanted source comb lines, e.g., due to finite FPC finesse, affects the astro-comb lineshape observed on a spectrograph. The lineshape can be modeled with knowledge of the FPC properties including mirror reflectivity and round trip phase delay. The intensity of a source comb line after the FPC is

where *I _{m}* is the intensity of the source comb line of optical frequency

*f*;

_{m}*T*is the resonant transmission of the FPC at optical frequency

_{m}*f*;

_{m}*F*is the finesse of the FPC near frequency

_{m}*f*; and

_{m}*ϕ*is the round trip phase delay. The phase delay may be expressed as [17]

_{m}where *L* is the length of the cavity, *c* is the speed of light in vacuum, *n _{m}* is the refractive index of the medium inside the cavity (air, vacuum, etc.) at optical frequency

*f*, and

_{m}*τ*(

*f*) is the frequency-dependent group delay of the mirrors. The first term in parentheses in Eq. (2) is the distance between the mirrors expressed in wavelengths; while the second term, the integral of

*τ*(

*f*), is the phase delay of the mirrors and represents the frequency-dependent penetration distance of light into the mirror. Maximum transmission occurs when a source comb line is resonant with the FPC, or equivalently

*ϕ*= 2

*πq*, with

*q*an integer [see Fig. 2(a)]. The FSR, Δ, is the frequency difference between two consecutive FPC resonant frequencies. Assuming that

*τ*varies slowly with frequency, the FSR can be approximated by

For the astro-comb shown in Fig. 1, which we deployed as a wavelength calibrator for the Tillinghast Reflector Echelle astrophysical Spectrograph (TRES) [16], typical operational parameters are Δ* _{m}* ≈ 31 GHz and

*F*≈ 180, with source comb lines nearest to the astro-comb peak suppressed by ≈ 22 dB [see Fig. 2(b) and Section 3]. This imperfect suppression leads to systematic inaccuracies in astro-comb line centers at the 50 cm/s level, as observed on the TRES spectrograph across a 1000 Å bandwidth. The effects of these systematic shifts can be characterized by determining the FPC finesse (

_{m}*F*) and the frequency difference between astro-comb lines and the FPC resonance over the full spectral width. As described below, compensation can then be applied to determine wavelength solutions (conversions from spectrograph pixel number to wavelength) with accuracy at least an order of magnitude below the 50 cm/s level of these systematic shifts.

_{m}We determine *F _{m}* and the frequency offset of the FPC resonance from the astro-comb lines by adjusting the length of the FPC and measuring the intensities of astro-comb lines over the full spectral bandwidth. When we change the frequency of one FPC resonance from

*f*to

_{d}*f*+

_{d}*δf*the cavity length changes by

where Δ* _{d}* and

*τ*are the FSR of the FPC and the phase delay of the mirrors, respectively, at the frequency of the controlled resonance. The change of the cavity length leads to a change of the round trip phase delay of an astro-comb line

_{d}*f*by

_{m}As a result, the intensity of the astro-comb line varies with the change of the frequency of one FPC mode according to Eq. (1), from which we can derive the finesse and the phase errors at all astro-comb line frequencies.

## 3. Experiment

The astro-comb employed in our experimental demonstrations is shown schematically in Fig. 1. The octave-spanning source comb spectrum (Fig. 3) is generated by a mode-locked titaniumsapphire femtosecond laser (Octavius, Menlosystems, Inc.) with repetition rate *f _{r}* ≈ 1 GHz. The absolute frequencies of the comb lines can be expressed as

*f*=

_{m}*f*

_{0}+

*m f*, where

_{r}*m*is an integer and

*f*

_{0}is the carrier-envelope offset frequency. A PIN diode detects

*f*, which is then stabilized by adjusting the laser cavity length. The

_{r}*f*-2

*f*self-referencing method is used to produce a signal at

*f*

_{0}on an avalanche photodiode: comb lines around 11400 Å are frequency doubled in a 1 mm thick LBO (Lithium Triborate) crystal and beat with comb lines around 5700 Å.

*f*

_{0}is then stabilized by intensity modulation of the 7.6W, 5320 Å pump laser. Both

*f*and

_{r}*f*

_{0}are referenced to low-noise radio-frequency synthesizers, which are stabilized to a commercial rubidium frequency reference. For our source comb, typical values of these frequencies are:

*f*= 1.042048458 GHz and

_{r}*f*

_{0}= 0.09 GHz, and the resulting linewidth of individual source-comb spectral lines is < 1 MHz.

The source comb beam passes through the FPC, filtering out unwanted comb lines and increasing the transmitted line spacing. Two flat mirrors with ≈ 98.5% reflectivity and minimal group delay dispersion (<10 fs^{2} from 7500 Å to 9000 Å) comprise the FPC. Due to the low dispersion mirrors, the FSR of the FPC is almost constant within the 1000 Å bandwidth: typically one in thirty comb lines is resonant with the FPC [Fig. 2(a)]. The measured finesse of the FPC is ≈ 180, which is consistent with the theoretical limit estimated from the mirror reflectivity.

We stabilize the FPC by locking one transmission resonance to an injected diode laser. The diode laser is modulated at 30 MHz with an electro-optic modulator (EOM) and injected into the FPC with a polarization orthogonal to the comb light. After the FPC, the diode beam is separated from the comb, demodulated, and used to derive a feedback signal that stabilizes the FPC to the diode laser. The diode laser itself is offset-locked to one of the source comb lines at a wavelength ≈ 7950 Å. The frequency of the offset lock [see Fig. 2(a)] is adjusted to optimize the astro-comb bandwidth, and compensates for FPC dispersion between the diode laser wavelength and the central astro-comb wavelength due to mirror properties and air in the cavity as discussed in Secs. 2 and 4. The diode laser wavelength is close to the Rb D1 transition. Absorption spectroscopy of the diode laser light using a thermal rubidium vapor cell thus identifies the diode frequency to within one source comb line. The offset lock then determines the absolute frequency of the diode laser with the accuracy of the underlying atomic frequency reference.

In its current incarnation as a wavelength calibrator for an astrophysical spectrograph, the astro-comb spectrum is analyzed by the Tillinghast Reflector Echelle Spectrograph (TRES) [16], a fiber-fed multi-order echelle spectrograph at the 1.5 m telescope of the Whipple Observatory at Mt. Hopkins in Arizona. Spectral dispersion is provided by an echelle grating with cross-dispersion from a prism operated in double-pass mode. TRES covers a spectral bandwidth from 3700 Å to 9100 Å with a resolving power of *R* ≡ *λ*/*δλ* ≈ 50000 [16] corresponding to a resolution of ≈ 7 GHz at 8000 Å. Astro-comb light is injected into an integrating sphere and then sent to TRES on a 100-*µ*m multimode fiber. The integrating sphere reduces lineshape fluctuations due to input-dependent illumination of the spectrograph optics. Each astro-comb spectral measurement by TRES is typically integrated for 60 s and recorded on a cooled (100 K) two dimensional CCD [see Fig. 2(b)]. Flat-field correction is applied to the astro-comb spectrum in order to remove artifacts caused by variations in the pixel-to-pixel sensitivity of the detector and by distortions in the optical path. The implementation of both the integrating sphere and the flat-field correction are essential for high accuracy astrophysical spectroscopy but is not crucial to the work described in this paper. The separation of the FPC mirrors is set such that astro-comb lines are separated by ≈ 4 resolution elements of the spectrograph. In the measurement of the intensity of each astro-comb lines, typically the counts from 18 CCD pixels around the central pixel are binned from the one dimensional extracted spectrograph spectrum to obtain each measurement points in Fig. 4. The 18 pixel integrations capture all the intensity from each astro-comb line (spectrograph resolution ≈ 6 pixels) while avoiding cross-talk between neighboring lines (line separation ≈ 24 pixels). We have also varied the binning and found no appreciable change in the result. The length of the FPC was swept (by changing the diode frequency) back and forth through the FPC resonance: we took steps of 80 MHz in one direction and then steps in the reverse direction at intermediate positions. Different step sizes or schemes have led to consistent results.

We also characterized the astro-comb spectrum with a commercial optical spectrum analyzer (OSA, ANDO AQ6315), which has much lower resolution (≈ 300 GHz) when operated in broadband (1000 Å) mode. There are then ≈ 10 astro-comb lines in one resolution element. Extraction of astro-comb intensities from the OSA is significantly simpler than the two dimensional TRES spectrograph, and the OSA, therefore, provides a useful cross-check. We find that all parameters measured with the OSA are consistent with those measured with TRES. Additionally, the calibration procedure has been performed more than 10 times with both the TRES spectrograph and the OSA over 10 days. All measurements are found to be consistent over this time period.

## 4. Results and Analysis

We find good agreement between TRES measurements of the astro-comb spectrum and the model presented above. For example, Figure 4 shows the measured variation of the peak intensity of one astro-comb line vs diode laser frequency (and thus cavity length), as well as a fit of Eq. (1) to the data. From the fit, we derive the astro-comb resonant intensity *T _{m}I_{m}*; the FPC finesse

*F*; and the phase deviation

_{m}*δϕ*. The phase deviation is the offset of the round trip phase from an integer multiple of 2

_{m}*π*at an astro-comb line, and is determined from the detuning of the astro-comb peak at frequency

*f*from center of the FPC resonance at nominal diode laser setting

_{m}*f*

_{FPC}. The phase deviation is then given by

where Δ* _{m}* is the free spectral range near frequency

*f*. The finesse,

_{m}*F*, of the FPC can be derived from the fit as

where Δ*ϕ*
_{FWHM} is the full width at half maximum in phase of an FPC resonance. The uncertainties in measured FPC resonance centers and widths are approximately 1 MHz, corresponding to phase uncertainties of 0.2 mrad. The variation of these parameters as a function of wavelength provides the information needed to characterize the astro-comb spectrum as a wavelength calibrator.

As shown in Fig. 5, the FPC finesse and phase deviation are found to vary slowly across the entire astro-comb spectrum, as measured with TRES. The FPC finesse is approximately 180 [Fig. 5(a)] and the phase deviation of the FPC relative to astro-comb lines is < 20 mrad [Fig. 5(b)]. Reflections from the back surfaces of the FPC mirrors used in these measurements produce the small, rapid phase variations observed in the figure. Despite antireflection-coatings the mirror substrates reflected 0.1% of the incident power. A simple model with realistic parameters reproduces the observed phase variations and allows systematic effects associated with these reflections to be eliminated. In future work, slight wedging (for example, 0.5°) of the substrates will eliminate these variations. After fitting the phase deviation with a sixth order polynomial, we compared our result with the phase deviation derived from a model of air dispersion [18] and the mirror phase delay measured with a white light interferometer (Fig. 6). The two different methods agree well, given systematic limitations to the white light interferometer measurement and resultant effect on the model fit.

Using these *in-situ* measurements of the FPC finesse and the phase deviation, along with the source comb intensity variation and resonant FPC transmission, we can determine the suppression of unwanted source comb lines and the resultant frequency shifts of astro-comb line centers. The variation with optical frequency of source comb line intensity and resonant FPC transmission [*T _{m}* ·

*I*in Eq. (1)] is measured by offsetting the diode laser frequency from its nominal value by the source comb repetition rate, and thereby tuning the FPC to resonance with a neighboring (and normally suppressed) source comb mode. The maximal transmission of each transmitted source comb line and thus

_{m}*T*·

_{m}*I*is measured; and then the procedure is repeated for all normally suppressed source comb modes. We find that the variation of

_{m}*T*·

_{m}*I*is within our measurement uncertainties (< 1%). Evaluating Eq. (1) for a single astro-comb line and all of its neighboring (suppressed) source comb lines within one spectrograph resolution element, and determining the net transmitted spectrum weighted by the frequency offset from the central astro-comb frequency, allows us to evaluate the effective line center shift in the astro-comb spectrum. For moderate finesse and

_{m}*f*≫ Δ

_{r}*/*

_{m}*F*, the systematic shift

*δf*of an astro-comb line of frequency

_{m}*f*, as measured on a spectrograph, can be approximated as

_{m}where the term in parentheses parameterizes the mean suppression of nearest source comb side modes, the first term in square brackets results from the difference in source comb intensity between the upper and lower side modes, and the second term in square brackets results from FPC phase deviation leading to asymmetry of the comb lines relative to the FPC mode. The frequency shift in an astro-comb line centroid due to variations of neighboring source comb line intensities at the 1% level is thus ≈ 70 kHz or 7 cm/s. Note that we do not expect neighboring source comb lines to differ as much as 1%; also a Gaussian distribution of comb-line intensities will average away much of this astro-comb centroid frequency shift. Nonetheless, we expect that a double pass configuration through the same FPC [19, 20], moderately higher FPC finesse (*F* ≈ 500), a higher repetition rate for the source comb (*f _{r}* ≈ 5 GHz), or more precise measurements of line to line intensity variations will be required to assure 1 cm/s accuracy In practical astro-comb calibrators. In Fig. 7, we show the shift of the estimated systematic error in the wavelength calibration of TRES as a function of wavelength, when using our current astro-comb as the calibration reference, if the effect of nearest suppressed source-comb lines are not included in the fit model. With proper inclusion of such effects, residual uncertainty in the wavelength calibration is at the 1 cm/s level. Thus we conclude that while astro-comb line frequency shifts caused by mirror dispersion and resultant wavelength dependence of FPC finesse and phase deviation can be absorbed into the wavelength calibration if reproducibility is all that is desired, at the few cm/s level these corrections must be applied to the spectrograph calibration to achieve accurate stellar radial velocity measurements.

## 5. Conclusion

In summary, we have demonstrated an *in-situ* method to determine systematic shifts of astro-comb spectral lines due to imperfect suppression of source comb lines by the Fabry-Perót filter cavity (FPC). This method involves measurement of FPC finesse as well as the phase deviation of all astro-comb lines over a bandwidth of 1000 Å. Such measurements can be performed at a telescope with either an astrophysical spectrograph or a commercial optical spectrum analyzer of lower resolution. From the measured phase deviation, the dispersion of the FPC is derived and found to be consistent with that calculated from the intra-cavity air dispersion and the group delay dispersion of the FPC mirrors. Applying resultant corrections for shifts in the centroid of astro-comb spectral lines allows us to model the astro-comb spectrum as measured on an astrophysical spectrograph with accuracy to better than 0.1 MHz, i.e., 10 cm/s for measurement of a stellar radial velocity. Such high accuracy is important for many applications of astro-comb wavelength calibrators, such as the search for habitable exoplanets, direct measurement of the expansion of the universe, and searches for a temporal variation of physical constants.

## References and links

**1. **C. Lovis, M. Mayor, F. Pepe, Y. Alibert, W. Benz, F. Bouchy, A. C. M. Correia, J. Laskar, C. Mordasini, D. Queloz, N. C. Santos, S. Udry, J.-L. Bertaux, and J.-P. Sivan, “An extrasolar planetary system with three neptune-mass planets,” Nature **441**, 305–309 (2006). [CrossRef] [PubMed]

**2. **A. Sandage, “The change of redshift and apparent luminosity of galaxies due to the deceleration of selected expanding universes,” Astrophys. J. **136**, 319–333 (1962). [CrossRef]

**3. **A. Loeb, “Direct measurement of cosmological parameters from the cosmic deceleration of selected expanding universes,” Astrophys. J. **499**, L111–L114 (1998). [CrossRef]

**4. **T. Udem, R. Holzwarth, and T. W. Hansch, “Optical frequency metrology,” Nature **416**, 233–237 (2002). [CrossRef] [PubMed]

**5. **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**, 839–847 (2007). [CrossRef]

**6. **P. O. Schmidt, S. Kimeswenger, and H. U. Kaeufl, “A new generation of spectrometer calibration techniques based on optical frequency combs,” in Proc. 2007 ESO Instrument Calibration Workshop (ESO Astrophysics Symposia series, Springer, in the press)

**7. **C. Araujo-Hauck, L. Pasquini, A. Manescau, T. Udem, T. W. Hänsch, R. Holzwarth, A. Sizmann, H. Dekker, S. D’Odorico, and M. T. Murphy, “Future wavelength calibration standards at ESO: the laser frequency comb,” ESO Messenger **129**, 24–26 (2007).

**8. **S. Osterman, S. Diddmas, M. Beasley, C. Froning, L. Hollberg, P. MacQueen, V. Mbele, and A. Weiner, “proposed laser frequency comb-based wavelength reference for high-resolution spectroscopy,” in Proc. SPIE, 6693, G1 (2007).

**9. **C.-H. Li, A. Benedick, P. Fendel, A. Glenday, F. Kärtner, D. Phillips, D. Sasselov, A. Szentgyorgyi, and R. Walsworth, “A laser frequency comb that enables radial velocity measurements with a precision of 1 cm s^{−1},” Nature **208**, 610–612 (2008). [CrossRef]

**10. **T. Steinmetz, T. Wilken, C. Araujo-Hauck, R. Holzwarth, T. W. Hansch, 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**, 1335–1337 (2008). [CrossRef] [PubMed]

**11. **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**, 57 (2008). [CrossRef]

**12. **C. Lovis, F. Pepe, F. Bouchy, G. L. Curto, M. Mayor, L. Pasquini, D. Queloz, G. Rupprecht, S. Udry, and S. Zucker, “The exoplanet hunter HARPS: unequalled accuracy and perspectives toward 1 cm s^{−1} precision,” Proc. SPIE **6269**, 62690P1–62690P23 (2006).

**13. **S. Udry, X. Bonfils, X. Delfosse, T. Forveille, M. Mayor, C. Perrier, F. Bouchy, C. Lovis, F. Pepe, D. Queloz, and J.-L. Bertaux, “The HARPS search for southern extra-solar planets. xi. super-earths (5 and 8 M⊕) in a 3-planet system,” Astron. Astrophys. **469**, L43–L47 (2007). [CrossRef]

**14. **M. J. Thorpe, R. J. Jones, K. D. Moll, J. Ye, and R. Lalezari, “Precise measurements of optical cavity dispersion and mirror coating properties via femtosecond combs,” Opt. Express **13**, 882–888 (2005). [CrossRef] [PubMed]

**15. **A. Schliesser, C. Gohle, T. Udem, and T. W. Hänsch, “Complete characterization of a broadband high-nesse cavity using an optical frequency comb,” Opt. Express **14**, 5975–5983 (2006). [CrossRef] [PubMed]

**16. **G. Furesz, “Design and application of high resolution and multiobject spectrographs: Dynamical studies of open clusters,” Ph.D. thesis, University of Szeged, Hungary (2008).

**17. **A. Yariv and P. Yeh, *Photonics* (Oxford University Press, Oxford, 2006).

**18. **P. E. Ciddor, “Refractive index of air: new equations for the visible and near infrared,” Appl. Opt. **35**, 1566–1573 (1996). [CrossRef] [PubMed]

**19. **M. S. Kirchner, D. A. Braje, T. M. Fortier, A. M. Weiner, L. Hollberg, and S. A. Diddams, “Generation of 20 GHz, sub-40 fs pulses at 960 nm via repetition-rate multiplication,” Opt. Lett. **34**, 872 (2009). [CrossRef] [PubMed]

**20. **T. Steinmetz, T. Wilken, C. Araujo-Hauck, R. Holzwarth, T. W. Hänsch, and T. Udem, “Fabry-Perót filter cavities for wide-spaced frequency combs with large spectral bandwidth,” Appl. Phys. B **96**, 251 (2009). [CrossRef]