Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Stability analysis of mode-coupling-assisted microcombs in normal dispersion

Open Access Open Access

Abstract

We theoretically study the stability of mode-coupling-assisted frequency comb generation in normal-dispersion microresonators. With the aid of mode coupling, quantitative analysis of the modulational instability is explored in the parameter space of pump power and detuning. By exploring the coupled mode number, dispersion, and coupling strength in the normalized Lugiato-Lefever model, the modulational stability gain exists and yields extended spatial structures within the regime of eigenvalue bifurcations. Moreover, the dynamics and efficiency of microcombs are discussed, providing the accessibility of high-efficient, stable, and controllable combs. This work offers universal guidelines for operating mode-coupling-assisted combs in a normal-dispersion system.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Kerr frequency comb generation consisting of a serial of equidistant frequency components in high quality factor (Q) microresonators has been widely studied within the past decades [16]. This phenomenon is initiated by pumping the cavity mode with a continuous-wave (CW) laser and transferring the photon energy to different cavity modes through the cascaded four-wave-mixing (FWM) process. Models for Kerr comb generation have been developed from coupled mode equations with limited mode numbers [7,8] to the spatiotemporal formulation known as Lugiato-Lefever equation (LLE), deriving from a driven, detuned and damped nonlinear Schrödinger equation [911]. Early theoretical models predict microcombs in the anomalous-dispersion regime through the modulation instability (MI) in microresonators [8,12] while hard excitation offers the accessibility of microcomb generation in both normal- and anomalous- dispersion regimes [1315]. Although Turing patterns and dark cavity solitons are possible nontrivial solutions in the normal-dispersion regime, it requires either large frequency detunings (from cavity resonance) or a suitable initial condition of cavity fields [16]. Lately, comb generation is also identified in normal dispersion with soft excitation [1720]. This can be theoretically explained by the local dispersion change modified by backward-scattering, multiple spheroidal modes, or waveguide modes [21,22]. Since different mode families have morphologically dependent group index, the free-spectral range (FSR) and waveguide dispersion are significantly different, limiting comb excitation only from one mode family [21]. As two mode families approach each other around a specific frequency, mutual coupling induces significant frequency shift in the resonant frequencies, explaining by the mode-coupling model [23], and then the shift results in the avoided mode crossings (AMX) [5,19,21]. This effect not only changes the comb spectra locally [24,25] but also paves the way for microcomb initiation [5,25] in normal-dispersion microresonators. In addition, by controlling the position of AMX, this mechanism dominates the microcomb generation in the normal-dispersion regime, in which the initial comb lines (or primary combs) can be selectively generated in the form of Turing rolls or dark solitons [26]. Even with weak mode coupling instead of observing significant AMX, microcomb can be still generated in the large normal-dispersion regime [27]. Theoretically, by utilizing the coupled mode equations and LLE, model shows good agreement with the experimental results [5,28]. However, although these reports based on the LLE provide significant insight into the generation dynamics for Kerr microcombs in normal dispersion, the discussion on stability regime of threshold power and comb dynamics is still limited. Meanwhile, as for the conversion efficiency, previous literatures have analytically shown that the comb efficiency in anomalous dispersion can be > 40% as the pattern of Turing rolls [29]; while in normal dispersion, it was shown that a high-efficient, dark pulse is generated [30]. However, the conversion efficiency is mainly studied in the dark-soliton regime requiring high-power operation or a nontrivial initial condition of intracavity fields. There are limited works on comb efficiency with mode coupling, especially in the low-power, controllable Turing-roll regimes [26,28]. It is now essential to address these issues for practical operation in normal dispersion.

In this letter, we theoretically investigate the relation between mode coupling, MI gain, and conversion efficiency of combs by solving the normalized mean-field LLE. Similar to the previous identification in anomalous dispersion [31], the stability regime of coupling-mode-assisted microcomb generation in normal dispersion is explored in the space of detunings and pump powers. This work delivers several new findings. First, the stability regime of microcomb generation in normal dispersion is shown by introducing the mode coupling to the normalized LLE. We observe that the MI gain and stability regime are strongly correlated to the coupling strength between different modes. In addition, unlike that in anomalous dispersion, an upper bound on power for MI gain is identified. Second, we study the accessibility of microcomb. Through the normalized LLE, the dynamics of the mode-coupling-assisted combs is explored in the similar parameter space of detuning and pump power [16,29,31,32]. Due to the cubic dependency of CW steady-state solutions in a nonlinear cavity, hysteresis behavior is prominently identified by shifting the detuning in opposite directions. In addition, the threshold boundary is clarified with mode coupling. Last, the conversion efficiency of mode-coupling-assisted combs is discussed in the normal-dispersion regime. Stable, high efficient (≈ 50%) Turing rolls are observed with low operating power (down to a few tens milliwatt). This analytical work allows us to have more insights into practical experiments of low-power, stable, and controllable microcomb generation.

2. Theory

The generalized mean-field LLE can be written as [29,31,32]:

$${T_R}\frac{{\partial E(t,\theta )}}{{\partial t}} = [ - (\alpha + i{\delta _0}) - iL\frac{{{\beta _2}}}{2}{(\frac{{2\pi }}{{{T_R}}})^2}\frac{{{\partial ^2}}}{{\partial {\theta ^2}}} + i\gamma L|E{|^2}]E + \sqrt \kappa {E_{in}}$$
where TR is the roundtrip time, α is the loss per roundtrip of the cavity, δ0 is the pump detuning, respectively. β2 is the 2nd-order dispersion coefficient of Taylor series expansion of propagation constant at the pump frequency while the higher order dispersion is ignored, γ is the nonlinearity coefficient, and L is the circumference of the resonator. t is the slow time describing the evolution of the intracavity field E (normalized such that the cavity power Pcavity=|E|2) while θ is the azimuthal coordinate describing the wave in the resonator. Ein is the input pump field with the coupling coefficient ${\boldsymbol \kappa }$. The normalized LLE is then given by setting $A({t^{\prime},\; \theta } )= E({t,\; \theta } )\sqrt {\gamma L/\alpha } $, $t^{\prime} = \alpha t/{T_R}$, $\Delta = {\delta _0}/\alpha $, $S = {E_{in}}\sqrt {\gamma L\kappa /{\alpha ^3}} $, and $\beta = 4{\pi ^2}L{\beta _2}/\alpha T_R^2$:
$$\frac{{\partial A(t^{\prime},\theta )}}{{\partial t^{\prime}}} = [ - (1 + i\Delta ) - i\frac{\beta }{2}\frac{{{\partial ^2}}}{{\partial {\theta ^2}}} + i|A{|^2}]A + S$$

The normalized field envelope of a single mode family can be then expressed as

$$A({t^{\prime},\; \theta } )= \mathop \sum \nolimits_l {A_l}({t^{\prime}} ){e^{il\textrm{FSR}\frac{\theta }{{2\pi }}{T_R}}} = \mathop \sum \nolimits_l {A_l}({t^{\prime}} ){e^{il\theta }}$$
where ${A_l}({t^{\prime}} )$ is the complex field envelopes for individual mode number l of microcombs with respect to the pump line (l = 0). By inserting this expansion into Eq. (2), we can show the dynamics of complex field envelopes of mode number l [16]:
$$\frac{{\partial {A_l}(t^{\prime})}}{{\partial t^{\prime}}} = [ - (1 + i\Delta ) + i\frac{\beta }{2}{l^2}]{A_l} + \delta (l)S + i\sum\limits_{m,n,p} {\delta (m - n + p - l){A_m}} A_n^\ast {A_p}$$
where $\delta (x)$ is the Kronecker function. Next, we consider the mode coupling by utilizing additional phase shift Δshift from the resonance shift of mode coupling [5].
$$\frac{{\partial {A_l}(t^{\prime})}}{{\partial t^{\prime}}} = [ - (1 + i(\Delta + {\Delta _{shift}})) + i\frac{\beta }{2}{l^2}]{A_l} + \delta (l)S + i\sum\limits_{m,n,p} {\delta (m - n + p - l){A_m}} A_n^\ast {A_p}$$

We should note this shift will be only applied to the field envelope (Al) of coupled mode l while other modes are remained unchanged as in Eq. (4). By introducing small perturbation $\delta {A_l}({t^{\prime}} )$ in ${A_l}({t^{\prime}} )$, the matrix of the linearized system can be then expressed as [16,31]:

$$\left[ {\begin{array}{{c}} {\frac{{d\delta {A_l}({t^{\prime}} )}}{{dt^{\prime}}}}\\ {\frac{{d\delta {A_{ - l}}^\ast ({t^{\prime}} )}}{{dt^{\prime}}}} \end{array}} \right] = \left[ {\begin{array}{{cc}} {{M_l}}&N\\ {{N^\ast }}&{{M_{ - l}}^\ast } \end{array}} \right]\left[ {\begin{array}{{c}} {\delta {A_l}({t^{\prime}} )}\\ {\delta {A_{ - l}}^\ast ({t^{\prime}} )} \end{array}} \right]$$
where ${M_l} ={-} ({1 + i({\mathrm{\Delta } + {\mathrm{\Delta }_{shift}}} )+ 2\textrm{i}{{|{{A_0}} |}^2} + i\beta {l^2}/2} )$, ${M_{ - l}} ={-} ({1 + i\mathrm{\Delta }} )+ 2\textrm{i}{|{{A_0}} |^2} + i\beta {l^2}/2)$, and $N = iA_0^2$. An additional phase shift is applied to mode l in the presence of mode coupling while ${M_{ - l}}$ for the corresponding mode –l remains unchanged. It should be noted that here we explore the comb initiation of mode l by introducing small perturbation and other modes are negligible in comparing to the pump (${A_0}$). The eigenvalues λ of the above linearized system can be written as:
$$\lambda (l )={-} 1 - i\frac{{{\Delta _{shift}}}}{2} \pm \sqrt { - {{\left( {\Delta - 2{{|{{A_0}} |}^2} + \frac{{{\Delta _{shift}}}}{2} - \frac{{\beta {l^2}}}{2}} \right)}^2} + {{|{{A_0}} |}^4}} $$
When Re(λ) > 0, it results in an exponentially growing perturbation solution, similar to the modulational instability identified in anomalous dispersion. In this paper, we will focus on eigenvalue solutions in the normal-dispersion regime (β > 0). Traditionally in normal dispersion without considering mode coupling and detunings, Re(λ) is always negative and no MI gain can be obtained. With the aid of mode coupling, positive Re(λ) can exist and achieve nontrivial solutions besides equilibria (flat solutions). We will then study the stability boundary in the parameter space of S (in related to the cavity pump field ${A_0}$) and $\Delta $ by exploring different dispersion parameters β, mode numbers l, and strength of mode coupling ${\Delta _{shift}}$. The simulation parameters are set with intrinsic Q = 3 million (α= 0.0017 for critical coupling), γ = 1.09 W−1·m−1, and L = 2π·(100 µm) with group index ng = 2, unless mentioned otherwise. The FSR for the resonator is around 239 GHz. These parameters are based on the mature SiN-based platform [32].

3. Gain spectrum and stability analysis of mode-coupling-assisted Kerr microcombs

3.1 Mapping of MI gain with mode coupling

By solving the eigenvalues from Eq. (7), we show the maps of Re(λ) and the corresponding power spectra for different coupling strength in Fig. 1. The parameter space is in the range of mode number 0 < l < 15 and pump power 0 < S2 < 50 (0 < |Ein|2 < 210 mW). The pump detuning is set at $\Delta $ = 0, assuming pumping at the resonance of the cold cavity. First, for normal dispersion without mode coupling (${\Delta _{shift}} = 0$) as shown in Fig. 1(a), Re(λ) is always negative. This means that small perturbation $\delta {A_{ {\pm} l}}({t^{\prime}} )$ does not grow and the intracavity power is stable as a CW. When the mode coupling (${\Delta _{shift}}$) is applied, positive Re(λ) exists as shown in Fig. 1(b)–(c). The MI gain results in stable Turing patterns and the corresponding spectra in the cavity are shown in Fig. 1(e). This identification is in consistent with the previous experimental findings [5,26,28]. The central idea of the periodic Turing rolls is that they emerge from noise while a flat solution becomes unstable. One interesting behavior is that, although the mode coupling (${\Delta _{shift}}$) is applied to a specific mode l and MI gain is possible, the gain spectrum has a gain cutoff lmax. Both the available gain and gain bandwidth grow with the strength of mode coupling in analogy with the case of anomalous dispersion [33]. For instance, lmax increases from 7 to 10 by raising the strength ${\Delta _{shift}}$ from 6 to 8. This effect can be also identified from the relative microcomb power of the generated spectra in Fig. 1(e) and Fig. 1(f). It suggests stronger coupling strength is needed for the mode number farer away from the pump line.

 figure: Fig. 1.

Fig. 1. Mapping of MI gain versus the mode number and pump power for (a) ${\Delta _{shift}}$ = 0, (b) ${\Delta _{shift}}$ = 6, and (c) ${\Delta _{shift}}$ = 8. The corresponding spectra at S2 = 10 and l = 5 are shown in (d)-(f), respectively.

Download Full Size | PDF

Next, for mode coupling at a particular mode, we see that the gain spectra only exist within a certain range of pump power. Unlike that a stronger pump power is favored in anomalous dispersion, providing more gain and broader bandwidth for microcomb generation, an upper boundary of pump power is now identified in normal dispersion. It implies that the pump power needs not only to satisfy the threshold condition of cavity loss [29] but also to work within the boundary of MI gain. For instance, the exemplary spectra are shown in Fig. 2 at S2 = 2, 10, and 50 by setting Δ= 0, Δshift= 8, and l = 5. Microcomb in the form of Turing-roll patterns can be initiated from the random noise fluctuation for S2 = 10, but no oscillation can be observed for the pump below threshold (e.g. S2 = 2) or above the upper boundary (e.g. S2 = 50). This poses a limitation on the required pump for microcomb formation in the normal-dispersion regime. Also, as stated above, the mode spacing is 5 FSRs, determining by the position of mode coupling.

 figure: Fig. 2.

Fig. 2. Spectra for different pump powers S2 = 2, 10, and 50. The microcomb is initiated from the random noise at Δ = 0, ${\Delta _{shift}}$ = 8, l = 5, and S2 = 10 but no oscillation can be found at S2 = 2 and 50.

Download Full Size | PDF

In Fig. 3, we compare the MI gain without and with mode coupling in the parameter space of pump power and detuning. As shown in Fig. 3(a)-(c), for the gain spectrum without mode coupling (${\Delta _{shift}}$ = 0), microcombs can only be initiated with a large detuning and pump power (e.g. Δ= 3 and S2 = 6) [1316], suggesting a much narrower gain regime in related to that in anomalous cases. An example of the microcomb spectrum is shown in Fig. 3(d), evidencing a periodic stable solution (Turing rolls). As for the one with mode coupling (${\Delta _{shift}}$ = 8, Fig. 3(c)), we can observe the positive MI gain is now accessible with a small detuning. Again, similar Turing-roll patterns can be identified. The corresponding spectrum is shown in Fig. 3(e). One significant difference is that the periodicity of microcombs (or the position of primary combs) is determined by the coupled mode number. This is in consistent with the previous findings [26,28] due to the fact that positive gain is only available from the mode coupling. It should be noted that, for microcomb generation without mode coupling, the oscillation induced by MI gain may come close to the unstable branch of the intracavity power (in-between the two stable equilibria, black-dashed lines). Therefore, unless well-split stable branches of pump with a relative large detuning are fulfilled, the oscillation may exhibit unstable Turing patterns and could be repelled to other stable equilibrium (flat solutions) [16]. As for the cases with mode coupling, MI gain can be generated very far away from the unstable point, providing a stable Turing-roll patterns.

 figure: Fig. 3.

Fig. 3. MI gain spectra without (${\Delta _{shift}}$ = 0) for (a) l = 5 and (b) l = 8 and (c) with (${\Delta _{shift}}$ = 8) mode coupling for l = 5. The corresponding spectra (d) without and (e) with mode coupling.

Download Full Size | PDF

3.2 Stability regime analysis

To further study the effect on mode coupling, the maps of eigenvalue bifurcations are shown in Fig. 4, again in $\Delta $-S2 space. By introducing mode coupling and solving Eq. (7), we show the maps at different effective resonance shift $y^{\prime} = \; \frac{{{\Delta _{shift}}}}{2} - \frac{{\beta {l^2}}}{2}$, normalized dispersion parameters $\beta $, mode numbers l, and coupling strength ${\Delta _{shift}}$. We aim to find the MI gain regime in the presence of mode coupling. The curves in Eq. (7) suggest Re(λ) = 0 solutions in $\Delta $-S2 space. Since both mode coupling and dispersion shift the resonant frequencies, we first analyze the effect on effective resonance shift ($y^{\prime}$) in Fig. 4(a). Similar to that observed in Ref. [16], the normalized LLE yields a cubic equation in which multiple states emerge as the detuning is increasing. For instance, if we consider the effective resonance shift $y^{\prime} = 0$, only single solution of the eigenvalue exist and there is no net gain for $\Delta < \sqrt 3 $. For $\Delta > \sqrt 3 $, two boundary lines, corresponding to two critical pump powers, can be identified. However, still no gain is found and the cavity field remains as a stable CW. As increasing the effective resonance shift (e.g. $y^{\prime} = 4$), bifurcation now emerges at a smaller detuning. Meanwhile, for the eigenvalue solutions in-between the bifurcation, the intermediate Re(λ) is positive and therefore results in MI gain for microcomb generation (unstable perturbation $\delta {A_l}$).

 figure: Fig. 4.

Fig. 4. Maps of eigenvalue bifurcations in normal dispersion in related to different (a) effective coupling strength y’, (b) normalized dispersion $\beta $, (c) mode number l, and (d) coupling strength ${\Delta _{shift}}$.

Download Full Size | PDF

Next, we study the effect on the magnitude of dispersion parameter β with fixed coupling strength (${\Delta _{shift}} = 6$) at l = 5 in Fig. 4(b). Significantly, as normal dispersion gets stronger, this results in the increase of critical detuning (${\Delta _{crit}}$) to achieve multiple solutions in the bifurcation maps, narrowing down the MI gain regime. For example, ${\Delta _{crit}}$ moves from −0.82 to 0.24 for β increasing from 0.042 to 0.126, suggesting that stronger detuning is needed to compensate the deep normal dispersion. In contrast to Fig. 4(a), the right boundaries are almost identical to different dispersion values. This can be explained by that the solution on the right boundary is dependent on $\Delta + \frac{{{\Delta _{shift}}}}{2} - \frac{\beta }{2}{l^2}$. Since $\Delta + \frac{{{\Delta _{shift}}}}{2}$ is order of magnitude stronger than $\frac{\beta }{2}{l^2}$ within the operation regime, the difference in the right boundary is negligible. As for the effect on the mode number l, similar trends can be identified as those in Fig. 4(c). Increasing the mode number effectively results in stronger normal dispersion. The quadratic dependency of mode number exhibits more significance on the power of the right boundary.

Lastly, we evaluate the effect on coupling strength. For weak mode coupling (e.g. ${\Delta _{shift}} = 3$), only a small intermediate regime ($\Delta $ from 0.73 to 3.57) of Re(λ) is positive. The area of MI gain is small; this regime expands as increasing the coupling strength. At the same time, stronger coupling contributes to lower ${\Delta _{crit}}$. For instance, ${\Delta _{crit}}$ shifts from 0.73 to −0.77 as doubling the strength ${\Delta _{shift}}$ from 3 to 6.

4. Accessibility, threshold power, generation dynamics, and efficiency of microcombs in normal dispersion

4.1 Comb accessibility and generation dynamics

The analysis of MI gain in Section 3 allows us to obtain the possible regimes for mode-coupling-assisted microcomb generation. However, beside the corresponding detuning and pump power in $\Delta $-S2 space, the normalized LLE in Eq. (2) with cubic nonlinearity exhibits the well-known hysteresis dynamics as adjusting either the detuning [16] or pump power [29]. The nontrivial equilibria are critical to the existence of microcombs in the bifurcation map [16]. We now study the accessibility of microcombs with mode coupling in the normal-dispersion regime. Figure 5 shows maps of comb efficiency in $\Delta $-S2 space for different modes l and coupling strength ${\Delta _{shift}}$, following the paths of red- (Fig. 5(a)–(d)) and blue-detuning (Fig. 5(e)–(h)). The dashed line shows the eigenvalue bifurcations as that in Figs. 4. Clearly, no matter which path we choose, microcomb is generally initiated as passing the MI gain regime, in the form of periodic stable solution (Turing rolls). This periodic wave remains stable into the regime where the CW is stable (no MI gain). This identification is in analogy to the soliton case in anomalous dispersion [31]. One significant difference is that, since the MI gain will only occur within the regime where mode coupling is introduced, the periodicity of Turing rolls is determined by the coupled mode number l for the entire path. With the determined periodicity by mode-coupling, the Turing patterns may later evolve into the soliton crystals (SCs) with a series of equidistant dissipative Kerr solitons (DKSs) while a stronger pump power and large detunings are studied [34,35]. We can also notice that, due to the higher available comb power on the right boundary of bifurcation curves, the hysteresis is more significant for red-detuning than blue-detuning. In addition, as stated in the previous section, strong mode coupling with mode number close to pump line (e.g. l = 1) leads to larger gain regime for comb generation, evidencing the trend in Fig. 5; as for weak mode coupling with large mode number (e.g. l = 5), the gain regime becomes relatively narrow.

 figure: Fig. 5.

Fig. 5. Maps of comb efficiency in $\Delta $-S2 space under red-detuning for (a) l = 1 and ${\Delta _{shift}}$ = 3, (b) l = 1 and ${\Delta _{shift}}$ = 6, (c) l = 5 and ${\Delta _{shift}}$ = 3, and (d) l = 5 and ${\Delta _{shift}}$ = 6, and (e)(h) the corresponding maps under blue-detuning.

Download Full Size | PDF

To further address the hysteresis dynamics, we discuss two examples in more detail for the blue-detuned cases. Although we have explicitly explored the MI gain regime with mode coupling, this scheme does not give the full explanation as operating the system, especially with large detunings. Figure 6(a) shows the blue-detuned case for l = 1 and ${\Delta _{shift}}$ = 6 while the white-dashed line indicates the eigenvalue bifurcation as that in Fig. 4, the red-dashed line shows the normalized intracavity pump ${|{{A_0}} |^2} = 1$, and the black-dashed line shows the two equilibria of the intracavity pump. The significant difference is that, when the pump is blue-detuned crossing the gain regime, the oscillation induced by MI is determined not by the right boundary of the eigenvalue bifurcations but close to the line of ${|{{A_0}} |^2} = 1$. This can be explained by that, although we show the gain regime with mode coupling, the eigenvalue solutions of pump (white-dashed line) in-between the two branches of pump equilibria (black-dashed lines) correspond to the unstable branch of the three equilibria. The negative slope of intracavity response is unstable with respect to CW perturbation [5]. Therefore, even with MI gain in-between the pump equilibria as shown in the zoom-in spectrum of Fig. 6(a), the intracavity pump may move to the stable equilibria and no stable microcomb can be identified. The oscillation regime is then determined by increasing the detuning and intracavity power as that in traditional normal-dispersion cavities [5]. Figure 6(b) shows another blue-detuned example for l = 5 and ${\Delta _{shift}}$ = 6. Since the mode number increases, it results in a smaller gain regime as explained in the previous section. We can see that again oscillation regime matches with the right side of the line of ${|{{A_0}} |^2} = 1$ when eigenvalue solution is in the unstable regime. We should note here that the tuning direction (blue-detuning) is opposite to the resonance red-shifting induced by the thermal effect. It generally results in unstable operation within this MI regime.

 figure: Fig. 6.

Fig. 6. Maps of comb efficiency in $\Delta $-S2 space under blue-detuning for (a) l = 1 and ${\Delta _{shift}}$ = 6 and (b) l = 5 and ${\Delta _{shift}}$ = 6. The corresponding zoom-in spectra are also shown on the right.

Download Full Size | PDF

4.2 Comb efficiency in mode-coupling-assisted Kerr microcombs

Last, we discuss the efficiency maps in normal dispersion with mode coupling. First, unlike the frequency-sensitive MI gain in anomalous dispersion [33], the periodicity of Turing rolls is determined by the mode coupling. Therefore, as shown in Fig. 5 and Fig. 6, this provides a smooth map instead of showing abrupt boundaries in efficiency which is induced by the evolution of mode number in the anomalous-dispersion case [29,31]. Second, the microcomb efficiency can be up to ≈ 25% for the regime of interests. By red-detuning the pump, higher efficiency is obtained. This can be explained by that, while the detuning is increasing and sweeping across the cavity resonance, higher cavity power and enhancement factor are available, as shown in Figs. 7. Since the phase is red-detuned, the enhancement factor is determined by the upper branch of pump. When the comb threshold is reached, pump starts to saturate and transfer additional power above threshold to the microcombs [29,36]. Therefore, it gives better efficiency for microcomb generation. Third, the efficiency is also strongly dependent to the coupling condition 529]. Although a critical coupling scheme is applied to the above discussion, higher comb efficiency is available with an over-coupled device [29]. As shown in Fig. 8, by tripling the coupling coefficient ($\kappa $ = 0.0051) for over-coupling, ≈ 50% efficiency is achieved with the input power down to a few tens milliwatt (e.g. S2 = 3).

 figure: Fig. 7.

Fig. 7. (a) Normalized intracavity power versus phase detuning at S2 = 2, 5, and 8. (b) The enhancement factor of the intracavity power when the phase detuning is continuously increasing (under red-detuning). The red-dashed and black-dashed lines represent the normalized intracavity pump ${|{{A_0}} |^2} = 1$ and the two equilibria of the intracavity pump, respectively.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. Maps of comb efficiency in $\Delta $-S2 space under red-detuning for (a) l = 1 and ${\Delta _{shift}}$ = 1.5.

Download Full Size | PDF

5. Conclusions

To conclude, by solving the normalized LLE, we theoretically study the MI gain, stability, accessibility, and dynamics for comb generation in normal-dispersion resonators with mode coupling. Considering the phase shift of mode coupling at a specific mode, parametric gain allows MI in the normal-dispersion resonators and shows strong correlation between the coupling strength, mode number, pump power, and detuning. We solve the bifurcation maps of MI gain to model the accessibility of microcomb generation. In comparison with the previous methods, our analytic work provides universal guidelines for operating microcomb in the normal-dispersion regime, especially with small detuning and low pump power. In addition, we investigate the comb dynamics under different detuning directions, showing good agreement with the bifurcation maps. With the aid of mode coupling, high efficient, stable microcombs can be generated in the form of Turing rolls with a controllable repetition rate.

Funding

National Science and Technology Council (NSTC), Taiwan (109-2221-E-008-091-MY2, 111-2221-E-008 -026).

Acknowledgment

We thank research financial support from the National Science and Technology Council, Taiwan.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. P. Del’Haye, A. Schliesser, O. Arcizet, T. Wilken, R. Holzwarth, and T. J. Kippenberg, “Optical frequency comb generation from a monolithic microresonator,” Nature 450(7173), 1214–1217 (2007). [CrossRef]  

2. J. S. Levy, A. Gondarenko, M. A. Foster, A. C. Turner-Foster, A. L. Gaeta, and M. Lipson, “CMOS-compatible multiple-wavelength oscillator for on-chip optical interconnects,” Nat. Photonics 4(1), 37–40 (2010). [CrossRef]  

3. L. Razzari, D. Duchesne, M. Ferrera, R. Morandotti, S. Chu, B. E. Little, and D. J. Moss, “CMOS-compatible integrated optical hyper-parametric oscillator,” Nat. Photonics 4(1), 41–45 (2010). [CrossRef]  

4. W. Liang, A. A. Savchenkov, A. B. Matsko, V. S. Ilchenko, D. Seidel, and L. Maleki, “Generation of near-infrared frequency combs from a MgF2 whispering gallery mode resonator,” Opt. Lett. 36(12), 2290–2292 (2011). [CrossRef]  

5. X. Xue, Y. Xuan, Y. Liu, P. H. Wang, S. Chen, J. Wang, D. E. Leaird, M. Qi, and A. M. Weiner, “Mode-locked dark pulse Kerr combs in normal-dispersion microresonators,” Nat. Photonics 9(9), 594–600 (2015). [CrossRef]  

6. A. L. Gaeta, M. Lipson, and T. J. Kippenberg, “Photonic-chip-based frequency combs,” Nat. Photonics 13(3), 158–169 (2019). [CrossRef]  

7. A. A. Savchenkov, A. B. Matsko, D. Strekalov, M. Mohageg, V. S. Ilchenko, and L. Maleki, “Low threshold optical oscillations in a whispering gallery mode CaF2 resonator,” Phys. Rev. Lett. 93(24), 243905 (2004). [CrossRef]  

8. Y. K. Chembo, D. V. Strekalov, and N. Yu, “Spectrum and dynamics of optical frequency combs generated with monolithic whispering gallery mode resonators,” Phys. Rev. Lett. 104(10), 103902 (2010). [CrossRef]  

9. T. Kippenberg, S. Spillane, and K. Vahala, “Kerr-nonlinearity optical parametric oscillation in an ultrahigh-Q toroid microcavity,” Phys. Rev. Lett. 93(8), 083904 (2004). [CrossRef]  

10. L. A. Lugiato, F. Prati, M. L. Gorodetsky, and T. J. Kippenberg, “From the Lugiato–Lefever equation to microresonator-based soliton Kerr frequency combs,” Phil. Trans. R. Soc. A. 376(2135), 20180113 (2018). [CrossRef]  

11. T. Herr, V. Brasch, J. D. Jost, C. Y. Wang, N. M. Kondratiev, M. L. Gorodetsky, and T. J. Kippenberg, “Temporal solitons in optical microresonators,” Nat. Photonics 8(2), 145–152 (2014). [CrossRef]  

12. I. H. Agha, Y. Okawachi, M. A. Foster, J. E. Sharping, and A. L. Gaeta, “Four-wave-mixing parametric oscillations in dispersion-compensated high-Q optical microspheres,” Phys. Rev. A 76(4), 043837 (2007). [CrossRef]  

13. A. B. Matsko, A. A. Savchenkov, and L. Maleki, “Normal group-velocity dispersion Kerr frequency comb,” Opt. Lett. 37(1), 43–45 (2012). [CrossRef]  

14. A. B. Matsko, A. A. Savchenkov, V. S. Ilchenko, D. Seidel, and L. Maleki, “Hard and soft excitation regimes of Kerr frequency combs,” Phys. Rev. A 85(2), 023830 (2012). [CrossRef]  

15. T. Hansson, D. Modotto, and S. Wabnitz, “Dynamics of the modulational instability in microresonator frequency combs,” Phys. Rev. A 88(2), 023819 (2013). [CrossRef]  

16. C. Godey, I. V. Balakireva, A. Coillet, and Y. K. Chembo, “Stability analysis of the spatiotemporal Lugiato-Lefever model for Kerr optical frequency combs in the anomalous and normal dispersion regimes,” Phys. Rev. A 89(6), 063814 (2014). [CrossRef]  

17. P. H. Wang, F. Ferdous, H. X. Miao, J. Wang, D. E. Leaird, K. Srinivasan, L. Chen, V. Aksyuk, and A. M. Weiner, “Observation of correlation between route to formation, coherence, noise, and communication performance of Kerr combs,” Opt. Express 20(28), 29284–29295 (2012). [CrossRef]  

18. A. A. Savchenkov, A. B. Matsko, W. Liang, V. S. Ilchenko, D. Seidel, and L. Maleki, “Kerr combs with selectable central frequency,” Nat. Photonics 5(5), 293–296 (2011). [CrossRef]  

19. Y. Liu, Y. Xuan, X. Xue, P.-H. Wang, S. Chen, A. J. Metcalf, J. Wang, D. E. Leaird, M. Qi, and A. M. Weiner, “Investigation of mode coupling in normal-dispersion silicon nitride microresonators for Kerr frequency comb generation,” Optica 1(3), 137–144 (2014). [CrossRef]  

20. N. M. Kondratiev and V. E. Lobanov, “Modulational instability and frequency combs in whispering-gallery-mode microresonators with backscattering,” Phys. Rev. A 101(1), 013816 (2020). [CrossRef]  

21. A. A. Savchenkov, A. B. Matsko, W. Liang, V. S. Ilchenko, D. Seidel, and L. Maleki, “Kerr frequency comb generation in overmoded resonators,” Opt. Express 20(24), 27290–27298 (2012). [CrossRef]  

22. N. M. Kondratiev, V. E. Lobanov, and D. V. Skryabin, “Backward-wave induced modulational instability in normal dispersion,” in Laser Congress 2019 (ASSL, LAC, LS&C), OSA Technical Digest (Optica Publishing Group, 2019), paper JTu3A.32.

23. H. A. Haus and W. Huang, “Coupled-mode theory,” Proc. of the IEEE 79, 10 (1991).

24. S. Ramelow, A. Farsi, S. Clemmen, J. S. Levy, A. R. Johnson, Y. Okawachi, M. R. E. Lamont, M. Lipson, and A. L. Gaeta, “Strong polarization mode coupling in microresonators,” Opt. Lett. 39(17), 5134–5137 (2014). [CrossRef]  

25. C. Bao, Y. Xuan, D. E. Leaird, S. Wabnitz, M. Qi, and A. M. Weiner, “Spatial mode-interaction induced single soliton generation in microresonators,” Optica 4(9), 1011–1015 (2017). [CrossRef]  

26. X. Xue, Y. Xuan, P.-H. Wang, Y. Liu, D. E. Leaird, M. Qi, and A. M. Weiner, “Normal-dispersion microcombs enabled by controllable mode interactions,” Laser Photon. Rev. 9(4), L23–L28 (2015). [CrossRef]  

27. G. Lin and T. Sun, “Mode crossing induced soliton frequency comb generation in high-Q yttria-stabilized zirconia crystalline optical microresonators,” Photonics Res. 10(3), 731–739 (2022). [CrossRef]  

28. S. Fujii, Y. Okabe, R. Suzuki, T. Kato, A. Hori, Y. Honda, and T. Tanabe, “Analysis of Mode Coupling Assisted Kerr Comb Generation in Normal Dispersion System,” IEEE Photonics J. 10(5), 1–11 (2018). [CrossRef]  

29. P.-H. Wang, K.-L. Chiang, and Z.-R. Yang, “Study of microcomb threshold power with coupling scaling,” Sci Rep 11(1), 9935 (2021). [CrossRef]  

30. ÓB Helgason, A. Fülöp, J. Schröder, P. A. Andrekson, A. M. Weiner, and V. Torres-Company, “Superchannel engineering of microcombs for optical communications,” J. Opt. Soc. Am. B 36(8), 2013–2022 (2019). [CrossRef]  

31. Z. Qi, S. Wang, J. Jaramillo-Villegas, M. Qi, A. M. Weiner, G. D’Aguanno, T. F. Carruthers, and C. R. Menyuk, “Dissipative cnoidal waves (Turing rolls) and the soliton limit in microring resonators,” Optica 6(9), 1220–1232 (2019). [CrossRef]  

32. J. A. Jaramillo-Villegas, X. Xue, P.-H. Wang, D. E. Leaird, and A. M. Weiner, “Deterministic single soliton generation and compression in microring resonators avoiding the chaotic region,” Opt. Express 23(8), 9618–9626 (2015). [CrossRef]  

33. A. M. Weiner, Ultrafast Optics (Wiley, 2009), Chap. 6.

34. M. Karpov, M. H. P. Pfeiffer, H. Guo, W. Weng, J. Liu, and T. J. Kippenberg, “Dynamics of soliton crystals in optical microresonators,” Nat. Phys. 15(10), 1071–1077 (2019). [CrossRef]  

35. T. Huang, J. Pan, Z. Cheng, G. Xu, Z. Wu, T. Du, S. Zeng, and P. P. Shum, “Nonlinear-mode-coupling-induced soliton crystal dynamics in optical microresonators,” Phys. Rev. A 103(2), 023502 (2021). [CrossRef]  

36. P. H. Wang, Y. Xuan, L. Fan, L. T. Varghese, J. Wang, Y. Liu, X. Xue, D. E. Leaird, M. Qi, and A. M. Weiner, “Drop-port study of microresonator frequency combs: power transfer, spectra and time-domain characterization,” Opt. Express 21(19), 22441–22452 (2013). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1.
Fig. 1. Mapping of MI gain versus the mode number and pump power for (a) ${\Delta _{shift}}$  = 0, (b) ${\Delta _{shift}}$  = 6, and (c) ${\Delta _{shift}}$  = 8. The corresponding spectra at S2 = 10 and l = 5 are shown in (d)-(f), respectively.
Fig. 2.
Fig. 2. Spectra for different pump powers S2 = 2, 10, and 50. The microcomb is initiated from the random noise at Δ = 0, ${\Delta _{shift}}$  = 8, l = 5, and S2 = 10 but no oscillation can be found at S2 = 2 and 50.
Fig. 3.
Fig. 3. MI gain spectra without ( ${\Delta _{shift}}$  = 0) for (a) l = 5 and (b) l = 8 and (c) with ( ${\Delta _{shift}}$  = 8) mode coupling for l = 5. The corresponding spectra (d) without and (e) with mode coupling.
Fig. 4.
Fig. 4. Maps of eigenvalue bifurcations in normal dispersion in related to different (a) effective coupling strength y’, (b) normalized dispersion $\beta $ , (c) mode number l, and (d) coupling strength ${\Delta _{shift}}$ .
Fig. 5.
Fig. 5. Maps of comb efficiency in $\Delta $ -S2 space under red-detuning for (a) l = 1 and ${\Delta _{shift}}$  = 3, (b) l = 1 and ${\Delta _{shift}}$  = 6, (c) l = 5 and ${\Delta _{shift}}$  = 3, and (d) l = 5 and ${\Delta _{shift}}$  = 6, and (e)(h) the corresponding maps under blue-detuning.
Fig. 6.
Fig. 6. Maps of comb efficiency in $\Delta $ -S2 space under blue-detuning for (a) l = 1 and ${\Delta _{shift}}$  = 6 and (b) l = 5 and ${\Delta _{shift}}$  = 6. The corresponding zoom-in spectra are also shown on the right.
Fig. 7.
Fig. 7. (a) Normalized intracavity power versus phase detuning at S2 = 2, 5, and 8. (b) The enhancement factor of the intracavity power when the phase detuning is continuously increasing (under red-detuning). The red-dashed and black-dashed lines represent the normalized intracavity pump ${|{{A_0}} |^2} = 1$ and the two equilibria of the intracavity pump, respectively.
Fig. 8.
Fig. 8. Maps of comb efficiency in $\Delta $ -S2 space under red-detuning for (a) l = 1 and ${\Delta _{shift}}$  = 1.5.

Equations (7)

Equations on this page are rendered with MathJax. Learn more.

T R E ( t , θ ) t = [ ( α + i δ 0 ) i L β 2 2 ( 2 π T R ) 2 2 θ 2 + i γ L | E | 2 ] E + κ E i n
A ( t , θ ) t = [ ( 1 + i Δ ) i β 2 2 θ 2 + i | A | 2 ] A + S
A ( t , θ ) = l A l ( t ) e i l FSR θ 2 π T R = l A l ( t ) e i l θ
A l ( t ) t = [ ( 1 + i Δ ) + i β 2 l 2 ] A l + δ ( l ) S + i m , n , p δ ( m n + p l ) A m A n A p
A l ( t ) t = [ ( 1 + i ( Δ + Δ s h i f t ) ) + i β 2 l 2 ] A l + δ ( l ) S + i m , n , p δ ( m n + p l ) A m A n A p
[ d δ A l ( t ) d t d δ A l ( t ) d t ] = [ M l N N M l ] [ δ A l ( t ) δ A l ( t ) ]
λ ( l ) = 1 i Δ s h i f t 2 ± ( Δ 2 | A 0 | 2 + Δ s h i f t 2 β l 2 2 ) 2 + | A 0 | 4
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.