## Dissipative discrete time crystals in a pump-modulated Kerr microcavity

### Model and system symmetry

If the CW laser driving a Kerr resonator is phase-modulated at frequency *ω*_{M}, sidebands will be generated symmetrically with respect to this pump. Depending on *ω*_{M}/*D*_{1}, the ratio of the modulation frequency to cavity FSR, and *δ*_{M}, the modulation depth or intensity, modulation-induced sidebands will have different effects on the frequency comb generation process. It has been shown that seeding microresonator-based Kerr comb (microcomb) generation^{20} by a modulated CW laser can improve comb stability^{21,22} and deterministically create a train of temporal dissipative Kerr solitons^{23,24,25,26,27,28}. Pump modulation at the cavity FSR (*M* = 1) with a stable “clock” can enhance the long-term stability of the microcomb and improve the spectral purity and phase noise of the microwave signal^{22,23}. Additionally, modulation depth provides added control over the beatnote signal power while system dynamics and attractors are different from those in a cavity pumped by a free-running laser.

In the current study, we focus on the case of the modulation frequency being an integer multiple *M* > 1 of the FSR, i.e., *ω*_{M} = *M**D*_{1}. In this notation, *ω*_{M} is an angular frequency in rad/s and *D*_{1}/2*π* = 1/*T*_{R} denotes the FSR in Hz and the reciprocal of the round-trip time *T*_{R}^{29,30}. Similar to DTC formation in a dichromatically pumped Kerr cavity^{12,31,32,33}, drive periodicity defines a discrete temporal symmetry exhibited in the modulation of the intra-cavity CW background (Fig. 1b), while the spontaneous appearance of one or more temporal dissipative solitons imprints a larger time period upon the system output (Fig. 1c, d). Therefore, at any fixed position on the output waveform path, photon count probability will demonstrate a periodicity that is an integer multiple of that of the same observable monitored before the resonator. As a result, discrete TTSB accompanied by characteristic *subharmonic* generation in the frequency domain, occurs.

The time-domain description of the system in the laboratory reference frame takes the form

$$\frac{\partial A}{\partial t}= \,\left(\frac{-\kappa }{2}-{{{{{{{\rm{i}}}}}}}}\sigma +\mathop{\sum }\limits_{m = 1}^{N}{{{{{{{{\rm{i}}}}}}}}}^{m+1}\frac{{D}_{m}}{m!}\frac{{\partial }^{m}}{\partial {\theta }^{m}}-{{{{{{{\rm{i}}}}}}}}g| A{| }^{2}\right)A\\ +\sqrt{{\kappa }_{{{{{{{{\rm{c}}}}}}}}}}\,{F}_{0}\exp \left[{{{{{{{\rm{i}}}}}}}}{\delta }_{{{{{{{{\rm{M}}}}}}}}}\sin ({\omega }_{{{{{{{{\rm{M}}}}}}}}}t-M\theta )\right],$$

(1)

where the complex waveform *A*(*t*, *θ*) is normalized such that ∫_{2π}∣*A*∣^{2} d*θ* renders the number of intra-cavity photons, *κ* denotes the power decay rate, *κ*_{c} is the coupling coefficient, \(\sigma ={{\Omega }}-{\omega }_{{j}_{0}}\) indicates the detuning between the CW laser and pumped resonance frequencies Ω and \({\omega }_{{j}_{0}}\), *D*_{m} for 1 ≤ *m* ≤ *N* are dispersion coefficients (*m* and *N* both integers), *θ* is the azimuthal angle around the resonator, related to the fast time *τ* via *θ* = 2*π**τ*/*T*_{R} (modulo 2*π*), *g* is the four-wave mixing (FWM) gain or nonlinear coupling coefficient, and *F*_{0} represents the laser pump power^{17,20} (see the “Methods” section). Equation (1) is invariant under *t* → *t* + 2*π*/*ω*_{M} or *θ* → *θ* + 2*π*/*M* (equivalently, *τ* → *τ* + 2*π*/*ω*_{M}), and therefore possesses time-translation symmetry defined by the drive. Soliton formation can cause discrete TTSB and DTC formation for *M* > 1, because the drive is periodic by *T*_{Drive} = 2*π*/*ω*_{M} = *T*_{R}/*M* while the pulse train exiting the resonator has a larger periodicity *T*_{Response} = *n**T*_{Drive}, the integer *n* satisfying 1 < *n* ≤ *M*.

Equation (1) can readily be recast by the change of variable *θ* → *θ*−*ω*_{M}*t*/*M*, equivalent to transitioning to a reference frame rotating with angular velocity *ω*_{M}/*M*. The resulting equation is better conducive to numerical modeling and upon normalization reads

$$\frac{\partial \psi }{\partial \bar{t}}= \,\left[-1-{{{{{{{\rm{i}}}}}}}}\alpha +{{{{{{{\rm{i}}}}}}}}| \psi {| }^{2}+\frac{2}{\kappa }\left(\frac{{\omega }_{{{{{{{{\rm{M}}}}}}}}}}{M}-{D}_{1}\right)\frac{\partial }{\partial \theta }\right.\\ -\left.\mathop{\sum }\limits_{m=2}^{N}{\left(-{{{{{{{\rm{i}}}}}}}}\right)}^{m+1}\frac{{d}_{m}}{m!}\frac{{\partial }^{m}}{\partial {\theta }^{m}}\right]\psi +{f}_{0}\exp \left[{{{{{{{\rm{i}}}}}}}}{\delta }_{{{{{{{{\rm{M}}}}}}}}}\sin (M\theta )\right].$$

(2)

The field envelope *ψ* = *A*^{*}/*A*_{th} and external pump amplitude \({f}_{0}=\left[\sqrt{{\kappa }_{{{{{{{{\rm{c}}}}}}}}}}/(\kappa /2)\right]\,{F}_{0}^{* }/{A}_{{{{{{{{\rm{th}}}}}}}}}\) have been normalized with respect to the modulus of the comb generation threshold \(| {A}_{{{{{{{{\rm{th}}}}}}}}}| =\sqrt{\kappa /(2g)}\)^{34}; note the complex conjugation in going from Eq. (1) to Eq. (2)^{35}. The slow time \(\bar{t}=t/(2{t}_{{{{{{{{\rm{ph}}}}}}}}})\) has been normalized by twice the pumped resonance photon life-time *t*_{ph} = *κ*^{−1}, while detuning *α* = −2*σ*/*κ* and dispersion coefficients *d*_{m} = −2*D*_{m}/*κ* have been re-scaled through division by the pumped mode half-width at half-maximum (HWHM). We will base our subsequent numerical modeling and stability analysis using the split-step Fourier transform method in Eq. (2), which is a damped and driven nonlinear Schrödinger or a modified Lugiato–Lefever equation (NLSE or LLE, respectively)^{36}. In the context of Bose–Einstein condensates, the NLSE is often called the Gross–Pitaevskii equation. The split-step Fourier transform method has been shown to offer excellent accuracy in modeling Kerr microcomb and fiber cavity soliton systems, and to agree well with experiments^{29,35,37}.

Before proceeding, we estimate parameters of a resonator suitable for the demonstration of the proposed phenomena. Consider a whispering gallery mode (WGM) resonator made out of magnesium fluoride (MgF_{2}), interrogated with a 1.5 μm laser. Accounting for geometrical and material dispersion, a resonator with a 10 GHz FSR is characterized with the dispersion parameter \({D}_{2}={\omega }_{{j}_{0}+1}+{\omega }_{{j}_{0}-1}-2{\omega }_{{j}_{0}}\approx 0.7\) kHz. With a HWHM modal linewidth of 150 kHz, this translates into *d*_{2} = −5 × 10^{−3}, and comb generation power threshold falls in the few mW range. Loading of the cavity can be modified by a few orders of magnitude through controlling the gap between the evanescent field coupler and the resonator surface as well as by changing the coupling element material (the prism coupler). Besides MgF_{2}, fused silica or another transparent optical material can be utilized. These estimations show that the experiments proposed here can be readily achieved in practice. Additionally, pump modulation at the FSR frequency (*M* = 1) for MgF_{2} resonators with FSRs as low as 12 and 14 GHz have been reported in recent years^{22,25}. Subharmonic generation, one of the criteria for DTCs, can in principle be realized also in normal dispersion^{32,38}.

### DTC waveforms

Two sets of examples are shown in Fig. 2 in a resonator with *d*_{2} = −6.3 × 10^{−3} (anomalous dispersion), pump intensity *f*_{0} = 1.275 and pump-resonance detuning *α* = 2 (red detuned), modulation frequency *ω*_{M} = 4*D*_{1} (i.e., *M* = 4), and two different modulation depths *δ*_{M} = 0.7 (Fig. 2a–d) and *δ*_{M} = 0.3 (Fig. 2e–h). The latter case of *δ*_{M} = 0.3 corresponds to the same set of numerical values also used in producing Fig. 1b–d. It is clear that without soliton formation, the temporal waveform leaving the resonator is periodic with *T*_{Drive} = *T*_{R}/4 (see Fig. 2c, d and g, h). However, when either one or three solitons are trapped in the lattice created by pump modulation, the periodicity of the soliton train changes to *T*_{Response} = *T*_{R} = 4*T*_{Drive}, thereby realizing 4-DTCs. The examples shown in Fig. 1d corresponds to the case of *T*_{Response} = *T*_{R}/2 = 2*T*_{Drive}. If four pulses emerge and all lattice spots per round-trip time are occupied with soliton peaks the discrete symmetry of the system will be respected and TTSB will not happen. With larger *M*, other DTC sizes (response to drive ratios) are possible. For instance, using the same parameters at *δ*_{M} = 0.7, we observed discrete TTSB with *M* up to 19. The pulses emerge from a random high-energy initial state (i.e., through hard excitation^{39}) and propagate stably for as long as the drive is on; the integration time in our numerical simulations have been hundreds of cavity photon lifetimes.

As the modulation frequency is increased successively, lattice sites created on the modulated CW background will be more packed and shorter pulses ought to be generated. While solitons tend to become shorter through increasing the pump-resonance detuning, packing more solitons in one round-trip naturally translates into a need for cavity dispersion engineering to reduced the magnitude of *D*_{2}^{40}. Furthermore, from an experimental perspective, faster modulators will be required. As low-*V*_{π} commercial modulator frequencies are currently bound to tens of GHz (typically below 100 GHz), demonstration of large dissipative DTCs in this system requires resorting to longer resonators. Increasing the resonator size in WGM crystalline resonators entails competing with overmodedness^{41} and avoided mode crossings^{42} which are known to challenge soliton formation. This issue can be more easily addressed in integrated platforms^{17} where racetrack resonators or those with wedge-shaped rims have been shown to support repetition rates down to a few GHz^{43}. We have shown elsewhere that pumping with two independent lasers can also create dissipative DTCs without suffering from the constraints of off-the-shelf high-frequency modulators^{12}. While integrated, fast, and low-*V*_{π} opto-electronic modulators are being developed in the integrated photonics community^{44}, the principles developed here can also apply to fiber resonators^{24}, where the cavity round-trip time is much longer (typically 3–4 orders of magnitude) compared to microresonators.

### Stability analysis

What is the relationship between the modulation frequency and the number and temporal separation of stable solitons leaving the resonator per round-trip time? Using the method of moments and in analogy with the notion of particle momentum in quantum mechanics, a formalism revealing the force exerted by pump phase modulation on dissipative solitons and the potential trapping them can be developed, as detailed in the “Methods” section. The equation of motion for the waveform momentum \({{{{{{{\mathcal{P}}}}}}}}\) is \({{{{{{{\rm{d}}}}}}}}{{{{{{{\mathcal{P}}}}}}}}/{{{{{{{\rm{d}}}}}}}}\bar{t}=-2{{{{{{{\mathcal{P}}}}}}}}+{{{{{{{\mathcal{F}}}}}}}}\), in which pulse and modulation parameters determine the force \({{{{{{{\mathcal{F}}}}}}}}\)^{31}. For modulation depth *δ*_{M} of order unity, this force will vanish at fixed points \({\theta }_{k}^{* }=(2k+1)\pi /(2M)\), with integer *k*. The potential, on the other hand, will assume the form \({{{{{{{\mathcal{V}}}}}}}}(\theta )=-a\sin (M\theta )+b\), where *a* and *b* are constants (see the “Methods” section). Hence, the potential is essentially the phase modulation profile with a negative sign, \(-\sin (M\theta )\), plus a constant offset not affecting stability. Potential minima (half of the fixed points) are the stable equilibria of the system.

These predictions are substantiated by the numerical integration of Eq. (2); at steady state, the soliton peaks always settle at the potential profile minima, as seen in Fig. 3a. The soliton peaks in Fig. 3a illustrate the final state of the pulse train in Fig. 3b, where Eq. (2) has been integrated with a random high-energy initial state and at *ω*_{M} = *M**D*_{1}, leading to the appearance of multiple peaks, all of which are dragged to the stable equilibria. Non-ideal locking of the modulation frequency to the pumped FSR has been assumed in Fig. 3c by considering (*ω*_{M}/*M*−*D*_{1}) to be a noisy signal with a uniform distribution of zero mean and 2% HWHM standard deviation. As the plot shows, it may take longer for the solitons to merge, but the steady state remains settling to the potential minima. Also, with sudden large-enough disruptions of the pump conditions, other soliton peaks may appear in the cavity, but they will converge under the drag imposed by the modulation, again to the expected stable states (see Fig. 4). Note that in the final states in Figs. 3 and 4, discrete TTSB does not occur (both the drive and response have the same periodicity of *T*_{R}/4).

### Absence of DTCs for M = 1

Before proceeding to state control, we digress briefly to emphasize the significance of the foregoing analysis. Besides underpinning temporal long-range order in the proposed DTCs, this analysis highlights the central role of phase modulation at an integer multiple of the FSR in the realization of discrete TTSB in this platform. Without modulating the pump, photon detection probability in the microcomb drive possesses continuous time-translation symmetry and cannot demonstrate DTCs^{12}. However, when the modulation frequency matches the FSR (*M* = 1), the framework presented above illustrates the existence of two stable equilibria in the system, one with no solitons and the other with one soliton per round-trip time^{23,26}. Therefore, the drive introduces a discrete temporal symmetry which is not violated by soliton formation in steady state; even if multiple pulse peaks emerge initially, only one may ultimately survive such that response periodicity equals *T*_{R} both with or without soliton generation^{23}. Finally, when the pump is modulated at a frequency that is a larger-than-unity integer multiple of the FSR (*M* > 1), the discrete symmetry introduced by the drive can be violated through some distributions of soliton peaks over potential lattice sites to create DTCs. This dichotomy is well manifested in the intra-cavity energy ∫_{2π}∣*ψ*∣^{2} d*θ* versus detuning curves of Fig. 5. In this region, corresponding to normalized detuning roughly between 3 and 8 for the considered parameters in this example, each step of higher energy indicates a stable state with one soliton peak more per round-trip time than the adjacent lower-energy one, whereas the lowest-energy step signifies soliton-free states. In the absence of modulation, different numbers of solitons can be realized per round-trip time and the relative temporal delay between them is not dictated by the pump^{29} (Fig. 5a). With modulation at an integer multiple of the FSR, *ω*_{M} = *M**D*_{1}, a discrete time-translation symmetry is defined by the drive and the maximum number of soliton peaks as well as their relative time delay becomes constrained by *M*. For *M* > 1, TTSB is possible, as seen in Fig. 5b, where *M* = 4 and frequency comb states lying on middle steps (i.e., steps other than the lowest- and highest-energy) correspond to DTCs (cf. Fig 2a, e). In contrast, only two steady states are possible when *M* = 1 (i.e., there is no middle step between the top and bottom ones) and neither of the pertinent waveforms leads to TTSB (see Fig. 5c). DTCs can thus not form when *ω*_{M} = *D*_{1}.

### DTC generation and controlled phase switching

As the multistability embodied in the soliton steps in Fig. 5 signifies, the number of solitons per round-trip time forming spontaneously in the cavity cannot be accurately controlled. Based on the understanding furnished by the foregoing analysis, one can devise ways to realize desirable DTC states in a controlled manner. Two example scenarios are illustrated in Fig. 6, where *ω*_{M} = 4*D*_{1}. In Fig. 6a–c, the initial state is a 4-DTC. The modulation is shortly turned off and then back on, now at *ω*_{M} = 2*D*_{1}; see the modulation profile in Fig. 6c, at which stage the solitons are dragged towards the stable equilibria of the new phase-modulated system and the two previously centered at *θ*_{1} = *π*/8 and *θ*_{2} = 5*π*/8 merge into one. This procedure realizes a 2-DTC when the modulation frequency goes back to *ω*_{M} = 4*D*_{1}. The waveforms and frequency spectra of the initial (blue) and final (red) DTC states are plotted in Fig. 7a, b. Panels d–f of Fig. 6 depict another example for transitioning from a 2-DTC to another 4-DTC (this time one with one soliton per round-trip time). To achieve this, the modulation frequency is changed to *ω*_{M} = *D*_{1}, again resulting in the merging of the pulses into one before it resumes the original value of *ω*_{M} = 4*D*_{1}. Figure 7 shows the waveforms (c) and frequency spectra (d) of the initial (blue) and final (red) DTC states.

It should be stressed that for a system to host DTCs, first it ought to possess discrete time translation symmetry. This symmetry is absent in monochromatically driven Kerr microcombs, even in the presence of the 3rd- and higher-order dispersion or avoided mode crossings; this has been discussed in detail elsewhere^{12}. As a result, dispersive wave emission and soliton crystals in CW-pumped Kerr microresonators do not constitute DTCs. We note, however, that if effects such as higher-order dispersion or Raman-induced frequency shift change the repetition rate with respect to the FSR, the modulation frequency in the current DTC generation scheme can readily be adjusted to compensate for soliton recoil and dominate the dynamics, and this effect is well captured by Eq. (2). An example is depicted in Fig. 8. Non-negligible *D*_{3} can create dispersive waves on one side of the pump and cause soliton recoil to the opposite side^{40,45}. Therefore, the soliton train will start moving in the rotating reference frame of Eq. (2), such that it will generally be more difficult to trap and merge solitons with *ω*_{M} = *M**D*_{1} (see Fig. 8a–c), where all parameters are the same as those in Fig. 6 except for a non-zero third-order dispersion coefficient. Phase modulation tends to reduce the temporal separation of the solitons when they are in the vicinity of a potential minimum but increase it as one of the pulses drifts away to a neighboring one, as observed in Fig. 8a. The successive increase and decrease of temporal separations, and the accompanying changes of soliton interaction energy, imprints a periodic modulation on the intra-cavity energy, seen in the region marked by *ω*_{M} = 2FSR in Fig. 8c. The repetition rate shift can be offset by adjusting the modulation frequency: if the soliton spectral center moves from the pumped mode *j*_{0} to *j*_{0} + *η*_{r}, the repetition rate varies by Δ*D*_{1} ≈ *η*_{r}*D*_{2} and modifying the modulation frequency accordingly will balance soliton recoil, as illustrated in Fig. 8d–f.

In the state control approach resulting in Figs. 6–8, we have interjected intervals of no modulation between those with different modulation frequencies because the formulation presented herein assumes integer *M*. Although this step can be implemented experimentally, it is nonetheless not necessary and a continuous sweep of the modulation frequency is expected to produce the same effect. Adiabatic increase of the modulation frequency while preserving the number of solitons (e.g., one soliton per round-trip time) can transition the system through various DTC phases (e.g., successive single-soliton *n*-DTCs). Yet, as this procedure constitutes new TTSB states while altering the symmetry of the underlying system, it does not accurately parallel phase transitions in the context of condensed matter crystals where new phases under the same symmetry are studied.

It has been shown that the number of solitons in a CW-pumped Kerr cavity can be controlled by the forward or backward sweep of the pump-resonance detuning^{29,46}. Through the prudent choice of the sweeping rate, direction, and final detuning, this procedure can alter the number of solitons per round-trip time and in principle guarantee achieving any number of them, thereby enabling DTC state control. Moreover, driving the microresonator by two pulse trains^{27} with judiciously chosen repetition rates and relative delays can effectuate switching between different DTC states when one pulse train locks a subset of the microcomb solitons while the other tweezes the rest through changing the delay^{24}. In both cases of detuning sweep or driving by delayed pulse trains, the modulation can adiabatically be turned off prior to state switching and then back on (as in Figs. 6–8).

### Phase noise reduction of the RF beatnote

Analysis of DTC stability in the current system naturally connects to the notion of subharmonic phase noise. Therefore, we consider here the phase noise of the generated subharmonic and how it relates to that of the oscillator driving the modulator and a Kerr microcomb driven by a free-running laser.

Oscillator phase noise reduction by frequency division is frequently used to generate RF and microwave signals of high spectral purity from an optical frequency comb^{47,48}. The comb should first be stabilized by optical means, e.g., two-point locking or *f*−2*f* self-referencing accompanied with single-point locking to an external reference. The phase noise of the lasers utilized for comb stabilization is often much higher than standard high-quality RF sources. In contrast, the phase noise of the photonic RF signal demodulated on a fast photodiode can be better than high-quality RF sources. This improvement is achieved through frequency division.

Generally, when a frequency *f*_{0} = 〈*f*_{0}〉 + Δ*f*_{0} characterized with an average value 〈*f*_{0}〉 and random fluctuations Δ*f*_{0}(*t*) is divided by *n*, fluctuations of the smaller frequency *f*_{RF} = *f*_{0}/*n* is also reduced by the same factor to Δ*f*_{RF}(*t*) = Δ*f*_{0}(*t*)/*n*. It can readily be shown that the power spectral density of the two signals will be related through \({S}_{{{\Delta }}{f}_{{{{{{{{\rm{RF}}}}}}}}}}={S}_{{{\Delta }}{f}_{0}}/{n}^{2}\). Therefore, the phase noise \({{{{{{{{\mathcal{L}}}}}}}}}_{{\phi }_{{{{{{{{\rm{RF}}}}}}}}}}={S}_{{{\Delta }}{f}_{{{{{{{{\rm{RF}}}}}}}}}}/(2{f}^{2})\) of the divided frequency will also drop as *n*^{2}. In this expression *f* is the offset from the average frequency and the power spectral densities \({S}_{{{\Delta }}{f}_{0}}\) and \({S}_{{{\Delta }}{f}_{{{{{{{{\rm{RF}}}}}}}}}}\) indicate the random fluctuations of frequency rather than phase^{49,50}.

In view of the crucial role of subharmonic generation in DTC formation, TTSB in the platform discussed herein is naturally accompanied with phase noise reduction. A prominent feature of optical DTCs predicated on a pump-modulated or dichromatically pumped Kerr microcomb is indeed the possibility of decreasing the phase noise of the RF beatnote in comparison with the drive phase noise. We have recently reported this effect in DTCs in a dually pumped Kerr microcomb^{12,51}, where the subharmonic phase noise is reduced compared with the beatnote of the two driving lasers. In pump-modulated Kerr microcomb DTCs, subharmonic phase noise is reduced with respect to the signal modulating the laser. In both cases, the reduction factor is *n*^{2}, the DTC size squared.

In addition to the phase noise reduction due to frequency division, formation of DTCs in a parametrically seeded microcomb results in lowering the fundamental phase noise associated with Kerr frequency comb generation. Solitons in a Kerr comb pumped by a free-running laser experience random walk and the associated phase noise of the RF signal displays an *f*^{−2} dependence^{50}—the noise increases without limit at small frequency offsets, hinting at the unbound soliton timing jitter in the cavity^{52}. Subdued by modulation, this jitter in soliton arrival time can be constrained in a pump-modulated comb. The phase diffusion is suppressed and the phase noise gains an \({({f}^{2}+{\kappa }_{{{{{{{{\rm{Loop}}}}}}}}}^{2})}^{-1}\) dependence instead, where the nonzero relaxation factor *κ*_{Loop} is the locking loop bandwidth, which depends on the parameters of the modulated light, duration of the pulse, and finesse of the optical cavity. Owing to this coefficient, the DTC subharmonic phase noise can be smaller than the fundamental quantum noise of a free-running Kerr comb photonic oscillator at small offsets.

To better understand this phenomenon, we introduce the following expression, based on the theory of injection-locked oscillators and phase-locked loops^{53,54}. The phase noise of the slave oscillator (here, the DTC) improves when using a pure master oscillator (modulator) because

$${{{{{{{{\mathcal{L}}}}}}}}}_{{\phi }_{{{{{{{{\rm{RF}}}}}}}}}}=\frac{{f}^{2}}{{f}^{2}+{\kappa }_{{{{{{{{\rm{Loop}}}}}}}}}^{2}}{{{{{{{{\mathcal{L}}}}}}}}}_{{{{{{{{\rm{Kerr}}}}}}}}}+\frac{{\kappa }_{{{{{{{{\rm{Loop}}}}}}}}}^{2}}{{f}^{2}+{\kappa }_{{{{{{{{\rm{Loop}}}}}}}}}^{2}}\frac{{{{{{{{{\mathcal{L}}}}}}}}}_{{{{{{{{\rm{Mod}}}}}}}}}}{{n}^{2}},$$

(3)

where \({{{{{{{{\mathcal{L}}}}}}}}}_{{{{{{{{\rm{Kerr}}}}}}}}}\) is the phase noise of the comb-based oscillator, and \({{{{{{{{\mathcal{L}}}}}}}}}_{{{{{{{{\rm{Mod}}}}}}}}}\) is that of the oscillator driving the modulator. At offset frequencies much smaller than the loop bandwidth, the pump phase noise divided by *n*^{2} dictates the stability, while at larger offset frequencies the noise of the free-running Kerr frequency comb dominates. Therefore, the resultant phase noise can in some spectral bands be smaller than the fundamental noise of the Kerr comb. A detailed analysis of this phenomenon is beyond the scope of this work.