1 Introduction
Precise timing synchronization is tremendously important in the large-scale scientific facilities[ Reference Xin, Şafak and Kärtner 1 ]. For instance, in high-power laser facilities (HPLFs) such as the SG-II project[ Reference Lin, Deng, Fan, Wang, Chen, Zhu, Qian, Shen, Xu, Zhu, Ma, Xie, Zheng, Zhang, Chen, Ling, Huang and Zhang 2 ] in China, the HERCULES laser[ Reference Bahk, Rousseau, Planchon, Chvykov, Kalintchenko, Maksimchuk, Mourou and Yanovsky 3 , Reference Yanovsky, Chvykov, Kalinchenko, Rousseau, Planchon, Matsuoka, Maksimchuk, Nees, Cheriaux, Mourou and Krushelnick 4 ], the SULF laser[ Reference Zhang, Wu, Hu, Yang, Gui, Ji, Liu, Wang, Liu, Lu, Xu, Leng, Li and Xu 5 ], the J-KAREN-P laser[ Reference Pirozhkov, Fukuda, Nishiuchi, Kiriyama, Sagisaka, Ogura, Mori, Kishimoto, Sakaki, Dover, Kondo, Nakanii, Huang, Kanasaki, Kondo and Kando 6 ] and the Texas Petawatt Laser[ Reference Tiwari, Gaul, Martinez, Dyer, Gordon, Spinks, Toncian, Bowers, Jiao, Kupfer, Lisi, McCary, Roycroft, Yandow, Glenn, Donovan, Ditmire and Hegelich 7 ], coherent beam combining (CBC) is regarded as a promising approach to increase the extreme intensity[ Reference Baker, Homoelle, Utterback and Jones 8 – Reference Bashinov, Gonoskov, Kim, Mourou and Sergeev 10 ]. The key to achieving high-efficiency CBC lies in ensuring that each individual laser is precisely spatially co-aligned and timing-synchronized with exceptional precision. To achieve even greater laser intensities, ultrahigh-power laser facilities, including Apollon[ Reference Zou, Le Blanc, Papadopoulos, Chériaux, Georges, Mennerat, Druon, Lecherbourg, Pellegrina, Ramirez, Giambruno, Fréneaux, Leconte, Badarau, Boudenne, Fournet, Valloton, Paillard, Veray, Pina, Monot, Chambaret, Martin, Mathieu, Audebert and Amiranoff 11 ], ELI-NP[ Reference Gales, Tanaka, Balabanski, Negoita, Stutman, Tesileanu, Ur, Ursescu, Andrei, Ataman, Cernaianu, D’Alessi, Dancus, Diaconescu, Djourelov, Filipescu, Ghenuche, Ghita, Matei, Seto, Zeng and Zamfir 12 ], EP-OPAL[ Reference Bromage, Bahk, Begishev, Dorrer, Guardalben, Hoffman, Oliver, Roides, Schiesser, Shoup, Spilatro, Webb, Weiner and Zuegel 13 ] and SEL[ Reference Wang, Liu, Lu, Chen, Long, Li, Chen, Chen, Bai, Li, Peng, Liu, Wu, Wang, Li, Xu, Liang, Leng and Li 14 ], capable of delivering outputs of 10 PW and above, are currently being developed, in which timing synchronization plays an even more important role. In addition to HPLFs, the next-generation X-ray free electron lasers (XFELs) have also garnered significant global attention[ Reference Abela, Aghababyan, Altarelli and Altucci 15 ]; these include the European XFEL in Germany, FERMI[ Reference Allaria, Appio, Badano, Barletta, Bassanese, Biedron, Borga, Busetto, Castronovo, Cinquegrana, Cleva, Cocco, Cornacchia, Craievich, Cudin, D’Auria, Dal Forno, Danailov, De Monte, De Ninno, Delgiusto, Demidovich, Di Mitri, Diviacco, Fabris, Fabris, Fawley, Ferianis, Ferrari, Ferry, Froehlich, Furlan, Gaio, Gelmetti, Giannessi, Giannini, Gobessi, Ivanov, Karantzoulis, Lonza, Lutman, Mahieu, Milloch, Milton, Musardo, Nikolov, Noe, Parmigiani, Penco, Petronio, Pivetta, Predonzani, Rossi, Rumiz, Salom, Scafuri, Serpico, Sigalotti, Spampinati, Spezzani, Svandrlik, Svetina, Tazzari, Trovo, Umer, Vascotto, Veronese, Visintini, Zaccaria, Zangrando and Zangrando 16 ] in Italy, SwissFEL in Switzerland and the Linac Coherent Light Source (LCLS)[ Reference Emma, Akre, Arthur, Bionta, Bostedt, Bozek, Brachmann, Bucksbaum, Coffee, Decker, Ding, Dowell, Edstrom, Fisher, Frisch, Gilevich, Hastings, Hays, Hering, Huang, Iverson, Loos, Messerschmidt, Miahnahri, Moeller, Nuhn, Pile, Ratner, Rzepiela, Schultz, Smith, Stefan, Tompkins, Turner, Welch, White, Wu, Yocky and Galayda 17 ] along with LCLS II in the US[ Reference Stohr 18 ]. To achieve the goal of generating attosecond (as, 10–18 s) X-ray pulses[ Reference Prat and Reiche 19 ] of unparalleled brightness and filming ultrafast physical and chemical processes with the spatial-temporal resolution of atoms, attosecond-precision timing synchronization is exigently required.
A typical synchronization system comprises a timing reference that produces exceptionally stable timing signals, a target signal that requires alignment, a timing detector that assesses the timing difference between the target and reference signals and a control unit that adjusts the timing of the target to match the reference. When the target device is located far from the reference, a timing link is essential to convey the timing signal from the reference to the target[ Reference Xin, Şafak, Peng, Kalaydzhyan, Wang, Mücke and Kärtner 20 ].
The passively mode-locked laser (MLL), which shows ultra-low timing jitter at high frequencies (e.g., <1 fs above 10 kHz)[ Reference Haus and Mecozzi 21 ], is considered as an ideal timing reference for synchronizing multiple microwave and optical sources. To break the limitation of amplitude modulation (AM)–phase modulation (PM) noise and provide precise timing resolution for the synchronization system, the balanced optical cross-correlator (BOC)[ Reference Jiang, Wong, Grein, Ippen and Haus 22 , Reference Schibli, Kim, Kuzucu, Gopinath, Tandon, Petrich, Kolodziejski, Fujimoto, Ippen and Kärtner 23 ] scheme was introduced. Over the past decades, the pulse-based BOC timing synchronization method has benefited several extreme laser facilities, such as SULF, ELI, the European XFEL, FERMI, SwissFEL and LCLS II.
However, as science facilities become steadily more complex, the limitations of the BOC timing synchronization system are gradually emerging. Currently, the primary challenge facing BOC-based timing synchronization systems lies in the inherent contradiction between optical power and synchronization precision. To ensure sub-fs synchronization precision while simultaneously suppressing additional timing jitter induced by fiber nonlinear effects[ Reference Xin, Şafak, Peng, Kalaydzhyan, Wang, Mücke and Kärtner 20 ], the optical power transmitted through the link must be constrained within the range of 10–20 mW. Taking transmission losses into account, this constraint limits the link length to under 5 km. In addition, since the output power of an MLL oscillator is usually 100–200 mW, the maximum number of synchronization terminals is also limited to approximately 10. Such limitations significantly impede the scalability of femtosecond-pulse-based timing synchronization technology toward greater distances and larger-scale implementations, thereby failing to meet the growing demands of large scientific facilities in terms of both physical dimensions and the number of terminal nodes.
To improve the scalability of these large facilities and ensure their continued operation with remarkable stability over the forthcoming decades, an ultra-precise timing and synchronization system that efficiently utilizes optical power and maintains a compact layout is urgently desired. In this paper, we have proposed a novel timing synchronization technology based on linear-optics timing detectors (LOTDs)[ Reference Wang, Ren, Şafak, Kärtner and Xin 24 ]. Each timing link can be stabilized to sub-fs precision with less than 1 mW optical power; therefore, the additional timing jitter induced by fiber nonlinear effects can be neglected. This work would potentially extend synchronization distances from several kilometers to tens of kilometers, and enable a single MLL oscillator to synchronize hundreds of terminals rather than just a few, paving the way for more sophisticated and powerful scientific facilities.
In the following, a detailed theoretical and experimental demonstration of the LOTD is first given in Sections 2 and 3. Using this timing detector, the residual dispersion of a 5.2 km fiber link is precisely compensated in Section 4. A detailed feedback model is then built to analyze the factors influencing the out-of-loop timing jitter in Section 5. Finally, a simulation of the complete timing link stabilization system with our new timing method is presented, which can satisfy the critical and demanding synchronization needs of next-generation photon science facilities and will be experimentally carried out in our future work.
2 Working principle of the linear-optics timing detector
Figure 1(a) illustrates the basic principle of the LOTD. The input from the laser source is injected into an acousto-optic modulator (AOM), which diffracts the pulse trains into different directions. A radio-frequency (RF) signal is then applied to the AOM to generate the original zeroth-order light and the modulated first-order light. The timing target under measurement is placed in the path of the first-order light, introducing a perturbative timing error Δt to the optical signal. A 3 dB fiber coupler integrates the two diffracted light paths, generating beat note signals that are subsequently detected by a photodiode (PD).
(a) Basic concept of the LOTD. FC, fiber collimator. (b) Simulated timing characterization curve, α = 1, τ = 170 fs. (c) Normalized timing sensitivity with different timing errors.
For the input optical pulse train, we can describe its electric field by the following expression:
$$\begin{align}{E}_0(t)=\sum \limits_{k=-\infty}^{+\infty }A\left(t- kT\right){e}^{-j{\omega}_0t},\end{align}$$
In Equation (1), A(t) characterizes the temporal envelope function, ω 0 defines the optical carrier frequency and T represents the temporal separation between pulses. The optical fields entering the collimators from the zeroth- and first-order diffraction light can be expressed as follows:
$$\begin{align}{E}_0(t)&=\sum \limits_{k=-\infty}^{+\infty }A\left(t- kT\right){e}^{-j{\omega}_0t}=\sum \limits_{n=-\infty}^{+\infty }{A}_n{e}^{- jn{\omega}_{\mathrm{rep}}t}{e}^{-j{\omega}_0t}, \end{align}$$
$$\begin{align}{E}_1(t)&=\sum \limits_{k=-\infty}^{+\infty }A\left(t-\Delta t- kT\right){e}^{-j\left({\omega}_0+{\omega}_{\mathrm{RF}}\right)\left(t-\Delta t\right)}\nonumber\\&=\sum \limits_{n=-\infty}^{+\infty }{B}_n{e}^{- jn{\omega}_{\mathrm{rep}}t}{e}^{-j\left({\omega}_0+{\omega}_{\mathrm{RF}}\right)\left(t-\Delta t\right)},\end{align}$$
where f rep = 1/T is the pulse train repetition frequency, f RF is the frequency of the RF driving signal, ω rep = 2πf rep, ω RF = 2πf RF and An , Bn are the Fourier series of the zeroth- and first-order pulse train profiles. After passing through the coupler, the optical power reaching the PD is given by the following:
$$\begin{align}P\left(t,\Delta t\right)\propto {\left|{E}_0-{jE}_1\right|}^2={\left|{E}_0\right|}^2+{\left|{E}_1\right|}^2+j\left({E}_0{E}_1^{\ast }-{E}_0^{\ast }{E}_1\right).\end{align}$$
To filter out the ω RF frequency component P(t, Δt), we use a band pass filter (BPF). Thus, the output signal is given by the following:
$$\begin{align}{V}_{\mathrm{BPF}}&\propto j\sum \limits_{n=-\infty}^{+\infty}\left[{A}_n{B}_n^{\ast }{e}^{j\left[{\omega}_{\mathrm{RF}}t-\left({\omega}_0+{\omega}_{\mathrm{RF}}\right)\Delta t\right]}\right.\nonumber\\&\nonumber\qquad\qquad\left.-{A}_n^{\ast }{B}_n{e}^{-j\left[{\omega}_{\mathrm{RF}}t-\left({\omega}_0+{\omega}_{\mathrm{RF}}\right)\Delta t\right]}\right]\\&=X\cos \left[\theta +{\omega}_{\mathrm{RF}}t-\left({\omega}_0+{\omega}_{\mathrm{RF}}\right)\Delta t\right],\end{align}$$
where
$X=\left|\sum \limits_{n=-\infty}^{+\infty }{A}_n{B}_n^{\ast}\right|$
, and
$j\sum \limits_{n=-\infty}^{+\infty }{A}_n{B}_n^{\ast}={Xe}^{j\theta}$
. It should be noted that the timing information Δt is reflected in both the amplitude and phase part of V
BPF. To eliminate the distraction of the phase component, a zero-bias Schottky diode (ZSD) is employed to extract the power fluctuation of V
BPF, thereby simplifying the resolution of Δt. Using the principle of Parseval’s theorem, the ZSD response can be formulated as follows:
$$\begin{align}{V}_{\mathrm{ZSD}}\left(\Delta t\right)\propto {\left|\sum \limits_{n=-\infty}^{+\infty }{A}_n{B}_n^{\ast}\right|}^{\alpha}\propto {\left|\underset{-\infty }{\overset{+\infty }{\int }}A(t)A\left(t-\Delta t\right)\mathrm{d}t\right|}^{\alpha },\end{align}$$
where α denotes a parameter linked to the nonlinear behavior of the ZSD. Hence, the timing characterization curve that indicates the relationship between timing error Δt and V ZSD can be derived. Suppose A(t) is a hyperbolic secant function[ Reference Haus 25 ]:
$$\begin{align}A(t)={A}_0\operatorname{sech}\left(t/\tau \right),\end{align}$$
where τ is a parameter that characterizes the pulse width. Substituting Equation (7) into Equation (6), the analytical solution of V ZSD can be further derived as follows:
$$\begin{align}V\left(\Delta t\right)\propto \left\{\begin{array}{l}{\left(\frac{4\Delta {te}^{\Delta t/\tau }}{e^{2\Delta t/\tau }-1}\right)}^{\alpha },\kern0.6em \Delta t\ne 0,\\ {}{\left(2\tau \right)}^{\alpha },\kern1.8em \Delta t=0.\end{array}\right.\end{align}$$
To eliminate all the trivial coefficients, the normalized voltage can be used:
$$\begin{align}{V}_{\mathrm{norm}}\left(\Delta t\right)=\frac{V\left(\Delta t\right)}{V(0)}={\left(\frac{2\frac{\Delta t}{\tau }{e}^{\Delta t/\tau }}{e^{2\Delta t/\tau }-1}\right)}^{\alpha }.\end{align}$$
A simulated normalized timing curve is given in Figure 1(b). On either side of this curve, V norm shows a consistent and monotonic change with Δt, which means it can be used to characterize Δt. In the magnified portion, the maximum slope of 0.0018 fs–1 is obtained, giving the maximum timing sensitivity. Based on Equation (9), the normalized timing sensitivity is given by the following:
$$\begin{align}\frac{\mathrm{d}{V}_{\mathrm{norm}}}{\mathrm{d}\Delta t}=\frac{1}{\tau }f(x),\end{align}$$
where
$$\begin{align}x=\frac{\Delta t}{\tau },\end{align}$$
$$\begin{align}f(x)=2\alpha {\left(\frac{2{xe}^x}{e^{2x}-1}\right)}^{\alpha -1}{e}^x\frac{e^{2x}\left(1-x\right)-\left(1+x\right)}{{\left({e}^{2x}-1\right)}^2}.\end{align}$$
Figure 1(c) shows the relationship between the normalized timing sensitivity and Δt. At the two peaks, the maximum sensitivity can be achieved when Δt = ± 1.6061τ.
The balanced detection is commonly employed to effectively lower the timing detection floor. A typical balanced operation consists of two identical timing detection circuits with a relative time delay, as shown in Figure 2(a). The zeroth-order light from the AOM is split into two paths, one of which incorporates a fixed delay T D. Simultaneously, the first-order diffraction light is divided via a 50:50 coupler: one branch combines with the undelayed zeroth-order light, while the other combines with the T D-delayed zeroth-order light through two additional 50:50 couplers. Both combined beams are processed through the PD, BPF and ZSD to generate independent timing signals. By differentially analyzing these two outputs, common-mode noises such as laser amplitude fluctuations and environmental disturbances are theoretically canceled, thereby enhancing the timing precision.
(a) Basic concept of the LOTD using balanced structure. (b) Balanced timing characterization curve with different delay T D. (c) Balanced normalized timing sensitivity with different delay T D.
The balanced output voltage is given by the following:
$$\begin{align}{V}_{\mathrm{B}}\left(\Delta t\right)={V}_{\mathrm{norm}}\left(\Delta t-{T}_{\mathrm{D}}\right)-{V}_{\mathrm{norm}}\left(\Delta t\right),\end{align}$$
and the corresponding timing curve becomes an ‘s’ shape, as depicted in Figure 2(b). At the zero-crossing Δt
0 of this curve, V
B(Δt
0) = 0, that is, V
norm(Δt
0 – T
D) = V
norm(Δt
0). Based on the symmetry of V
norm(Δt), Δt
0 = T
D/2. In order to attain maximum timing sensitivity at the crossing point, the slope of V
norm(Δt – T
D) and V
norm(Δt) at Δt = Δt
0 needs to be maximum and minimum, respectively. Therefore, the optimum delay time for balanced detection is
$T_\mathrm{D}=\pm3.2212\cdot \tau$
the detailed relation between dV
B/dΔt and T
D is given in Figure 2(c)). In Figure 2(b), several balanced timing curves with different T
D values are plotted; the one with T
D= 550 fs (closest to 3.2122τ) gives the maximum timing sensitivity around the zero-crossing point.
The LOTD has two key inherent advantages that make it well-suited for various applications. Firstly, if the zeroth- and first-order pulse trains have the same chirp, the dispersion phase coefficients can be canceled by the product
$A_n B_n^*$
in Equation (6). Therefore, transform-limited pulse widths are not essential for maintaining high timing detection sensitivity. This characteristic is particularly beneficial in situations where it is difficult or cost-prohibitive to compensate for high-order dispersion, yet identical chirp between two pulse trains is readily achievable, as exemplified in ultra-long fiber loop-based timing sensors. In addition, it offers enhanced timing resolution at low power levels when compared to nonlinear-optics-based detectors such as the BOC. Taking advantage of this merit, it is possible to perform timing link stabilization with ultra-low optical power, which will be discussed in Section 5.
3 Balanced linear-optics timing detector experimental characteristics
The detailed experimental setup of the balanced linear-optics timing detector (BLOTD) is presented in Figure 3(a). A stable MLL source (MENHIR-1550) generates pulse trains centered at 1555 nm with 216.667 MHz repetition frequency and 170 fs duration, and subsequently injected into an AOM. An 80 MHz sinusoidal waveform from a function generator, initially at 200 mV root mean square (RMS), is amplified to 36 dBm to drive the AOM. The power distribution between the zeroth- and first-order paths of diffraction light should be near-equal.
(a) The experimental setup for the BLOTD. λ/2, half-wave plate; λ/4, quarter-wave plate; PBS, polarization beam splitter; RFL, retroreflector; 50:50, 3 dB coupler; PC, personal computer; DA, differential amplifier. (b) Measured balanced timing characterization curve. (c) Measured results: TJSD and corresponding IDTJ.
A motorized delay line (MDL) is positioned in the first-order path to introduce a time delay Δt between the two beams. Balanced detection requires two independent timing measurement circuits for the zeroth- and first-order light. In the first circuit, the relative time delay between the two beams is Δt, while the second circuit introduces an additional delay T D (total: Δt+T D). To implement this, the zeroth-order light from the AOM is split into two paths via a polarization beam splitter (PBS). The light in one path is reflected by a silver mirror (SM) mounted on a manual stage, so as to pass twice through a quarter-wave plate. By finely adjusting the distance between the SM and the PBS using the manual stage, the relative delay T D for the optimal timing sensitivity near the zero-crossing point of the ‘s’-shaped curve can be obtained. With the help of three 50:50 fiber couplers, the zeroth- and first-order pulse trains beat on two avalanche photodiode detectors (APDs), and then two 80 MHz beat note signals can be obtained. The 100 MHz 3 dB bandwidth APD is customized with low noise equivalent power (2 pW/√Hz) and high responsivity (0.9 A/W). The beat note from each APD is initially extracted using a BPF and subsequently amplified by a low-noise amplifier (LNA). This process ensures that the power fluctuations, which contain the timing information of the beat note, are adequately amplified for further extraction by the ZSD. The balanced detection is finally achieved by differentially amplifying the two ZSDs’ outputs. The differential amplifier is driven by a matched low-noise power supply. The timing characterization curve, which illustrates the relationship between Δt and the ZSD output, is generated by recording the output voltage through a data acquisition card (DAQ) while concurrently moving the computer-controlled MDL. With 1 mW power for each pulse train at each APD, the timing sensitivity around the zero-crossing is about 43.537 mV/fs, as shown in Figure 3(b). Due to the nonuniformity of the two ZSDs, the s-like curve in Figure 3(b) is a little asymmetric.
The timing jitter spectral density (TJSD) of the BLOTD is assessed by routing the output from the differential amplifier to the baseband port of a signal source analyzer (SSA). The findings are illustrated with the dark blue curve in Figure 3(c). The lowest measured timing detection noise floor reaches approximately 1×10–10 fs²/Hz, which is nearly at the shot-noise limit (marked by the red dashed line). By opting for photodetectors with enhanced responsivity, it is feasible to further decrease the detection floor and approach the standard quantum limit for timing jitter in pulse trains, as predicted in Ref. [Reference Giovannetti, Lloyd and Maccone26] (represented by the blue dashed line), which is only 10 dB higher than the current detection threshold. Displayed at the bottom of Figure 3(c) is the integrated timing jitter (IDTJ). The total IDTJ within the range of 1 Hz to 1 MHz is found to be only 26.57 as.
4 Residual dispersion compensation with the linear-optics timing detector
In a typical pulse-based large-scale timing synchronization system[ Reference Xin, Şafak, Peng, Kalaydzhyan, Wang, Mücke and Kärtner 20 ], the timing link is stabilized by comparing the relative timing error between the round-trip link pulse and the new incoming pulse from the master laser. The dispersion of the timing link needs to be carefully compensated to guarantee high-precision timing detection as well as to restrict the Gordon–Haus timing jitter[ Reference Gordon and Haus 27 ] generated during link transmission. In this section, with the help of the LOTD, the residual dispersion of a timing link can be easily measured, and then accurate dispersion compensation can be performed afterwards.
In Equation (6), if the zeroth- and first-order pulse trains do not have the same chirp, the dispersion phase coefficients can no longer be canceled by the product
$A_n B_n^*$
, which will broaden the timing curve. Therefore, the width of the timing curve can be used to characterize the link residual dispersion. The experimental setup is implemented as shown in Figure 4(a). The output of the MLL is separated into two paths. An AOM is placed in one path to shift the optical pulses’ carrier frequency before they enter the timing link. The relative timing errors between the link pulse and the original laser pulse are measured by the LOTD electronics (PD, BPF and ZSD), and the timing curve can be obtained by moving the MDL. The timing link is composed of 4.5 km single-mode fiber (SMF) and 0.7 km dispersion compensating fiber (DCF). The DCF can compensate both the second- and third-order dispersion of the SMF simultaneously. To allow room for fine adjustments, the DCF overcompensates for the dispersion of about 10 m SMF.
Experimental setup for residual dispersion compensation and results. (a) Experimental setups. (b) Timing curve with different fiber patch cord lengths. (c) Measured timing curves with and without the timing link. (d) Timing curve measured with MDL movements in the forward and backward directions. (e) Normalized timing sensitivity with different fiber patch cord lengths.
The measured timing curves with and without the timing link are compared in Figure 4(c), in which an obvious broadening effect is observed due to the links’ residual dispersion. Different lengths of fiber patch cord are then introduced to finely tune the residual dispersion. With 4.8 m fiber patch cord, the two measured timing curves by moving the MDL in the forward and backward directions are compared in Figure 4(d). The two curves exhibit near-perfect overlap, indicating that the temperature drift of the timing link does not affect the measurement results. The timing curves with different lengths of fiber patch cord are given in Figure 4(b), and the maximum normalized timing sensitivity for each case is summarized in Figure 4(e). With 9.8 m patch cord, the narrowest timing curve is obtained with the highest timing sensitivity, which means the residual dispersion of the timing link is almost compensated.
5 Feedback model of the balanced linear-optics timing detector-based timing link stabilization
After the timing link’s dispersion is precisely compensated, it can be used to realize large-scale remote timing synchronization. Figure 5 gives a setup for timing link stabilization using the BLOTD. In the link stabilization block, the input optical pulse train is divided into two pulse trains in the reference path and in the link path. The link path consists of an MDL, an AOM, a fiber stretcher, a dispersion-compensated fiber link, a 45° Faraday rotator and a 40/60 transmission/reflection fiber mirror. The 45° Faraday rotator can make sure that the polarization directions of forward and backward transmitted light in the link are orthogonal, so as to minimize the timing jitter induced by polarization mode dispersion. The relative timing error J IL between the pulses in the reference and link path is measured by the in-loop BLOTD and fed back to control the fiber stretcher through a proportional–integral (PI) controller, and thus can stabilize the length of the timing link. In the out-of-loop measurement block, the residual timing jitter J O between the link output pulse and the original pulse from the MLL is measured by the out-of-loop BLOTD to evaluate the performance of the link stabilization.
Detailed scheme of the link network stabilization. HVA, high-voltage amplifier; T40/R60, transmission 40%/reflection 60% fiber mirror; FR, Faraday rotator.
Since the BLOTD is a linear-optics device, it is possible to operate this link stabilization with very low power levels. To investigate the minimum required optical power for sub-fs precision link stabilization, based on the setup in Figure 5, a feedback model using experimental parameters is constructed. The flow diagram of the feedback model is shown in Figure 6.
Feedback flow diagrams of (a) in-loop timing link stabilization and (b) out-of-loop jitter measurement. Here, ω is the complex frequency and s = jω; J I, inherent jitter of the mode-locked laser; J E, environmental jitter imposed on the link for single-trip link transmission; J IL, detected timing jitter by the in-loop BLOTD; H BLD, transfer function of the in-loop BLOTD; E N, electronic noise in the in-loop BLOTD electronics; H DA, transfer function of the differential amplifier in the BLOTD; H PI, transfer function of the PI controller; E PI, electronic noise of the PI controller; H FS, transfer function of the fiber stretcher; J C, equivalent timing delay generated by the control loop for compensation; J S, shot-noise-equivalent timing jitter; T L, single-trip link transmission time; J O, relative timing jitter between the link output pulses and the original pulses from the mode-locked laser.
The in-loop BLOTD detected timing jitter J IL and the shot-noise-equivalent timing jitter J S are converted to voltage signals with the BLOTD transfer function H BLD, amplified by the differential amplifier H DA together with the electronic noise E N(V/√Hz), and then sent to the PI controller H PI. The output of the PI controller with additive noise E PI then undergoes amplification and subsequent conversion to timing delays J C through the fiber stretcher transfer function H FS, where J C acts as a compensating jitter to maintain the lock. Due to the round-trip link propagation, J C is added to the link delay twice, and the round-trip link propagation time 2T L is accounted for in both J C and the MLL’s inherent jitter J I; whereas single-trip link propagation time T L is utilized for environmental jitter J E, as the latter integrates all the single-trip environmental jitter in its definition. For the out-of-loop measurement, J C and J I experience a single-pass link propagation delay T L. All the noise sources J I, J E, J S, E N, E PI are characterized with the method in Appendix A. The transfer functions can be calculated by the equations given in Appendix B.
Based on Figure 6, we have
$$\begin{align}\kern-10pt {J}_{\mathrm{IL}}&={J}_{\mathrm{I}}\left[\exp \left(-2{sT}_{\mathrm{L}}\right)-1\right]+{J}_{\mathrm{E}}\left[\exp \left(-{sT}_{\mathrm{L}}\right)+1\right]\nonumber\\& \quad +{J}_{\mathrm{C}}\left[\exp \left(-2{sT}_{\mathrm{L}}\right)+1\right], \end{align}$$
$$\begin{align}{J}_{\mathrm{C}}&=\left[-\left({J}_{\mathrm{IL}}{H}_{\mathrm{BLD}}+{J}_{\mathrm{S}}{H}_{\mathrm{BLD}}+{E}_{\mathrm{N}}\right){H}_{\mathrm{DA}}{H}_{\mathrm{PI}}+{E}_{\mathrm{PI}}\right]{H}_{\mathrm{FS}},\end{align}$$
$$\begin{align}{J}_{\mathrm{O}}&={J}_{\mathrm{I}}\left[\exp \left(-{sT}_{\mathrm{L}}\right)-1\right]+{J}_{\mathrm{E}}+{J}_{\mathrm{C}}\exp \left(-{sT}_{\mathrm{L}}\right).\end{align}$$
Substituting Equations (14) and (15) into Equation (16), the explicit expression of J O can be derived as follows:
$$\begin{align}{J}_{\mathrm{O}}={C}_{\mathrm{I}}{J}_{\mathrm{I}}+{C}_{\mathrm{E}}{J}_{\mathrm{E}}+{C}_{\mathrm{S}}{J}_{\mathrm{S}}+{C}_{\mathrm{N}}{E}_{\mathrm{N}}+{C}_{\mathrm{PI}}{E}_{\mathrm{PI}},\end{align}$$
where
$$\begin{align}{C}_{\mathrm{I}}&=\left[\exp \left(-{sT}_{\mathrm{L}}\right)-1\right]\frac{1+H\left[1-\exp \left(-{sT}_{\mathrm{L}}\right)\right]}{1+H\left[1+\exp \left(-2{sT}_{\mathrm{L}}\right)\right]}, \end{align}$$
$$\begin{align}{C}_{\mathrm{E}}&=\frac{1+H\left[1-\exp \left(-{sT}_{\mathrm{L}}\right)\right]}{1+H\left[1+\exp \left(-2{sT}_{\mathrm{L}}\right)\right]},\end{align}$$
$$\begin{align}{C}_{\mathrm{S}}&=-\frac{H\exp \left(-{sT}_{\mathrm{L}}\right)}{1+H\left[1+\exp \left(-2{sT}_{\mathrm{L}}\right)\right]},\end{align}$$
$$\begin{align}{C}_{\mathrm{N}}&=-\frac{H_{\mathrm{DA}}{H}_{\mathrm{PI}}{H}_{\mathrm{FS}}\exp \left(-{sT}_{\mathrm{L}}\right)}{1+H\left[1+\exp \left(-2{sT}_{\mathrm{L}}\right)\right]},\end{align}$$
$$\begin{align}{C}_{\mathrm{PI}}&=\frac{H_{\mathrm{FS}}\exp \left(-{sT}_{\mathrm{L}}\right)}{1+H\left[1+\exp \left(-2{sT}_{\mathrm{L}}\right)\right]},\end{align}$$
$$\begin{align}H&={H}_{\mathrm{BLD}}{H}_{\mathrm{DA}}{H}_{\mathrm{PI}}{H}_{\mathrm{FS}}.\end{align}$$
All the noise contributions CiJi
(i = E, S, I) and
$C_jE_j (j=\mathrm{N},\ \mathrm{PI})$
on the right-hand side of Equation (17) are mutually uncorrelated, and the mean values of their cross-products are cancelled. Consequently, the average out-of-loop jitter spectral density (s²/Hz) can be represented as follows:
$$\begin{align}{J_{\mathrm{O}}}_2={\left|{C}_{\mathrm{I}}\right|}^2\overline{{J_{\mathrm{I}}^2}}+{\left|{C}_{\mathrm{E}}\right|}^2\overline{{J_{\mathrm{E}}^2}}+{\left|{C}_{\mathrm{S}}\right|}^2\overline{{J_{\mathrm{S}}^2}}+{\left|{C}_{\mathrm{N}}\right|}^2\overline{E_{\mathrm{N}}^2}+{\left|{C}_{\mathrm{PI}}\right|}^2\overline{E_{\mathrm{PI}}^2}\\[-18pt]\nonumber\end{align}$$
As Equation (24) indicates, J O has five main contributions: the environmental noise imposed on the link; the electronic noise of the BLOTD; the noise introduced by the PI controller; the shot noise; and the inherent jitter of the MLL, with coefficients C E, C N, C PI, C S and C I, respectively. Figures 7(a)–7(e) present the calculated coefficients |Ci| (i = E, N, PI, S, I) with a 10.4 km link length for several distinct PI controller gain settings. To effectively reduce environmental noise and the laser’s inherent noise at frequencies lower than 5 kHz, it is essential to maintain high gain settings, as demonstrated in Figures 7(a) and 7(e). Nevertheless, elevating the gain leads to increased BLOTD electronic noise and shot noise, as illustrated in Figures 7(b) and 7(d). The large gain peaks can be observed near frequencies n/4T L (n = 1, 3, 5, …) as well as the fiber stretcher’s resonant frequency (approximately 18 kHz) in Figures 7(a)–7(e).
Simulation results for a 10.4 km timing link stabilization system. (a) The coefficient for the environmental noise imposed on the link. (b) The coefficient for the electronic noise. (c) The coefficient for the noise from the PI controller. (d) The coefficient for the shot noise. (e) The coefficient for the laser’s inherent jitter. (f) Out-of-loop integrated timing jitter with different PI gains and corner frequencies (1 mW optical power at each PD of the in-loop BLOTD).
Using the calculated coefficients above, each noise contribution in Equation (24), as well as J O 2, can be calculated. Typical experimental parameters are used in the simulation as given in Table 1 in Appendix C. Firstly, the optical power at each PD of the in-loop BLOTD is set as 1 mW. Through continuous optimization of the PI gain and corner frequency, we can achieve optimal link stabilization with minimal IDTJ of J O 2, as shown in Figure 7(f). With the optimized PI controller settings, the TJSD of J O 2 and each noise contribution are given in Figure 8(a). The total IDTJ from 1 Hz to 1 MHz is only 0.36 fs (black curve). The environmental jitter of the link (pink curve) is the predominant noise factor in the low-frequency range. The inherent noise contribution of the MLL becomes more noticeable for offset frequencies above 100 Hz.
Out-of-loop jitter simulation results of individual jitter contributions with the optical power of (a) 1 mW and (b) 1 μW at each PD of the in-loop BLOTD. The timing link length is 10.4 km. The legend on the right shows the color code of the jitter spectral densities. Integrated jitter in this graph is shown on a logarithmic scale.
To investigate the link stabilization performance at ultra-low power levels, the optical power at each PD of the in-loop BLOTD is reduced to 1 μW. After re-optimizing the PI controller parameters, the TJSD and IDTJ results are as shown in Figure 8(b), and the timing link can still maintain sub-fs precision with 0.48 fs IDTJ at [1 Hz, 1 MHz]. To compensate the low timing sensitivity of the BLOTD at the 1 μW power level, a larger feedback gain is necessary to help suppress the environmental noise of the timing link. Compared with Figures 8(a) and 8(b), it can be seen that as the link operational power decreases, both the shot noise and the electronic noise from the PI controller are elevated, while the electronic noise of the BLOTD increases the most dramatically and becomes the major contributor to the increase of IDTJ. Nevertheless, the total out-of-loop timing jitter is still dominated by J E and J I.
Furthermore, timing link stabilization with different fiber lengths is also simulated. Keeping the in-loop PD power as 1 μW, an IDTJ of 0.79 fs can still be achieved for 20 km link stabilization (Figure 9(a)). While if the link length is reduced to 5.2 km, an unprecedented low optical power of 240 nW at each in-loop PD is sufficient to provide sub-fs link stabilization precision (with 0.93 fs integrated jitter, as shown in Figure 9(b)). Leveraging the advantages of the BLOTD, we theoretically demonstrated the capability to extend the link length to nearly five times that of previous works using the BOC, and the required operational power is also reduced by more than three orders of magnitude. This breakthrough overcomes the current power-imposed limitations on transmission distance and the number of synchronization terminals, holding significant implications for the expansion and construction of various large-scale laser facilities.
Out-of-loop jitter simulation results of individual jitter contributions at the in-loop PD power of (a) 1 μW with transmission length of 20 km and (b) 240 nW with transmission length of 5.2 km.
6 Conclusions
In summary, we have experimentally demonstrated an attosecond-precision balanced linear-optics-based timing detector. With 1 mW power per pulse train per photodetector, a shot-noise-limited timing detection floor of 1×10–10 fs2/Hz is achieved. The balanced configuration effectively mitigates environmental and laser amplitude noise. The total integrated jitter from 1 Hz to 1 MHz is merely 26.57 as. As one of the applications, we have conducted the residual dispersion evaluation with this powerful tool. Finally, we have presented a comprehensive jitter analysis for the timing link stabilization. At 1 mW and 1 μW optical power at each PD of the in-loop BLOTD, the 10.4 km link is stabilized with total integrated jitter of 0.359 fs and 0.48 fs, respectively. The system can still maintain sub-fs precision with lower optical power or longer link length. For instance, with 5.2 km link length, an integrated jitter of 0.93 fs can be achieved with only 240 nW power, while with 1 μW power, after 20 km link transmission, the residual timing jitter is only 0.79 fs.
We believe that the timing link stabilization at ultra-low optical power levels will not only greatly expand the distance and number of terminals of the current synchronization systems, but also be beneficial in many other cutting-edge fields such as ultra-long-range laser radar and gravitational wave detection.
Appendix
A Noise source characterization
The inherent jitter J I is obtained from our previous BOC characterization results of a commercial MLL[ Reference Şafak, Xin, Callahan, Peng and Kärtner 28 ]. The shot-noise-equivalent timing jitter J S can be calculated by the Equation (9) of Ref. [Reference Wang, Ren, Şafak, Kärtner and Xin24]. The electronic noise E N of the BLOTD (generated from the photodetector, the BPF, the RF amplifiers, the ZSD and the differential amplifier) is measured when there is no input light in the system and then converted to timing jitter (in units of s/√Hz) by multiplying with the timing sensitivity. The environmental noise jitter J E is divided into two parts: at low frequencies (<100 Hz), the noise is modeled by a combination of 1/f and 1/f² noise, reflecting the typical properties of environmental noise; for frequencies of more than 100 Hz, the measured data from Ref. [Reference Şafak, Xin, Zhang, Chia, Mücke and Kärtner29] are utilized. The electronic noise of the PI controller E PI is evaluated utilizing a noise model predicated with the characteristics of the amplifier circuits presented in Figure 10(a), and can be estimated by Equations (25) and (26) [ Reference Şafak, Xin, Zhang, Chia, Mücke and Kärtner 29 ]:
$$\begin{align}{E}_{\mathrm{PI}}^2&=\frac{4{Z}_{\mathrm{f}}^2{k}_{\mathrm{B}}T}{R_{\mathrm{in}}}+4{k}_{\mathrm{B}}{TR}_{\mathrm{f}}+{\left(1+\frac{Z_{\mathrm{f}}}{R_{\mathrm{in}}}\right)}^2{V}_{\mathrm{n}}^2, \end{align}$$
$$\begin{align}{Z}_{\mathrm{f}}&=\frac{1/{R}_{\mathrm{f}}- j\omega C}{{\left(1/{R}_{\mathrm{f}}\right)}^2+{\omega}^2{C}^2}.\end{align}$$
(a) Noise model of the PI controller. V in and V out, input and output voltage; V n, input voltage noise of the PI controller; R in and i th,in, input impedance and its thermal noise current; R f and i th,f, feedback gain resistance and its thermal noise current; C, integrating capacitance. (b) Equivalent circuit of the fiber stretcher.
B Transfer function calculation
The transfer functions H BLD, H DA, H PI and H FS in Figure 6 can be obtained by the following equations:
$$\begin{align}{H}_{\mathrm{BLD}}&={k}_{\mathrm{BLD}}, \end{align}$$
$$\begin{align}{H}_{\mathrm{PI}}&={k}_{\mathrm{PI}}\left(\frac{s+2\pi {f}_{\mathrm{PI}}}{s}\right),\end{align}$$
$$\begin{align}{H}_{\mathrm{DA}}&={\left(1+\frac{s}{2\pi {f}_{\mathrm{BW}}}\right)}^{-1},\end{align}$$
$$\begin{align}{H}_{\mathrm{FS}}&={H}_{\mathrm{PZT}}{H}_{\mathrm{F}}.\end{align}$$
where k BLD is the sensitivity of the BLOTD (in units of mV/fs), k PI is the PI gain, f PI is the PI corner frequency and f BW is the 3 dB bandwidth of the differential amplifier. The transfer function H FS is the product of two parts: the lead zirconate titanate (PZT) response function H PZT fitted from manufacturer data and the filter function H F by the equivalent driven circuits (shown in Figure 10(b))[ 30 ].
C Typical parameters of the link stabilization used for simulation
Parameters used for simulation of the timing link stabilization.
Acknowledgements
This work was supported by the National Key R&D Program of China (Grant No. 2021YFC2201902) and the National Natural Science Foundation of China (Grant No. 61975149).
Open access
