How long wavelengths can one extract from silica-core fibers?

Similar documents
International Conference on Information Sciences, Machinery, Materials and Energy (ICISMME 2015)

Highly Nonlinear Fibers and Their Applications

Generation of infrared supercontinuum radiation: spatial mode dispersion and higher-order mode propagation in ZBLAN step-index fibers

Nonlinear effects and pulse propagation in PCFs

Photonic crystal fiber with a hybrid honeycomb cladding

Generation of supercontinuum light in photonic crystal bers

Tailoring Nonlinearity and Dispersion of Photonic Crystal Fibers Using Hybrid Cladding

Estimating modal instability threshold for photonic crystal rod fiber amplifiers

Vectorial analysis of dielectric photonic crystal VCSEL

Observation of spectral enhancement in a soliton fiber laser with fiber Bragg grating

Empirical formulae for hollow-core antiresonant fibers: dispersion and effective mode area

Maximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers

Numerical modeling of the conduction and radiation heating in precision glass moulding

Increasing the blue-shift of a supercontinuum by modifying the fiber glass composition

Effective area of photonic crystal fibers

Research Article Nonlinear Phenomena of Ultra-Wide-Band Radiation in a Photonic Crystal Fibre

Sources supercontinuum visibles à base de fibres optiques microstructurées

Frequency-selective self-trapping and supercontinuum generation in arrays of coupled nonlinear waveguides

Fiber-Optics Group Highlights of Micronova Department of Electrical and Communications Engineering Helsinki University of Technology

Effect of cross-phase modulation on supercontinuum generated in microstructured fibers with sub-30 fs pulses

Berøringsløs optisk måling af gassammensætning ved UV og IR spektroskopi

Dmitriy Churin. Designing high power single frequency fiber lasers

The dipole moment of a wall-charged void in a bulk dielectric

Generalized Nonlinear Wave Equation in Frequency Domain

Mid-infrared supercontinuum covering the µm molecular fingerprint region using ultra-high NA chalcogenide step-index fibre

Ion exchange model for phase proton exchange waveguide in LiNbO3.

Molecular orientation via a dynamically induced pulse-train: Wave packet dynamics of NaI in a static electric field

Deep-blue supercontinnum sources with optimum taper profiles verification of GAM

Validity of PEC Approximation for On-Body Propagation

Demonstration of ultra-flattened dispersion in photonic crystal fibers

Advanced Vitreous State The Physical Properties of Glass

Experimental studies of the coherence of microstructure-fiber supercontinuum

Accelerating interior point methods with GPUs for smart grid systems

Title. Author(s)Nagasaki, Akira; Saitoh, Kunimasa; Koshiba, Masanori. CitationOptics Express, 19(4): Issue Date Doc URL.

Reflection of sound from finite-size plane and curved surfaces

Modelling Salmonella transfer during grinding of meat

Self-Phase Modulation in Optical Fiber Communications: Good or Bad?

Exposure Buildup Factors for Heavy Metal Oxide Glass: A Radiation Shield

On numerical evaluation of two-dimensional phase integrals

Reducing Uncertainty of Near-shore wind resource Estimates (RUNE) using wind lidars and mesoscale models

On the Einstein-Stern model of rotational heat capacities

Supercontinuum generation in dispersion engineered highly nonlinear (γ = 10 /W/m) As 2 S 3 chalcogenide planar waveguide

Optics Communications

Coherent Raman imaging with fibers: From sources to endoscopes

Optical solitons and its applications

Optical Cherenkov radiation in an As 2 S 3 slot waveguide with four zero-dispersion wavelengths

Threshold for strong thermal dephasing in periodically poled KTP in external cavity frequency doubling

Fiber Gratings p. 1 Basic Concepts p. 1 Bragg Diffraction p. 2 Photosensitivity p. 3 Fabrication Techniques p. 4 Single-Beam Internal Technique p.

An Equivalent Source Method for Modelling the Lithospheric Magnetic Field Using Satellite and Airborne Magnetic Data

A Note on Powers in Finite Fields

Harnessing and control of optical rogue waves in. supercontinuum generation

Copyright 2008 IEEE. Reprinted from IEEE Journal of Selected Topics in Quantum Electronics, 2007; 13 (3):

Sound absorption properties of a perforated plate and membrane ceiling element Nilsson, Anders C.; Rasmussen, Birgit

Blue-enhanced Supercontinuum Generation in a Fluorine-doped Graded-index Multimode Fiber

A short tutorial on optical rogue waves

Crack Tip Flipping: A New Phenomenon yet to be Resolved in Ductile Plate Tearing

Principle of photonic crystal fibers

Measurements of the angular distribution of diffuse irradiance

The colpitts oscillator family

Plane Wave Medical Ultrasound Imaging Using Adaptive Beamforming

Realtime 3D stress measurement in curing epoxy packaging

Feasibility of non-linear simulation for Field II using an angular spectrum approach

Nonlinear Optics (WiSe 2016/17) Lecture 9: December 16, 2016 Continue 9 Optical Parametric Amplifiers and Oscillators

Square Lattice Elliptical-Core Photonic Crystal Fiber Soliton-Effect Compressor at 1550nm

Single Shooting and ESDIRK Methods for adjoint-based optimization of an oil reservoir

University of Bath. Publication date: Document Version Peer reviewed version. Link to publication

Break up of the azimuthal symmetry of higher order fiber modes

Multi-cycle THz pulse generation in poled lithium niobate crystals

Reflection and absorption coefficients for use in room acoustic simulations

Design of Seven-core Photonic Crystal Fiber with Flat In-phase Mode for Yb: Fiber Laser Pumping

Diffusion of sulfuric acid in protective organic coatings

Breakup of Ring Beams Carrying Orbital Angular Momentum in Sodium Vapor

Four-wave mixing in PCF s and tapered fibers

Dynamic Effects of Diabatization in Distillation Columns

Validation of Boundary Layer Winds from WRF Mesoscale Forecasts over Denmark

The Glass Ceiling: Limits of Silica. PCF: Holey Silica Cladding

Alexander Gaeta Department of Applied Physics and Applied Mathematics Michal Lipson Department of Electrical Engineering

The Wien Bridge Oscillator Family

Supplementary Figure 1 Schematics of an optical pulse in a nonlinear medium. A Gaussian optical pulse propagates along z-axis in a nonlinear medium

Single-Polarization Single-Mode Photonic Band Gap Fiber

Matrix representation of a Neural Network

Group interactions of dissipative solitons in a laser cavity: the case of 2+1

Polarization Properties of Photonic Crystal Fibers Considering Thermal and External Stress Effects

Radioactivity in the Risø District July-December 2013

COMPRESSION of soliton pulses propagating in conventional

Design of Robust AMB Controllers for Rotors Subjected to Varying and Uncertain Seal Forces

2D PIM Simulation Based on COMSOL

Chalcogenide glass Photonic Crystal Fiber with flattened dispersion and high nonlinearity at telecommunication wavelength

Linear pulse propagation

Dispersion Properties of Photonic Crystal Fiber with Four cusped Hypocycloidal Air Holes in Cladding

Multiple Scenario Inversion of Reflection Seismic Prestack Data

Benchmarking of optimization methods for topology optimization problems

Detailed comparison of two approximate methods for the solution of the scalar wave equation for a rectangular optical waveguide

Lecture 4 Fiber Optical Communication Lecture 4, Slide 1

Determination of packaging induced 3D stress utilizing a piezocoefficient mapping device

Efficient surface second-harmonic generation in slot micro/nano-fibers

APPLICATION NOTE. Supercontinuum Generation in SCG-800 Photonic Crystal Fiber. Technology and Applications Center Newport Corporation

A practical implementation of microphone free-field comparison calibration according to the standard IEC

Three octave spanning supercontinuum by red-shifted dispersive wave in photonic crystal fibers

Transcription:

Downloaded from orbit.dtu.dk on: Dec 27, 28 How long wavelengths can one extract from silica-core fibers? Lægsgaard, Jesper; Tu, Haohua Published in: Optics Letters Link to article, DOI:.364/OL.38.458 Publication date: 23 Document Version Publisher's PDF, also known as Version of record Link back to DTU Orbit Citation (APA): Lægsgaard, J., & Tu, H. (23). How long wavelengths can one extract from silica-core fibers? Optics Letters, 38(2), 458-452. DOI:.364/OL.38.458 General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

458 OPTICS LETTERS / Vol. 38, No. 2 / November, 23 How long wavelengths can one extract from silica-core fibers? Jesper Lægsgaard, * and Haohua Tu 2 DTU Fotonik, Department of Photonics Engineering, Technical University of Denmark, Ørsteds Plads 343, DK-28 Kongens Lyngby, Denmark 2 Biophotonic Imaging Laboratory, University of Illinois at Urbana-Champaign, Urbana, Illinois 68, USA *Corresponding author: jlag@fotonik.dtu.dk Received August, 23; revised September 27, 23; accepted October, 23; posted October 2, 23 (Doc. ID 95); published October 3, 23 The generation of wavelengths above 3 μm by nonlinear processes in short silica photonic crystal fibers is investigated numerically. It was found that wavelengths in the 3 3.5 μm range may be generated quite efficiently in centimeter-long fiber pieces when pumping with femtosecond pulses in the.55 2 μm range. Wavelengths in the range of 3.5 4 μm can in principle be generated, but these require shorter fiber lengths for efficient extraction. The results indicate that useful 3 μm sources may be fabricated with existing silica-based fiber technology. 23 Optical Society of America OCIS codes: (6.239) Fiber optics, infrared; (6.437) Nonlinear optics, fibers; (6.5295) Photonic crystal fibers; (6.74) Ultrafast processes in fibers; (9.553) Pulse propagation and temporal solitons; (32.225) Femtosecond phenomena. http://dx.doi.org/.364/ol.38.458 Fiber-based light sources for the mid-infrared wavelength range above 3 μm are of high interest because fiberized systems are expected to have advantages in robustness and ease of operation compared to solutions based on wavelength conversion in nonlinear crystals. A number of soft-glass material platforms are currently being investigated for this purpose []. Pure silica fibers, such as photonic crystal fibers (PCFs), seem less attractive in this wavelength range because of the strongly increasing losses beyond 2 2.5 μm[2]. On the other hand, silica-based fiber technology has developed to a very mature state compared to most material alternatives, which allows for easy monolithic integration with silica-based fiber lasers [3]. It is therefore of interest to investigate how far this material platform can be pushed regarding mid-ir generation. Conventional fiber-based supercontinuum sources that utilize picosecond pump pulses and fiber lengths of several meters rapidly lose spectral density at wavelengths beyond 2.3 μm [4]. Similarly, four-wave mixing based approaches suffer low conversion efficiencies to long wavelengths because of loss limitations [5,6]. Xia et al. demonstrated an extension of the supercontinuum edge toward 3 μm using 5 cm fiber pieces pumped with nanosecond.55 μm pulses [7]. In our work, the generation of long wavelengths from femtosecond (fs) pump pulses at wavelengths of.55 and 2 μm in still shorter fibers was studied numerically. PCFs with tailored dispersion properties served as the nonlinear medium. It is shown that existing silica PCF technology allows for the construction of sources in the 3 4 μm range provided one is able to work with sufficiently short fiber pieces. The phenomenon of resonant energy transfer from solitons to dispersive waves in nonlinear fiber optics is well known [8], and the dispersive waves generated are often termed Cherenkov radiation (CR) because of formal analogies to Cherenkov radiation in one dimension. The intensity of the CR is controlled by the spectral density of the soliton at the CR frequency, ω CR, which in turn is given by a phase matching condition whose linear (intensity-independent) part reads [8] ~β ω CR ~β ω p ; ~β ω β ω β ω p ω ω : () Here, β is the frequency-dependent propagation constant of the fiber, β dβ dω and ω p is the center frequency of the pump pulse. It is important to note that CR emission is not critically dependent on the presence of a stable soliton [9]. For instance, when intense fs pulses are coupled into optical fibers with anomalous dispersion, one initially observes strong spectral broadening by self-phase modulation (SPM), which in turn leads to dispersive pulse compression. With sufficient power, this process can generate a few-cycle pulse that rapidly breaks up into one or more fundamental solitons and residual dispersive waves around ω p. However, around the point of maximal compression, one can observe strong transfer of optical power to CR frequencies fulfilling Eq. (), and occurring over a mm or sub-mm length scale. Such rapid spectral conversion relaxes the constraints imposed by fiber loss, which is the basic idea behind the mid-ir generation scheme proposed here. A standard PCF design, in which a triangular array of airholes make up the cladding region of the fiber, with one missing airhole defining the core, is capable of satisfying Eq. () over wide wavelength ranges with appropriate fiber parameters. This fact has been widely utilized for supercontinuum generation in the visible and near-ir wavelength regions []. Direct use of the CR mechanism for generating specific visible wavelengths from near-ir pump sources has also been demonstrated [3,]. For PCFs having two zero-dispersion wavelengths, CR on the long-wavelength side of the soliton may be observed [,2]. In Fig., ~β curves are shown for four fiber designs having d Λ.8, where d is the airhole diameter and Λ is the center-to-center distance between neighboring airholes. The curves were calculated by finding fully vectorial eigenmodes of Maxwell s equations with periodic 46-9592/3/2458-4$5./ 23 Optical Society of America

November, 23 / Vol. 38, No. 2 / OPTICS LETTERS 459 β β β (ω ω ) (µm - ).2. -. -.2 Λ=.5 µm, λ p =.55 µm Λ=.8 µm, λ p =.55 µm Λ=.75 µm, λ p =2 µm Λ=.95 µm, λ p =2 µm.5.5 2 2.5 3 3.5 4 Fig.. ~β β ω p curves for a PCF with d Λ.8, and pump wavelengths λ p characteristic of Er- or Tm-doped fiber lasers. boundary conditions in a planewave basis using the freely available software package MIT Photonic-Bands [3]. Silica material dispersion was described by a perturbative scheme [4] using the Sellmeier parameters given by Okamoto [5]. These are very similar to the parameters found by Malitson [6] by fitting to refractiveindex data ranging from 24 nm to 3.7 μm, thus providing some confidence in their use for the purposes of the present work. The pump wavelength was chosen to be either.55 or 2 μm, which is characteristic of either Er- or Tm-doped fiber lasers. The CR phase matching to wavelengths in the 3 4 μm range was readily obtained by scaling Λ in the range between.5 2 μm. There was also a phase matched ultraviolet (UV) CR frequency around 3 nm. Similar phase matching could also be obtained for lower d Λ values. The high d Λ value of.8 was chosen to ensure good confinement of the guided mode at the long CR wavelength. The numerical simulations of nonlinear propagation were done using the generalized nonlinear Schrödinger equation (GNLSE) formulated in [7], although only nontapered fibers were considered here. The frequency integrals involved in the evaluation of the nonlinear operator were done over the full frequency axis i.e., no restriction to positive frequencies was done. This makes the formalism fairly general, the main approximations being the restriction to propagation in one polarization state of the fundamental mode. Restricting the integrations to positive frequencies was found to have a significant impact on the calculated spectra, although the qualitative predictions were similar. Frequency grids spanning wavelengths from 2 25 nm to 6 μm were employed to fully capture the effects of long- and short-wavelength CR as well as non-phase-matched third-harmonic generation (THG). The main challenge for mid-ir generation is the silica material loss, so a reasonable model of the loss function is essential. In this work, the parametrization found by Frosz et al. was used [8,9] α A UV e λ UV λ A IR e λ IR λ ; (2) where A UV 3 db km, A IR 6 db km, λ UV 4.67 μm, and λ IR 47.8 μm. The IR constants were fitted to measurements in the.5 2.2 μm range. Loss curves for longer wavelengths have been reported by Izawa et al. [2]. The curve beyond 3 μm is bumpy and difficult to reproduce with a simple fit. Equation (2) compares reasonably with the experimental data, which include absorptions from a 3 ppm OH content, at 3 and 5 μm, while somewhat overestimating the losses in between. The predictions of the present paper should therefore be conservative. Note, however, that OH-induced absorptions will depend on the preparation and ambient conditions of the fiber. Also, confinement loss of the guided mode at the CR wavelength is not taken into account, since this can in principle be suppressed to any desired level by making the microstructured cladding sufficiently large. Figure 2 shows simulated output spectra at wavelengths of.55 and 2. μm for fibers with different Λ-values when pumping with Gaussian pulses having 5 fs full width at half-maximum (FWHM). The input pulse energy was 4 nj, which corresponds to a peak power of 75 kw. The propagation distance was set to cm. The long CR wavelength was easily tuned in the 3 4 μm range by changing Λ. The conversion efficiency decreased strongly with increasing CR wavelengths mainly because of the attenuation rising from.7 db cm at 3 μmto 4 db cm at 4 μm. As the main CR generation was found to occur within the first 2 4 mm of the fiber, the long-wavelength radiation was severely attenuated before the fiber exit. The spectra in Fig. 2 also show peaks from nonresonant THG, and UV CR in the case of.55 μm pumping. The UV CR peak also appeared for 2 μm pumping at higher energies. It is undesirable for a long-wavelength source, since the presence of intense UV light is likely to degrade the fiber over time. -6-6 Λ=.5 µm Λ=.6 µm Λ=.7 µm Λ=.8 µm 2 3 4 Λ=.75 µm Λ=.85 µm Λ=.95 µm.5.5 2 2.5 3 3.5 4 Fig. 2. Output spectra after cm of propagation for a 4 nj 5 fs input pulse at a wavelength of.55 μm and 2 μm. The PCF has d Λ.8.

452 OPTICS LETTERS / Vol. 38, No. 2 / November, 23 - -3 =5.3 nj =8 nj =.6 nj Filtered pulse energy (nj).75.5.25 λ>3.5 µm λ>4 µm 2 3 4 2 3 4 5 6 7 8 9 Propagation distance (mm) Filtered pulse energy (nj).... e-5 λ>3 µm λ>3.5 µm λ>4 µm λ>4 µm, Λ=.95 µm 5 5 Integrated pulse energy (nj) e-6 4 6 8 Input pulse energy (nj) Fig. 3. Output spectra after cm propagation for a PCF with Λ.75 μm, pumped by 5 fs pulses at 2 μm. Longpass filtered pulse energies versus input pulse energy for a fiber with Λ.75 μm and one with Λ.95 μm..5.5 2 2.5 3 3.5 4 4.5 Fig. 4. Long-pass filtered pulse energy versus propagation distance for a.6 nj, 5 fs, Gaussian input pulse at 2 μm wavelength in a fiber with Λ.85 μm. Spectrum and cumulative integral after 3.3 mm propagation. The conversion efficiency to the longer wavelengths can be increased by increasing the pump power. In Fig. 3, spectra for 2 μm pump pulses in a -cm PCF with Λ.75 μm are shown at different pulse energies with fluences well below expected damage thresholds [2]. Also shown are the integrated pulse energy above cutoff wavelengths of 3, 3.5, and 4 μm. The energy above 3 μm reached.5 nj for an input pulse energy of.6 nj. However, above 3.5 μm, the maximal extracted energy was around pj, and it appeared at an input energy of 8 nj. The decrease of the long-wavelength output energy with increasing pump power was due to the shift of the CR radiation toward longer wavelengths with correspondingly higher losses. Above 4 μm, the design with Λ.75 μm delivered very little energy. Increasing Λ to.95 μm gave an output of 3 pj for.6 nj input energy, which still represents a very poor conversion factor. The conversion efficiency to longer wavelengths may in principle be improved by further shortening of the fiber so that the desired long-wavelength radiation is generated just before the output facet. In this connection it should be noted that working with cm or sub-cm fiber pieces poses a substantial practical problem. In Fig. 4, the integrated energy beyond 3.5 and 4 μm for a.6 nj input pulse is shown as a function of propagation distance for a fiber with Λ.85 μm. The pulse energy above 3.5 μm peaks around 2.3 mm reaching a maximum of.6 nj i.e., a % conversion efficiency. The energy beyond 4 μm reaches a maximum of 42 pj after 3.3 mm. Given the assumption of ultrafast long-wavelength pump pulses and very short fiber lengths, these figures should be considered optimistic estimates of the conversion efficiencies obtainable in silica fibers. The spectrum after 3.3 mm is shown in the bottom panel of the figure along with the integrated (from the long-wavelength side) pulse energy. Note the appearance of the UV CR peak below 3 nm, which in this case carried around 5 pj of energy. The use of 5 fs input pulses improved conversion efficiency, but this is not a necessary requirement for generating long-wavelength CR. In Fig. 5, the spectrum of a 2.8 nj input pulse with 3 fs FWHM is shown after 2 cm propagation in a fiber with Λ.75 μm. The energy beyond 3 μm was 985 pj, which corresponded to a conversion efficiency of 7.7%. Pulse breakup and CR generation initiates after about cm propagation, a distance which decreases with increasing input power. Thus, the use of longer pulses relaxes the need for extremely short fiber pieces a bit. Although noise and shot-to-shot coherence was not modeled here, it must be expected that the use -.5 2 2.5 3 3.5 Fig. 5. Spectrum and its cumulative integral of a 3 fs 2.8 nj input pulse after 2 cm propagation in a fiber with Λ.75 μm. - Integrated pulse energy (nj)

November, 23 / Vol. 38, No. 2 / OPTICS LETTERS 452 of longer pulses will lead to a noisier CR signal in line with what one finds in fs-pumped continuum generation. In summary, the generation of mid-ir light in the 3 4 μm wavelength range was shown to be feasible with existing silica-based PCFs. Accessing the 3.5 4 μm range in an efficient manner was found to require very short (sub-cm) fiber pieces, which represents a major practical challenge. On the other hand, wavelengths around 3 μm seem readily accessible with fiber lengths of a few cm, considering material losses on the order of db cm. This wavelength range is of high interest for spectroscopy, imaging, and femtochemistry because of the presence of important fundamental stretching vibrations e.g., from O-H, N-H, and C-H chemical bonds. The main limitations with the proposed generation scheme are the need to work with very short fiber pieces, the small core dimensions needed to obtain CR phase matching, which make power scaling difficult, and the possible appearance of unwanted UV radiation in the fiber. None of these issues are of a fundamental nature. Phase matching in larger cores without simultaneous phase matching to UV wavelengths might be obtained by more advanced fiber designs. The problem of working with short fiber lengths could possibly be avoided by chirping the input pulse and/or using tapered fibers. Research along these lines is currently in progress. J. Lægsgaard acknowledges financial support by the Danish Council for Independent Research Technology and Production Sciences (FTP). H. Tu acknowledges grant support from the National Institutes of Health (R CA6639). References. J. H. V. Price, T. Monro, H. Ebendorff-Heidepriem, F. Poletti, P. Horak, V. Finazzi, J. Y. Y. Leong, P. Petropoulos, J. Flanagan, G. Brambilla, X. Feng, and D. Richardson, IEEE J. Sel. Top. Quantum Electron. 3, 738 (27). 2. T. Izawa, N. Shibata, and A. Takeda, Appl. Phys. Lett. 3,33 (977). 3. X. Liu, J. Lægsgaard, U. Møller, H. Tu, S. A. Boppart, and D. Turchinovich, Opt. Lett. 37, 2769 (22). 4. U. Møller, S. T. Sørensen, C. Jakobsen, J. Johansen, P. M. Moselund, C. L. Thomsen, and O. Bang, Opt. Express 2, 285 (22). 5. C. Jauregui, A. Steinmetz, J. Limpert, and A. Tünnermann, Opt. Express 2, 24957 (22). 6. A. Herzog, A. Shamir, and A. A. Ishaaya, Opt. Lett. 37, 82 (22). 7. C. Xia, M. Kumar, M.-Y. Cheng, O. P. Kulkarni, M. N. Islam, A. Galvanauskas, F. L. Terry, Jr., M. J. Freeman, D. A. Nolan, and W. A. Wood, IEEE J. Sel. Top. Quantum Electron. 3, 789 (27). 8. N. Akhmediev and M. Karlsson, Phys. Rev. A 5, 262 (995). 9. K. E. Webb, Y. Q. Xu, M. Erkintalo, and S. G. Murdoch, Opt. Lett. 38, 5 (23).. J. M. Dudley, G. Genty, and S. Coen, Rev. Mod. Phys. 78, 35 (26).. H. Tu and S. A. Boppart, Opt. Express 7, 9858 (29). 2. D. V. Skryabin, F. Luan, J. C. Knight, and P. S. J. Russell, Science 3, 75 (23). 3. S. G. Johnson and J. D. Joannopoulos, Opt. Express 8, 73 (2). 4. J. Lægsgaard, A. Bjarklev, and S. E. B. Libori, J. Opt. Soc. Am. B 2, 443 (23). 5. K. Okamoto, Fundamentals of Optical Waveguides (Academic, 2). 6. I. H. Malitson, J. Opt. Soc. Am. 55, 25 (965). 7. J. Lægsgaard, J. Opt. Soc. Am. B 29, 383 (22). 8. M. H. Frosz, P. M. Moselund, P. D. Rasmussen, C. L. Thomsen, and O. Bang, Opt. Express 6, 276 (28). 9. S. T. Sørensen, Deep-blue supercontinuum light sources based on tapered photonic crystal fibers, Ph.D. dissertation (Technical University of Denmark, 23). 2. R. Cherif, M. Zghal, I. Nikolov, and M. Danailov, Opt. Commun. 283, 4378 (2).