Abstract
We consider the interaction of counter-propagating waves in a bi-directionally pumped ring microresonator with Kerr nonlinearity. We introduce a hierarchy of the mode expansions and envelope functions evolving on different time scales set by the cavity linewidth and nonlinearity, dispersion, and repetition rate, and provide a detailed derivation of the corresponding hierarchy of the coupled mode and of the Lugiato-Lefever-like equations. An effect of the washout of the repetition rate frequencies from the equations governing the dynamics of the counter-propagating waves is elaborated in details.
Published by The Optical Society under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
1. Introduction
Microresonator frequency combs have been attracting a significant recent attention with their numerous practical applications and as an experimental setting to study fundamental physics of dissipative optical solitons, see [1–3] for recent reviews. So called Lugiato-Lefever (LL) model has become a paradigm in this research area [3–6]. Its soliton solutions have some pre- and post- Lugiato-Lefever history in and outside the optics context, see, e.g., [7–14]. However, the area has exploded after a breakthrough experimental demonstration of Ref. [5]. In terms of the first principle approach to the Kerr microresonator model development, the decade old work [6] has remained a main reference. However, together with experimental progress in the area of Kerr microresonators, the underpinning theory deserves a refreshed outlook. One of the recent challenges has emerged after a series of experiments with birectionally pumped microresonators, where combs and solitons have been observed in counter-propagating waves [15–19]. Bi-directionally pumped and, related dual-ring, microresonators have also been recently studied for symmetry breaking [20–24] and gyroscope [25–30] related effects, including idealised PT-symmetric cases [31–33].
A variety of models has been reported in the context of experiments dealing with a single mode operation in each direction [20–24,34]. We note, here that studies into single mode bidirectional lasers, laser gyroscopes and symmetry breaking in them have history going back to 1980’s, see, e.g., [35–37]. To interpret recent soliton experiments, [15,18] have used models without nonlinear cross-coupling, while [16] has accounted for it. As we will see below, neglecting by the nonlinear cross-coupling was probably a better approach to analyse the experimental measurements under the circumstances, when modelling in neither of [15,16,18] included the effect of opposing group velocities, i.e., opposite signs of the resonator repetition rates for counter-propagating waves.
Due to complexity of the problem and diversity of equations both met in literature and the ones that are encountered during first principle analysis of the problem, it appears beneficial to have a detailed reference derivation that can be followed and tailored by a reader. Such mathematically transparent and physically motivated derivation that can be readily mapped onto a variety of experimental setups is present below. Focus of our work is to identify a hierarchy of the mode expansions and envelope functions evolving on different time scales set by the cavity linewidth, nonlinearity and 2nd order dispersion (slow time scales), and by the repetition rate (fast time scale), which can be used to derive a hierarchy of the coupled mode and envelope equations. We pay particular attention to comprehensive explanations of our derivation steps and interpretation of the results.
2. Hierarchies of mode expansions and envelope functions
This Section introduces physical system and discusses a hierarchy of mode amplitudes and envelope functions accounting for different time-scales. It also outlines plan of work for the rest of the paper.
Maxwell equations written for the electric field components ${{\cal E}}_\alpha$ using Einstein’s notations read as
Electric field vector ${{\cal E}}_\alpha$ inside a ring resonator is expressed as a superposition of its linear modes $F_{\alpha j}(r,z)e^{ij\theta \pm i\omega _jt}$, which are solutions of Eq. (1) with ${{\cal N}}_\alpha =0$:
In order to cut notational complexity and drop the $\alpha$ index, we consider TM family, so that from now on $F_{z j}\to F_j$, and
Results of our derivations would be the same for TE modes, $F_{r j}\to F_j$, and therefore, what we are loosing is only formal consideration of the nonlinear coupling between the TE and TM families.We assume that inhomogeneities of the resonator surfaces result in scattering in general and in backscattering, in particular, and hence lead to the linear coupling between the modes. We account for these effects assuming
Here $\varepsilon _{id}$ is the dispersive dielectric function of the ideal (no backscattering) geometry, that does not depend on $\theta$, while relatively small $\varepsilon _{in}(\theta )$ accounts for inhomogeneities along the ring. Mode profiles $F_{j}(r,z)$ are calculated for $\varepsilon _{in}=0$.${{\cal E}}$ is measured in V/m, hence normalising linear modes as $\max _{r,z} |F_{ j}|=1$ makes units of $b_jB^{\pm }_j$ to be V/m. Real field amplitude of a CCW mode is $2b_{j}|B^{+}_{j}|$, so that its intensity is $I^{+}_j=2c\epsilon _{vac}n_{j}b_{j}^{2}|B^{+}_{j}|^{2}$ and power is $I^{+}_j/S_j$. $S_j$ is the effective transverse mode area, $S_j={\big (}\iint |F_j^{\prime }|^{2} dxdz{\big )}^{2}/\iint |F_j^{\prime }|^{4} dxdz$ and $F_j^{\prime }=F_j(r,z)\vert _{y=0}$. $n_j$ is the linear refractive index, $n_j^{2}=\int _{-\infty }^{\infty }\varepsilon _{id}(\tau ,r=r_0,z=0)e^{i\omega _j\tau }d\tau$, where $r_0$ is the distance between the $z$ axis and a point of maximum of $|F_j|$. We define scaling factors $b_j$ as
so that the $|B^{\pm }_j|^{2}$ are measured in Watts. $\epsilon _{vac}$ is the vacuum susceptibility.We assume that the resonator is pumped into its $j_p$ mode and introduce the mode index offset $\mu =j-j_p$. The real field expression, Eq. (3), is then
$D_1$ is the resonator repetition rate (or free spectral range (FSR)) and $D_2$ is its group velocity dispersion. $D_2>0$ implies anomalous and $D_2<0$ normal dispersion. For example, the work [15] deals with a bi-directionally pumped silica ring with radius $1.5$mm and it has $D_1= 2\pi \times 22$GHz, $D_2=2\pi \times 16$kHz. The linewidth of this resonator is $\kappa =2\pi \times 1.5$MHz, and hence the corresponding finesse ${{\cal F}}=D_1/\kappa \simeq 13000$. The mode area estimate is $S_{j_p}\simeq 30\mu$m$^{2}$, which gives $b_j^{2}\simeq 4\times 10^{12}$V$^{2}$W$^{-1}$m$^{-2}$. Pump laser wavelength was $\simeq 1550$nm ($\omega _0\simeq 2\pi \times 193$THz) and the comb spectra observed there were relatively narrow and span over $\sim 20$nm bandwidth, corresponding to about $300$ modes, and the momentum of a mode nearest to the pump is estimated as $j_p=8700$.
In order to introduce a new set of mode amplitudes important in what follows, we transform Eq. (7) further:
In order to take control of $D_1$ in our future calculations, we define yet another set of slow amplitudes
Here $D_1$ is moved away from the exponential factors defining our third and final set of amplitudes $A_\mu ^{\pm }$. Instead, exponents with $D_1$ appear explicitly in the total field equation that uses $A_\mu ^{\pm }$,To summarize this section: $B_\mu ^{\pm }$ amplitudes absorb only the slowest time scales associated with the nonlinear effects and resonator losses. $A_\mu ^{\pm }$ absorb time scales associated with the second and higher order dispersions, in addition to the ones already inside $B_\mu ^{\pm }$. $Q_\mu ^{\pm }$ amplitudes evolve with the highest in our hierarchy frequency determined by the resonator repetition rate. To see how these different time scales and mode amplitudes are used to express the total field, ${{\cal E}}$, one should compare Eqs. (10a), (10b) and (15).
The rest of this work is structured as follows: In Section 3, we first derive a system of equations for $B^{\pm }_\mu$ and perform its exact reduction to the equations for $Q^{\pm }_\mu$. In Section 4, we come back to the equations for $B^{\pm }_\mu$, make the $D_1$ role explicit, eliminate the associated fast oscillations and derive a simpler system for $A^{\pm }_\mu$. Corresponding mean-field equations for the envelope functions $Q_{\pm }$ and $A_{\pm }$ and their counter parts with the reflected spatial coordinates are derived in Section 5.
3. Coupled mode equations
3.1 Separating equations for CW and CCW amplitudes
Substituting $t'=t-\tau$ in Eqs. (1), (3) we then assume that material response is fast so that $C_j(t-\tau )\simeq C_j(t)-\tau \partial _t C_j+\dots$. Neglecting all the 2nd and higher order time derivatives of $B^{\pm }_\mu$ we find that Eq. (1) transforms to
We now multiply the left and right hand-sides of Eq. (17) by $b_{j_p}F_{j_p}\exp ^{-ij_{\mu '}\theta }$, integrate in $r,z$ and $\theta$, and approximate $\omega _\mu \simeq \omega _0=\omega _{j_p}$, $n_\mu \simeq n_0$ inside all the pre-factors, but not in the powers of the exponents. The resulting model, see Eqs. (22), makes use of the two scattering matrices having dimensions of angular frequencies. One characterises scattering induced coupling between the co-propagating modes
Equation (22) can now be split, as per rotating wave approximation, into the parts proportional to $e^{\pm i\omega _{j_\mu }t}$ exponents, so that we have two equations defined on the slow, $D_2$ related, time scales:
In order to be used to describe laboratory experiments with microresonators, Eqs. (23) have to be amended with the single mode pump term and losses accounting for the finite linewidth. We take, for the laser frequency at the exact cavity resonance $\Omega =\omega _{\mu =0}=\omega _{j=j_p}$ and for the low pump levels, i.e., linear regime, the intracavity powers of CW and CCW waves to be $|{{\cal H}}_\pm |^{2}=|B^{\pm }_\mu |^{2}$. This is achieved via a phenomenological substitution
If pump is absent, then the field power would decay with the rate $\kappa$ (full width of the resonance). An expression linking ${{\cal H}}_\pm$ with the laser powers ${{\cal W}}_\pm$ is
where ${{\cal W}}_\pm$ are the laser powers pumping, respectively, CW and CCW waves. $\eta <1$ is the coupling efficiency via, e.g., a prism or a waveguide, into a resonator mode. $\eta =\kappa _{c}/\kappa$, where $\kappa _{c}$ is the coupling pump rate (equals coupling loss rate). ${{\cal F}}/\pi$ is the cavity induced power enhancement. Detailed theoretical and experimental studies of the power enhancement effect and coupling in and out considerations for ring cavities can be found in, e.g., [39,40].$\widehat {{\cal R}}_{\mu \mu }\sim 2\pi \times 4$ kHz in Ref. [15]. In this regime, it is safe to assume that $\kappa$ dominates over $\widehat \Gamma$ and $\widehat {{\cal R}}$ terms. Using this we disregard $\widehat \Gamma _{\mu '\mu }$ in what follows, and retain only the dominant diagonal terms in $\widehat {{\cal R}}_{\mu '\mu }$, i.e., $\widehat {{\cal R}}_{\mu '\mu \ne \mu '}\approx 0$. Dispersion of the diagonal terms is also disregarded, $\widehat {{\cal R}}_{\mu \mu }\simeq \widehat {{\cal R}}_{00}= R$. Accounting for all of the above and complex conjugating second of Eqs. (23) we conclude this subsection with
3.2 Opening up nonlinearity
andUsing Eqs. (11) to express amplitudes $B_\mu ^{\pm }$ via $Q_\mu ^{\pm }$ we find that all the time dependent exponents cancel out and the resulting coupled mode equations for $Q_\mu ^{\pm }$ amplitudes are
4. Washout of the repetition rate timescales from the coupled mode equations
Systems of Eqs. (30), (11), (12) on one side, and Eqs. (34), (12) on the other, are mathematically and physically equivalent. However, there are important observations to be made here. If one could assume that $|Q_+|^{2}+2|Q_-|^{2}$ under the integrals in the right hand sides of Eqs. (30) and Eqs. (34) is a slow function of time, then these integrals would be approximately equal to $Q_{\mu }^{\pm }e^{-i\delta _\mu t}$, see Eqs. (11). Balancing these with the $e^{i\delta _\mu t}$ exponents before the integrals in Eqs. (30), one would end up with equations involving time scales determined only by the linewidths, pump detuning and nonlinear resonance shifts, which are all order of MHz. MHz frequencies would be far simpler to resolve numerically, compare to GHz-THz frequencies associated with $D_1$, that are directly implicated inside $\delta _\mu$ in the linear parts of Eqs. (34).
In this Section, we demonstrate that there are both slow and fast time scales inside the nonlinear terms in Eqs. (30), and that the latter can be eliminated resulting in a simpler and better balanced system of equations for the $A_\mu ^{\pm }$ amplitudes, see Eqs. (14), (39).
We proceed by taking Eqs. (30a), express $Q_{\pm }$ via $B_\mu ^{\pm }$, see Eqs. (11), (12), calculate integrals in the nonlinear terms, see Eq. (35a), and perform the two step transformation, see Eqs. (35b), (35c),
Using $A_{\mu }^{\pm }$ amplitudes and detunings $\delta '_\mu$, which are both $D_1$ free, see Eqs. (14), allows to hide $e^{iD_2\mu ^{2}t/2}$ exponents in Eqs. (35). Adding the CCW equation, we have
5. Envelope models
Connection of the coupled mode equations to the wave dynamics becomes more intuitive, if one now derives the envelope, Lugiato-Lefever like, equations. First, we take the $Q_\mu ^{\pm }$ model, see Eqs. (34), and multiply Eq. (34a) with $e^{i\mu \theta }$ and (34b) with $e^{-i\mu \theta }$. We then sum up each of the equations in $\mu$ and use Eqs. (12), (13) connecting the envelopes $Q_{\pm }$ and the reflected envelopes $Q_{\pm }^{(r)}$ to their mode amplitudes. This procedure is free from approximations and it leads to a system of partial differential equations for $Q_\pm$ and $Q_\pm ^{(r)}$,
Starting from the equations for $A_\mu ^{\pm }$, see Eqs. (39), we follow a modified procedure. Namely, we multiply both CW and CCW equations by the same exponent $e^{i\mu \theta }$, use the envelope definitions in Eqs. (16), observe that $\int _0^{2\pi }|A_-|^{2}A_-e^{i\mu \theta }d\theta =\int ^{0}_{-2\pi }|A_-^{(r)}|^{2}A_-^{(r)}e^{-i\mu \theta }d\theta$ and, due to periodicity, $=\int _0^{2\pi }|A_-^{(r)}|^{2}A_-^{(r)}e^{-i\mu \theta }d\theta$, sum up in $\mu$, and derive the following envelope equations
Eqs. (40) (not Eqs. (41)) could in fact, be written without a rigorous derivation, by simply relying on common knowledge, let aside reflected envelopes in the backscattering terms. These equations include traditional cross-phase modulation, and also repetition rates terms and other odd order dispersion terms with the opposite signs. Contrary, Eqs. (41) have no repetition rate terms, i.e., $D_1$-terms, and the remaining odd dispersions, i.e., $D_3$, $D_5$, etc., come with the same signs. Simultaneously, phase sensitive nonlinear wave mixing effects induced by CW-CCW interaction have been washed out. The only nonlinear cross-interaction left comes from the integrated power, which merely shifts the detuning parameters. Thus, in the absence of backscattering a nonlinear bi-directional resonator operates as a uni-directional one, but with the detuning parameter altered by the total power of the counter-propagating wave.
6. Summary
We have derived coupled mode equations describing nonlinear wave mixing processes in Kerr microresonator with counter-propagating waves. Features of the first two coupled mode formulations given by Eqs. (30) and Eqs. (34) are that they fully account for the repetition rate effects and that nonlinear terms are taken in the real space, and can be evaluated via Fourier transforms, see also [42]. We then proceeded to present simplified multi-mode equations that neglect the repetition rate dynamics driving the phase sensitive terms responsible for nonlinear interaction between the counter-propagating fields (washout effect, Section 4), and again deal with the nonlinearity in the real space, see Eqs. (39).
Finally, we demonstrated that coupled mode equations ( 34) and Eqs. (39) are equivalent to two different, Lugiato-Lefever-like, envelope models. The one that involves the repetition rate dynamics, see Eqs. (40), links two usual envelopes for the CW and CCW fields, with two of their space reflections. While the one with the repetition rate averaged out, see Eqs. (41), makes a closed system already for two envelopes, $A_\pm$, one of which is reflected. We note, that $Q_\pm$ can be used directly to reconstruct total electric field, see Eq. (10c), while $A_\pm$ can not, but their respective mode amplitudes can, see Eqs. (15), (16a).
We have taken care to reveal all mathematical transformations, that allow a reader to verify our derivation steps and apply modifications if required. Opportunities for future theoretical and numerical studies offered by the models presented here are numerous, as well as their potential to guide and interpret experimental work.
Funding
Horizon 2020 Framework Programme (812818); Engineering and Physical Sciences Research Council (EP/R008159); Russian Science Foundation (17-12-01413).
Acknowledgement
Discussions with authors of [19] and with greatly missed M.L. Gorodetsky are gratefully acknowledged.
Disclosures
The author declares no conflicts of interest.
References
1. A. L. Gaeta, M. Lipson, and T. J. Kippenberg, “Photonic-chip-based frequency combs,” Nat. Photonics 13(3), 158–169 (2019). [CrossRef]
2. A. Pasquazi, M. Peccianti, L. Razzari, D. J. Moss, S. Coen, M. Erkintalo, Y. K. Chembo, T. Hansson, S. Wabnitz, P. Del’Haye, X. Xue, A. M. Weiner, and R. Morandotti, “Micro-combs: A novel generation of optical sources,” Phys. Rep. 729, 1–81 (2018). [CrossRef]
3. T. J. Kippenberg, A. L. Gaeta, M. Lipson, and M. Gorodetsky, “Dissipative Kerr solitons in optical microresonators,” Science 361(6402), eaan8083 (2018). [CrossRef]
4. L. A. Lugiato and R. Lefever, “Spatial Dissipative Structures in Passive Optical Systems,” Phys. Rev. Lett. 58(21), 2209–2211 (1987). [CrossRef]
5. 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]
6. Y. K. Chembo and N. Yu, “Modal expansion approach to optical-frequency-comb generation with monolithic whispering-gallery-mode resonators,” Phys. Rev. A 82(3), 033801 (2010). [CrossRef]
7. D. J. Kaup and A. C. Newell, “Theory of nonlinear oscillating dipolar excitations in one-dimensional condensates,” Phys. Rev. B 18(10), 5162–5167 (1978). [CrossRef]
8. D. W. Mc Laughlin, J. V. Moloney, and A. C. Newell, “Solitary waves as fixed points of infinite-dimensional maps in an optical bistable ring cavity,” Phys. Rev. Lett. 51(2), 75–78 (1983). [CrossRef]
9. K. Nozaki and N. Bekki, “Solitons as attractors of a forced dissipative nonlinear Schrödinger equation,” Phys. Lett. A 102(9), 383–386 (1984). [CrossRef]
10. I. V. Barashenkov, M. M. Bogdan, and T. Zhanlav, in Nonlinear World: XX International Workshop on Nonlinear and Turbulent Processes in Physics, V. G. Bariakhtar, ed. (World Scientific, 1990).
11. W. J. Firth and A. Lord, “Two-dimensional solitons in a Kerr cavity,” J. Mod. Opt. 43(5), 1071–1077 (1996). [CrossRef]
12. I. V. Barashenkov and Y. S. Smirnov, “Existence and stability chart for the ac-driven, damped nonlinear Schrödinger solitons,” Phys. Rev. E 54(5), 5707–5725 (1996). [CrossRef]
13. D. V. Skryabin, “Energy of internal modes of nonlinear waves and complex frequencies due to symmetry breaking,” Phys. Rev. E 64(5), 055601 (2001). [CrossRef]
14. D. V. Skryabin, “Energy of the soliton internal modes and broken symmetries in nonlinear optics,” J. Opt. Soc. Am. B 19(3), 529–536 (2002). [CrossRef]
15. Q. F. Yang, X. Yi, K. Y. Yang, and K. Vahala, “Counter-propagating solitons in microresonators,” Nat. Photon. 11(9), 560–564 (2017). [CrossRef]
16. S. Fujii, A. Hori, T. Kato, R. Suzuki, Y. Okabe, W. Yoshiki, A. C. Jinnai, and T. Tanabe, “Effect on Kerr comb generation in a clockwise and counter-clockwise mode coupled microcavity,” Opt. Express 25(23), 28969 (2017). [CrossRef]
17. C. Joshi, A. Klenner, Y. Okawachi, M. Yu, K. Luke, X. Ji, M. Lipson, and A. L. Gaeta, “Counter-rotating cavity solitons in a silicon nitride microresonator,” Opt. Lett. 43(3), 547 (2018). [CrossRef]
18. W. Weng, R. Bouchand, E. Lucas, and T. J. Kippenberg, “Polychromatic Cherenkov Radiation Induced Group Velocity Symmetry Breaking in Counterpropagating Dissipative Kerr Solitons,” Phys. Rev. Lett. 123(25), 253902 (2019). [CrossRef]
19. 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]
20. W. Yoshiki, A. Chen-Jinnai, T. Tetsumoto, and T. Tanabe, “Observation of energy oscillation between strongly-coupled counter-propagating ultra-high Q whispering gallery modes,” Opt. Express 23(24), 30851 (2015). [CrossRef]
21. Q. Cao, H. Wang, C. Dong, H. Jing, R. Liu, X. Chen, L. Ge, Q. Gong, and Y. Xiao, “Experimental Demonstration of Spontaneous Chirality in a Nonlinear Microresonator,” Phys. Rev. Lett. 118(3), 033901 (2017). [CrossRef]
22. M. T. M. Woodley, J. M. Silver, L. Hill, F. Copie, L. D. Bino, S. Zhang, G. Oppo, and P. Del’Haye, “Universal symmetry-breaking dynamics for the Kerr interaction of counterpropagating light in dielectric ring resonators,” Phys. Rev. A 98(5), 053863 (2018). [CrossRef]
23. L. D. Bino, J. M. Silver, S. L. Stebbings, and P. Del’Haye, “Symmetry Breaking of Counter-Propagating Light in a Nonlinear Resonator,” Sci. Rep. 7(1), 43142 (2017). [CrossRef]
24. J. M. Silver, K. T. V. Grattan, and P. Del’Haye, “Critical Dynamics of an Asymmetrically Bidirectionally Pumped Optical Microresonator,” arXiv:1912.08262 (2019).
25. J. Li, M.-G. Suh, and K. Vahala, “Microresonator Brillouin gyroscope,” Optica 4(3), 346–348 (2017). [CrossRef]
26. Y. H. Lai, Y. K. Lu, M. G. Suh, Z. Q. Yuan, and K. Vahala, “Observation of the exceptional-point-enhanced Sagnac effect,” Nature 576(7785), 65–69 (2019). [CrossRef]
27. A. B. Matsko and L. Maleki, “Bose-Hubbard hopping due to resonant Rayleigh scattering,” Opt. Lett. 42(22), 4764–4767 (2017). [CrossRef]
28. A. B. Matsko, W. Liang, A. A. Savchenkov, V. S. Ilchenko, and L. Maleki, “Fundamental limitations of sensitivity of whispering gallery mode gyroscopes,” Phys. Lett. A 382(33), 2289–2295 (2018). [CrossRef]
29. W. Liang, V. S. Ilchenko, A. A. Savchenkov, E. Dale, D. Eliyahu, A. B. Matsko, and L. Maleki, “Resonant microphotonic gyroscope,” Optica 4(1), 114–117 (2017). [CrossRef]
30. M. Dong and H. G. Winful, “Unified approach to cascaded stimulated Brillouin scattering and frequency-comb generation,” Phys. Rev. A 93(4), 043851 (2016). [CrossRef]
31. M. P. Hokmabadi, A. Schumer, D. N. Christodoulides, and M. Khajavikhan, “Non-Hermitian ring laser gyroscopes with enhanced Sagnac sensitivity,” Nature 576(7785), 70–74 (2019). [CrossRef]
32. J. Ren, H. Hodaei, G. Harari, A. U. Hassan, W. Chow, M. Soltani, D. Christodoulides, and M. Khajavikhan, “Ultrasensitive micro-scale parity-time-symmetric ring laser gyroscope,” Opt. Lett. 42(8), 1556–1559 (2017). [CrossRef]
33. C. Milián, Y. V. Kartashov, D. V. Skryabin, and L. Torner, “Cavity solitons in a microring dimer with gain and loss,” Opt. Lett. 43(5), 979–982 (2018). [CrossRef]
34. G. Lin, S. Diallo, J. M. Dudley, and Y. K. Chembo, “Universal nonlinear scattering in ultra-high Q whispering gallery-mode resonators,” Opt. Express 24(13), 14880–14894 (2016). [CrossRef]
35. W. R. Christian and L. Mandel, “Frequency dependence of a ring laser with backscattering,” Phys. Rev. A 34(5), 3932–3939 (1986). [CrossRef]
36. R. J. C. Spreeuw, R. Centeno Neelen, N. J. van Druten, E. R. Eliel, and J. P. Woerdman, “Mode coupling in a He-Ne ring laser with backscattering,” Phys. Rev. A 42(7), 4315–4324 (1990). [CrossRef]
37. D. V. Skryabin, A. G. Vladimirov, and A. M. Radin, “Spontaneous phase symmetry breaking due to cavity detuning in a class-A bidirectional ring laser,” Opt. Commun. 116(1-3), 109–115 (1995). [CrossRef]
38. A. R. Aftabizadeh, Y. K. Huang, and J. Wiener, “Bounded solutions for differential equations with reflection of the argument,” J. Math. Anal. Appl. 135(1), 31–37 (1988). [CrossRef]
39. D. Romanini, I. Ventrillard, G. Méjean, J. Morville, and E. Kerstel, “Introduction to Cavity Enhanced Absorption Spectroscopy,” in Cavity Enhanced Spectroscopy and Sensing, G. Gagliardi and H.-P. Loock, eds. (Springer, 2013).
40. I. Pupeza, Power Scaling of Enhancement Cavities for Nonlinear Optics (Springer, 2012).
41. D. C. Cole, A. Gatti, S. B. Papp, F. Prati, and L. Lugiato, “Theory of Kerr frequency combs in Fabry-Perot resonators,” Phys. Rev. A 98(1), 013831 (2018). [CrossRef]
42. T. Hansson, D. Modotto, and S. Wabnitz, “On the numerical simulation of Kerr frequency combs using coupled mode equations,” Opt. Commun. 312, 134–136 (2014). [CrossRef]