Journal of Electromagnetic Analysis and Applications
Vol.07 No.05(2015), Article ID:56130,10 pages

Numerical Study of Highly Nonlinear Photonic Crystal Fiber with Tunable Zero Dispersion Wavelengths

Alejandro Barrientos García1*, Igor A. Sukhoivanov1, José Amparo Andrade Lucio1, Oscar Gerardo Ibarra Manzano1, Igor Guryev2, Juan Carlos Hernández García1,3, Gabriel Ramos Ortiz4

1Department of Electronics Engineering, Universidad de Guanajuato, Salamanca, Mexico

2Department of Multidisciplinary Studies, Universidad de Guanajuato, Yuriria, Mexico

3CONACYT, Cátedras CONACYT, Ciudad de Mexico, Mexico

4Centro de Investigaciones en Óptica, A. C., GPOM, León, Mexico

Email: *

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 31 January 2015; accepted 5 May 2015; published 6 May 2015


Solid-core silica photonic crystal fiber is proposed borrowing the concept of golden ratio (1.618) and keeping it between pitch and air hole diameter in a subset of six rings of air-holes with hexagonal arrangement. In the case when we have a pitch equal to one micron (Λ = 1 μm), we need air-holes diameters d = 0.618 μm in order to achieve two zero dispersion wavelength (ZDW) points at 725 nm and 1055 nm; this gives us the possibility to use the fiber for supercontinuum generation, by pumping close to that points, pulse compression or reshaping. We analyzed a series of fibers using this relation and showed the possibilities of tunable ZDW in a wide range of wavelengths from 725 nm to 2000 nm, with low losses and small effective area. In agreement with the ZDW point needed, the geometry of the structure can be modified to the point of having only three rings of air holes that surround the solid core with low losses and good confinement mode. The design proposed here is analyzed using the finite element method with perfectly matched layers, including the material dispersion directly into the model applying the Sellmeier’s equation.


Photonic Crystal Fibers, Supercontinuum Generation, Zero Dispersion Wavelength

1. Introduction

Today the necessities for different applications in science and engineering of photonic crystal fibers (PCF) have boosted proposals on novel designs that satisfy specific conditions of linear and non-linear properties. Because of this, the control of dispersion, losses, mode confinement and conditions for generation of supercontinuum or pulse reshaping is a wide field of research. A series of techniques implying modification of some parameters of the fiber for various applications have been published. Since the invention of PCFs [1] [2] , their designs and applications have been increased. However, some structures suffer from complex geometry. Traditionally optical fibers are used in various applications, for instance, to transport laser energy for material machining, as a gain medium in fiber lasers or optical telecommunication [3] . However, new possibilities of applications arise when nonlinear properties of PCFs are manifested. Specifically, nonlinearities in PCFs allow new light sources with very broad spectra, i.e., supercontinuum generation (SCG) [4] [5] . The possibility to control group velocity dispersion (GVD) in PCFs, their small effective area and elevated nonlinearity [6] has allowed the observation of SCG in a much wider range of source parameters than has been possible with bulk media or conventional fibers [7] . It is known that in the case of normal GVD is present in the PCFs, the Raman scattering effect and self-phase modulation (SPM) dominate the spectrum broadening [8] , and soliton dynamics govern spectral broadening when central wavelength of the propagating pulse is in anomalous GVD. Finite element and finite-difference methods as well as the multipole method and the plane wave expansion method [9] [10] were used to characterize the electromagnetic field distribution in PCFs. The accurate control of the dispersion characteristics is obligatory in variety of applications, such as pulse recompression [11] mode locked fiber lasers [12] , nonlinear effects, including Raman scattering, four-wave mixing, SPM or soliton generation [13] , and SCG [5] [7] . The flexible geometry of PCFs offers many unique and novel dispersion properties in comparison to conventional single mode fibers such as nearly zero ultraflattened dispersion [14] -[16] , high negative chromatic dispersion [17] , and anomalous dispersion in short-wavelength regime [18] . As is known, the chromatic dispersion coefficient of an optical fiber is a key quantity in various analysis and design issues in fiber transmission systems [19] [20] .

According to the guidance mechanism, PCFs can be divided into two categories: first, total internal reflection (TIR), where air holes surround a high index solid core, and second, photonic band gap (PBG), where the light is confined to a central air hole [21] . In PCFs, the cladding region is usually formed by array of air holes in a regular triangular lattice with separation Λ [22] . Although there is a rich nomenclature of PCFs at this moment, we will restrict ourselves to index-guiding PCFs. In this kind of fibers, the core index is greater than the average index of the cladding because of the presence of air holes, and the fiber can guide the light by total internal reflection as a standard fiber does. That is, the guided mode has effective index that satisfies the condition


where is the mode propagation constant, is the core index, and is the cladding effective index, which is termed also the effective index of the fundamental space-filling mode (FSM). In case of a PCF made of pure silica, is the refractive index of silica. The losses in PCFs occur due to intrinsic material absorption, structural imperfection, Rayleigh scattering, and imperfect confinement. Fabrication-related losses can be reduced by carefully optimizing the fabrication process [23] . Confinement losses are an additional form of energy leakage that occurs in single-material fibers because of inherently leaky nature of the guided modes in these fibers [24] . This is because the core index is the same as the index of the outer cladding without air holes. Confinement losses can be reduced by increasing the number of rings of air holes which surround the solid core, and thus are determined by the geometry of the structure. It is important to know how many rings of air holes are required to reduce the confinement loss under the Rayleigh scattering limit. At the same time, geometry of the structure allows the chromatic dispersion to be easily controlled by varying the air hole diameter and the pitch of PCF. Here we show that some properties such as losses, ZDW and intermodal dispersion can be controlled keeping the ratio between fiber pitch and hole diameter of the microstructured cladding equal to the called number golden ratio, i.e., 1.618.

Controllability of chromatic dispersion in PCFs is a very important problem for practical applications in optical communication systems, and nonlinear optics. The chromatic dispersion profile can be easily controlled by varying the values of any or all of the following parameters: air hole diameter, pitch, number of rings (Nr) and the number of missing holes (MNr). The aim of this work is to demonstrate how the dispersion profile may be modified by variation of the pitch and the hole diameter, keeping the golden ratio between them, i.e., keeping a ratio between them equal to 1.618, providing a tunable ZWD in the range of 725 nm to 2000 nm.

2. Characteristics of the Fiber

Figure 1 shows the geometry of the solid core PCF proposed for SCG. The PCF consists of an arrangement of six rings with hexagonal structure and pitch size from 1 to 2 μm (Λ) and air hole diameter from 0.618 to 1.236 μm. This PCF has small effective area and tunable zero wavelength dispersion points in the range from 725 to 2000 nm. Although PCFs can be characterized by variety of methods including finite-difference time-domain method [25] we have used finite-element (FEM) approximation of the Helmholtz equation applying Sellmeier’s equation to take the material dispersion into account when calculating effective mode index. The curve of total chromatic dispersion of the PCFs is obtained using following expression


where is the speed of light in vacuum, is the effective mode index, and is the radiation wave- length.

This structure, proposed to tune the ZDW coinciding with the pump laser wavelengths to facilitate SCG or pulse reshaping applications, is inspired by the nature where the golden ratio can be found in perfect structures of the leaves, animals or the human body. Thus, we propose to use such a golden ratio in the design of a PCF and evaluate its impact on the optical properties. When this proportion is satisfied in the nature, the efficiency of the structures in macro- and micro-world is increased a lot, independently of the eye appeal. Talking about micro-world, we can found the golden ratio in the structure of some viruses, living cells, deoxyribonucleic acid, and bodies of honey bee, leaves or animals. In the macro-world, golden ratio can be found in the structure of galaxies. Due to this we probed to use the golden ratio for synthesis of photonic microstructures just like it happened in architecture to create some design of buildings, for example, keeping this relation. We found that this relation can improve the properties of the PCF; when we keep pitch-to-hole-diameter ratio equal to the golden ratio, it helps us to obtain simpler structures with better characteristics corresponding to the application of need.

In Figure 2 and Figure 3 we show the effective area and nonlinear coefficient of the designed PCF with the pitch Λ = 1, 1.5, 2 μm and air-holes diameter d = 0.618, 0.927, 1.236 μm. Using only six rings of air holes surrounding the solid core we have good confinement of the fundamental mode, small effective area, and high nonlinear coefficient. The structure keeps the golden ratio between pitch and diameter, and with this, we achieve

Figure 1. Cross-section of the PCF proposed where we keep the golden ratio between pitch and air hole diameter (Λ/d).

Figure 2. Effective area with Λ = 1, 1.5, 2 μm and d = 0.618, 0.927, 1.236 μm.

Figure 3. Nonlinear coefficient with Λ = 1, 1.5, 2 μm and d = 0.618, 0.927, 1.236 μm.

an improvement of the characteristics of PCF correspondingly to an application needed. The effective area and the nonlinear coefficient have been calculated applying the concepts of conventional fiber theory [26] with modifications necessary to include photonic-crystal cladding into account [26] -[28] .





where is the effective core radius, is the operating wavelength, is the core index, is the effective index of the fundamental space-filling mode, is the half of mode field diameter (MFD) called also the effective modal spot size, is the effective mode area, is the correction factor [29] , is the nonlinear coefficient [30] and is the nonlinear index of silica. In Section 4 we show an application of the supercontinuum generation using this specific value of nonlinear coefficient shown in Figure 3.

3. Tunable Zero Dispersion Wavelength

In Figure 4 we show the calculated total chromatic dispersion obtained using submicron air hole diameters in order to keep the golden ratio parameter. We can see that varying the pitch from 1 to 1.5 μm two different ZDW points are obtained which are tunable depending on the value of the pitch. Moreover, the curve becomes more gently sloped with increasing the value of pitch, giving the possibility to use this kind of fiber in different applications. In Section 4 we use this design in simulation of super continuum generation. The curves show a maximal chromatic dispersion from 38 to 80 ps/(nm*km) and always we obtain two ZWD points which can be used for pumping by two different lasers to generate a more wide and flattened SC spectrum [31] -[33] .

In Figure 5 we show chromatic dispersion curves obtained when air hole diameters exceed 1 μm. Therefore, the pitch shall be from 1.618 to 2 μm in order to maintain the golden ratio proportion. The curves shows a maximal chromatic dispersion from 85 to 90 ps/(nm*km), and, as in the previous case, we always have two ZWD points in the range of 725 - 2000 nm. We believe that fabrication of this kind of fiber is simple due to micrometric scale of the constituents.

4. Simulation of Super Continuum Generation, Pulse Reshaping and Compression

The PCFs proposed here were analyzed using commercial software COMSOL Multiphysics which uses the proven finite element method (FEM). This software runs the finite element analysis together with adaptive meshing and error control using a variety of numerical solvers. The RF Module of COMSOL provides the capabilities to design and is used to simulate electromagnetic wave propagation in PCFs. The mode analysis is made on a cross-section in the xy-plane of the fiber, by expressing the propagating wave in z direction in the following form

Figure 4. A family of the chromatic dispersion curves calculated by keeping the golden ratio. The ZDW points are moved from 725 to 1800 nm when pitch and air hole diameter are varied.

Figure 5. The behavior of chromatic dispersion with variation of geometry parameters, when diameter of holes is greater than 1 μm.


and, thereupon, solving the homogeneous Helmholtz equation


for the eigenvalue applying circular perfectly matched layer, with being the propagation constant.

The pulse propagation in the PCFs is numerically investigated within the Generalized Nonlinear Schrödinger Equation (GNLSE) [34] [35]


where is the attenuation constant, is the complex field envelope, is the distance, is the retarded time traveling at the envelope group velocity, is the nonlinear coefficient and are the dispersion coefficients obtained by a Taylor series expansion of the propagation constant around the center frequency. The GNLSE is solved by Fourth-Order Runge-Kutta in the Interaction Picture Method (RK4IP) [36] [37] .

The Raman response function for silica fibers is described by [38]


To simulate the super continuum generation we consider the waveform with the full width at half maximum (FWHM), TFWHM = 50 fs, and 10 kW of peak power and 873 nm, propagating into a 5-mm length of a PCF with maximal chromatic dispersion at 873 nm of 39.78 ps/(nm*km). The PCF designed has low losses, small effective area (Aeff = 1.76 μm2 at 873 nm), and large nonlinear coefficient in the range of 800 to 1600 nm. Pulse evolution has been simulated using RK4IP method for solution of the generalized nonlinear Schrödinger equation The fiber parameters used in the calculation are given as, , , , , in which we achieve a more wide spectrum that reported using similar parameters.

Figure 6 shows the super continuum spectrum obtained by simulation of ultra short pulse propagation in the designed PCF, where we use the property of anomalous dispersion with high nonlinear coefficient of the PCF proposed, for creation a wide spectrum using femtosecond pump laser and shorter length of PCF. We can see a relatively smooth spectrum and wide SC spectrum of 2100 nm bandwidth at the level of −40-dB. In Figure 7 we show temporal evolution of the pump pulse. We observe that soliton-like pump pulse has broken onto seven waveforms which can be interpreted as pulses; it gives the output pulse a form of the train of very short pulses. Taken this observation into account together with the fact that wavelength of the pump pulse is in anomalous dispersion region of the PCF designed, we conclude that Figure 7 show us manifestation of the soliton fission, which is known as a primary nonlinear effect under these propagation conditions [7] . In the spectral domain the

Figure 6. Output spectrum obtained at fiber length of 5 mm. The curve show SC spectrum with the 40-dB bandwidth of 2100 nm using a pump laser centered at 873 nm.

Figure 7. The red curve shows the initial pulse (sech2 waveform) and the blue curve shows the waveform after propagation in 5 mm-long PCF.

soliton fission manifests itself as an extreme broadening of the spectrum of incident pulse. It is exactly what we observe in Figure 6. Thus, the main effect of spectral broadening in our case is the soliton fission.

The structure of the PCF designed has also been used for investigation of a possibility of pulse compression and transformation of the pulse shape supposing a sech2 initial pulse envelope (IPE) waveform. Each of these applications requires a particular combination of dispersion and nonlinearity. Transformation of pulse to another high-quality waveform in a fiber requires usually a normal dispersion at the central wavelength of the incident pulse [8] [39] , as the primary phenomenon in this case is an interplay between the self-phase modulation and normal dispersion [8] [40] -[42] . Therefore, for simulation of pulse shape transformation we use following parameters and central wavelength of the pulse to be 600 nm: TFWHM = 140 fs, peak power P = 8 kW, fiber length 13.5 mm, , , , , . The wavelength of 600 nm has been chosen because of all dispersion curves we obtained using golden ratio, demonstrate normal dispersion at this wavelength; and this is exactly what we need to achieve pulse transformation with sufficiently high quality of the output pulse.

All of golden-ratio combinations of pitch and hole diameter demonstrate highly nonlinear dispersion profile; therefore, those fibers cannot be applied directly to pulse compression like it made in case of chirped mirrors, for example [43] [44] . However, these combinations provide wavelength ranges with anomalous group velocity dispersion where formation of higher-order solitons is possible. Thus, within these ranges one can compress the pulse using soliton-effect compression [45] . This has been implemented for the case of Λ = 1.618 μm and air-holes diameter with five rings of holes in the cladding. In this case, zero dispersion wavelengths are 800 nm and 1956 nm, and maximum of the chromatic dispersion equals to 99.27 ps/(nm*km) @1303 nm, as shown in Figure 4. Other parameters are following: TFWHM = 250 fs @1303 nm, peak power P = 10 kW, , , , ,. The optimal length of the fiber when maximum compression is achieved, equals to 15 mm. The compression factor in this case equals to 100. Figure 8 shows the incident pulses for cases of shape transformation and compression, compressed pulse, and reshaped pulse with a flat-top waveform. The shape of the compressed pulse demonstrates a pedestal, which is typical for pulses obtained via soliton-effect compression. Spectra of corresponding output pulses are shown in Figure 9. Spectrum of the flat-top pulse demonstrates signatures of the self-phase modulation. This nonlinear effect together with normal group velocity dispersion at 600 nm produces flat top pulses. A closer look to the reshaped pulse in Figure 8 allows noting oscillations at the trailing edge of the pulse, which is manifestation of the optical wave breaking [46] . These results indicate a possibility to use the PCF proposed in variety of engineering problems of nonlinear optics.

Figure 8. Initial pulse envelope (green and red curve), compressed pulse envelope (blue curve) and flat-top pulse (black curve).

Figure 9. Spectrum for a compressed pulse and flat-top pulse. The dash line show the spectrum of compressed pulse and solid line show the spectrum of the flat-top pulse.

5. Conclusion

In this work, we proposed to evaluate how the optical properties of a PCF were influenced when the proportion between pitch and diameter of holes was equal to the golden ratio. In particular, when this ratio is maintained and the pitch and diameter holes are varied, the ZDW could be tunable between the range of 725 nm to 2000 nm. Additionally, with a low number of rings, we obtain a small loss confinement with small effective area and high nonlinear coefficient. Finally, we obtain a flattened chromatic dispersion curve between of two ZDW points. We present an example of applications of this kind of designs for supercontinuum generation, pulse recompression and pulse reshaping. A good quality of flattened and wide spectrum is achieved. The spectrum possesses 2100 nm bandwidth at 40 dB level. The compression factor achieved is 100 producing a compressed pulse of 2.5 fs from an initial pulse of 250 fs. It is believed that the design proposed could be fabricated easily with techniques available nowadays.


This work was supported in part by the National Council for the Science and technology (CONACYT) of Mexico, and by Universidad de Guanajuato (project DAIP-192/13 and DAIP-430/2014), FOMIX (CONACyT & CONCyTEG) under grant GTO-2012-C04-195229 and DAIP 458/2014. A. Barrientos-García and J. Carlos Hernandez-Garcia are grateful to CONACYT for a student research grant No. 329409/336847 and CONACYT project no. 3155.


  1. Knight, J.C. (2003) Photonic Crystal Fibres. Nature, 424, 847-851.
  2. Russell, P. (2003) Photonic Crystal Fibers. Science, 299, 358-362.
  3. Böswetter, P., Baselt, T., Ebert, F., Basan, F. and Hartmann, P. (2011) Group-Velocity Dispersion in Multimode Photonic Crystal Fibers Measured Using Time-Domain White-Light Interferometry. Proceedings of SPIE, 7914, Article ID: 791422.
  4. Liao, M., Chaudhari, C., Qin, G., Yan, X., Suzuki, T. and Ohishi, Y. (2009) Tellurite Microstructure Fibers with Small Hexagonal Core for Supercontinuum Generation. Optics Express, 17, 12174-12182.
  5. Gaeta, A.L. (2002) Nonlinear Propagation and Continuum Generation in Microstructured Optical Fibers. Optics Letters, 27, 924-926.
  6. Reeves, W.H., Skryabin, D.V., Biancalana, F., Knight, J.C., Russell, P.St.J., Omenetto, F.G., Efimov, A. and Taylor, A.J. (2003) Transformation and Control of Ultra-Short Pulses in Dispersion-Engineered Photonic Crystal Fibres. Nature, 424, 511-515.
  7. Dudley, J.M., Genty, G. and Coen, S. (2006) Supercontinuum Generation in Photonic Crystal Fibers. Reviews of Modern Physics, 78, 1135-1184.
  8. Sukhoivanov, I.A., Iakushev, S.O., Shulika, O.V., Díez, A. and Andrés, M. (2013) Femtosecond Parabolic Pulse Shaping in Normally Dispersive Optical Fibers. Optical Express, 21, 17770-17785.
  9. Johnson, S.G. and Joannopoulos, J.D. (2001) Block-Iterative Frequency-Domain Methods for Maxwell’s Equations in a Planewave Basis. Optics Express, 8, 173-190.
  10. White, T.P., Kuhlmey, B.T., McPhedran, R.C., Maystre, D., Renversez, G., de Sterke, C.M. and Botten, L.C. (2002) Multipole Method for Microstructured Optical Fibers. I. Formulation. Journal of the Optical Society of America B, 19, 2322-2330.
  11. Laegsgaard, J. and Roberts, P.J. (2008) Dispersive Pulse Compression in Hollow-Core Photonic Band-Gap Fibers. Optics Express, 16, 9628-9644.
  12. Lim, H. and Wise, F. (2004) Control of Dispersion in a Femtosecond Ytterbium Laser by Use of Hollow-Core Photonic Band-Gap Fiber. Optics Express, 12, 2231-2235.
  13. Dudley, J.M. and Taylor, J.R. (2010) Supercontinuum Generation in Optical Fibres. Cambridge University Press, Cambridge, 32-96.
  14. Ferrando, A., Silvestre, E., Miret, J.J., Monsoriu, J.A., Andres, M.V. and Russell, P.St.J. (1999) Designing a Photonic Crystal Fiber with Flattened Chromatic Dispersion. Electronics Letters, 35, 325-326.
  15. Ferrando, A., Silvestre, E., Miret, J.J. and Andres, P. (2000) Nearly Zero Ultraflattened Dispersion in Photonic Crystal Fibers. Optics Letters, 25, 790-792.
  16. Reeves, W.H., Knight, J.C., Russell, P.St.J. and Roberts, P.J. (2002) Demonstration of Ultra-Flattened Dispersion in Photonic Crystal Fibers. Optics Express, 10, 609-613.
  17. Birks, T.A., Mogilevstev, D., Knight, J.C. and Russell, P.St.J. (1999) Dispersion Compression Using Single Material Fibers. IEEE Photonics Technology Letters, 11, 674-676.
  18. Knight, J.C., Arriaga, J., Birks, T.A., Ortigosa-Blanch, A., Wadsworth, W.J. and Russell, P.St.J. (2000) Anomalous Dispersion in Photonic Crystal Fiber. IEEE Photonics Technology Letters, 12, 807-809.
  19. Mogilevtsev, D., Birks, T.A. and Russell, P.S.J. (1998) Group-Velocity Dispersion in Photonic Crystal Fibers. Optics Letters, 23, 1662-1664.
  20. Knight, J.C., Arriaga, J., Birks, T.A., Ortigosa-Blanch, A., Wadsworth, W.J. and Russell, P.S.J. (1998) Anomalous Dispersion in Photonic Crystal Fiber. IEEE Photonics Technology Letters, 12, 817-819.
  21. Cheng, H., Crutchfield, W.Y., Doerg, M. and Greengard, L. (2004) Fast, Accurate Integral Equation Methods for the Analysis of Photonic Crystal Fibers I: Theory. Optics Express, 12, 3791-3805.
  22. Saitoh, K., Koshiba, M., Hasegawa, T. and Sasaoka, E. (2003) Chromatic Dispersion Control in Photonic Crystal Fibers: Application to Ultra Flattened Dispersion. Optics Express, 11, 843-852.
  23. Tajima, K., Zhou, J., Kurokawa, K. and Nakajima, K. (2011) Low Water Peak Photonic Crystal Fibers. European Conference on Optical Communication (ECOC), Rimini, September 2003, Th4.1.6.
  24. White, T.P., McPhedran, R.C., de Sterk, C.M., Botten, L.C. and Steel, M.J. (2001) Confinement Losses in Microstructured Optical Fibers. Optics Letters, 26, 1660-1662.
  25. Shulika, O., Guryev, I., Gurieva, N., Sukhoivanov, I. and Andrade Lucio, J.A. (2013) Graphics Processing Unit-Acce- lerated Finite-Difference Time-Domain Method for Characterization of Photonic Crystal Fibers. Optical Engineering, 52, Article ID: 126109.
  26. Marcuse, D. (1977) Loss Analysis of Single-Mode Fiber Splices. Bell System Technical Journal, 56, 703-718.
  27. Koshiba, M. and Saitoh, K. (2004) Applicability of Classical Optical Fiber Theories to Holey Fibers. Optics Letters, 29, 1739-1741.
  28. Mortensen, N.A., Folkenberg, J.R., Nielsen, M.D. and Hansen, K.P. (2003) Modal Cut-Off and the V Parameter in Photonic Crystal Fibers. Optics Letters, 28, 1879-1881.
  29. Miyagi, K., Namihira, Y., Razzak, S.M.A., Kaijage, S.F. and Begum, F. (2010) Measurements of Mode Field Diameter and Effective Area of Photonic Crystal Fibers by Far-Field Scanning Technique. Optical Review, 17, 388-392.
  30. Liu, Z.L., Hou, L.T. and Wang, W. (2009) Tailoring Nonlinearity and Dispersion of Photonic Crystal Fibers Using Hybrid Cladding. Brazilian Journal of Physics, 39, 50-54.
  31. Frosz, M., Falk, P. and Bang, O. (2005) The Role of the Second Zero-Dispersion Wavelength in Generation of Supercontinua and Bright-Bright Soliton-Pairs across the Zero-Dispersion Wavelength. Optics Express, 13, 6181-6192.
  32. Hilligsøe, K.M., Andersen, T., Paulsen, H., Nielsen, C., Mølmer, K., Keiding, S., Kristiansen, R., Hansen, K. and Larsen, J. (2004) Supercontinuum Generation in a Photonic Crystal Fiber with Two Zero Dispersion Wavelengths. Optics Express, 12, 1045-1054.
  33. Tse, M.L.V., Horak, P., Poletti, F., Broderick, N.G., Price, J.H., Hayes, J.R. and Richardson, D.J. (2006) Supercontinuum Generation at 1.06µm in Holey Fibers with Dispersion Flattened Profiles. Optics Express, 14, 4445-4451.
  34. Shulika, O.V., Sukhoivanov, I.A., Guryev, I.V., Andrade-Lucio, J.A., Ibarra-Manzano, O., Barrientos-Garcia, A. and Ramos-Ortiz, G. (2013) Numerical Study of Few-Cycle Pulses Generation from Supercontinuum in ANDi-PCF. Acta Universitaria, 23, 40-44.
  35. Shulika, O.V., Sukhoivanov, I.A., Ramos-Ortiz, G., Guryev, I.V., Iakushev, S.O. and Andrade Lucio, J.A. (2013) Characterization of All-Normal Dispersion Microstructured Optical Fiber via Numerical Simulation of Passive Nonlinear Pulse Reshaping and Single-Pulse Flat-Top Supercontinuum. Journal of Nanophotonics, 8, Article ID: 083890.
  36. Hult, J. (2007) A Fourth-Order Runge-Kutta in the Interaction Picture Method for Simulating Supercontinuum Generation in Optical Fibers. Journal of Lightwave Technology, 25, 3770-3775.
  37. Hernandez-Garcia, J.C., Pottiez, O., Estudillo-Ayala, J.M. and Rojas-Laguna, R. (2012) Numerical Analysis of a Broadband Spectrum Generated in a Standard Fiber by Noise-Like Pulses from a Passively Mode-Locked Fiber Laser. Optics Communications, 285, 1915-1919.
  38. Blow, K.J. and Wood, D. (1989) Theoretical Description of Transient Stimulated Raman Scattering in Optical Fibers. IEEE Journal of Quantum Electronics, 25, 2665-2673.
  39. Iakushev, S.O., Shulika, O.V. and Sukhoivanov, I.A. (2012) Passive Nonlinear Reshaping towards Parabolic Pulses in the Steady-State Regime in Optical Fibers. Optics Communications, 285, 4493-4499.
  40. Agrawal, G. (2012) Nonlinear Fiber Optics. 5th Edition, Academic Press, Waltham, 648 p.
  41. Iakushev, S.O., Shulika, O.V., Sukhoivanov, I.A., Fesenko, V.I., Andrés, M.V. and Sayinc, H. (2014) Formation of Ultrashort Triangular Pulses in Optical Fibers. Optics Express, 22, 29119-29134.
  42. Sukhoivanov, I.A., Iakushev, S.O., Shulika, O.V., Andrade Lucio, J.A., Díez, A. and Andrés, M. (2014) Supercontinuum Generation at 800 nm in All-Normal Dispersion Photonic Crystal Fiber. Optics Express, 22, 30234-30250.
  43. Yakushev, S.O., Shulika, O.V., Petrov, S.I. and Sukhoivanov, I.A. (2008) Chirp Compression with Single Chirped Mirrors and Its Assembly. Microelectronics Journal, 39, 690-695.
  44. Iakushev, S.O., Shulika, O.V., Lysak, V.V. and Sukhoivanov, I.A. (2008) Air-Gap Silicon Nitride Chirped Mirror for Few-Cycle Pulse Compression. Optoelectronics and Advanced Materials―Rapid Communications, 2, 686-688.
  45. Agrawal, G. (2008) Application of Nonlinear Fiber Optics. 2nd Edition, Academic Press, Waltham.
  46. Iakushev, S.O., Shulika, O.V. and Garcia, J.J.R. (2010) Optical Wave Breaking Cancelation and Formation of Quasi- Parabolic Ultrashort Pulses. 10th International Conference on Laser and Fiber-Optical Networks Modeling, Sevastopol, 12-14 September 2010, 176-179.


*Corresponding author.