Broadband Resonator-Waveguide Coupling for Efficient Extraction of Octave Spanning Microcombs

TL;DR

This paper presents a novel pulley coupling scheme for microresonators that significantly improves the extraction efficiency of octave-spanning microcombs, especially at short wavelengths, by optimizing waveguide-resonator coupling.

Contribution

It introduces a pulley waveguide design to enhance coupling uniformity across the octave, validated by coupled mode theory and experimental results.

Findings

20 dB improvement in short-wavelength extraction

Good agreement between theory and experiment

Enhanced microcomb extraction efficiency

Abstract

Frequency combs spanning over an octave have been successfully demonstrated in Kerr nonlinear microresonators on-chip. These micro-combs rely on both engineered dispersion, to enable generation of frequency components across the octave, and on engineered coupling, to efficiently extract the generated light into an access waveguide while maintaining a close to critically-coupled pump. The latter is challenging as the spatial overlap between the access waveguide and the ring modes decays with frequency. This leads to strong coupling variation across the octave, with poor extraction at short wavelengths. Here, we investigate how a waveguide wrapped around a portion of the resonator, in a pulley scheme, can improve extraction of octave-spanning microcombs, in particular at short wavelengths. We use coupled mode theory to predict the performance of the pulley couplers, and demonstrate good agreement with experimental measurements. Using an optimal pulley coupling design, we demonstrate a 20~dB improvement in extraction at short wavelengths compared to straight waveguide coupling.

Figures

Equations

$${\kappa }_{\mathrm{r}\to \mathrm{wg}}={\text{bigintss}}_L\Gamma (\omega ,l){\mathrm{e}}^{i\varphi }\mathrm{d}l$$

$$\Gamma (\omega ,l)=\frac{i\omega }{4}{\text{bigintss}}_S({\epsilon }_{\mathrm{w}g}-{\epsilon }_{\mathrm{R}}){\mathbf{E}}_{\mathrm{R}}^{\ast }\cdot {\mathbf{E}}_{\mathrm{wg}}\mathrm{d}r\mathrm{d}z$$

$$Q_{\mathrm{c}}=\omega \frac{n_{\mathrm{g}}^{\mathrm{R}}}{c_0}\frac{2\pi R}{{\left|{\kappa }_{\mathrm{r}\to wg}\right|}^2}$$

$${\kappa }_{\mathrm{r}\to wg}=\Gamma (\omega ){\text{bigintss}}_{L_{\mathrm{c}}}{\mathrm{e}}^{i\varphi }\mathrm{d}L=\Gamma L_{\mathrm{c}}\mathrm{sinc}\left(L_{\mathrm{c}}\sqrt{{\left(\text{nicefrac}\Delta \beta 2\right)}^2+{\Gamma }^2}\text{lower}0.4pt\text{vrule}height=9.30444pt,depth=-7.44359pt\right)$$

$$Q_{\mathrm{c}}=\omega \frac{n_{\mathrm{g}}^{\mathrm{R}}}{c_0}2\pi R{\left[\Gamma L_{\mathrm{c}}\mathrm{sinc}\left(L_{\mathrm{c}}\sqrt{{\left(\text{nicefrac}\Delta \beta 2\right)}^2+{\Gamma }^2}\text{lower}0.4pt\text{vrule}height=9.30444pt,depth=-7.44359pt\right)\right]}^{-2}$$

$$\Delta n_{\mathrm{eff}}^{\mathrm{wg}}=2\frac{c_0}{\omega }\sqrt{{\left(\text{nicefrac}k\pi \widetilde{L_{\mathrm{c}}}\right)}^2-{\Gamma }^2}\text{lower}0.4pt\text{vrule}height=9.30444pt,depth=-7.44359pt;k\in \mathbb{N}.$$

Full text

††journal: ol

Kerr solitons generated in nonlinear microresonators [1] are promising for many applications in telecommunications [2], range measurement [3], and optical frequency metrology [4]. However, for frequency metrology in particular, many applications require octave-spanning bandwidth for full stabilization through $f$-2$f$ technique [5]. Suitable engineering of the resonator dispersion profile for octave bandwidth [6, 7], or even super-octave bandwidth [8], has been widely reported and octave-spanning soliton frequency combs have been demonstrated [9, 10], along with $f$-2$f$ stabilization [4, 11]. Such stabilization requires sufficient power at the frequencies of interest. This ultimately depends not only on the generated intracavity field and ability to take advantage of effects like targeted dispersive wave (DW) emission [7, 12, 13], but also on the extraction of the intracavity field, usually through evanescent coupling to an in-plane waveguide (or waveguides) for microring resonators. Efficient extraction over an octave of bandwidth is particularly non-trivial due to the wavelength dependence of both the phase matching and the spatial mode overlap between the resonator and waveguide modes.

In this letter, we characterize an approach to overcome this challenge, particularly at short wavelengths, based on a pulley configuration in which a portion of the access waveguide is wrapped around the microring [14]. Though utilized in our recent octave-spanning microcomb works [9, 4, 11], this approach was not studied in detail. Here, we present a basic coupled mode theory (CMT) formalism to design the pulley to improve resonator-waveguide coupling, thus comb extraction at short wavelengths, while maintaining desirable coupling in the pump and long wavelength bands. One consequence of pulley coupling is the introduction of narrow spectral windows in which essentially no coupling occurs, due to complete phase-mismatch between the ring and waveguide modes. Consequently, it is important to control the spectral position of these windows in which no coupling occurs, which we refer to as anti-phase-matched frequencies, so that they are separated from the regions of interest, namely the pump and DW frequencies. Experimentally, we validate both the control of the pulley anti-phase-matched frequencies and the enhancement of short wavelength extraction, by $\approx$ 20 dB relative to conventional point coupling using straight waveguides.

A number of computational approaches have been used to model coupling between ring resonators and waveguides [15, 14, 16]. Here, we model resonator-waveguide coupling in an integrated planar geometry by considering only the region over which their fields interact (Fig. 1(a)), with microring outer radius $R$ and ring width $RW$ separated by a gap $G$ from the coupling waveguide of width $W$. Using the spatial CMT formalism common to waveguide directional couplers [17, 18], we determine the per-pass coupling coefficient from resonator to waveguide $\kappa_{\mathrm{r\rightarrow wg}}$, based on the overlap of the ring and waveguide fields, integrated over the coupling portion. From there, the coupling quality factor $Q_{\mathrm{c}}$ is obtained, and compared to a typical resonator intrinsic quality factor ($Q_{\mathrm{i}}$) to gauge whether the coupling is at an appropriate level, namely close to critical coupling at the pump $Q_{\mathrm{c}}\approx Q_{\mathrm{i}}$.

The coupling coefficient between the ring and the waveguide is:

[TABLE]

with $L$ being the optical path and $\Gamma$ the overlap of the ring mode projected onto the waveguide mode as:

[TABLE]

with $\omega$ the angular frequency, $\bm{E_{\mathrm{R,wg}}}$ the electric field of ring and waveguide mode respectively, normalized such that $P=\frac{1}{2}\iint\left(\bm{E}\times\bm{H^{*}}\right)\cdot\hat{\theta}\mathrm{d}r\mathrm{d}z=1$ [18], with $r$, $\theta$, and $z$ being the radial, azimuthal, and vertical directions as taken from the center of the ring (Fig. 1(a)). $\varepsilon_{\mathrm{R},wg}$ is the dielectric permittivity considering only the ring and waveguide, respectively. The accumulated phase term in Eq. 1 corresponds to \phi=l\mathchoice{{\hbox{\displaystyle\sqrt{\left(\nicefrac{{\Delta\beta}}{{2}}\right)^{2}+\Gamma^{2},}}\lower 0.4pt\hbox{\vrule height=9.30444pt,depth=-7.44359pt}}}{{\hbox{\textstyle\sqrt{\left(\nicefrac{{\Delta\beta}}{{2}}\right)^{2}+\Gamma^{2},}}\lower 0.4pt\hbox{\vrule height=9.30444pt,depth=-7.44359pt}}}{{\hbox{\scriptstyle\sqrt{\left(\nicefrac{{\Delta\beta}}{{2}}\right)^{2}+\Gamma^{2},}}\lower 0.4pt\hbox{\vrule height=6.53888pt,depth=-5.23112pt}}}{{\hbox{\scriptscriptstyle\sqrt{\left(\nicefrac{{\Delta\beta}}{{2}}\right)^{2}+\Gamma^{2},}}\lower 0.4pt\hbox{\vrule height=5.03888pt,depth=-4.03113pt}}}, where $\Delta\beta=\frac{m}{R_{\mathrm{wg}}}-n_{\mathrm{eff}}^{\mathrm{wg}}\frac{\omega}{c_{\mathrm{0}}}=\frac{\omega}{c_{\mathrm{0}}}\Delta n_{\mathrm{eff}}^{\mathrm{wg}}$ is the difference of propagation constant between the ring and the waveguide mode, within the waveguide, with $m$ the azimuthal mode number of the ring for a given resonance frequency, $R_{\mathrm{wg}}=R+G+\nicefrac{{W}}{{2}}$, and $c_{0}$ the speed of light in vacuum.

The effects of the coupling coefficient, $\kappa_{\mathrm{r\rightarrow wg}}$, the phase-mismatch $\Delta\beta$ and the mode overlap $\Gamma$ can be combined into a single quantity describing the coupling strength known as the coupling quality factor $Q_{\mathrm{c}}$, defined as:

[TABLE]

where $n^{\mathrm{R}}_{\mathrm{g}}$ is the group index of the ring resonator. It is convenient to compare $Q_{\mathrm{c}}$ with the intrisic quality factor $Q_{\mathrm{i}}$, and to define the extraction efficiency as $\eta=\left(1+\nicefrac{{Q_{\mathrm{c}}}}{{Q_{\mathrm{i}}}}\right)^{-1}$.

The basic challenge that we address is conceptually illustrated in Fig. 1, and is easy to explain using this CMT framework. Straight waveguide coupling to a ring resonator involves a limited interaction length over which the waveguide and resonator modes spatially overlap, leading to close to a point-like coupling region, particularly for small diameter rings. At long wavelengths (low frequencies) and for a carefully chosen gap size $G$, the overlap $\Gamma$ can be appreciable enough that a short interaction length is adequate, yielding $Q_{\textrm{c}}$ comparable to $Q_{\textrm{i}}$ (i.e., critical coupling). However, as seen in Fig. 1(a,b), as the wavelength decreases (frequency increases), each mode is more confined, leading to a reduction in $\Gamma$, and $Q_{\textrm{c}}$ increases exponentially with frequency. This results in poor coupling at short wavelengths, with $Q_{\textrm{c}}$ orders of magnitude higher than $Q_{\textrm{i}}$ (Fig. 1(b)). This is problematic for octave-spanning combs, as illustrated in Fig. 1(c). Here, the spectrum of the intracavity field is simulated by solving the Lugiato-Lefever equation with the open-source pyLLE package [20], for a geometry appropriate for supporting octave-span operation. Though the dispersion has been engineered to support nearly harmonic dispersive waves at 280 THz and 155 THz, the $>$100 $\times$ difference in $Q_{\textrm{c}}$ will lead to very different out-coupled powers (given a $Q_{\textrm{i}}\leavevmode\nobreak\ \approx\leavevmode\nobreak\ 3\times 10^{6}$ that is not expected to significantly vary with wavelength), as seen in the plot of $\eta$ in Fig. 1(c). This will be a major impediment to direct self-referencing.

To overcome this issue, it is possible to increase the interaction length between the waveguide and the ring by wrapping the former around the latter, resulting in a pulley coupling design [14] shown schematically in Fig. 2(c). We note that the overlap coefficient $\Gamma$ in Eq. 2 is independent of position along the optical path $L$, and the accumulated phase has to be accounted for across the pulley length $L_{\mathrm{c}}$, length for which the gap is constant between the ring and the waveguide. Thus Eq. 1 becomes:

[TABLE]

with $\mathrm{sinc}(x)=\nicefrac{{\mathrm{sin}(x)}}{{x}}$. Hence Eq. 3 can be rewritten as:

[TABLE]

The above only accounts for the region where the resonator-waveguide gap is constant, and not where the waveguide bends towards and away from the ring, as seen in Fig. 2(c). To account for this, we evaluate Eq. 1 in these regions, resulting in a coupling coefficient $\kappa_{1}$. The ratio $\nicefrac{{\kappa_{1}(\omega)}}{{\Gamma_{\mathrm{pulley}}(\omega)}}$ gives the effective length of the curved portion (Fig. 2(c)). Hence, one can then introduce an effective pulley length $\widetilde{L_{\mathrm{c}}}(\omega)=L_{\mathrm{c}}+2\nicefrac{{\kappa_{1}(\omega)}}{{\Gamma_{\mathrm{pulley}}(\omega)}}$ that replaces the pulley length in Eq. 4.

Interestingly, Eqs. 4 and 5 suggest that resonances in the coupling will happen according to:

[TABLE]

Physically, these resonances correspond to locations where the access waveguide and the ring waveguide are anti-phase-matched, and are not observed for a straight waveguide due to the limited interaction length (over which the gap is continuously varying). To investigate further, we calculate $Q_{\mathrm{c}}$ using parameters that correspond to the experimental system studied, that is, 780 nm thick silicon nitride (Si3N4) microrings that are symmetrically clad in silica (SiO2), with $R=$23.3 $\mu$m. We pick $RW=1600$ nm (resulting in the simulated frequency comb shown in Fig. 1(b)), along with $W=550$ nm and $G=700$ nm. As shown in Fig. 2, the anti-phase-matching condition results in sharp peaks in $Q_{\mathrm{c}}$ for frequencies that vary with $L_{c}$. At these frequencies, regardless of the overlap between the ring and the waveguide modes, no transfer of energy occurs. The behavior of $Q_{\mathrm{c}}$ on either side of the resonances is important for octave-spanning comb applications. On the blue side (short wavelengths), the overall increase in interaction length results in smaller $Q_{\mathrm{c}}$ (improved coupling) than in the straight waveguide case. On the red side (longer wavelengths), the difference in $\Delta n_{\mathrm{eff}}^{\mathrm{wg}}$, accumulated over $L_{\mathrm{c}}$, keeps $Q_{\mathrm{c}}$ larger than in the straight waveguide case, where the rings are generally overcoupled. The net result is a reduced wavelength-dependence in $Q_{\textrm{c}}$ (outside of the anti-phase-matched window) than for a straight waveguide. We also note that when the pulley is sufficiently long, higher orders of anti-phase-matching can be satisfied (i.e., $k>1$), leading to multiple resonances in $Q_{\mathrm{c}}$.

To validate the CMT modeling, we first verify the pulley resonance behavior through linear transmission measurements of devices designed to show a $Q_{\mathrm{c}}$ resonance within the $182$ THz to $207$ THz tuning range of our laser source. We keep the pulley parameters fixed, namely gap $G=800$ nm, waveguide width $W=750$ nm, and pulley coupling length $L_{\mathrm{c}}=40$ $\mu$m, while the ring width $RW$ is varied. This results in variation of the effective index of the microring $n_{\mathrm{eff}}^{\mathrm{ring}}$, leading to a modification of the anti-phase-matching condition and hence the spectral position of the corresponding frequencies. By fitting approximately 25 resonances of the first order transverse electric (TE) mode family that appear within the laser scan range, we extract the spectral dependence of $Q_{\mathrm{c}}$ (Fig. 3(a)) for each $RW$, taking into account internal losses, coupling, and backscattering [19].

Simulations of $Q_{\textrm{c}}$ through CMT match both the values and the trend of the experimental $Q_{\textrm{c}}$, including the divergence at resonance. Moreover, one can reproduce the variation in $Q_{\textrm{c}}$ anti-phase-matching spectral position with $RW$ due to the change of effective index. The simulations also show the difference of sensitivity in dimension between the ring and the waveguide. As the waveguide is narrower with a mode less confined than the ring, a small variation of its width results in a significant change in its effective index. Hence, by only changing the waveguide width by 10 nm, the anti-phase matching frequency shifts by about 5 THz. To achieve the same shift, one needs to modify the $RW$ by 50 nm. This gives the ability to tune the position of the pulley anti-phase-matched frequency while keeping the microring resonator at a fixed geometrical dimension that is likely already dispersion-optimized.

Outside of the tunable laser range, it is possible to extract the position of the pulley anti-phase-matched frequency by measuring the spectra of modulation instability (MI) combs generated through strong pumping ($P_{\mathrm{pmp}}=200$ mW) of the resonators at $\approx$ 1550 nm, on the blue-detuned side of a cavity resonance. The spectral components in these MI combs are not phase-locked, and the overall comb acts as a quasi-continuous-wave, spectrally broadband source. Thus, the pulley coupling can be studied using these states over a spectral range as broad as the comb bandwidth. To confirm this, we measure the MI comb spectra of the devices characterized linearly (Fig. 3(b)). We observe that the position of the anti-phase-matched frequency obtained through linear characterization of the device (by extracting $Q_{\textrm{c}}$) and through measuring the position of the dip in the MI comb are consistent. The latter method also agrees with the CMT simulations, and the position of the anti-phase-matched frequency is within the uncertainty in the fabricated geometry.

We now compare a pulley coupling design optimized for extraction of an octave-spanning microcomb, namely, with a coupling anti-phase-matched frequency in-between the pump and the short wavelength DW, against straight waveguide coupling for the same ring parameters. We first characterized $Q_{\textrm{c}}$ (through linear transmission measurements) in both the pump band and near the short wavelength DW, around 193 THz and 280 THz, respectively (Fig. 4(a)) using two continuous tunable laser centered around 1550 nm and 1050nm. We were unable to measure any resonance of the first order TE mode in the 280 THz range for the straight waveguide, as expected from simulations where $Q_{\textrm{c}}$ is orders of magnitude higher than the expected $Q_{\textrm{i}}\approx 3\times 10^{6}$ (as measured in other bands). In contrast, the pulley devices, for both $L_{c}$=15 $\mu$m and $L_{c}$=17 $\mu$m, exhibit a difference in $Q_{\textrm{c}}$ of only one order of magnitude between the two bands, and show good agreement with the values predicted by the CMT. Finally, from MI comb spectra (Fig. 4(b)), the advantage of using the pulley coupling approach for extraction is apparent. Pumping both the straight waveguide and the $L_{c}$=15 $\mu$m pulley devices such that the long DW and overall comb shape are the same, the pulley coupling shows a clear advantage in extracting the short DW with $>20$ dB increase in power obtained. This enhancement of short DW extraction has recently been applied in studies of octave-spanning soliton microcombs [4, 9, 11, 21].

In conclusion, we have presented a CMT formalism to design pulley couplers to help with the extraction of octave-spanning spectra from chip-integrated, microring-based frequency combs. We use the CMT to elucidate the roles of phase-mismatch and spatial overlap in the wavelength-dependent coupling spectrum. Finally, we show that using such pulley coupling increases by $\approx$20 dB the extraction of the short wavelength part of an octave-spanning frequency comb compared to the same resonator with a straight waveguide coupling.

Acknowledgement

This work is supported by the DARPA DODOS, ACES, and NIST-on-a-chip programs. G.M., X.L., Q.L. and A.R. acknowledge support under the Cooperative Research Agreement between UMD and NIST-PML, Award no. 70NANB10H193.