## Abstract

We investigate compression of ultrashort laser pulses by nonlinear propagation in gas-filled planar hollow waveguides, using (3+1)-dimensional numerical simulations. In this geometry, the laser beam is guided with a fixed size in one transverse dimension, generating significant spectral broadening, while it propagates freely in the other, allowing for energy up-scalability. In this respect the concept outperforms compression techniques based on hollow core fibers or filamentation. Small-scale self-focusing is a crucial consideration, which introduces mode deterioration and finally break-up in multiple filaments. The simulation results, which match well with initial experiments, provide important guidelines for scaling the few-cycle pulse generation to higher energies. Pulse compression down to few-cycle duration with energies up to 100mJ levels should be possible.

© 2009 Optical Society of America

## 1. Introduction

One of the main goals of today’s ultrafast nonlinear optics is the compression of ultrashort laser pulses down to the few-cycle regime. Intense few-cycle infrared (IR) laser pulses with reproducible carrier envelope offset are key ingredients in attosecond physics, serving both as driving pulses for the generation of intense coherent sources in the extreme ultraviolet (XUV) wavelengths [1] and as short intense streaking pulses in the near infrared [2, 3]. Near single-cycle pulses are commonly generated through nonlinear propagation in gas-filled hollow fibers [4] or through optical filamentation in noble gases [5, 6, 7, 8]. Due to potential damage to the fiber, the hollow fiber scheme supports pulse energies typically limited to sub-millijoules, in particular when driven at kHz-repetition rates. At low repetition rate and by using a pressure gradient along the fiber the output was recently increased to some millijoules of pulse energy [9]. Pulse compression through filamentation is limited to the same energy levels, caused by intensity clamping, inherent chirp and angular dispersion within the filament [10, 11, 12] or by the onset of multiple filaments. To date, few-cycle pulses with higher energies (up to ~100mJ) are mainly generated through optical parametric chirped pulse amplification (OPCPA) [13, 14]. These setups involve considerable experimental complexity and require the use of additional pump lasers.

Recently, the theoretical work of Nurhuda *et al*. [15] suggested a compression scheme using gas-filled planar hollow waveguides, which can address the issue of energy upscalability with significantly less experimental complexity. As the beam is guided in only one transverse direction (short axis), the perpendicular direction (long axis) can be utilized to freely adjust the beam size and intensity with the objective of limiting photo ionization on the one hand, and achieving strong self-phase modulation (SPM) and spectral broadening, essential to reach the few-cycle regime, on the other hand. While the simulation results presented by Nurhuda *et al*. [15] are promising, the initial experiments by Chen *et al*. and Akturk *et al*. showed that the transverse spatial propagation dynamics imposes severe restrictions on the energy up-scalability [16, 17]. Imperfections and noise in the input spatial mode are amplified by small-scale self-focusing, resulting in output beams with complicated spatial structure and multiple filaments. Nonetheless, the planar waveguide scheme holds its promise for high energy pulse compression, but a better understanding and more practical and complete analysis of nonlinear pulse propagation inside the waveguide is needed.

In this work, we present a comprehensive analysis of nonlinear propagation and pulse compression via planar hollow waveguides. In practice a planar hollow waveguide can be composed of two polished glass surfaces facing each other and separated by thin spacers. It should be noted that this structure strictly speaking is not a planar waveguide, since the modes are lossy. However, the losses are small for the thickness of the spacers considered. For simplicity we keep the term waveguide throughout the paper. We present the results of detailed (3+1)-dimensional simulations of the spatio-temporal pulse dynamics in the waveguide, which provides understanding of the energy and compressibility limits of this technique, from which a practical criterion to reach the best pulse compression without compromising the spatial mode is established. Careful validation of our modeling is achieved by comparing simulated compressed pulse temporal shape and spatial mode with experimental results [17]. The model is then applied to define the conditions for the compression of very energetic pulses (>100mJ) down to the few-cycle regime (<10fs) as well as for the possible compression of few-cycle pulses with more than 15mJ energy to the single-cycle regime.

## 2. Numerical simulation results and experimental comparison

Nonlinear propagation inside the gas-filled planar hollow waveguide is modeled by the Nonlinear Envelope Equation (NEE) [18], known to be valid to the single-cycle regime [19]. The description of the interaction of the ultrashort pulse with the argon gas inside the waveguide comprises a Kerr-type nonlinearity for self-focusing and self-phase modulation (nonlinear index of refraction *ñ*
_{2}=*n*
_{2}/*P*=0.98×10^{−19} cm^{2}W^{−1}bar^{−1} [20]), plasma generation via optical field ionization with rates calculated from the generalized Keldysh-PPT formula [18, 21], plasma defocusing, as well as higher order dispersion, self-steepening, and space-time focusing. The NEE is numerically integrated in a split-step method, treating the nonlinear part via the field dependent nonlinear polarization and the linear part in mode space, allowing us to properly describe exact coupling between high-order modes, dispersion and losses induced by the waveguide as well as guided propagation along the short waveguide axis and free propagation in the perpendicular transverse direction. For this purpose, the electric field is decomposed into a basis of leaky transverse electric (TE) waveguide even modes *V _{p}*=2

*m*(

*x*)=sin (

*pπx*/2

*a*) and odd modes

*V*=2

_{p}*m*−1 (

*x*)=cos (

*pπx*/2

*a*). The transformation from one representation to the other reads:

where *E* is the complex envelope of the field and *A _{p}* of the leaky mode with index

*p*; 2

*a*=127

*µ*m is the separation of the glass slabs (refractive index

*n*=1.45) constituting the waveguide. The dependence upon frequency of propagation and attenuation constants

*β*(

_{p}*ω*)=

*k*(1-(

*pπ*/2

*ka*)

^{2})

^{1/2}with

*β*

_{p,0}=

*β*(

_{p}*ω*

_{0}) and ${\alpha}_{p}\left(\omega \right)={p}^{2}{\pi}^{2}/4{k}^{2}{a}^{3}\sqrt{{n}^{2}-1}$ is obtained from the wave number

*k*(

*ω*) calculated from the dispersion data of argon [22], which also depends on pressure. In the linear part, the 2-dimensional Fourier transform

*Ã*(

_{p}*k*) is propagated, thus no expansion around ω0 is made in the dispersion, absorption, and group velocity for each mode. Higher order modes decay much faster compared to the fundamental. Ten modes are sufficient to numerically propagate ultrashort pulses inside the waveguide.

_{y},ω,zThe parameters for the calculations presented in Fig. 1 and 2 were chosen exactly as in the experiments of reference [17], namely a pulse duration of 43fs (FWHM), a beam size in the guided direction of *w _{x}*=0.7

*a*close to the optimum coupling conditions, a beam size of

*w*=7.2mm along the long waveguide axis, a pressure of argon of 1.5atm, and a waveguide length of 21cm (see also Tab. 1). Figure 1 shows the normalized power inside the waveguide compared to the power of the fundamental mode as well as higher order modes. Nonlinear coupling mainly occurs between the fundamental (

_{y}*p*=1) and the first excited mode (

*p*=3), whereas the length scale of the power oscillations corresponds to the beating length

*L*≡2

_{b}*π*/(

*β*

_{1,0}−

*β*

_{3,0})~1cm. Only a small portion of power is exchanged between higher order modes. In the simulations, the total throughput of the input energy is about 93% compared to 65% in the experiment. Different effects, which are not considered in the modeling, may cause this discrepancy, such as coupling losses due to inhomogeneities in the beam, imperfections of the waveguide, reflections from the windows of the gas cell, and beam clipping inside the cell. To better match the experimental conditions, the input pulse energy for the simulations was corrected for the additional loss; furthermore the measured input temporal intensity and phase was used. To mimic the experimental compression by negatively dispersive mirrors, the output pulse obtained from the simulations is compressed by applying second order phase delay only, leaving higher order chirp, picked up during nonlinear interaction, uncorrected. To also diagnose the output spatial mode quality, which is a crucial parameter for energy up-scalability, a measured beam profile with 4% intensity modulations compared to an ideal Gaussian profile was used as input along the long waveguide axis.

Figure 2 compares the experimental results of reference [17] with the modeling, showing convincing overall agreement. Note in particular that the total spectral and temporal widths coincide well in the spectrograms (2A,B), generated by means of frequency resolved optical gating (FROG) [10]. The phase of simulated broadened spectra (data not shown) features a strong second-order component, explaining the good compressibility with chirped mirrors. Finally, also the transverse structure of the experimental output mode is numerically reproduced, provided a measured input beam profile is considered in the long transverse direction (2E,F). However, the depth of the transverse modulations is somewhat deeper experimentally. This could be due to additional imperfections of the beam and the waveguide, not accounted for in the simulations. Asymmetries in the guided direction are filtered by the waveguide (2E). The intensity of *I*
_{0}=3.35×10^{13}Wcm^{−2} keeps the density of generated free carriers well below 10^{15} cm^{−3} (data not shown), minimizing propagation losses and limiting the nonlinear interaction to mainly self-phase modulation.

#### 2.1. Analysis of the beam stability

The simulations clearly reveal that pulse compression in planar hollow waveguides is compromised by the deterioration of the output mode due to small-scale self-focusing (see Fig. 3, lower plot). To identify a regime of stable propagation, one must evaluate the length scales for instabilities to be avoided, namely *L _{sf}* for self-focusing of a beam of width wy along the long transverse axis and

*L*for small-scale self-focusing (modulational instability) [23]. An analytical expression for

_{mi}*L*(dashed line in Fig. 3) is obtained by deriving an evolution equation for the second order moment

_{sf}*w*

^{2}(

*z*)≡∫

^{+∞}

_{−∞}

*y*

^{2}|

*E*|

^{2}

*dy*/∫

^{+∞}

_{−∞}|

*E*|

^{2}

*dy*[18]. Assuming whole beam Kerr self-focusing with a frozen beam profile,

*e.g*.

*E*∝exp(−

*y*

^{2}/

*w*

^{2}(

*z*)), a direct integration between the entrance of the waveguide and the nonlinear focus can be performed which yields:

*L*=2

_{sf}*πz*(

_{Ry}*P**/

*P*)×[4(

_{c}*P**/

*P*−1)]

_{c}^{−3/2}, where

*z*≡

_{Ry}*k*

_{0}

*w*

^{2}

_{y}/2 is the Rayleigh length associated with

*w*,

_{y}*P**=

*πI*

_{0}

*w*

^{2}

_{y}/2 the power of a 2D beam with the same intensity and width and

*P*≡2

_{c}*πn*

_{0}/

*k*

^{2}

_{0}

*n*

_{2}is the critical power for 2D collapse. For

*P**/

*P*≫1 and constant intensity,

_{c}*L*is roughly proportional to wy. Small-scale self-focusing follows the theory of modulational instability [24] from which an evaluation of

_{sf}*L*≡log(

_{mi}*G*)

*k*

_{0}

*P*/2

_{c}*πI*

_{0}is found, where G denotes a gain that represents the amplification factor of small-scale inhomogeneities (both beam imperfections and generic noise) before saturation of the modulational instability in the form of filaments.

In Fig. 3 the upper plot shows the lengths after which either a nonlinear focus or beam break-up is numerically observed, as we scale the transverse beam size *w _{y}* (long axis) and pulse energy, respectively at constant input intensity. Simulations with a perfect Gaussian transverse input beam profile along the long axis (indicated as triangles) greatly overestimate the distance of stable propagation as compared to simulations with a realistic (measured) profile used as an input condition (indicated as circles). The shaded area identifies the regime of stable propagation determined by the measured beam profile. For small

*w*the propagation distance is limited by self-focusing of either the whole beam or beam imperfections; for

_{y}*w*≳20mm modulational instability defines the limit. The solid star corresponds to the experimental conditions of reference [17] and the simulations presented in Fig. 2, showing that the length of the waveguide (dotted line) was at the border of the permissible propagation distance. For these exact conditions the fluence along the waveguide is shown in the lower plot. Indeed, shortly after the length of the experimental waveguide (dotted line) the profile breaks into individual filaments, unusable for compression. Note that the length scales and the stability regime in Fig. 3 are valid for a particular combination of parameters only (

_{y}*I*

_{0}

*n*

_{2}=const). However, when modulational instability and not self-focusing determines the stability region, the distance of stable propagation can also be associated with a specific value of the B-integral

*B*~

*k*

_{0}

*I*

_{0}

*n*

_{2}

*z*(upper plot, solid line), where z is the propagation distance. Its value

*B*~15.5 is solely determined by the amount of imperfections and noise contained in the measured input profile. This corresponds to a gain

*G*~4.5×10

^{4}. Experimentally, due to additional sources of noise, modulational instability may saturate below

*B*~15.5. For the parameters of reference [17] the B-integral corresponds to a value of

*B*~8.

#### 2.2. Energy upscalability

Having defined the stability region for hollow planar waveguide pulse compression, we are now in a position to determine the waveguide and laser beam parameters, allowing to obtain higher pulse energies and shorter pulse durations. We particularly examine two cases: First, the compression of very energetic pulses (>100mJ) to durations of sub 10fs; and second, the compression of pulses, which already are in the few-cycle regime, to further advance to the single-cycle regime. The parameters are listed in Tab. 1. For the sake of simplicity we assume Gaussian spatial and temporal beam profiles and choose conditions, which can also be indicated in Fig. 3 (open stars). However, due to the different pulse durations only the length scale (lower axis) but not the energy scale (upper axis) is applicable. The compressed temporal output profiles for both cases are shown in Fig. 4. For the high energy pulse we calculate a compression from 40fs duration to below 10fs at 100mJ output energy (left column). The short pulse can be compressed from 10fs to below 3fs (almost a single-cycle pulse) at energy of 15.8mJ (right column), provided chirped mirrors of sufficient bandwidth (500–1200nm) are available. The 2D plots of the compressed pulse profiles (upper row) reveal a homogeneous compression close to the axis, whereas the pulse becomes longer towards the edge of the profile due to less self-phase modulation. For the case of the 100mJ pulse (Fig. 4, left column) about 70% of the pulse energy falls within a duration of 1.5 times the duration at the center. Simultaneously, the compressed pulse features a phase and pulse front delay in the center. Therefore, to achieve good focusability, the spatial phase has to be either pre or post compensated by, *e.g*. adaptive optics.

## 3. Conclusion

In conclusion, we have investigated the potential of planar hollow waveguides for the compression of ultrashort laser pulses in practical experimental conditions. Stability analysis revealed that the limit of possible compression is given by a B-integral threshold determined by the input spatial beam profile. Nonetheless, within the regime of stable propagation the scheme scales up to high energies and can as well be applied to advance from the few-cycle to the single-cycle regime. In this respect a two-stage compression seems feasible, with the additional benefit of also employing the waveguide in the second stage to filter spatial noise generated within the first stage.

## Acknowledgements

We gratefully acknowledge the funding of C.L. Arnold by the Deutsche Akademie der Natur-forscher Leopoldina, Grant No. BMBF-LPD 9901/8-181.

## References and links

**1. **M. Hentschel, R. Kienberger, C. Spielmann, G. A. Reider, N. Milosevic, T. Brabec, P. Corkum, U. Heinzmann, M. Drescher, and F. Krausz, “Attosecond metrology,” Nature **414**, 509–513 (2001).
[CrossRef] [PubMed]

**2. **R. Kienberger, E. Goulielmakis, M. Uiberacker, A. Baltuska, V. Yakovlev, F. Bammer, A. Scrinzi, T. Wester-walbesloh, U. Kleineberg, U. Heinzmann, M. Drescher, and F. Krausz, “Atomic transient recorder,” Nature **427**, 817–821 (2004).
[CrossRef] [PubMed]

**3. **P. Eckle, M. Smolarski, P. Schlup, J. Biegert, A. Staudte, M. Schöffler, H. G. Muller, R. Dörner, and U. Keller, “Attosecond angular streaking,” Nature Physics **4**, 565–570 (2008).
[CrossRef]

**4. **M. Nisoli, S. de Silvestri, O. Svelto, R. Szipöcz, K. Ferencz, C. Spielmann, S. Sartania, and F. Krausz, “Compression of high-energy laser pulses below 5 fs,” Opt. Lett. **22**, 522–524 (1997).
[CrossRef] [PubMed]

**5. **C. P. Hauri, W. Kornelis, F. W. Helbing, A. Heinrich, A. Couairon, A. Mysyrowicz, J. Biegert, and U. Keller, “Generation of intense, carrier-envelope phase-locked few-cycle laser pulses through filamentation,” Appl. Phys. B **79**, 673–677 (2004).
[CrossRef]

**6. **A. Couairon, M. Franco, A. Mysyrowicz, J. Biegert, and U. Keller, “Pulse-compression to the single cycle limit by filamentation in a gaz with a pressure gradient,” Opt. Lett. **30**, 2657–2659 (2005).
[CrossRef] [PubMed]

**7. **A. Couairon, J. Biegert, C. P. Hauri, W. Kornelis, F. W. Helbing, U. Keller, and A. Mysyrowicz, “Self-compression of ultrashort laser pulses down to one optical cycle by filamentation,” J. Mod. Opt. **53**, 75–85 (2006).
[CrossRef]

**8. **A. Mysyrowicz, A. Couairon, and U. Keller, “Self-compression of optical laser pulses by filamentation,” New J. Phys. **10**, 025023 (2008).
[CrossRef]

**9. **A. Suda, M. Hatayama, K. Nagasaka, and K. Midorikawa, “Generation of sub-10-fs, 5-mJ-optical pulses using a hollow fiber with a pressure gradient,” Appl. Phys. Lett. **86**, 111116 (2005).
[CrossRef]

**10. **S. Akturk, A. Couairon, M. Franco, and A. Mysyrowicz, “Spectrogram representation of pulse self compression by filamentation,” Opt. Express 17626–17636 (2008).
[CrossRef] [PubMed]

**11. **H. S. Chakraborty, M. B. Gaarde, and A. Couairon, “Single attosecond pulses from high harmonics driven by self-compressed filaments,” Opt. Lett. **31**, 3662–3664 (2006).
[CrossRef] [PubMed]

**12. **A. Couairon, H. S. Chakraborty, and M. B. Gaarde, “From single-cycle self-compressed filaments to isolated attosecond pulses in noble gases,” Phys. Rev. A **77**, 053814–053824 (2008).
[CrossRef]

**13. **F. Tavella, A. Marcinkevicius, and F. Krausz, “90 mJ parametric chirped pulse amplification of 10 fs pulses,” Opt. Express **14**, 12822–12827 (2006).
[CrossRef] [PubMed]

**14. **A. Dubietis, R. Butkus, and A. P. Piskarskas, “Trends in Chirped Pulse Optical Parametric Amplification,” IEEE J. SEL. Top. Quantum Electron. **12**, 163–172 (2006).
[CrossRef]

**15. **M. Nurhuda, A. Suda, S. Bohman, S. Yamaguchi, and K. Midorikawa, “Optical Pulse Compression of Ultrashort Laser Pulses in an Argon-Filled Planar Waveguide,” Phys. Rev. Lett. **97**, 153902 (2006).
[CrossRef] [PubMed]

**16. **J. Chen, A. Suda, E. J. Takahashi, M. Nurhuda, and K. Midorikawa, “Compression of intense ultrashort laser pulses in a gas-filled planar waveguide,” Opt. Lett. **33**, 2992–2994 (2008).
[CrossRef] [PubMed]

**17. **S. Akturk, C. L. Arnold, B. Zhou, and A. Mysyrowicz, “High energy ultrashort laser pulse compression in hollow planar waveguides,” Opt. Lett. **34**, 1462–1464 (2009).
[CrossRef] [PubMed]

**18. **A. Couairon and A. Mysyrowicz, “Femtosecond filamentation in transparent media,” Phys. Rep. **441**, 47–189 (2007).
[CrossRef]

**19. **T. Brabec and F. Krausz, “Nonlinear optical pulse propagation in the single-cycle regime,” Phys. Rev. Lett. **78**, 3282–3285 (1997).
[CrossRef]

**20. **G. Tempea and T. Brabec, “Theory of self-focusing in a hollow waveguide,” Opt. Lett. **23**, 762–764 (1998).
[CrossRef]

**21. **A. M. Perelomov, V. S. Popov, and M. V. Terent’ev, “Ionization of atoms in an alternating electric field,” Sov. Phys. JETP **23**, 924–934 (1966).

**22. **Dalgarno and Kingston, “The refractive indices and Verdet constants of the inert gases,” Proc. Roy. Soc. London Ser. A **259**, 424–429 (1966).

**23. **G. Fibich, S. Eisenmann, B. Ilan, Y. Erlich, M. Fraenkel, Z. Henis, A. L. Gaeta, and A. Zigler, “Self-focusing distance of very high power laser pulses,” Opt. Express **13**, 5897–5903 (2005).
[CrossRef] [PubMed]

**24. **V. I. Bespalov and V. I. Talanov, “Filamentary structure of light beams in nonlinear liquids,” Zh. Eksper. Teor. Fiz. Pis’ma 3, 471–476 (1966). [JETP Lett. 3 (1966) 307–310].