Towards an analytical framework for tailoring supercontinuum generation

: A fully analytical toolbox for supercontinuum generation relying on scenarios without pulse splitting is presented. Furthermore, starting from the new insights provided by this formalism about the physical nature of direct and cascaded dispersive wave emission, a uniﬁed description of this radiation in both normal and anomalous dispersion regimes is derived. Previously unidentiﬁed physics of broadband spectra reported in earlier works is successfully explained on this basis. Finally, a foundry-compatible few-millimeters-long silicon waveguide allowing octave-spanning supercontinuum generation pumped at telecom wavelengths in the normal dispersion regime is designed, hence showcasing the potential of this new analytical approach.


Introduction
Supercontinuum (SC) generation is a complex nonlinear phenomenon that can rely on different mechanisms leading to a rich phenomenology.It can be described by generalizations of the standard nonlinear Schrödinger equation (NLSE).Particularly, high-order dispersion (HOD) plays a key role for SC generation [1].Moreover, the theoretical analysis of SC is often challenging because it must address nonlinear pulse evolution away from stationary states such as solitons.As a result, most of the studies look for approximate solutions to the NLSE based on a suitable ansatz and different perturbative methods [2,3].They give rise to effective theories that are extremely valuable to deal with complex scenarios such as soliton fission and subsequent processes [4][5][6].Although such theories are often needed, it cannot be assumed that this fully mathematical approach is always the most convenient.
In nonlinear photonics, the NLSE can be derived from Maxwell's equations when backscattered waves induced by nonlinearities are neglected (this corresponds to the slowly-varying envelope approximation in the z coordinate [7] and converts the second-order wave equation into a firstorder equation) [8] and when diffraction is addressed by means of modal methods (this allows separating transverse and longitudinal coordinates and reduces the model to a single equation if monomode propagation can be assumed) [9].These approximations do not affect the nonlinear processes allowed in the medium [7,10].If only the third-order polarization that accounts for ω 1 + ω 2 → ω 3 + ω 4 , i.e., the four-wave mixing (FWM) process, is retained, then the standard NLSE arises [11].From this point of view, although new frequency generation can be analyzed based on a mathematical solution of NLSE -a high-level concept -, it necessarily relies on FWM processes -a low-level concept.The level of detail of a description will depend on the nature of the problem itself, but also on the purposes of the description.In scenarios where many FWM processes take place or intricate interplays between dispersion and nonlinearities over the whole pulse are produced (giving rise to, e.g., solitons), low-level descriptions are often not useful.However, if the pulse spectrum develops well-defined resonances through, e.g., dispersive wave (DW) emission, one might ask if the dominant FWM processes could be identified in this case and so, a low-level description could be derived.
DW emission refers to the formation of spectral resonances due to HOD [12].This mechanism plays a crucial role for SC generation [1,2,14,15] and particularly, it was recently exploited for efficient wavelength conversion into the visible range in photonic crystal fiber (PCF) [16,17].It was pointed out early that optical solitons perturbed by HOD [e.g., propagating in the neighborhood of the zero-dispersion frequency (ZDF)] radiated at specific frequencies [12].An effective radiative process involving the soliton [i.e., a nonlinear wave-packet] and a smallamplitude signal [i.e., a linear, often narrow-band, wave-packet] was considered to analyze this phenomenon [13].With this approach one succeeded in finding the phase-matching condition that leads to the resonant frequency.
Spectral resonances observed in SC in PCFs and integrated waveguides were usually interpreted as DWs radiated by the solitons emerged after soliton fission [1].According to [18], a higher-order soliton is a bound state of fundamental solitons.However, it becomes unstable in the presence of small perturbations [e.g., HOD or stimulated Raman scattering (SRS)] and tends to split into its constituent solitons after propagation.Consequently, the pulse splitting that precedes radiative processes can be analyzed as soliton fission, and the subsequent generation of spectral resonances (also called nonsolitonic radiation) as DWs radiated by the solitons emerged after fission [14].After these stages, intricate scattering processes involving the ejected solitons and DWs can take place [2].These processes are also of key importance to explain the SC features [4][5][6].
Although SC generation was initially explained in these terms [1,2,14,15], Austin et al. pointed out that DWs are radiated before the soliton fission [19].This is also in accordance with recent results in integrated waveguides [20,21].Interestingly, DW emission by pulses that propagate in the normal dispersion regime (called generalized DWs) has also been demonstrated [22,23].Despite the importance of new frequency generation processes triggered by soliton fission, these observations on DWs reveal that SC assisted by DW generation does not need pulse splitting, even when pumped in the anomalous dispersion regime.In these cases, the soliton dynamics is not essential and low-level strategies could be useful, in line with early interpretations of some SC spectra based on intrapulse FWM processes [24][25][26].Here our aim is to obtain a unified analytical description of such DW emission, in both normal and anomalous dispersion regimes, through a low-level approach.Since numerical solutions of the NLSE are not required, our formalism provides, in addition, straightforward guidelines for designing SC sources.
An illustrative example of low-level descriptions has been the interpretation of DW emission by nearly nondispersive pulses in terms of cascaded FWM processes by Erkintalo et al. [27].It did not only provide a deeper understanding of the radiation emitted by solitons, but also demonstrated that spectral resonances can be predicted without a complete solution for the NLSE.This work also reported the appearance of spectral resonances from pulses pumped in the normal dispersion regime close to the ZDF.Another very recent example has been the low-level analysis of the effective reflection of a weak linear wave induced by a soliton based on certain FWM processes [28,29].
If frequency dispersion in time, i.e., the dispersion of the instantaneous pulse frequency along the pulse duration, is significant, then new mechanisms can also arise [22].In the normal dispersion regime, the optical wave-breaking (OWB) process can take place [30] provided the system is sufficiently nonlinear [31].Conforti et al. studied this regime based on the shock-wave solution [3].Although several cases can be understood from this high-level concept and some particular cases can be solved analytically [32], the group-velocity of the shock-wave must be numerically determined in general [3], which limits its applicability in design tasks.Recently, we proposed a combined approach where self-phase modulation (SPM) leading to OWB and DW emission mechanisms are studied based on the envelope evolution -a high-level concept -and discrete wave mixing -a low-level concept -, respectively [23].One FWM process enabled by frequency dispersion in time and favored by a suitable HOD was identified for the DW emission, leading to an analytical expression that relates the position of the spectral resonance with the dispersion and the nonlinear coefficient of the waveguide and input pulse features.In this paper we extend our approach to new scenarios in both normal and anomalous dispersion regimes, where DWs are produced according to several FWM processes induced by the dispersion of frequencies in time.
The approach that we present in this paper will be applicable to both normal and anomalous dispersion regimes.This is in contrast with high-level formalisms which are fundamentally different for the normal and anomalous dispersion regimes because they rely on nonlinear envelope states.Particularly, these nonlinear states correspond to solitons in the anomalous dispersion case [1,2,14] and a shock-wave in the normal dispersion case [3,32].However, there are scenarios where a simple picture of the pulse as a coherent superposition of monochromatic waves, that can disperse in time, can be suitable.In such cases, a unified approach for both the normal and anomalous dispersion regimes should be feasible [23,27] and we here present such an approach for SC produced in scenarios that do not involve pulse splitting (thus before soliton fission occurs in the anomalous dispersion case [33,34]).
In Section 2, we present a fully analytical approach to define key stages in the SC generation in both normal and anomalous dispersion regimes.We construct a picture of the frequency distribution in time based on these stages.In Section 3, several FWM processes are identified to produce direct and cascaded DW emission.This allows the derivation of analytical expressions that characterize the different scenarios and provide the resonance frequencies.Finally, Section 4 will be devoted to the application of this framework in actual (nonideal) waveguides.On the one side, an octave-spanning SC generation in a 2 mm-long foundry-compatible (220-nm-thick) silicon-on-insulator (SOI) strip waveguide [35] pumped at telecom wavelengths in the normal dispersion regime will be presented.On the other side, a recently reported SC spectrum obtained in a silicon nitride waveguide pumped in the anomalous dispersion regime [21] will be analyzed according to our framework.

A propagation equation for generalized lengths
A suitable definition of the different stages of the pulse evolution is extremely useful to depict a simplified, although consistent, model of the SC dynamics.The most extended scenario valid in the anomalous dispersion regime is composed of an initial pulse compression stage driven by SPM and anomalous group-velocity dispersion (GVD), followed by the soliton fission process due to higher-order effects (HOD, SRS or even noise) and subsequent radiation by the emerging solitons, and scattering processes involving such radiation and solitons [1,2].Usually these stages are defined in terms of the changes experienced by the pulse (in both the time and frequency domains), thus requiring the numerical solution of the propagation equation.Nevertheless, such complete solutions contain information about the pulse details that is not essential to define the stages.This approach is not particularly efficient for determining the different SC stages, because the goal of defining these stages is precisely to provide a simpler analysis of the pulse evolution.
Such a simpler analysis can be enabled by looking for some properties of the pulse and its spectrum that allow the characterization of the propagation stages without solving numerically the NLSE.An intuitive example is given by the classical dispersive length, L D = T 2 0 /| β 2 |, where T 0 represents the pulse duration and β 2 is the GVD; and the nonlinear length, L NL = 1/(γ 0 P 0 ), where γ 0 is the nonlinear coefficient and P 0 corresponds to the input pulse peak power.These parameters allow comparing the strength of the processes at work, namely, GVD and SPM (provided the pulse shape is not notably altered).Similarly, we recently proposed the generalized dispersive and nonlinear length functions, L D (z) and L NL (z), respectively, to extend this kind of analysis over the entire propagation distance [36].For the convenience of the reader, we briefly review their meaning here.For the sake of clarity, let us consider the standard NLSE (although our remarks remain valid when HOD is included), where ω 0 is the carrier frequency, f (ω is the Fourier transform centered at ω 0 , and Ã is the envelope of the analytic signal of the electric field [11].Next we take averages of the right-hand side terms in Eq. ( 1) and introduce the following inverse length functions, See details of the formula for ξ GWB in Table 1.
On the one hand, L D and L NL represent the length scales over which GVD and SPM, respectively, act at any z distance.On the other hand, the Hamiltonian conservation law [37] requires that [36].This property allows the nonlinear pulse propagation to be interpreted as a competition between the generalized inverse lengths.Therefore, these functions can provide useful insight to identify the dynamic stages in the SC generation.Since Eq. ( 1) has actually one degree of freedom [37,38], parametrized by N 2 = L D /L NL = T 2 0 γ 0 P 0 / β 2 , we from now on work with a normalized propagation distance, ξ = z/L D .
Figure 1(a) represents the typical evolution of . The crosspoint at ξ GWB (the -normalizedgeneralized wave-breaking distance) is a notable feature of this regime because it can define two stages without a full knowledge of the pulse envelope: Since L −1 NL (ξ) > L −1 D (ξ) for ξ < ξ GWB , the spectral broadening (pointed out by the growth of L −1 D , which is proportional to the variance of the spectral distribution) should be mainly driven by SPM in this stage.Moreover, L −1 NL captures the pulse broadening induced by β 2 > 0 through its decrease.For ξ > ξ GWB , dispersive effects should dominate the pulse propagation.Figure 1(b) corresponds to the anomalous dispersion regime [i.e., sign( Although no crosspoint appears in this case, again two different stages can also be identified in a natural way.The generalized lengths initially experience a relatively slowly varying monotonic evolution along ξ and evolve into an oscillatory behavior afterwards.Analogously to the previous case, the spectral broadening should be mainly induced by SPM initially.We point out that the pulse compression due to β 2 < 0 explains the growth of L −1 NL .From this approach, no fundamental difference between the normal and anomalous dispersion regimes is observed.Similarly, the self-compression and self-decompression distances were also proposed without any particular distinction between the anomalous and normal regimes [34].Therefore, we define the (normalized) generalized wavebreaking distance, ξ GWB , for both sign( β 2 ) = ±1, as L −1 D (ξ GWB ) = sign( β 2 )L −1 NL (0)/2.We will return later to the motivation of this definition (see next page).We now present a propagation equation for L −1 NL and, correspondingly, for L −1 D .We take only two effects into account to determine the evolution of the generalized lengths for ξ < ξ GWB .Firstly, we assume that SPM rules the spectral broadening in this stage.This assumption is expected to be completely valid when L −1 NL (0) L −1 D (0).Secondly, we neglect any pulse reshaping and consider an effective broadening or compression in time (see Appendix).Based on Eq. ( 1) and keeping in mind these assumptions, we derive in Appendix where σ 2 is an input-pulse form factor (see its value in the caption of Table 1 and details in Appendix), and s 2 = sign( β 2 ).To check the validity of Eq. ( 4), first we numerically solve Eq. ( 1) to obtain ξ GWB according to its definition.Then we analytically solve Eq. ( 4) and impose , so that we obtain In Fig. 1(c) and Fig. 1(d) we compare the analytical results based on Eq. ( 5) with those evaluated numerically.Despite the assumptions made in Eq. ( 4), good agreement is observed.This indicates that our approximations (see Appendix) are in accordance with the conditions of the scenarios illustrated in Fig. 1.For ξ > ξ GWB , the dynamics become notably different in Figs.1(a) and 1(b).In both regimes, L −1 D strengthens as the pulse propagates.It also indicates that, besides the spectral broadening, frequency dispersion starts to play a more important role.It is well-known that in the normal dispersion regime, the dispersion can produce frequency overtaking in the pulse tails that leads to strong oscillations in the pulse intensity.This phenomenon was interpreted early on as shockwave formation and OWB [30,31,39,40].Furthermore, the observation of these oscillations has been greatly improved very recently by means of an optical sampling oscilloscope [41].Interestingly, although shock-wave formation enables nonlinear mixing involving frequencies in the pulse tails [3,23,31,32], these mixing processes could also appear due to a partial overlapping.From this point of view, shock-wave and OWB would be extreme cases of frequency overlapping in time.The latter suggests that the GVD-induced spreading in time of frequencies generated initially by SPM can enable new nonlinear mixing between these frequencies in both normal and anomalous dispersion regimes, the only difference between the two cases being the frequencies involved: in the normal dispersion regime, it are those contained in the pulse tails, (ω − SPM , ω 0 ) or (ω 0 , ω + SPM ); and, in the anomalous dispersion regime, it are those contained in the central part of the pulse, (ω − SPM , ω + SPM ), where ω ± SPM = ω 0 ± δω SPM with δω SPM representing the maximum SPM-induced chirp (see Fig. 2).This similar treatment for both dispersion regimes motivates the definition of one generalized wave-breaking distance.When using this definition, one has to Fig. 2. Illustration of the frequency overlapping and the resulting FWM processes in (a) the normal dispersion regime, s 2 = 1, and (b) the anomalous dispersion regime, s 2 = −1.The schematic plots include the absolute instantaneous frequency, ω(t) = ω 0 +δω(t) (continuous lines), and instantaneous power, P(t) (dashed lines).Thick lines highlight the frequencies that can overlap.The time shifting induced by dispersion can induce the FWM processes that are represented.take into account that in the normal case the spectral broadening due to SPM is mainly produced before ξ GWB , whereas in the anomalous case, due to the pulse compression, L −1 NL , and therefore SPM, experiences the most important growth just after ξ GWB .Clearly, this feature cannot be neglected in a reasonable model.
From Eq. ( 4), we can estimate the frequencies that could interact in the second stage (i.e., ξ > ξ GWB ), provided the distance at which frequency overlapping takes place, ξ ol , is known.In the anomalous dispersion scenario [see Fig. 1(b)] the distance at which L −1 NL attains its first maximum, L −1 NL (ξ ol ) = max(L −1 NL ), could be an appropriate choice because it is related to the maximum compression distance in the time domain and hence maximum frequency overlapping.In contrast, it is not straightforward to define such a distance in the normal dispersion case.Based on the above qualitative description of ξ ol , the pulse broadening at ξ ol can be estimated to be roughly twice that attained at ξ GWB .Consequently, in the normal dispersion scenario, we define L −1 NL (ξ ol ) = L −1 NL (ξ GWB )/2.Note that the aim of defining ξ ol is to facilitate evaluating the frequency ranges that will be involved in the nonlinear mixing in the second stage.
Let us first analyze the anomalous dispersion case.Equation (4) takes into account SPM and pulse compression induced by GVD.As a result, it predicts that L −1 NL diverges for ξ > ξ GWB .Obviously, it is an artificial divergence as we can observe in Fig. 1(b).In fact, for ξ > ξ ol , the pulse experiences both spectral and temporal narrowing and broadening, hence avoiding any divergence or pulse collapse.The physical reasons explaining that behavior are analogous to those accounting for the maximum pulse compression of a positively chirped pulse in a linear waveguide with anomalous dispersion [38].Indeed, a clear similarity with the case we are dealing with arises, the SPM-induced chirp playing the role of the initial positive chirp in the linear case.Furthermore, this analogy suggests the following strategy to evaluate the spectral broadening at ξ ol : if we define a chirped Gaussian pulse at ξ GWB satisfying L −1 D (ξ GWB ) = sign( β 2 )L −1 NL (0)/2, then we can estimate L −1 NL (ξ ol ) as defined above from the analytical results of the linear case (see Appendix).Once L −1 NL (ξ ol ) is known, both ξ ol and the range of frequencies eventually overlapping can be estimated through an extrapolation of Eq. ( 4).
A similar approximation can be applied in the normal dispersion regime, but then considering that L −1 NL (ξ ol ) = L −1 NL (0)/4.A synopsis of the most important analytical results presented here Normal dispersion regime, β 2 > 0 Anomalous dispersion regime, and derived in Appendix is shown in Table 1.We use these results in the next section, where several scenarios that exploit DW emission to generate SC are identified.

Direct and cascaded dispersive wave emission
The analysis of the generalized lengths allows identifying the frequencies that can overlap, thus can generate new frequencies through FWM in the second stage, after ξ GWB (see Table 1).We now study new frequency production through direct processes, i.e., one FWM interaction fed by SPM (both pump and signal waves generated in the first stage by SPM), but also through cascaded processes, i.e., FWM interactions fed by other FWM processes also produced in the second stage.As will be shown soon, direct DW emission can take place in waveguides with, in addition to GVD, third-order dispersion, while cascaded processes require third and fourth-order dispersion.Equation (1) must then be extended, in which case it is called generalized NLSE (GNLSE), where The evolution of the generalized lengths is assumed to be ruled by SPM and GVD and the most important effect of HOD is to favor specific FWM processes in the second stage.In this section we pursue to determine the HOD, and therefore the waveguide dispersion curve, that induces DW emission at targeted frequencies to further enhance the pulse spectral broadening in the second propagation stage.Consequently, at this point, the FWM processes to be exploited must be selected a priori to derive the HOD.First, we address direct DW emission in the normal dispersion regime.Among the two ranges of frequencies that get closer in the time domain due to the GVD [e.g., the frequency range present in the trailing pulse edge (ω 0 , ω + SPM )], the FWM process 2 ω + SPM → ω 0 +ω DW with ω + SPM as pump and ω 0 as signal [see Fig. 2(a)] is more favorable in terms of power [see Fig. 2(a), P(t 1 ) P(t 2 )] than the reciprocal process where ω 0 acts as pump wave and ω + SPM as signal wave.Consequently, new frequency generation is expected to be more pronounced through the former channel.[The same considerations can be applied to the range (ω − SPM , ω 0 ) and the process 2 ω − SPM → ω 0 +ω DW .]In addition, it can provide the highest (lowest) new frequency.If phase matching is imposed in this process, then β(ω DW ) + β(ω 0 ) − 2 β(ω + SPM ) = β 2 (ω + SPM )[δω SPM ] 2 = 0 (the nonlinear phase mismatch is assumed to be negligible) when dispersion up to the third-order ( β 3 ) is considered [38].As a result, the phase matching condition can be rewritten as β 3 = − β 2 /δω SPM , as we already reported in [23].Note that the FWM pump frequency, ω + SPM (or ω − SPM ), corresponds, in this case, to the ZDF, in line with previous experimental works where spectral peaks were observed when using input pump pulses at normal dispersion regime [26,42].Interestingly, this condition also leads to a group-velocity matching between the signal and idler waves involved in this FWM process, β 1 (ω 0 ) = β 1 (ω DW ).Therefore, in this case, both phase matching and group-velocity matching can be attained simultaneously by means of just one additional degree of freedom, namely, that related to β 3 .If group-velocity dispersion allows space-time overlapping as described above and phase matching permits gain in the nonlinear FWM process, group-velocity matching becomes important to sustain the intrapulse frequency conversion [24,42].In Fig. 3(a), we show the numerical solutions of Eq. (6) for the system described in Fig. 1(a) with β 3 given by the corresponding expression collected in Table 2 and δω SPM = δω ol (see Table 1).An excellent agreement is observed between the idler wave frequency produced through the above mentioned FWM process, ν th DW , and the frequency of the spectral resonance appearing in the numerical simulation.
Second, we study direct DW emission in the anomalous dispersion regime, where GVD induces frequency overlapping between the frequencies (ω − SPM , ω + SPM ) because β 2 < 0. In this case, we focus on 2 ω , and again it leads to β 3 = − β 2 /δω SPM and to group-velocity matching between waves that take part in this process, β 1 (ω − SPM ) = β 1 (ω DW ).Analogously to the direct DW emission in the normal dispersion regime, the FWM pump frequency corresponds to the ZDF, similarly to [24].
When comparing the idler frequency analytically determined, indicated by ν th DW , with the spectral resonance obtained numerically in Fig. 3(b), we again observe a good agreement.
Next we study some cascaded scenarios keeping in mind the idea of maximizing the spectral broadening.It is known that a cascade can be induced through nonresonant intermediate processes [27].In the normal dispersion regime, we aim to get 2 ω + SPM → ω 0 + ω i (ω i is not necessarily resonant), and subsequently 2 ω i → ω 0 + ω CDW .Therefore, if only phase matching and groupvelocity matching on the last process of our cascade, 2 [38] and β 1 (ω 0 ) = β 1 (ω i ).Note that, in this case, the group-velocity matching involves the pump and signal waves of the second process of the cascade.These conditions lead to β 3 = −5 β 2 /3δω SPM and β 4 = β 2 /δω 2 SPM .We want to exploit this process with as few requirements as possible.As such, a GNLSE with one additional degree of freedom, that linked to β 4 , is used.The numerical output spectrum plotted in Fig. 3(c) shows a resonance in line with the theoretical position of the cascaded DW, ν th CDW .The cascade studied in the anomalous dispersion regime includes the following processes: 2 ω + SPM → ω − SPM + ω i and 2 ω − SPM → ω i + ω CDW .Similarly to the previous case, we impose In this case, we derive β 3 = −2 β 2 /9δω SPM and β 4 = −2 β 2 /3δω 2 SPM .The numerically simulated resonance observed in Fig. 3(d) again corresponds well with the theoretical ν th CDW .This agreement in several and notably different scenarios strongly supports our approach as a valid tool to address SC generation relying on SPM and DW emission.
Note that the HOD used here served to stimulate a resonance at ω CDW and a spectral peak has appeared there according to our theory.Nevertheless, other resonant processes can also occur; in fact, an additional resonance is excited in this scenario, contributing to the spectral broadening.From this point of view, the above results indicate that the expressions in Table 2 (including the information provided by Table 1) represent sufficient conditions to efficiently induce DWs.Note that, regarding design tasks, this is the most useful information.Furthermore, the analysis of phase and group-velocity matching, i.e., the conditions leading to the results in Table 2, can also improve the understanding of any DW emission.To illustrate this, the following section will be devoted to the analysis of all the spectral resonances appearing in the numerically calculated spectra in actual waveguides.That strategy is also valid to interpret the high-frequency resonance in Fig. 3(d).
It is worth remarking that, in the anomalous dispersion regime, DW formation is often interpreted based on the coupling between a soliton and linear waves [2,14,15,19].This highlevel theory assumes an ansatz to Eq. ( 6) that consists of two parts: a solution to Eq. ( 1) that belongs to the family of fundamental solitons and a small-amplitude linear wave [2,4,13,38].
After neglecting the nonlinear coupling terms in the resulting equation as indicated in [4,13], the following phase matching condition is derived, β(ω DW ) = β 0 + β 1 (ω DW − ω 0 ) (in case one neglects the nonlinear contribution to the soliton wavenumber).Interestingly, our framework can explain this ansatz because the distance at which L −1 NL attains its first maximum is related to the distance where the pulse chirp becomes zero.As a result, a generalized fundamental soliton can be identified around this distance (the self-compression point) due to its straight-line dispersion curve.Moreover, our approach also identifies FWM (low-level) processes capable of emitting DWs and conditions allowing efficient radiation.(Note that we do not impose group-velocity matching between the soliton and the DW, but it affects waves involved in the FWM process that is selected.) We would like to stress that the new physical insights obtained from our approach also pave the way to an inverse nonlinear engineering.The procedure hereto can be outlined in three steps.In the first one, the inverse nonlinear design step, the results above can be used to obtain the optimal dispersive features of a waveguide exhibiting a given nonlinear response: From the available input pulse (pulse shape, N and center wavelength), Table 1 and Table 2 determine the dispersion curves and the characteristic distances associated to several possible scenarios.The desired output spectral bandwidth, the features of the available pump source, and the degrees of freedom for the waveguide design will determine the most convenient scenario to be selected.In the second step, the inverse linear design step, a waveguide cross-section is obtained from the target dispersion curve defined in the previous step [45].Finally, once a suitable waveguide is found, and β(ω) and γ(ω) have been calculated, the GNLSE, including higher-order effects, is numerically solved in order to check the output spectrum of the realistic waveguide design [23].
Most of the times, the target dispersion could not be exactly realized with realistic waveguide cross-sections.For example, the waveguide thickness is often constrained in photonic fabrication foundries.Then dispersion control is restricted to the optimization of a few number of parameters (e.g., the waveguide width [43]) and consequently, dispersion engineering is limited.Even in these cases, Table 1 and Table 2 provide valuable information since it sets a target that should be realized as closely as possible with the dispersion engineering possibilities at hand.What is more, our framework can also be helpful to understand the origin of spectral resonances produced in a waveguide that does not feature exactly the HOD presented in Table 2. To illustrate the latter, we apply our framework to concrete SC examples based on waveguides that have been already fabricated.

Supercontinuum relying on generalized wave-breaking: case studies
Designing a new octave-spanning SC source in a 220 nm-thick silicon waveguide pumped at 1550 nm As a first case study, we use our framework for designing a SC source based on a SOI strip waveguide with a fixed thickness of 220 nm, in line with the fabrication rules of multi-project wafer runs of photonic foundries [46].From the dispersion profile in [47] for such a foundrycompatible (i.e., 220 nm-thick) SOI waveguide with a width of 900 nm, we can derive that the type of scenario corresponding to Fig. 3(c) could be feasible in that waveguide [see inset in Fig. 4(a)].Guided by Table 1, we consider 50 fs-long sech pump pulses with 150 W of peak power at 1550 nm. Figure 4(a) shows the octave-spanning output spectrum after propagation through this waveguide along 2 mm.The solid blue curve is the simulated spectrum taking into account linear losses of 2 dB cm −1 [43], two-photon absorption (TPA), the dispersion of the nonlinear coefficient [23,48,49], free-carrier-related effects and SRS [43] and the spectrum without these higher-order effects but with an effective nonlinear coefficient [23] is indicated by the green dashed curve.Common features between Fig. 4(a) and Fig. 3(c) become apparent.Consequently, we analyze the output spectrum in Fig. 4(a) keeping in mind the dynamics corresponding to Fig. 3(c).Table 2. Four different scenarios for SC relying on direct and cascaded DW.The optimal HOD parameters that induces group-velocity matching and phase matching for those FWM and the spectral broadening produced through each mechanism (ω DW or ω CDW ) are also included.δω SPM can be estimated by δω ol in Table 1.

sign( β 2 )
Favored FWM processes The inverse generalized lengths, L −1 NL and L −1 D , in Fig. 4(b) evolve similarly as in Fig. 1(a).They also provides valuable information about propagation distances of interest.In the first stage (z < z GWB ), SPM broadens the spectrum.The analysis of the relative group velocity, β 1 (ω) − β 1 (ω 0 ), in Fig. 4(c) allows identifying which frequencies generated through SPM in the first stage could overlap due to their GVD in the second stage (z > z GWB ).In this case, red-shifted frequencies, ω < ω 0 with β 1 (ω) < β 1 (ω 0 ) (indicated in the green background), tend to approach ω 0 in the leading pulse edge [see also Fig. 2(a)].Among the FWM processes favored by the frequency-dispersion-induced overlapping, those allowed by energy conservation and phase matching [see Fig. 4(d)] can extend the spectral broadening initiated by SPM.Note that the spectrum spans an octave through this mechanism also in the realistic simulation that includes all effects in the silicon waveguide [see solid blue line in Fig. 4(a)].
Moreover, based on Section 3, spectral resonances arise if both phase matching and groupvelocity matching are fulfilled.Phase matching can be evaluated through the linear phase mismatch, ∆ β = β(ω i ) + β(ω s ) − 2 β(ω p ), of the degenerate FWM processes that can take place, 2 ω p → ω s + ω i .Unlike the group-velocity matching, the phase-matching condition depends on linear and nonlinear contributions.Although the linear phase mismatch can be readily evaluated, the nonlinear term would require details of the pulse due to its dynamical nature.For these reasons, we choose to restrict our study to ∆ β, while keeping in mind an uncertainty related to the nonlinear mismatch.This analysis is done in Fig. 4(d).In addition, the solid green line that intersects with "S" in this figure relates the pump and signal waves involved in processes with group-velocity matching between pump and signal waves.Accordingly, the spectral resonance ω CDW can be interpreted as the idler of a FWM process involving the pump and signal waves depicted as "P" and "S", respectively.This is in agreement with our analysis in Fig. 3(c).
It is worth reminding that group-velocity and phase matching are conditions allowing efficient FWM.Notwithstanding, new frequency generation can also take place if these requirements are not fulfilled provided pump and signal waves can interact.These processes, besides broadening the spectrum, can also feed other FWM interactions, the second step in cascaded processes [e.g., 2ω i → ω 0 + ω CDW in Fig. 3(c) and Table 2].Indeed, these processes can contribute to the pump wave "P" production according to Fig. 4(d).Of course, higher-order effects, e.g., nonlinear losses, can affect these processes.The comparison between the green dashed and blue solid lines in Fig. 4(a) indicates the position of ω CDW is not altered, but its efficiency is reduced (as expected due to TPA) and new resonant processes become more efficient (see the spectral peak around 140 THz).
This example shows the (nearly unexplored) potential of pumping in the normal dispersion regime to produce broadband light in silicon waveguides compared to previous approaches pumping in the anomalous dispersion regime [20,43,50].It can represent an important breakthrough towards the experimental demonstration of an octave-spanning SC spectrum in silicon pumping in the short-wavelength infrared since, according to our knowledge, it has only been achieved in SOI waveguides pumping beyond 2.2 µm, i.e., above the TPA threshold of silicon [44].

Analyzing SC results in literature
As a second example of how to use our framework in practice, we analyze the dynamics of the experimental SC reported recently in [21].There a silicon nitride waveguide was used, pumped in the anomalous dispersion regime with 92 fs-long Gaussian pulses with 260 W of peak power at 1030 nm.In Fig. 5(a), we show the spectra at different distances, including the output spectrum corresponding to Fig. 2 of [21].We numerically calculate the evolution of the generalized lengths for this system and plot the result in Fig. 5  scenario represented in Fig. 3(d) with cascaded DW emission.In Fig. 5(c), we highlight the range of frequencies that could overlap due to GVD by means of a green background.SPM generates blue-shifted frequencies (ω > ω 0 ) at positive times and red-shifted frequencies (ω < ω 0 ) at negative times around the pulse center.As a result, blue-shifted frequencies that move faster, β 1 (ω) < β 1 (ω 0 ), and red-shifted frequencies that move more slowly, β 1 (ω) > β 1 (ω 0 ), will tend to get closer to each other [see Fig. 2(b)].Among the potential interactions, those with suitable phase matching and group-velocity matching conditions will generate new frequencies more efficiently.Analogously to the first case studied in this section, FWM processes that satisfy both conditions [see Figs.5(c) and 5(d)] yield spectral resonances observed in Fig. 5(a) for the green solid curve corresponding to a propagation distance of 7.5 mm ( i.e., the distance where L −1 NL attains its first maximum).This physical explanation for the observed spectral resonances has not been presented before and provides new insights in the SC dynamics.
If the pulse propagates longer distances [see spectrum depicted by the blue dotted curve in Fig. 5(a)], additional processes can also take place, including the blue-shift in the resonance around 400 THz [19] and soliton fission.Since conservation of the initial pulse train structure is important for several practical applications, we have focused in this paper on spectral broadening mechanisms before pulse splitting [33].Finally, we point out that, in addition to the cases studied here, we have verified that our framework also allows analyzing a wide panoply of other results reported in the literature [15,19,22,23,26,42,43].

Fig. 1 .
Fig. 1. (a)-(b): Evolution of the generalized lengths for N = 10 in (a) the normal dispersion and (b) the anomalous dispersion regimes.(c)-(d): Comparison of the analytically calculated (normalized) generalized wave-breaking distance (red squares) with the numerical results (blue circles) in (c) the normal dispersion and (d) the anomalous dispersion regimes.The parameters considered in theses cases are included in the insets of (c) and (d).See details of the formula for ξ GWB in Table1.

Fig. 3 .
Fig. 3. Direct DW emission in (a) the normal regime and (b) the anomalous regime.Cascaded DW emission in (c) the normal regime and (d) the anomalous regime.HOD has been calculated according to the analytical results in Table2using δω SPM = δω ol and| β 2 | = 1 ps 2 m −1 , γ 0 = 4 W −1 m −1 , T 0 = 0.5 ps and ν 0 corresponding to 1550 nm.The position of the resonances are successfully predicted in our framework in all these cases.The propagations have been stopped at the distances ξ, where the resonances achieve their maximum power levels while additional processes not included in our model have not impacted the dynamics yet.Yellow dashed arrows represent the spectral broadening relying on SPM and green solid arrows indicate the FWM processes that have been considered.The dispersion profiles, including arrows to indicate the pumping frequency, have been added as insets.

Fig. 4 .
Fig. 4. (a) Output spectra simulated through Eq. (6) including higher order effects (blue solid curve) and without them (green dashed line).(b) Evolution of the generalized lengths [cf.Fig. 1(a)].(c) Plot of the relative β 1 , i.e., the inverse of the group velocity.The green window includes the frequencies that can overlap in the leading pulse edge.A and N indicate anomalous and normal dispersion, respectively.(d) Linear phase mismatch (only negative values are represented).The green solid curve points out processes with group-velocity matching (see details in the text).
(b).The behavior of the generalized lengths resembles that in Fig.1(b) and the total waveguide length (z WG = 8 mm) is slightly beyond the distance where L −1 NL attains its first maximum.Whereas the generalized lengths in Fig.5(b) and Fig.1(b) are similar, the features of the spectrum in Fig.5(a) strongly suggest that it is related to the

Fig. 5 .
Fig. 5. (a) Output spectra at several distances between z GWB and the total length of the waveguide, z WG , according to [21].(b) Evolution of the generalized lengths [cf.Fig. 1(b)].(c) Plot of the relative β 1 , i.e., the inverse of the group velocity.The green window includes the frequencies that can overlap.A and N indicate anomalous and normal dispersion, respectively.(d) Linear phase mismatch.The green solid curve points out the processes with group-velocity matching (see details in the text).

Table 1 .
Characteristic lengths, ξ GWB and ξ ol , and the maximum chirps generated by SPM at such lengths.Here a can be considered an auxiliary variable, σ 2 , and Υ are input-pulse form factors.In particular, they equal to 16/35 and 4/ Gaussian pulse.σ2 is an auxiliary parameter equals to 1/ √ 2 (see details in Appendix).Remember that N 2 = L D /L NL .