Transverse-mode dependence of femtosecond filamentation

doi:10.1364/OE.17.012217

« Show journal navigation

Transverse-mode dependence of femtosecond filamentation

Zhenming Song, Zhigang Zhang, and Takashi Nakajima »View Author Affiliations

Optics Express, Vol. 17, Issue 15, pp. 12217-12229 (2009)
http://dx.doi.org/10.1364/OE.17.012217

View Full Text Article

Acrobat PDF (422 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Model
3. Results and discussions
4. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (38)
Cited By
Figures (12)
Metrics

Abstract

We theoretically investigate the transverse-mode dependence of femtosecond filamentation in Ar gas. Three different transverse modes, Bessel, Gaussian, and Laguerre modes, are considered for incident laser pulses. By solving the extended nonlinear Schrödinger equation coupled with the electron density equation, we find that the lengths of the filament and the plasma channel induced by the Bessel incident beam is much longer than the other transverse modes with the same peak intensity, pulse duration, and beam diameter. Moreover we find that the temporal profile of the pulse with the Bessel incident mode is nearly undistorted during the propagation. Since the pulse energy that the Bessel beam can carry is more than one order of magnitude larger than the other modes for the same peak intensity, pulse duration, and beam diameter, the Bessel beam can be a very powerful tool in ultrafast nonlinear optics involving propagation in a Kerr medium.

1. Introduction

It is well known that the propagation of intense laser pulses in any kinds of Kerr media (gases, liquids, and transparent solid materials) can result in the formation of filamentation under certain conditions. Recently various aspects of filamentation, such as long range propagation [1

G. Méchain, C. D’Amico, Y.-B. André, S. Tzortzakis, M. Franco, B. Prade, A. Mysyrowicz, A. Couairon, E. Salmon, and R. Sauerbrey, “Length of plasma filaments created in air by a multiterawatt femtosecond laser,” Opt. Commun. 247, 171–180 (2005). [CrossRef]

], light detection and ranging techniques (LIDAR) [2

P. Rairoux, H. Schillinger, S. Neirdeimer, M. Rodriguez, F. Ronneberger, R. Sauerbrey, B. Stein, D. Waite, C. Wedekind, H. Wille, L. Wöste, and C. Ziener, “Remote sensing of the atmosphere using ultrashort laser pulses,” Appl. Phys. B 71, 573–580 (2000). [CrossRef]

], virtual antenna [3

G. Méchain, A. Mysyrowicz, M. Depiesse, and M. Pellet, “A virtual antenna produced in air by intense femtosecond laser pulses,” In: D.H. Titterton, S.M. Kirkpatrick, R. Stoian, R. Appleby, J.M. Chamberlain, and K.A. Krapels, (Eds.), Technologies for Optical Countermeasures II; Femtosecond Phenomena II; and Passive Millimetre-Wave and Terahertz Imaging II, vol. 5989. SPIE , p. 5989 (2005).

], spectral broadening [4

J. Kasparian, R. Sauerbrey, D. Mondelain, S. Niedermeier, J. Yu, J. -P. Wolf, Y. -B. André, M. Franco, B. Prade, S. Tzortzakis, A. Mysyrowicz, M. Rodriguez, H. Wille, and L. Wöste, “Infrared extension of the supercontinuum generated by femtosecond terawatt laser pulses propagating in the atmosphere,” Opt. Lett. 25, 1397–1399 (2000). [CrossRef]

], pulse self-compression [5

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]

] and THz generation [6

S. Tzortzakis, G. Méchain, G. Patalano, Y.-B. André, B. Prade, M. Franco, A. Mysyrowicz, J.-M. Munier, M. Gheudin, G. Beaudin, and P. Encrenaz, “Coherent subterahertz radiation from femtosecond infrared filaments in air,” Opt. Lett. 27, 1944–1946 (2002). [CrossRef]

], have attracted lots of attention. Interestingly, recent experimental as well as theoretical studies have shown that the filamentation can be formed without self-guiding [7

A. Dubietis, E. Gaižauskas, G. Tamošauska,s, and P. Di. Trapani, “Light filaments without self-channeling,” Phys. Rev. Lett. 92, 253903 (2004). [CrossRef] [PubMed]

By recalling that the filamentation is an outcome of the balance of the diffraction, self-focusing due to the nonlinear index refraction and self-defocusing due to ionization, it is quite natural to come to the idea that the use of different transverse modes for incident laser pulses may result in different propagation dynamics, which may affect the filamentation formation, etc. Indeed, because the Bessel beam [8

J. Durnin, J. J. Miceli Jr, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987). [CrossRef] [PubMed]

, 9

J. Durnin, “Exact solutions for nondiffracting beams. I. The scalar theory,” J. Opt. Soc. Am. A 4, 651–654 (1987). [CrossRef]

] is known to be diffraction-free when it propagates in vacuum, the use of the Bessel beam draws lots of interests in recent years [10

P. Polesana, M. Franco, A. Couairon, D. Faccio, and P. Di Trapani, “Filamentation in Kerr media from pulsed Bessel beams,” Phys. Rev. A 77, 043814 (2008). [CrossRef]

-12

A. Dubietis, P. Polesana, G. Valiulis, A. Stabinis, P. Di Trapani, and A. Piskarskas, “Axial emission and spectral broadening in self-focusing of femtosecond Bessel beams,” Opt. Express 15, 4168–4175 (2007). [CrossRef] [PubMed]

]. In a laboratory, however, it is not exactly an ideal Bessel beam but rather a quasi-Bessel beam that is produced. A quasi-Bessel beam is an apertured Bessel beam which carries large but still finite energy. Compared with a Gaussian beam with the same beam diameter, even a quasi-Bessel beam exhibits a remarkable resistant to diffraction during the propagation [13

B. Lü, W. Huang, B. Zhang, F. Kong, and Q. Zhai, “Focusing properties of Bessel beams,” Opt. Commun. 131, 223–228 (1996). [CrossRef]

-15

J. Durnin, J. J. Miceli Jr, and J. H. Eberly, “Comparison of Bessel and Gaussian beams,” Opt. Lett. 13, 79–80 (1988). [CrossRef] [PubMed]

Comparisons between the Bessel and Gaussian beams for incident pulses have been reported both theoretically and experimentally from different aspects such as the beam divergence, power-transport efficiency [15

J. Durnin, J. J. Miceli Jr, and J. H. Eberly, “Comparison of Bessel and Gaussian beams,” Opt. Lett. 13, 79–80 (1988). [CrossRef] [PubMed]

], second harmonic generation efficiency [16

K. Shinozaki, C. Xu, H. Sasaki, and T. Kamijoh, “A comparison of optical second-harmonic generation efficiency using Bessel and Gaussian beams in bulk crystals,” Opt. Commun. 133, 300–304 (1997). [CrossRef]

], propagation diffracted by a circular aperture [17

P. L. Overfelt and C. S. Kenney, “Comparison of the propagation characteristics of Bessel, Bessel-Gauss, and Gaussian beams diffracted by a circular aperture,” J. Opt. Soc. Am. A 8, 732–745 (1991). [CrossRef]

], and plasma channels [18

P. Polynkin, M. Kolesik, A. Roberts, D. Faccio, P. D. Trapani, and J. Moloney, “Generation of extended plasma channels in air using femtosecond Bessel beams,” Opt. Express 16, 15733–15740 (2008). [CrossRef] [PubMed]

, 19

S. Akturk, B. Zhou, M. Franco, A. Couairon, and A. Mysyrowicz, “Generation of long plasma channels in air by focusing ultrashort laser pulses with an axicon,” Opt. Commun. 282, 129–134 (2009).

]. To our knowledge, however, a detailed study on the transverse-mode dependence of filamentation has not been reported in the literature.

The purpose of this paper is to theoretically investigate how the transverse-mode of the incident femtosecond pulse affects the propagation in the Kerr medium, in particular in terms of the formation of filament and plasma channel. We assume an Ar gas as a Kerr medium and numerically solve the extended nonlinear Schrödinger equation coupled with the electron density equation. Through the comparisons of the propagation dynamics for three different transverse-modes for incident pulses, i.e., Bessel, Gaussian, and Laguerre modes, we will show that the use of the Bessel mode significantly improves the propagation characteristics in the Kerr medium, for instance, in the formation of filament and plasma channel. Perhaps more importantly we will also show that the temporal distortion of the Bessel incident beam integrated over the beam radius is extremely small compared with the other modes. Since the pulse energy that the Bessel beam can carry is more than one order of magnitude larger than the other modes for the same peak intensity and pulse duration, these findings clearly indicate that the Bessel beam can be a very powerful tool in ultrafast nonlinear optics involving propagation in a Kerr medium.

2. Model

In this paper, we theoretically investigate the transverse-mode dependence of femtosecond filamentation in the Ar gas. Three different transverse modes are considered for incident pulses: Bessel, Gaussian, and Laguerre modes which have all axial symmetry with respect to the radial coordinate r. The extended nonlinear Schrödinger equation coupled with the electron density equation is used to describe the pulse propagation in the Ar gas.

2.1 Transverse modes of the incident beams

The transverse modes of incident pulses we employ in this paper are the Bessel

ε_{b} (r, t, 0) = ε_{0} J_{0} (α r) \exp (- t^{2} ⁄ t_{0}^{2})

(1)

Gaussian

ε_{b} (r, t, 0) = ε_{0} \exp ({- r}^{2} ⁄ w_{0}^{2} - t^{2} ⁄ t_{0}^{2})

(2)

and Laguerre modes

ε_{l} (r, t, 0) = ε_{0} (1 - 2 {(β r)}^{2}) \exp (- {(β r)}^{2} - t^{2} / t_{0}^{2})

(3)

where ε₀ is the peak amplitude of the incident electric field, t ₀ and w ₀ are the half temporal width and the spatial radius of the intensity (for 1/e²), respectively, and J ₀ is the zero order Bessel function. The parameters α and β are chosen in such a way that the full widths at half maximum of the beam diameter,

(w_{WFHM} = \sqrt{2 \ln 2} w_{0})

, of different transverse modes becomes equal. The spatial profiles of the three transverse modes are shown in Fig. 1. In practice, a Bessel beam can be obtained from a Gaussian beam through the mode-conversion with an axicon lens [10

P. Polesana, M. Franco, A. Couairon, D. Faccio, and P. Di Trapani, “Filamentation in Kerr media from pulsed Bessel beams,” Phys. Rev. A 77, 043814 (2008). [CrossRef]

, 18

, 19

S. Akturk, B. Zhou, M. Franco, A. Couairon, and A. Mysyrowicz, “Generation of long plasma channels in air by focusing ultrashort laser pulses with an axicon,” Opt. Commun. 282, 129–134 (2009).

], while a Laguerre beam can be obtained in a similar manner from a Gaussian beam with spiral phase plates [20

G.A. Turnbull, D.A. Robertson, G.M. Smith, L. Allen, and M.J. Padgett, “The generation of free-space Laguerre-Gaussian modes at millimetre-wave frequencies by use of a spiral phaseplate,” Opt. Commun. 127, 183–188 (1996). [CrossRef]

, 21

M. W. Beijersbergen, R. P .C. Coerwinkel, M. Kristensen, and J. P. Woerdman, “Helical-wavefront laser beams produced with a spiral phaseplate,” Opt. Commun. 112, 321–327 (1994). [CrossRef]

] or computer-generated holograms [22

N. R. Heckenberg, R. McDuff, C. P. Smith, H. Rubinsztein-Dunlop, and M. J. Wegener, “Laser beams with phase singularities,” Opt. Quantum Electron. 24, S951–S962 (1992). [CrossRef]

, 23

J. Arlt, K. Dholakia, L. Allen, and M. J. Padgett, “The production of multiringed Laguerre-Gaussian modes by computer-generated holograms,” J. Mod. Opt. 45, 1231–1237 (1998). [CrossRef]

Fig. 1. Spatial intensity distribution of the (a) Bessel, (b) Gaussian, and (c) Laguerre beams for the same beam diameter defined by the full width at the half maximum as shown by red arrows on the right panels.

2.2 Equations

In order to describe the pulse propagation in the Ar gas, we employ the extended nonlinear Schrödinger equation coupled with the electron density equation. The extended nonlinear Schrödinger equation for the electric field envelop ε(r, t, z) in a reference frame (t=t_lab -z/v_g ) moving at the group velocity v_g can be written as,

\partial ε ⁄ \partial z = i (\partial^{2} ε ⁄ \partial r^{2} + \partial ε ⁄ \partial r ⁄ r) 2 k_{0} n_{0} - i k " (\partial^{2} ε ⁄ \partial t^{2}) ⁄ 2

+ i k_{0} n_{2} {∣ ε ∣}^{2} ε - i k_{0} ρ ε ⁄ 2 n_{0} ρ_{c} - σ ρ ε ⁄ 2 - U_{i} W (I) (ρ_{n t} -) ε ⁄ 2 I

(4)

where the electron density ρ can be computed from the following equation,

\partial ρ ⁄ \partial t = W (I) (ρ_{n t} - ρ) + σ ρ {∣ ε ∣}^{2} ⁄ U_{i}

(5)

In Eqs. (4) and (5), k ₀=2π/λ ₀ and k″=20 fs²/m are the wave number in vacuum and the second order dispersive coefficient, respectively. n ₀ and n ₂=5×10^-19 cm²/W are the linear refractive index and the nonlinear refractive index, respectively. ρ_c =1.73×10²¹ cm^-3 is the critical plasma density above which the plasma becomes opaque, while ρ _nt=2.7×10¹⁹ cm^-3 represents the neutral atom density. σ=1×10^-19 cm² stands for the cross-section for inverse bremsstrahlung, and U_i =15.76 eV is the ionization potential. I=εε* denotes the laser intensity and W(I) is the photoionization rate. Now Eqs. (4) and (5) are solved simultaneously by the split-step method [24

G. P. Agrawal, Nonlinear Fiber Optics (Academic, San Diego, 1995).

, 25

M. Mlejnek, E.M. Wright, and J. V. Moloney, “Femtosecond pulse propagation in argon: A pressure dependence study,” Phys. Rev. E 58, 4903–4910 (1998) [CrossRef]

]; the diffraction part is calculated by the Crank-Nicholson method and the dispersion part is calculated in the frequency domain by Fourier transform, while the nonlinear part is solved in the time domain by the Runge-Kutta method.

2.3 Description of ionization

Filamentation takes place as a result of the balance between the diffraction, self-focusing due to the nonlinear index of refraction, and self-defocusing due to ionization. Clearly it is very important to choose the most appropriate model to describe the ionization processes. We note that different models for ionization may result in one or even a few orders of magnitude differences in the ionization rate at the same intensity, and naturally a filament with a wrong filament diameter and length will be numerically observed at a wrong intensity. This is so, since a filament originates from the dynamic balance between self-focusing and self-defocusing where the latter is directly related to the electron density in the medium. To describe the ionization process, it is often convenient to introduce the Keldysh parameter, γ, which is a measure of the adiabaticity [26

L.V. Keldysh, “Ionization in the field of a strong electromagnetic wave,” Sov. Phys. JETP 20, 1307–1314 (1965).

]. For an atom with the ionization potential U_i exposed to the linearly polarized laser field with the field amplitude E and the frequency ω ₀, the Keldysh parameter, γ, is defined as,

γ = ω_{0} \sqrt{2 m U_{i}} ⁄ e E

(6)

where e is the electron charge and m the electron mass.

In this paper we employ the PPT (Perelemov-Popov-Terent’ev) model [27

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

] instead of more commonly used tunneling or multiphoton ionization models, since the former is an improved version of the latters and valid for the wider intensity range.

2.3.1 PPT model

As mentioned above, the PPT model [27

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

] is a kind of improved version of the Keldysh theory [26

L.V. Keldysh, “Ionization in the field of a strong electromagnetic wave,” Sov. Phys. JETP 20, 1307–1314 (1965).

] where the ionization rate for any atom or ion with the orbital quantum number l and its projection m onto the quantization axis and the charge state Z is given by,

W (ω_{0}, E) = ω_{a . u .} \sqrt{6 ⁄ π} ∣ C_{n^{*}, l^{*}} ∣^{2} f (l, m) (U_{i} ⁄ 2 U_{H})

\times A_{m} (ω_{0}, γ) {(2 E_{0} ⁄ E \sqrt{1 + γ^{2}})}^{{2 n}^{*} - ∣ m ∣ - 3 ⁄ 2} \exp (- 2 E_{0} g (γ) ⁄ 3 E)

(7)

Naturally the PPT model can be reduced to the either tunneling or multiphoton ionization limits under certain conditions as we will briefly describe below. For detailed description, we refer to Refs. [26

L.V. Keldysh, “Ionization in the field of a strong electromagnetic wave,” Sov. Phys. JETP 20, 1307–1314 (1965).

-28

M. V. Ammosov, N. B. Delone, and V. P. Krainov, “Tunnel ionization of complex atoms and of atomic ions in an alternating electromagnetic field,” Sov. Phys. JETP 64, 1191–1194 (1986).

2.3.2 Tunneling and multiphoton ionization limits

Consider ionization of atoms in the initial state with l=m=0, n*=(U_i /U_H )^-1/2 and |C _n*,l*|²=2²ⁿ*/n*Γ(2n*). If the laser field is very strong, we have γ≪1 and tunneling ionization takes place. Under this condition we obtain A_m (ω ₀,γ)→1, and accordingly Eq. (7) can be reduced to

W (E) = ω_{a . u .} \sqrt{6 ⁄ π} ∣ C_{n^{*}, l^{*}} ∣^{2} (U_{i} ⁄ 2 U_{H}) {(2 E_{0} ⁄ E)}^{2 n^{*} - 3 ⁄ 2} \exp (- 2 E_{0} ⁄ 3 E)

(8)

which is the ionization rate in the tunneling ionization regime (ADK model) [28

M. V. Ammosov, N. B. Delone, and V. P. Krainov, “Tunnel ionization of complex atoms and of atomic ions in an alternating electromagnetic field,” Sov. Phys. JETP 64, 1191–1194 (1986).

]. From Eq. (8), we can see that the tunneling ionization rate is independent of the laser frequency. Similarly, if the laser field is sufficiently weak we have γ≫1 and multiphoton ionization takes place. Under this condition, Eq. (7) can be reduced to

W (ω_{0}, I) = σ_{K} I^{K}

(9)

where K=11 and σ_K =5.06×10^-140 s^-1 cm²²/W¹¹ for Ar at the central laser wavelength of 800 nm.

Fig. 2. Ionization probabilities of Ar calculated by the three different formulas, PPT model (solid), multiphoton ionization model (dashed), and ADK model (dotted), as a function of laser intensity. The duration and the central wavelength of the pulse are assumed to be 30 fs (FWHM) with a Gaussian temporal profile and 800 nm, respectively.

In Fig. 2 we compare the ionization probabilities of Ar for the 30 fs pulse (FWHM with a Gaussian temporal profile) at the central wavelength of 800 nm calculated by the PPT model, and the tunneling (i.e., the ADK model) and multiphoton ionization limits, respectively. Figure 2 suggests that the PPT model is most appropriate to describe ionization for the intensity range of 10¹³ W/cm²~10¹⁴ W/cm² which is the intensity range we assume in this paper. Note that γ is about equal to or smaller than unity in this intensity range.

3. Results and discussions

To start with, we solve Eqs. (4) and (5) in vacuum to see the effect of diffraction. Figure 3 shows the variation of the beam diameter for the three different transverse modes of incident pulses as a function of propagation distance, z, which is normalized with respect to the Rayleigh length, z ₀, where z ₀=πw ² ₀/λ₀ with λ ₀ being the central laser wavelength in vacuum. We can see that the pulse with the Bessel incident mode is free of diffraction over a long distance while the pulse with the Gaussian incident mode has the largest diffraction compared with the other two modes. For the Gaussian beam with w _FWHM=100 µm and λ ₀=800 nm, the Rayleigh length is z ₀=2.83 cm.

Fig. 3. Variation of the beam diameters in vacuum for the Bessel (solid), Gaussian (dashed), and Laguerre (dotted) incident beams for the pulse with a 30 fs duration, which is nothing but the diffraction.

Having checked the behavior of diffraction in vacuum, we now proceed to perform the propagation calculation in Ar at the atmospheric pressure with the incident pulse of 30 fs (FWHM) with a Gaussian temporal profile and a beam diameter of w _FWHM=100 µm for the three different transverse modes and the peak intensity of 3.2×10¹³ W/cm². Note that all the results shown below have been performed under these conditions except parameters mentioned specially, and the step sizes and the maximum beam radius we have employed to numerically solve Eqs. (4) and (5) are typically dt=0.5 fs, dr=5~10 µm, dz=0.01 cm and 2r _max=100w _FWHM.

It is very important to mention that each lobe of the Bessel beam (see Fig. 1(a)) contains approximately the same energy as the central part. This means that the pulse energy that the Bessel beam can carry is more than one order of magnitude larger than those of the Gaussian and Laguerre beams at the same peak intensity, pulse duration, and beam diameter. Indeed this is a very nice feature of the Bessel beam and we will make further remarks in terms of filamentation at the end of section 3.1.

3.1 Formation of filament

Now we perform propagation calculations in Ar gas at the atmospheric pressure and investigate the formation of filament for the three different transverse modes. Since the concept of filament is rather general as we have described in Sec. 1, we can always refer to the formation of filament for the incident pulse with any transverse mode. The length of the filament [29

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

] can be defined as a distance of the propagation where the pulse spatially keeps a narrow and nearly constant diameter with a clamping intensity [30

A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and G. Mourou, “Self-channeling of high-peak-power femtosecond laser pulses in air,” Opt. Lett. 20, 73–75 (1995). [CrossRef] [PubMed]

-32

A. Becker, N. Aközbek, K. Vijayalakshmi, E. Oral, C. M. Bowden, and S. L. Chin, “Intensity clamping and re-focusing of intense femtosecond laser pulses in nitrogen molecular gas,” Appl. Phys. B 73, 287–290 (2001).

] which is sufficiently high to induce ionization in the medium. An alternative way to define the filament length is to refer to the length of the plasma channel. In this paper we have adopted the former definition for the filament length.

Fig. 4. Variation of the beam diameter as a function of propagation distance in the Ar gas for the (a) Bessel, (b) Gaussian, and (c) Laguerre incident beams.

Fig. 5. Color-coded plot of the time-integrated propagating pulse as a function of propagation distance and beam radius for the (a) Bessel, (b) Gaussian, and (c) Laguerre incident beams.

Figures 4(a)-4(c) show the variation of the beam diameter, w _FWHM, during the propagation for the Bessel, Gaussian, and Laguerre incident beams, respectively. Due to the self-focusing the beam diameter becomes smaller by ~20 % for all transverse modes during the first few cm of propagation in Figs. 4(a)-4(c). Accordingly the ionization probability at r=0 significantly increases from 10^-8 (at z=0 cm) to 10^-2~10^-3 (at z=3 cm) as shown in Figs. 10(a)-10(c), which now implies that self-defocusing starts into play. From the fact that both self-focusing and self-defocusing play important roles it is obvious that we are certainly not in the linear propagation regime but in the nonlinear propagation regime. Hence we can interpret the propagation of the pulse with a nearly constant beam diameter (which is smaller than the initial beam diameter as we have mentioned above) in Fig. 4 as a formation of filament. From Fig. 4 it is clear that the length of the filament is the longest for the Bessel incident beam, while that of the Gaussian incident beam is the shortest. We can also see that, for the same peak intensity (we have repeated calculations at intensities other than 3.2×10¹³ W/cm² which are not shown here), the diameters of the filaments are nearly the same for all three modes. It is interesting to point out that Fig. 4(c) for the Laguerre incident beam exhibits two separated filaments located at 0<z<10 cm and z>12 cm. This can be understood as a result of refocusing, as already discussed for the filament by the Gaussian incident beam [32

-34

L. T. Vuong, R.B. Lopez-Martens, C. P. Hauri, and A.L. Gaeta, “Spectral reshaping and pulse compression via sequential filamentation in gases,” Opt. Express 16, 390–401 (2008). [CrossRef] [PubMed]

]. Indeed refocusing takes place at different intensities for different transverse modes. For instance, similar to the case of Laguerre incident beam in Fig. 4(c), two separated filaments for the Bessel incident beam can be observed at a lower intensity (1.5×10¹³ W/cm² at which the Gaussian beam cannot form a filament).

Figures 5(a)-5(c) show the color-coded plot of the variation of the time-integrated pulse, i.e., I (r,z)=∫^∞ ^-∞|ε(r,t,z)|² dt for the Bessel, Gaussian, and Laguerre incident beams. We noted that the intensity distribution of the Bessel incident beam at small r is more uniform than the other two modes. It is also clear that the length of the filament for the Gaussian incident beam is the shortest among the three transverse modes. We would like to point out that the results shown in Fig. 4 are good measures about the detailed propagation characteristics as shown in Fig. 5. Clearly, provided the same peak intensity, pulse duration, and beam diameter, the Bessel incident beam has a remarkable superiority over the Gaussian and Laguerre incident beams in terms of the persistence of filamentation for a long propagation distance. This means that, even if the pulse energy is very large we can obtain a long filament by using the Bessel incident beam, since the Bessel beam can carry much more pulse energy than the Gaussian beam at the same peak intensity, pulse duration, and beam diameter as mentioned before. In contrast, the use of the Gaussian incident beam results in multi-filament [35

Y.R. Shen, The Principles of Nonlinear Optics (Wiley-Interscience, NewYork, 1984).

-38

A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and G. Mourou, “Self-channeling of high-peak-power femtosecond laser pulses in air,” Opt. Lett. 20, 73–75 (1995). [CrossRef] [PubMed]

] if the pulse energy is too large.

It is also interesting to compare the propagation of the Bessel and Gaussian incident beams with the same pulse energy and pulse duration. To investigate this, we perform similar calculations for the Bessel incident beam at the same peak intensity (3.2×10¹³ W/cm²) and pulse duration (30 fs) but with a much smaller beam diameter, w _FWHM=10 µm, which means that the pulse energy of the w _FWHM=10 µm Bessel incident beam is comparable to that of the w _FWHM=100 µm Gaussian incident beam. We find (results not shown here) that the beam diameter of the w _FWHM=10 µm Bessel incident beam is nearly unchanged during the propagation while the ionization probability is negligible, i.e., <10^-5 at r=0, which should be compared with Fig. 10(a) for the w _FWHM=100 µm Bessel incident beam where the ionization probability is much larger, i.e., ~10^-2. Such a small ionization probability for the w _FWHM=10 µm Bessel incident beam implies that the self-defocusing does not play any role during the propagation, indicating that the nearly unchanged beam diameter for the w _FWHM=10 µm Bessel incident beam during propagation has nothing to do with the filament. It simply comes from the non-diffractive nature of the Bessel beam in the linear propagation regime. As for the w _FWHM=10 µm Gaussian incident beam it suffers from huge diffraction. Thus, provided the comparable pulse energy and pulse duration (and naturally with different beam diameters), the use of the Bessel incident beam again leads to the significant improvement in propagation characteristics compared with the Gaussian incident beam, although it may not be called a filament if the beam diameter is chosen to be too small (such as 10 µm).

3.2 Change of the maximum intensity

In Fig. 6, we show the maximum intensity of the pulse (I_M (z)=max(|ε(r,t,z)|²)) as a function of propagation distance for the three transverse modes. We can see that, at the beginning of propagation, i.e., z<3 cm where self-focusing is the dominant process during propagation, the behavior of I_M (z) is almost the same for the three transverse modes. However, at z>3 cm where self-defocusing also set in due to significant ionization by the self-focused beam and hence self-focusing and self-defocusing compete to each other, quite different behaviors are seen for the three transverse modes: The Bessel incident beam has the best ability to maintain the beam diameter at the higher peak intensity during the propagation, which is followed by the Laguerre incident beam. The most commonly used Gaussian incident beam shows the worst feature.

Fig. 6. Maximum intensity of the pulse as a function of propagation distance for the Bessel (solid), Gaussian (dashed), and Laguerre (dotted) incident beams.

3.3 Evolution of the temporal profile during the propagation

We further investigate the time-evolution of the three different modes during the propagation. As shown in Fig. 7(a) we have assumed that the temporal profile of the incident pulse is Gaussian and identical for the three transverse modes. In Figs. 7(b)-7(d) we show the temporal profile of the propagating pulse integrated over the beam radius, i.e., I(t, z)=∫^∞ ^-∞|ε(r,t,z)|² 2πrdr, at z=25 cm. We can see that the temporal profile, I(t, z=25 cm), for the Bessel incident beam is nearly undistorted, while those of the Gaussian and Laguerre incident beams exhibit significant distortions. The similar result at z=50 cm for the Bessel incident beam is shown in Fig. 7(e). Clearly the temporal profile still remains almost unchanged. Figures 8(a)-8(c) are the color-coded plot of the temporal profile, I(t, z), for the three different transverse modes as a function of propagation distance. Note that Figs. 7(a)-7(e) are nothing but the representative temporal profiles, I(t, z), at z=0, 25, and 50 cm in Fig. 8, respectively. As Fig. 8 clearly shows the Bessel incident beam is nearly undistorted during the propagation. This suggests that the use of the Bessel incident beam results in the significantly improved propagation characteristics in the Kerr medium.

Fig. 7. Temporal profile of the pulse integrated over the beam radius (a) at z=0 (note that this is exactly the same for the Bessel, Gaussian, and Laguerre incident beams), (b) at z=25 cm for the Bessel, (c) at z=25 cm for the Gaussian, (d) at z=25 cm for the Laguerre, and (e) at z=50 cm for the Bessel incident beams.

Fig. 8. Color-coded plot of the temporal profiles integrated over the beam radius as a function of propagation distance, for the (a) Bessel, (b) Gaussian, and (c) Laguerre incident beams.

Fig. 9. Temporal profiles of the pulse on axis (r=0) for the Bessel incident beam at different distances, (a) z=25 and (b) z=50 cm.

In Figs. 9(a) and 9(b) we show the temporal profile of the pulse on axis, i.e., I(r=0, t, z)=|ε(r=0,t, z)|² for the Bessel incident beam at z=25 and 50 cm, respectively. It is interesting to point out that, although the temporal profile of the Bessel beam integrated over the beam radius remains nearly unchanged during the propagation (see Figs. 7(a), 7(b), and 7(e)), the on-axis temporal profile itself is distorted during the propagation. In this sense the behavior of the Bessel incident beam is still similar to those of the Gaussian and Laguerre incident beams. The results shown in Figs. 7-9 for the Bessel incident beam seem to suggest that the redistribution of pulse energy for the Bessel incident beam takes place in space only, while that for the Gaussian and Laguerre incident beams takes place in both space and time. To our knowledge, the peculiarity of the energy transfer for the Bessel incident beam has never been reported before.

3.4 Plasma channel

Having shown representative results for the pulse from the different viewpoints, we now present some results for the medium, i.e., Ar gas. Plasma channel is a spatial area with a nearly uniform plasma density along the propagation direction in the medium after significant ionization. In Fig. 10 we show the ionization probability at r=0 produced by the three different transverse modes for incident pulses. Clearly the plasma channel produced by the Bessel incident beam is much longer than that by the Gaussian incident beam. For the Gaussian incident beam, the plasma channel is usually limited to the short length due to the severe distortion of the pulse during the propagation. Increasing the peak intensity or beam diameter does not help to extend the length of the plasma channel for the Gaussian incident beam, since it simply results in the formation of multi-filaments [35

Y.R. Shen, The Principles of Nonlinear Optics (Wiley-Interscience, NewYork, 1984).

-38

A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and G. Mourou, “Self-channeling of high-peak-power femtosecond laser pulses in air,” Opt. Lett. 20, 73–75 (1995). [CrossRef] [PubMed]

], indicating that the plasma channel is lost. In contrast the Bessel incident beam is far superior in terms of the length of the plasma channel. The dips appearing in Figs. 10(a) and 10(c) are due to the reduction of the intensity at r=0 originating from the temporal distortion of the pulse. Figures 11(a)-11(c) are the color-coded plot of the ionization probability as a function of propagation distance and beam radius for the Bessel, Gaussian and Laguerre incident beams, respectively. We can see that the plasma channel formed by the Bessel incident beam is the longest and most uniform among the three transverse modes. A comparison between Fig. 5 and Fig. 11 shows that the diameter of the plasma channel is smaller than that of the pulse. This is simply because the formation of the plasma channel requires highly nonlinear processes in terms of ionization (see Eqs. (5), (7)–(9)). That is, a small change of intensity can result in a large change in the plasma density.

Fig. 10. Ionization probability at r=0 as a function of propagation distance for the (a) Bessel, (b) Gaussian, and (c) Laguerre incident beams.

Fig. 11. Color-coded plot of the ionization probability as a function of propagation distance and beam radius for the (a) Bessel, (b) Gaussian, and (c) Laguerre incident beams.

3.5 Robustness of the superiority of the Bessel incident beam

So far we have shown numerical results for the 100 µm beam diameter and 3.2×10¹³ W/cm² peak intensity for the incident beam. Since the formation of filament and plasma channel is based on the rather delicate dynamic balance of self-focusing and self-defocusing, it is very important to show that the superiority of the Bessel incident beam to the others is robust and not accidental. For that purpose we have undertaken further calculations for different beam diameters and peak intensities. Some representative results are presented in Fig. 12 for different beam diameters and peak intensities for the incident pulse. For the Gaussian incident beam, the length of the filament may become longer if we employ a larger beam diameter and higher intensity (compare Fig. 12 and Fig. 4(b)). But it is still shorter than that for the Bessel incident beam at the same beam diameter and peak intensity. We would like to point out that the filament lengths for the Bessel incident beams are beyond the distance we have calculated, 50 cm, for all five cases in Figs. 4(a), and Figs. 12(a), 12(b), 12(e) and 12(f). Moreover, increasing the beam diameter and peak intensity for the Gaussian incident beam to have a longer filament does not work forever, since multi-filament will take place at a certain intensity. Having the results presented in Fig. 12, we can now claim, with much more confidence, that the propagation characteristics of the Bessel incident beam is remarkably nicer compared with those of the Gaussian and Laguerre incident beams in terms of the lengths of filament and plasma channel and its potential to hold much more (possibly more than 1-2 orders of magnitude) energy in the pulse. Another results (not shown here) similar to Fig. 7 also show that the temporal profile of the pulse integrated over the beam radius is nearly undistorted for the Bessel incident beam. This also confirms the robustness of our findings.

Fig. 12. Variation of the beam diameter as a function of propagation distance in the Ar gas for the incident beams with (a) 150 µm, 3.2×10¹³ W/cm², (b) 200 µm, 3.2×10¹³ W/cm², (e) 100 µm, 2.0×10¹³ W/cm², and (f) 100 µm, 4.0×10¹³ W/cm². Ionization probability at r=0 corresponding to (a), (b), (e), and (f) are shown in (c), (d), (g), and (h), respectively.

4. Conclusions

We have theoretically investigated the transverse-mode dependence of femtosecond pulse propagation in a Kerr medium (Ar gas) for the Bessel, Gaussian, and Laguerre incident beams by solving the extended nonlinear Schrodinger equation. Due to the different spatial intensity distributions and diffraction behaviors, we have found that the propagation characteristics represented by the filamentation and plasma channel are significantly different for the different transverse modes of the incident beams. Among others the Bessel incident beam turned out to have huge advantages over the others, since the length of the filament and also the plasma channel is much longer for the Bessel incident beam. Perhaps more importantly, we have also found that the temporal distortion of the pulse (integrated over the beam radius) is very small for the Bessel incident beam, while this is not the case for the Gaussian and Laguerre incident beams. A very similar trend has been found for different beam diameters and peak intensities, indicating that our findings are robust. This suggests that the femtosecond Bessel beam can be very useful for the high intensity signal transmission over the long distance. All these findings clearly indicate that the Bessel beam can be a very powerful tool in ultrafast nonlinear optics involving propagation in a Kerr medium.

Acknowledgments

This work was partially supported by a Grant-in-Aid for Scientific Research from the Ministry of Education and Science of Japan.

References and links

1.	G. Méchain, C. D’Amico, Y.-B. André, S. Tzortzakis, M. Franco, B. Prade, A. Mysyrowicz, A. Couairon, E. Salmon, and R. Sauerbrey, “Length of plasma filaments created in air by a multiterawatt femtosecond laser,” Opt. Commun. 247, 171–180 (2005). [CrossRef]
2.	P. Rairoux, H. Schillinger, S. Neirdeimer, M. Rodriguez, F. Ronneberger, R. Sauerbrey, B. Stein, D. Waite, C. Wedekind, H. Wille, L. Wöste, and C. Ziener, “Remote sensing of the atmosphere using ultrashort laser pulses,” Appl. Phys. B 71, 573–580 (2000). [CrossRef]
3.	G. Méchain, A. Mysyrowicz, M. Depiesse, and M. Pellet, “A virtual antenna produced in air by intense femtosecond laser pulses,” In: D.H. Titterton, S.M. Kirkpatrick, R. Stoian, R. Appleby, J.M. Chamberlain, and K.A. Krapels, (Eds.), Technologies for Optical Countermeasures II; Femtosecond Phenomena II; and Passive Millimetre-Wave and Terahertz Imaging II, vol. 5989. SPIE , p. 5989 (2005).
4.	J. Kasparian, R. Sauerbrey, D. Mondelain, S. Niedermeier, J. Yu, J. -P. Wolf, Y. -B. André, M. Franco, B. Prade, S. Tzortzakis, A. Mysyrowicz, M. Rodriguez, H. Wille, and L. Wöste, “Infrared extension of the supercontinuum generated by femtosecond terawatt laser pulses propagating in the atmosphere,” Opt. Lett. 25, 1397–1399 (2000). [CrossRef]
5.	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]
6.	S. Tzortzakis, G. Méchain, G. Patalano, Y.-B. André, B. Prade, M. Franco, A. Mysyrowicz, J.-M. Munier, M. Gheudin, G. Beaudin, and P. Encrenaz, “Coherent subterahertz radiation from femtosecond infrared filaments in air,” Opt. Lett. 27, 1944–1946 (2002). [CrossRef]
7.	A. Dubietis, E. Gaižauskas, G. Tamošauska,s, and P. Di. Trapani, “Light filaments without self-channeling,” Phys. Rev. Lett. 92, 253903 (2004). [CrossRef] [PubMed]
8.	J. Durnin, J. J. Miceli Jr, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987). [CrossRef] [PubMed]
9.	J. Durnin, “Exact solutions for nondiffracting beams. I. The scalar theory,” J. Opt. Soc. Am. A 4, 651–654 (1987). [CrossRef]
10.	P. Polesana, M. Franco, A. Couairon, D. Faccio, and P. Di Trapani, “Filamentation in Kerr media from pulsed Bessel beams,” Phys. Rev. A 77, 043814 (2008). [CrossRef]
11.	E. Gaižauskas, E. Vanagas, V. Jarutis, S. Juodkazis, V. Mizeikis, and H. Misawa, “Discrete damage traces from filamentation of Gauss-Bessel pulses,” Opt. Lett. 31, 80–82 (2006). [CrossRef] [PubMed]
12.	A. Dubietis, P. Polesana, G. Valiulis, A. Stabinis, P. Di Trapani, and A. Piskarskas, “Axial emission and spectral broadening in self-focusing of femtosecond Bessel beams,” Opt. Express 15, 4168–4175 (2007). [CrossRef] [PubMed]
13.	B. Lü, W. Huang, B. Zhang, F. Kong, and Q. Zhai, “Focusing properties of Bessel beams,” Opt. Commun. 131, 223–228 (1996). [CrossRef]
14.	S. Chávez-Cerda and G. H. C. New, “Evolution of focused Hankel waves and Bessel beams,” Opt. Commun. 181, 369–377 (2000). [CrossRef]
15.	J. Durnin, J. J. Miceli Jr, and J. H. Eberly, “Comparison of Bessel and Gaussian beams,” Opt. Lett. 13, 79–80 (1988). [CrossRef] [PubMed]
16.	K. Shinozaki, C. Xu, H. Sasaki, and T. Kamijoh, “A comparison of optical second-harmonic generation efficiency using Bessel and Gaussian beams in bulk crystals,” Opt. Commun. 133, 300–304 (1997). [CrossRef]
17.	P. L. Overfelt and C. S. Kenney, “Comparison of the propagation characteristics of Bessel, Bessel-Gauss, and Gaussian beams diffracted by a circular aperture,” J. Opt. Soc. Am. A 8, 732–745 (1991). [CrossRef]
18.	P. Polynkin, M. Kolesik, A. Roberts, D. Faccio, P. D. Trapani, and J. Moloney, “Generation of extended plasma channels in air using femtosecond Bessel beams,” Opt. Express 16, 15733–15740 (2008). [CrossRef] [PubMed]
19.	S. Akturk, B. Zhou, M. Franco, A. Couairon, and A. Mysyrowicz, “Generation of long plasma channels in air by focusing ultrashort laser pulses with an axicon,” Opt. Commun. 282, 129–134 (2009).
20.	G.A. Turnbull, D.A. Robertson, G.M. Smith, L. Allen, and M.J. Padgett, “The generation of free-space Laguerre-Gaussian modes at millimetre-wave frequencies by use of a spiral phaseplate,” Opt. Commun. 127, 183–188 (1996). [CrossRef]
21.	M. W. Beijersbergen, R. P .C. Coerwinkel, M. Kristensen, and J. P. Woerdman, “Helical-wavefront laser beams produced with a spiral phaseplate,” Opt. Commun. 112, 321–327 (1994). [CrossRef]
22.	N. R. Heckenberg, R. McDuff, C. P. Smith, H. Rubinsztein-Dunlop, and M. J. Wegener, “Laser beams with phase singularities,” Opt. Quantum Electron. 24, S951–S962 (1992). [CrossRef]
23.	J. Arlt, K. Dholakia, L. Allen, and M. J. Padgett, “The production of multiringed Laguerre-Gaussian modes by computer-generated holograms,” J. Mod. Opt. 45, 1231–1237 (1998). [CrossRef]
24.	G. P. Agrawal, Nonlinear Fiber Optics (Academic, San Diego, 1995).
25.	M. Mlejnek, E.M. Wright, and J. V. Moloney, “Femtosecond pulse propagation in argon: A pressure dependence study,” Phys. Rev. E 58, 4903–4910 (1998) [CrossRef]
26.	L.V. Keldysh, “Ionization in the field of a strong electromagnetic wave,” Sov. Phys. JETP 20, 1307–1314 (1965).
27.	A.M. Perelemov, V. S. Popov, and M. V. Terent’ev, “Ionization of atoms in an alternating electric field,” Sov. Phys. JETP 23, 924–933 (1966).
28.	M. V. Ammosov, N. B. Delone, and V. P. Krainov, “Tunnel ionization of complex atoms and of atomic ions in an alternating electromagnetic field,” Sov. Phys. JETP 64, 1191–1194 (1986).
29.	A. Couairon and A. Mysyrowicz, “Femtosecond filamentation in transparent media,” Phys. Rep. 441, 47–189 (2007) [CrossRef]
30.	A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and G. Mourou, “Self-channeling of high-peak-power femtosecond laser pulses in air,” Opt. Lett. 20, 73–75 (1995). [CrossRef] [PubMed]
31.	J. Kasparian, R. Sauerbrey, and S. L. Chin, “The critical laser intensity of self-guided light filaments in air,” Appl. Phys. B 71, 877–879 (2000).
32.	A. Becker, N. Aközbek, K. Vijayalakshmi, E. Oral, C. M. Bowden, and S. L. Chin, “Intensity clamping and re-focusing of intense femtosecond laser pulses in nitrogen molecular gas,” Appl. Phys. B 73, 287–290 (2001).
33.	M. Mlejnek, E. M. Wright, and J. V. Moloney, “Dynamic spatial replenishment of femtosecond pulses propagating in air,” Opt. Lett. 23, 382–384 (1998). [CrossRef]
34.	L. T. Vuong, R.B. Lopez-Martens, C. P. Hauri, and A.L. Gaeta, “Spectral reshaping and pulse compression via sequential filamentation in gases,” Opt. Express 16, 390–401 (2008). [CrossRef] [PubMed]
35.	Y.R. Shen, The Principles of Nonlinear Optics (Wiley-Interscience, NewYork, 1984).
36.	L. Bergé, S. Skupin, F. Lederer, G. Méjean, J. Yu, J. Kasparian, E. Salmon, J. P. Wolf, M. Rodriguez, L. Wöste, R. Bourayou, and R. Sauerbrey, “Multiple filamentation of terawatt laser pulses in air,” Phys. Rev. Lett. 92, 225002 (2004). [CrossRef] [PubMed]
37.	F. Vidal and T.W. Johnston, “Electromagnetic beam breakup: multiple filaments, single beam equilibria, and radiation,” Phys. Rev. Lett. 77, 1282–1285 (1996). [CrossRef] [PubMed]
38.	A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and G. Mourou, “Self-channeling of high-peak-power femtosecond laser pulses in air,” Opt. Lett. 20, 73–75 (1995). [CrossRef] [PubMed]

OCIS Codes
(320.2250) Ultrafast optics : Femtosecond phenomena
(320.7110) Ultrafast optics : Ultrafast nonlinear optics

ToC Category:
Ultrafast Optics

History
Original Manuscript: April 16, 2009
Revised Manuscript: June 10, 2009
Manuscript Accepted: June 11, 2009
Published: July 6, 2009

Citation
Zhenming Song, Zhigang Zhang, and Takashi Nakajima, "Transverse-mode dependence of femtosecond filamentation," Opt. Express 17, 12217-12229 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-15-12217

Sort: Author | Year | Journal | Reset

References

G. Méchain, C. D’Amico, Y.-B. André, S. Tzortzakis, M. Franco, B. Prade, A. Mysyrowicz, A. Couairon, E. Salmon, and R. Sauerbrey, "Length of plasma filaments created in air by a multiterawatt femtosecond laser," Opt. Commun. 247, 171-180 (2005). [CrossRef]
P. Rairoux, H. Schillinger, S. Neirdeimer, M. Rodriguez, F. Ronneberger, R. Sauerbrey, B. Stein, D. Waite, C. Wedekind, H. Wille, L. Wöste, and C. Ziener, "Remote sensing of the atmosphere using ultrashort laser pulses," Appl. Phys. B 71, 573-580 (2000). [CrossRef]
G. Méchain, A. Mysyrowicz, M. Depiesse, and M. Pellet, "A virtual antenna produced in air by intense femtosecond laser pulses," In: D.H. Titterton, S.M. Kirkpatrick, R. Stoian, R. Appleby, J.M. Chamberlain, K.A. Krapels, (Eds.), Technologies for Optical Countermeasures II; Femtosecond Phenomena II; and Passive Millimetre-Wave and Terahertz Imaging II, vol. 5989. SPIE, p. 5989 (2005).
J. Kasparian, R. Sauerbrey, D. Mondelain, S. Niedermeier, J. Yu, J. -P. Wolf, Y. -B. André, M. Franco, B. Prade, S. Tzortzakis, A. Mysyrowicz, M. Rodriguez, H. Wille, and L. Wöste, "Infrared extension of the supercontinuum generated by femtosecond terawatt laser pulses propagating in the atmosphere," Opt. Lett. 25, 1397-1399 (2000). [CrossRef]
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]
S. Tzortzakis, G. Méchain, G. Patalano, Y.-B. André, B. Prade, M. Franco, A. Mysyrowicz, J.-M. Munier, M. Gheudin, G. Beaudin, and P. Encrenaz, "Coherent subterahertz radiation from femtosecond infrared filaments in air," Opt. Lett. 27, 1944-1946 (2002). [CrossRef]
A. Dubietis, E. Gaižauskas, G. Tamošauskas and P. Di. Trapani, "Light filaments without self-channeling," Phys. Rev. Lett. 92, 253903 (2004). [CrossRef] [PubMed]
J. Durnin, J. J. MiceliJr, and J. H. Eberly, "Diffraction-free beams," Phys. Rev. Lett. 58, 1499-1501 (1987). [CrossRef] [PubMed]
J. Durnin, "Exact solutions for nondiffracting beams. I. The scalar theory," J. Opt. Soc. Am. A 4, 651-654 (1987). [CrossRef]
P. Polesana, M. Franco, A. Couairon, D. Faccio and P. Di Trapani, "Filamentation in Kerr media from pulsed Bessel beams," Phys. Rev. A 77, 043814 (2008). [CrossRef]
E. Gaižauskas, E. Vanagas, V. Jarutis, S. Juodkazis, V. Mizeikis, and H. Misawa, "Discrete damage traces from filamentation of Gauss-Bessel pulses," Opt. Lett. 31, 80-82 (2006). [CrossRef] [PubMed]
A. Dubietis, P. Polesana, G. Valiulis, A. Stabinis, P. Di Trapani, and A. Piskarskas, "Axial emission and spectral broadening in self-focusing of femtosecond Bessel beams," Opt. Express 15, 4168-4175 (2007). [CrossRef] [PubMed]
B. Lü, W. Huang, B. Zhang, F. Kong, and Q. Zhai, "Focusing properties of Bessel beams," Opt. Commun. 131, 223-228 (1996). [CrossRef]
S. Chávez-Cerda and G. H. C. New, "Evolution of focused Hankel waves and Bessel beams," Opt. Commun. 181, 369-377 (2000). [CrossRef]
J. Durnin, J. J. Miceli, Jr, and J. H. Eberly, "Comparison of Bessel and Gaussian beams," Opt. Lett. 13, 79-80 (1988). [CrossRef] [PubMed]
K. Shinozaki, C. Xu, H. Sasaki, and T. Kamijoh, "A comparison of optical second-harmonic generation efficiency using Bessel and Gaussian beams in bulk crystals," Opt. Commun. 133, 300-304 (1997). [CrossRef]
P. L. Overfelt and C. S. Kenney, "Comparison of the propagation characteristics of Bessel, Bessel-Gauss, and Gaussian beams diffracted by a circular aperture," J. Opt. Soc. Am. A 8, 732-745 (1991). [CrossRef]
P. Polynkin, M. Kolesik, A. Roberts, D. Faccio, P. D. Trapani, and J. Moloney, "Generation of extended plasma channels in air using femtosecond Bessel beams," Opt. Express 16, 15733-15740 (2008). [CrossRef] [PubMed]
S. Akturk, B. Zhou, M. Franco, A. Couairon, and A. Mysyrowicz, "Generation of long plasma channels in air by focusing ultrashort laser pulses with an axicon," Opt. Commun. 282, 129-134 (2009).
G.A. Turnbull, D.A. Robertson, G.M. Smith, L. Allen, and M.J. Padgett, "The generation of free-space Laguerre-Gaussian modes at millimetre-wave frequencies by use of a spiral phaseplate," Opt. Commun. 127,183-188 (1996). [CrossRef]
M. W. Beijersbergen, R. P.C. Coerwinkel, M. Kristensen, and J. P. Woerdman, "Helical-wavefront laser beams produced with a spiral phaseplate," Opt. Commun. 112,321-327 (1994). [CrossRef]
N. R. Heckenberg, R. McDuff, C. P. Smith, H. Rubinsztein-Dunlop, and M. J. Wegener, "Laser beams with phase singularities," Opt. Quantum Electron. 24,S951-S962 (1992). [CrossRef]
J. Arlt, K. Dholakia, L. Allen, and M. J. Padgett, "The production of multiringed Laguerre-Gaussian modes by computer-generated holograms," J. Mod. Opt. 45,1231-1237 (1998). [CrossRef]
G. P. Agrawal, Nonlinear Fiber Optics (Academic, San Diego, 1995).
M. Mlejnek, E.M. Wright, and J. V. Moloney, "Femtosecond pulse propagation in argon: A pressure dependence study," Phys. Rev. E 58, 4903-4910 (1998) [CrossRef]
L.V. Keldysh, "Ionization in the field of a strong electromagnetic wave," Sov. Phys. JETP 20, 1307-1314 (1965).
A.M. Perelemov, V. S. Popov, and M. V. Terent’ev, "Ionization of atoms in an alternating electric field," Sov. Phys. JETP 23, 924-933 (1966).
M. V. Ammosov, N. B. Delone, and V. P. Krainov, "Tunnel ionization of complex atoms and of atomic ions in an alternating electromagnetic field," Sov. Phys. JETP 64, 1191-1194 (1986).
A. Couairon and A. Mysyrowicz, "Femtosecond filamentation in transparent media," Phys. Rep. 441,47 - 189 (2007) [CrossRef]
A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and G. Mourou, "Self-channeling of high-peak-power femtosecond laser pulses in air," Opt. Lett. 20, 73-75 (1995). [CrossRef] [PubMed]
J. Kasparian, R. Sauerbrey, and S. L. Chin, "The critical laser intensity of self-guided light filaments in air," Appl. Phys. B 71, 877-879 (2000).
A. Becker, N. Aközbek, K. Vijayalakshmi, E. Oral, C. M. Bowden, and S. L. Chin, "Intensity clamping and re-focusing of intense femtosecond laser pulses in nitrogen molecular gas," Appl. Phys. B 73, 287-290 (2001).
M. Mlejnek, E. M. Wright, and J. V. Moloney, "Dynamic spatial replenishment of femtosecond pulses propagating in air," Opt. Lett. 23, 382-384 (1998). [CrossRef]
L. T. Vuong, R.B. Lopez-Martens, C. P. Hauri, and A.L. Gaeta, "Spectral reshaping and pulse compression via sequential filamentation in gases," Opt. Express 16, 390-401 (2008). [CrossRef] [PubMed]
Y.R. Shen, The Principles of Nonlinear Optics (Wiley-Interscience, NewYork, 1984).
L. Bergé, S. Skupin, F. Lederer, G. Méjean, J. Yu, J. Kasparian, E. Salmon, J. P. Wolf, M. Rodriguez, L. Wöste, R. Bourayou, and R. Sauerbrey, "Multiple filamentation of terawatt laser pulses in air," Phys. Rev. Lett. 92, 225002 (2004). [CrossRef] [PubMed]
F. Vidal and T.W. Johnston, "Electromagnetic beam breakup: multiple filaments, single beam equilibria, and radiation," Phys. Rev. Lett. 77, 1282-1285 (1996). [CrossRef] [PubMed]
A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and G. Mourou, "Self-channeling of high-peak-power femtosecond laser pulses in air," Opt. Lett. 20, 73-75 (1995). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper