# Design of solid core photonic bandgap fibers for visible supercontinuum generation

V. Pureur, J.M. Dudley

► To cite this version:

V. Pureur, J.M. Dudley. Design of solid core photonic bandgap fibers for visible supercontinuum generation. Optics Communications, 2011, 284 (6), pp.1661-1668. 10.1016/j.optcom.2010.11.037. hal-00582921

HAL Id: hal-00582921

https://hal.science/hal-00582921v1

Submitted on 6 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Design of solid core photonic bandgap fibers for visible supercontinuum generation

Vincent Pureur*, John M. Dudley

Institut FEMTO-ST, Département d'Optique P. M. Duffieux, CNRS UMR 6174, Université de Franche-Comté, 25030 Besançon, France

We present a numerical study of two dimensional solid core photonic bandgap fiber design criteria for their particular application to blue/visible supercontinuum generation. By exploiting their strong frequency-dependent dispersion when compared to index guiding micro-structured fibers, we highlight the design of solid core photonic bandgap fibers to fulfill group index matching conditions between the first ejected optical soliton and the trapped dispersive wave generated in the visible wavelength range. We study how these matching conditions depend on the opto-geometrical parameters of the micro-structured cladding, and we use frequency-domain numerical simulations to determine the expected supercontinuum spectral characteristics for selected cases. We investigate design criteria to generate short wavelengths by pumping in such photonic bandgap fibers in different pulse duration regimes and we identify a novel class of short wavelength (blue/visible) supercontinuum generation in the 3rd bandgap of a typical structure by pumping into the 2nd bandgap through a high attenuation spectral region.

## 1. Introduction

The generation of broadband supercontinuum spectra has attracted much interest over the last decade since the advent of photonic crystal fiber (PCF) [1–5]. Taking advantage of the novel dispersive and non linear properties [6], the generation of new frequencies through supercontinuum generation in PCF has found many applications in source development, optical coherence tomography, spectroscopy, confocal microscopy and optical frequency metrology.

An important area of the current activity is the study of techniques to extend the non linear spectral broadening further into the visible or the blue spectral range. One approach that has been extensively studied is the mechanism whereby a dispersive wave (DW) generated in the normal group velocity dispersion (GVD) regime is trapped by a soliton propagating in the infrared in the anomalous GVD regime such that there is enhanced frequency shifting to shorter wavelengths [7–10]. In order to achieve this, group index matching conditions must be fulfilled between the dispersive wave and the soliton, and this effect has been observed experimentally in the visible range in index guiding PCFs both in the picosecond and CW pumping regime pumping around 1060–1080 nm [11–13].

The extent of the supercontinuum in the infrared is limited in silica PCFs by the increase of OH⁻ ions absorption beyond 2.4 µm with a strong peak around 1.38 µm, limiting then the coupled energy to the trapped dispersive wave in the visible [13]. In order to improve short

wavelength component generation in PCF, one way is to decrease water absorption around 1.38 µm from typically 600 dB/km to about 90 dB/km by reducing surface contamination during the fabrication process through chemical cleaning [14]. However, another way is to modify the group index curve of the fundamental core guided mode by tailoring the geometry of the photonic crystal cladding and the core size and by keeping the zero dispersion wavelength (ZDW) close enough to the pump. Because the group index on the short wavelength side is mainly dominated by the material dispersion (silica), efforts have been recently made to optimize the waveguide group index in the infrared of index guiding PCFs to push the supercontinuum further into the near ultraviolet [13]. However in this case the optimization also shifts the ZDW further away from the pump, limiting the efficiency of the short wavelength generation in the quasi-CW regime [15].

In this work we consider how the novel guidance characteristics of two dimensional solid core photonic bandgap (PBG) fibers can overcome this limitation, and be used to generate short wavelengths via supercontinuum generation. In fact, recent experimental work has reported spectral broadening in PBG fibers in the fundamental bandgap using sub-nanosecond pulses at 1064 nm [16] but in the infrared from 1 µm to 1.6 µm and not in the visible. We show here that group index matching conditions can be found between dispersive waves in the blue/visible and optical solitons propagating in the infrared without reaching wavelength above 1 µm. In addition, motivated by recent numerical studies showing energy transfer between two different transmission bands across a spectral region of high attenuation [17], we also simulate trapped dispersive wave generation in the blue edge of the 3rd bandgap of a typical structure by pumping into the anomalous GVD regime of the 2nd transmission

<footnote>* Corresponding author.
E-mail address: Vincent.Pureur@femto-st.fr (V. Pureur).</footnote>
band in order to highlight the potential of PBG fibers to generate blue/ visible components.

## 2. Solid core photonic bandgap fibers: group index matching

Two dimensional solid core PBG fibers propagate light in a low index defect through a PBG guidance mechanism [18]. In order to obtain such an effect, high refractive index inclusions (Germanium/ fluorine doped-silica or oil-infiltrated air holes) are periodically embedded into a low index background made generally of pure silica as the core (see Fig. 1) [18]. This guidance mechanism involves propagation of the light through different spectral windows of different order (different PBGs) each associated with its own linear GVD and confinement losses and non linear properties; both the linear and non linear properties are strongly frequency-dependent [18].

According to the Anti Resonant Reflecting Optical Waveguide (ARROW) approach [19], the spectral position of the different transmission windows depend on the well-known normalized frequency (V-parameter) of an isolated high index rod, and depends on its diameter *d* and the refractive index contrast $\Delta n$ with the low index medium. Theoretically the PBGs can be shifted spectrally by changing only one of these two opto-geometrical parameters, but the confinement and bending losses can be limiting factors [20] as they also depend strongly on the *d/λ* ratio and $\Delta n$, where *λ* is the rod-to-rod space and on the micro-structured cladding size relative to the wavelength (proportional to the ratio *λ/λ* [18]. The opto-geometrical parameters of the micro-structure need then to be carefully optimised in order to obtain simultaneously low loss and a desired spectral position of the different transmission bands.

To highlight the potential of solid core PBG fibers for super- continuum (SC) generation [21] toward short wavelengths, we firstly consider a PBG structure with a low index core surrounded by four confinement rings of high index rods (Cargille oil, *n* = 1.6 at *λ* = 589.3 nm, step refractive index profile). Pitch and diameter are 3.10 µm and 1.23 µm respectively. Using a multipole method [22], we calculate the dispersive and non linear guidance properties of the fundamental core guided mode. For both the silica and inclusion the material dispersion is included; we neglect material absorption over the wavelength range considered here. Fig. 2 plots the PBG guidance characteristics obtained as a function of the normalized frequency *V* of an isolated high index rod. The top panel shows the linear guidance properties: confinement losses (left axis) and group index curves (right axis). The bottom panel plots the effective area (left axis) and the associated non linear parameter γ (right axis). In order to illustrate the interest of solid core PBG fibers compared to index guiding PCF, we also plot on Fig. 2 the group index (top panel, dashed line) evolution of the

fundamental mode guided into the index guiding PCF with a silica core surrounded by air holes (dimensions as above, with material dispersion).

Due to energy coupling from the core guided mode to cladding modes close to each PBG, the confinement losses, group index and effective area of the core guided mode all undergo dramatic increase contrary to index guiding PCFs. These strong variations of group index curves in solid core PBG fibers close to each bandgap edges can be used to fulfill matching conditions between the trapped dispersive wave and the optical soliton in the infrared.

## 3. Group index curve optimization and supercontinuum generation

In this section we use a finite element method (COMSOL software) to calculate the linear and non linear properties of the fundamental core guided mode for a wider variety of micro-structure configurations considered as mono-modes [23]. The confinement losses are computed using a circular Perfectly Matched Layer (PML) and the material dispersion is always taken into account (silica or doped silica). For the purposes of our simulations, we choose to work with a step refractive index profile of the high index inclusions in the cladding, the core being one missing rod. In order to give general guidelines of group index optimization we work in the first PBG order which is spectrally the widest but also characterized by the highest confinement losses. Section 3.1 considers the effect of varying $\Delta n$ and Section 3.2 considers the effect of varying *d/λ*. Note that the dispersion characteristics depend only weakly on the number of confinement rings considered. On the other hand, the dispersion and non linearity are very strongly affected by the particular opto-geometric parameters of pitch, diameter and $\Delta n$ which are the focus of our study here. In each of these sections we present explicit simulations of the supercontinuum generation dynamics to link the spectral broadening mechanisms with the particular choice of fiber parameters.

### 3.1. Varying $\Delta n$

We plot in Fig. 3 the dependence on wavelength of the (a) effective index, (b) confinement losses, (c) group index ($n_g$) and (d) non linear parameter $\gamma$ of the fundamental core guided mode for different micro- structure configurations with an increase of $\Delta n$ from 0.01 to 0.1. The ratio *d/λ* is fixed to 0.4 in all cases and six confinement rings are simulated. Even if such high refractive index contrasts are impossible to achieve in silica (all solid structures), the infiltration of an air hole PCF with high refractive index fluids can lead to such high $\Delta n$ values [24]. In order to keep the same spectral position of PBG1 (ARROW

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** (a) Schematic representation of a two dimensional solid core PBG fiber. The dark regions represent the high refractive index inclusions embedded in a pure silica background (light regions). (b) Transversal cut showing a step refractive index profile in the inclusions.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Top axes: confinement losses (red, left) and group index curves (black, right) as a function of normalized frequency *V*. Bottom axes: effective area (blue, left) and non linear parameter *γ* (green, right). The dashed line represents group index for an index-guiding PCF. The vertical lines show lowest-order guided mode cut-offs delimiting approximately the bandgaps.</figcaption></figure>

model [19]), the rod diameter (and thus $\Lambda$) is decreased for each $\Delta n$
value. For example $\Lambda$ goes from 3.9 µm to 1.22 µm.

In Fig. 3(a) and (c), we firstly observe that increasing $\Delta n$ involves a decreasing of the fundamental mode effective index, due to deeper PBG (in effective index) [18]. Consequently, the group index value rises with $\Delta n$ in the center of the PBG, and $n_g$ increases quicker to the long wavelength edge with $\Delta n$. This steep increase of the group index in the infrared is a characteristic of the PBG guidance mechanism and is not observed in the modified total internal reflection PCFs. Moreover we observe that the minimum group index related to the ZDW shifts toward the short wavelength side when increasing $\Delta n$ (from roughly 845 nm to 570 nm).

Fig. 3(b) and (d) show that an increase of $\Delta n$ is followed by a rise of the non linear coefficient and confinement losses over all the PBG. Indeed, when $\Delta n$ increases, $\Lambda$ must be decreased to keep the PBG at the same wavelength. Then the ratio $\Lambda/\lambda$ (size of the micro-structured cladding relative to the wavelength) decreases and leads to a rise of confinement losses and $\gamma$ with $\Delta n$. However (and in contrast to the group index curves variation with $\Delta n$ in Fig. 3(c)), Fig. 3(b) and (d) show that the spectral minima of confinement losses and maxima of non linear parameter in the PBG depend slightly on $\Delta n$. Then for this $d/\Lambda$ value, the shift of the ZDW highlighted in Fig. 3(c) toward the infrared for weak $\Delta n$ is followed by a strong decrease of the light confinement and non linearity. For example at the ZDW the confinement losses go

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Spectral evolution of the fundamental mode effective (a) and group index (c), confinement losses (b) and non linear parameter $\gamma$ (d) in PBG1 for different $\Delta n$ values.</figcaption></figure>
from $2 \times 10^4$ dB/km to $4.5 \times 10^5$ dB/km when $\Delta n$ decreases from 0.1 to 0.01 and $\gamma$ goes from $100 (\text{W km})^{-1}$ to $3.8 (\text{W km})^{-1}$.

### 3.1.1. Supercontinuum generation for fixed soliton order at constant $\beta_2$

To model non linear propagation of femtosecond pulses in PBG fibers we use the rigorous frequency-domain approach described by Laesggaard [25] using Raman parametrization as in Ref. [26]. This approach allows an approximation-free treatment of the strong frequency-dependent loss, dispersion and non linearity of PBG fibers [17]. To highlight the $\Delta n$ influence on the short wavelength generation, we simulate non linear spectral broadening in PBG1 in the different structures presented in Fig. 3 for a constant soliton order $N^2 = \gamma P_0 T_0^2 / |\beta_2|$. For each set of parameter, the pump wavelength is changed so that the GVD is constant at $\beta_2 = -10 \text{ ps}^2/\text{km}$ (anomalous regime). The temporal pulse width being held constant at 50 fs, the input peak power $P_0$ is changed to keep the product $P_0 \gamma$ constant at the pump wavelength ($N \approx 10$). The computed non linear length at the pump is close to 0.7 mm for all the cases.

For the lowest value of $\Delta n$ (0.01) and for a $d/\Lambda$ of 0.4, Fig. 3(c) shows that the increase of $n_g$ to the infrared is not sufficient to generate a trapped dispersive wave on the blue side before reaching high confinement losses close to the PBG edge. Fig. 4 shows then simulated SC evolutions for the three different $\Delta n$ values 0.1 (a), 0.05 (b) and 0.03 (c). We observe in Fig. 4(a) ($P_0 \sim 16.5$ kW) that a coupling between the soliton and the phase matched dispersive wave (around 450 nm) occurs around the first millimetres of propagation while the first ejected optical soliton is largely attenuated after a propagation distance of 10 cm. When $\Delta n$ decreases, Fig. 4(b) and (c) show a shift of the whole SC toward the infrared where losses are higher. Consequently the first ejected optical soliton is attenuated over a shorter distance (around 4 cm and 1 cm for

$\Delta n = 0.05$ and 0.03 respectively) and the coupling efficiency with the red shifted dispersive wave is limited. Moreover small $\Delta n$ values involve higher input peak power ($P_0 \sim 100$ kW for $\Delta n = 0.03$) due to reduced non linearity.

### 3.1.2. Supercontinuum generation for a fixed input pulse wavelength

In this section we propose to highlight the $\Delta n$ influence on the non linear spectral broadening when the input pulse wavelength is fixed. As in section 3.1.1 the soliton order at the pump wavelength is kept constant through the ratio $P_0 \gamma / |\beta_2|$ ($N \approx 8$). Fig. 5 shows then the supercontinuum evolution over 20 cm of fiber when pumped at 650 nm for the different $\Delta n$ values: 0.03 (a), 0.05 (b) and 0.1 (c).

By changing $\Delta n$ and keeping the same input pulse wavelength, we modify $\beta_2$ and the GVD regime. For $\Delta n = 0.1$, $\beta_2$ is negative at 650 nm ($\sim -20 \text{ mps}^2/\text{km}$, anomalous regime). Fig. 5(a) shows the non linear spectral broadening occurring for an input peak power of 20 kW. We clearly observe after a few millimetres of propagation the generation of a dispersive wave around 400 nm phase matched with optical soliton propagating in the infrared. When $\Delta n = 0.05$, $\beta_2$ is very weak and close to the ZDW ($\sim -1 \text{ ps}^2/\text{km}$). Because the required input peak power is low at $\sim 2.3$ kW, Fig. 5(b) does not show non linear spectral broadening. When $\Delta n$ decreases to 0.03, the GVD regime becomes normal ($\beta_2 \sim +5.5 \text{ ps}^2/\text{km}$ at 650 nm and $P_0 \sim 18$ kW). Fig. 5(c) shows the limited spectral broadening due to self phase modulation. In this regime the non linear processes are not dominated by soliton fission and dispersive waves are not generated on the short wavelength side. Note that the conversion efficiency from pump to short wavelength dispersive wave is here comparable to that seen in the standard index guiding PCF ($< 1\%$) but with the strong pump powers used, useful intensities of short wavelength light can nonetheless be generated.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Simulated SC evolutions in PBG1 over a fiber length of 20 cm for $\Delta n = 0.1$ (a), 0.05 (b) and 0.03 (c). $\beta_2$ is fixed at $-10 \text{ ps}^2/\text{km}$.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Simulated SC evolutions in PBG1 over a fiber length of 20 cm for $\Delta n = 0.1$ (a), 0.05 (b) and 0.03 (c). The input pulse wavelength is fixed at 650 nm.</figcaption></figure>

### 3.2. Varying $d/\Lambda$

We plot in Fig. 6 the effective (a) and group index (c) variations of the fundamental core guided mode, its confinement losses (b) and its non linear parameter $\gamma$ (d) versus wavelength in the first PBG for the different micro-structure configurations with $d/\Lambda$ ratios of 0.2, 0.4 and 0.6. The refractive index contrast $\Delta n$ is set to 0.03 in all cases and seven confinement rings are simulated. As the spectral position of the transmission windows depends only on $d$ and $\Delta n$ [19], only the pitch

of the micro-structure is changed from 4.5 µm to 1.5 µm. In Fig. 6(a) and (c), we observe that increasing the ratio $d/\Lambda$ leads to a decrease of the fundamental mode effective index, due to deeper PBG [18]. Consequently the group index value rise with $d/\Lambda$ around the center of the PBG and $n_g$ increases quicker with $d/\Lambda$ in the infrared. Moreover, the ZDW shifts toward the short wavelength when increasing $d/\Lambda$ (from roughly 895 nm to 600 nm).

Fig. 6(b) and (d) show that an increase of $d/\Lambda$ is followed by an increase of $\gamma$ (around the center of the PBG) and of the confinement

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Spectral evolution of the fundamental mode effective (a) and group index (c), confinement losses (b) and non linear parameter $\gamma$ (d) in PBG1 for different $d/\Lambda$ ratio.</figcaption></figure>
losses. Indeed, the rise of the *d*/Λ ratio involves a decrease of Λ (*d* being constant). The ratio Λ/Λ decreases and leads to an increase of confinement losses and γ with *d*/Λ. For the smallest *d*/Λ ratio, the high index inclusions are widely spaced apart. The confinement losses characteristics of the core guided mode correlate with the scattering properties of a single cylinder which are Fano resonance [27]. Due to these properties the confinement losses of the fundamental mode in the first PBG can be very low in the center of the transmission band even with a few rings. Fig. 6(b) and (d) highlight moreover that the loss minimum and non linear parameter maximum shift slightly toward long wavelength when increasing *d*/Λ (comparing 0.4 and 0.6) whereas the *d*/Λ=0.2 case is not comparable regarding the loss.

3.2.1. Supercontinuum generation for a fixed β₂

We simulate supercontinuum generation in PBG1 for the different *d*/Λ presented in Fig. 6 and for a constant soliton order. Again the pump wavelength is changed such that β₂ is fixed to −10 ps²/km (anomalous regime) and the input peak power P₀ is changed to keep the product P₀γ constant at the pump wavelength (N≈5).

Fig. 7 shows the simulated SC evolution over 20 cm. We observe
that the shift in the pump wavelength associated with maintaining
constant β₂ and with increasing d/Λ leads to a shift of the generated SC
toward the short wavelength side as the ZDW decreases (P₀ = 100 kW
for d/Λ= 0.2). Consequently the generated DW phase matched with
the optical soliton reaches shorter wavelength with a high d/Λ value,
but the SC becomes spectrally narrower according to the “U” shaped
form of the group index presented in Fig. 6(c). Moreover the shift of
the ZDW to short wavelength when d/Λ increases enables propagation

in a regime of low confinement losses as the loss minimum is closer to
the short wavelength side in PBG1. The required input peak power is
also reduced at ~10 kW for d/Λ=0.6 because of the higher non linear
parameter. However for d/Λ=0.6 (Fig. 7(c)) the confinement losses
are very high both at the dispersive wave and the optical soliton
wavelengths and thus the coupled energy to the visible decreases
rapidly with fiber length.

3.2.2. Supercontinuum generation for a fixed input pulse wavelength

In this section we investigate the *d*/Λ influence on the non linear
spectral broadening when the input wavelength is fixed. As in the
previous section the soliton order is constant at the pump wavelength
through the ratio P₀γ/|β₂| (N≈5). The simulated structures shown in
Figs. 6 and 8 represent the SC evolutions when pumped at 750 nm for
the different *d*/Λ values: 0.2 (a), 0.4 (b) and 0.6 (c).

As the ZDW changes with the ratio d/Λ by keeping the input pulse
wavelength constant we modify the GVD regime. For d/Λ=0.2, β₂ is
positive at 750 nm (~+26 ps²/km). In the normal GVD regime Fig. 8
(a) shows the initial and symmetric non linear spectral broadening
mainly due to self phase modulation. Optical solitons are not
generated in the normal regime and neither DW in short wavelengths.
For d/Λ=0.4, β₂ is close to the ZDW and in the anomalous GVD regime
(~−19 ps²/km). Fig. 8(b) put in evidence generation of optical soliton
phase matched with DW around 550 nm. When d/Λ=0.6, β₂
decreases to ~−116 ps²/km. Whereas DW generation at shorter
wavelengths is expected, Fig. 8(c) shows that the level of losses
around 750 nm is too high and prevent non linear spectral broadening
and soliton fission even after a few millimetres.

<figure><img src="image_8.png" /><figcaption>Fig. 7. Simulated SC evolutions in PBG1 over a fiber length of 20 cm for d/Λ=0.2 (a), 0.4 (b) and 0.6 (c). β₂ is fixed at −10 ps²/km.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>Fig. 8. Simulated SC evolutions in PBG1 over a fiber length of 20 cm for $d/λ=0.2$ (a), 0.4 (b) and 0.6 (c). The input pulse wavelength is fixed at 750 nm.</figcaption></figure>

### 3.3. Discussion

Numerical results of Sections 3.1 and 3.2 show that when the PBG become deeper in effective index, increasing for example the refractive index contrast $\Delta n$ and/or the ratio $d/λ$ of the micro-structured cladding, the group index curve of the fundamental core guided mode rises mainly on the long wavelength side, while the short wavelength remains dominated by material dispersion in the spectral range of interest. Figs. 3(c) and 6 (c) show that the frequency dependence of group index can be stronger than in index guiding PCFs and highlight for the first time the possibilities of group index matching without reaching a soliton wavelength above 1 µm in the infrared.

Additionally Figs. 3(b) and 6 (c) show that such high values of $\Delta n$ and $d/λ$ involve lower ZDW and high confinement losses that can represent a limiting factor. Indeed the results highlighted in Sections 3.1 and 3.2 show that the “U” form of the group index can improve the shift of the generated dispersive wave toward the blue for high $\Delta n$ and/or $d/λ$ values, but the losses prevent the optical soliton propagation in the infrared and thus limit the coupled energy to the dispersive wave. Consequently a compromise must be found between rising of the group index in the infrared and weak losses. A careful design must be undertaken to decrease the losses in the fundamental bandgap by increasing for example the number of confinement rings. The adding of an external ring of air holes around the micro-structured cladding [28] or interstitial air holes between each rod [29] is well known to reduce strongly both confinement and bending losses. With sufficient low losses the use of solid core PBG fibers can lead to the generation of short wavelength components without reaching wavelength above 1 µm in the infrared.

### 4. Energy transfer across different PBG orders

In this section we propose to demonstrate numerically short wavelengths generation in the blue pumping into the anomalous GVD

of a typical solid core PBG structure in the femtosecond regime. In order to highlight the potential of such PBG fibers on the non linear spectral broadening of ultra-short pulses, we choose to generate a dispersive wave across a spectral region of high attenuation pumping into an adjacent PBG [17].

We consider a PBG structure with a low index core (pure silica) surrounded by four confinement rings of high index rods (Cargille oil, $n=1.6$ at $λ=589.3$ nm). Pitch and $d/λ$ ratio are 2.44 µm and ~0.4 respectively. All the linear and non linear optical properties of the fundamental core guided mode are simulated using the multipole method (with material dispersion) [22]. We plot in Fig. 9 the group

<figure><img src="image_10.png" /><figcaption>Fig. 9. Spectral evolution of the group index (a) and confinement losses (b) of the fundamental core mode within the 2nd and 3rd PBG order.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>Fig. 10. (a) SC evolution in PBG2 and PBG3 over fiber length of 50 cm ($P_{in}$ = 100 kW, $\lambda_{pump}$ = 700 nm). (b) Normalized output spectrum after 50 cm of propagation.</figcaption></figure>

index and confinement losses spectral variations of this mode confined within the second and third PBG. We observe potential group index matching (dashed lines) between the long wavelength side of the 2nd PBG (where the optical soliton propagates when pumping into its anomalous GVD regime) and the blue edge of PBG3 around 420 nm where a dispersive wave can be generated.

Pumping at 700 nm with an input peak power of 100 kW and 50 fs pulse duration (anomalous GVD regime), Fig. 10 highlights the energy transfer from the first ejected optical soliton propagating in PBG2 to a dispersive wave generated in the blue side of PBG3. Even if the fundamental soliton in PBG2 is attenuated after a propagation distance of 20 cm, the energy transfer occurs over the first millimetres of propagation because the computed non linear length at the pump is close to 0.4 mm. Then once generated in the short wavelength in low loss regions of PBG3 the dispersive wave see much lower attenuation and can propagate over a longer distance. Fig. 10(b) plots the normalized output spectrum simulated after a propagation distance of 50 cm. To the best of our knowledge this is the first demonstration of the possibility of such short wavelength generation in PBG fibers. As the spectral position of the different transmission bands can be shifted by varying the rod diameter or the refractive index contrast $\Delta n$ of the micro-structured cladding, a slight reduction of the rod-to-rod space ($d/\Lambda$fixed) or $\Delta n$ (using another oil) can shift the generated dispersive wave below 400 nm. This shift of the bands toward the blue is also associated with an increase of fundamental mode confinement losses which can be decreased by increasing the number of confinement rings for example.

## 5. Conclusion

In this paper, we have studied numerically the use of solid core photonic bandgap fibers for visible supercontinuum generation. The strong frequency-dependent group index curve of the fundamental mode close to the bandgap edge allows dispersive wave generation in short wavelengths without the need to reach soliton wavelengths above 1 µm as is the case with modified total internal reflection PCF. We numerically investigate the impact of opto-geometrical parameters of the micro-structured cladding on supercontinuum generation in the visible, such as the rod-to-rod space and the refractive index contrast between the inclusion and the low index medium. Design criteria are highlighted to generate efficiently visible wavelength which have to be associated with sufficient low confinement losses. Moreover taking into account all the strong frequency-dependent optical properties of such structures we demonstrate the generation of short wavelengths on the

blue edge of a bandgap pumping into its adjacent bandgap through a spectral region of high attenuation.

## Acknowledgment

Authors thank M. Abdilahi for the helpful discussions. We acknowledge support from the French Agence Nationale de la Recherche projects MANUREVA ANR-08-SYSC-019 and IMFINI ANR-09-BLAN-0065.

## References

[1] J.K. Ranka, R.S. Windeler, A.J. Stentz, Opt. Lett. 25 (2000) 25.

[2] S. Coen, A.H.L. Chau, R. Leonhardt, J.D. Harvey, J.C. Knight, W.J. Wadsworth, P.St.J. Russell, Opt. Lett. 26 (2001) 1356.

[3] T.A. Birks, W.J. Wadsworth, P.St.J. Russell, Opt. Lett. 25 (2000) 1415.

[4] J.M. Dudley, G. Genty, S. Coen, Rev. Mod. Phys. 78 (2006) 1135.

[5] J.M. Dudley, J.R. Taylor, Nat. Photonics 3 (2009) 85.

[6] J.C. Knight, T. Birks, P. Russell, D. Atkin, Opt. Lett. 21 (1996) 1547.

[7] G. Genty, M. Lehtonen, H. Ludvigsen, Opt. Express 12 (2004) 4614.

[8] A.V. Gorbach, D.V. Skryabin, J.M. Stone, J.C. Knight, Opt. Express 14 (2006) 9854.

[9] A.V. Gorbach, D.V. Skryabin, Nat. Photonics 1 (2007) 653.

[10] J.C. Travers, A.B. Rulkov, B.A. Cumberland, S.V. Popov, J.R. Taylor, Opt. Express 16 (2008) 14435.

[11] A. Kudlinski, A.K. George, J.C. Knight, J.C. Travers, A.B. Rulkov, S.V. Popov, J.R. Taylor, Opt. Express 14 (2006) 5715.

[12] A. Kudlinski, A. Mussot, Opt. Lett. 33 (2008) 2407.

[13] J.M. Stone, J.C. Knight, Opt. Express 16 (2008) 2670.

[14] A. Kudlinski, G. Bouwmans, M. Douay, M. Taki, A. Mussot, J. Lightwave Technol. 27 (2009) 1556.

[15] W.J. Wadsworth, N. Joly, J.C. Knight, T.A. Birks, F. Biancalana, P.St.J. Russell, Opt. Express 12 (2004) 299.

[16] A. Bétourné, A. Kudlinski, G. Bouwmans, O. Vancincq, A. Mussot, Y. Quiquempois, Opt. Lett. 34 (2009) 3083.

[17] V. Pureur, J.M. Dudley, Opt. Lett. 15 (2010) 2813.

[18] V. Pureur, A. Bétourné, G. Bouwmans, L. Bigot, A. Kudlinski, K. Delplace, A. Le Rouge, Y. Quiquempois, M. Douay, Fiber Integr. Opt. 28 (2009) 27.

[19] N.M. Litchinitser, S.C. Dunn, B. Usner, B.J. Eggleton, T.P. White, R.C. McPhedran, C.M. De Sterke, Opt. Express 11 (2003) 1243.

[20] T.A. Birks, F. Luan, G.J. Pearce, A. Wang, J.C. Knight, D.M. Bird, Opt. Express 14 (2006) 5688.

[21] A. Fuerbach, P. Steinvurzel, J.A. Bolger, A. Nulsen, B.J. Eggleton, Opt. Lett. 30 (2005) 830.

[22] T.P. White, B.T. Kuhlmey, R.C. McPhedran, D. Maystre, G. Renversez, C.M. De Sterke, L.C. Botten, JOSA B 19 (2002) 2322.

[23] V. Pureur, J.C. Knight, B.T. Kuhlmey, Opt. Express 18 (2010) 8906.

[24] B.T. Kuhlmey, B.J. Eggleton, D.K.C. Wu, J. Lightwave Technol. 27 (2009) 1617.

[25] J. Laegsgaard, Opt. Express 15 (2007) 16110.

[26] G.P. Agrawal, Nonlinear Fiber Optics, Academic Press, San Diego, 2007.

[27] P. Steinvurzel, C. Martijn de Sterke, M.J. Steel, B.T. Kuhlmey, B.J. Eggleton, Opt. Express 14 (2006) 8797.

[28] A. Bétourné, V. Pureur, G. Bouwmans, Y. Quiquempois, L. Bigot, M. Perrin, M. Douay, Opt. Express 15 (2007) 316.

[29] A. Bétourné, G. Bouwmans, Y. Quiquempois, M. Perrin, M. Douay, Opt. Lett. 32 (2007) 1719.