Spatial-temporal structure of the femtosecond third harmonic generation in photonic-crystal fibers

We analyze the third harmonic generation by the supercontinuum at 1.5μm in nonlinear-glass microstructured fibers. The numerical model includes the real field dynamics and allows analyzing t he spectral as well as the temporal structure of the generated field. The cal culated third harmonic and supercontinuum spectra are compared with the e xperiment in the SF6 glass PCF pumped by a femtosecond Cr 4+:YAG oscillator. The spectral structure of the third harmonic is composed of t he 10–20 high-order modes excited at the phase matching points aroun d 500–550 nm. The individual third harmonic modes have spectral width s of 2–4 nm. In the time domain, the third harmonic signal is a ∼10-ps pulse with psscale slow modulation, containing subpicosecond ripples a t it trailing edge. © 2007 Optical Society of America OCIS codes:(190.5530) Pulse propagation and solitons; (190.4380) Non linear optics, fourwave mixing. References and links 1. J. M. Dudley, G. Genty, S. Coen, “Supercontinuum generati on in photonic crystal fiber,” Rev. Mod. Phys. 78, 1135–1184 (2006). 2. V. L. Kalashnikov, E. Sorokin, I. T. Sorokina, “Raman effec ts in the infrared supercontinuum generation in softglass PCFs,” Appl. Phys. B. 87, 37–44 (2007) 3. R. H. Stolen, “Phase-matched-stimulated four-wave mixing i n sil ca-fiber waveguides,” IEEE J. Quantum Electron.QE-11,100–103 (1975). 4. R. H. Stolen and W. H. Leibolt, “Optical fiber modes using sti mulated four photon mixing,” Appl. Opt. 11, 239–243 (1976). 5. J. M. Dudley, L. Provino, N. Grossard, H. Mailotte, R. S. Wi ndeler, B. J. Eggleton, S. Coen, “Supercontinuum generation in air-silica microstructured fibers with nanose c nd and femtosecond pulse pumping,” J. Opt. Soc. Am. B 19,765–771 (2002). 6. J. H. V. Price, T. M. Monro, K. Furusawa, W. Belardi, J. C. Ba ggett, S. Coyle, C. Netti, J. J. Baumberg, R. Paschotta, D. J. Richardson, “UV generation in a pure-silic a holey fiber,” Appl. Phys. B77,291–298 (2003). 7. C. G. Durfee III, S. Backus, M. M. Murnane, and H. C. Kapteyn , “Ultrabroadband phase-matched optical parametric generation in the ultraviolet by use of guided waves,” Opt. Lett.22,1565–1567 (1997). 8. A. M. Zheltikov, N. I. Koroteev, A. N. Naumov, “Selfand cro ss-phase modulation accompanyning thirdharmonic generation in a hollow waveguide,” J. Exp. Theor. Ph ys.88,857–867 (1999). 9. J. K. Ranka, R. S. Windeler, A. J. Stentz, “Optical propert ies of high-delta air-silica microstructure optical fibers, ” Opt. Lett.25,796–798 (2000). 10. A. N. Naumov, D. A. Sidorov-Biryukov, F. Giammanco, A. B. Fed otov, P. Marsili, A. Ruffini, O. A. Kolevatova, and A. M. Zheltikov, “Four-wave mixing of picosecond pulses i n hollow fibers: phase matching and the influence of high-order waveguide modes,” J. Exp. Theor. Phys. 93,247–255 (2001). #83343 $15.00 USD Received 24 May 2007; revised 19 Jul 2007; accepted 15 Aug 2007; published 22 Aug 2007 (C) 2007 OSA 3 September 2007 / Vol. 15, No. 18 / OPTICS EXPRESS 11301 11. S. O. Konorov, A. B. Fedotov, and A. M. Zheltikov, “Enhanc ed four-wave mixing in a hollow-core photonic crystal fiber,” Opt. Lett.28,1448–1450 (2003). 12. K. D. Moll, D. Homoelle, and A. L. Gaeta, “Conical harmonic g eneration in isotropic matherials,” Phys. Rev. Lett. 88,153901 (2002). 13. T. Schreiber, J. Limpert, H. Zellmer, A. T ünnermann, K. P. Hansen, “High average power supercontinuum generation in photonic crystal fibers,” Opt. Commun. 228,71–78 (2003). 14. A. N. Naumov, A. B. Fedotov, A. M. Zheltikov, V. V. Yakovlev , L. A. Mel’nikov, V. I. Beloglazov, N. B. Skibina, A. V. Shcherbakov, “Enhanced χ(3) interactions of unamplified femtosecond Cr:forsterite laser pulses in photonic-crystal fibers,” J. Opt. Soc. Am. B 19,2183–2190 (2002). 15. F. G. Omenetto, A. J. Taylor, M. D. Moores, J. Arriaga, J. C. Knight, W. J. Wadsworth, P. S. Jt. Russell, “Simultaneous generation of spectrally distinct third harmonics i n a photonic crystal fiber,” Opt. Lett. 26, 1158–1160 (2001). 16. N. Y. Joly, T. A. Birks, A. Yulin, J. C. Knight, and P. St. J. Russell, “Linear and nonlinear guidance in an ultralow loss planar glass membrane,” Opt. Lett. 30,2469–2471 (2005). 17. A. Efimov, A. J. Taylor, F. G. Omenetto, J. C. Knight, W. J. Wad sworth, P. St. J. Russell, “Phase-matched third harmonic generation in microstructured fibers,” Opt. Express 11,2567–2576 (2003). 18. A. A. Ivanov, D. Lorenc, F. Uherek, E. E. Serebrynnikov, S . O. Konorov, M. V. Alfimov, D. Chorvat, and A. M. Zheltikov, “Multimode anharmonic third-order harmonic gener ation in a photonic-crystal fiber,” Phys. Rev. E 73, 016610 (2006). 19. A. B. Fedotov, Ping Zhou, Yu. N. Kondrat’ev, S. N. Bagayev , V. S. Shevandin, K. V. Dukel’ski ı̂, V. B. Smirnov, A. P. Tarasevitch, D. von der Linde, and A. M. Zheltikov, “The mode structure and spectral properties of supercontinuum emission from microstructured fibers,” J. Exp. Theo r. Phys.95,851–860 (2002). 20. V. Grubsky, A. Savchenko, “Glass micro-fibers for efficien t third harmonic generation,” Opt. Express 13, 6798– 6806 (2005). 21. R. A. Sammut, A. V. Buryak, Y. S. Kivshar, “Modification of so litary waves by third-harmonic generation,” Opt. Lett. 22, 1385–1387 (1997). 22. R. A. Sammut, A. V. Buryak, Y. S. Kivshar, “Bright and dark so litary waves in the presence of third-harmonic generation,” J. Opt. Soc. Am. B 15, 1488–1496 (1998). 23. F. G. Omenetto, A. J. Taylor, M. D. Moores, J. Arriaga, J. C. Knight, W. J. Wadsworth, and P. St. J. Russell, “Simulataneous generation of spectrally distinct third har monics in a photonic fiber,” Opt. Lett. 26,1158–1160 (2001). 24. F. G. Omenetto, A. Efimov, J. Taylor, J. C. Knight, W. J. Wadsw orth, and P. St. J. Russell, “Polarization dependent harmonic generation in microstructured fibers,” Opt. Express 11,61–67 (2003). 25. A. Efimov, A. J. Taylor, “Spectral-temporal dynamics of ultr ashort Raman solitons and their role in third-order harmonic generation in photonic crystal fibers,” Appl. Phys. B 80,721–725 (2005). 26. E. E. Serebrynnikov, A. B. Fedotov, A. M. Zheltikov, A. A. Ivanov, M. V. Alfimov, V. I. Beloglazov, N. B. Skibina, D. V. Skryabin, A. V. Yulin, J. C. Knight, “Third-ha rmonic generation by Raman-shifted solitons in a photonic-crystal fiber,” J. Opt. Soc. Am. B 23,1975–1980 (2006). 27. A. N. Naumov, and A. M. Zheltikov, “Asymmetric spectral broa dening and temporal evolution of cross-phasemodulated third-harmonic pulses,” Opt. Express 10,122–127 (2002). 28. D. A. Akimov, A. A. Ivanov, A. N. Naumov, O. A. Kolevatova, M. V Alfimov, T. A. Birks, W. J. Wadsworth, P. St. J. Russell, A. A. Podshivalov, A. M. Zheltikov, “Generat ion of a spectrally asymmetric third harmonic with unamplified 30-fs Cr:forsterite laser pulses in a tapered fibe r,” Appl. Phys. B76,515–519 (2003). 29. A. Zheltikov, “Multimode guided non-3 ω third-harmonic generation by ultrashort laser pulses,” J. O pt. Soc. Am. B 22,2263–2269 (2005). 30. A. M. Zheltikov, “Third-harmonic generation with no sign al at 3ω,” Phys. Rev. A72,043812 (2005). 31. G. Genty, P. Kinsler, B. Kibler, J.M. Dudley, “Non linear envelope equation modeling of sub-cycle dynamics and harmonic generation in nonlinear waveguides,” Opt. Express 15, 5382–5387 (2007). 32. A.N. Berkovsky, S.A. Kozlov, Yu.A. Shpolyanskiy, “Self -focusing of few-cycle light pulses in dielectric media,” Phys. Rev. A72, 043821 (2005). 33. V. L. Kalashnikov, E. Sorokin, S. Naumov, I. T. Sorokina, V . V. Ravi Kanth Kumar, A. K. George, “Low-threshold supercontinuum generation from an extruded SF6 PCF using a c ompact Cr 4+:YAG laser,” Appl. Phys. B79, 591– 596 (2004). 34. E. Sorokin, S. Naumov, I.T. Sorokina, “Ultrabroadband in frared solid-state lasers,” IEEE J. Sel. Top. Quantum Electron.11, 690–712 (2005). 35. P. Xie, Zh.-Q. Zhang, “Large enhancement of third-harmoni c generation induced by coupled gap solitons in χ(3) nonlinear photonic crystals,” Phys. Rev. E 71, 026610 (2005). 36. G.P. Agrawal: Nonlinear fiber optics3rd ed., (Academic Press, New York, 2001), p 263. #83343 $15.00 USD Received 24 May 2007; revised 19 Jul 2007; accepted 15 Aug 2007; published 22 Aug 2007 (C) 2007 OSA 3 September 2007 / Vol. 15, No. 18 / OPTICS EXPRESS 11302


Introduction
Generation of the spectral supercontinuum (SC) in microstructured fibers (photonic crystal fibers, PCFs) has opened new possibilities for the applied nonlinear optics (see [1] and references therein).SC is of interest for a number of applications ranging from spectroscopy and metrology to the medicine.Possibility of the SC generation at unprecedentedly low level of the input pulse energy (sub-nanojoules) in the microstructured fibers results from the enhancement of nonlinearity and long interaction length due to the controllable fiber group-delay dispersion (GDD).The enhancement of the nonlinearity is caused by the strong mode confinement while the GDD control is made possible through the waveguide contribution to the fiber dispersion.As a result of this favorable combination, the spectrum evolution includes, besides the self-phase modulation, also the stimulated Raman scattering (SRS) [2] and generation of new components due to the four-wave mixing (FWM) [3,4], including the third-harmonic (TH) generation [8,9].
The waveguide contribution to GDD provides the phase-matched FWM and thereby enhances the efficiency of the frequency conversion [3].In particular, the required phase matching can result from a stimulated excitation of the high-order waveguide modes in conventional fibers [4] or PCFs [5,6].Such a process allows generating the third harmonic in a gas-filled hollow fiber [7,8] and a PCF [9] by both the direct conversion 3ω = ω + ω + ω and the cascading process 3ω = 2ω + 2ω − ω [10,11].It was found that the surface effects, the symmetry breaking and the Cherenkov's type phase-matching (like that in a bulk material [12]) can contribute to the TH generation in a PCF [13,14].However, the waveguide contribution to TH in the fibers is a dominant factor in contrast to the bulks.Such a contribution results in a "discreteness" of the TH spectra [15,16], meaning that the isolated TH spectral lines reproduce the discreteness of the higher-order fiber modes providing a phase-matching with the low-order mode of the pump pulse [17,18].
In spite of progress in the investigations of the TH generation in the PCFs, the theoretical analysis concerned mainly the phase-matching effects [17,14,18], while the field dynamics remained unstudied.Both analytical [14,18] and numerical [20] analysis of the TH generation by a non-dispersing pump pulse soliton have been performed.The soliton dynamics (dispersion up to the second order) in the presence of TH generation was described in [21,22].Dynamics of the pulse propagation at high peak intensities of ∼ 10 13 W/cm 2 has been investigated numerically [31,32] beyond the slowly varying amplitude approximation and taking the high-order dispersion into account.For such pulses the nonlinear interaction is very strong and results in efficient TH generation without phase-matching already at distances of 100 µm [31] to 400 µm [32].In this work we perform the analysis for the real-world conditions of cm-long fibers with high modal and material dispersion and typical low-power oscillator pulse parameters (duration ∼10 light cycles and energy ≤ 2 nJ, corresponding to peak intensities ≤ 10 12 W/cm 2 ).Such pulses are conveniently produced by a number of sources and their dynamics substantially differs from high-intensity cases of amplified pulses.We present the results of the numerical analysis of the SC dynamics in the presence of the TH generation provided by the multimode phase-matched FWM.The spatial, spectral and temporal structure of TH is studied theoretically and found to agree with the experimental data for a SF6 glass PCF pumped by a Cr 4+ :YAG oscillator operating at 60-70 fs pulse duration [33,34] As demonstrated in Ref. [33], an extruded SF6 PCF with approximately 4.5 µm rectangular core can generate a low-threshold octave-spanning SC when excited by a femtosecond pulse at λ 0 ≈1.5 µm.This is due to the proximity of the zero-GDD wavelength to 1.5 µm and very high nonlinear refraction in the SF6 material.Large nonlinearity in the combination with the strong mode confinement provides both, efficient FWM and multimode phase matching.As a result, 2ω, 3ω and mixed frequencies can be observed [33].In this work, we theoretically and numerically analyse the FWM 3ω = ω + ω + ω process in a SF6 PCF and compare it with the experiment.
A typical PCF used for the infrared SC generation is shown in Fig. 1.The mode structure of PCF was analyzed on the basis of the full-vector model using the real core geometry [33].The fiber is multimode in the infrared, but only the lowest modes M=1-4 (Fig. 1) can contribute to the SC and TH generation because of large shift of the zero-dispersion wavelength for the higher-order modes [33] from λ 0 .The rectangular slightly asymmetric PCF core geome- try gives rise to slightly different mode parameters depending on polarization [33].This should result in the polarization effects, but we were unable to observe any noticeable birefringence in the experiment, neither in the SC nor in the TH signal.We suppose that the core inhomogeneity along the fiber (evidenced by the fluctuations of the visible light emission) mixes the polarizations within the first millimeters of propagation.As a result of the waveguiding contribution to the effective refractive index, it is possible to satisfy the phase matching condition: n eff,M (ω) = n eff,N (3ω), where M and N are the mode indices and M < N. Fig. 2 shows the effective refractive indices of the lowest-order modes in the vicinity of λ 0 (M =1, 2) and the high-order modes in the vicinity of λ 0 /3 (N =69-80) as a function of wavelength for the SF6 PCF with a 2.5 µm core.The mode enumeration rule corresponds to Ref. [33], so that the two lowest-order modes result from the two orthogonal polarizations of the zero-order PCF mode.The phase matching points are shown by the red circles.2) and at the third harmonic (≈ λ 0 /3, M = 69-80).Red circles show the phase matching points.SF6 PCF with 2.5 µm core.

M=1 M=3 M=74
The experimental TH spectra were obtained from a 4 cm section of the SF6 PCF pumped by the Cr 4+ :YAG oscillator generating 60-70 fs pulses at 100 MHz repetition rate with the average output power of up to 170 mW [33].The typical observed infrared and TH spectra are shown in Figs.3,4.Fig. 4 demonstrates the TH spectrum together with its decomposition according to the multimode phase matching scheme presented in Fig. 2. The TH spectrum can be interpreted as a multimode composite, excited by the FWM process, which is phase-matched with the infrared SC generated within the lowest modes.
Since the multimode phase matching mechanism provides a good description of the TH spectral structure, one can suppose that the intra-fiber power is not too high or/and that the TH radiation losses are sufficiently strong in order to suppress the contribution of the crossphase modulation to the mode matching.In the general case, such a conclusion would not be valid [27] and one can see the effect of the cross-phase modulation on the TH spectrum in the numerical simulation below.
Intensities of the TH spikes are defined by several factors including the SC dynamics at the fundamental wavelength as discussed in the next section.Another important contribution comes from the mode spatial structure, because various modes will have different overlapping with the pump mode (affects FWM efficiency) and different waveguide radiation losses (propagation attenuation).Both overlapping factor and radiation losses were found numerically on the basis of the generalized eigenvalue method [33].The overlapping factors of the most intense modes m 74 (Fig. 1) and m 75 are 12% and 14%, and the radiation losses are 9 and 20 dB/m, respectively.As another example, the radiation loss for the weaker mode m 71 with comparable overlapping factor is ≈2000 dB/m (leaking mode).Another leaking mode m 73 contributes in the right maximum spike together with m 74 .One can conclude that the TH spectrum shape depends on both overlapping factor and radiation losses for the excited modes.
PCF with 2.5 µm core provides a higher effective nonlinearity due to its small effective mode area (A eff =3.7 µm 2 for the lowest-order mode).As a result, the FWM efficiency is higher than in Fig. 4. Experimental TH spectrum from the SF6 PCF with a 2.5 µm core and its mode decomposition according to from Fig. 2. The spectrometer resolution was 0.5 nm the case of the larger core sizes (e.g.A eff =9.4 µm 2 for the lowest-order mode of a PCF with 4.5 µm core).This fact is easily observable in the experiment.However, since the zero-dispersion wavelength is λ ≈1.2 µm for the small core PCF, the SC spectrum is smooth and comparatively narrow [33].The 4.5 µm core PCF produces TH with lower intensity, but i) its density of highorder modes is higher around λ 0 /3 (Fig. 5) and ii) its SC is wider and more modulated due to the proximity of the zero-dispersion wavelength to λ 0 .The latter is most interesting for an analysis because it affects the TH temporal structure and dynamics.

Spatial, spectral and temporal structures of TH
For the 4.5 µm core PCF the zero-dispersion wavelength is close to λ 0 ≈ 1.5 µm.This provides a comparatively wide SC shown in Fig. 6.Its spectral structure can be interpreted as consisting of i) the soliton spike centered at 1.6 µm, ii) the auto-modulated remainder of the pump pulse at 1.5 µm (narrow line is the cw component from the oscillator), iii) the dispersion wave at ≈1.26 µm, and iv) the TH spikes in the vicinity of 505 and 540 nm.
Since the TH signal is generated by the broadband SC pump pulse, many modes are generated simultaneously.This results in a broad almost featureless spectrum (Fig. 6), as well as Wavelength, nm Spectral power, arb.un.smooth structureless mode pattern at the output (Fig. 7), which is an overlap of many mode patterns (Fig. 1) with close wavelengths.Such an overlap develops already at the initial propagation stage and the individual modes cannot be resolved.In similar experiments with fused silica-based PCFs and tunable narrowband pump pulses one can observe and resolve the individual phase-matched high-order modes [17], which are also conveniently separated from each other due to the lower mode density in the 2.5 µm-core fused silica as opposed to the 4.5 µm-core SF6.
To analyze the SC dynamics numerically, we reduce the full 3+1 dimensionality to 1+1dimensional problem for each pair of modes (one lowest-order pump mode and one high-order mode at λ 0 /3).We also approximate the slightly birefringent PCF core by a core with cylin- drical symmetry and 5.2 µm diameter, providing the same zero-dispersion wavelength at ≈ 1.5 µm.Spatial profiles of the three high-order modes providing the phase-matching with the infrared SC are shown in Fig. 8.The mode structure of the cylindrical core differs from that of the experimental PCF.Therefore, the TH spectral structure can be reproduced only approximately.In particular, the birefringence effects disappear and there are 4 times less phase-matching points (only three in the wavelength interval 450-600 nm).In the rest of the paper we concentrate on the TH dynamics and its co-evolution with the supercontinuum.
The model which describes an interaction between the fundamental (SC, a F ) and third harmonic (a T H ) fields is based on the slowly varying amplitude approximation, which is valid for our pulse width (the case of the extremely short pulses is considered in [31,32]).The equations consider the fields separately but with the coupling through FWM [8,21,22,20,35]: Here z is the propagation distance, t is the local time ( frame of reference is co-moving with the group velocity v F of the fundamental field at ω , n 2 is the nonlinear refractive index, c is the velocity of light.A F eff and A T H eff are the effective mode areas at ω F and ω T H , respectively.Parameter β 2 describes the frequency-dependent GDD in the vicinity of ω F (first equation) or ω T H (second equation).Parameters ν F and ν T H are the group velocities at ω F and ω T H , respectively.The overlapping of modes [36] at ω F and ω T H defines the parameter ξ .TH initial filed is assumed null a T H (z = 0) = 0.
The spectral profile of β 2 (ω) reproduces the fiber dispersion and includes the contribution of high-order dispersions up to the 12 th order.Taking into account of the higher-order dispersions is necessary to describe the SC structure and its dynamics, which affects the TH.
Numerical simulations were performed pairwise between the fundamental and each higherorder modes.Such an approach is valid in our case because the energy conversion into TH is small (Fig. 11).The simulations confirm that the signal in the higher-order modes appears in the vicinity of the tripled pump wavelength at λ ≈ λ 0 /3 and ∆k ≈ 0 [11].As new wavelengthes are be- ing generated in the infrared SC along the propagation coordinate (within both solitonic and dispersion-wave parts of SC), the phase-matched TH higher-order modes grow in the form of isolated spectral spikes [18].Fig. 9 shows three modes with the highest simulated intensity.They correspond to initial ∆k close to 0. The spikes have widths of the order of few nanometers.There is a spectral overlap between some TH modes.The spectra have a fine structure and asymmetry [27,28], resulting from the cross-phase modulation with the fundamental wave.The cross-phase modulation is hardly visible in the experiment because of the radiation losses.
In the numerical simulation, where we neglect the radiation losses, the cross-phase modulation effect is stronger.Fig. 10 shows the instantaneous snapshot of the co-propagating signals at a fiber length of 4 cm.Because of the large group velocity difference, the TH wave is constantly retarded from the SC pulse.In the calculation, the local time frame is moving with a constant velocity, which is equal to the group velocity of the input pulse at 1.5 µm, and is close to the group velocity of the SC signal.Plotting the temporal intensity profile of the TH signal in the local time frame thus shows a generation history of the TH signal during the propagation.
It is interesting to note, that at the beginning of the TH generation there are ripples with a period of the order of ∼ 100 fs (Fig. 10b).The reason for this can be understood, if we take into account, that at the initial stage, the frequency of the TH signal is exactly ω T H = 3ω F .For most of the high-order modes, this is not the phase-matched frequency, and the signal shows typical oscillations of the non-phase-matched harmonic generation with ∆k = 0. Later on, in the first few millimeters of propagation [34], the spectrum expansion takes place, and the TH generation frequency shifts to the perfect phase-matching ∆k = 0.
The long-term modulation of the TH temporal profile (≈ 10 ps, Fig. 8a) results from the variation of the SC spectral power at the phase-matched wavelengths along the propagation distance z.A combination of the short-scale (due to ∆k = 0) and long-scale (due to the SC evolution) modulations creates the unique temporal structure of the TH signal.Note that the long-scale evolution is mode-dependent (Fig. 8), because it is connected with the SC intensity evolution at the tripled wavelength of the mode-specific phase-matching point.
Finally, Fig. 11 shows the evolution of the overall energy in the TH modes as a fraction of the input pulse energy.At 4 cm propagation length, the maximum conversion efficiency is ∼0.8%.This value is likely the upper estimation of the conversion efficiency, because the model does not account for the radiation losses of the high-order TH modes.

Conclusion
Multimode phase-matched four-wave-mixing providing the third-harmonic generation in microstructured fibers has been explored numerically taking into account the propagation dynamics within the framework of slowly-varying amplitude approximation.It has been found, that the the TH signal is a superposition of spectral peaks, which are excited in the high-order spatial modes providing the phase matching with the supercontinuum in the PCF.The calculated TH peaks and SC spectra are in a good agreement with the experimental ones, which have been obtained from the SF6 glass PCF pumped by a femtosecond pulse Cr 4+ :YAG oscillator.The structure of the experimental TH spectrum is completely described by the multimode phasematching conditions in the actual geometry of the PCF core.
Temporal structure of the TH has been analyzed in connection with the SC dynamics.The TH appears in the vicinity of the tripled pump frequency in the form of isolated spectral lines propagating in the modes with ∆k = 0. Imperfect phase-matching causes the TH envelope ripples with the sub-picosecond period.Such ripples decay when the SC approaches the perfect phase-matching frequency, where ∆k → 0. After that the temporal structure of TH is affected mainly by the dynamics of SC.The resulting TH signal has the form of a intensity-modulated pulse with approximately 10-ps duration, which stretches out of the SC pulse due to the differ-

Fig. 1 .
Fig. 1.Cross section of the SF6 PCF (upper graph) and field intensity distributions for the two lowest infrared modes and one visible mode (lower graph, see text for mode designations).

Fig. 6 .Fig. 7 .
Fig.6.Experimental spectra of the SC (red) and TH (black) from SF6 PCF with 4.5 µm core.Blue curve shows the input pulse spectrum.All three spectra are normalized separately for better visibility.

Fig. 8 .
Fig. 8. Spatial profiles of the three TH modes in the cylindrical 5.2 µm core of the SF6 PCF.

Fig. 9 .
Fig. 9. Calculated: the SC (a) and the TH (b) spectra for the modes of Fig. 8.

Fig. 10 .
Fig. 10.Calculated: SC (fundamental) and TH temporal profiles for the modes of Fig. 6 at 4 cm propagation length.(a) Full time window of the TH signal (log scale).(b) Trailing edge of the TH pulse, corresponding to the beginning of the TH generation (linear scale).

Fig. 11 .
Fig. 11.Calculated evolution of the total TH energy fraction for the modes of Fig. 8.

. Multimode phase matching in a SF6 glass PCF
.