On the robustness of phase locking in Kerr optical frequency combs

We theoretically investigate the phase locking phenomena between the spectral components of Kerr optical frequency combs in the dynamical regime of Turing patterns.We show that these Turing patterns display a particularly strong and robust phase locking, originating from a cascade of phase locked triplets which asymptotically lead to a global phase locking between the modes. The local and global phase locking relationships defining the shape of the comb are analytically determined. Our analysis also shows that solitons display a much weaker phase locking that can be destroyed more easily than in the Turing pattern regime. Our results indicate that Turing patterns are generally the most suitable for applications requiring the highest stability. Experimental generation of such combs is also discussed in detail, and is in excellent agreement with the numerical simulations. © 2014 Optical Society of America OCIS codes: (140.3945) Microcavities; (140.4050) Mode-locked lasers; (190.4380) Nonlinear optics, four-wave mixing; (190.5530) Pulse propagation and temporal solitons. http://dx.doi.org/10.1364/OL.39.001529

Kerr optical frequency combs have been the focus of a strong interest from the scientific community since their first generation a few years ago [1].These combs result from a four-wave mixing process where new frequencies are generated through the nonlinear interaction between the light and the bulk material of a monolithic whispering gallery mode (WGM) resonator.Depending on the frequency and power of the pump laser, this phenomenology leads to the excitation of equally spaced spectral lines, thereby forming the so-called Kerr comb.This multiwavelength optical source is very promising for various applications such as high resolution spectroscopy and ultrastable microwave generation.In the latter case, the generated microwave corresponds to the beating between the different spectral modes, and its phase noise is directly related to the phase relationship between the comb's lines [2][3][4].The understanding of the phase locking process and its robustness with regard to external perturbations is therefore of the highest importance for ultrastable microwave generation, and the objective of this Letter is to unveil the mechanisms ensuring a stable and robust phase locking between the modes.
The experimental setup used for Kerr comb generation is shown on Fig. 1.The output of a continuous-wave 1550 nm laser is amplified and coupled to an ultrahigh Q WGM resonator using a tapered fiber.Above a given pump power threshold, a Kerr comb is generated and can be monitored using an optical spectrum analyzer.
A spatiotemporal formalism based on the Lugiato-Lefever equation (LLE) has recently been proposed to describe Kerr comb generation in WGM resonators [5][6][7].In our case, the LLE can be written as [6] where ψθ; τ is the complex amplitude of the overall field in the rotating frame of reference of the cavity, θ ∈ −π; π is the azimuthal angle along the cavity circumference, τ t∕2τ ph is the normalized time, with τ ph 1∕Δω tot being the photon lifetime, Δω tot ω 0 ∕Q tot the loaded linewidth, and Q tot the loaded quality factor.The parameter α −2Ω − ω 0 ∕Δω is the detuning between the pump laser and the cold resonance frequencies Ω and ω 0 , respectively.The parameter β represents the second-order dispersion of the resonator, while F stands for the amplitude of the pump laser field.
Various nontrivial solutions can emerge from the LLE, depending on the parameters α, β and F [8,9].In the regime of anomalous dispersion, the time-independent stationary solutions are either Turing patterns or bright solitons.Turing patterns correspond to an integer number of intensity extrema along the azimuthal direction of the cavity.In the Fourier domain, they correspond to the well-known Kerr combs with multiple free spectral range (FSR) spacing [10][11][12][13][14].
A few examples of such combs are presented in Fig. 2. They have been generated in a magnesium fluoride (MgF 2 ) resonator with intrinsic Q-factor equal to 1.5 × 10 9 , and a FSR equal to Ω FSR ∕2π 5.86 GHz.By choosing different resonances, pump powers, and frequencies we were able to obtain different mode spacings, each corresponding to a different integer number of FSR [11].For each one of these spectra, numerical simulations were performed using the split-step Fourier method on the LLE, and a good agreement with the experiments was found, hence validating the theoretical description.In the case of Fig.To investigate theoretically the phase locking of Turing patterns and their robustness relatively to solitons, an intermediate set point was chosen, with β −0.03, α 0, F 2 2.5, and noisy initial conditions, leading to the generation of the 12th order Turing pattern shown in Fig. 3(a).Contrary to Turing patterns that entirely fill the cavity, the bright cavity soliton is a unique structure spinning in the cavity [14].The spectrum of the soliton is therefore composed of spectral lines separated by only one FSR, as shown on Fig. 3 Turing patterns and cavity solitons are the two most promising Kerr comb regimes for ultrapure microwave generation.They are, indeed, structurally stable solutions of the LLE and can be obtained with parameter values far from chaotic regimes.In this application, the microwave signal corresponds to the beating between the different modes of the comb detected by a fast photodiode.To evaluate the phase properties of this beat note one needs to investigate the relative phases of the different modal components of the comb.From an analytical point of view this task can hardly be performed using the LLE.Instead, the following modal decomposition, allows direct access to the phase, φτ, of the slowly varying modal amplitudes, Ψ l jΨ l je iφ , with l l − l 0 being the azimuthal eigennumber of the photons with respect to the pumped mode (the pumped mode is therefore l 0, while the sidemodes correspond to l 1; 2; …).By plugging the expansion of Eq. ( 2) inside the LLE, the following ordinary differential equations ruling the dynamics of each mode are obtained: where δx is the usual Kronecker function, while m, n, p, and l are eigennumbers labeling the interacting modes following the interaction ℏω m ℏω p ↔ ℏω n ℏω l .Note that strict equivalence between the modal and spatiotemporal models had been demonstrated [6], so that the theoretical analysis can be performed with either formalism depending on their specific advantages.
The regime of the Turing pattern is essentially originating from modulational instability, which was analyzed in detail in Refs.[12][13][14][15][16][17].Turing patterns actually arise from the following cascade of interactions.When the pump is set above the threshold F 2 th 1 1 − α 2 , the modes around l th , with l th 2α − 2∕β 1∕2 , experience positive gain and are thereby excited.After transient dynamics characterized by modal competition, two symmetric sidemodes, L ≃ l th , are stabilized around the pump and the net flux of photons follows the degenerate interaction 2ℏω 0 → ℏω L ℏω −L .From then, the sidemodes L provide, in their turn, a net flux of photons to the sidemodes 2L using interactions ℏω 0 ℏω L → ℏω ∓L ℏω 2L .This cascaded process is further generalized later through the interactions ℏω 0 ℏω L → ℏω ∓kL ℏω k1L where the sidemodes kL provide a net flux of photons to the sidemodes k 1L.This cascaded fountain effect asymptotically yields a Kerr comb characterized by a multiple-FSR spacing of multiplicity L. The effect of this cascading mechanism can be observed in Fig. 4, displaying the phase dynamics of the 12th order Turing patterns of Fig. 3(a).It can be clearly seen that,  after transient dynamics, the various modal phases, ϕ l φ l − φ 0 , with respect to the pump are reaching a time-independent value.The phase stabilization occurs first for the mode l 12, and then to l 24, l 36, and so on.
It is important to note that the phase locking value depends on the initial condition, as illustrated on Fig. 5(a), where the phase difference, ϕ l φ l − φ 0 for l 12, has been plotted for different initial conditions.This phenomenon is true for all of the excited modes of the Turing pattern.However, it can be shown that an intrinsic phase locking relationship takes place between the modes of the comb.According to the scenario described earlier, a given mode, kL, mainly interacts with its two neighboring modes, k 1L.In the steady state (dΨ kL ∕dτ ≡ 0), this three-mode interaction can be rewritten as [15] κ kL Ψ kL δkF jΨ kL j 2 Ψ kL with κ l −1 iα iβl 2 ∕2.This approximate relationship, obtained from a photon flux balance standpoint, can be used to show that consecutive triplets of modes are necessarily phase locked according to Hence, the phases, φ kL , might individually have an asymptotic value that depends on the initial conditions but, for the same Turing pattern, the combination φ k−1L − 2φ kL φ k1L will be independent of initial conditions and will converge to the same value, ξ k , as can be seen in Fig. 5(b).The value of ξ k depends only on the parameters α, β, F, and k.The sequence is, therefore, the following.Initially, the triplet f−L; 0; Lg is the first to phase lock.Then, the two triplets f0; L; 2Lg phase lock and, later on, all the cascade of triplets fk − 1L; kL; k 1Lg as k → ∞ successively phase lock as well.This sequence of triplet phase locking ultimately leads to a global phase coupling through a cascaded process, and the global phase locking relationship can be obtained by considering that since any ξ k is constant so is any linear combination of them.If we consider 2N 1 oscillating modes in a Turing pattern comb (the pump and N excited modes on each side), we can show that If we consider, for example, the 7 oscillating modes that are visible in Fig. 3 .Hence, when the permanent state is reached, the linear combination, Ξ N , which involves all the modal phases is a constant that is independent of the initial conditions, as evidenced by Fig. 5(c).These intrinsic phase relations, ξ k and Ξ N , are, therefore, responsible for the shape of the Kerr comb in the Turing pattern regime.
In order to test the robustness of the phase locking in both the Turing pattern and the soliton regimes, we introduce a perturbation after the stationary state is reached.The phase of the pump signal is abruptly shifted by constant offset, and we monitor the difference, Δϕ l t ϕ l t − ϕ st l , with ϕ l φ l − φ 0 .When the offset in the pump phase is set to π∕5, the relative phases of the different modes oscillate in response to this perturbation, but the phase locking is preserved.In the case of the Turing patterns, the phase oscillations are of small amplitude with respect to the perturbation and decrease when the mode number increases, as if the perturbation was attenuated during its propagation to higher order modes.In contrast, the amplitudes of the phase oscillations of the soliton are larger than the initial stimulus and do not substantially decrease for higher order modes.This result indicates that Turing patterns are more robust to external phase perturbations and could lead to better phase noise performances.To strengthen this claim, the phase shift of the pump signal was increased to π∕2 and the resulting phase evolution is plotted in Figs.6(c) and 6(d).The Turing pattern behavior remains very stable and identically phase locked, while the soliton is destroyed (as well as the phase locking).
In this work, we have given an insight into the dynamics of phase locking in a Kerr frequency comb.The fact that the photonic interactions are strongly dominated by the exchange of microwave photons of energy L × ℏΩ FSR is a direct consequence of the modulational instability, which generates the highest gain for photons whose offset frequency relative to the pump is precisely L × Ω FSR .We have analytically demonstrated that phase locking in Turing patterns is achieved through a cascaded phase locking of adjacent triplets, which ultimately leads to a global phase locking relationship that we have determined.We have also shown that Turing patterns are more robust than solitons with regard to external perturbations, thereby demonstrating their stronger potential for ultrastable microwave generation.However, further investigations are needed in order to generate Turing patterns with frequency spacings below 50 GHz, using for example WGMs with larger anomalous dispersion.Y. C. K. acknowledges financial support from the European Research Council (ERC) through the project NextPhase.