Photonics Research, Volume. 12, Issue 9, 2033(2024)

Dynamic counterpropagating all-normal dispersion (DCANDi) fiber laser

Neeraj Prakash1, Jonathan Musgrave1, Bowen Li1,2,3, and Shu-Wei Huang1、*
Author Affiliations
  • 1Department of Electrical, Computer, and Energy Engineering, University of Colorado Boulder, Boulder, Colorado 80309, USA
  • 2Key Laboratory of Optical Fiber Sensing and Communications (Ministry of Education), University of Electronic Science and Technology of China, Chengdu 611731, China
  • 3e-mail: bowen.li@uestc.edu.cn
  • show less

    The fiber single-cavity dual-comb laser (SCDCL) is an emerging light-source architecture that opens up the possibility for low-complexity dual-comb pump-probe measurements. However, the fundamental trade-off between measurement speed and time resolution remains a hurdle for the widespread use of fiber SCDCLs in dual-comb pump-probe measurements. In this paper, we break this fundamental trade-off by devising an all-optical dynamic repetition rate difference (Δfrep) modulation technique. We demonstrate the dynamic Δfrep modulation in a modified version of the recently developed counterpropagating all-normal dispersion (CANDi) fiber laser. We verify that our all-optical dynamic Δfrep modulation technique does not introduce excessive relative timing jitter. In addition, the dynamic modulation mechanism is studied and validated both theoretically and experimentally. As a proof-of-principle experiment, we apply this so-called dynamic CANDi (DCANDi) fiber laser to measure the relaxation time of a semiconductor saturable absorber mirror, achieving a measurement speed and duty cycle enhancement factor of 143. DCANDi fiber laser is a promising light source for low-complexity, high-speed, high-sensitivity ultrafast dual-comb pump-probe measurements.

    1. INTRODUCTION

    Ultrafast pump-probe measurement has played a pivotal role in the study of femtosecond light–matter interaction and carrier dynamics as well as the characterization of chemical, solid-state, and biological materials [1]. Among all the variants, the dual-comb pump-probe technique offers a unique nonmechanical delay scanning mechanism that solves the measurement speed bottleneck and enables high-speed applications including terahertz time-domain spectroscopy [24], ultrafast photoacoustic characterization [5,6], high-sensitivity photothermal spectroscopy [7], multidimensional coherent spectroscopy [8,9], and coherent Raman spectroscopy [1013]. In the dual-comb pump-probe technique, the delay scanning is controlled by setting the repetition rate difference (Δfrep) between the two combs such that the two pulse trains will temporally walk off each other by steps of Δfrep/frep2, where frep is the average repetition rate. The total delay scan is 1/frep, and the measurement speed is limited by 1/Δfrep.

    The single-cavity dual-comb laser (SCDCL) is an emerging dual-comb architecture that has been studied over the past few years in various platforms including on-chip microresonators [14], Kerr-lens mode-locked lasers [15], semiconductor disk lasers [16], and fiber lasers [1722]. Therein, fiber SCDCLs are a particularly attractive architecture owing to their compactness and environmental robustness. Fiber SCDCL opens up the possibility for low-complexity dual-comb systems and can potentially lead to a paradigm shift in dual-comb light-source development [23,24]. In particular, we recently demonstrated the counterpropagating all-normal dispersion (CANDi) fiber laser that pushed the energy limit in fiber SCDCLs by 2 orders of magnitude to 8-nJ pulse energy in each comb [25,26]. 190-fs compressed pulse duration and 39-fs integrated relative timing jitter (1 kHz, 20 MHz) have also been demonstrated and characterized. All these features render CANDi fiber laser an ideal source candidate for ultrafast dual-comb pump-probe measurements.

    Another obstacle to the widespread adoption of the dual-comb pump-probe technique is the fundamental trade-off between the time resolution and measurement speed, due to their opposite scaling laws with respect to the Δfrep [27]. In addition, for applications such as coherent Raman spectroscopy, where the vibrational coherence lifetime of most molecules is only a few picoseconds, the dual-comb pump-probe technique’s typical >1-ns total delay scan is much longer than necessary, resulting in a <1% measurement duty cycle and a significant waste of dual-comb power that effectively lowers the measurement sensitivity [12]. To break the fundamental trade-off and enhance the time and energy efficiency, a few elegant ways have been devised recently by nonmechanical cavity length control of one of the dual combs [11,13,28]. Impressive results such as orders-of-magnitude measurement speed enhancement without sacrificing resolution and sensitivity and near 100% duty cycle coherent Raman spectroscopy have been demonstrated [11,13]. However, none of these methods are compatible with the SCDCL architecture, let alone fiber SCDCL.

    In this paper, we devise an all-optical dynamic Δfrep modulation technique that is applicable to SCDCL to overcome the fundamental trade-off between time resolution and measurement speed. We demonstrate a dynamic CANDi (DCANDi) fiber laser architecture and illustrate a dynamic Δfrep switching between >+10  Hz and <10  Hz at kilohertz-level modulation frequency. We characterize the relative timing jitter of DCANDi fiber laser to be the same as that of CANDi fiber laser, confirming that our all-optical dynamic Δfrep modulation technique does not introduce excessive noise. Furthermore, the dynamic modulation mechanism is studied and validated both theoretically and experimentally. As a proof-of-principle experiment, we apply DCANDi fiber laser to measure the relaxation time of a semiconductor saturable absorber mirror (SESAM), achieving a measurement speed and duty cycle enhancement of 143 times compared to traditional dual-comb pump-probe measurement. DCANDi fiber laser is a promising light source for low-complexity, high-speed, high-sensitivity ultrafast dual-comb pump-probe measurements.

    2. DCANDI: ARCHITECTURE AND PRINCIPLE

    The simplified DCANDi laser configuration is shown in Fig. 1(a). The cavity of DCANDi laser is constructed similar to that in Ref. [25]. However, in order to bring Δfrep closer to zero, the lengths of the passive fiber on the two sides of the ytterbium-doped fiber (YDF) are deliberately kept the same, which makes the cavity symmetric for both directions. The fiber section consists of 2  m long 6-μm-core YDF (Thorlabs YB1200-6/125DC) and 1  m of HI1060 passive fiber on each side of the YDF. The length of gain fiber is chosen to be long enough to achieve the gain asymmetry for dynamic modulation (mechanism discussed in Section 3 and Appendix A) but short enough to avoid reabsorption loss for efficient lasing. The mode locking is achieved using nonlinear polarization rotation. Moreover, the YDF is pumped bidirectionally by two 980-nm laser diodes with same power level to make the gain distribution inside the YDF symmetric for both directions. In this way, Δfrep can be easily tuned to close to zero by fine-tuning the intracavity wave plates.

    Principle of operation. (a) Simplified experimental setup of the DCANDi laser; (b) (top) schematic diagram of the Δfrep modulation; (middle) temporal walk-off of the two combs under modulation; (bottom) enhancement of interferogram repetition rate under modulation.

    Figure 1.Principle of operation. (a) Simplified experimental setup of the DCANDi laser; (b) (top) schematic diagram of the Δfrep modulation; (middle) temporal walk-off of the two combs under modulation; (bottom) enhancement of interferogram repetition rate under modulation.

    Under this condition, by modulating the current of either of the pumps (pump2 is modulated here) with a square waveform, the sign of Δfrep can be dynamically flipped, as schematically shown in Fig. 1(b). In order to make sure the absolute value of Δfrep is the same before and after flipping, a slow phase-locked loop (PLL) with 100 Hz control bandwidth is applied to pump1 to lock the average Δfrep to zero. More details on the PLL can be found in Appendix E. Since the control bandwidth is much smaller than the modulation frequency (kilohertz), the PLL does not counteract the modulation. When the average value of Δfrep is successfully locked to zero, the pulses from the two combs will temporally walk off each other alternatively as the sign of Δfrep flips [Fig. 1(b)] without needing to walk off for a whole round trip. As a result, the temporal interferogram frequency (i.e., the DC detection frequency) will be twice the modulation frequency, instead of Δfrep, as in a conventional DC system. In this way, the DCANDi system can dramatically increase its DC detection speed without sacrificing its temporal sampling resolution and pulse energy, thus providing an ideal tool for high-speed nonlinear DC applications.

    3. MECHANISM OF DYNAMIC MODULATION IN DCANDI

    The mechanism of pump power change-induced dynamic Δfrep modulation is mainly attributed to the nonidentical center frequency change coupled with intracavity group delay dispersion (β2), which is identified here through numerical simulation. The dependence of frep on pump power has been extensively studied in Refs. [25,26]. The change in frep due to change in pump power can be due to several factors, including spectral shift (ωΔ) coupled with β2, change in root mean square spectral width coupled with third-order dispersion, gain, and self-steepening terms. In Ref. [26], we identified the ωΔ term coupled with β2 to be the main coupling mechanism, resulting in pump-power-induced frep and Δfrep change in the CANDi laser. Considering only this term, the total pump-power-(P)-induced changes in frep and Δfrep are governed by the following equations, as explained in Ref. [26]: dfrepdP=frep2×(β2dωΔdP),dΔfrepdP=frep2×(β2dωΔCWdPβ2dωΔCCWdP)=frep2×(β2dδωΔdP).

    In the above equation, ωΔCW and ωΔCCW are the spectral shift in clockwise (CW) and counterclockwise (CCW) directions, while δωΔ is the relative spectral shift between the two directions. Since both combs in the DCANDi laser share the same pump, the mechanism for pump-power-induced change on Δfrep can be attributed to the asymmetric response of the two directions to pump-power change. When only one of the pumps is modulated, the gain distribution along the length of the YDF becomes asymmetric for the two directions. This results in unequal gain center shift, resulting in different spectral shift for the two directions for a given pump-power change. This causes unequal change in frep in two directions and hence a change in Δfrep.

    To confirm this hypothesis, we developed a simulation routine to model DCANDi based on solving the generalized complex Ginzburg–Landau equation coupled with rate equations. Since the gain inside the YDF dynamically changes as the pump power is modulated, a key requirement of the simulation is to calculate the true gain by solving the rate equations and generalized complex Ginzburg–Landau equation simultaneously. Note that a steady state condition of upper state population and a perfect symmetric cavity is assumed in the simulation. The results of simulation are shown in the Figs. 2(a) and 2(b). Figure 2(a) illustrates the δωΔ as a function of round trip for a particular mode-locked state as the pump power is modulated. Note that the simulation is displayed for 600 round trips, with the pump power modulated every 100 round trips. Figure 2(b) shows the comparison of the simulated Δfrep deduced from the difference in pulse timing against Δfrep calculated from the values in Fig. 2(a) using Eq. (2). The simulated Δfrep matches well with the Δfrep calculated using β2 coupled with the δωΔ. This confirms that unequal gain center shift, resulting in different ωΔ for the two directions for a given pump-power change is the main reason behind dynamic modulation of Δfrep. More insights on the importance of longitudinal gain profile in causing pump-power-dependent δωΔ can be found in the Appendix section. Later, in Section 4, we experimentally confirm the findings of the simulation to cement our understanding on the mechanism. More details on the simulation and simulation results are included in Appendix A.

    DCANDi simulation routine results. (a) Modulation of the δωΔ in DCANDi due to modulation of pump2 power; (b) comparison of the expected evolution of Δfrep calculated from cavity β2 coupled with ωΔ (red circle) with the Δfrep retrieved from pulse timing (black).

    Figure 2.DCANDi simulation routine results. (a) Modulation of the δωΔ in DCANDi due to modulation of pump2 power; (b) comparison of the expected evolution of Δfrep calculated from cavity β2 coupled with ωΔ (red circle) with the Δfrep retrieved from pulse timing (black).

    These simulation results help us understand the mechanism of Δfrep modulation in DCANDi laser and pave the way for custom designing DCANDis in the future. Beyond unraveling the mechanism of dynamic modulation in DCANDi, the simulation routine developed in this work contributes significantly to the broader field of SCDCL. The simulation can model complex interplay of parameters governing bidirectional SCDCL performance and help achieve precision controlled bidirectional SCDCL. The theoretical foundation discussed in this work can not only help optimize existing SCDCL designs but also lead to novel designs and applications.

    4. RESULTS AND DISCUSSION

    Figure 3 shows the basic performance of the DCANDi laser. Under static state (no current modulation), the optical spectra of the two combs are shown in Fig. 3(a). The two spectra show high spectral overlapping, and both represent a “batman” spectrum typical for ANDi lasers. The DCANDi laser has a repetition rate of 48 MHz and pulse energies of about 2 nJ for both directions. Higher pulse energies can be obtained by replacing the current 6-μm-core fibers with 10-μm-core large-mode area fibers [26]. A single-pulse compressor made of transmission grating pairs is used to compress pulses from both directions to around 180 fs, with total transmission of over 85% (see Appendix B for more details). Therefore, the peak power for compressed pulses is about 11 kW, which directly enables many nonlinear DC applications. To intuitively show the operation of dynamic Δfrep modulation, we first modulate the current of pump2 with about 1 Hz frequency and simultaneously monitor the repetition rate of both combs using two frequency counters referenced to the same clock. The corresponding current modulation depth is around 3%, which does not destroy the mode-locking state, thanks to the wide mode-locking range of ANDi lasers. The resulting Δfrep is shown in Fig. 3(b), which periodically and deterministically switches from +11 to 11  Hz. The local fluctuation results from the limited resolution of frequency counter under 0.03 s arm time. To demonstrate the real potential of dynamic operation of the DCANDi laser, the modulation frequency is then increased to 1 kHz and the PLL is turned on. The two combs are then combined through a 50/50 coupler and launched together into an amplified photodetector. Since the pulses from the two combs do not walk off monotonically when the average Δfrep is locked to zero, a variable optical delay line (ODL) is deployed on the path of one comb to make sure the pulses can temporally overlap each other at the PD. The output of the PD is low-pass filtered with 28 MHz; the corresponding temporal waveform is shown in Fig. 3(c). As expected, the temporal interferogram between the two combs shows a repeating frequency of 1 kHz, while two interferograms appear in each period, which results from forward and backward (shadowed with red) scanning. The time separations between a forward-scanning interferogram and two neighboring backward-scanning interferograms can be tuned to be identical by adjusting the relative delay between the two combs, leading to an effective DC detection speed of 2 kHz. The zoomed-in temporal waveforms of a forward-scanning and a backward-scanning interferograms are shown in Figs. 3(d) and 3(f), respectively, and their corresponding spectra are obtained by performing fast Fourier transform (FFT), as shown in Figs. 3(e) and 3(g). As observed, the two spectra are located at different center RF frequencies. This is due to the change of carrier-envelope-offset frequency (Δfceo) in response to the modulation of pump power. Nevertheless, both Figs. 3(e) and 3(g) show an RF spectrum that matches well with the overlapping optical spectrum in Fig. 3(a), which indicates the feasibility of fast and accurate DC metrology using the DCANDi laser.

    Basic performance of DCANDi laser. (a) Optical spectra of the two combs without modulation; (b) evolution of repetition rate difference under slow pump modulation; (c) DC interferogram under 1 kHz modulation; interferogram during (d) forward and (f) backward temporal scanning; RF spectrum obtained from (e) forward and (g) backward interferogram.

    Figure 3.Basic performance of DCANDi laser. (a) Optical spectra of the two combs without modulation; (b) evolution of repetition rate difference under slow pump modulation; (c) DC interferogram under 1 kHz modulation; interferogram during (d) forward and (f) backward temporal scanning; RF spectrum obtained from (e) forward and (g) backward interferogram.

    In the previous section, using simulation, we identified the unequal ωΔ in the two directions coupled with β2 as the main reason behind the dependence of Δfrep on pump power change. Here, we experimentally confirm this by measuring the evolution of the pulse spectrum in both directions during the pump-power modulation using the dispersive Fourier transform (DFT) technique. More information on measurement setup can be found in Appendix F. The modulation frequency was set at 5 kHz. The measured spectral evolution is shown in the appendix, while the δωΔ calculated from the DFT measurements are shown in Fig. 4(a). Figure 4(b) shows the expected Δfrep calculated based on Eq. (2) and data from Fig. 4(a). The β2 coupled with δωΔ results in Δfrep change from 11 to +11  Hz. This matches with magnitude of dynamic modulation in Δfrep achieved experimentally, as shown in the frequency counter measurement in Fig. 3(b). This validates the results from simulation and confirms that the unequal ωΔ change due to pump-power change contributes the most in dynamic modulation of Δfrep.

    (a) Evolution of the δωΔ in DCANDi calculated from the spectral evolution measured using DFT; (b) evolution of Δfrep due to the β2 coupled with δωΔ calculated using Eq. (2).

    Figure 4.(a) Evolution of the δωΔ in DCANDi calculated from the spectral evolution measured using DFT; (b) evolution of Δfrep due to the β2 coupled with δωΔ calculated using Eq. (2).

    It should be noted that the difference in the modulation depth of Δfrep in simulation and experiments is due to the difference in the mode-locking state in the experiment and simulation and also due to the fact that the simulation assumes a perfectly symmetric cavity, while it could be slightly asymmetric in the real experiments due to human error. However, the fact that dynamic modulation can be achieved even in a perfectly symmetric cavity, as shown through simulation results, rules out the necessity of cavity asymmetry (asymmetry in passive fiber length) for creating unequal response of the two directions to pump power and strengthens our understanding of spectral shifting playing the major role in this process. Since the asymmetric gain distribution along the length of the YDF in the two directions is the primary cause for realizing the dynamic tunable DC presented in this study, such dynamic modulation can be achieved in any symmetric bidirectional SCDCL systems where such an imbalanced longitudinal gain in the gain medium can be achieved. This is another advantage of bidirectional geometry compared to dual-wavelength and dual-polarization SCDCL systems, along with lesser cross talk and overlapping spectra. Hence, this technique is widely applicable in bidirectional SCDCL geometries, such as soliton fiber DCs and free-space DCs, and is not limited to CANDi lasers.

    5. NOISE CHARACTERIZATION OF DCANDI

    To ensure the enhanced frame rate in DCANDi does not come at the expense of CANDi’s timing stability, we use DFT-based spectral interferometry to measure the relative timing jitter during dynamic modulation with femtosecond resolution. Experimental details on the measurement technique can be found in the appendix and in Ref. [26]. The spectral interferogram formed by the two separation-evolving comb pulses is shown in Fig. 5(a). The two positions around 140 and 290 μs, where the fringe density drops to zero, correspond to the times when the two comb pulses meet each other during the forward and backward scanning. By performing an FFT on the spectral fringes in each round trip, the evolution of pulse separation is reconstructed, as shown in Fig. 5(b). Furthermore, the corresponding Δfrep can be obtained by taking the derivative of separation evolution in neighboring round trips, as shown in Fig. 5(c). The local fluctuation at 140 and 290 μs was caused by the separation-retrieving error when the two pulses overlap with each other. The near-zero spectral fringe density makes the separation calculation inaccurate. Nevertheless, Figs. 5(b) and 5(c) clearly show the separation and Δfrep evolution during the switching event. As observed, the two pulses first walk off with each other with constant speed before 195 μs. The corresponding Δfrep was about 11 Hz. Then, the pump power dropped, which drove the Δfrep to the negative value. The transition process took around 90 μs, after which the Δfrep appeared stable at 11  Hz. The duration of transition process is attributed to the laser response bandwidth, which is around 10 kHz. This also sets the upper speed limit for the dynamic Δfrep modulation. To verify the laser stability during dynamic operation, we calculate the relative timing jitter during 1000 consecutive round trips (48 μs) before and after the transition stage [marked by the red and blue dashed areas in Fig. 5(b)]. The corresponding power spectral densities (PSDs) and the corresponding integrated jitters are shown in Figs. 5(d) and 5(e), respectively. As observed, the PSDs for both time spans appear almost identical, and the integrated jitter (from 24 MHz to 20 kHz) is only around 5 fs, which is as stable as a conventional CANDi laser without modulation [26]. This result confidently proves that the dynamic modulation does not sacrifice the relative timing stability, thus providing a powerful technique for significantly enhancing the DC detection speed using SCDCL.

    Precise characterization of DC time-delay evolution during dynamic modulation. (a) Spectrogram during one Δfrep switching; (b) pulse separation evolution obtained by performing FFT on (a); (c) Δfrep evolution obtained from (b); (d) PSD of relative timing jitter calculated from 1000 consecutive round trips before (red) and after (blue) the transition, as marked in (b); (e) integrated timing jitter calculated from (d).

    Figure 5.Precise characterization of DC time-delay evolution during dynamic modulation. (a) Spectrogram during one Δfrep switching; (b) pulse separation evolution obtained by performing FFT on (a); (c) Δfrep evolution obtained from (b); (d) PSD of relative timing jitter calculated from 1000 consecutive round trips before (red) and after (blue) the transition, as marked in (b); (e) integrated timing jitter calculated from (d).

    6. FAST AND HIGH-FIDELITY NONLINEAR METROLOGY

    To demonstrate the application of high frame rate and temporal resolution of the proposed dynamic Δfrep modulation technique, we characterize the relaxation time constant of an off-the-shelf SESAM sample (from BATOP) using the DCANDi laser. The pump fluence incident on the SESAM is 58  μJ  cm2, while the probe fluence is 0.5  μJcm2. The pump-probe experiment was done in a reflection geometry, as shown in Fig. 6(a). The reflected probe was collected using a Si-switchable gain detector (Thorlabs PDA100A). The relaxation time constant is characterized in both conventional DC pump-probe and the dynamic Δfrep modulation DC pump-probe methods. The results are shown in Figs. 6(b)–6(e). In the conventional DC method, where the Δfrep was fixed at 14 Hz (temporal sampling time step is 6  fs), Fig. 6(b) shows the collected pump-probe signal for 600 ms (nine pump-probe curves) and Fig. 6(c) shows the averaged SESAM response, which reveals a time constant of 5.18±0.18  ps. The same experiment was repeated using the dynamic modulation of Δfrep at 1 kHz. Figure 6(d) shows the measured pump-probe signal using dynamic Δfrep modulation technique for 5 ms. The temporal sampling time step is maintained at 6  fs. Figure 6(e) shows the average SESAM response, revealing a fitted time constant of 5.19±0.22  ps. The number of pump-probe curves used for calculating the average response is kept at nine for both the techniques to ensure the same SNR.

    Characterizing the relaxation time constant of SESAM. (a) Experimental setup. Pump-probe technique is employed in a reflection geometry. (b) Measured pump-probe signal using conventional DC technique for 600 ms; (c) mean response of the SESAM calculated by averaging nine pump-probe signals in (b); (d) measured pump-probe signal using proposed dynamic Δfrep modulation technique for 5 ms; (e) mean response of the SESAM calculated by averaging nine pump-probe signals in (d).

    Figure 6.Characterizing the relaxation time constant of SESAM. (a) Experimental setup. Pump-probe technique is employed in a reflection geometry. (b) Measured pump-probe signal using conventional DC technique for 600 ms; (c) mean response of the SESAM calculated by averaging nine pump-probe signals in (b); (d) measured pump-probe signal using proposed dynamic Δfrep modulation technique for 5 ms; (e) mean response of the SESAM calculated by averaging nine pump-probe signals in (d).

    Figure 6 confirms a frame rate enhancement of 143 (2000 Hz/14 Hz) without sacrificing resolution, sensitivity, and accuracy. Of note, while similar enhancement has been demonstrated in two-laser dual-comb systems [11,13], this is the first dynamic Δfrep modulation demonstration using an SCDCL. The inherent common-mode-noise cancellation of SCDCLs helps in passively reducing relative timing jitter, which opens up the possibility for low-complexity dual-comb pump-probe measurements. Apart from frame rate enhancement, as a byproduct, assuming the noise is uncorrelated, one can also expect to obtain a higher SNR in the proposed dynamic Δfrep modulation technique for same averaging time. From the standard deviation of the fitted response in both cases, it is also clear that dynamic modulation measurement has a similar error (or SNR) compared to the conventional case. Hence, the suggested technique can be used for employing higher frame rate or higher SNR depending on application.

    7. CONCLUSION AND OUTLOOK

    In summary, we demonstrate an all-optical technique to break the measurement speed-time resolution trade-off in SCDCL dual-comb pump-probe measurements. We achieve the all-optical dynamic Δfrep modulation in a dual-pumped CANDi fiber laser by applying pump-power modulation to one of the pumps and showcase its basic operation. We present a simulation routine to accurately model the complex interplay of parameters governing bidirectional SCDCLs and use it to identify the unequal change in center frequency due to the asymmetric longitudinal gain distribution in the YDF as the mechanism of dynamic Δfrep modulation of the DCANDi fiber laser. In addition, we confirm that the demonstrated all-optical dynamic Δfrep modulation technique does not introduce excessive relative timing jitter and identify that the measurement speed is limited by the laser response bandwidth to 10 kHz. As a proof-of-principle experiment, we apply the DCANDi fiber laser to measure the relaxation time of an off-the-shelf SESAM and achieve a 143-fold enhancement in measurement speed and duty cycle. The DCANDi fiber laser is a promising light source for low-complexity, high-speed, high-sensitivity ultrafast DC pump-probe measurements. The working principle can be generalized to other SCDCL platforms where gain asymmetry in the two directions dominates the Δfrep, thus representing a universal solution for breaking the measurement speed-time resolution trade-off in bidirectional SCDCLs. In particular, gain media with a shorter upper-state lifetime, such as highly doped Yb multicomponent glass [29], phosphate glass fiber [30], and Ti:sapphire crystal [31] can be utilized to further enhance the dynamic modulation speed above 10 kHz.

    Acknowledgment

    Acknowledgment. The authors acknowledge the Optics and Photonics Research Group at the University of Colorado Boulder supervised by Dr. Juliet Gopinath for providing the SESAM used in this work.

    APPENDIX A: NUMERICAL SIMULATIONS

    Simulation Routine

    The numerical simulation for the DCANDi laser with gain dynamics is based on a model that jointly solves the rate equations (REs) and a coupled set of Ginzburg–Landau equations (CGLE), using a split-step Fourier method (SSFM). Limiting our model to include dispersion up to the third order and ignoring Raman effects, our set of coupled equations describing propagation takes the form δUδz=Δβ1δUδtiβ22δ2Uδt2iβ36δ3Uδt3αg2U+iγ[(|U|2+23|V|2)+13U*V2],δVδz=Δβ1δVδtiβ22δ2Vδt2iβ36δ3Vδt3αg2V+iγ[(|V|2+23|U|2)+13V*U2].U and V represent the electric fields in orthogonal polarizations, Δβ1 is the polarization mode dispersion, β2 is the group velocity dispersion (GVD), and β3 is the third-order dispersion (TOD). α is the intrinsic fiber loss. γ=2πn2/λAeff is the nonlinear coefficient calculated using the effective mode area, Aeff, and the nonlinear index n2. g represents the gain of the active fiber. For propagation through the segments of passive fiber, Eq. (A1) is sufficient to model the CW and CCW propagation directions. Without loss of generality, the following simulations ignored polarization dispersion as well as the intrinsic fiber loss (α=Δβ=0).

    To capture the gain dynamics in the active fiber, the frequency-dependent gain is calculated during the linear step of the SSFM using the pulse spectra and employing the following set of REs valid for high repetition pulse amplification [32,33]: δPpδz=Γp[σeN2(z)σaN1(z)]ρPp,δPsδz=Γs[σeN2(z)σaN1(z)]ρPs,N2(z)=R12+W12R12+R21+W12+W21+1τ21.

    The spectroscopic data from Ref. [34] were used to calculate the emission (σe) and absorption (σa) cross sections. Γp(s)=1exp(2rcore2/w2) is the modal overlap of the signal, where rcore is the core radius and w is the mode field radius at 1/e2 approximated by the Whitely model [35]. The population inversion is calculated with R12 (W12) and R21 (W21) representing the stimulated absorption and emission rate for the pump (signal). For a fiber core-area A, Rij=Γpσe,a(Pp++Pp)/ωpA and Wij=Γsσe,aPs±/ωsA, where Pp+ represents the copropagating pump and Pp is the counterpropagating pump, while Ps± represents signal power in CW and CCW directions. Finally, Eqs. (A1) and (A2) are coupled through the gain coefficient g˜=g˜0/[1+(Ps++Ps)/Psat], where Psat is the saturation power and g˜0 is the small signal gain calculated at each step, δz, of the fiber using the following relation: g˜0=1δzln[Ps(z+δz)Ps(z)].

    The four wave plates in the cavity are modeled with two sets of polarization controllers (PCs). The transfer function of both PCs for the CW direction is represented as TPC=(cosθsinθsinθcosθ)(eiϕ/200eiϕ/2)(cosθsinθsinθcosθ).θ is the angle between the coordinates of the PC and the polarization beam splitter, and ϕ is the phase delay between the two orthogonal polarization components induced by the PC.

    The logic of the simulation is illustrated in the flow chart in Fig. 7, and Table 1 shows the parameters and their values used in the simulation.

    Logic of the DCANDi simulation.

    Figure 7.Logic of the DCANDi simulation.

    (a) Evolution of center wavelength for direction 1 (black) and direction 2 (red) for the case of Δfrep=−2.5 Hz; (b) comparison of the saturated gain spectrum of direction 1 (black) and direction 2 (red) at different points along the YDF marked in (a), which are the beginning, middle, and end of the YDF. The insets are zoomed-in to the gain peaks, and the dotted lines in the inset denote the center of the gain peak for each direction.

    Figure 8.(a) Evolution of center wavelength for direction 1 (black) and direction 2 (red) for the case of Δfrep=2.5  Hz; (b) comparison of the saturated gain spectrum of direction 1 (black) and direction 2 (red) at different points along the YDF marked in (a), which are the beginning, middle, and end of the YDF. The insets are zoomed-in to the gain peaks, and the dotted lines in the inset denote the center of the gain peak for each direction.

    (a) Evolution of center wavelength for direction 1 (black) and direction 2 (red) for the case of Δfrep=+2.5 Hz; (b) comparison of the saturated gain spectrum of direction 1 (black) and direction 2 (red) at different points along the YDF marked in (a), which are the beginning, middle, and end of the YDF. The insets are zoomed-in to the gain peaks, and the dotted lines in the inset denote the center of the gain peak for each direction.

    Figure 9.(a) Evolution of center wavelength for direction 1 (black) and direction 2 (red) for the case of Δfrep=+2.5  Hz; (b) comparison of the saturated gain spectrum of direction 1 (black) and direction 2 (red) at different points along the YDF marked in (a), which are the beginning, middle, and end of the YDF. The insets are zoomed-in to the gain peaks, and the dotted lines in the inset denote the center of the gain peak for each direction.

    Figure 9 shows the center wavelength evolution and gain curve for the case of positive Δfrep (+2.5  Hz). Figure 9(a) shows the evolution of the center wavelength inside the cavity where the evolution behavior is different from Fig. 8(a). In this case, the center wavelength of direction 1 is larger than that of direction 2 both inside the passive fibers and after the Gaussian filter. However, within the YDF, although the center wavelengths of both directions are blueshifted in the initial length of the fiber, the center wavelengths have two crossing points where the center wavelength of direction 1 becomes smaller than that of direction 2 after the first crossing and becomes larger again once the center wavelengths start to redshift toward the back end of the YDF. The saturated gain profiles of both directions at beginning (left, blue); middle (center, green); and end (right, magenta) of YDF are shown in Fig. 9(b). Similar to the earlier case, for both directions, the gain peak is on the shorter wavelength side in the beginning part of the YDF and gets pushed to the red side (closer to the DCANDi spectral region) due to gain saturation as the pulses are amplified inside the YDF. Even though the gain saturation results in lesser gain in direction 1 compared to direction 2, the peak gain wavelength [dotted lines in the insets of middle and right plots in Fig. 9(b)] for the two directions is almost similar at the middle of the YDF (closer to the second crossing). However, by the end of the YDF, the peak gain wavelength of direction 1 is redder compared to direction 2.

    APPENDIX B: PULSE COMPRESSION

    Both the combs from DCANDi are compressed using a common grating compressor in the transmission geometry. The grating compressor has a transmission of 85%. The compressed pulses are characterized using a conventional FROG setup; the results are shown in Fig. 10. The retrieved pulse duration, assuming a Gaussian pulse shape, of comb1 is about 171 fs, while that of comb2 is 185 fs. The slight difference in pulse duration could be attributed to the slight difference in spectral width and accumulated nonlinear phase along both directions. Note that the error in reconstructed spectrogram is less than 1% for both the cases.

    Characterization of compressed pulse using a conventional FROG setup. In the figure, (a) and (b) are the measured and reconstructed spectrogram of comb 1, with (c) showing the retrieved pulse shape of comb 1. The Gaussian fitted pulse duration is 171 fs. Similarly, (d) and (e) are the measured and reconstructed spectrogram of comb 2, with (f) depicting the retrieved pulse shape of comb 2. The Gaussian fitted pulse duration is 185 fs.

    Figure 10.Characterization of compressed pulse using a conventional FROG setup. In the figure, (a) and (b) are the measured and reconstructed spectrogram of comb 1, with (c) showing the retrieved pulse shape of comb 1. The Gaussian fitted pulse duration is 171 fs. Similarly, (d) and (e) are the measured and reconstructed spectrogram of comb 2, with (f) depicting the retrieved pulse shape of comb 2. The Gaussian fitted pulse duration is 185 fs.

    APPENDIX C: DCANDI RESPONSE TO PUMP POWER

    The dynamic modulation of Δfrep in DCANDi is realized by modulation of one of the pump laser powers (pump2) while keeping the other pump power constant. Figure 11(a) shows the response of the Δfrep to change in pump2 power. Note that the pump1 power is kept constant at 1.26 W, and measurement of Δfrep is done using frequency counters. The response is almost linear, with a slope of 0.2  Hz/mW and a tuning range of 34 Hz for a power change of 168 mW. These values depend on the mode-locking position, as different mode-locking positions can result in different slopes and tuning ranges. Note that, for the depicted mode-locking state, Δfrep can be further increased by increasing the pump power up to 1.54 W without disturbing the mode locking. Above 1.54 W and below 1.26 W, the bidirectional mode-locking state is disturbed. Figure 11(b) shows the corresponding response of the output power in the two directions to the change in the pump2 power. Direction 2 (red) has smaller slope (0.036 mW/mW) compared to direction 1 (black) (0.062 mW/mW), with an overall change of 6.5% and 11.5% in output power for a change of 160 mW pump power. Note that these are the values for the whole tuning range. For a change of Δfrep=±10  Hz, the change in output power is lower (3.3% and 5.2%, respectively).

    (a) Response of Δfrep to change in pump2 power. (b) Change in output power for direction 1 (black) and direction 2 (red) as a function of change in pump2 power. Note that the pump1 power is fixed at 1.26 W for this measurement.

    Figure 11.(a) Response of Δfrep to change in pump2 power. (b) Change in output power for direction 1 (black) and direction 2 (red) as a function of change in pump2 power. Note that the pump1 power is fixed at 1.26 W for this measurement.

    APPENDIX D: RESPONSE OF PUMP POWER TO SQUARE-WAVE MODULATION

    Figure 12 shows the response of the pump power (measured using a photodiode) to square modulation that is applied to DCANDi for dynamic modulation of Δfrep. Figures 12(a) and 12(b) show the response for 1 kHz modulation frequency and 5 kHz modulation frequency, respectively. The insets in both figures show a zoomed-in picture of the first rising edge of the square wave. In both cases, the modulation amplitude is kept constant at 250 mV. The response shows that the pump power follows the modulation wave well, at least until 5 kHz. However, the zoomed insets reveal that the pump power undergoes an oscillation relaxation-like behavior until 50  μs before settling down. This implies that the maximum modulation frequency required to achieve a stable pump power is limited to <10  kHz (considering the full period of both high and low of the square wave). This implies a maximum pump-limited frame rate (or half-cycle) of 20 kHz in a DCANDi. However, as discussed in the main text, currently, the maximum frame rate of DCANDi is limited by the laser response to 10  kHz.

    Modulation of pump power by a square wave of (a) 1 kHz and (b) 5 kHz frequency. The insets show the zoomed-in picture of the first rising edge of the square wave. The amplitude of modulation is 250 mV.

    Figure 12.Modulation of pump power by a square wave of (a) 1 kHz and (b) 5 kHz frequency. The insets show the zoomed-in picture of the first rising edge of the square wave. The amplitude of modulation is 250 mV.

    APPENDIX E: LOW-FREQUENCY PID LOOP

    In order to make sure the absolute value of Δfrep is the same before and after flipping, a slow PLL with 100 Hz control bandwidth is applied to pump1 to lock the average Δfrep to zero. Since the control bandwidth is much smaller than the modulation frequency (a few kilohertz), the PLL does not counteract the modulation. To form the PLL, 1% of output power from both combs is fed into two photodetectors, respectively, to obtain the fundamental repetition frequency, which are then frequency-mixed. The output of the mixer serves as the error signal, which is seeded into a PID controller that generates the control signal for controlling the current of pump1. The PLL on the pump1 power is obtained by using only the proportional and integral options of the PID control in the Moku:Lab, with an overall gain of +21  dB, proportion gainof12  dB, and integral gainof50  dB. This results in a frequency bandwidth of 100  Hz. Figure 13(a) shows the unlocked and locked error signal collected using an oscilloscope. With the locking, the slow frequency is suppressed and the standard deviation of error signal is reduced from 0.032 to 0.00006 V. The FFTs of the error signals are shown in Fig. 13(b). This shows that the noise until 100 Hz is well suppressed by the PLL, while the 2 kHz signal is not suppressed, implying that PLL does not affect the pump-power modulation operation.

    (a) Error signal of the unlocked (blue) and locked (red) case measured with an oscilloscope; (b) FFT of the error signals in (a).

    Figure 13.(a) Error signal of the unlocked (blue) and locked (red) case measured with an oscilloscope; (b) FFT of the error signals in (a).

    APPENDIX F: DFT MEASUREMENT

    The DFT measurement setup comprises a 10 km HI1060 fiber spool providing a group delay dispersion of 230  ps2, followed by an ytterbium-doped fiber amplifier (YDFA) stage for amplification, a fast photodetector with bandwidth of 12.5 GHz, and a high-speed real-time oscilloscope (10 GHz). The net pulse energy launched into the setup is limited to 25  pJ to avoid nonlinearity in the fiber. For measuring the evolution of the pulse spectrum in both directions during the pump power modulation, a pulse from each individual direction was sent to the DFT fiber spool alternatively while pump2 was modulated. For measuring the relative timing jitter of DCANDi, both directions are combined using a 50:50 coupler that has a polarization controller installed in one of the arms to align the polarization of the two pulses. One of the output arms of the 50:50 coupler is sent to the DFT fiber spool.

    In order to experimentally validate the findings from simulation on the origin of dependence of Δfrep on pump-power change, we measured the evolution of the pulse spectrum in both directions during the pump-power modulation using the DFT technique. The spectral evolution of the two combs when pump2 power is modulated at 5 kHz frequency is shown in Figs. 14(a) and 14(b). Center wavelengths are calculated for each comb at each round trip, and these values are used in Eq. (2) in the main text to calculate the plots in Fig. 4(a) and subsequently Fig. 4(b). As observed, comb1 is more sensitive to pump-power change compared to comb2.

    Normalized spectral evolution along (a) direction 1 and (b) direction 2 while the pump2 is modulated, measured using the DFT technique.

    Figure 14.Normalized spectral evolution along (a) direction 1 and (b) direction 2 while the pump2 is modulated, measured using the DFT technique.

    [1] A. M. Weiner. Ultrafast Optics(2011).

    [31] H. Jelínková. Lasers for Medical Applications: Diagnostics, Therapy and Surgery(2013).

    Tools

    Get Citation

    Copy Citation Text

    Neeraj Prakash, Jonathan Musgrave, Bowen Li, Shu-Wei Huang, "Dynamic counterpropagating all-normal dispersion (DCANDi) fiber laser," Photonics Res. 12, 2033 (2024)

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Lasers and Laser Optics

    Received: Apr. 29, 2024

    Accepted: Jul. 4, 2024

    Published Online: Sep. 2, 2024

    The Author Email: Shu-Wei Huang (ShuWei.Huang@colorado.edu)

    DOI:10.1364/PRJ.528873

    Topics