• Advanced Photonics
  • Vol. 2, Issue 4, 046001 (2020)
Jianhao Zhang1、*, Vincent Pelgrin1, Carlos Alonso-Ramos1, Laurent Vivien1, Sailing He2, and Eric Cassan1、*
Author Affiliations
  • 1Université Paris-Saclay, CNRS, Centre de Nanosciences et de Nanotechnologies, Palaiseau, France
  • 2Zhejiang University, Centre for Optical and Electromagnetic Research, State Key Laboratory for Modern Optical Instrumentation, Hangzhou, China
  • show less
    DOI: 10.1117/1.AP.2.4.046001 Cite this Article Set citation alerts
    Jianhao Zhang, Vincent Pelgrin, Carlos Alonso-Ramos, Laurent Vivien, Sailing He, Eric Cassan. Stretching the spectra of Kerr frequency combs with self-adaptive boundary silicon waveguides[J]. Advanced Photonics, 2020, 2(4): 046001 Copy Citation Text show less

    Abstract

    Dispersion engineering of optical waveguides is among the most important steps in enabling the realization of Kerr optical frequency combs. A recurring problem is the limited bandwidth in which the nonlinear phase matching condition is satisfied, due to the dispersion of the waveguide. This limitation is particularly stringent in high-index-contrast technologies such as silicon-on-insulator. We propose a general approach to stretch the bandwidth of Kerr frequency combs based on subwavelength engineering of single-mode waveguides with self-adaptive boundaries. The wideband flattened dispersion operation comes from the particular property of the waveguide optical mode that automatically self-adapts its spatial profile at different wavelengths to slightly different effective spatial spans determined by its effective index values. This flattened dispersion relies on the squeezing of small normal-dispersion regions between two anomalous spectral zones, which enables it to achieve two Cherenkov radiation points and substantially broaden the comb, achieving a bandwidth between 2.2 and 3.4 μm wavelength. This strategy opens up a design space for trimming the spectra of Kerr frequency combs using high-index-contrast platforms and can provide benefits to various nonlinear applications in which the manipulation of energy spacing and phase matching are pivotal.

    1 Introduction

    Third-order nonlinear optical processes in photonic waveguides have raised increasing interest in recent years, due to their unique capabilities for on-chip light generation and for the manipulation of guided light in spectral/temporal domains, with an immense potential for the implementation of light sources based on optical parametric amplification, supercontinuum, or frequency comb generation processes.19 Kerr frequency comb generation has great potential for applications in a wide range of fields, from metrology-on-a-chip to the synthesis of radiofrequency signals or integrated spectroscopy.1019 Comprising discrete and equally spaced frequency lines, the Kerr frequency combs require precise control of the waveguide dispersion and nonlinearity, the gain and loss of optical waveguides, especially when temporal patterns with few solitons are expected.20

    Controlling the chromatic dispersion of optical waveguides is a key preliminary step for the exploitation of Kerr frequency combs. As such, the problem is quite simple in principle: it is essentially a question of optimizing the geometrical dimensions of the waveguides in order to compensate for the material natural dispersion as well as that provided by the nonlinear effect itself, in order to satisfy, as a whole, the energy and momentum conservation conditions. To balance the nonlinearity-induced wavevector phase mismatch, anomalous dispersion is generally envisaged;10,11,20 it can be supported directly from the materials (e.g., silica at telecom wavelengths) or induced by the waveguide dispersion with well-designed waveguide cross-sections. As a result, toroidal-shape cavities using silica10 or MgF212,13 are frequently used for frequency comb generation due to their ultrahigh Q factors of up to a few millions and to the low dispersion in these structures. Another classical material platform used for frequency comb demonstrations is based on silicon nitride (SiN) waveguides. Due to the high quality-factors of SiN microring resonators (up to 106) and a nonlinear Kerr index around 2.4×1019  m2/W, high-performance frequency combs extending two octaves can be achieved on-chip.1518 Despite the merits and inherent advantages of this platform, the SiN waveguides are less attractive for compact photonic integration than silicon ones. The realization of Kerr frequency combs directly in Si waveguides thus presents great interest. However, achieving broadband phase-matching in Si waveguides is difficult due to the high-index contrast between Si and SiO2 (Δn2) that makes the dispersion engineering condition highly sensitive to wavelength.5,21 It appears that, by moving from silicon nitride to silicon waveguides, an advantage (a lower pump power) and a drawback (a reduced frequency comb spectral width) result simultaneously. Consequently, addressing the problem of the spectral bandwidth of microcombs in silicon waveguides is an important point. Broadband phase-matching was shown based on the optimization of high order dispersion terms, e.g., by relying on fourth-order dispersive waveguides,2226 or by implementing slightly etched rib geometries.6 Yet, the proposed solutions require complex fabrication processes, with deposition of different materials or a tight control of the rib and slab thicknesses, respectively. Though mode-locked pulse generation has been achieved using normal-dispersion,27 a flexible method to flatten the silicon waveguide dispersion and trim the spectra for Kerr frequency comb generation is still expected. The feasibility of spectrum engineering28 could further enhance optical frequency combs as strong candidates for the realization of on-chip silicon light sources29 and could provide a strong basis for many applications in on-chip spectroscopy or metrology,30 as well as for emerging research domains such as time-space-frequency mapping.31

    We propose here a method to trim the dispersion of strong index contrast waveguides using a self-adaptive boundary condition to significantly extend the bandwidth of Kerr comb spectra. This boundary condition was originally introduced in Ref. 32 for multimode optical waveguides in which an equivalent potential well3335 in the waveguide transverse direction was achieved by the use of a lateral index profile, enabling an automatic phase-matching condition of degenerate four-wave-mixing processes. We explore here the properties of such waveguides, now operating in the single mode regime, and investigate their use for frequency comb generation. Our study shows that, in fact, such waveguides have an immense potential in terms of chromatic dispersion engineering, which can be made minimal or slightly abnormal in a wide spectral band, while being controllable by the geometrical parameters of the structures. The interest lies in the general nature of the proposed method, which can be implemented on various photonic platforms and can also be applied to waveguides with high core-cladding index contrasts. The metric used in this study is the spectral bandwidth of the Kerr frequency combs achieved from this approach and the comparison with standard silicon waveguides.

    2 Comb Bandwidth of Silicon Step-Index Waveguide

    The generation of Kerr frequency combs is traditionally achieved in a Kerr nonlinear ring/disk resonator, and the overall process of comb generation has been discussed in previous works.12,3651 A narrowband continuous optical source (laser) is centered close to an optical resonance of the ring resonator. If the chromatic dispersion of the ring is well chosen, it can result in a cascade of FWM processes in the frequency domain, which leads to the generation of a frequency comb. From a temporal point of view, the generation of such a broad-spectrum comb corresponds to a regime of short solitonic pulses or even, in the best of cases, to a regime of monosoliton propagation in the ring. The full dynamics of the comb generation process in which a Kerr microring resonator is fed by a side-coupled waveguide can be quantitatively described by solving the nonlinear Schrödinger equation, in the form of the Lugiato–Lefever equation (LLE)3945,47,51where A(t,τ) is the pulse amplitude described by the time variable τ in a time frame corresponding to each round trip, t stands for the accumulated calculation time while τR is the roundtrip span, βk is the k’th derivative of the wavevector with respect to the angular frequency βk=kβωk (e.g., β2 is the group velocity dispersion), γ is the Kerr nonlinear parameter, α is the linear waveguide loss, β2PA, and β3PA stand for the two-photon5153 and three-photon absorption coefficients,54,55 respectively, and Aeff is the waveguide effective mode area. The free-carrier effects56,57 induced by multiphoton absorption are considered here through the αFCA coefficient, which is governed by the carrier density N(t,τ). The carrier density N(t,τ) is dynamically related to the intracavity power and the carrier lifetime tFC through τN(t,τ)=β2PA2ωAeff2|A(t,τ)|4+β3PA3ωAeff3|A(t,τ)|6N(t,τ)tFC. The pump is characterized by its amplitude Ain and δ0, which describes the phase detuning between the pump frequency ω0 and the closest n’th-order cold-cavity resonant frequency ωn of the ring resonator.

    For the purpose of making reasonable comparison with previously proposed architectures and sticking to the fact that a PIN structure is sometimes needed to obviate free carrier accumulation issues,14,53 a rib-shaped waveguide structure can be reasonably considered in a first approach, as shown in the inset of Fig. 1(a). The dispersion curve from Ref. 14 is reproduced in Fig. 1(a) (gray solid line), quantitatively through estimating D=2πcλ2·d2βdω2. Other compressed rib waveguides geometries (340-nm-thick silicon film, without SiN cladding) with rather similar dispersion properties but centered at smaller wavelengths are also chosen (orange and red curves), as shown in Fig. 1(a). Two configurations based on this thickness are shown in the inset, with only a difference in the waveguide width. Compared with a fully etched strip waveguide with the same or even a thinner thickness, these rib waveguides present a lower dispersion in the considered wavelength range. Four working points, including that of Ref. 14 (point 0), one for the 950-nm rib waveguide width configuration (point 1), and two for the 1050-nm wide rib case (points 2 and 3), are selected for a further investigation of the dispersion-governed comb spectrum properties.

    (a) Dispersion curves of silicon-on-insulator (SOI) waveguides with different dimensions, including the one (gray line) from Ref. 14. Inset: Cross-section of a silicon ring resonator used for Kerr frequency comb modeling. (b), (d) The wavevector mismatch and associated spectra of frequency combs, with the working points indicated in (a). The near-bottom spectra are the corresponding primary comb lines. (c) The wavelength-dependent Kerr coefficient and its interpretation in the time domain in the single-soliton comb state, indicated by smooth and modulated curves in black. The pattern of the temporal soliton is also given in blue, corresponding to the upper (τ) and right (power) axis.

    Figure 1.(a) Dispersion curves of silicon-on-insulator (SOI) waveguides with different dimensions, including the one (gray line) from Ref. 14. Inset: Cross-section of a silicon ring resonator used for Kerr frequency comb modeling. (b), (d) The wavevector mismatch and associated spectra of frequency combs, with the working points indicated in (a). The near-bottom spectra are the corresponding primary comb lines. (c) The wavelength-dependent Kerr coefficient and its interpretation in the time domain in the single-soliton comb state, indicated by smooth and modulated curves in black. The pattern of the temporal soliton is also given in blue, corresponding to the upper (τ) and right (power) axis.

    Due to the chromatic dispersion, light propagating along the ring resonator with different frequencies features different group velocities. The wavevector β at the frequency ω, will deviate from the wavevector of the pump (β0 at frequency ω0) by β0+β1(ωω0). Using Taylor’s expansion, this mismatch is noted as Δβ=[ββ0β1(ωω0)]=βn2(ωω0)2++βnn!(ωω0)n, which is considered an operator ik=2nβkk!(ikτk) in Eq. (1). This classical figure of merit for evaluating a Kerr comb bandwidth13,15 is usually called “integrated dispersion” and is expressed in units of the microring free spectral range number. We prefer here to consider it as the dispersion “wavevector mismatch” consistent with its native expression. In Fig. 1(b), we show the wavevector mismatch Δβ for the different waveguide configurations mentioned above. Points 0, 1, and 2 are very typical operating points where Δβ rapidly evolves from Δβ=0 with near-parabolic-shape local traces and increases back again with a second Δβ=0 intersection where Cherenkov radiations are generated15 [see also the secondary spectral peaks in Fig. 1(d) for several waveguide configurations close to the secondary dispersion cancellation points]. Compared with point 1, points 0 and 2 show smaller dispersion values (and possibly narrower anomalous regions) and therefore smaller deviations from the Δβ=0 condition, which leads to larger spectrum bandwidth, as shown in Fig. 1(d). Parameters used here have been directly extracted from Ref. 14, i.e., the pump power is fixed to 150 mW, a ring radius of 100  μm is considered, and a free-carrier lifetime of 10 ps is adopted, while propagation losses of 0.7  dB/cm are assumed, and critical coupling conditions are retained. The wavelength-dependent Kerr nonlinear coefficient46,21,28 was obtained using γ=2πn2ηλAneff·(ngnSi)2, where n2, η, and Aneff are the nonlinear index of silicon, the interaction overlap between the electric field and the silicon core, and the effective mode area, respectively. ng and nSi are the group index of the optical mode and the material index, respectively. The interpretation of the nonlinear coefficient in the time domain of the single-soliton state, labeled here by a blue line, is also shown in Fig. 1(c) in black. We clearly see that even though the nonlinear coefficient decreases rapidly with increasing wavelength due to an enlarged mode area, the nonlinear coefficient is more or less anchored to the value obtained at the pump wavelength, being modulated only near the soliton position where energy from all frequencies comes to make a contribution.

    One of the possible ways to expand the frequency comb spectrum is to shift the working point to a position closer to Δβ=0 on the same curve. Comparing point 3 to point 2, it can be observed that this action leads to the flattening of the wavevector mistmach Δβ [Fig. 1(b)] and therefore to a bandwidth improvement in Fig. 1(d). Configuration using working point 3 presents a comb bandwidth that is comparable (more precisely, a bit better) to that reported in Ref. 14 and is thus used hereafter as the reference case. The primary comb lines12,13,15 that are generated in the early stage of the comb generation process provide information that is useful for evaluating the bandwidth of parametric process;2426 it is also plotted for point 3 at the bottom of Fig. 1(d). This is considered by the phase matching condition with the nonlinearity-induced phase and the group velocity dispersion taken into account and is approximately described using23,24Here, Pintra is the intracavity power, while Δω is the frequency distance between the generated waves (idle/signal). These primary frequency lines are critical to the following subcombs and to the overall comb shape (and its bandwidth) around the pump line. It is thus very interesting to reduce and engineer the dispersion coefficient D globally while keeping the other parameters fixed (γ, Pin, etc.).2226 This matter is discussed in the next section of the article.

    3 Silicon Frequency Comb with Engineered Dispersion in Self-Adaptive Boundary Waveguides

    In a standard step-index single-mode waveguide like the one described in Fig. 2(a), several wavelengths actually have the same physical boundaries. When the light wavelength increases, a more and more significant fraction of the mode energy is pushed out of the waveguide core into the cladding material. To illustrate this point, in Fig. 2(b), we consider the dispersion curve of a 750  nm×300  nm SOI waveguide (single-mode for λ>2.1  μm). In view of this curve, it is quite remarkable that the variation of the chromatic dispersion of the mode remains moderate up to about 2.7  μm. Then, its evolution becomes very fast. This illustrates the notable point that will guide our strategy to optimize the geometry of optical waveguides for the generation of Kerr frequency combs: the chromatic dispersion of the waveguides is actually strongly related to the degree of spatial confinement of the guided mode. In situations in which the mode cannot extend laterally due to spatial confinement from fixed boundary conditions, its chromatic dispersion amplitude tends to increase. One possible way to limit the waveguide mode dispersion is therefore to ensure a smooth spatial lateral confinement and to minimize the effect of the lateral waveguide boundaries by allowing the electromagnetic guided mode to spatially spread as freely as possible. Investigating this kind of waveguide modes has led us to identify a so-called self-adaptive boundary (SAB) condition.32

    (a) Schematics of a step-index waveguide. (b) The chromatic dispersion parameter (D) of a silicon waveguide with a width and a height of 750 and 300 nm, respectively [n=3.48 for the waveguide core, with an air upper cladding and silica (n=1.44) for the bottom cladding]. (c) Schematics of a graded-index waveguide ensuring a self-adapting (SAB) condition (as defined in the text) and operating in a single-mode regime in the explored wavelength range. (d) Chromatic dispersion parameter D of the waveguide mode. The width and height of the waveguide are 750 and 340 nm, respectively.

    Figure 2.(a) Schematics of a step-index waveguide. (b) The chromatic dispersion parameter (D) of a silicon waveguide with a width and a height of 750 and 300 nm, respectively [n=3.48 for the waveguide core, with an air upper cladding and silica (n=1.44) for the bottom cladding]. (c) Schematics of a graded-index waveguide ensuring a self-adapting (SAB) condition (as defined in the text) and operating in a single-mode regime in the explored wavelength range. (d) Chromatic dispersion parameter D of the waveguide mode. The width and height of the waveguide are 750 and 340 nm, respectively.

    With this objective in mind, we decided to explore the use of SAB waveguides32 within the perspective of using their unusual dispersion properties to generate Kerr frequency combs. Figure 2(c) shows a typical lateral index profile of an SAB waveguide, where the longitudinal axis of the upper diagram and the intensity of the gray region of the lower diagram both represent the material index of the waveguide core. A graded-index profile along the transverse direction of the waveguide is used here, with ncent, nb, and nc representing the material index at waveguide center, at the waveguide inner edge, and the one of the cladding region. In a traditional strip waveguide (i.e., ncent=nb), the standing-wave region is confined by the boundary at the nb/nc interface since nb>neff. The standing-wave region is then the same for all wavelengths. In contrast, in a waveguide like the one represented in Fig. 2(c), if nb<neff, then the standing wave region is, instead, confined by the points where the material and the mode effective index values are equal. Since the mode effective index depends on the light wavelength, this leads to confining regions of different widths for different wavelengths, as illustrated in Fig. 2(c). In this situation, each confined field at a given wavelength automatically adapts its profile to a given spatial extension. To interrogate the effect of this SAB on the waveguide dispersion, the dispersion curves of a strip and two gradual SAB waveguides are put together for comparison in Fig. 2(d). Such graded-index profiles can be practically achieved by relying on the toolbox of subwavelength waveguide engineering,5860 as we mentioned in Ref. 32. As shown in the inset of Fig. 2(d), the waveguide is composed by a periodical corrugation in which the individual period is itself characterized by a unitary silicon filling factor at the center and a reduced one at the edge. Using such a typical patterned structure, an equivalent material index can be flexibly designed.32 In the precise study case reported in Fig. 2(d), the period, thickness, and width of the stair-case graded index waveguide were chosen as 150, 340, and 750 nm, respectively.

    The dispersion of the strip waveguide displays a trend that is similar to the one of Fig. 2(b), highlighting a rapid change at long wavelengths. When the SAB condition is applied, the dispersion in Fig. 2(d) can actually be separated in two regions: the left part (blue region) corresponding to the nb<neff condition and the right one (red region) for nb>neff. Dispersion of a strip waveguide and SAB waveguide with nb=1.5 and nb=2.8 is shown in Fig. 2(d) for comparison. The SAB allows us to trim the dispersion of the “long” wavelength where nb<neff and the wave is confined by the index contrast of nb/nc with the phase integral strongly dependent on the index nb. In contrast, since the short wavelength (nb>neff) is confined by the effective width of the waveguide where the electric field passes smoothly, the dispersion does not vary a lot. The self-adaptive response can be achieved in a continuously varying index profile, or in a discrete scheme, as shown in Fig. 3(a), where a multilevel silicon segment waveguide is described. This structure is, on the one hand, compatible with the liberation of the silicon waveguide from the silica substrate for eliminating the absorption of SiO2 for wavelengths higher than 3.6  μm; on the other hand, it provides room for the possible realization of a lateral PIN structure if needed for sweeping out free-carriers through a reverse applied bias, as shown in the inset.

    (a) Schematic of a graded-index waveguide ensuring a self-adapting (SAB) condition. Inset: Possible a configuration comprising doped regions for free-carrier extraction. (b) Chromatic dispersion parameter D of the four possible configurations of the proposed bilevel geometry. The silicon height is t=340 nm. (c) The calculated wavevector mismatch and (d) the single-soliton comb spectrum for our proposed SAB waveguide, at the working point in (b), compared with the reference rib microring described in Fig. 1 (red color). The two predicted Cherenkov radiation points are marked by the dashed circles in (c).

    Figure 3.(a) Schematic of a graded-index waveguide ensuring a self-adapting (SAB) condition. Inset: Possible a configuration comprising doped regions for free-carrier extraction. (b) Chromatic dispersion parameter D of the four possible configurations of the proposed bilevel geometry. The silicon height is t=340  nm. (c) The calculated wavevector mismatch and (d) the single-soliton comb spectrum for our proposed SAB waveguide, at the working point in (b), compared with the reference rib microring described in Fig. 1 (red color). The two predicted Cherenkov radiation points are marked by the dashed circles in (c).

    To further expand the previous ideas, we now illustrate the proposed principle with a few typical examples chosen for their applicative relevance. The silicon thickness is chosen first at 340 nm as proposed in the reference structure in Fig. 1, which represents an appropriate value for optimizing the waveguide dispersion in the 2 to 4  μm wavelength range. The length and width of the nanoarms are fixed at Larm=50  nm and warm=2.5  μm, which are values that are compatible with deep ultraviolet immersion lithography fabrication techniques. Let us first consider bilevel cases in which the value of w2 in Fig. 3(a) is zero. Four configurations are plotted for comparison. The corresponding parameters are listed in the figure and figure caption. As the SAB condition can flatten the mode anomalous dispersion level for wavelengths not confined by the outer boundary, ripples containing two local peaks of anomalous dispersion can be expected using a bilevel SAB waveguide geometry, as observed in Fig. 3(b). The period of Bragg reflection is around 430 nm, which is much larger than the longitudinal period that we consider (240 nm). Accordingly, this small period creates a bandgap only for wavelengths below 1.47  μm, which is far from our pump wavelength. This enables the periodic structure to operate as a homogeneous material5860 in which the frequency–wavevector curve is much less dispersive than those near the band edge. Despite the fact that photonic bandgap effects from this periodic structure do not spoil the predicted dispersion and the overall validity of the proposed design, previous studies61,62 have indicated that splitting of some resonant peaks can be observed under a certain grating phase condition. However, this high-order grating effect is much slighter than the first one and is not included in our discussion. In these cases, a moderate index contrast with a filling factor varying from 1 to 0.66 (i.e., L1/p=1) is chosen. As observed in Fig. 3(b), the whole wavelength range extending from 2.2 to 3.4  μm is greatly flattened. Compared with those classical strip/rib waveguides where the dispersion curve normally displays a parabolic shape, there is more than one local peak (e.g., one at 2.3  μm and one at 3.3  μm) existing in an anomalous region, which gives us a new degree of freedom for dispersion engineering purpose. In particular, a small normal-dispersion region (gray zone) is squeezed out between these two wide anomalous regions (highlighted by blue and purple colors, from left to right). We also emphasize that the equivalent index of the outermost subwavelength cladding is 1.84 compared with the core index 3.48, which promises a high index contrast and tight mode confinement. The dispersion distortion associated with the bending radius could be seamlessly addressed by increasing the ring radius or by engineering the filling factor of the waveguide.

    Two working points are selected for the calculation of the wavevector mismatch Δβ, as shown in Fig. 3(c). The curve with point 1 shows a typical parabolic shape due to the high dispersion value and the rapid change of dispersion. However, by pushing it to point 2, two intersection points (therefore two Cherenkov radiation regions) with Δβ=0 are now achieved; they are highlighted by the black dashed circles. The consecutiveness of these two intersections, only spaced by 300  nm, forces the whole curve to be flattened. Consequently, even better results than the reference one (red curve) are then obtained. Though fourth-order dispersion15,21 and even higher-order dispersion63 have been adopted for extending the zero-GVD working point, an improvement in comb bandwidth has not yet been witnessed. Flat dispersion has been achieved in a similar way in a vertical-slot waveguide,64 but the heterogeneous implementation was less attractive compared with our monolithic approach. The corresponding single-soliton comb spectra obtained by solving the LLE are given in Fig. 3(d). Since experimental demonstration6568 already indicated loss levels smaller than 3  dB/cm and according to previous experiments carried out on similar sub-wavelength waveguide structures, a linear waveguide loss level of 2  dB/cm was considered here. The considered radius is the same as the one of the reference rib waveguide, i.e., 100  μm. To overcome the effect of the loss and enable the generation of a comb state, the pump power is slightly increased to 200 mW. Note that critical coupling is always considered. The frequency detuning for achieving the final single-soliton comb is 1.513 GHz from the initialization. In all of the cases considered here, the two-photon-induced and three-photon-induced absorption lead to extra loss of not more than 0.001 and 0.01  dB/cm, in which the additional loss penalty is much smaller than the linear loss level itself. The free carrier lifetime considered in this case is 50 ps. From our model, we estimated the influence of a free-carrier lifetime ranging from 10 to 100 ps and observed, in the considered conditions, a weak influence of its value on the comb dynamics and the required pumping power. In a complementary fashion, we noticed that free-carriers exert a strong influence on the whole process only when their lifetime is larger than 200 ps, in which case the carrier density can reach a few 1016  cm3, which in return puts strong demand on the pump. Since the three-photon absorption process increases more quickly than the Kerr gain, at some point, comb generation is then totally prohibited for any pump power level. To sum up, we observed a lifetime cut-off with respect to each waveguide configuration. Regarding the main case corresponding to a propagation loss of 2  dB/cm, the cut-off lifetime was estimated to be around 220 ps. The situation could obviously be eased by lowering the linear loss. For instance, for propagation losses of 0.5  dB/cm, the cut-off lifetime increases to 420 ps.

    As shown in Fig. 3(b), the accumulated wavevector mismatch from the deep-blue anomalous region is soon compensated for by the gray normal-dispersion region, which gives the first Cherenkov radiation. The remaining part of this gray region, however, is compensated for again by the right anomalous region (purple), which leads to the generation of the second Cherenkov radiation. These results evidence the potential of SAB waveguides for flattening the overall waveguide dispersion profiles and forcing the wavevector mismatch to intersect the Δβ=0 axis twice without much deviation from it, which is difficult to achieve with traditional step-index waveguides of strip or rib shapes.

    To improve the performance, a trilevel case (w20) is considered to locally engineer the dispersion with an additional degree of freedom. As shown in Fig. 4(a), by varying the width of these three sections, we open a window of low-anomalous dispersion in a large wavelength range of up to 1000  nm, while small areas of each zone are achieved. The filling factor of Secs. 1 and 2 are fixed at L1/p=0.8 and L2/p=0.5. The frequency detuning for achieving the final single-soliton comb is 1.272 GHz from the initialization. Other parameters are the same as those previously mentioned, i.e., the loss, the carrier lifetime, the radius, and the pump power are 2  dB/cm, 50 ps, 100  μm, and 200 mW, respectively. As shown in Fig. 4(b), the wavevector mismatch Δβ is tremendously flattened near the axis from 2.4 to 3.2  μm wavelengths, with again two intersection points and small ripples across and near the zero Δβ axis. This flattened dispersion profile assists us in achieving, as shown in Fig. 4(c), a comb spectrum that gives a 20/30-dB bandwidth of 44/48 THz, which represents a great improvement compared with the optimal reference microring structure (red color). The two Cherenkov radiations can be distinguished from the smooth background. As visible, the primary comb lines of the SAB case (blue color), on the bottom part of Fig. 4(c), are much wider than that of the reference case. The evolution of the temporal and frequency patterns of the comb generation of the SAB waveguide comb with respect to point 2 is also presented in Figs. 4(d) and 4(e), respectively.

    (a) Chromatic dispersion parameter D of four bilevel configurations with different parameters and based on the subwavelength silicon waveguides shown in Fig. 3(a). The silicon height is t=340 nm. (b) and (c) The calculated wavevector mismatch and single-soliton comb spectra corresponding to the working point in (a) and that from Fig. 1 (red color/reference Si rib waveguide). The near-bottom spectra are the corresponding primary comb lines. (d) and (e) The corresponding comb generation process of the comb in deep blue. The two predicted Cherenkov radiation points are marked by the dashed circles or the white labels. The frequency detuning for achieving the final single-soliton comb is 1.272 GHz from the initialization.

    Figure 4.(a) Chromatic dispersion parameter D of four bilevel configurations with different parameters and based on the subwavelength silicon waveguides shown in Fig. 3(a). The silicon height is t=340  nm. (b) and (c) The calculated wavevector mismatch and single-soliton comb spectra corresponding to the working point in (a) and that from Fig. 1 (red color/reference Si rib waveguide). The near-bottom spectra are the corresponding primary comb lines. (d) and (e) The corresponding comb generation process of the comb in deep blue. The two predicted Cherenkov radiation points are marked by the dashed circles or the white labels. The frequency detuning for achieving the final single-soliton comb is 1.272 GHz from the initialization.

    As already stated before, the loss, the carrier lifetime, and the pump power of the reference comb case were fixed at 0.7  dB/cm, 10 ps, and 150 mW, respectively, and while raising the loss and carrier lifetime to 2  dB/cm and 50 ps in our SAB case, respectively, the required pump power is only raised to 200 mW. With this moderate pump power increase, the observed bandwidths have gained a large improvement, with a 20-dB bandwidth raised from 32 to 44 THz and a 30-dB bandwidth from 38 to 48 THz. Globally, these results prove that the nature of SAB waveguides implemented through a multisection lateral geometry opens up a new degree of freedom for engineering the needed waveguide dispersion and widening the Kerr frequency comb spectra.

    4 Conclusion

    We propose a general approach to trim and expand the spectral bandwidth of Kerr frequency combs in high-index-contrast optical platforms using self-adaptive boundary waveguides enabled by subwavelength index engineering. As a matter of illustration and a prime example, this approach enabled us to flatten the dispersion of silicon (on-insulator) waveguides and squeeze a small normal-dispersion region sandwiched between two anomalous zones. This improved dispersion shape enables us to achieve two Cherenkov radiation points and substantially broaden the bandwidth of comb spectra between 2.2 and 3.4  μm wavelengths. This strategy addresses the dispersion issue in high-index-contrast platforms and can be used to further explore versatile nonlinear applications in which the manipulation of energy spacing and phase matching is of primary importance. We believe that this new approach will have broad applications in integrated nonlinear photonics for the realization of frequency combs, supercontinuum sources, and other nonlinear effects applicable to signal processing, metrology, and spectroscopy on-a-chip.

    5 Appendix

    The dispersion of periodically segmented silicon waveguides and their mode profiles are obtained from the band analysis of FDTD simulation (Lumerical Inc.), while the comb generation is modeled with a self-built Python code. The Lugiato–Lefever equation is solved by applying a split-step Fourier integration scheme to Eq. (1), and multistep calculation is considered in every round trip. The influence of waveguide dispersion on the light pulse is considered in the Fourier frequency domain while the loss and parametric amplification are accounted for in the time domain. The initial power level in the cavity is set at 200  dBm. The pump wave is continuously injected into the input waveguide; it is first set at a wavelength slightly blue-detuned from a certain cold-cavity resonance and then slowly scanned across the resonance to longer wavelengths with increasing times. During this scanning process, the modulation instability, primary frequency lines, and cascaded four-wave-mixing processes are observed successively with the increase of the intracavity power. The overall time range is considered from nanoseconds to microseconds to confirm the stability of the generation of the Kerr frequency combs.

    References

    [1] G. P. Agrawal. Nonlinear Fiber Optics(1989).

    [2] J. Hansryd et al. Fiber-based optical parametric amplifiers and their applications. IEEE J. Sel. Top. Quantum Electron., 8, 506-520(2002).

    [3] J. Leuthold, C. Koos, W. Freude. Nonlinear silicon photonics. Nat. Photonics, 4, 535-544(2010).

    [4] M. C. Borghi et al. Nonlinear silicon photonics. J. Opt., 19, 093002(2017).

    [5] M. A. Foster et al. Broad-band optical parametric gain on a silicon photonic chip. Nature, 441, 960-963(2006).

    [6] A. C. Turner-Foster et al. Frequency conversion over two-thirds of an octave in silicon nanowaveguides.. Opt. Express, 18, 1904-1908(2010).

    [7] X. Liu et al. Mid-infrared optical parametric amplifier using silicon nanophotonic waveguides. Nat. Photonics, 4, 557-560(2010).

    [8] S. Zlatanovic et al. Mid-infrared wavelength conversion in silicon waveguides using ultracompact telecom-band-derived pump source. Nat. Photonics, 4, 561-564(2010).

    [9] X. Liu et al. Bridging the mid-infrared-to-telecom gap with silicon nanophotonic spectral translation. Nat. Photonics, 6, 667-671(2012).

    [10] P. A. Del’Haye et al. Optical frequency comb generation from a monolithic microresonator. Nature, 450, 1214-1217(2007).

    [11] Y. Okawachi et al. Octave-spanning frequency comb generation in a silicon nitride chip. Opt. Lett., 36, 3398-3400(2011).

    [12] T. Herr et al. Universal formation dynamics and noise of Kerr-frequency combs in microresonators. Nat. Photonics, 6, 480-487(2012).

    [13] T. Herr et al. Temporal solitons in optical microresonators. Nat. Photonics, 8, 145-152(2014).

    [14] A. G. Griffith et al. Silicon-chip mid-infrared frequency comb generation. Nat. Commun., 6, 6299(2014).

    [15] V. Brasch et al. Photonic chip-based optical frequency comb using soliton Cherenkov radiation. Science, 351, 357-360(2016).

    [16] X. Xue et al. Second-harmonic-assisted four-wave mixing in chip-based microresonator frequency comb generation. Light: Sci. Appl., 6, e16253(2017).

    [17] H. Guo et al. Mid-infrared frequency comb via coherent dispersive wave generation in silicon nitride nanophotonic waveguides. Nat. Photonics, 12, 330-335(2018).

    [18] N. Singh et al. Octave-spanning coherent supercontinuum generation in silicon on insulator from 1.06 μm to beyond 2.4 μm. Light: Sci. Appl., 7, 17131(2018).

    [19] B. Stern et al. Battery-operated integrated frequency comb generator. Nature, 562, 401-405(2018).

    [20] T. J. Kippenberg et al. Dissipative Kerr solitons in optical microresonators. Science, 361, eaan8083(2018).

    [21] M. A. Foster et al. Broad-band continuous-wave parametric wavelength conversion in silicon nanowaveguides. Opt. Express, 15, 12949-12958(2007).

    [22] A. B. Matsko et al. Clustered frequency comb. Opt. Lett., 41, 5102-5105(2016).

    [23] M. Haelterman, S. Trillo, S. Wabnitz. Additive-modulation-instability ring laser in the normal dispersion regime of a fiber. Opt. Lett., 17, 745-747(1992).

    [24] N. L. B. Sayson et al. Octave-spanning tunable parametric oscillation in crystalline Kerr microresonators. Nat. Photonics, 13, 701-706(2019).

    [25] S. Fujii et al. Octave-wide phase-matched four-wave mixing in dispersion-engineered crystalline microresonators. Opt. Lett., 44, 3146-3149(2019).

    [26] J. Yang et al. Coherent satellites in multispectral regenerative frequency microcombs. Commun. Phys., 3, 27(2020).

    [27] S.-W. Huang et al. Mode-locked ultrashort pulse generation from on-chip normal dispersion microresonators. Phys. Rev. Lett, 114, 053901(2015).

    [28] B. Yao et al. Gate-tunable frequency combs in graphene–nitride microresonators. Nature, 558, 410-414(2018).

    [29] A. Pasquazi et al. Micro-combs: a novel generation of optical sources. Phys. Rep., 729, 1-81(2018).

    [30] N. Picqué et al. Frequency comb spectroscopy.. Nat. Photonics, 13, 146-157(2019).

    [31] P. Feng et al. Dual-comb spectrally encoded confocal microscopy by electro-optic modulators. Opt. Lett., 44, 2919-2922(2019).

    [32] J. Zhang et al. Self-adaptive waveguide boundary for inter-mode four-wave mixing. IEEE J. Sel. Top. Quantum Electron., 26, 5100108(2020).

    [33] M. Eichenfield et al. Optomechanical crystals. Nature, 462, 78-82(2009).

    [34] F. Alpeggiani et al. Effective bichromatic potential for ultra-high Q-factor photonic crystal slab cavities. Appl. Phys. Lett., 107, 261110(2015).

    [35] A. Simbula et al. Realization of high-Q/V photonic crystal cavities defined by an effective Aubry-André-Harper bichromatic potential. Appl. Phys. Lett., 2, 056102(2017).

    [36] J. M. Dudley et al. Supercontinuum generation in photonic crystal fiber. Rev. Mod. Phys., 78, 1135(2006).

    [37] V. E. Zakharov et al. Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media. Sov. Phys. JETP, 34, 62-69(1972).

    [38] I. V. Barashenkov et al. Existence and stability chart for the ac-driven, damped nonlinear Schrödinger solitons. Phys. Rev. E, 54, 5707-5725(1996).

    [39] K. Ikeda et al. Multiple-valued stationary state and its instability of the transmitted light by a ring cavity system. Opt. Commun., 30, 257-261(1979).

    [40] L. Lugiato et al. Spatial dissipative structures in passive optical systems. Phys. Rev. Lett., 58, 2209-2211(1987).

    [41] L. Lugiato et al. From the Lugiato–Lefever equation to microresonator based soliton Kerr frequency combs. Trans. R. Soc. A, 376, 20180113(2018).

    [42] M. Haelterman et al. Dissipative modulation instability in a nonlinear dispersive ring cavity. Opt. Commun., 91, 401-407(1992).

    [43] A. B. Matsko et al. Mode-locked Kerr frequency combs. Opt. Commun., 36, 2845-2847(2011).

    [44] S. Coen et al. Universal scaling laws of Kerr frequency combs. Opt. Lett., 38, 1790-1792(2013).

    [45] S. Coen et al. Modeling of octave-spanning Kerr frequency combs using a generalized mean-field Lugiato–Lefever model. Opt. Lett., 38, 37-39(2013).

    [46] T. Hansson et al. Dynamics of microresonator frequency comb generation: models and stability. Nanophotonics, 5, 231-243(2016).

    [47] M. R. E. Lamont et al. Route to stabilized ultrabroadband microresonator-based frequency combs. Opt. Express, 38, 3478-3481(2013).

    [48] H. Guo et al. Universal dynamics and deterministic switching of dissipative Kerr solitons in optical microresonators. Nat. Phys., 13, 94-102(2017).

    [49] Y. K. Chembo et al. Spectrum and dynamics of optical frequency combs generated with monolithic whispering gallery mode resonators. Phys. Rev. Lett., 104, 103902(2010).

    [50] Y. K. Chembo et al. Modal expansion approach to optical-frequency-comb generation with monolithic whispering-gallery-mode resonators. Phys. Rev. A, 82, 033801(2010).

    [51] R. K. W. Lau et al. Effects of multiphoton absorption on parametric comb generation in silicon microresonators. Opt. Lett., 40, 2778-2781(2015).

    [52] A. D. Bristow et al. Two-photon absorption and Kerr coefficients of silicon for 850–2200 nm. Appl. Phys. Lett., 90, 191104(2007).

    [53] M. Först et al. High-speed all-optical switching in ion-implanted silicon-on-insulator microring resonators. Opt. Lett., 32, 2046-2048(2015).

    [54] S. Pearl et al. Three photon absorption in silicon for 2300–3300 nm. Appl. Phys. Lett., 93, 131102(2008).

    [55] F. Gholami et al. Third-order nonlinearity in silicon beyond 2350 nm. Appl. Phys. Lett., 99, 081102(2011).

    [56] D. K. Schroder et al. Free carrier absorption in silicon. IEEE J. Solid-State Circuits, 13, 180-187(1978).

    [57] R. R. Vardanyana et al. Modeling free carrier absorption in silicon. J. Contemp. Phys., 47, 73-79(2012).

    [58] D. Ortega et al. Analysis of “quasi-modes” in periodic segmented waveguides. J. Lightwave Technol., 17, 369-375(1999).

    [59] R. Halir et al. Waveguide sub-wavelength structures: a review of principles and applications. Laser Photonics Rev., 9, 25-49(2015).

    [60] P. Cheben et al. Subwavelength integrated photonics. Nature, 560, 565-572(2018).

    [61] X. Lu et al. Selective engineering of cavity resonance for frequency matching in optical parametric processes. Appl. Phys. Lett., 105, 151104(2014).

    [62] K. McGarvey-Lechable, P. Bianucci. Bloch-Floquet waves in optical ring resonators. Phys. Rev. B., 97, 214204(2018).

    [63] G. Moille et al. Phased-locked two-color single soliton microcombs in dispersion-engineered Si3N4 resonators. Opt. Lett., 43, 2772-2775(2018). https://doi.org/10.1364/OL.43.002772

    [64] Y. Guo et al. Ultra-flat dispersion in an integrated waveguide with five and six zero-dispersion wavelengths for mid-infrared photonics. Photonics Res., 7, 1279-1286(2019).

    [65] J. S. Penadés et al. Suspended SOI waveguide with sub-wavelength grating cladding for mid-infrared. Opt. Lett., 39, 5661-5664(2014).

    [66] J. S. Penadés et al. Suspended silicon mid-infrared waveguide devices with subwavelength grating metamaterial cladding. Opt. Express, 24, 22908-22918(2016).

    [67] H. Lin et al. Mid-infrared integrated photonics on silicon: a perspective. Nanophotonics, 7, 393-420(2018).

    [68] S. A. Miller et al. Low-loss silicon platform for broadband mid-infrared photonics. Optica, 4, 707-712(2017).

    Jianhao Zhang, Vincent Pelgrin, Carlos Alonso-Ramos, Laurent Vivien, Sailing He, Eric Cassan. Stretching the spectra of Kerr frequency combs with self-adaptive boundary silicon waveguides[J]. Advanced Photonics, 2020, 2(4): 046001
    Download Citation