Soliton and quasi-soliton frequency combs due to second harmonic generation in microresonators

TL;DR

This paper explores how second harmonic generation in microring resonators with quadratic nonlinearity leads to various soliton and quasi-soliton frequency combs, revealing new nonlinear resonance behaviors and soliton states.

Contribution

It demonstrates the formation of different soliton types in quadratic nonlinear microresonators and links their transition to spectral gap closure, advancing understanding of nonlinear optical combs.

Findings

Bright quasi-solitons on low intensity background

Exponential localized solitons with same dispersion signs

Transition linked to spectral gap closure

Abstract

We report how a doublet of the symmetric oppositely tilted bistable resonance peaks in a microring resonator with quadratic nonlinearity set for generation of the second harmonic can transform into a Kerr-like peak on one side of the linear cavity resonance and into a closed loop structure disconnected from the quasi-linear resonance on the other. Both types of the nonlinear resonances are associated with the formation of the soliton combs for dispersion profiles of a typical LiNbO$_3$ microring. We report bright quasi-solitons propagating on a weakly modulated low intensity background when the group velocity dispersions have the opposite signs for the fundamental and second harmonic. We also show exponentially localized solitons when the dispersion signs are the same. Finally, we demonstrate that the transition between these two types of soliton states is associated with the closure of the forbidden gap in the spectrum of quasi-linear waves.

Figures

Equations

$$i{\partial }_T{\Psi }_f=\left({\omega }_f-i{D}_{1f}{\partial }_{\varphi }-\frac{1}{2}{D}_{2f}{\partial }_{\varphi }^2\right){\Psi }_f-i{\kappa }_f{\Psi }_f-\alpha {\Psi }_f^{\ast }{\Psi }_s-\alpha h{e}^{-i{\omega }_{p}T},$$

$$i{\partial }_T{\Psi }_s=\left({\omega }_s-i{D}_{1s}{\partial }_{\varphi }-\frac{1}{2}{D}_{2s}{\partial }_{\varphi }^2\right){\Psi }_s-i{\kappa }_s{\Psi }_s-\alpha {\Psi }_f^2.$$

$$i{\partial }_t{\psi }_f=-2\pi \left(iV{\partial }_{\theta }+\frac{1}{2}{d}_{2f}{\partial }_{\theta }^2\right){\psi }_f+({\delta }_f-i{\gamma }_f){\psi }_f-{\psi }_f^{\ast }{\psi }_s-h,$$

$$i{\partial }_t{\psi }_s=-2\pi \left(i(V-U){\partial }_{\theta }+\frac{1}{2}{d}_{2s}{\partial }_{\theta }^2\right){\psi }_s+({\delta }_s-i{\gamma }_s){\psi }_s-{\psi }_f^2.$$

$${\psi }_{f,s}(\theta )={\psi }_{f0,s0}+{\zeta }_{f,s}{e}^{i\mu \theta }+{\xi }_{f,s}^{\ast }{e}^{-i\mu \theta }.$$

Full text

OCIS codes: (190.4380) Nonlinear optics, four-wave mixing; (140.3945) Microcavities; (190.1450) Bistability

1 Introduction

Frequency comb generation in microring resonators with Kerr nonlinearity has been intensely studied in the last decade. These studies have demonstrated a plethora of novel solitonic effects with immediate applications in the precision measurements and optical signal processing [1, 2]. Using nonlinear effects other than Kerr nonlinearity for comb generation is an active research area. One can mention here examples of the comb generation due to cascaded Raman and Brillouin effects [3, 4, 5, 6], polariton [7, 8] and graphene nonlinearities [9], and quadratic, $\chi^{(2)}$, nonlinear effects [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. The latter hold a particular promise. $\chi^{(2)}$ materials are widely available and well studied for frequency conversion and soliton applications. They provide high levels of nonlinear response allowing to reduce pump power requirements. Importantly, phase matching for the three wave mixing in $\chi^{(2)}$ materials and the soliton formation conditions do not require anomalous dispersion, as it takes place in Kerr, $\chi^{(3)}$, materials [23]. This significantly extends the playing field for the comb research in terms of enhancing flexibility of the pump source choices and generating combs over wider spectral ranges between mid-infrared and ultraviolet.

Following pioneering observations of 2nd harmonic generation in a whispering gallery mode resonator [24], research into frequency conversion for those devices was limited by the difficulty in achieving a wide mode spectra generated in the proximity of both pump and signal waves over some period of time, see [25] for a review. Only recently the first convincing experimental observation of a broad frequency combs and a cascaded multi-wave mixing has been reported in a micro-resonator with $\chi^{(2)}$ nonlinearity [18, 19]. Since most practical low noise and, at the same time, broad combs are provided by bright dissipative optical solitons [1] in microresonators, theoretical understanding of cavity solitons supported by quadratic nonlinearity is essential. We note that first theories on $\chi^{(2)}$ cavity solitons date two decades back, we recall for example early studies on diffractive solitons in planar microcavities, see, e.g., [26, 27]. However, $\chi^{(2)}$ cavity solitons had not received much of the experimental attention so far, in contrast to their Kerr couterparts [1]. Nowadays, due to the growing interest in comb generation in $\chi^{(2)}$ microresonators [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22], the theoretical studies of quadratic solitons are particularly relevant to stimulate and guide growing experimental efforts in this area. We recall, for example, the recent experiment in a bow-tie cavity, on the observation of $\chi^{(2)}$ simultons (coupled bright and dark pulses) under conditions where the 2nd order dispersion is negligible [16]. The theoretical problems that are relevant at present and that are connected to the study below are the frequency comb generation under conditions of the normal and mixed group velocity dispersion signs at the fundamental and second harmonics and dispersive wave emission by solitons, which is expected to have a different effect in a ring geometry compare to the waveguides.

The paper is structured as follows: Section 2 introduces the reader into the Lugiato-Lefever model and maps its parameters on the dispersion profiles corresponding to a LiNbO3 microring wire resonator. Section 3 covers single mode solutions and their nonlinear resonances. In Section 4 theoretical considerations concerning the criteria for the existence of quadratic solitons are discussed. In particular, we explain why these type of solutions belong to a class of gap solitons. Numerical results on the soliton and quasi-soliton states are described in Section 5. Our work has been performed and results have been obtained independently from [20] on comb solitons due to second harmonic generation published in December 2018. While finalising our material for publication we have not included data overlapping with [20] related to the case of the symmetric resonance splitting as in our Fig. 2(a).

2 Lugiato-Lefever model for a microring with $\chi^{(2)}$ nonlinearity

To simulate the spatio-temporal evolution of the fundamental and 2nd harmonic fields in a microring resonator we use the following model equations

[TABLE]

Here $T$ is the physical time and $\phi$ is the polar angle varying along the resonator circumference. The amplitudes of the dimensionless electric fields are given by $\Psi_{f}$ and $\Psi_{s}$ for the fundamental and for the 2nd harmonic, respectively. $\omega_{p}$ is the frequency of the CW pump laser while $\omega_{f}$ corresponds to its closest resonant frequency in the linear spectrum of the microring. Similarly, $\omega_{s}$ is the nearest resonance to the 2nd harmonic of the pump, $2\omega_{p}$. We have assumed the exact angular phase matching conditions here, that can be achieved through the quasi-phase matching technique [25]. Accounting for the time derivatives and the bracketed expressions in both Eqs. (1) and (2) allows to approximate the microring dispersion around the pump and its 2nd harmonic. The coefficients $D_{1f,1s}$ correspond to the local free spectral ranges (FSRs) of the microresonator and $D_{2f,2s}$ set the FSR dispersion (group velocity dispersion). Finally, $\kappa_{f}$ and $\kappa_{s}$ represent the photon loss rates, $h$ is the dimensionless amplitude of the CW pump and $\alpha$ is a scaling parameter with units of frequency.

Our primary interest here is in finding the soliton combs, therefore we look for shape preserving pulses rotating around the microring with the constant frequency $\Omega$: $\Psi_{f}(\phi,T)=\psi_{f}(\theta,t)e^{-i\omega_{p}T}$, $\Psi_{s}(\phi,T)=\psi_{s}(\theta,t)e^{-i2\omega_{p}T}$, where $\phi=\theta-\Omega T$ and $t=T/T_{f}$ is the dimensionless time normalized to the pulse round trip time $T_{f}=2\pi/D_{1f}$ at the fundamental frequency $2\pi D_{1f}$. Then $\alpha$ is conveniently fixed as $\alpha=1/T_{f}$. Let us stress that $\Omega$ is the same for the fundamental and 2nd harmonic fields, so that if $\psi_{s,f}$ are assumed to be $t$ independent, then both fields will be permanently locked one to another. We also note that this assumption can be practically satisfied only if $D_{1f}$ and $D_{1s}$ are sufficiently close, as will be discussed in detail below. It is therefore convenient to assume that $\Omega=D_{1f}(1-V)$, where $V$ is proportional to the FSR shift induced by the nonlinearity. The above substitutions result in the dimensionless system of equations convenient for our numerical and theoretical studies:

[TABLE]

Here $\gamma_{f,s}=2\pi\kappa_{f,s}/D_{1f}$ and $d_{2f,2s}=D_{2f,2s}/D_{1f}$ are the normalized losses and dispersion parameters. We also introduce $\delta_{f}=2\pi(\omega_{f}-\omega_{p})/D_{1f}$ and $\delta_{s}=2\pi(\omega_{s}-2\omega_{p})/D_{1f}=2\delta_{f}+\epsilon$ which represent the normalized detuning parameters. Due to the discrete nature of the spectrum of a microresonator and dispersion, even if we tune $\omega_{p}$ to perfectly match $\omega_{f}$, the 2nd harmonic $\omega_{s}$ will in general be off-resonance. To carefully account for that we define $\epsilon=2\pi(\omega_{s}-2\omega_{f})/D_{1f}$ to be a constant offset between the resonance frequency $\omega_{s}$ nearest to $2\omega_{p}$ and the doubled frequency $\omega_{p}$ of the fundamental resonance. The importance of $\varepsilon$ is in shaping the nonlinear response of the microresonator, and its effect will be discussed in detail in the following section. Finally, $U=(D_{1f}-D_{1s})/D_{1f}$ is the normalized difference between the FSRs at the fundamental and at 2nd harmonic frequencies. This parameter characterizes the difference of the pulse rotation frequencies of the fundamental and 2nd harmonic fields. In the soliton regime $U$ is compensated through nonlinear effects [23], so that the two fields lock and propagate together. $U$ can also be minimised or eliminated by adjusting the pump frequency and cavity dispersion, see below. We note, that in the later case $U=0$ implies the selection of a particular value of $\epsilon$.

Considering a toroidal microring with the large radius $100\,\mu$m and the small one $1\,\mu$m made of LiNbO3 and using the theory in [28] we estimate values of the dispersion and detuning parameters to be used in our numerical analysis below. These estimates can also be used as a qualitative guide for the planar on-chip wire microresonators, see, e.g., [12, 22]. Considering the same polarisation states for the fundamental and 2nd harmonic fields we found that it is possible to match FSRs at both frequencies across the zero dispersion point, see Fig. 1(a). Using an ordinary wave for the fundamental and the extraordinary for the 2nd harmonic it is possible to match $D_{1f}$ to $D_{1s}$ maintaining $D_{2f}$ and $D_{2s}$ negative, see Fig. 1(b). Note, that the negative/positive $D_{2}$ corresponds here to the normal/anomalous dispersion.

3 Nonlinear resonances and bistability

We start our analysis of the nonlinear microring states from the CW-solutions that correspond to homogeneous solutions of Eqs. (3) and (4) ($\partial_{t}=\partial_{\theta}=0$) and we check how quadratic nonlinearity distorts the shape of the linear resonance response. Figure 2 shows how the amplitude of CW solutions vary with $\delta_{f}$ for several values of $\varepsilon$. Note that changing $\delta_{f}$ corresponds experimentally to tuning the pump laser frequency $\omega_{p}$ around $\omega_{f}$. We also stress that during this process the value of $\varepsilon$ remains constant. In Fig. 2(a) we consider the ideal case when $\varepsilon=0$. For such a value of $\varepsilon$ the linear resonance splits symmetrically and two tilted bistable resonances are formed  [20]. To capture all relevant properties of the system one needs to account for $\varepsilon\neq 0$, which leads to a shift and reshaping of the structure of nonlinear resonances. When $\varepsilon>0$, see Figs. 2(b)-2(d), the most of reshaping appears in the negative range of values of $\delta_{f}$, while the situation is reversed when $\epsilon<0$. By further increasing the value of $\epsilon$, the negatively tilted resonance undergoes a distortion which narrows down its bistability range. This process ultimately leads to the detachment of the tilted resonance from the CW low amplitude background, see Fig. 2(c). We stress that the latter scenario can have an important implication concerning the experimental technique used to achieve solitons in microresonators. It is well known that the Kerr soliton combs can be obtained by adiabatic scanning of the laser pump across a resonance [1]. In this way, unstable high amplitude CW solutions tend to decay into nearby soliton states. However, the presence of a resonance loop disconnected from low amplitude state may compromise this experimental technique, which suggests that hard excitation methods might be a more suitable choice to achieve comb-soliton states. While stability properties of the upper branches of the CW-solutions involve many details reporting on which goes beyond our present scope, the low amplitude solution are linearly stable through ranges of the soliton existence shown with full lines in Fig. 5 below.

4 Soliton existence conditions

In this Section we analyse the comb soliton existence conditions with two different approaches. For sufficiently large $\varepsilon$, see Fig. 2(d), the nonlinear resonance effects for negative detunings become suppressed. The system enters into the Kerr-like regime, also called the cascading nonlinearity regime [23]. Here the 2nd harmonic simply adiabatically follows the fundamental one, such that $\psi_{s}\simeq\psi_{f}^{2}/\delta_{s}$. The latter approximation reduces Eqs. (3) and (4) to the standard Lugiato-Lefever model for stationary comb solitons [1]: $\pi d_{2f}\partial^{2}_{\theta}\psi_{f}\simeq\delta_{f}\psi_{f}-|\psi_{f}|^{2}\psi_{f}/\delta_{s}-(h+2i\pi V\partial_{\theta}+i\gamma_{f})\psi_{f}$. One of the advantages of working with a cascading $\chi^{(2)}$ nonlinearity, over a pure Kerr system, is that the sign of the cascading nonlinearity is controlled by the value of $\delta_{s}=2\delta_{f}+\varepsilon$. Hence bright soliton combs can be expected for both normal ($D_{2f}<0$, $\delta_{f,s}<0$) and anomalous ($D_{2f}>0$, $\delta_{f,s}>0$) dispersion.

Since the role played by the 2nd harmonic field in the soliton existence condition remains obscure in the cascading limit analysis, we are analysing the soliton existence conditions using a more formal, but physically grounded approach. For high quality microresonators the loss term can be disregarded in first approximation. This approximation allows to analyse qualitatively and quantitatively the role played by the velocity parameter $V$ and by the 2nd harmonic field in the soliton formation. The low amplitude state of the nonlinear resonances is known to serve as the background state for the bright comb solitons. Thus one can say that the necessary condition for the soliton existence is that this state does not support the extended modes representing quasi-linear small amplitude perturbations distorting the soliton. In particular, we are anticipating that for some parameter values there should exists a gap in $V$-values such that for all modal numbers $\mu$ one can not find a $V$ value inside the gap giving a physically realizable quasi-linear periodic wave forms. We assume that such quasi-linear waves can be sought as $t$-independent solutions of Eqs. (3) and (4) in the form

[TABLE]

Here $\psi_{f0}$ and $\psi_{s0}$ are the fundamental and 2nd harmonic amplitudes of the lowest amplitude state within the bistability ranges of the nonlinear resonances, while $\zeta_{f,s}$ and $\xi_{f,s}$ are small amplitudes of the sought waves. We note that $\mu$ is one of the modes of the microresonator counted from zero. Here $\mu=0$ corresponds to the mode associated to $\omega_{f}$. Substituting Eqs. (5) into Eqs. (3) and (4) and linearising for small $\zeta_{f,s}$ and $\xi_{f,s}$ we find that, for the solutions to exist, the velocity parameter has to satisfy the fourth order equation $a_{4}V^{4}+a_{3}V^{3}+a_{2}V^{2}+a_{1}V+a_{0}=0$, where all $a_{n}$’s are parameterized by $\mu$.

Figure 3 shows typical plots of four roots of $V$ vs $\mu$. Figure 3(a) demonstrates the situation when dispersion signs are the same and the mismatch of group velocities $U$ is sufficiently large to overcome the nonlinear effects to pull two pulses apart. This is the physical argument prohibiting the existence of the mutually locked fundamental and 2nd harmonic pulses. This argument matches with our gap analysis results showing that the small amplitude periodic waves are allowed for all $V$’s for these system parameters. However, when $U$ is sufficiently small or zero, then an interval (gap) of $V$ values, where the small amplitude periodic waves are forbidden, is emerging, see Fig. 3(b) and cf. Fig. 1(b). We expect that the bright comb solitons can exist under these conditions since a pulsed excitation with, e.g., $V=0$, will not couple to a continuum of quasi-linear waves. When $D_{2f}>0$ and $D_{2s}<0$, see Fig. 3(c) and cf. Fig. 1(a), the gap between the two roots closes for some relatively large $\mu=\mu_{cr}$ at $V=0$ and even when $U=0$ is chosen. This implies, in the spirit of the theory of dispersive wave emission by solitons [29], that the quasi-soliton excitation with $V=0$ is going to develop tails with a small amplitude dispersive wave at $\mu=\mu_{cr}$. If $\mu_{cr}$ is sufficiently large ($\mu_{cr}$ tends to infinity with $D_{2f}\to 0$), then an impact of the dispersive waves on the soliton core is expected to be negligible, see previous theories for the dispersive wave emission by Kerr solitons [29].

5 Bright soliton combs

We will now numerically investigate existence and dynamics of bright soliton solutions, which spectrally correspond to the frequency comb states due to 2nd harmonic generation. Throughout this study, the parameters $D_{2f,2s}$ and $\varepsilon$ are chosen using the dispersion data in Fig. 1. Numerical simulations to find soliton families in the time-independent version of Eqs. 3 and 4 are performed using a Newton method in combination with a Bi-conjugate gradient method implemented in the Fourier space, such that periodic boundary conditions are accounted for automatically. To simulate soliton dynamics we use standard pseudo spectral approach. All our results are obtained considering 2048 modes.

We start from the case when dispersions are normal ($D_{2f,2s}<0$) for both fundamental and 2nd harmonic, which is the case shown in Fig. 1(b). Bright soliton combs in a material with positive Kerr do not exist under these common conditions. However, $\chi^{(2)}$ effects allow their existence. The dimensionless parameters in Eqs. (3) and (4) which allow for the FSR matching are estimated to be $\epsilon=0.51$, $d_{2f}=-1.83\times 10^{-5}$ and $d_{2s}=-6.59\times 10^{-5}$. We note that these parameters are the same as the ones used in Fig. 3(b). This suggests the presence of an open gap in the soliton velocity parameter $V$. More specifically, both the $V$-gap analysis and the cascading approximation predict the existence of solitons when negative values of detuning are considered. Solitons are indeed found when $\delta_{f}<0$. Figures 4(a)-4(c) show the spatial profiles of the two soliton components and their corresponding frequency combs. The structure of both nonlinear resonances and soliton families for different values of $\delta_{f}$ are shown in Fig. 5(a). We note that, continuing the soliton family solutions in $\delta_{f}$, we found a closed loop structure similar to the one for the cw state, see the zoomed inset in Fig. 5(a). This is true for both fundamental and 2nd harmonic solitons family. Performing linear stability analysis we have confirmed that the low amplitude soliton state in the loop is unstable. Its dynamics, in fact, leads either to the excitation of an attractor in the proximity of the high amplitude state or to the decay into the lower amplitude stable background. On the other hand, the soliton state with the larger amplitude can be stable. We note the presence of some instabilities for small values of $\delta_{f}$, see Fig. 5(a), but reporting and studying peculiarities of those goes beyond the aim of this work.

We will now discuss the case of different signs of dispersions: anomalous ($D_{2f}>0$) for the fundamental harmonic and normal ($D_{2s}<0$) for the 2nd harmonic. In particular, we will consider the scenario when both fields are extraordinary polarized, see Fig. 1(a). The dimensionless parameters in Eqs. (3) and (4) which allow for the FSR matching are given by $\epsilon=0.49$, $d_{2f}=6.89\times 10^{-5}$ and $d_{2s}=-2.94\times 10^{-5}$. Similarly to the previous case, these parameters have already been used in Fig. 3(c), which does not show any open $V$-gap. Despite the absence of an open gap, the Lugiato-Lefever model with Kerr type of nonlinearity, is approximately valid in the cascading limit, and still predicts the existence of bright solitons for positive detunings. Solitons are indeed found when $\delta_{f}>0$ and their profiles and the corresponding spectral combs are shown in Figs. 4(d)-4(f). In order to explain this contradicting result one needs to account for the finite size of the microresonator. The theory of dispersive waves emission by solitons [29] implies that solitons in this situation are not truly localized objects (quasi-solitons), but they will transfer energy to small amplitude dispersive waves with the modal index closely matching the mode number corresponding to the crossing of different $V$-branches within the former gap. In our case the critical mode number is given by $\mu_{cr}=\pm 82$. However, being in a confined system the transfer of energy from the soliton to the dispersive (Cherenkov) radiation does not necessarily causes the decay of the soliton. This is due to the fact that the radiation travels along the ring and transfers energy back to the soliton. A family of these quasi-solitons was numerically traced in $\delta_{f}$ and is shown in Fig. 5(b) together with the corresponding CW solutions. We note that here the absence of the left tilted resonance is due to the relatively large choice of $\gamma_{f,s}=0.005$.

A feature of the quasi-solitons in a ring geometry is that the radiation creates an effective potential landscape for the soliton which makes it more robust with respect to the translational shifts. This has been reported for Kerr comb solitons with the 4th order dispersion [30] and is expected here. The radiation amplitude is significantly stronger for the second harmonic, which corroborates the hypothesis that the gap closure in Fig. 3(c) involves the dispersion branches $V(\mu)$ associated with the 2nd harmonic field. In fact, if the cascading limit is considered, the dispersive branches for the 2nd harmonic are eliminated and exponentially damped tails are recovered. We stress that these strong dispersive resonances substantially amplify spectral tails of the comb around the 2nd harmonic, which can be utilised for practical purposes.

Finally, we confirm the comb generation using dynamical simulations. More specifically, we found that the solitons emerge from modulational instabilities of the upper branch of the CW solutions, see Figs. 6(a)-6(b) and Fig. 6(e)-(f). Using the same parameters of Fig. 5, plots in Figs. 6(a)-6(d) represent the case with the same signs of dispersions, while plots in Figs. 6(e)-6(h) correspond to the opposite signs. These simulations not only reproduce combs similar to the ones in Fig. 4, but also confirm the robustness of $\chi^{(2)}$ comb solitons. We note, that solitons also remain stable against perturbations coming from the background radiation waves, which are particularly strong during the initial stages of the soliton evolution.

6 Conclusions

We have developed a theory that demonstrates that comb solitons in microring resonators with quadratic nonlinearity belong to a class of gap-solitons. We demonstrated the existence of the families of comb solitons for normal dispersion at both the fundamental and second harmonic fields. Some of these solitons exist under conditions of the strongly modified nonlinear resonances, that are not smoothly connected to the linear resonance tails. When the solitonic gap is closed through the action of the opposite dispersion signs at the fundamental and second harmonic, we have found that the spectral tails of the second harmonic comb acquire strong resonance peaks corresponding to the dispersive wave emitted by the soliton. The associated spectral broadening may find practical applications. Dispersion values and group velocity matching conditions in our work have been calculated for the resonator made of the $1\mu$m wide LiNbO3 waveguide.

Funding

The Leverhulme Trust (RPG-2015-456); H2020 (691011, Soliring); Russian Science Foundation (17-12-01413).

Acknowledgments

We thank I. Breunig and J. Szabados for helpful suggestions and discussions.

References and links

• [1]

T.J.  Kippenberg, A.L.  Gaeta, M.  Lipson, and M.L.  Gorodetsky, “Dissipative Kerr solitons in optical microresonators,” Science 361, eaan8083 (2018).

• [2]

D. V Strekalov, C. Marquardt, A. B. Matsko, H. G. L. Schwefel, and G. Leuchs, “Nonlinear and quantum optics with whispering gallery resonators,” J. Opt. 18, 123002 (2016).

• [3]

Y. Okawachi, M. J. Yu, V. Venkataraman, P. M. Latawiec, A. G. Griffith, M. Lipson, M. Loncar, and A. L. Gaeta, “Competition between Raman and Kerr effects in microresonator comb generation,” Opt. Lett. 42, 2786-2789 (2017).

• [4]

A. V. Cherenkov, N. M. Kondratiev, V. E. Lobanov, A. E. Shitikov, D. V. Skryabin, and M. L. Gorodetsky, “Raman-Kerr frequency combs in microresonators with normal dispersion,” Opt. Express 25, 31148-31158 (2017).

• [5]

S. Gundavarapu, G.M. Brodnik, M. Puckett, T. Huffman, D. Bose, R. Behunin, J. F. Wu, T. Q. Qiu, C. Pinho, N. Chauhan, J. Nohava, P. T. Rakich, K. D. Nelson, M. Salit, and D. J. Blumenthal, “Sub-hertz fundamental linewidth photonic integrated Brillouin laser,” Nature Photon. 13, 60-67 (2019).

• [6]

M. Sumetsky, “Optical frequency combs generated mechanically,” Opt. Lett. 42, 3197-3200 (2017).

• [7]

R. Castro-Beltran, V.M. Diep, S. Soltani, E. Gungor, and A.M. Armani, “Plasmonically Enhanced Kerr Frequency Combs,” ACS Photonics 4, 2828-2834 (2017).

• [8]

O. A. Egorov, and D. V. Skryabin, “Frequency comb generation in a resonantly pumped exciton-polariton microring resonator,” Opt. Express 26, 24003-24009 (2018).

• [9]

B.C. Yao, S.W. Huang, Y. Liu, A.K. Vinod, C. Choi, M. Hoff, Y.N. Li, M.B. Yu, Z.Y. Feng, D.L. Kwong, Y. Huang, Y.J. Rao, X.F. Duan, and C.W. Wong, “Gate-tunable frequency combs in graphene-nitride microresonators,” Nature 558, 410-414 (2018).

• [10]

X. Xue, X. Zheng, and A. M. Weiner,

“Soliton trapping and comb self-referencing in a single microresonator with $\chi^{(2)}$ and $\chi^{(3)}$ nonlinearities,”

Opt. Lett. 42, 4147-4150 (2017).

• [11]

X. Xue, F. Leo, Y. Xuan, J. A Jaramillo-Villegas, P.-H. Wang, D. E. Leaird, M.  Erkintalo, M.  Qi, and A. M. Weiner.

“Second-harmonic-assisted four-wave mixing in chip-based microresonator frequency comb generation,”

Light Sci. Appl. 6, e16253 (2017).

• [12]

R.  Wolf, I.  Breunig, H.  Zappe, and K.  Buse,

“Cascaded second-order optical nonlinearities in on-chip micro rings,”

Opt. Express 25, 29927-29933 (2017).

• [13]

S. Mosca, M. Parisi, I. Ricciardi, F. Leo, T. Hansson, M. Erkintalo, P. Maddaloni, P. De Natale, S. Wabnitz, and M. De Rosa,

“Modulation instability induced frequency comb generation in a continuously pumped optical parametric oscillator,”

Phys. Rev. Lett. 121, 093903 (2018).

• [14]

S. K.  Meisenheimer, J. U. Fürst, K.  Buse, and I.  Breunig,

“Continuous-wave optical parametric oscillation tunable up to an 8$\mu$m wavelength,”

Optica 4, 189-192 (2017).

• [15]

F. Leo, T. Hansson, I. Ricciardi, M. De Rosa, S. Coen, S. Wabnitz, and M. Erkintalo, “Walk-off-induced modulation instability, temporal pattern formation, and frequency comb generation in cavity-enhanced second-harmonic generation,”

Phys. Rev. Lett. 116, 033901 (2016).

• [16]

M.  Jankowski, A.  Marandi, C. R. Phillips, R.  Hamerly, K. A. Ingold, R. L. Byer, and M. M. Fejer,

“Temporal simultons in optical parametric oscillators,”

Phys. Rev. Lett. 120, 053904 (2018).

• [17]

R.  Ikuta, M.  Asano, R.  Tani, T.  Yamamoto, and N.  Imoto.

“Frequency comb generation in a quadratic nonlinear waveguide resonator,”

Opt. Express 26, 15551-15558 (2018).

• [18]

S. J. Herr, V.  Brasch, J.  Szabados, E.  Obrzud, Y. Jia, S.  Lecomte, K.  Buse, I.  Breunig, and T.  Herr,

“Frequency comb up- and down-conversion in synchronously driven $\chi^{(2)}$ optical microresonators,” Opt. Lett. 43, 5745-5748 (2018).

• [19]

D.J. Wilson, K. Schneider, S. Hoenl, M. Anderson, T.J. Kippenberg, and P. Seidler, “Integrated gallium phosphide nonlinear photonics,” arXiv e-prints, arXiv:1808.03554 (2018).

• [20]

T. Hansson, P. Parra-Rivas, M. Bernard, F. Leo, L. Gelens, and S. Wabnitz. “Quadratic Soliton Combs in Doubly-Resonant Second-Harmonic Generation,” Opt. Lett. 43, 6033-6036 (2018).

• [21]

T. Hansson, F. Leo, M. Erkintalo, S. Coen, I. Ricciardi, M. De Rosa, and S. Wabnitz, “Singly resonant second-harmonic-generation frequency combs,” Phys. Rev. A 95, 013805 (2017).

• [22]

M. Zhang, C. Wang, R. Cheng, A. Shams-Ansari, and M. Loncar, “Monolithic ultra-high-Q lithium niobate microring resonator,” Optica 4, 1536-1537 (2017).

• [23]

A.V. Buryak, P. Di Trapani, D.V. Skryabin, and S. Trillo, “Optical solitons due to quadratic nonlinearities: From basic physics to futuristic applications,” Phys. Rep. 370, 63-235 (2002).

• [24]

J.U. Fürst, D.V. Strekalov, D. Elser, M. Lassen, U. L. Andersen, C. Marquardt, and G. Leuchs, “Naturally phase-matched secondharmonic generation in a whispering-gallery-mode resonator,” Phys. Rev. Lett. 104, 153901 (2010).

• [25]

I. Breunig, “Three-wave mixing in whispering gallery resonators,” Laser Photon. Rev. 10, 569-587 (2016).

• [26]

C. Etrich, U. Peschel, and F. Lederer, “Solitary waves in quadratically nonlinear resonators,” Phys. Rev. Lett. 79, 2454-2457 (1997).

• [27]

D.V. Skryabin, A.R. Champneys, and W.J. Firth, “Frequency selection by soliton excitation in nondegenerate intracavity down-conversion,” Phys. Rev. Lett. 84, 463-466 (2000).

• [28]

M.L. Gorodetsky, and A.E. Fomin, “Geometrical theory of whispering-gallery modes,” IEEE J. Sel. Top. Quantum Electron. 12, 33-39 (2006).

• [29]

D.V. Skryabin, and A. V. Gorbach, “Colloquium: Looking at a soliton through the prism of optical supercontinuum,” Rev. Mod. Phys. 82, 1287-1299 (2010).

• [30]

D.V. Skryabin, and Y. V. Kartashov, “Self-locking of the frequency comb repetition rate in microring resonators with higher order dispersions,” Opt. Express 25, 27442-27451 (2017).