Review Article  Open Access
Manipulating the Propagation of Solitons with SolidCore Photonic Bandgap Fibers
Abstract
We review the dynamics of soliton selffrequency shift induced by Raman gain in special solidcore photonic bandgap fibers and its consequences in terms of supercontinuum generation. These photonic bandgap fibers have been designed to allow nonlinear experiments in the first bandgap without suffering from significant loss even when working close to the photonic bandgap edge. We studied experimentally, numerically, and analytically the extreme deceleration of the soliton selffrequency shift at the longwavelength edge of the first transmission window. This phenomenon is interpreted as being due to a large variation of the groupvelocity dispersion in this spectral range and has been obtained with no significant power loss. Then, we investigated experimentally and numerically the generation of supercontinuum in this kind of fibers, in both spectral and temporal domains. In particular, we demonstrated an efficient tailoring of the supercontinuum spectral extension as well as a strong noise reduction at its longwavelength edge.
1. Introduction
Solidcore photonic bandgap (PBG) fibers are one class of microstructured optical fibers in which light is confined in a lowindex solid core by the PBGs of the cladding [1]. In the case of twodimensional PBGs, the periodic microstructured cladding is usually composed of periodic highindex inclusions embedded in a lowindex background [1]. The core region then corresponds to a defect (lack of inclusion) in the center of the periodic structure. Because of the intrinsic nature of the PBG waveguidance, the fiber transmission properties are characterized by discrete spectral bands [1, 2] in which the groupvelocity dispersion (GVD), attenuation, and effective mode area of the fundamental core mode are strongly wavelength dependent [3], especially near the PBG edges [1, 4]. These singular characteristics are thus of particular interest in the field of nonlinear (NL) fiber optics. The potential of solidcore PBG fibers for NL propagation experiments has already been pointed out by the report of soliton propagation and phasematched dispersive wave generation [5, 6], even across different PBG [7]. Degenerate fourwave mixing [8], frequency doubling and tripling [9], and supercontinuum (SC) generation across adjacent PBGs have also been investigated numerically [10].
In this paper, we present an overview of our recent work showing how the particular properties of solidcore PBG fibers can influence soliton propagation and SC generation. First, we detail the linear properties of the fabricated solidcore PBG fibers under investigation through experiments and numerical modelling. Then, we demonstrate experimentally that the soliton selffrequency shift (SSFS) effect can be strongly mitigated [11]. A simple analytical model indicates that this is caused by the strong thirdorder dispersion experienced by solitons as they approach the PBG edge. Finally, we show that this SSFS quasicancellation allows to control the spectral extent of infrared SC in the long pulse regime [12] and, simultaneously, leads to a significant reduction of spectral power fluctuations at the SC longwavelength edge [13].
2. Linear Properties of the SolidCore Photonic Bandgap Fibers
The solidcore PBG fibers used in various experiments reported in this paper are all based on the same double periodicity structure [11, 12], although their optogeometrical characteristics (summarized in Table 1) slightly vary from one fiber to another. The scanning electron microscope (SEM) image shown in Figure 1 illustrates such a double periodicity solidcore PBG fiber. The periodical highindex resonators (light gray regions in Figure 1) present in the cladding are responsible for the PBG guidance and are made of germaniumdoped silica with a parabolic refractive index profile and a maximal index difference of (as compared to pure silica). They are embedded in a pure silica background (dark gray). Because we decided to work in the first PBG (mainly to minimize the effective area in order to maximize the NL coefficient), air holes of diameter (depicted in black in Figure 1) were periodically added into the cladding structure for reducing the otherwise relatively high confinement losses in this PBG [14, 15]. The core region corresponds to a defect (lack of air hole) in the center of the periodic structure. The geometrical parameters of the three solidcore PBG fibers investigated here are summarized in Table 1, where the pitch corresponds to the distance between two neighboring germaniumdoped inclusions.

Although all fibers were characterized the same way, we only present here the linear properties of fiber 1, whose SEM image is depicted in Figure 1. Other fibers present similar characteristics, with a slight redshift of the PBG transmission windows due to a slightly greater pitch . Figure 2 shows the experimental transmission spectrum of a white light SC source through a 2 m long sample of the solidcore PBG fiber. The first PBG is very broad and extends from 740 to 1535 nm, whereas the second one is much sharper and extends from 510 to 650 nm. As for usual solidcore PBG fibers, typical resonances can be observed in the first ring of highindex resonators [4], with the number of intensity maxima indicating the PBG order (insets of Figure 2). The red curve (right axis) shows confinement losses across the first PBG calculated with a finite element mode solver. They are less than 1 dB/km across most of the PBG, and as expected they drastically increase to the dB/m range near the PBG edge. Figure 3 shows the calculated GVD curve of the solidcore PBG fiber 1 (solid black line), as well as measurements (markers) performed with a low coherence interferometric method [16]. The GVD slightly increases from the zerodispersion wavelength (ZDW), located around 900 nm, to about 200 ps/nm/km at 1510 nm. For higher wavelengths, the GVD increase suddenly becomes much more important because of an anticrossing point between the fundamental core mode and cladding modes at 1535 nm, delimiting the practical PBG transmission window [12]. In this region near the PBG edge, the GVD reaches extreme values well higher than 1000 ps/nm/km: measurements (not displayed in Figure 3 for the sake of clarity) indicates a GVD value 4760 ps/nm/km at 1535 nm, in excellent agreement with the calculated curve. This leads the thirdorder dispersion to reach extremely high values in the order of s^{3}/m in this spectral region. Finally, we also studied numerically the spectral dependance of the NL coefficient defined as [17], with the NL refractive index of pure silica and the mode effective area. Because the effective area also increases strongly near the PBG edge (due to the anticrossing phenomenon), the NL coefficient decreases rapidly near the PBG edge (see red curve in Figure 3) from about 6 W^{−1}km^{−1} at 1200 nm to less than 1.3 W^{−1}km^{−1} at 1535 nm.
In the following, we will see that the strong and unusual spectral dependencies of GVD and NL coefficient near the PBG edge can be used to tailor and manipulate the propagation of solitons approaching this spectral region through SSFS, included in the context of longpulse SC generation.
3. Quasicancellation of the Soliton SelfFrequency Shift
3.1. Experimental Evidence
To illustrate this possible manipulation of the SSFS, NL propagation experiments were first performed by launching 270 fs pulses with a repetition rate of 1 kHz at 1200 nm into a 5 m long sample of the solidcore PBG fiber 1. This pump wavelength has been chosen from a compromise between a wavelength not too close to the longwavelength edge of the first PBG (in order to observe a significant SSFS) and not too close to the ZDW (in order to avoid SC generation). Similarly, the input peak power was adjusted so that it was high enough to observe SSFS for several solitons while being low enough to avoid SC generation. The pump peak power effectively launched into the fiber core was estimated to about 10–20 kW. Figure 4 shows a representation of the output spectrum (in linear scale) as a function of fiber length from 0.4 to 5 m. It was obtained from a cutback measurement in which the output spectrum was recorded every 20 cm while keeping the launch conditions unchanged. The input Gaussian pulse breaks up into two fundamental solitons (labeled and ) through a fission mechanism [18, 19] caused by Raman scattering and higherorder dispersion effects [19, 20]. As expected, both solitons then redshift through Ramaninduced SSFS [21, 22]. The most striking feature in the soliton dynamics is that, after 2.5 m, the first ejected soliton, , becomes frequency locked at a central wavelength of about 1510 nm, that is, when its spectrum starts to overlap with the PBG located at 1535 nm (depicted in white dotted line in Figure 4). This is even more clearly illustrated in Figure 5(a) which shows a plot of the soliton central wavelength (black squares) as a function of the propagation distance. A clear slope change in the SSFS rate is observed at the fiber length of 2.5 m. Red squares in Figure 5(b) represent the evolution of soliton full width at half maximum (FWHM) spectral width versus propagation length. It is shown that, as soliton redshifts, it initially experiences a spectral compression from about 40 nm to 15 nm within the first 3 m. For longer propagation distances, that is, when the soliton is frequency locked, its spectral width stabilizes to about 15 nm, without any significant change in its spectral power density, as observed in Figure 4. Note that the same phenomenon was observed for the soliton (at a longer fiber length) and even for more than five ejected solitons by increasing the pump power. These experiments clearly show that the SSFS is extremely decelerated (and even quasicancelled) when the solitons approach the PBG edge, this effect being accompanied by a spectral compression.
3.2. Confirmation with Numerical Simulations
In order to get further insight into this phenomenon, we performed numerical simulations by integrating the generalized NL Schrödinger equation (GNLSE) taking into account the mode profile dispersion [23]. In this case, it writes [23, 24]: where the NL coefficient takes the following form [24]: and , being the complex spectral envelope. The coefficient corresponds to the mode propagation constant at the center frequency of the input pulse spectrum , and is the first derivative of (i.e., ). The symbol denotes the Fourier transform of the function in argument, and is the fiber linear power attenuation at . The silica response function is defined by where and is the Raman response function which is approximated by the analytical expression given in [25].
Figure 6 shows a propagation simulation over 5 m of a 270 fs Gaussian pulse with 15 kW peak power centered at 1200 nm, corresponding to the input experimental pulses. The frequency dependence of all parameters, that is, , , and , associated with the experimental fiber was taken into account. This simulation shows excellent qualitative agreement with the experimental results of Figure 4. The input pulse experiences fission leading to the ejection of two fundamental solitons, labeled to , and the SSSF associated to soliton is strongly decelerated near the PBG edge. This effect is accompanied by a spectral compression and an increase of its spectral power density. Black and red solid lines in Figure 5, respectively, correspond to the simulated soliton central wavelength and FWHM spectral width. They show good quantitative agreement with experiments, especially when the solitons approach the PBG edge (for propagation distances greater than 2.5 m). The blue solid curve in the inset of Figure 5 represents the soliton peak power as a function of propagation distance. Although its spectral power density increases (as shown in Figure 6), its peak power (in the temporal domain) decreases along propagation, because the soliton temporally broadens. The variation of total photon number as a function of fiber length is displayed in the green dashed curve in the inset of Figure 5. It decreases (due to linear attenuation) by only less than 1.5% within the first 5 m, which indicates that the soliton does not significantly suffer from the high attenuation at the PBG edge. This suggests that the quasicancellation of the SSFS does not occur at the expense of significant attenuation. At this point, however, the physics at the origin of this phenomenon has not been elucidated yet.
3.3. Discussion
First, let us recall that, as a soliton reaches the PBG edge, the values of the GVD evaluated at its central frequency strongly increase while its NL coefficient is largely reduced (see Figure 3). All these spectral evolutions are expected to be detrimental to the efficiency of the SSFS as they lead, directly or indirectly, to a significant increase of the pulse duration which is known to strongly affect the SSFS efficiency [17, 22]. In the following, we will present an analytical approach that allows to quantify the relative weight of the spectral dependence of the GVD () and NL coefficient () on the SSFS rate. For the sake of clarity, we only consider hereafter fundamental solitons as input pulses. This is justified by the fact that experimental input pulses experience fission which leads to pulse breakup and ejection of actual fundamental solitons [18–20]. The characteristics of these input pulses (hyperbolic secant pulses of 21 fs duration FWHM and 100 kW peak power) were estimated from the measurement of the soliton spectral width in Figure 5(b), assuming a soliton order .
The starting point of our analysis is the wellknown Gordon formula allowing to predict the SSFS rate, , of a fundamental soliton of mean frequency and duration as a function of propagation distance [22]. It writes where and is the imaginary part of the Fourier transform of the Raman response, which is related to the Raman gain spectrum [22]. Then, assuming that (i) the soliton remains fundamental, (ii) confinement losses are negligible, and (iii) the photon number is conserved, it is possible to find the following evolution of the soliton duration as a function of propagation distance [11]: where is the initial soliton temporal width. Figure 7 shows the evolution of the soliton central wavelength along propagation calculated by integrating (3) with help of (5) (red line) and by the numerical resolution of the GNLSE given by (1) (black line). A reasonable agreement is obtained between both curves, especially for the longest propagation distances (i.e., when the soliton appears to be frequency locked), the relative discrepancy on the frequency shift being less than 3% for values of larger than 1 m. Thus, this Figure 7 shows that (i) the SSFS deceleration near the PBG edge is directly linked to the fiber properties, as claimed above and that (ii) the soliton central frequency can be predicted easily and very accurately for the longest propagation distances, without requiring the resolution of the GNLSE.
Then, in order to distinguish the effect of the spectral dependence of and the one of the spectral dependence of , we now calculate the variation of the SSFS rate with propagation distance . Assuming that , that is, , this writes [11] where
From (6), it can be seen that the variation of SSFS rate during the propagation depends on the variation of and with . Since the fiber is assumed to be invariant along its length, the variations of and with are in fact only due to their spectral evolution as the soliton redshifts, that is, and . Consequently, the derivative of and with respect to in (6) can mathematically be expressed in terms of their derivative with respect to . We then introduce the ratio defined as so that values of higher than 1 indicate that the NL coefficient variation dominates the evolution of SSFS rate while values of smaller than 1 mean that the dynamics of the SSFS is mainly driven by the GVD variation , that is, by the thirdorder dispersion term .
The blue curve of Figure 7 (right axis) shows a plot of the ratio as a function of soliton wavelength calculated with (8). From this plot, it can be seen that the spectral variation of dominates the beginning of the SSFS dynamics, whereas for wavelengths higher than 1475 nm the effect induced by the spectral variation of becomes predominant. More importantly, the value just below the PBG edge (located at 1535 nm) is about 0.1, which clearly means that the extreme deceleration of the SSFS is due to the large GVD slope near the PBG edge rather than the large NL coefficient diminution.
Finally, let us mention that, in the absence of loss, a real cancellation of the SSFS, that is, , implies that and/or are infinite in order to get also equal to infinite. Thus, the strong decrease of the SSFS efficiency near the PBG edge does not lead, strictly speaking, to a true suppression of the SSFS for a finite propagation distance. However, the extreme deceleration observed all along this paper is enough to consider that the soliton frequency can be considered, at least from a practical point of view, as frequency locked.
4. Tailoring the SC Extension
The SSFS is one of the main mechanisms involved in the generation of octavespanning SC [26]. Indeed, in both the short and longpulse pumping regimes [27], the longwavelength part of the SC spectrum is composed of redshifting fundamental solitons initially created, respectively, by pulse fission and modulation instability (MI) [26]. Given that the spectral extent in the visible is intimately linked to the infrared one through a group index matching between solitons and trapped dispersive waves [28, 29], the SSFS dynamics lies at the heart of the SC generation process. There have been many attempts to optimize the spectral extension of SC sources, but, in most applications, only a small part of the SC spectrum is useful. There is thus an increasing interest in tailoring the spectral extent so that SC generation only occurs over the required spectral region. Fibers with two ZDWs have been employed in this way [30, 31]. Because the region of anomalous GVD can be suitably adjusted through a control of the microstructured cladding parameters, the SSFS can be cancelled near the longwavelength ZDW, which results in a stop of the infrared spectral broadening. However, the SSFS cancellation is accompanied by the emission of a dispersive wave at higher wavelengths, so that the energy of frequency locked solitons exponentially decreases. In this section, we show that the cancellation of the SSFS observed in solidcore PBG fibers described above allows to control the SC spectral extent without any significant loss of energy.
4.1. Dynamics of the SC Formation
SC experiments were performed using the solidcore PBG fiber 2 of Table 1, in which longwavelength PBG edge is located at 1580 nm. The pump laser used to generate the SC was a linearly polarized microchip laser emitting 0.6 ns pulses at 1064 nm with a repetition rate of 7 kHz. The coupling efficiency in the fibre core is estimated to be 50%. In the following, input power always refers to the power effectively launched into the fibre core.
At first, we investigated the dynamics of the SC formation with propagation distance by cutting the fiber back from 8 to 0.1 m and recording the output spectrum every 0.5 m, for a fixed pump power of 2.3 kW. Figure 8 shows the evolution of the output spectrum as a function of propagation distance. The spectrum initially presents two symmetric side lobes characteristic of MI classically observed when pumping in the anomalous GVD regime. In the temporal domain, this corresponds to the generation of a solitonic train, in which pulses strongly differ from each other in terms of peak power and duration. The spectrum then mainly broadens towards the infrared (because we pump relatively far away from the ZDW in the anomalous dispersion regime) until a fiber length of about 4 m is reached. At this propagation distance, the SC ranges from 1050 to 1580 nm. The infrared spectral broadening is due to the SSFS of each individual soliton generated from MI. Since they are characterized by very different properties, their SSFS efficiency can be very different, which leads to the SC. For further increasing fiber lengths to 8 m, the longwavelength spectral broadening completely stops, due to the SSFS cancellation of each soliton approaching the PBG edge, following the mechanism described in Section 3. This leads to an accumulation of solitons just below the PBG edge (located at 1580 nm), which can be seen as a spectral power peak centered around 1550 nm. Consequently, the SC longwavelength spectral extent can be tailored over a given spectral range through a suitable design of the fiber linear properties and, in particular, of the longwavelength PBG edge.
To further illustrate this, we performed an additional experiment in which the dynamics of the SC formation was investigated as a function of pump peak power, for a fixed fiber length of 10 m. To do that, the output spectrum was recorded with 0.2 kW increments for pump peak powers ranging from 0.2 to 2.6 kW. The result of these measurements is displayed in Figure 9. Similar to previous experiments of Figure 8, the spectrum mainly broadens towards the infrared when the pump peak power is increased from 0 up to about 1.7 kW. For higher pump powers, the spectral broadening suddenly stops at the longwavelength PBG edge. This is even more clearly illustrated in Figure 10(a) which represents the longwavelength SC edge (measured for a 10 dB drop with respect to the spectral power at 1100 nm) as a function of input peak power. A clear cancellation of the spectral broadening at the PBG edge (depicted in gray) is observed for peak powers higher than 1.7 kW. As discussed above, this effect is accompanied by an accumulation of solitons which can be seen as a spectral power peak at 1550 nm, that is, just below the PBG edge located at 1580 nm.
(a)
(b)
4.2. Discussion
As shown by Figure 4 and discussed in Section 3, the SSFS cancellation occurs without any reduction of the spectral power density. This thus suggests, in the frame of SC generation, that the SC spectrum width can be limited by the PBG effect without any significant power loss. In order to check this experimentally, we plotted in Figure 10(b) the output average power measured over the whole SC spectrum (red squares) as a function of the input peak power. These measurements can be linearly fitted (red line), and the slope does not change once the SSFS cancellation has occurred and the SC spectral broadening has stopped (corresponding to the green region). In order to check that this was not due to the residual pump peak in the SC, we performed the same kind of measurements after filtering it with a notch filter at 1064 nm at the fiber output. The results are displayed with blue markers in Figure 10 and show that the output SC power (without the contribution of the 1064 nm pump) still increases with the input peak power. This confirms that the tailoring of the SC extent is not done through a simple longwavelength filtering due to the PBG effect, which would lead to a significant output power reduction. In our case, the fact that solitons are frequency locked rather than attenuated near the PBG edge allows to keep most of the launched power into the SC and thus to optimize the pump power budget.
5. Reduction of PulsetoPulse Power Fluctuations at the SC LongWavelength Edge
In this section, we focus our attention on pulsetopulse fluctuations of the SC spectral power near the longwavelength edge of the spectrum. Indeed, in the longpulse regime, the spectral broadening originates from the MI process which is seeded by noise and the longwavelength SC edge is characterized by statistically rare optical rogue waves [32], originating from soliton dynamics [32–35]. Consequently, the SC is characterized by a low coherence, and pulsetopulse fluctuations of the spectral power at the longwavelength SC edge are very important [32, 36]. Here, we investigate the impact of the SSFS cancellation observed in solidcore PBG fibers on pulsetopulse fluctuations of the longwavelength SC edge.
5.1. Experiments
Experiments were performed using the solidcore PBG fiber 3 of Table 1, in which longwavelength PBG edge is located at 1590 nm. Following the results of Section 4, a SC was generated using a 7 m long fiber sample and a pump power of 3 kW. Figure 11(a) shows a measurement of the spectral broadening as a function of propagation distance (made using a cutback technique). As expected, the infrared spectral broadening stops near the PBG edge, depicted with the dashed white line. Figures 11(b1) and 11(b2) shows the spectra obtained at 7 m and 2.5 m, respectively. Although the spectral extent is comparable in both cases, the shape of the spectrum is different. Indeed, while the longwavelength edge is relatively smooth for a fiber length of 2.5 m, it becomes much sharper at 7 m, where the cancellation of the SSFS imposes a saturation of the spectral broadening. This phenomenon causes a soliton accumulation whose spectral signature is the peak located around 1560 nm, as discussed in the previous section.
(a)
(b)
(c)
(d)
Pulsetopulse power fluctuations of the corresponding SC obtained at 2.5 and 7 m were then studied experimentally. To do this, the SC output was collimated and spectrally filtered with a 10 nm bandpass filter centered at 1550 nm, that is, on the soliton accumulation peak. Pulsetopulse fluctuations in this 10 nm wide spectral regions were measured using the method firstly proposed in [37]. The energy of spectrally filtered pulses is proportional to their peak power, so that rogue events and their associated characteristics can be captured through a simplified measurement of shottoshot pulse energy [37]. This was done in our experiments using a photodiode and an analog oscilloscope, with the same photodiode maximum signal amplitude and trigger level for both fiber lengths. Note that, in our experiments, this setup provides a measurement of the average energy of each pulse, but it does not allow to quantify noisesensitive instabilities related to the actual temporal coherence of the SC, as explained in [38] for instance. The percentage of shottoshot variations, , was evaluated using where and are, respectively, the maximum and minimum photodiode signal amplitude measured at least for 10 of the 10,000 recorded pulses.
Figures 11(c1) and 11(c2) show samples of pulse trains measured for each investigated fiber length after the 10 nm bandpass filter at 1550 nm, from which pulsetopulse fluctuations can be evaluated. The pulse amplitude has been normalized to its average value over the 10,000 measurements for this plot. Pulsetopulse fluctuations calculated at 1550 nm are as high as 80% after 2.5 m while they drop down to 10% after 7 m. Besides, the corresponding statistical distributions at 1550 nm have very different shapes, as attested by Figures 11(d1) and 11(d2) which represent histograms calculated over 10,000 measurements for fiber lengths of 2.5 m and 7 m, respectively. The distribution is highly asymmetric after 2.5 m but becomes Gaussianlike at 7 m, that is, once solitons have reached the PBG edge. These observations thus suggest that there is a correlation between the cancellation of the SSFS occurring near the PBG edge (leading to a saturation of the SC spectrum for propagation distances longer than about 3 m), and the enhanced pulsetopulse stability and statistical distributions in the same spectral region.
5.2. Confirmation with Simulations
In order to get insight into these results, we performed numerical simulations using the GNLSE given by (1). At first, we simply attempted to reproduce experimental results of Figure 11. To do this, all available fiber characteristics and experimental conditions have been taken into account with no free parameter, with the exception of the pulse duration that has been reduced to 50 ps (against 600 ps in experiments). Indeed, shortening the pulse duration in simulations (while staying in the longpulse regime though) allows to significantly reduce the computation time without significantly affecting the SC dynamics. The input pump peak power was 3 kW, and quantum noise was modeled by adding one photon per mode with random phase on each spectral discretization bin of the input field [26]. In order to study the statistical behavior of the pulse train, we performed 200 simulations with random initial noise conditions.
Simulations results are displayed in Figure 12, with the same organization as Figure 11 for easy comparison. Figures 12(a), 12(b1), and 12(b2), which show the dynamics of the SC formation, correspond to the averaged spectra over 200 simulations shots. The agreement with experiments is excellent, and the typical spectral features discussed above about the spectrum shape at its longwavelength edge are accurately reproduced. The simulated temporal properties shown in Figures 12(c1), 12(c2), 12(d1), and 12(d2) have been calculated using a numerical gaussian filter of 10 nm (FWHM) centered at 1562 nm, which corresponds to the maximum of the spectral power peak observed near the PBG edge (see plot 12(b1)). Taking the inverse Fourier transform gives the temporal profile of the SC filtered at the corresponding wavelength, whose average power has been calculated. Samples of the modeled pulse trains corresponding to 20 simulation shots are displayed in Figures 12(c1) and 12(c2), for, respectively, fiber lengths of 7 and 2.5 m. Corresponding statistical distributions calculated from all the 200 simulations are shown in Figures 12(d1) and 12(d2) and also show excellent agreement with experiments.
(a)
(b)
(c)
(d)
All these simulation results therefore allow to reproduce not only the main spectral features observed in the experiments (dynamics of the SC formation, typical spectral shape at the longwavelength edge), but also the measured statistical features of the pulse train (pulsetopulse fluctuations as a function of wavelength and propagation distance, shape of the distribution). At this point, however, numerical simulations do not allow to discuss on the physical mechanisms responsible for the enhanced pulsetopulse stability near the PBG edge.
5.3. Discussion
Numerical simulations using the GNLSE accurately reproduce our experiments and can consequently be used to discuss the physical mechanisms causing the enhanced stability of the SC longwavelength edge. To this end, the spectrotemporal representation has proved to be a powerful tool in studying the dynamics of SC generation [26]. Figures 13(a)–13(d) show the evolution of a simulated spectrogram for fiber lengths of 1, 2.5, 4, and 7 m.
(a)
(b)
(c)
(d)
Figure 13(a) shows how MI initially generates solitons that subsequently redshifts through Ramaninduced SSFS as can be seen from Figure 13(b). In such a longpulse pumping regime, the presence of statistically rare temporal events in SC experiments can find two complementary physical origins. Firstly, it can be due to a single high peak power soliton generated from MI for particular initial noise conditions [32, 33]. Because of its higher peak power, it experiences a more efficient SSFS than other solitons and, consequently, becomes statistically rare at highest wavelengths due to the low probability for these particular noise conditions to happen [32]. Note that, in this case, the requirements for a soliton to be statistically rare at the longwavelength SC edge, are only to experience a slightly more efficient SSFS than the other ones [39]. Secondly, following these early interpretations, it has been suggested that rare and brief events can arise from the collision of two or more solitons travelling with different group velocities [34, 35, 39, 40] because of the convective nature of the system [41]. Note that these two explanations are complementary for explaining longtail statistical distributions usually observed at the SC longwavelength edge [41].
In the case of the present experiments, however, this dynamics is strongly affected by the SSFS cancellation occurring near the PBG edge, as can be seen from Figures 13(c) and 13(d). Firstly, these figures show that, as long as solitons are frequency locked around 1562 nm, they all have very close group velocities at this location in the spectrum. Consequently, once they have reached the PBG edge, they cannot collide anymore, which prevents the formation of brief spikes [34] and leads to an enhanced pulsetopulse stability at the SC edge. Secondly, since all solitons whose SSFS has been cancelled near the PBG edge have very close characteristics, power fluctuations are reduced in this spectral region as compared to the usual case in which only a few solitons with slightly higher peak power are rare at the longwavelength edge [32, 33].
As a consequence, the reduction of pulsetopulse fluctuations at the longwavelength edge as a function of propagation distance is intimately linked to the cancellation of the SFSS discussed in Section 3. It is thus in fact due to the specific linear properties inherent to solidcore PBG fibers and, in particular, to the strong thirdorder dispersion near the PBG edge, as shown in Figure 3.
6. Conclusion
In conclusion, we showed that the SSFS and its dynamics can be strongly and easily manipulated thanks to the specific properties of the solidcore PBG fibers used in our various experimental and numerical studies. Indeed, we demonstrated that a quasicancellation of the SSFS could be obtained in this kind of fibers thanks to the large increase of observed on the longwavelength edge of the first PBG. We established, trough experimental and numerical studies, that this phenomenon could be achieved without significant power loss. Thanks to these results, we demonstrated and clearly interpreted that these fibers could be used to tailor efficiently the spectral extension of SC at no power cost. Moreover, we revealed and explained that the longwavelength edge of SC generated in this kind of fibers can present strong quantitative and qualitative differences in terms of temporal stability compared to SC generated in more conventional fibres: the power fluctuations in this spectral domain are indeed strongly reduced and have a Gaussianlike statistical distribution and not the more usual shaped statistics associated with the presence of optical rogue waves generally observed at the longwavelength edge of more conventional SC. We believe these results to be helpful for building SC sources for the next generation of timeresolved coherent antistokes Raman scattering (CARS) microscopes for instance [42].
Acknowledgments
The authors acknowledge Rémi Habert for experimental assistance and Karen Delplace for assistance in fiber fabrication. This work was partly supported by the French Ministry of Higher Education and Research, the NordPas de Calais Regional Council, and FEDER through the “Contrat de Projets Etat Région (CPER) 20072013" and the “Campus Intelligence Ambiante" (CIA).
References
 F. Luan, A. K. George, T. D. Hedley et al., “Allsolid photonic bandgap fiber,” Optics Letters, vol. 29, no. 20, pp. 2369–2371, 2004. View at: Publisher Site  Google Scholar
 A. Argyros, T. A. Birks, S. G. LeonSaval, C. M. B. Cordeiro, and P. S. J. Russell, “Guidance properties of lowcontrast photonic bandgap fibres,” Optics Express, vol. 13, no. 7, pp. 2503–2511, 2005. View at: Publisher Site  Google Scholar
 V. Pureur, A. Betourne, G. Bouwmans et al., “Overview on solid core photonic bandgap fibers,” Fiber and Integrated Optics, vol. 28, no. 1, pp. 27–50, 2009. View at: Publisher Site  Google Scholar
 G. Bouwmans, L. Bigot, Y. Quiquempois, F. Lopez, L. Provino, and M. Douay, “Fabrication and characterization of an allsolid 2D photonic bandgap fiber with a lowloss region ($<20$ dB/km) around 1550 nm,” Optics Express, vol. 13, pp. 8452–8459, 2005. View at: Google Scholar
 A. Fuerbach, P. Steinvurzel, J. A. Bolger, A. Nulsen, and B. J. Eggleton, “Nonlinear propagation effects in antiresonant highindex inclusion photonic crystal fibers,” Optics Letters, vol. 30, no. 8, pp. 830–832, 2005. View at: Publisher Site  Google Scholar
 B. Kibler, T. Martynkien, M. Szpulak et al., “Nonlinear femtosecond pulse propagation in an allsolid photonic bandgap fiber,” Optics Express, vol. 17, no. 12, pp. 10393–10398, 2009. View at: Publisher Site  Google Scholar
 A. S. Cerqueira Jr., C. M. B. Cordeiro, F. Biancalana et al., “Nonlinear interaction between two different photonic bandgaps of a hybrid photonic crystal fiber,” Optics Letters, vol. 33, no. 18, pp. 2080–2082, 2008. View at: Google Scholar
 P. D. Rasmussen, J. Laegsgaard, and O. Bang, “Degenerate four wave mixing in solid core photonic bandgap fibers,” Optics Express, vol. 16, no. 6, pp. 4059–4068, 2008. View at: Google Scholar
 A. Bétourné, Y. Quiquempois, G. Bouwmans, and M. Douay, “Design of a photonic crystal fiber for phasematched frequency doubling or tripling,” Optics Express, vol. 16, no. 18, pp. 14255–14262, 2008. View at: Publisher Site  Google Scholar
 V. Pureur and J. M. Dudley, “Nonlinear spectral broadening of femtosecond pulses in solidcore photonic bandgap fibers,” Optics Letters, vol. 35, no. 16, pp. 2813–2815, 2010. View at: Publisher Site  Google Scholar
 O. Vanvincq, A. Kudlinski, A. Bétourné, Y. Quiquempois, and G. Bouwmans, “Extreme deceleration of the soliton selffrequency shift by the thirdorder dispersion in solidcore photonic bandgap fibers,” Journal of the Optical Society of America B, vol. 27, no. 11, pp. 2328–2335, 2010. View at: Publisher Site  Google Scholar
 A. Bétourné, A. Kudlinski, G. Bouwmans, O. Vanvincq, A. Mussot, and Y. Quiquempois, “Control of supercontinuum generation and soliton selffrequency shift in solidcore photonic bandgap fibers,” Optics Letters, vol. 34, no. 20, pp. 3083–3085, 2009. View at: Publisher Site  Google Scholar
 O. Vanvincq, B. Barviau, A. Mussot, G. Bouwmans, Y. Quiquempois, and A. Kudlinski, “Significant reduction of power fluctuations at the longwavelength edge of a supercontinuum generated in solidcore photonic bandgap fibers,” Optics Express, vol. 18, no. 23, pp. 24352–24360, 2010. View at: Publisher Site  Google Scholar
 A. Bétourné, G. Bouwmans, Y. Quiquempois, M. Perrin, and M. Douay, “Improvements of solidcore photonic bandgap fibers by means of interstitial air holes,” Optics Letters, vol. 32, no. 12, pp. 1719–1721, 2007. View at: Google Scholar
 M. Perrin, Y. Quiquempois, G. Bouwmans, and M. Douay, “Coexistence of total internal reflexion and bandgap modes in solid core photonic bandgap fibre with intersticial air holes,” Optics Express, vol. 15, no. 21, pp. 13783–13795, 2007. View at: Publisher Site  Google Scholar
 M. Tateda, N. Shibata, and S. Seikai, “Interferometric method for chromatic dispersion measurement in a singlemode optical fiber,” IEEE Journal of Quantum Electronics, vol. 17, no. 3, pp. 404–407, 1981. View at: Google Scholar
 G. P. Agrawal, Nonlinear Fiber Optics, Academic Press, New York, NY, USA, 4th edition, 2007.
 E. M. Dianov, A. Y. Karasik, P. V. Mamyshev et al., “StimulatedRaman conversion of multisoliton pulses in quartz optical fibers,” JETP Letters, vol. 41, pp. 294–297, 1985. View at: Google Scholar
 Y. Kodama and A. Hasegawa, “Nonlinear pulse propagation in a monomode dielectric guide,” IEEE Journal of Quantum Electronics, vol. 23, no. 5, pp. 510–524, 1987. View at: Google Scholar
 P. Beaud, W. Hodel, B. Zysset, and H. P. Weber, “Ultrashort pulse propagation, pulse breakup and fundamental soliton formation in a singlemode optical fiber,” IEEE Journal of Quantum Electronics, vol. 23, no. 11, pp. 1938–1946, 1987. View at: Google Scholar
 F. M. Mitschke and L. F. Mollenauer, “Discovery of the soliton selffrequency shift,” Optics Letters, vol. 11, no. 10, pp. 659–661, 1986. View at: Google Scholar
 J. P. Gordon, “Theory of the soliton selffrequency shift,” Optics Letters, vol. 11, no. 10, pp. 662–664, 1986. View at: Google Scholar
 J. Laegsgaard, “Mode profile dispersion in the generalized nonlinear Schrödinger equation,” Optics Express, vol. 15, no. 24, pp. 16110–16123, 2007. View at: Publisher Site  Google Scholar
 J. C. Travers, M. H. Frosz, and J. M. Dudley, “Nonlinear fibre optics overview,” in Supercontinuum Generation in Optical Fibers, J. M. Dudley and J. R. Taylor, Eds., chapter 3, Cambridge University Press, Cambridge, UK, 2010. View at: Google Scholar
 D. Hollenbeck and C. D. Cantrell, “Multiplevibrationalmode model for fiberoptic Raman gain spectrum and response function,” Journal of the Optical Society of America B, vol. 19, no. 12, pp. 2886–2892, 2002. View at: Google Scholar
 J. M. Dudley, G. Genty, and S. Coen, “Supercontinuum generation in photonic crystal fiber,” Reviews of Modern Physics, vol. 78, no. 4, pp. 1135–1184, 2006. View at: Publisher Site  Google Scholar
 The longpulse pumping regime refers to cases in which the pump pulse duration $\u25b3T$ is much longer than the MI oscillation period $\u25b3{T}_{\text{MI}}$, given by $\u25b3{T}_{\text{MI}}=2\pi \sqrt{{\beta}_{2}/(2\gamma P)}$, with ${\beta}_{2}$ the secondorder dispersion coefficient, $\gamma $ the NL coefficient and P the pump peak power. Cases in which $\u25b3T$ is in the order or less than $\u25b3{T}_{\text{MI}}$ correspond to the shortpulse pumping regime.
 A. V. Gorbach and D. V. Skryabin, “Light trapping in gravitylike potentials and expansion of supercontinuum spectra in photoniccrystal fibres,” Nature Photonics, vol. 1, no. 11, pp. 653–657, 2007. View at: Publisher Site  Google Scholar
 J. M. Stone and J. C. Knight, “Visibly “white” light generation in uniform photonic crystal fiber using a microchip laser,” Optics Express, vol. 16, no. 4, pp. 2670–2675, 2008. View at: Publisher Site  Google Scholar
 A. Kudlinski, G. Bouwmans, Y. Quiquempois, and A. Mussot, “Experimental demonstration of multiwatt continuouswave supercontinuum tailoring in photonic crystal fibers,” Applied Physics Letters, vol. 92, no. 14, Article ID 141103, 2008. View at: Publisher Site  Google Scholar
 A. Kudlinski, G. Bouwmans, M. Douay, M. Taki, and A. Mussot, “Dispersionengineered photonic crystal fibers for CWpumped supercontinuum sources,” Journal of Lightwave Technology, vol. 27, no. 11, pp. 1556–1564, 2009. View at: Publisher Site  Google Scholar
 D. R. Solli, C. Ropers, P. Koonath, and B. Jalali, “Optical rogue waves,” Nature, vol. 450, no. 7172, pp. 1054–1057, 2007. View at: Publisher Site  Google Scholar
 J. M. Dudley, G. Genty, and B. J. Eggleton, “Harnessing and control of optical rogue waves in supercontinuum generation,” Optics Express, vol. 16, no. 6, pp. 3644–3651, 2008. View at: Publisher Site  Google Scholar
 A. Mussot, A. Kudlinski, M. Kolobov, E. Louvergneaux, M. Douay, and M. Taki, “Observation of extreme temporal events in CWpumped supercontinuum,” Optics Express, vol. 17, no. 19, pp. 17010–17015, 2009. View at: Publisher Site  Google Scholar
 G. Genty, C. M. de Sterke, O. Bang, F. Dias, N. Akhmediev, and J. M. Dudley, “Collisions and turbulence in optical rogue wave formation,” Physics Letters, Section A, vol. 374, no. 7, pp. 989–996, 2010. View at: Publisher Site  Google Scholar
 F. Vanholsbeeck, S. MartinLopez, M. GonzálezHerráez, and S. Coen, “The role of pump incoherence in continuouswave supercontinuum generation,” Optics Express, vol. 13, no. 17, pp. 6615–6625, 2005. View at: Publisher Site  Google Scholar
 C. Lafargue, J. Bolger, G. Genty, F. Dias, J. M. Dudley, and B. J. Eggleton, “Direct detection of optical rogue wave energy statistics in supercontinuum generation,” Electronics Letters, vol. 45, no. 4, pp. 217–219, 2009. View at: Publisher Site  Google Scholar
 H. Kubota, K. R. Tamura, and M. Nakazawa, “Analyses of coherencemaintained ultrashort optical pulse trains and supercontinuum generation in the presence of solitonamplified spontaneousemission interaction,” Journal of the Optical Society of America B, vol. 16, no. 12, pp. 2223–2232, 1999. View at: Google Scholar
 M. Erkintalo, G. Genty, and J. M. Dudley, “On the statistical interpretation of optical rogue waves,” European Physical Journal, vol. 185, no. 1, pp. 135–144, 2010. View at: Publisher Site  Google Scholar
 K. Hammani, B. Kibler, C. Finot, and A. Picozzi, “Emergence of rogue waves from optical turbulence,” Physics Letters, Section A, vol. 374, no. 34, pp. 3585–3589, 2010. View at: Publisher Site  Google Scholar
 M. Taki, A. Mussot, A. Kudlinski, E. Louvergneaux, M. Kolobov, and M. Douay, “Thirdorder dispersion for generating optical rogue solitons,” Physics Letters, Section A, vol. 374, no. 4, pp. 691–695, 2010. View at: Publisher Site  Google Scholar
 D. A. SidorovBiryukov, E. E. Serebryannikov, and A. M. Zheltikov, “Timeresolved coherent antiStokes Raman scattering with a femtosecond soliton output of a photoniccrystal fiber,” Optics Letters, vol. 31, no. 15, pp. 2323–2325, 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 O. Vanvincq et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.