1State Key Laboratory of Precision Measurement Technology & Instruments, Department of Precision Instrument, Tsinghua University, Beijing 100084, China
2Department of Physics, 366 Le Conte Hall MS 7300, University of California, Berkeley, California 94720, USA
Resolution and bandwidth are critical for cavity-enhanced dual-comb spectroscopy (CE-DCS). Here, we pioneer an adaptive approach in CE-DCS to improve the broadband as well as the resolution. Postcorrections to dual-comb interferograms adaptively compensate the relative phase jitters of the optical frequency combs and result in both a mode-resolved spectral resolution and a signal-to-noise ratio of 440:1 in 1 s. Meanwhile, an adaptive comb-cavity locking scheme exploits more than 90% of the comb modes, covering (10 THz) at . For a single dual-comb interferogram, more than 40,000 comb teeth spaced by 250 MHz are measured in less than 7.5 ms, contributing to a noise equivalent absorption per spectral element of . This adaptive cavity-enhanced dual-comb spectroscopy technique provides an attractive spectroscopic tool that may be utilized in trace-gas sensing, breath and cancer analysis, and engine combustion diagnosis.
1. INTRODUCTION
Cavity-enhanced optical frequency comb spectroscopy (CE-OFCS) marries broadband spectrum with high sensitivity, and thus permits spectroscopy in cold molecules probing [1], chemical kinetics analysis [2], flame diagnosis [3], breath analysis [4], and remote sensing [5,6]. It requires not only the perfect coupling of the comb into a cavity but also an effective technique for resolving and detecting the comb elements. Even though their analogous equidistant modes provide the very intuition of their integrations, matching the two degrees of freedom (DOFs) of the comb, namely, the repetition rate () and the carrier envelope offset (), with the only parameter of the cavity, the free spectral range (FSR), can be challenging. A two-point locking scheme [7,8] intuitively fits the two DOFs of the comb to the cavity but has defects in increasing complexity and in limiting the coupled optical power and bandwidth, since, in some implementations [9], modulation is applied directly to the comb light. The one-point locking scheme [10–12], which inversely locks the cavity to the comb, circumvents these problems. But previous practices have imperfections in nonconstant transmitting [12], showing incompetent stability for the dual-comb technique. Steady coupling is realized in a recent work [10], yet a short length cavity is used for portability, suffering from a large mode-spacing of 2 GHz and multiple steps of sweeping of cavity length for coupling every adjacent comb tooth. For measuring the transmitted comb light, various configurations are attempted, such as the cavity ring-down [13], employing dispersers [14–16], the Vernier methods [17–19], the mechanical Fourier transform spectrometers (FTSs) [9,20], and dual-comb spectroscopy (DCS) [10,21,22]. Among these, the DCS potentially stands out in subgigahertz resolution and submicrosecond detecting speed [23–27]. The debut of the CE-DCS [22], however, fails to be superior in spectral resolution, which is 4.5 GHz and is still surpassed by the mechanical FTS ones. This is, to a large extent, due to the deteriorated mutual coherence caused by the phase fluctuations between the two combs [28], only permitting a resolution 10 times worse than the ideal value. Locking the combs to two optical frequencies can be effective solutions [10,27], but there is a sacrifice in cost and complexity. Another implementation uses the combs generated by electro-optical modulating of the continuous-wave (CW) lasers, revealing less complexity and even higher mutual coherence [21]. But these combs emit a narrow spectrum spanning no wider than (50 GHz), which is insufficient for the simultaneous analysis of multiple spectral bands and molecules. The dawn of the adaptive methods [29–33], which compensate the distorted interferograms (IGMs), enables simple configurations as well as high resolution. Up till now, no such an adaptive approach with the cavity-enhanced setup has happened.
In this paper, we demonstrate an adaptive cavity-enhanced dual-comb spectroscopy (ACE-DCS) that realizes high spectral resolution, high sensitivity, and broad spectral bandwidth. The adaptive postcorrections of the IGMs are implemented, achieving a signal-to-noise ratio (SNR) of 440:1 in 1 s and a mode-resolved resolution (250 MHz), and eliminating the tight-locking of the comb to state-of-the-art lasers. Meanwhile, an adaptive comb-cavity coupling scheme, including the optimizations of the two comb parameters and the Pound–Drever-Hall (PDH) technique, tightly locks the cavity length to adaptively track the variation of the comb. These bring about better coupling stability, a broader spectrum (), and, to the best of our knowledge, the largest number () of spectral elements among the CE-OFCS techniques. This ACE-DCS technique potentially facilitates precision spectroscopy, including atmospheric trace-gas observation, rapid breath analysis, and dynamic combustion diagnosis.
2. PRINCIPLE
In our ACE-DCS setup, we demonstrate an adaptive comb-cavity coupling. Apart from tuning the and of the comb for optimizing the coupling bandwidth, a CW laser is added as an intermediary that adaptively locks the cavity length and permits robust and constant comb-to-cavity coupling.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
The first step in this one-point coupling approach is to optimize the repetition rate of the comb. Notably, when sweeping the cavity length, for the two cases that the is greater [Fig. 1(a)] or lower [Fig. 1(b)] than the FSR, this envelope moves in opposite directions and transmits different portions of the comb spectrum. Such a characteristic provides a criterion for adjusting the FSR. When sweeping the cavity, if the current locates within the dithering range of the FSR, two adjacent transmitted peaks of the comb that have horizontally symmetrical shapes should be found, as shown in Fig. 1(d). Otherwise, as shown in Fig. 1(c), the needs further adjustment. Next, tuning the value of until observing three adjacent transmission peaks and the two symmetrical ones that lie on the sides have equal intensities, as shown in Fig. 1(e), optimizes the transmission power of the center peak.
Figure 1.Adjustment of and . Cavity filtering on comb when (a) and when (b) . CC, current filtering center; NC, next filtering center. Cavity transmission when (c) is not appropriate, (d) is appropriate and (e) both and are appropriate.
Then, the cavity length should be locked to permit the highest transmitted power of the comb, the in Fig. 1(e). Previous one-point coupling adopts the dither-locking method, which, unfortunately, is incapable of a tight-locking performance due to its large response latency. Also, unwanted phase shift to the coupled comb light is introduced because the locking point locates at the side of the transmission peak rather than at the highest point. To overcome these, we innovate the approach by using a CW laser to lock the cavity length adaptively through a sideband PDH technique. The principle is that, on the one hand, the output of the CW laser is phase-modulated by an electro-optical modulator (EOM). Then the PDH technique generates error signals using the carrier frequency and one sideband, which are reflected by the mirror, and locks the rest of the sidebands into the cavity. On the other hand, the light of the CW laser without the modulation is coupled with the comb, their beat note is detected, and frequency is discriminated and used for feedback controlling the cavity length, locking the beat note at the modulating frequency of the EOM. Analogous to the frequency-agile rapid scanning spectroscopy (FARS) [34], the locked comb tooth shares the same frequency with the locked sideband, and thus, the comb is coupled into the cavity with a spectral center close to the frequency of the CW laser. As a supplement to the optimization of the comb parameters in the previous steps, tuning the frequency of the CW laser to the spectral center of the comb can further widen the coupled bandwidth. Because the “equal intensities” used in previous implementations may be relatively imprecise, additional fine-tuning of the could further optimize the coupled spectrum. Once such adjustments are made, the steady coupling status can be maintained because the cavity length is constantly traced to the same comb tooth. That is to say, the length of the cavity is tuned accordingly to follow the jitters of the and the .
For obtaining a mode-resolved resolution, a postcorrecting method is adopted. In DCS, the jittering between each pair of comb teeth can be equivalent to the difference of the repetition rate, meaning insufficient stability for resolving each comb element. Our approach can remedy this through retrieving the phase jittering of the two combs using two CW lasers and applying adaptive corrections to the IGMs [32,33]. In our setup, the CW laser used for the comb-cavity coupling can also work for this adaptive correction, and only one additional CW laser with a different operating wavelength is needed. These two CW lasers beat with the two combs, and the four generated beat notes are used for the postcorrection algorithm. The details of this algorithm are described in the Appendix A. Through the mixing and resampling of the original IGMs, the corrected IGMs are free from the effects of both the jitters and the jitters.
Thanks to the stable transmission maintained by the adaptive comb-cavity coupling and the error-free IGMs produced by the postcorrections, high-quality dual-comb measurements can be continuously obtained. Subsequently, applying the direct time domain averaging using multiple IGMs, higher SNRs can result.
3. EXPERIMENTAL SETUP
The experimental setup, shown in Fig. 2, consists of a measurement part and an acquisition part. In Fig. 2(a), the output of an external cavity diode laser (ECDL) (Newport, TLB-6728) was divided into three equal power portions. The first portion was phase-modulated by an EOM (Newport, 4002) that generated sidebands for the PDH locking at a detuning frequency of 20 MHz. A phase shifter adjusted the polarity of the PDH error signal to move the locking point of the proportion-integration (PI) controller (Newport, LB1005) to one of the sidebands. Another output portion of the ECDL beat with Comb1 (Menlosystems, FC-250-1500). When the ECDL was locked to the cavity, the frequency of the ECDL followed with the cavity length, which in turn varied its beating frequency with the comb. Using a phase discriminator (Menlosystems, DXD 200) and a PI controller, we feedback-controlled the piezo actuator (PZT) mounted on the cavity mirror and locked the beating to the same frequency of the EOM (20 MHz). All the RF signals, including the modulating source of the EOM and the sampling clock of the data acquisition device, were synchronized to the same rubidium clock (Symmetricom, 8040C), which was also the reference of the combs. Thus, the beating comb tooth shared the same frequency with the ECDL sideband and could be coupled into the cavity. The last portion of the ECDL output was used for the acquisition part and will be explained later. The enhancement cavity, as shown in Fig. 2(a), is approximately 600 mm long, corresponding to an FSR of about 250 MHz, close to the repetition rate of the combs. The reflectivity of the cavity mirror is higher than 99.95% within the whole spectral range of the comb, denoting a finesse of around 3000. Four invar rods supported the cavity for stabilizing the cavity length, and a sandwich of heater and insulation cotton maintained the temperature with a thermostat (Wavelength Electronics, PTC-2.5 K). There are three inlets and outlets along the cavity, which were used for exchanging the gas and monitoring the inside pressure.
The acquisition part is shown in Fig. 2(b), wherein one of the two CW lasers that generated the reference signals is the rest portion of the ECDL, as noted above. The other CW laser [CW2, shown in Fig. 2(b)] is a narrow linewidth laser (Redfern Integrated Optics, 1534.233 nm). Because of the low output power of the comb around the wavelength of the CW2, we used two fiber Bragg gratings (FBGs) to improve the SNR of the beat notes. The comb power around the ECDL frequency was sufficient for omitting the use of FBGs. Since the beat note between Comb1 and the ECDL was tightly locked to 20 MHz with subhertz jittering (in the measurement part), this beat signal can be numerically constructed, and is needless for being recorded. Compared to Ref. [33], the required acquisition channels are reduced by one, and only a four-channel acquisition board (NI 5761, 14 bits, 250 MS/s) was used for synchronously digitizing the three beat notes and the dual-comb IGMs. Although in Ref. [32], two balanced detectors are used for simultaneously recording four beat notes using just two acquisition channels, only three channels are needed. But, in that setup, additional numerical filtering is required, and the effective dynamic range of the each digitizing beat note is reduced. We adopted a bandpass filter (17–32 MHz) for recording the DCS signal, which could help avoid frequency aliasing and increase the SNR. The maximum recording time is 1 s, limited by the size of the onboard memory. Since the acquisition card was populated with a field-programmable gate array (FPGA), it can be used for real-time processing of the IGMs and circumvent the limited acquisition time.
4. RESULTS
According to the adjustment procedure mentioned in Section 2, we recorded the transmitted optical power when sweeping the cavity length, as shown in Figs. 3(a)–3(c); the corresponding spectrum is shown in Fig. 3(e). In case that the was not fully matched with FSR, only a bandwidth of and an optical power of 0.5 mW could be coupled. After adjusting the two parameters of the comb until we got two peaks of identical height, as shown in Fig. 3(c), we further fine-tuned the to find the widest coupled bandwidth, shown as the light blue in Fig. 3(e). The determined optimal values of the and of the Comb1 were 250 and .
Figure 3.Optimization of and and coupled spectrum. Cavity transmission in the experiment when (a) is not appropriate, (b) is appropriate, and (c) and are appropriate. The nonuniform intervals between adjacent two peaks appearing in (b) and (c) were due to the nonlinearity of the piezo movement. (d) Cavity finesse and corresponding frequency detuning of FSR. (e) The spectra of the original comb and the transmitted light in cases of (a), (c), and fine-tuning of .
Thanks to the fine-tuning of comb parameters and the tight locking of the ECDL, the final optimized optical power was 3.6 mW, which is more than 17% of the incident comb power; the corresponding coupled bandwidth stably exceeded . Thus, a spectral bandwidth of over 10 THz is available for DCS measurement. Since the FSR of the cavity is the same as the of the combs, almost every comb tooth was coupled within the entire bandwidth of the comb output, with a total number of more than 40,000. The coupled bandwidth can be potentially further widened if the comb is spectrally expanded using nonlinear fiber optics. These results convincingly attest to the feasibility of this one-point cavity-comb coupling. It is fair to say the coupled comb intensity and the bandwidth were degraded compared to the original comb, the reason for which can be complicated. First, due to the dispersion of the cavity mirror, the modes of the cavity are not strictly equidistant, as indicated by the cavity finesse shown in Fig. 3(d). Thus, even if the comb modes around the frequency of the CW1 laser have perfect alignment with the cavity, modes at frequencies away from the CW1 laser may reveal mismatching with the cavity. Second, the response bandwidth of the PZT that is mounted on the cavity mirror is typically small and may raise instability of the PDH locking, causing the degradation of the coupling efficiency.
To demonstrate this ACE-DCS spectrometer, we measured rovibrational transitions of the gas sample around , which covers four spectral bands, namely, the (), the (), the (), and the (). The gas sample had a purity of more than 99.9% and was detected at a pressure of around 420 Pa and a temperature of 298 K. According to the passband of the filter and the determined parameters of Comb1, as mentioned above, we set the and of Comb2 to be 250.00014 and 20 MHz. Thus, according to Ref. [27], the frequency of the IGMs located within 20 to 25 MHz, and the measurement rate was 140 Hz. Figure 4(a) is the spectrum after applying Fourier transform to a measured IGM with a recording time of 1 s. As shown in the inset of Fig. 4, the different regions of the spectrum all reveal clear features of absorption transitions. The “overshoot” baselines, shown around and , are due to the misalignment between the cavity modes and the comb teeth caused by the mirror dispersion, which commonly happens for broadband cavity-enhanced frequency comb spectroscopy and has proved to be compensated when taking the dispersion parameters into the fitting process [10,35]. Figures 4(b)–4(d) and the zoomed-in spectrum, where each comb line can be distinguished, confirm the effectiveness of the adaptive postcorrection algorithm and a mode-resolved resolution of this ACE-DCS system.
Figure 4.Spectral measurement of ACE-DCS system. (a) The spectrum of sample around ; blue lines, the measured spectrum after Fourier transforming the IGM, exhibiting each individual comb mode; red lines, the fitted spectrum using the highest points of each comb mode, revealing clear absorption features. (b), (c), and (d) Mode-resolved spectrum in different scales, spanning , , and , respectively. (e) Transmittance and phase around ; (f) transmittance of the P(16) line of the () band and the Voigt fitted result and the line profile provided by HITRAN 2016.
After truncating the IGM of 1 s into 140 individual IGMs, we averaged them in the time domain and applied Fourier transform. Then we extracted the baseline and got the transmittance and the phase, both of which reveal good quality, as shown in Fig. 4(e). Figure 4(f) shows the observed transmittance for the P(16) line of the () band and its fitted curve using the Voigt line shape model. According to this plot, we calculated the SNR using the standard deviation of the transmittance, where no transition lines appear [10,22]. This calculation of SNR is intensity-normalized, and our system obtained a good SNR of higher than 440:1 (in 1 s). Compared to the theoretical line profile using the parameters provided by the HITRAN database [36], as shown in Fig. 4(f), we estimated the actual absorption length to be about 560 m, and the corresponding enhancement coefficient is higher than 900. Given the number of spectral elements of about 40,000 (span/resolution), the noise equivalent absorption (NEA) per spectral element [37] is determined to be , showing more than a tenfold improvement over the latest CE-DCS setup [10]. Compared to the cavity ringdown spectroscopy [38], this ACE-DCS system also reveals a commensurate sensitivity and outstands because of its broadband measurement ability.
5. CONCLUSION
In summary, we developed an adaptive cavity-enhanced dual-optical comb spectroscopy that enables CE-DCS with mode-resolved resolution, broad bandwidth, and high sensitivity. The adaptive corrections of the IGMs achieve a resolution of 250 MHz with a compact and simplified setup. The adaptive comb-cavity coupling allows tighter locking, wider bandwidth, and an increase in the number of the total spectral elements. Unlike the two-point locking scheme [7,8], our scheme locks the cavity to the comb. Thus, no strict vibration and thermal control are required for the cavity, and the comb could be locked to highly stable frequency standards, such as atomic clocks, which in turn permits excellent stability and traceability. Compared with the previous adaptive DCSs [30,32], the measurement sensitivity is improved more than 2 orders of magnitude because of the enhanced cavity. Given that no additional CW lasers are needed, since the lasers used for measuring the dual-comb jitters can also be used for coupling the cavity, very little increase in complexity is introduced. Considering the advantages of the DCS, which are the elimination of moving parts and the high measurement speed, this ACE-DCS technique is expected to benefit many precise spectrometry applications. It is especially true for the investigation of atmospheric observation and combustion diagnosis, where it requires not only high resolution and sensitivity but also the need for dynamic and multiplex analysis. By adopting cavity mirrors with higher reflectivity, this ACE-DCS system can enable broadband spectral measurements for the applications that require extremely high sensitivity, such as trace-gas detection and radiocarbon chronology. These advantages could mean they will surpass and replace CW laser-based spectral measuring techniques and show huge potential to promote fundamental science.
APPENDIX A: POSTCORRECTION ALGORITHM
The details of the postcorrection algorithm are illustrated in Fig. 5 and have two major steps: mixing and resampling [39].
Figure 5.Diagram for the postcorrection algorithm. Clk., clock; Sig., signal.
First, the mixing step is used to neutralize the mutual jitter of the difference of , represented by . Specifically, the two CW lasers produce beat notes at the frequencies around their adjacent comb modes: the numbers of and modes of Comb1 and those of and of Comb2, respectively. In our configuration, we chose a difference of repetition rate () to the one that makes the value of equal to and equal to , helping avoid frequency aliasing. Besides, we actually recorded only three beat notes (Beats 2, 3, and 4 in Fig. 5) and digitally constructed the last one (Beat 1 in Fig. 5). The recording of these beat notes was synchronized with that of the dual-comb IGMs. In terms of processing, we extracted the phase of the beat notes derived between the two combs and the same CW laser ( and for CW1, and for CW2) and calculated the difference between the two ( and ). These two phases represent the phase change between the two combs at the frequencies of the CW1 and CW2 lasers and can be expressed as Eqs. (A1) and (A2). Then, as shown in Fig. 5, we produced a using the value of and add (), where is an integer. This step is to get a so that the frequency of the lies outside the frequency of the IGMs. Subsequently, we mixed the IGMs with the and applied low-pass filtering to the product. Because the contains the jitter of , the shown in Eq. (A4), the mixed IGMs are free from the influence of and only consist of multiples of the repetition rate difference and its jitter ( and ): Second, the resampling could eliminate the jitter of the repetition rate on the IGMs. Specifically, we took the , shown in Eq. (A3), as a reference coordinate and resampled the IGMs by interpolation and at a specified interval on this coordinate. To avoid the loss of accuracy and the increase of the data size, the resampling interval was assigned to one that maintained the lengths of the resampled IGMs close to the original ones. After that, the resampled IGMs were free from the effects of both the jitters and the jitters.