Relative timing jitter compression in a Fabry-Perot cavity-2 assisted free-running dual-comb interferometry

9 Abstract . Dual-comb interferometric systems with high time accuracy have been realized for various applications. 10 The flourished ultra-low noise dual-comb system promotes the measurement and characterization of relative timing 11 jitter thus improving the time accuracy. With optical solutions, introducing an optical reference enables 10 5 harmonics 12 measurements therefore breaking the limit set by electrical methods, and nonlinear processes or spectral interference 13 schemes were also employed to track the relative timing jitter. However, such approaches operated in the time domain 14 either require additional continuous references or impose stringent requirements on the amount of timing jitter. This 15 paper proposes a scheme to correct the relative timing jitter of a free-running dual-comb interferometry assisted by a 16 Fabry-Perot cavity (F-P cavity) in the frequency domain. With high wavelength thermal stability provided by the F-P 17 cavity, the absolute wavelength deviation in the operating bandwidth is compressed to less than 0.4 pm, corresponding 18 to a sub-picosecond sensitivity of pulse-to-pulse relative timing jitter. Also, Allan deviation of 10 –10 is obtained under 19 multiple coherent averaging, which lays the foundation for mode-resolved molecular spectroscopic applications. The 20 spectral absorption features of hydrogen cyanide gas molecules at ambient temperature were measured and matched 21 to the HITRAN database. Our scheme promises to improve new ideas in sensitive measurements of relative timing


Introduction
Born from optical frequency comb (OFC) technology [1][2][3], the dual-comb system is made up of two OFCs with a slight difference in repetition frequency (f rep ) that inherits high stability and frequency accuracy better than 10 -18 supported by each OFC at optical frequency [4,5], which promotes the research of optical metrology applications [6][7][8][9][10][11][12].By establishing the mutual coherence between dual combs, the high level of temporal consistency promotes a range of applications, including absolute distance measurement based on time-of-flight methods [13], remote synchronization as well as calibration between combs [14], and on-chip clock distribution networks with femtosecond precision [15], to name a few.Specifically, two sequences of stable ultra-short pulses emitted from two independent OFCs walk off temporally with a fixed interval of δτ = 1/ f rep -1/(f rep +Δf rep ) ≈ Δf rep /f 2 rep from each other, where Δf rep is the offset repetition frequency between two combs, forming the regular asynchronous optical sampling (ASOPS) in the time domain.Such a down-sample approach enables optical carrier frequencies that would be difficult to detect directly to be captured by lower-cost, more mature radio frequency (RF) devices.
In such scenarios, the time axis is linearly magnified through the equivalent sampling by, for example, more than 10 5 times.This allows recording the optical signal in real-time at the RF domain, providing a powerful tool for measuring ultrafast phenomena within microseconds.In high-precision time-domain applications such as equivalent time sampling and dual-comb spectroscopy (DCS), the relative timing jitter requires rigorous measurement and subsequent reduction to guarantee a time synchronization between two pulse trains.Typically, a mode-lock fiber laser-based dual-comb ranging system with a timing jitter lower than the sub-femtosecond results in a high resolution of the sub-micrometer [13].For dual-comb-based spectroscopic applications, Shi et al.'s work [16] numerically derives an amplified spontaneous emissions (ASE) noise-induced quantum-limit relative timing jitter of 0.82 fs to introduce a relative linewidth of 650 Hz into the free-running single cavity DCS and finally set a 65-MHz spectral uncertainty at the measurement.
Timing jitter refers to the degree to which the optical pulse deviates from the perfect periodic position.For mode-locked lasers, the timing jitter was set by the intracavity ASE noise in two ways given by [17] : where f is the offset frequency, D T and D ω are the diffusion coefficients for ASE-induced group delay velocity and center wavelength shift of pulse trains, respectively, D is half the net dispersion of the cavity，    is the decay time of frequency fluctuations.For the first term in Eq. ( 1), ASE noise is coupled directly into the timing jitter; the second term is known as the Gordon-Haus jitter, where ASE noise is coupled to the timing jitter through dispersion-affected center frequency fluctuations [18].Besides the ASE noise-induced timing jitter, another factor lies in the coupling from intensity noise by several mechanisms, including slow saturable absorber, Kramers-Krönigrelated phase changes, and other nonlinear processes [19].
The measurement of relative timing jitter between dual combs can be roughly divided into two approaches: electrically and optically.In electrical schemes, fast photodetector (PD) is usually used to detect f rep and their higher harmonics n•f rep .Usually n > 100, then n•Δf rep is obtained by electrical mixing and finally analyzed by commercial signal analyzers [20].However, the optical signal is limited by the detection speed of the PD and electrical equipment, and the noise characteristics of this electrical reference set an inherent lower boundary for measurement sensitivity to satisfy the demands for phase noise characterization [21].Excess intensity-phase noise coupling also degrades the measuring accuracy of the noise of interest [22].GHz bandwidth limitations can be overcome by using the optical method, thus achieving n > 10 5 th harmonic measurements.To overcome the bottleneck of electrical approaches, the optical method either bypasses the noise-limited electrical reference or records the traces of relative timing jitter in realtime, employing interaction between dual combs, thereby improving the detection sensitivity to the femtosecond level [16,[23][24][25].For example, nonlinear ASOPS based on a sum-frequency generation process records the relative timing jitter to the idler component, which provides the ability to track the relative timing jitter directly and achieves the measurement capability of the attosecond level by equivalently amplifying the time axis [16].However, the occurrence of nonlinear processes sets a low threshold for power.In addition, such solutions inherently demand that the deviation of the timing jitter be less than the duration of the pulse.By introducing an optical reference, a multi-heterodyne detection method can derive a time-varying function of the offset repetition frequency while maintaining the minimal intensity-phase noise coupling, which can subsequently be analyzed for phase noise or calibrated in post-processing [23,24].The uncorrelated relative timing jitter of the free-running dual-comb system is measured to achieve the sub-attosecond detection sensitivity.Another reference-free solution borrows the idea of spectral interference to ingeniously transfer the undetectable timing jitter with fs-level in the time domain to the frequency domain [25].To this end, a dispersive Fourier transform (DFT) process [26] was carried out to map the interference patterns containing relative timing jitter back to the time domain with a mapping factor of 2πβL, where βL is the group delay dispersion of the system provided by the dispersive fiber [25].A high-speed PD and oscilloscope finally recorded the real-time optical signal stream.Through digital signal processing, the timing jitter between roundtrips can be calculated from interferograms (IGMs), and the relative timing jitter of a counterpropagating allnormal dispersion dual-comb fiber laser was measured to verify the scheme's feasibility.In this way, the measurement frame rate is the f rep of the laser, and the detection sensitivity was also increased by 5 times compared to the optical-referenced heterodyne method.However, the accuracy of this method is highly dependent on the dispersion linearity of the DFT process.
Specifically, the residual dispersion of the third or higher order makes the frequency-to-time mapping no longer linear, resulting in the distortion of the calculated relative timing jitter.Besides that, spectral interference patterns of fast oscillations in time-domain waveforms must be captured using large-bandwidth detection equipment, while smaller Fourier frequencies also rely on that.

Materials and Methods
In this work, we exploit the frequency-time mapping relationship to convert the study of relative timing jitter in the time domain to be performed in the frequency domain, inspirited by [25].The relative timing jitter is reflected as deviations in the RF spectrum for a dual-comb interferometry.
In turn, tiny deviations in the spectrum can be used to trace the relative timing jitter periodically.
The essential prerequisite for realizing the above time-frequency correspondence is that another degree of freedom, i.e., the relative carrier-envelope offset frequency (Δf ceo ) noise penalty is effectively suppressed or calibrated in real-time [27,28].In our previous work [29], it was demonstrated that a passive phase-shifted fiber Bragg grating (PS-FBG) can be used to calibrate spectral distortions caused by Δf ceo penalties owing to its ultra-narrow notch filter bandwidth (~100 MHz).Here, we compare the effect of locking mode on Δf rep under three different operating conditions, i.e., electrical-referenced locking, free-running, and Fabry-Perot (F-P) cavity-assisted locking, illustrated in Fig. 1.It should be emphasized that the v-axes in Fig. 1 have no real physical meaning.In reality, ASOPS implements pulse-to-pulse sampling of the probe light to the signal light, and the optical carrier signal at the effective time is linearly stretched into a laboratory time to yield the down-sampled RF waveform within the Nyquist bandwidth.The corresponding optical frequency can be inverted using up-conversion.For the case shown in Fig. 1(a), the reference signal for f rep stabilization is usually provided by a commercial RF reference source, and the degree of stability depends on the reference quartz crystal inside the facility.There is a definite frequency transfer relationship between f rep (MHz) and the optical comb frequency (THz), and the transfer quantity N is usually larger than 10 6 [3].Although f rep is accurately locked in the RF domain, the accumulation of frequency instability (i.e., σ(ν)) causes instability in the optical frequency, which significantly deteriorates the accuracy as known by the elastic tape theory [30,31]  Figures 1(c) and (e) compare the equivalent locking of the Δf rep when OFCs are completely free-running, where neither frequency degree of freedom is interfered with.Specifically, Δf rep jitter leads to telescopic changes in the RF comb lines, while Δf ceo jitter leads to an overall shift in the combs, both of which broaden the combs over long time scales, resulting in degradation of resolution and distortion of the measurements [27].Using a PS-FBG allows Δf ceo jitter to be tracked in real-time [29,32], setting the stage for subsequent equivalent Δf rep locking.Generally,

Free-running comb 1
Free-running comb 2 π mutation the optical field of the Gaussian pulse trains output by the two mode-locked fiber lasers can be approximated separately as: where p is the power intensity, T pw denotes the pulse width, Φ denotes the initial phase, and the subscripts 1 and 2 of all of the above items are used to differentiate between OFCs 1 and 2. After passing through the coupler, the two optical fields interfere to generate an IGM train, and the interference signal can be approximated as follows: The third term is the mutual interference term between each other.Bringing Eqs. ( 2) and (3) into Eq.( 4), the beat frequency signal can be approximated by considering only the beat frequency results between the most neighboring combs: where A(v) is a slow-varying envelope function, from Eq. ( 5) the IGM can be seen as a Gaussianshape envelope, whose oscillating frequency is the carrier frequency difference between dual combs.PS-FBGs typically have notch filter bandwidths in the order of 100 MHz, which can be indicated by the grey spectral pattern in Fig. 1(c).Thus, the beat frequency of the adjacent neighboring teeth between dual combs can be extracted as a calibration signal expressed as: where n c is the number of comb modes at the notch center of the PS-FBG, and p is a constant.
From Eq. ( 5) and ( 6), both contain Δf ceo (t) in the phase term, which can be eliminated by reconstructing the raw IGM, thereby equivalently locking the comb at v(n c ), as viewed by the red arrow in Fig. 1(c).From the frequency domain point of view, this process can be interpreted as "moving" the frequency center to v(n c ).However, the residual timing jitter still leads to a deviation ).This is attributed to the accumulation of σ(v) on both sides, with v(n c ) as a virtual "fixed point".Next, an F-P cavity is introduced, and two reflecting end surfaces form a standing wave field inside the cavity, only frequency components satisfying certain phase conditions can output, forming resonance peaks near each transmitted wavelength, resulting in discrete comb spectral features [33].These discretely distributed resonance peaks provide excellent wavelength references, especially when combined with the stable structural features of the F-P cavity.The wavelength deviation at the very edge (dark green arrow in Fig. 1(e)) can be corrected and then transfer the wavelength stability uniformly to each line in the working bandwidth (light green arrow in Fig. 1(e)), equivalent to performing an optical frequency division operation [34][35][36] on the wavelength axis.In this way, wavelength deviation under free running is effectively eliminated (evolution of the red to green curve in Fig. 1(f)) and the relative timing jitter is significantly compressed.The commonly used oven-controlled crystal oscillator (OCXO) provides 10 -9 stability and accumulates frequency deviations on the order of mHz when transferred to optical frequencies.A commercially available F-P cavity has a nominal wavelength thermal stability of ±6 pm over a temperature range of -5 to 70°C.Dividing the above parameters by the conversion factor M of a dual-comb system and the number of teeth N line in the operating bandwidth, one can obtain that each comb line shares a frequency stability of 12 mHz.
As evidenced by the following results, the locking accuracy of the relative timing jitter using the assisted by F-P cavity is superior to that of the electrically referenced, thus enabling a freerunning DCS.The flowchart of data processing is given, the frequency distribution of the phase calibration signal under free running is calculated, and it is characterized by power spectral densities to show that Δf ceo penalty correction is achieved in both free-running and electrical reference locking cases.Next, the results of the wavelength calibration show that the wavelength accuracy corrected by the calibration is superior to the case of electrical reference, even under free running.The absolute wavelength deviations below 0.4 pm over the entire band correspond to a sub-picosecond sensitivity of pulse-to-pulse relative timing jitter.Calculations of the Allan deviation show that stability of 10 -10 can be achieved after multiple coherent averaging.Finally, the gas absorption spectra of hydrogen cyanide gas molecules were characterized using the freerunning dual-comb interferometry and matched to a HITRAN database, enabling mode-resolved molecular spectroscopy applications.

Experimental setup
We carried out free-running DCS to demonstrate the superiority of the F-P cavity-assisted scheme compared with the conventional electrical-referenced one.The schematic diagram of the experimental setup is illustrated in Fig. 2. Two identical homemade Erbium-doped mode-locked fiber lasers served as the dual-comb source with a f rep of 100 MHz.The lasers work in the stretchedpulse state by optimizing the intracavity net dispersion to near zero in the negative dispersion region [37], which helps to reduce the Gordan-Haus jitter as well [38] according to Eq. ( 1).The working bandwidth of both combs exceeds 130 nm for broadband spectroscopic applications as shown in Fig. 2(b).The passive fiber devices-aided reference-free post-processing scheme was organized in Fig. 2(a) to establish the mutual coherence between two OFCs without introducing any bulky negative feedback control mechanisms.Two OFCs with a common path are divided into three paths by a combination of optical couplers (OCs), one of which is used for probing while the other two perform calibration operations.Among them, one PS-FBG with a 118.3-MHz notch bandwidth measured using a fast-sweeping method [29] is set to correct the Δf ceo penalty between dual combs.Although the PS-FBG is a temperature-sensitive device, the core is to filter one line in a dual-comb interferometer.Therefore, even if the temperature change causes a drift of the notch center, the ability does not change and does not lead to a failure of the phase calibration only to set a different reference for post-processing.To reduce the environmental impact, we sealed the PS-FBG in a metal box.Different from other paths, the optical signal passes through the F-P cavity are re-divided into two paths, simultaneously recording the transmission spectrum using a commercial optical spectral analyzer (OSA) with a 0.02-nm resolution for followed wavelength The hydrogen cyanide gas cell was chosen as the demonstration.16-frame synchronous raw data were recorded.Limited by the memory depth of the ADC card, the time window for each frame recorded is 160 ms.After interpolating the measured interference signal, the gravity of the envelope is determined to be the center burst of the IGM [39].The time difference between two adjacent central bursts in each frame represents the average repetition difference between the IGM and the two adjacent interference cycles, denoted as ∆f rep .The time window of 1/2∆f rep is extended to both sides around the center burst of each IGM, and it is divided into a single processing unit for the execution of online post-processing procedures.The individual IGM window split here is not the same as the actual interference period because the repetition difference changes randomly, and the interval determined by the peaks is only the average value.(g) The spectrum of the signal after performing phase calibration of (a) using the time-varying phase difference plot in (e).(h) Spectrum after performing wavelength calibration and spectral shaping on (g) taking (c) as the reference.
The block diagram shown in Fig. 3 illustrates the post-processing process of the raw data.The broadband F-P cavity used for wavelength calibration is considered the "sample under test", and the signal detected is I F-P (t).As the dual-comb passes through, the comb-like spectrum produces a signature of the periodic bulge near the center burst of each IGM (inserted figure in Fig. 3(a)).
Performing a fast Fourier transform (FFT) directly yields a chaotic spectrum in which the RF comb lines are overlapped.This is due to the deterioration of coherence between two combs, which causes the RF comb lines to be broadened and distorted over a long time.Concurrently, the signal detected from the PS-FBG path is recorded as a phase calibration signal and is denoted as I PS- FBG (t), presenting a single-frequency signal covering the entire interference period (inserted figure in Fig. 3(b)), which suggests its 118.3-MHz notch filtering bandwidth can only filter out a single comb tooth.Performing Hilbert transform on two signals to obtain their phases as a function of time as well as their time-varying phase difference, shown in Figs.3(d) and (e) respectively.The effective window for interference to occur is much smaller than the period of the IGM determined by Δf rep , so the 2π ambiguity can be avoided when unwrapping the phase [40].As confirmed by Eqs. ( 4) and ( 5), the Δf ceo jitter is contained in both time-varying phase signals φ F-P (t) and φ PS- FBG (t).Thus, reconstructing the time-domain signal utilizing the phase difference Δφ(t) (Fig. 3(e)) can effectively correct the Δf ceo penalty, mathematically expressed as |H[I F-P (t)] |•exp[iΔφ(t)], where H denotes for the Hilbert transform.This operation is known as the phase calibration, and the output spectrum v DCS F-P (λ) after phase calibration is attained by performing FFT on the reconstructed waveform, where the superscript "DCS" indicates the detect facility.As viewed in Fig. 3(g), comblike transmission spectra of F-P cavity features are presented instead of a chaotic one (Fig. 3

(f)).
A transmission spectrum with a central wavelength of 1550 nm, a 60-nm bandwidth spanning, and a 100-GHz free spectral range (FSR) was obtained.Since both signals (Fig. 3(a) and (b)) experience the same fluctuation, this temporal consistency makes phase calibration reliable.
However, the retrieved spectra still suffer from wavelength deviations due to residual relative timing jitter.To this end, the spectrum of the F-P cavity v OSA F-P (λ) measured by the commercial optical spectrum analyzer (OSA) with a 0.02-nm resolution was taken as a wavelength reference (Fig. 3(c)), whose spectral profile was removed by measuring the background spectra simultaneously.Similarly, the superscript "OSA" indicates the detection facility.In this way, wavelength deviations can be eliminated through a resampling process, equivalent to reconstructing the original wavelength axis λ in Fig.The validity of the two calibration processes needs to be further checked under the free-running state.Figure 4 illustrates the distribution of frequencies extracted from I PS-FBG (t).We select frames The waveform exhibits a single-frequency oscillation at each time point, consistent with the conclusion that a single-frequency signal was obtained from I PS-FBG (t), as described in the previous section, and the signal power is concentrated near the center burst (defined in each IGM as the moment t = 0) and covers the entire interference period.The single-frequency waveform shows an approximately sinusoidal frequency distribution in all three visual IGMs which oscillate at about 1 kHz despite different interference periods.This indicates that the Δf ceo -dominated frequency fluctuations vary with time at 1 kHz around their respective fundamental carriers.Besides the main component, it is still possible to observe "power fluctuations" superimposed on a sinusoidal waveform with an oscillatory degree (reflected in rapid frequency changes) much lower than the principal fluctuation.These rapid frequency changes are inferred to be attributed to Δf rep jitter.where 430 IGMs were selected and calculated subsequently, as plotted in the purple region in Figs.

5(a) and (b).
From the calculated PSD and the integrated frequency curves, the Δf ceo of the dualcomb system exhibits similar frequency noise characteristics in both states.In other words, the validity in the free-running state can be guaranteed when performing the phase calibration.The correction for wavelength deviation through wavelength calibration was then performed and displayed in Fig. 6.The wavelength at the peak position of each resonance peak of the F-P cavity within the operating bandwidth (defined here as 1520 to 1580 nm) was obtained and compared with the reference wavelength measured by OSA to obtain the absolute wavelength @ 1520.9 nm @ 1550.0 nm deviation σ(λ p ), where λ p is the set of all resonant peak wavelengths.The deviation of all resonance peaks is recorded for fitting to obtain the calibration function, and the original wavelength axis is resampled to obtain the calibrated one, which can be corrected for σ(λ p ).We compared the σ(λ p ) for three calibration conditions: electrical reference (i.e., DPLL), free-running, and free-running followed by wavelength calibration.For each IGM, since the wavelength deviation is accumulated uniformly from the center v(n c ) to both sides, the intermediate spare intervals can be approximately described using the calibration function, even though the calibration function is accurate only at the resonance peak wavelength position.A V-shaped distribution is predicted due to the uniform accumulation, as verified by three curves in Fig. 6(a).For the quantitative description of σ(λ p ), the down-conversion factor M of the dual-comb system is required to construct the wavelength axis.
In the case of DPLL, Δf rep is locked but inaccurate.When reconfiguring the wavelength axis with a fixed Δf rep , it will inevitably introduce a significant wavelength deviation, especially away from the center.On the other hand, although Δf rep is drifting for the free-running case, it is possible to calculate accurately by determining the interval of the center burst of the IGM with sub-Hz level (Figs.4(a) to (c)).This is sufficient to restore the wavelength axis more accurately even though there is a deviation, which leads to the reduction in σ(λ p ) (orange curve) compared with the DPLL one (cyan curve), as shown in Fig. 6(a).With wavelength calibration, the wavelength deviation is significantly compressed even in free running (magenta curve).Figure 6(b) illustrates the magnified results of the magenta curve in Fig. 6(a).Statistically, σ(λ p ) is reduced by about two orders of magnitude, and the deviation at all resonance peaks less than 0.4 pm, which is even lower than the theoretical spectral resolution (0.8 pm).This is meaningful because it sets the sensitivity for measuring the relative timing jitter.Thanks to the high-frequency fractionation due to the large bandwidth of the F-P cavity, pulse-to-pulse relative timing jitter measurements can be determined with a sensitivity on the sub-picosecond level.
In dual-comb applications, simple but effective coherent averaging [41,42] is performed to enhance the signal-to-noise ratio (SNR) of the signal, as the averaging process results in a significant reduction of the Gaussian white noise introduced by the detector.Two representative resonance peaks were chosen, 1520.9 nm away from and 1550 nm near the center wavelength, and the Allan variance was calculated under the three calibration conditions, as plotted in Figs.5(c) and (d), represented as triangles and squares, respectively.After multiple coherent averaging, the free-running case exhibits lower deviations than the DPLL case, and the instability of wavelength can be as low as 10 -10 .And for the wavelength calibration, the away-from-center one performs better (the last point on the magenta curve in Figs.5(c) and (d)).This is because the larger the wavelength deviation, the more accurate the calculation will be when calculating the calibration function.Conversely, when the wavelength deviation is very small, additional fluctuations in the calculated function values occur.

Gas molecular spectroscopy
Using this free-running dual-comb system, we demonstrated the spectral absorption characteristics of hydrogen cyanide (H 13 C 14 N, 25 Torr, 16.5 cm) gas cell at ambient temperature.The H 13 C 14 N molecule has two equally spaced absorption bands, the P and R branches, near 1.55 μm [43].
Physically, these absorption bands are generated by the multiplicative and combinatorial frequency absorption of low-energy electron leaps, and stretching vibrations of hydrogen-containing atomic groups (N-H bonds, C-H bonds) during the molecule's transition from the ground state to the excited state.A C-band wavelength division multiplexer is inserted in front of the gas cell to be measured and is used to concentrate the optical power in the range of the typical absorption peaks of H 13 C 14 N.  database [44] (dark blue curve in Fig. 7(a)), where the absorption lines are calculated using the standard Voigt profile since the line broadening is affected by both intermolecular collision and the Doppler effect.One should note the significant characteristic asymmetric spectral residue when comparing with the HITRAN database from Fig. 7(a), especially in the R-branch.Considering the symmetric DCS setup and the high peak power of the pulses, we believe that this deviation is dominated by the nonlinear processes, i.e., self-phase modulation (SPM) and cross-phase modulation (XPM), occurring in the common-fiber before passing through the gas sample [45,46].
Specifically, SPM alters each comb's initial spectral baseline, leading to intensity errors in the measurements.At the same time, XPM allows the pulses to attract with each other arising from the temporally intensity-dependent phase modulation.During asynchronous optical sampling, this leads to sampling errors in the timing and ultimately result in the unexpected asymmetric spectral residue.An effective solution to weaken the nonlinear penalty is to pre-chirp two lasers to reduce the peak power before the common fiber.Here we show four sets of spectral results without wavelength calibration, two of which show spectral linearity stretching and the other two correspond to spectral linearity compression, all of which set the PS-FBG's notch center as the "origin".Differences in elastic properties arise from inconsistencies in the pulse misalignment caused by relative timing jitter between dual combs during interference.Specifically, the greater the relative timing jitter, the more pronounced the wavelength deviation is.From the results, two conclusions need to be explained.Firstly, the Rbranch is further away from the "origin" than the P-branch, thus accumulating a more significant wavelength deviation in the same data set (conclusion of Fig. 6), whereas the P-branch presents more stable.Secondly, the wavelength dependence of the light source power determines the quality of the retrieved spectra.For our source, the average power in the R-branch is about 5 dB lower than that of the P-branch (Fig. 2(b) and Fig. 3(g)), leading to observed noise undulations between resonance peaks after coherent averaging.That would be effectively suppressed to achieve higherquality retrieved spectra if dual-comb spectra were pre-optimized.The R 4 and P 11 lines were selected as shown in Fig. 7 (d).The four curves show the spectral measurements performed at 100, 200, 300, and 428 times of coherent averaging, and all show agreement with the HITRAN reference, showing good potential for gas molecular spectroscopy applications.

Discussion
We demonstrated a free-running DCS in the Δf rep range of 100 to 300 Hz.However, the range is limited.For DCS applications, the effective spectral information within the optical detection bandwidth is mapped to be detected on the RF domain.Its spectral width can be expressed as ΔB For our system, with a λ 0 of 1550 nm and a f rep of 100 MHz, when Δν reaches 100 nm, the corresponding maximum detection frame rate can be calculated to be 400 Hz without being able to increase it since the aliasing of the RF spectrum.Using simultaneous in-phase and counterphase correction to extend the bandwidth to 3 Nyquist regions without comb lines aliasing [47], thus breaking the frame rate-bandwidth limit, is expected to increase the detection frame rate even further.In recent years, single-cavity dual-comb technology has flourished [48,49], perfectly solving the coherence problem of DCS on conventional dual-source platforms, as it produces sources with inherent excellence in inter-coherence without additional frequency locking.It has been realized in dual-wavelength lasers [50], polarization-multiplexing [51], directionalmultiplexing [52], and mechanical-sharing [53] scenarios, even in a dual-disk-laser system [54].
Our scheme can be ported to such single-cavity dual-comb systems for higher precision Δf replocking, where Δf ceo noise is suppressed by common mode instead of introducing a phase calibration process.
In summary, we present an F-P cavity-assisted scheme for enhancing relative timing jitter, i.e., Δf rep locking accuracy in a dual comb system.Thanks to the excellent wavelength thermal stability of the F-P cavity in the operating bandwidth, i.e. ±0.73 GHz in the range of -5 to 70 °C, the locking accuracy of Δf rep can be further improved.After correcting for Δf ceo penalty in the phase calibration process supported by a PS-FBG with a notch filter bandwidth of 118.3 MHz, the locking accuracy is even better than that of the electrically-referenced scheme.Thus, free-form DCS applications are permitted.Careful comparison of the frequency noise PSD curves of the signals from the phase calibration path under both free-running and electrically reference-locked shows that there is a similar frequency distribution in both states, thus guarantee the validity of performing phase calibration during free-running.Next, the validity of the wavelength calibration was also verified under free running, and the locking accuracy of the F-P cavity-assisted scheme was verified to be superior to the electrical reference scheme.Wavelength deviations below the theoretical optical resolution over the entire operating band allow precision spectroscopic measurements.At the same time, the lower bound of a 0.4-pm deviation provides the scheme with sub-picosecond pulse-topulse sensitivity in the measurement of relative timing jitter.As a practical spectroscopic measurement application, the spectral absorption characteristics of H 13 C 14 N gas molecules were measured at ambient temperature, which showed agreement with the HITRAN database in both the P and R branches.Our scheme promises for new ideas in sensitive measurements of relative timing jitter of dual dual-comb interferometry, such as absolute distance and precision spectral measurements.

Fig. 1 .
Fig. 1.Principle of the F-P cavity-assisted relative timing jitter locking scheme and comparison of offset repetition rate locking stabilities under three operating conditions.(a) The RF source provides a reference signal to lock f rep directly.(c) Schematic of equivalent locking using a phase-shifted fiber Bragg grating (PS-FBG) in free-running operation.(e) Schematic of F-P cavity assisted Δf rep locking in free-running operation.(b), (d), and (f) are the schematic representations of the frequency deviation corresponding to the three cases.
Fig. 2 The schematic diagram of the setup.(a) Two free-running combs pass through the sample, PS-FBG, and F-P cavity synchronously.Three optical streams are captured by a low-speed PD array and recorded through a 3-channel analog-to-digital converter (ADC) card.A commercial optical spectral analyzer (OSA) is employed to record transmission spectra of the F-P cavity simultaneously.All signals are processed on the computer.(b) The output spectra of the OFCs.OC, optical coupler.

F 2 3 .
Fig. 3. Block diagram of the on-line calibration process.(a), (b) are the raw data recorded by ADC card and (c) by OSA.The subscript of the signal indicates the fiber device that passed through.(d) and (e) plots the time-varying phase of (a), and (b) obtained by the Hilbert transform (HT) and their difference respectively.(f) The raw spectrum of (a) by performing the fast Fourier transform (FFT).
3(g), to a new one λ'.Meanwhile, the spectral profile in Fig. 3(g) was fitted and reshaped to obtain the normalized absorption spectrum v DCS F-P (λ') as an output version of the processing.From the vertical comparison of Figs.3(f) to (h), the raw, chaotic spectrum is calibrated for the Δf ceo penalty to obtain a secondary intermediate spectrum with a slight wavelength deviation, and the final normalized absorption spectrum is obtained by resampling with the reference spectrum to correct the wavelength deviation caused by the residual relative timing jitter.

Fig. 4 .
Fig. 4. The distribution of frequencies extracted from I PS-FBG (t).Calculation of ∆f rep evolution over time for (a) frame 2, (b) frame 7, and (c) frame 10 out of 16 frames under free-running operation.(d), (e) and (f) are the frequency distributions calculated by each IGM selected from the corresponding frame, respectively.

2 , 7 ,
and 10 from the 16-frame signal corresponding to Figs. 4(a) to (c), respectively.The scatter plots show the corresponding frames' computed ∆f rep evolution distribution over time.Under freerunning operation, the ∆f rep of the three frames are around 76.85, 165, and 246.9 Hz, respectively.As the ∆ rep increases, each frame contains more IGM cycles, i.e., corresponds to a denser scattering distribution with the same time window.By performing a short-time Fourier transform operation on the IGMs represented by each point in the scatter plot, the frequency-distributed evolution of I PS-FBG (t) can be calculated.Figures.4(d) to (f) are the results of calculating the frequency distributions of the 5th, 15th, and 22nd IGMs in the corresponding frames, respectively.

Fig. 5 .
Fig. 5. Comparison of (a) Frequency noise power spectral density (PSD) curves and (b) integral frequency jitter curves for two operating states calculated from the I PS-FBG (t).

Fig. 6 .
Fig. 6.Characterizations of wavelength calibration results for three operating states.(a) Statistical results of absolute wavelength deviations at the resonance peak of each F-P cavity for the three conditions.(b) Standard deviation after wavelength calibration in free-running operation.(c) and (d) are Allan deviations calculated after multiple coherent averaging of the resonance peaks of the F-P cavity away from and near the center wavelength, respectively.

Fig. 7 .
Fig. 7. Spectroscopic measurements of the P and R branches of the hydrogen cyanide (H 13 C 14 N) gas absorption cell in free-running operation.(a) The retrieved normalized spectral absorption features of H 13 C 14 N compared to the HITRAN database.(b) and (c) compare the difference before and after wavelength calibration in the R and P branches, respectively.(d) shows the spectral measurements of the R 4 and P 11 absorption lines of H 13 C 14 N molecules at different averaging times compared to the HITRAN database.

Figures 7 (
Figures 7(b) and (c) compare the retrieved transmission peaks from the R and P branches of

=
Δν/M = Δν•Δf rep /f rep , where ΔB and Δν are the electrical bandwidth and optical frequency bandwidth, respectively.The electrical bandwidth is also limited by Nyquist's law of sampling, i.e., ΔB ≤ f rep /2, and considering that there is a conversion relationship between the optical wavelength bandwidth and the optical frequency bandwidth Δν = c•Δλ/λ 2 0 , where c is the speed of light, λ 0, and Δλ are the central wavelengths of the light source and the effective bandwidth, respectively.From this, the constrained relationship of the detection frame rate Δf rep can be