Terahertz emission from a femtosecond laser focus in a two-color scheme

doi:10.1364/JOSAB.27.000016

« Show journal navigation

Terahertz emission from a femtosecond laser focus in a two-color scheme

Alexey V. Balakin, Alexander V. Borodin, Igor A. Kotelnikov, and Alexander P. Shkurinov »View Author Affiliations

JOSA B, Vol. 27, Issue 1, pp. 16-26 (2010)
http://dx.doi.org/10.1364/JOSAB.27.000016

View Full Text Article

Acrobat PDF (232 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (35)
Cited By
Figures (2)
Metrics

Abstract

We critically revise the theory of terahertz emission from a plasma filament induced in a gas media by one or two focusd femtosecond laser pulses. We distinguish a radiation pressure force (RPF) from a ponderomotive force (PF), discuss conditions for one of these forces to be the dominating contribution to the terahertz emission, and also show that the angular distribution of the emitted power critically depends on which of the two forces dominates in a particular experiment. We show that the experimentally observed periodic dependence of the emitted terahertz power on the gas pressure reveals the dominating role of the RPF over the PF, whereas the angular diagram of the emission allows us to determine the predominant direction of the force. We also emphasize that the terahertz emission originated by a transient photocurrent exhibits a different dependency from the phase difference between the first and the second harmonics of the optic laser field, which generally enables the experimental detection of the prevailing mechanism of the terahertz emission from the plasma filament.

1. INTRODUCTION

A particular subject of plasma science is to explain various kinds of electromagnetic radiation ranging from terahertz to x-ray frequencies induced by high-power laser pulses interacting with photoinduced plasma [1

P. Sprangle, E. Esarey, and A. Ting, “Nonlinear theory of intense laser-plasma interactions,” Phys. Rev. Lett. 64, 2011–2014 (1990). [CrossRef] [PubMed]

, 2

P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71, 1994–1997 (1993). [CrossRef] [PubMed]

]. It has been observed that a broadband electromagnetic pulse (EMP) of terahertz radiation can be emitted from a plasma filament produced by an ultrashort (less than 100 fs) high-intensity

(10^{12} - 10^{15} W / {cm}^{2})

laser pulse propagating in a gas medium [3

H. Hamster, A. Sullivan, S. Gordon, W. White, and R. W. Falcone, “Subpicosecond, electromagnetic pulses from intense laser-plasma interaction,” Phys. Rev. Lett. 71, 2725–2728 (1993). [CrossRef] [PubMed]

, 4

H. Hamster, A. Sullivan, S. Gordon, and R. W. Falcone, “Short-pulse terahertz radiation from high-intensity-laser-produced plasmas,” Phys. Rev. E 49, 671–677 (1994). [CrossRef]

, 5

P. Sprangle, J. R. Peñano, B. Hafizi, and C. A. Kapetanakos, “Ultrashort laser pulses and electromagnetic pulse generation in air and on dielectric surfaces,” Phys. Rev. E 69, 066415 (2004). [CrossRef]

]. Earlier published theories suggest that the basic mechanism behind the terahertz radiation relies on the effect of a ponderomotive force (PF), which separates electrons from ions spatially within the plasma filament [5

, 6

L. M. Gorbunov and A. A. Frolov, “Emission of low-frequency electromagnetic waves by a short laser pulse in stratified rarefied plasma,” J. Exp. Theor. Phys. 83, 967–973 (1996).

] and, consequently, generates a nonlinear electric current. Recently, it has been noted that the emitted terahertz power increases manifold if a two-color scheme is used where a fundamental optical wave with the frequency ω is mixed with its second harmonic (SH) satellite at the doubled frequency

2 ω

D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210–1212 (2000). [CrossRef]

, 8

T. Bartel, P. Gaal, K. Reimann, M. Woerner, and T. Elsaesser, “Generation of single-cycle THz transients with high electric-field amplitudes,” Opt. Lett. 30, 2805–2807 (2005). [CrossRef] [PubMed]

, 9

M. Kress, T. Löffler, M. D. Thomson, R. Dörner, H. Gimpel, K. Zrost, T. Ergler, R. Moshammer, U. Morgner, J. Ullrich, and H. G. Roskos, “Determination of the carrier-envelope phase of few-cycle laser pulses with terahertz-emission spectroscopy,” Opt. Lett. 29, 1120–1122 (2004). [CrossRef] [PubMed]

The nonlinear electric current gives rise to the plasma filament polarization. Phenomenological models, formulated in terms of nonlinear polarization susceptibilities, attribute the terahertz emission to a four-wave mixing (FWM) rectification process [7

D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210–1212 (2000). [CrossRef]

, 8

, 9

, 10

N. Bloembergen, “Recent progress in four-wave mixing spectroscopy,” in Laser Spectroscopy , H. Walther and K. W. Rothe, eds. (Springer, 1979), Vol. IV, pp. 340–348.

, 11

X. Xie, J. Dai, and X. -C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005 (2006). [CrossRef] [PubMed]

, 12

J. Reintjes, Nonlinear Optical Parametric Processes in Liquids and Gases , 2nd ed. (Academic, 1984).

]. The phenomenology describes polarization properties of the terahertz field reasonably well [9

, 13

A. Houard, Y. Liu, B. Prade, and A. Mysyrowicz, “Polarization analysis of terahertz radiation generated by four-wave mixing in air,” Opt. Lett. 33, 1195–1197 (2008). [CrossRef] [PubMed]

] and the fact that the terahertz field is subject to the phase difference φ between the fundamental and the SH waves. Specifically, it predicts that

E_{T H z} \propto E_{ω}^{2} E_{2 ω} cos φ

D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210–1212 (2000). [CrossRef]

, 11

X. Xie, J. Dai, and X. -C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005 (2006). [CrossRef] [PubMed]

, 14

J. Dai, X. Xie, and X. -C. Zhang, “Detection of broadband terahertz waves with a laser-induced plasma in gases,” Phys. Rev. Lett. 97, 103903 (2006). [CrossRef] [PubMed]

, 15

X. Lu, N. Karpowicz, and X. -C. Zhang, “Broadband terahertz detection with selected gases,” J. Opt. Soc. Am. B 26, A66–A73 (2009). [CrossRef]

] or

E_{T H z} \propto E_{ω}^{2} E_{2 ω} sin φ

]. Within the phenomenological approach it is hardly possible to distinguish whether the

cos φ

or the

sin φ

dependency is correct but the latter has been confirmed by experimental studies [9

, 16

K. Y. Kim, J. H. Glownia, A. J. Taylor, and G. Rodriguez, “Terahertz emission from ultrafast ionizing air in symmetry-broken laser fields,” Opt. Express 15, 4577–4584 (2007). [CrossRef] [PubMed]

]. Furthermore, it has been shown that

E_{T H z} \propto E_{2 ω}

and

E_{T H z} \propto E_{ω}^{2}

for sufficiently low fundamental wave energy (less than

100 μ J

) [8

, 13

A. Houard, Y. Liu, B. Prade, and A. Mysyrowicz, “Polarization analysis of terahertz radiation generated by four-wave mixing in air,” Opt. Lett. 33, 1195–1197 (2008). [CrossRef] [PubMed]

The understanding of plasma formation through optical breakdown has a paramount value for the unveiling of the physical mechanism of the terahertz emission from an optically generated plasma filament. Earlier dynamic theories of the emission assumed that the photoelectrons, pulled away from gas atoms by strong laser fields, are born with zero macroscopic (i.e., average) velocities [5

]. This natural contemplation was disputed in [16

, 17

]. It was suggested that the terahertz emission at a microscopic level originates from a transient photocurrent created by the photoelectrons produced in a symmetry-broken laser field composed of fundamental and phase-shifted SH optic waves.

In this paper we present a simple theory based on classical assumptions of the terahertz emission, which takes into account the FWM process. It has led us to the conclusion that the manifold increase of the radiated power, observed in a number of experiments, occurs mainly due to transient photocurrent effect and also due to enhanced plasma production in the process of photoionization, since the multiphoton ionization (MPI) at doubled frequency goes on faster than at the fundamental frequency. Additional enhancement occurs due to the effect of the plasma pressure force [18

C. -C. Cheng, E. M. Wright, and J. V. Moloney, “Generation of electromagnetic pulses from plasma channels induced by femtosecond light strings,” Phys. Rev. Lett. 87, 213001 (2001). [CrossRef] [PubMed]

]. For the first time, to our knowledge, we compute the radiation pressure force (RPF) for a laser pulse that contains both the first harmonic and the SH. We then argue that this newly computed force can produce a larger terahertz emission even if it is smaller than the PF. At a phenomenological level we also add to the theory the effect of the photocurrent by allowing a nonzero initial macroscopic velocity of the photoelectrons. Although we are not ready at the moment to compute the initial velocity from a strict quantum theory, we present a simple classical consideration of the photoionization process in a two-color laser field, which might hopefully elucidate major properties of the terahertz emission originated from the photocurrent. In particular, we note that the terahertz emission induced by the photocurrent scales as

sin φ

while that induced by the RPF or the PF varies in proportion to

cos φ

In our model the electrons and ions are initially created by the MPI within laser string, and a dipole moment is subsequently induced in the plasma filament via the RPF acting on the electrons and resulting from the velocity-dependent Lorentz force and from the PF resulting from spatial inhomogeneity of the laser fields. These forces separate the electrons from the heavy ions both longitudinally and transversely regarding the plasma filament axis on the short time scale of the laser pulse. Behind the laser pulse head the electron fluid oscillates at plasma frequency and generates an EMP propagating away from the plasma filament as first described in [18

]. Since plasma frequency is proportional to the square root of the local electron density, and the density varies from point to point within the plasma filament, the frequency spectrum of the terahertz emission ranges from a maximal value down to almost zero.

The paper is organized as follows. In Section 2 we compute the RPF and then in Section 3 we compute the PF. Slow motion of the plasma electrons under the action of these forces that spatially separates the electrons from the ions is analyzed in Section 4 and the resulting terahertz emission is computed in Section 5. Section 6 adds to our consideration of the effect of the transient photocurrent. Section 7 specifies experimental apparatus used to validate the theory and our conclusions are summarized in Section 8.

2. RADIATION PRESSURE FORCE

To calculate the RPF we extend a model first elaborated in [18

], where the force was computed for a single-color laser field. In contrast to the PF to be analyzed in Section 3, the RPF does not explicitly depend on the gradient sizes of the laser beam and, therefore, it dominates for a relatively wide beam. In a single-color scheme, the RPF is directed along the laser pulse path [18

], whereas a transverse component of the force emerges when a laser pulse at the fundamental frequency ω is mixed with that at the doubled frequency

2 ω

Following [3

, 4

H. Hamster, A. Sullivan, S. Gordon, and R. W. Falcone, “Short-pulse terahertz radiation from high-intensity-laser-produced plasmas,” Phys. Rev. E 49, 671–677 (1994). [CrossRef]

, 18

, 19

N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, “Optical second-harmonic generation in reflection from media with inversion symmetry,” Phys. Rev. 174, 813–822 (1968). [CrossRef]

], we employ an electron-fluid model of the plasma in the presence of the electric and magnetic fields E and B, which is governed by the equation

m \frac{\partial v}{\partial t} + m (v \cdot \nabla) v = e E + \frac{e}{c} v \times B - γ m v,

(1)

for the fluid velocity v of the electrons, where e and m are the charge and the mass of the electron, and

- γ m v

approximates a friction force due to electron scattering by heavy constituent particles (counting both ions and neutral atoms); throughout the paper, the heavy particles are assumed to remain unmovable.

The RPF indicated below as G formally originates from a velocity-dependent part of the Lorentz force, i.e., from the second term on the right-hand side (RHS) of Eq. (1). Since an oscillatory velocity of an electron in the field of a laser pulse is of the order of

v_{ω} \sim e E_{ω} / m ω

, and

B_{ω} \sim E_{ω}

, one formally finds that

G \sim m v_{ω}^{2} / ƛ

, where

ƛ = c / ω

is the reduced laser radiation wavelength.

The PF indicated below as F formally originates from the second term on the left-hand side of Eq. (1), which is evaluated as

m v_{ω}^{2} / a

, with a being a gradient length of the laser pulse envelope. Since

a ⪢ ƛ

, one might deduce that

G ⪢ F

in any circumstances, but this would be a rushed conclusion because the main term

m v_{ω}^{2} / ƛ

in G turns out to be zero due to phase cancellation. One can suggest although that the next term in an expansion of G over

v_{ω}

can be as large as F, alternatively, collisional friction can prevent complete phase cancellation.

At the first stage we ignore spatial derivatives in Eq. (1) and drop the term

m (v \cdot \nabla) v

; it will be retained over again in Section 3. Furthermore, spatial dimensions of the laser string are assumed in this section to be so enormous that the laser electromagnetic fields

E = E_{L}

and

B = B_{L}

can be approximated as spatially uniform over the transverse profile of the plasma filament.

In the rest of this section we sequentially consider a case of mutually parallel harmonics,

E_{2 ω} ∥ E_{ω}

, and that of perpendicular harmonics,

E_{2 ω} ⊥ E_{ω}

. After having completed the calculations, we shall see that the treatment of these two geometries is sufficient for a comprehensive characterization of the RPF.

2A. Case $E_{2 ω} ∥ E_{ω}$

Let a laser pulse be linearly polarized along the x axis while propagating in air along the positive direction of the z axis so that

E_{L} = E_{L} \hat{x}, B_{L} = E_{L} \hat{y},

(2)

where

E_{L} = \frac{1}{2} E_{ω} + \frac{1}{2} E_{2 ω} + c
                        .c
                        .

, and

E_{ω} = E_{1} (t - z / v_{g} (ω)) exp [i k (ω) z - i ω t],

E_{2 ω} = E_{2} (t - z / v_{g} (2 ω)) exp [i k (2 ω) z - i 2 ω t],

are the laser beam electric fields at the basic and doubled frequencies, with

t_{1}^{'} = t - z / v_{g} (ω)

and

t_{2}^{'} = t - z / v_{g} (2 ω)

being the retarded times in the reference frames moving with the laser pulse group velocities

v_{g} (ω)

and

v_{g} (2 ω)

, respectively.

The parameter

1 / c

in Eq. (1) formally tracks the order of perturbation theory and reflects the fact that the velocity-dependent Lorentz force is of the order of

v / c ⪡ 1

smaller than the dominant force due to the electric field under nonrelativistic conditions appropriate here. Then, to zeroth and first orders in the expansion parameter

1 / c

, Eq. (1) yields the triplet of equations,

m \frac{\partial v_{x}}{\partial t} = e E_{L} - e \frac{v_{z}}{c} E_{L} - γ m v_{x},

m \frac{\partial v_{y}}{\partial t} = - γ m v_{y},

m \frac{\partial v_{z}}{\partial t} = e \frac{v_{x}}{c} E_{L} - γ m v_{z} .

(3)

Under the assumption that the electric field envelope

E_{L}

varies slowly on the carrier time scale

ω^{- 1}

, the equation for

v_{x}

is approximately solved as

v_{x} \approx \frac{e}{2 m} [\frac{E_{ω}}{(γ - i ω)} + \frac{E_{ω}^{*}}{(γ + i ω)} + \frac{E_{2 ω}}{(γ - i 2 ω)} + \frac{E_{2 ω}^{*}}{(γ + i 2 ω)}] .

The equation for

v_{y}

means that

v_{y} = 0

, and the electron motion is limited to the

x z

plane.

Next, the inspection of the term

e (v_{x} / c) E_{L}

in the equation for

v_{z}

reveals that it has the low frequency part,

G_{z} = ⟨ e \frac{v_{x}}{c} E_{L} ⟩ = \frac{e^{2}}{2 m c} \frac{γ {| E_{ω} |}^{2}}{(γ^{2} + ω^{2})} + \frac{e^{2}}{2 m c} \frac{γ {| E_{2 ω} |}^{2}}{(γ^{2} + 4 ω^{2})},

(4)

representing the RPF along the direction of the laser pulse propagation.

Formula (4) has been derived in [18

]. It should not be confused with the light pressure force

f_{T} = σ_{T} ({| E_{ω} |}^{2} + {| E_{2 ω} |}^{2}) / 8 π

, where

σ_{T} = (8 π / 3) {(e^{2} / m c^{2})}^{2}

is the Thomson cross section of an electromagnetic wave scattering by a free electron. The light pressure force can be derived from Eq. (1) with the radiative reaction force

(2 e^{2} / 3 c^{3}) \ddot{v}

substituted instead of the friction force

- γ m v

. In most circumstances, the light pressure force is negligibly small,

f_{T} ⪡ G_{z}

Straightforward integration of Eq. (3) for

v_{z}

yields five harmonics, from zeroth to fourth, varying as

e^{\pm i 0 ω t}

e^{\pm i ω t}

e^{\pm i 2 ω t}

e^{\pm i 3 ω t}

, and

e^{\pm i 4 ω t}

, respectively,

\frac{v_{z}}{c} = {(\frac{e}{2 m c})}^{2} [\frac{E_{ω}}{(γ - i ω)} \frac{E_{ω}^{*}}{γ} + \frac{E_{ω}^{*}}{(γ + i ω)} \frac{E_{ω}}{γ} + \frac{E_{2 ω}}{(γ - i 2 ω)} \frac{E_{2 ω}^{*}}{γ} + \frac{E_{2 ω}^{*}}{(γ + i 2 ω)} \frac{E_{2 ω}}{γ}] + {(\frac{e}{2 m c})}^{2} [\frac{E_{ω}}{(γ - i ω)} \frac{E_{2 ω}^{*}}{(γ + i ω)} + \frac{E_{ω}^{*}}{(γ + i ω)} \frac{E_{2 ω}}{(γ - i ω)} + \frac{E_{2 ω}}{(γ - i 2 ω)} \frac{E_{ω}^{*}}{(γ - i ω)} + \frac{E_{2 ω}^{*}}{(γ + i 2 ω)} \frac{E_{ω}}{(γ + i ω)}] + {(\frac{e}{2 m c})}^{2} [\frac{E_{ω}}{(γ - i ω)} \frac{E_{ω}}{(γ - i 2 ω)} + \frac{E_{ω}^{*}}{(γ + i ω)} \frac{E_{ω}^{*}}{(γ + i 2 ω)}] + {(\frac{e}{2 m c})}^{2} [\frac{E_{ω}}{(γ - i ω)} \frac{E_{2 ω}}{(γ - i 3 ω)} + \frac{E_{ω}^{*}}{(γ + i ω)} \frac{E_{2 ω}^{*}}{(γ + i 3 ω)} + \frac{E_{2 ω}}{(γ - i 2 ω)} \frac{E_{ω}}{(γ - i 3 ω)} + \frac{E_{2 ω}^{*}}{(γ + i 2 ω)} \frac{E_{ω}^{*}}{(γ - i 3 ω)}] + {(\frac{e}{2 m c})}^{2} [\frac{E_{2 ω}}{(γ - i 2 ω)} \frac{E_{2 ω}}{(γ - i 4 ω)} + \frac{E_{2 ω}^{*}}{(γ + i 2 ω)} \frac{E_{2 ω}^{*}}{(γ + i 4 ω)}] .

The first pair of square brackets here comprises low (zero) frequency terms and the second pair contains the first harmonic, i.e., the terms varying as

e^{+ i ω t}

e^{- i ω t}

, etc. Putting the second and the third pairs into the term

e (v_{z} / c) E_{L}

of the equation for

v_{x}

gives the low frequency force,

G_{x} = ⟨ - e \frac{v_{z}}{c} E_{L} ⟩ = - \frac{e^{3}}{8 m^{2} c^{2}} \frac{3 γ^{2}}{(γ^{2} + ω^{2}) (γ^{2} + 4 ω^{2})} [E_{ω}^{2} E_{2 ω}^{*} + c
                              .c
                              .],

(5)

acting in the direction of the laser field polarization.

2B. Case $E_{2 ω} ⊥ E_{ω}$

Now let the SH be polarized perpendicularly to the first harmonic so that

E_{L} = \frac{1}{2} E_{ω} \hat{x} + \frac{1}{2} E_{2 ω} \hat{y} + c
                           .c
                           .,

B_{L} = - \frac{1}{2} E_{2 ω} \hat{x} + \frac{1}{2} E_{ω} \hat{y} + c
                           .c
                           .

(6)

The equations of motion now take the form

m \frac{\partial v_{x}}{\partial t} = \frac{e}{2} (E_{ω} + E_{ω}^{*}) - \frac{e}{2} \frac{v_{z}}{c} (E_{ω} + E_{ω}^{*}) - γ m v_{x},

m \frac{\partial v_{y}}{\partial t} = \frac{e}{2} (E_{2 ω} + E_{2 ω}^{*}) - \frac{e}{2} \frac{v_{z}}{c} (E_{2 ω} + E_{2 ω}^{*}) - γ m v_{y},

m \frac{\partial v_{z}}{\partial t} = \frac{e}{2} \frac{v_{x}}{c} (E_{ω} + E_{ω}^{*}) + \frac{e}{2} \frac{v_{y}}{c} (E_{2 ω} + E_{2 ω}^{*}) - γ m v_{z} .

(7)

Using the same method as in Subsection 2A for solving these equations recovers formula (4) for the longitudinal force

G_{z}

. The transverse force,

G_{y} = ⟨ - \frac{e}{2} \frac{v_{z}}{c} (E_{2 ω} + E_{2 ω}^{*}) ⟩ = - \frac{e^{3}}{8 m^{2} c^{2}} [\frac{E_{ω}^{2} E_{2 ω}^{*}}{(γ - i ω) (γ - i 2 ω)} + c
                              .c
                              .] \approx \frac{e^{3}}{16 m^{2} c^{2} ω^{2}} [E_{ω}^{2} E_{2 ω}^{*} + c
                              .c
                              .],

(8)

turns out to be parallel to

E_{2 ω}

. As shown in Subsection 2A, this is also the case for

E_{2 ω} ∥ E_{ω}

. However, in a general case the transverse force is not exactly parallel to

E_{2 ω}

since

G_{x} / G_{y} \neq E_{2 ω, x} / E_{2 ω, y}

. For

E_{2 ω, x} = E_{2 ω, y}

the ratio

G_{x} / G_{y} \approx 3 γ^{2} / 2 ω^{2}

is expected to be small. In other words, the effect of the SH is stronger when it is polarized differently from the first harmonic.

The longitudinal force

G_{z}

is small as compared with

G_{y}

provided that the SH is strong enough,

E_{2 ω, y} ⪢ 4 γ m c / e

. To evaluate γ we assume that the electron–atom collisions dominate over the electron–ion collisions, which is the case for sufficiently low ionization of the ambient gas in the plasma filament. Then

γ = n_{a} σ_{e a} v_{e}

, where

σ_{e a} \sim 10^{- 16} {cm}^{- 3}

is the cross section of the electron scattering by neutral atoms,

n_{a}

is the atoms' density, and the average velocity

v_{e}

of the electrons is of the order of their oscillatory velocity

v_{ω} \sim e E_{ω} / m ω

With these assumptions used, the condition

E_{2 ω} ⪢ 4 γ m c / e

is equivalent to the inequality

E_{2 ω} / E_{ω} ⪢ 4 c n_{a} σ_{e a} / ω \sim 0.2 n_{a} / n_{atm}

, where

n_{atm} \approx 3 \times 10^{19} {cm}^{- 3}

is the air density at normal conditions and

ω = 2 \times 10^{15} s^{- 1}

. The condition can be met in a rarefied gas. Furthermore, if it is not satisfied, the transverse force can still produce a more powerful radiation from the plasma filament as shown in Section 5.

Following [20

A. Proulx, A. Talebpour, S. Petit, and S. L. Chin, “Fast pulsed electric field created from the self-generated filament of a femtosecond Ti:sapphire laser pulse in air,” Opt. Commun. 174, 305–309 (2000). [CrossRef]

, 21

J. F. Ready, Effects of High Power Laser Radiation (Academic, 1971).

], many papers refer the estimation

γ \approx 10^{12} - 10^{13} s^{- 1}

to the collision rate in the air at atmospheric pressure. This is consistent with the estimation above for the electron energy in the range from 0.1 to 10 eV.

3. PONDEROMOTIVE FORCE

We can now retain the term

m (v \cdot \nabla) v

in Eq. (1) neglected earlier. Furthermore, we relax the assumption that the laser field is spatially uniform over the transverse profile of the plasma filament and retain spatial derivatives of the laser pulse envelope. Specifically, for a linearly polarized laser pulse with a given amplitude of the first harmonic,

E_{ω} = E_{1} (x, y, t - z / v_{g} (ω)) exp [i k (ω) z - i ω t],

exact expressions for the laser field have the form

E_{L} = \frac{1}{2} E_{ω} \hat{x} + c
                        .c
                        .,

B_{L} = \frac{1}{2} \frac{k c}{ω} E_{ω} \hat{y} - \frac{1}{2} \frac{i c}{ω} \nabla E_{ω} \times \hat{x} + c
                        .c
                        .

(9)

Similar formulas can be drawn for the SH. Approximate formulas (2, 6) are obtained from Eq. (9) by dropping

\nabla E_{ω}

and substituting

k c

with ω. The PF is identified as the sum of all terms with spatial derivatives of the pulse envelope, namely,

F = - \frac{1}{4} m (v_{ω}^{*} \cdot \nabla) v_{ω} + \frac{1}{4} \frac{e}{c} v_{ω}^{*} \times [- \frac{i c}{ω} \nabla E_{ω} \times \hat{x}] + c
                     .c
                     .

(10)

Straightforward calculations reveal a well known result [22

G. A. Askar’yan, “Cherenkov radiation from optical pulses,” Phys. Rev. Lett. 57, 2470 (1986). [CrossRef] [PubMed]

], which takes especially the simple form

F = - \frac{e^{2}}{4 m ω^{2}} \nabla {| E_{ω} |}^{2},

(11)

γ = 0

. In plasma physics formula (11) is sometimes called the Miller force. The standard derivation of the Miller force relies on the Lagrangian coordinates where an observer follows individual fluid particles as they move through space (see, e.g., [23

B. B. Kadomtsev, Collective Phenomena in Plasmas (Pergamon, 1982).

, 24

L. A. Artsimovich and R. Z. Sagdeev, Plasma Physics for Physicists (Atomizdat, 1979) (in Russian).

, 25

R. Fitzpatrick, The Physics of Plasmas (Lulu, 2008).

]). On the contrary, Eq. (10) gives the Miller force in the Eulerian coordinates that focus on specific locations in the space through which the fluid flows.

The PF can be neglected if

F_{z} ⪡ G_{z}

and

F_{y} ⪡ G_{y}

. The first condition is equivalent to the inequality

v_{ω}^{2} / l ⪡ e v_{ω} E_{ω} / m c

, where

l = c τ_{L}

is the coherence length of the laser pulse with the time duration

τ_{L}

. Putting here

v_{ω} \sim e E_{ω} / m ω

, one immediately concludes that the longitudinal PF

F_{z}

is small as compared with

G_{z}

l ⪢ c / γ

. The latter inequality is feasible even for short laser pulses.

The condition

F_{y} ⪡ G_{y}

leads to the inequality

c (E_{2 ω}^{2} / 8 π) π a^{2} ⪢ (c / 2 r_{e}) m c^{2},

(12)

where a is a characteristic radius of the plasma filament and

r_{e} = e^{2} / m c^{2}

is the classical electron radius. In other words, the laser power in the SH should exceed 4.4 GW for the PF to be smaller than the RPF.

The GW laser power is currently feasible, which provided the initial motivation for introducing into account the RPF and looking for a prospective physical mechanism that could create a dipole moment in the plasma filament [18

]. It is worth noting, however, that the transverse part

G_{⊥}

of the RPF can produce larger emission from the plasma filament than that of the PF

F_{⊥}

even if the ordering Eq. (12) fails. Since

G_{⊥}

has a predominant direction across the entire section of the plasma filament, it produces a dipole polarization and, hence, a dipole radiation. On the contrary,

F_{⊥}

is directed mainly toward the center of the filament. It produces therefore a dipole radiation only due to the unavoidable ellipticity of the filament cross section. Ceteris paribus, a quadrupole radiation is less powerful and it is therefore usually neglected in comparison with the dipole radiation caused by the longitudinal part of the PF

F_{z}

4. SLOW MOTION

The combined force

f = F + G

causes the plasma electrons to execute relatively slow motion that can be characterized by the electron–ion separation vector ξ . The latter is associated with the averaged velocity

⟨ v ⟩ = \partial ξ / \partial t

and obeys the equation

\frac{\partial^{2} ξ}{\partial t^{2}} + γ \frac{\partial ξ}{\partial t} = \frac{e}{m} E + \frac{f}{m} .

(13)

It can be derived from Eq. (1) by a straightforward averaging procedure and by subsequently dropping the nonlinear term

(\partial ξ / \partial t, \nabla) \partial ξ / \partial t

, which is justified by the smallness of ξ in comparison with a characteristic radius a of the plasma filament.

The polarization electric field E in Eq. (13) is related to the electron density perturbation

δ n_{e}

by the equation

div E = 4 π e δ n_{e} .

In turn,

δ n_{e}

obeys the linearized continuity equation

\frac{\partial δ n_{e}}{\partial t} + div (n_{e} \frac{\partial ξ}{\partial t}) = 0,

(14)

which yields

E = - 4 π e n_{e} ξ .

(15)

Putting Eq. (15) in Eq. (13) gives the final equation

\frac{\partial^{2} ξ}{\partial t^{2}} + γ \frac{\partial ξ}{\partial t} + ω_{p}^{2} ξ = \frac{f}{m},

(16)

for ξ , where

ω_{p} = \sqrt{4 π e^{2} n_{e} / m}

is the plasma frequency, and

n_{e}

is the local density of the electrons; in general, the density

n_{e} = n_{e} (r, t)

varies from point to point and also in time on a time scale different from the period

2 π / ω_{p}

of the plasma oscillations.

The driving force

f = f (r, t - z / v_{g})

is assumed below to be a function of the radius vector

r = (x, y, z)

and the time

t^{'} = t - z / v_{g}

in a comoving frame of reference. Within a plasma filament, a small difference between the group velocities of the fundamental harmonic,

v_{g} (ω) \approx c (1 - ω_{p}^{2} / 2 ω^{2})

, and that of the second one,

v_{g} (2 ω) \approx c (1 - ω_{p}^{2} / 8 ω^{2})

, can be neglected. Indeed, the mistiming length in the plasma,

Δ z_{p} = \frac{τ_{L}}{1 / v_{g} (2 ω) - 1 / v_{g} (ω)} \approx \frac{8 ω^{2}}{3 ω_{p}^{2}} c τ_{L},

(17)

for the harmonic envelopes

E_{1}

and

E_{2}

to become separated in space is approximately

10^{6}

times greater than the laser pulse head length

c τ_{L}

for a typical ordering

ω \sim 10^{3} ω_{p}

. As to the mistiming of the harmonics envelopes in the gas, it must be minimized by a proper experimental setup to a negligible amount since the effect of the SH cannot be observed otherwise. For this reason, we substitute

v_{g}

below with the speed of light in every occurrence of the retarded time,

t - z / v_{g} \approx t - z / c

On the contrary, in some circumstances the misphasing of the two harmonics should be taken into account. By misphasing we imply the phase factor

exp [i 2 k (ω) z - i k (2 ω) z]

that enters the term

E_{ω}^{2} E_{2 ω}^{*}

G_{⊥}

. Within the plasma filament,

Δ k_{p} = 2 k (ω) - k (2 ω) \approx - \frac{3}{4} \frac{ω_{p}^{2}}{c ω} .

(18)

In a neutral gas,

Δ k_{g} = \frac{2 ω}{c} [n_{ω} - n_{2 ω}]

(19)

is substantially smaller being proportional to the gas pressure through the refraction indices of the gas,

n_{ω}

and

n_{2 ω}

, at corresponding frequencies. Equation (18) is obtained from Eq. (19) by taking

n_{ω} = \sqrt{1 - ω_{p}^{2} / ω^{2}}

and

n_{2 ω} = \sqrt{1 - ω_{p}^{2} / {(2 ω)}^{2}}

Misphasing in the gas determines a relative phase difference φ between the first harmonic and the SH at the “entrance” to the plasma filament, whereas the misphasing in the plasma causes modulation of the driving force f within the filament. The modulation period,

ℓ_{p} = 2 π / | Δ k_{p} | = 8 π c ω / 3 ω_{p}^{2},

(20)

strongly depends on the density of the electron in the plasma filament. In what follows, the modulation is prescribed to the first argument r of the force

f (r, t - z / c)

. The plasma misphasing length

ℓ_{p}

substantially exceeds

c τ_{L}

but can be as short as a few centimeters.

Integrating Eq. (16) over time with the initial conditions

ξ = 0

and

\partial ξ / \partial t = 0

t = - \infty

and a constant electron density gives

ξ = \frac{1}{m \sqrt{ω_{p}^{2} - γ^{2} / 4}} \int_{0}^{\infty} exp [- γ τ / 2] sin [\sqrt{ω_{p}^{2} - γ^{2} / 4} τ] f (r, t - z / c - τ) d τ .

(21)

The electron density can be approximated by constant in time behind the ionization front of the laser pulse. We will relax this assumption in Section 6.

If force f varies slowly on the time scale

ω_{p}^{- 1}

, i.e.,

ω_{p} τ_{L} ⪢ π

, Eq. (21) reduces to

ξ \approx f (r, t - z / c) / m ω_{p}^{2},

within the laser pulse head,

| t - z / v_{g} | ≲ τ_{L}

, where the force f assumes a maximal magnitude. Behind the head,

t - z / v_{g} ⪢ τ_{L}

, the laser pulse excites a damping wake of standing electron plasma oscillations, and

ξ \approx - Im {f_{ω_{p}} exp [- (γ / 2 + i ω_{p}) (t - z / c)]} / 2 m ω_{p},

where

f_{ω_{p}}

is the time Fourier transform,

f_{Ω} (r) = \int_{- \infty}^{\infty} f (r, τ) e^{i Ω τ} d τ,

(22)

of the driving force at

Ω = ω_{p}

. The current density

j (r, t - z / c) = e n_{e} \partial ξ / \partial t

is calculated by differentiating Eq. (21) over time. Keeping only major terms yields

j \approx \frac{e n_{e}}{m} \int_{0}^{\infty} e^{- γ τ / 2} cos (ω_{p} τ) f (r, t - z / c - τ) d τ .

(23)

5. ELECTROMAGNETIC EMISSION FROM THE PLASMA FILAMENT

In a dipole approximation, the energy spectral intensity of the low frequency radiation generated by the slow motion of the plasma electrons per unit solid angle can be calculated directly from the well-known expression

\frac{d E}{d Ω d o} = \frac{Ω^{2}}{4 π c^{3}} {| n \times j_{Ω, K} |}^{2},

(24)

where n and

d o

are, respectively, a unit vector and an element of the solid angle in the direction of radiation,

K = n Ω / c

is the wave vector, Ω is the low frequency of the emitted radiation, and

j_{Ω, K} = \int d t d^{3} r e^{i Ω t - i K r} j (r, t)

is the space–time Fourier transform of the current density [see Eq. (66.9) in [26

L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields , Vol. II of Course of Theoretical Physics, 4th ed. (Butterworth & Heinemann, 1998).

]]. Note that Ω is assumed to be positive in Eq. (24), where it is already taken into account that

| j_{Ω, K} | = | j_{- Ω, - K} |

. Putting

f (r, t - z / c - τ) = \int \frac{d Ω}{2 π} f_{Ω} (r) e^{- i Ω (t - z / c - τ)}

in Eq. (23) and then integrating over τ, we obtain the time Fourier transform of the current density,

j_{Ω} = \frac{e}{2 m} n_{e} f_{Ω} (r) e^{i Ω z / c} [\frac{1}{γ / 2 + i ω_{p} - i Ω} + \frac{1}{γ / 2 - i ω_{p} - i Ω}] .

When taking the spatial Fourier transform of the current, we assume the radius of the plasma filament a to be smaller than the emitted wavelength

2 π c / Ω

so that

n_{e} \approx N_{e} δ (x) δ (y),

where

N_{e} = \int \int n_{e} d x d y

is the linear density of the electrons in the filament and δ stands for a delta function. In this approximation, only the driving force averaged over the filament cross section,

⟨ f ⟩ = \int \int d x d y n_{e} f / N_{e}, ⟨ f_{Ω} ⟩ = \int \int d x d y n_{e} f_{Ω} / N_{e},

determines radiated power and the Fourier transform of the current density,

j_{Ω, K} = \frac{e N_{e}}{2 m} \int ⟨ f_{Ω} ⟩ e^{i (Ω / c - K_{z}) z} d z [\frac{1}{γ / 2 + i ω_{p} - i Ω} + \frac{1}{γ / 2 - i ω_{p} - i Ω}],

(25)

does not depend on the transverse components

K_{x}, K_{y}

of the wave vector. The longitudinal component of the wave vector

K_{z}

is related to the angle θ between K and the axis z of the plasma filament by the equation

K_{z} = (Ω / c) cos θ

so that

Ω / c - K_{z} = 2 (Ω / c) {sin}^{2} θ / 2

To proceed further, we approximate the dependency of the force

f (r, t^{'})

on the retarded time

t^{'} = t - z / c

in the head of the laser pulse by a tablelike peak with the duration

τ_{L}

so that

f (r, t^{'}) = f (r) H (t^{'} + τ_{L} / 2) H (τ_{L} / 2 - t^{'})

, where H stands for the Heaviside step function. Then,

⟨ f_{Ω} ⟩ = ⟨ f (r) ⟩ τ_{L} sinc [Ω τ_{L} / 2],

(26)

where

sinc (ξ) = sin (ξ) / ξ

. We also take into account that the plasma filament has a finite length L, and some components of the driving force (namely,

G_{⊥}

) can be modulated due to the effect of the misphasing as discussed above. We simulate both these facts by writing

⟨ f (r) ⟩ = f_{0} H (z + L / 2) H (L / 2 - z) cos [φ + Δ k_{p} z],

(27)

where

f_{0}

is a constant vector, φ is the phase shift between the two harmonics due to optical path difference in the gas, and

Δ k_{p} z

is the phase modulation within the plasma filament. Both φ and

Δ k_{p}

should be set to zero for the PF F and the longitudinal part of the RPF

G_{z}

Performing the integration in Eq. (25) and putting the result in Eq. (24) gives an expression that is a bit cumbersome expression,

\frac{d E}{d Ω d o} = \frac{Ω^{2}}{4 π c^{3}} {(\frac{e N_{e} L τ_{L}}{2 m})}^{2} \frac{(γ^{2} / 4 + Ω^{2}) {sinc}^{2} (Ω τ_{L} / 2)}{(γ^{2} / 4 + {(ω_{p} - Ω)}^{2}) (γ^{2} / 4 + {(ω_{p} + Ω)}^{2})} {{[sinc (\frac{Ω L}{c} {sin}^{2} \frac{θ}{2} + \frac{Δ k_{p} L}{2}) + sinc (\frac{Ω L}{c} {sin}^{2} \frac{θ}{2} - \frac{Δ k_{p} L}{2})]}^{2} {cos}^{2} φ + {[sinc (\frac{Ω L}{c} {sin}^{2} \frac{θ}{2} + \frac{Δ k_{p} L}{2}) - sinc (\frac{Ω L}{c} {sin}^{2} \frac{θ}{2} - \frac{Δ k_{p} L}{2})]}^{2} {sin}^{2} φ} {| n \times f_{0} |}^{2} .

(28)

It simplifies to

\frac{d E}{d Ω d o} = \frac{Ω^{2}}{4 π c^{3}} {(\frac{e N_{e} L τ_{L}}{m})}^{2} \frac{(γ^{2} / 4 + Ω^{2}) {sinc}^{2} (Ω τ_{L} / 2)}{(γ^{2} / 4 + {(ω_{p} - Ω)}^{2}) (γ^{2} / 4 + {(ω_{p} + Ω)}^{2})} {sinc}^{2} (\frac{Ω L}{c} {sin}^{2} \frac{θ}{2}) {cos}^{2} φ {| n \times f_{0} |}^{2},

(29)

Δ k_{p} L / 2 ⪡ 1

, that is, if the plasma filament length is substantially smaller than the modulation length Eq. (20),

L ⪡ ℓ_{p}

This constitutes our principal result describing the spectral and angular distribution of the radiation. The frequency spectrum has a maximum at the plasma frequency, corresponding to

ω_{p} / 2 π \sim 1 THz

for

n_{e} = 10^{16} {cm}^{- 3}

. Experimental data exhibit a broadband frequency spectrum rather than a resonantlike one. It would mean either that the friction coefficient γ is comparable with

ω_{p}

or that the variation

ω_{p}

across the plasma filament cross section should be taken into account. This can be done by averaging

f / [γ / 2 \pm i ω_{p} - i Ω]

in Eq. (25) rather than just the force f alone. In our opinion, both these effects are equally important.

Although wide, the frequency spectrum in most experiments with femtosecond lasers is narrower than one could evaluate from the spectral width

π / τ_{L}

of the seed laser pulse. It means that the exact shape of the seed pulse is not very important, and the factor

{sinc}^{2} (Ω τ_{L} / 2)

can be dropped from Eqs. (28, 29) as it is approximately equal to 1 within the main part of the spectrum. We utilize this observation in Section 6.

As it is seen from Eqs. (28, 29), the emitted terahertz power is affected by the phase difference φ between the fundamental and the SH pulses caused by the dispersion of the refractive index along the path of these optical pulses in the gas to the foci, as reported by Cook and Hochstrasser [7

D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210–1212 (2000). [CrossRef]

] and later by Kress et al. [9

] and Chen et al. [27

Y. Chen, M. Yamaguchi, M. Wang, and X. -C. Zhang, “Terahertz pulse generation from noble gases,” Appl. Phys. Lett. 91, 251116 (2007). [CrossRef]

]. Furthermore, Eq. (28) also describes the effect of the dispersion inside the plasma filament, and the latter limits the effective length of the filament that contributes to the terahertz emission. In particular, Eq. (29) predicts that the total emitted power oscillates from zero to a maximal value as the phase shift φ varies with the gas pressure. The more accurate Eq. (28) excludes lowering the emitted power down to zero as the result of the phase modulation within the plasma filament. Instead of being zero, the total radiated power falls down to a minimal value, which depends on the ratio of the plasma filament length L to the modulation length

ℓ_{p}

. Since

φ = 0

both for the PF and longitudinal RPF, detecting the power oscillations would certainly indicate that the transverse RPF plays a dominating role in the emission of the terahertz waves from the plasma filament. It is remarkable that the maximal emitted power is proportional to the square of the total number

N_{e} L

of the free electrons in the plasma filament.

The angular distribution also provides crucial information regarding which of the two forces plays a dominating role. According to Eq. (29), the angular distribution contains multiple lobes, the angular positions of which are defined by the condition

2 L {sin}^{2} θ / 2 = N λ

for the emitted wavelength

λ = 2 π c / Ω

and integer number

N = 0, 1, 2, \dots

The lobe

N = 0

does not exist if the force

f_{0}

is directed along the plasma filament. This occurs when the longitudinal force is greater than the transverse one, i.e.,

f_{z} ⪢ f_{⊥}

. In this case, the angular distribution is proportional to

{sin}^{2} θ {sinc}^{2} (Ω L {sin}^{2} (θ / 2) / c)

, and the main lobe,

N = 0

, is empty inside resembling an inverse Cherenkov cone as shown in Fig. 1a for the case

L / λ = 50

. The lobe

N = 1

is the strongest and corresponds to the cone opening angle

θ \approx \sqrt{2 λ / L} .

(30)

On the contrary, if the transverse force is dominant, i.e.,

f_{⊥} ⪢ f_{z}

, the angular distribution is proportional to

(1 - {sin}^{2} α {sin}^{2} θ) {sinc}^{2} (Ω L {sin}^{2} (θ / 2) / c)

, where α is the azimuthal angle (around the filament axis z) counted from the direction

f_{⊥}

. In this case, the main lobe is filled. It is also slightly distorted in the azimuthal direction but this distortion is very weak since

(1 - {sin}^{2} α {sin}^{2} θ) \approx 1

for small θ. In addition, the contribution of the longitudinal force is suppressed by the small factor

{sin}^{2} θ \sim 2 λ / L

. It means that the domination of the transverse force over the longitudinal one requires a less severe condition

f_{⊥} ⪢ f_{z} (λ / L)

to be satisfied instead of the inequality

f_{⊥} ⪢ f_{z}

specified above from superficial consideration.

Comparing Figs. 1a, 1b shows that measuring the angular distribution of the radiated power allows determining the dominant direction of the driving force f. However, it is important to note that the contribution of the transverse PF

F_{⊥}

to the emitted power is substantially suppressed and might seem to be lower than that of the RPF

G_{⊥}

of the same magnitude. It follows from the fact that

\int \int n_{e} \nabla_{⊥} {| E_{ω} |}^{2} d x d y = 0

for an axisymmetric distribution of both

n_{e}

and

{| E_{ω} |}^{2}

The angular distribution becomes even more complicated if

ℓ_{p} \sim L

, when the exact formula (28) should be used instead of Eq. (29). In this case, not only the total emitted power varies with the gas pressure (through the phase shift φ) but also the angular distribution. A more certain conclusion that follows from Eq. (28) is that it predicts a drastic decrease in the emitted power when the plasma filament length greatly exceeds the modulation length,

L ⪢ ℓ_{p}

. In this case, the terms in the curly brackets drop as

{(ℓ_{p} / L)}^{2}

, and the total power occurs to be proportional to

{(N_{e} ℓ_{p})}^{2}

instead of

{(N_{e} L)}^{2}

. Thus,

ℓ_{p}

is an effective length of the plasma filament that effectively contributes to the radiation.

6. EFFECT OF THE TRANSIENT PHOTOCURRENT

We now relax the assumption that photoelectrons are born with zero average velocity. Although the MPI is essentially a quantum process, we employ a simple classical model in the hope that it catches such a crucial feature of the photoionization as the dependency on the phase shift φ between the first harmonic and the SH of the laser field.

Let us consider the motion of a free electron in the two-color laser field,

E (t) = E_{ω} cos (ω t) + E_{2 ω} cos (2 ω t + φ) .

(31)

We assume that the electron is born in the result of the MPI at the instant of time

t_{N}

and that its velocity at that moment is equal to zero. Then, we have

v (t) = v_{ω} sin (ω t) + α v_{ω} sin (2 ω t + φ) + v_{N},

(32)

for

t > t_{N}

, where

v_{ω} = e E_{ω} / m ω

and

α = E_{2 ω} / E_{ω}

. The last term in Eq. (32), which does not depend on time, is set by the condition

v (t_{N}) = 0

and turns out to be the average velocity of the electron in the laser field,

v_{N} = - v_{ω} sin (ω t_{N}) - α v_{ω} sin (2 ω t_{N} + φ) .

(33)

Employing the classical model of the photocurrent generation [16

, 17

], we make one more proposition, the validity of which we will discuss later. Namely, we assume that due to the extremely large nonlinearity of the MPI process tearing-off a bound electron from an atom occurs exclusively at the instants of time, corresponding to the maxima of the absolute magnitude of the electric field strength. Making use of the smallness of the parameter

α = E_{2 ω} / E_{ω}

in realistic experimental conditions, we find the phase of the fundamental harmonic at these instants of time,

ω t_{N} = π N - 2 α sin φ cos π N,

(34)

where N stands for an integer. Having calculated the average velocity of an electron,

v_{N} = \frac{3}{2} α v_{ω} sin φ,

(35)

born at these moments, we note that it does not depend on N and, hence, is the same both for odd and even maxima. However, these maxima have different magnitudes,

E_{N} \equiv | E (t_{N}) | = | E_{ω} | [1 + α cos π N] .

(36)

The numbers of the electrons produced at odd and even maxima are different.

Generally speaking, this difference might be very huge even for a very small ratio α since a photon with doubled frequency has doubled energy and, thus, a two times smaller number of such photons is required to complete the act of ionization. However, for the sake of simplicity we assume that the parameter α is so small that the photoionization is produced exclusively by the fundamental harmonic of the laser field. It is then reasonable to think that the number of the electrons produced in the odd maxima differs from that produced in the even ones by the value of the order of α. If it is not so and the ionization by the SH dominates over the fundamental harmonic, one needs to think that the electrons are predominantly born at the maxima of the SH but we will not consider this case. Since the number of the electrons produced at different moments differs by a small amount of the order of α according to our assumption, the difference can be neglected when calculating the total (macroscopic) initial momentum of the electrons because their average velocity Eq. (35) itself contains the small parameter α. We note also that any initial spread of the electron velocities does not appear in our model. It means that a hydrodynamic approximation with zero temperature can be used to describe further the motion of the electron fluid.

To include the effect of the initial electron velocity, we extend the equations of slow motion from Section 4. First, we need to add a source producing

\dot{n}

electrons per unit of time in a unit of volume and take into account that the electrons acquire an initial velocity

v_{0}

when they are born. Respectively, we add the term

\dot{n}

to the RHS of the continuity equation,

\frac{\partial n_{e}}{\partial t} + div (n_{e} v) = \dot{n},

(37)

and the summand

\dot{n} v_{0}

to the RHS of the equation of motion,

\frac{\partial}{\partial t} n_{e} v + \nabla \cdot (n_{e} v v) = \frac{e n_{e}}{m} E - γ n_{e} v + n_{e} \frac{f}{m} + \dot{n} v_{0} .

Recall that E stands here for a slow polarization field appearing due to the spatial separation of the electrons from the neutralizing ion background. Excluding

\partial n_{e} / \partial t

from the second equation with the aid of the first one, we obtain

\frac{\partial}{\partial t} v + (v \cdot \nabla) v = \frac{e}{m} E - γ v + \frac{f}{m} + \frac{\dot{n}}{n_{e}} [v_{0} - v] .

(38)

To proceed further, we linearize Eqs. (37, 38), assuming the velocity v and the perturbation of the electron density

δ n_{e}

to be equally small quantities. By an unperturbed electron density we shall imply the quantity

n = \int_{- \infty}^{t} \dot{n} (t) d t,

(39)

which is the formal solution of the continuity equation for

v \equiv 0

. In a focus of a femtosecond laser pulse, secondary ionization has no time to enter the scene and the source of the electrons

\dot{n}

can therefore be considered as a given function of time. Taking this fact into account when linearizing Eqs. (37, 38), putting there

n_{e} = n + δ n_{e}

, and expressing the velocity

v = \partial ξ / \partial t

through the vector ξ , we obtain

\frac{\partial δ n_{e}}{\partial t} + div (n_{e} \frac{\partial ξ}{\partial t}) = 0,

(40)

\frac{\partial^{2} ξ}{\partial t^{2}} + γ \frac{\partial ξ}{\partial t} = \frac{e}{m} E + \frac{f}{m} + \frac{\dot{n}}{n} [v_{0} - \frac{\partial ξ}{\partial t}] .

(41)

Since Eq. (40) exactly coincides with Eq. (14) from Section 4, expression (15) for the polarization electric field E also remains valid. Putting it in Eq. (41) gives the final equation

\frac{\partial^{2} ξ}{\partial t^{2}} + γ \frac{\partial ξ}{\partial t} + ω_{p}^{2} ξ = \frac{f}{m} + \frac{\dot{n}}{n} [v_{0} - \frac{\partial ξ}{\partial t}],

(42)

for ξ , where

ω_{p} = \sqrt{4 π e^{2} n / m}

is again the plasma frequency. This equation differs from Eq. (16) by the last term in the RHS. The driving force

f = f (r, t - z / v_{g})

again can be interpreted as a function of the radius vector

r = (x, y, z)

and the time

t^{'} = t - z / v_{g}

in a comoving frame of reference. The same assumption is also applicable to the quantities n and

v_{0}

Integrating Eq. (42) over time should be done with the initial conditions

ξ = 0

and

\partial ξ / \partial t = 0

t = - \infty

. Since the plasma density n, the plasma frequency

ω_{p}

, the driving force f, and the initial velocity

v_{0}

depend on time, the integration can be generally done only numerically. There is, however, a practically important case, when one can present an approximate solution.

As has been noted at the end of Section 4, the period

2 π / ω_{p}

of the plasma oscillations in present-day experiments with femtosecond lasers does not exceed the duration of the laser pulse. Restricting ourselves to the solution of Eq. (42) at larger times, we can say that the separation vector ξ in an arbitrary point with the coordinate z starts to differ from zero only after the laser pulse head passes through the point, i.e., at

t > z / v_{g}

. Behind the pulse head the unperturbed electron density remains constant as well as the plasma frequency. As a result, Eq. (42) is transformed into the ordinary differential equation

\frac{\partial^{2} ξ}{\partial t^{2}} + γ \frac{\partial ξ}{\partial t} + ω_{p}^{2} ξ = 0,

(43)

with zero RHS and constant coefficients, which depend on the point coordinates only as on the parameters. This equation should be solved with the initial conditions

ξ = 0,

(44)

\frac{\partial ξ}{\partial t} = \frac{1}{n} \int [n f / m + \dot{n} v_{0}] d t,

(45)

t = z / v_{g} + 0

. It is assumed here that n before the integral sign in the RHS of Eq. (45) is the electron density established after the laser pulse head passing, and the range of the integrations over the time comprises the entire interval of the pulse passing (e.g., from

t = z / v_{g} - τ_{L}

t = z / v_{g} + τ_{L}

). The solution of the reduced problem [Eqs. (43, 44, 45)] is trivial,

ξ = \frac{e^{- γ (t - z / v_{g}) / 2}}{\sqrt{ω_{p}^{2} - γ^{2} / 4}} sin [\sqrt{ω_{p}^{2} - \frac{γ^{2}}{4}} (t - \frac{z}{v_{g}})] {\frac{\partial ξ}{\partial t} |}_{t = z / v_{g} + 0} .

(46)

Elementary comparison of the summands in the integrand in Eq. (45) shows that the effect of the initial electron momentum acquired by the electrons at the act of the ionization causes a stronger plasma polarization if the momentum exceeds the momentum gained by the electrons from the driving force for the entire duration of the laser pulse, i.e.,

m v_{0} ⪢ f τ_{L}

. A formal estimation

v_{0} \sim α v_{ω}

, following from the classical treatment in the beginning of this section, leads to the conclusion that the condition

m v_{0} ⪢ f τ_{L}

is satisfied practically for any laser power feasible to the date. However, available quantum theories predict a substantially lower value of

v_{0}

, especially in the case of the laser pulse with the duration significantly exceeding the period of the laser field,

ω τ_{L} ⪢ 1

Indeed, the notion of the instant of time of the ionization has no evident sense in quantum physics. According to the uncertainty principle, the instant of time of the ionization is determined at most with the accuracy of the order of the period of the bound electron gyration around the atomic core

2 π / ω_{a}

(if the frequency

ω_{a} = J / ℏ

, corresponding to the ionization potential J can be interpreted so) or with the accuracy of the ionizing field period

2 π / ω

. As a quantum theory of the photoionization [28

V. S. Popov, “Tunnel and multiphoton ionization of atoms and ions in a strong laser field (Keldysh theory),” Phys. Usp. 47, 855–885 (2004). [CrossRef]

, 29

V. S. Popov, “Multiphoton ionization of atoms by an ultrashort laser pulse,” JETP Lett. 73, 1–5 (2001). [CrossRef]

] says, the ionization occurs at times close to the absolute maximum of the laser field (taking into account the field envelope time profile) in a so called limit of low frequency, when

v_{ω} ⪢ v_{e} = \sqrt{2 J / m}

[30

This condition can be rewritten as $c E_{ω}^{2} / 8 π ⪢ (α^{8} m c^{3} / 32 π r_{e}^{3}) {(J / J_{H})}^{3} / {(J / ℏ ω)}^{2}$ , where α is the fine structure constant, $r_{e}$ is the classical electron radius, $J / J_{H}$ is the ratio of the potential of ionization to that of the atom of hydrogen, and $J / ℏ ω$ is the number of photons required for the ionization. Numerically, this condition means that the laser power flux density should exceed $8.8 \times 10^{15} {(J / J_{H})}^{3} / {(J / ℏ ω)}^{2} W / {cm}^{2}$ .

]. According to the results of [29

V. S. Popov, “Multiphoton ionization of atoms by an ultrashort laser pulse,” JETP Lett. 73, 1–5 (2001). [CrossRef]

], the velocity distribution of the produced electrons has a maximum near

v_{max} = \int_{0}^{\infty} d t e E (t) / m

. Although the computation in [29

V. S. Popov, “Multiphoton ionization of atoms by an ultrashort laser pulse,” JETP Lett. 73, 1–5 (2001). [CrossRef]

] was performed for a special case, where

E (- t) = E (t)

, which corresponds to the phase shift φ, multiple of π, its result assumes that

v_{max}

drops quickly as the pulse duration increases since the integral

\int_{0}^{\infty} d t e E (t) / m

decreases in proportion to

1 / ω τ_{L}

for a pulse with a tablelike time envelope and even more quickly for smoother envelope profiles. For example, it gives

v_{max} \propto exp (- 1 / ω^{2} τ_{L}^{2})

for the Gaussian profile. For a monochromatic wave, the theory yields

v_{max} = 0

On the other hand, a numerical solution of the Schrödinger and experimental data, reported in [31

N. Karpowicz and X. -C. Zhang, “Coherent terahertz echo of tunnel ionization in gases,” Phys. Rev. Lett. 102, 093001 (2009). [CrossRef] [PubMed]

], supports the dependency of the average velocity on the phase shift proportional to

sin φ

, as predicted by our simple classical model. However, the same result contradicts to the results of [29

V. S. Popov, “Multiphoton ionization of atoms by an ultrashort laser pulse,” JETP Lett. 73, 1–5 (2001). [CrossRef]

] as the latter gives nonzero

v_{max}

for

φ = π N

7. EXPERIMENTAL RESULTS AND DISCUSSION

Experimental evidence is needed to establish the veracity of the model presented above. We used an experimental scheme analogous to that described in [7

D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210–1212 (2000). [CrossRef]

We employed an amplified Ti:sapphire laser system providing fundamental 120 fs and 800 nm pulses at the repetition rate of 1 kHz. The SH is generated in a beta barium borate crystal (BBO, type-I,

300 μ m

thick). Both the intense fundamental laser pulses and their SH descendants are focused into a gas cell by a lens with a focal length of

f = 18 cm

. Radiated terahertz waves are collected by a set of crystal quartz lenses and then focused by an off-axis parabola

(f = 25 mm)

onto

300 μ m

thick

⟨ 110 ⟩

oriented ZnTe crystal for electro-optical detection. The metal gas cell of 100 mm length has a window with a diameter of 25 mm. The front window is made of a

300 μ m

thick fused quartz plate and a rear window is fabricated from Teflon. The gas cell is filled with air or with a noble gas (argon, neon, xenon, or krypton) at adjustable pressure in the range from 0 to 1.5 bar. Pure gas can be pumped into or evacuated from the cell trough a gas diluting system. To separate and reduce the contribution of the terahertz wave from the BBO crystal, we used an adjustable special filter placed inside the gas cell in the focus area.

According to basic principles of the terahertz time domain measurements, the temporal profile of the emitted terahertz wave is measured by scanning the time delay of a probe beam at the fundamental frequency ω relative to the terahertz pulse. The phase shift φ between the fundamental and the SH pulses can be changed [7

D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210–1212 (2000). [CrossRef]

, 9

, 11

X. Xie, J. Dai, and X. -C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005 (2006). [CrossRef] [PubMed]

] by varying the distance

Δ l

between the SH generation BBO crystal and the laser focus. In addition, φ varies with the gas pressure in the cell since both

n_{ω} - 1

and

n_{2 ω} - 1

are proportional to the pressure but with different coefficients of proportionality. In both cases, the modulation of the terahertz wave amplitude has been detected in accordance with Eq. (28). Similar behavior has also been observed by Kress et al. [9

] who moved a SH nonlinear crystal (BBO) away from to the focus. In this case the phase shift can be written as

φ = (2 ω / c) (n_{ω} - n_{ω}) Δ l

The variation in the phase shift leads also to the changes in the temporal profile of the terahertz field and its polarity [11

X. Xie, J. Dai, and X. -C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005 (2006). [CrossRef] [PubMed]

, 16

] in agreement with our discussion at the end of Section 5. In our experiments, the phase shift changes by π as the gas (Xe) pressure changes by an amount in the range from 0.6 to 0.7 bar. Figure 2 illustrates this statement.

Using indices of refraction from [32

C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of krypton and xenon and their relation to those of helium and argon,” Proc. R. Soc. London, Ser. A 84, 2805–2807 (1910).

, 33

C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of neon,” Proc. R. Soc. London, Ser. A 83, 149–151 (1910). [CrossRef]

, 34

S. A. Korff and G. Breit, “Optical dispersion,” Rev. Mod. Phys. 4, 471–503 (1932). [CrossRef]

, 35

E. R. Peck and D. J. Fisher, “Dispersion of argon,” J. Opt. Soc. Am. A 54, 1362–1364 (1964). [CrossRef]

], we have estimated the pressure periods to be 1.1, 0.6, 0.25, and 10 atm for Ar, Kr, Xe, and Ne, respectively. These estimates agree fairly well with experimentally observed values, especially at lower pressures.

8. CONCLUSION

In this paper we discussed basic principles of the terahertz wave generation initiated by the optical breakdown in air of a noble gas. We focused mainly on the two-color scheme of the terahertz wave generation. We presented a simple classical model of the multiple photoionization in a two-color laser field that hopefully catches some peculiar features of the process. In particular, we emphasize that the contribution of the radiation pressure force (RPF) to the emitted terahertz power varies with the phase shift φ due to the optical phase difference in the air to the laser focus in proportion to

cos φ

, whereas the contribution of the transient photocurrent is proportional to

sin φ

We have distinguished the RPF from the ponderomotive force (PF) and, for the first time, we have calculated the transverse component of the RPF. We have shown that the effect of the second harmonic (SH) is stronger when it is polarized differently from the first harmonic. In addition, the effect of the transverse PF is suppressed as compared with the effect of the transverse RPF because the PF varies across the section of the laser beam. Furthermore, the effect of the longitudinal driving force is suppressed as compared with the transverse force by the factor

λ / L

Both empty and filled angular diagrams of the terahertz emission (see Fig. 1) are allowed by our theory and both types of the diagram were observed experimentally although no special care has been taken so far to quantitatively validate experimental data against our theory (as it was not available earlier). By inspecting the angular diagram of the terahertz radiation, one can deduce which component—the longitudinal or the transverse—makes a dominating contribution to the emitted terahertz power.

We showed that the effective length

ℓ_{p}

[given by Eq. (20)] of the emitting part of the plasma filament is related to the angle spread of the terahertz radiation angular diagram if

ℓ_{p} ⪡ L

. We explained that the effective length shortens due to the optical dispersion within the filament in inverse proportion to the electron density. We attributed the broadband character of the terahertz radiation to the spatial inhomogeneity of the electron density within the plasma filament.

If we compare our theoretical speculations with the complete experimental data available so far, we conclude that they do not allow validating our theory without doubts since some experimental results are mutually exclusive. We hope that future experiments will resolve this discrepancy.

ACKNOWLEDGMENTS

The authors are grateful to V. Davydenko for consulting us regarding properties of the laser-induced plasma, X.-C. Zhang and A. B. Savelev for useful discussions of various aspects of the laser-induced terahertz radiation, and M. Nazarov for the assistance in the experiments. This work was supported by the Russian Foundation for Basic Research (RFBR) in the frameworks of grants RFBR 08-02-00869 and 09-02-12198 and by the Russian Federal Agency for Science and Innovation (Rosnauka), the state contract 02.740.11.0223.

References and links

1.	P. Sprangle, E. Esarey, and A. Ting, “Nonlinear theory of intense laser-plasma interactions,” Phys. Rev. Lett. 64, 2011–2014 (1990). [CrossRef] [PubMed]
2.	P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71, 1994–1997 (1993). [CrossRef] [PubMed]
3.	H. Hamster, A. Sullivan, S. Gordon, W. White, and R. W. Falcone, “Subpicosecond, electromagnetic pulses from intense laser-plasma interaction,” Phys. Rev. Lett. 71, 2725–2728 (1993). [CrossRef] [PubMed]
4.	H. Hamster, A. Sullivan, S. Gordon, and R. W. Falcone, “Short-pulse terahertz radiation from high-intensity-laser-produced plasmas,” Phys. Rev. E 49, 671–677 (1994). [CrossRef]
5.	P. Sprangle, J. R. Peñano, B. Hafizi, and C. A. Kapetanakos, “Ultrashort laser pulses and electromagnetic pulse generation in air and on dielectric surfaces,” Phys. Rev. E 69, 066415 (2004). [CrossRef]
6.	L. M. Gorbunov and A. A. Frolov, “Emission of low-frequency electromagnetic waves by a short laser pulse in stratified rarefied plasma,” J. Exp. Theor. Phys. 83, 967–973 (1996).
7.	D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210–1212 (2000). [CrossRef]
8.	T. Bartel, P. Gaal, K. Reimann, M. Woerner, and T. Elsaesser, “Generation of single-cycle THz transients with high electric-field amplitudes,” Opt. Lett. 30, 2805–2807 (2005). [CrossRef] [PubMed]
9.	M. Kress, T. Löffler, M. D. Thomson, R. Dörner, H. Gimpel, K. Zrost, T. Ergler, R. Moshammer, U. Morgner, J. Ullrich, and H. G. Roskos, “Determination of the carrier-envelope phase of few-cycle laser pulses with terahertz-emission spectroscopy,” Opt. Lett. 29, 1120–1122 (2004). [CrossRef] [PubMed]
10.	N. Bloembergen, “Recent progress in four-wave mixing spectroscopy,” in Laser Spectroscopy , H. Walther and K. W. Rothe, eds. (Springer, 1979), Vol. IV, pp. 340–348.
11.	X. Xie, J. Dai, and X. -C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005 (2006). [CrossRef] [PubMed]
12.	J. Reintjes, Nonlinear Optical Parametric Processes in Liquids and Gases , 2nd ed. (Academic, 1984).
13.	A. Houard, Y. Liu, B. Prade, and A. Mysyrowicz, “Polarization analysis of terahertz radiation generated by four-wave mixing in air,” Opt. Lett. 33, 1195–1197 (2008). [CrossRef] [PubMed]
14.	J. Dai, X. Xie, and X. -C. Zhang, “Detection of broadband terahertz waves with a laser-induced plasma in gases,” Phys. Rev. Lett. 97, 103903 (2006). [CrossRef] [PubMed]
15.	X. Lu, N. Karpowicz, and X. -C. Zhang, “Broadband terahertz detection with selected gases,” J. Opt. Soc. Am. B 26, A66–A73 (2009). [CrossRef]
16.	K. Y. Kim, J. H. Glownia, A. J. Taylor, and G. Rodriguez, “Terahertz emission from ultrafast ionizing air in symmetry-broken laser fields,” Opt. Express 15, 4577–4584 (2007). [CrossRef] [PubMed]
17.	M. Kress, T. Löffler, M. D. Thomson, R. Dörner, H. Gimpel, K. Zrost, T. Ergler, R. Moshammer, U. Morgner, J. Ullrich, and H. G. Roskos, “Determination of the carrier-envelope phase of few-cycle laser pulses with terahertz-emission spectroscopy,” Nat. Phys. 2, 327–331 (2006). [CrossRef]
18.	C. -C. Cheng, E. M. Wright, and J. V. Moloney, “Generation of electromagnetic pulses from plasma channels induced by femtosecond light strings,” Phys. Rev. Lett. 87, 213001 (2001). [CrossRef] [PubMed]
19.	N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, “Optical second-harmonic generation in reflection from media with inversion symmetry,” Phys. Rev. 174, 813–822 (1968). [CrossRef]
20.	A. Proulx, A. Talebpour, S. Petit, and S. L. Chin, “Fast pulsed electric field created from the self-generated filament of a femtosecond Ti:sapphire laser pulse in air,” Opt. Commun. 174, 305–309 (2000). [CrossRef]
21.	J. F. Ready, Effects of High Power Laser Radiation (Academic, 1971).
22.	G. A. Askar’yan, “Cherenkov radiation from optical pulses,” Phys. Rev. Lett. 57, 2470 (1986). [CrossRef] [PubMed]
23.	B. B. Kadomtsev, Collective Phenomena in Plasmas (Pergamon, 1982).
24.	L. A. Artsimovich and R. Z. Sagdeev, Plasma Physics for Physicists (Atomizdat, 1979) (in Russian).
25.	R. Fitzpatrick, The Physics of Plasmas (Lulu, 2008).
26.	L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields , Vol. II of Course of Theoretical Physics, 4th ed. (Butterworth & Heinemann, 1998).
27.	Y. Chen, M. Yamaguchi, M. Wang, and X. -C. Zhang, “Terahertz pulse generation from noble gases,” Appl. Phys. Lett. 91, 251116 (2007). [CrossRef]
28.	V. S. Popov, “Tunnel and multiphoton ionization of atoms and ions in a strong laser field (Keldysh theory),” Phys. Usp. 47, 855–885 (2004). [CrossRef]
29.	V. S. Popov, “Multiphoton ionization of atoms by an ultrashort laser pulse,” JETP Lett. 73, 1–5 (2001). [CrossRef]
30.	This condition can be rewritten as $c E_{ω}^{2} / 8 π ⪢ (α^{8} m c^{3} / 32 π r_{e}^{3}) {(J / J_{H})}^{3} / {(J / ℏ ω)}^{2}$ , where α is the fine structure constant, $r_{e}$ is the classical electron radius, $J / J_{H}$ is the ratio of the potential of ionization to that of the atom of hydrogen, and $J / ℏ ω$ is the number of photons required for the ionization. Numerically, this condition means that the laser power flux density should exceed $8.8 \times 10^{15} {(J / J_{H})}^{3} / {(J / ℏ ω)}^{2} W / {cm}^{2}$ .
31.	N. Karpowicz and X. -C. Zhang, “Coherent terahertz echo of tunnel ionization in gases,” Phys. Rev. Lett. 102, 093001 (2009). [CrossRef] [PubMed]
32.	C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of krypton and xenon and their relation to those of helium and argon,” Proc. R. Soc. London, Ser. A 84, 2805–2807 (1910).
33.	C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of neon,” Proc. R. Soc. London, Ser. A 83, 149–151 (1910). [CrossRef]
34.	S. A. Korff and G. Breit, “Optical dispersion,” Rev. Mod. Phys. 4, 471–503 (1932). [CrossRef]
35.	E. R. Peck and D. J. Fisher, “Dispersion of argon,” J. Opt. Soc. Am. A 54, 1362–1364 (1964). [CrossRef]

Fig. 1 Angular distribution of the terahertz radiation computed by Eq. (29) for

L / λ = 50

: (a) driving force is directed along the plasma filament; (b) driving force is perpendicular to the plasma filament.

Fig. 2 Dependence of the terahertz field amplitude generated in the Xe versus the gas pressure. In the figure, open squares correspond to 980 mW and the dark triangles correspond to the 600 mW of the overage power of the ω fundamental field.

OCIS Codes
(190.2620) Nonlinear optics : Harmonic generation and mixing
(260.3090) Physical optics : Infrared, far
(300.2570) Spectroscopy : Four-wave mixing
(350.5400) Other areas of optics : Plasmas

ToC Category:
Nonlinear Optics

History
Original Manuscript: September 3, 2009
Manuscript Accepted: October 22, 2009
Published: December 11, 2009

Virtual Issues
January 8, 2010 Spotlight on Optics

Citation
Alexey V. Balakin, Alexander V. Borodin, Igor A. Kotelnikov, and Alexander P. Shkurinov, "Terahertz emission from a femtosecond laser focus in a two-color scheme," J. Opt. Soc. Am. B 27, 16-26 (2010)
http://www.opticsinfobase.org/josab/abstract.cfm?URI=josab-27-1-16

Sort: Author | Year | Journal | Reset

References

P. Sprangle, E. Esarey, and A. Ting, “Nonlinear theory of intense laser-plasma interactions,” Phys. Rev. Lett. 64, 2011-2014 (1990). [CrossRef] [PubMed]
P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71, 1994-1997 (1993). [CrossRef] [PubMed]
H. Hamster, A. Sullivan, S. Gordon, W. White, and R. W. Falcone, “Subpicosecond, electromagnetic pulses from intense laser-plasma interaction,” Phys. Rev. Lett. 71, 2725-2728 (1993). [CrossRef] [PubMed]
H. Hamster, A. Sullivan, S. Gordon, and R. W. Falcone, “Short-pulse terahertz radiation from high-intensity-laser-produced plasmas,” Phys. Rev. E 49, 671-677 (1994). [CrossRef]
P. Sprangle, J. R. Peñano, B. Hafizi, and C. A. Kapetanakos, “Ultrashort laser pulses and electromagnetic pulse generation in air and on dielectric surfaces,” Phys. Rev. E 69, 066415 (2004). [CrossRef]
L. M. Gorbunov and A. A. Frolov, “Emission of low-frequency electromagnetic waves by a short laser pulse in stratified rarefied plasma,” J. Exp. Theor. Phys. 83, 967-973 (1996).
D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210-1212 (2000). [CrossRef]
T. Bartel, P. Gaal, K. Reimann, M. Woerner, and T. Elsaesser, “Generation of single-cycle THz transients with high electric-field amplitudes,” Opt. Lett. 30, 2805-2807 (2005). [CrossRef] [PubMed]
M. Kress, T. Löffler, M. D. Thomson, R. Dörner, H. Gimpel, K. Zrost, T. Ergler, R. Moshammer, U. Morgner, J. Ullrich, and H. G. Roskos, “Determination of the carrier-envelope phase of few-cycle laser pulses with terahertz-emission spectroscopy,” Opt. Lett. 29, 1120-1122 (2004). [CrossRef] [PubMed]
N. Bloembergen, “Recent progress in four-wave mixing spectroscopy,” in Laser Spectroscopy, H.Walther and K.W.Rothe, eds. (Springer, 1979), Vol. IV, pp. 340-348.
X. Xie, J. Dai, and X.-C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005 (2006). [CrossRef] [PubMed]
J. Reintjes, Nonlinear Optical Parametric Processes in Liquids and Gases, 2nd ed. (Academic, 1984).
A. Houard, Y. Liu, B. Prade, and A. Mysyrowicz, “Polarization analysis of terahertz radiation generated by four-wave mixing in air,” Opt. Lett. 33, 1195-1197 (2008). [CrossRef] [PubMed]
J. Dai, X. Xie, and X.-C. Zhang, “Detection of broadband terahertz waves with a laser-induced plasma in gases,” Phys. Rev. Lett. 97, 103903 (2006). [CrossRef] [PubMed]
X. Lu, N. Karpowicz, and X.-C. Zhang, “Broadband terahertz detection with selected gases,” J. Opt. Soc. Am. B 26, A66-A73 (2009). [CrossRef]
K. Y. Kim, J. H. Glownia, A. J. Taylor, and G. Rodriguez, “Terahertz emission from ultrafast ionizing air in symmetry-broken laser fields,” Opt. Express 15, 4577-4584 (2007). [CrossRef] [PubMed]
M. Kress, T. Löffler, M. D. Thomson, R. Dörner, H. Gimpel, K. Zrost, T. Ergler, R. Moshammer, U. Morgner, J. Ullrich, and H. G. Roskos, “Determination of the carrier-envelope phase of few-cycle laser pulses with terahertz-emission spectroscopy,” Nat. Phys. 2, 327-331 (2006). [CrossRef]
C.-C. Cheng, E. M. Wright, and J. V. Moloney, “Generation of electromagnetic pulses from plasma channels induced by femtosecond light strings,” Phys. Rev. Lett. 87, 213001 (2001). [CrossRef] [PubMed]
N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, “Optical second-harmonic generation in reflection from media with inversion symmetry,” Phys. Rev. 174, 813-822 (1968). [CrossRef]
A. Proulx, A. Talebpour, S. Petit, and S. L. Chin, “Fast pulsed electric field created from the self-generated filament of a femtosecond Ti:sapphire laser pulse in air,” Opt. Commun. 174, 305-309 (2000). [CrossRef]
J. F. Ready, Effects of High Power Laser Radiation (Academic, 1971).
G. A. Askar'yan, “Cherenkov radiation from optical pulses,” Phys. Rev. Lett. 57, 2470 (1986). [CrossRef] [PubMed]
B. B. Kadomtsev, Collective Phenomena in Plasmas (Pergamon, 1982).
L. A. Artsimovich and R. Z. Sagdeev, Plasma Physics for Physicists (Atomizdat, 1979) (in Russian).
R. Fitzpatrick, The Physics of Plasmas (Lulu, 2008).
L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields, Vol. II of Course of Theoretical Physics, 4th ed. (Butterworth & Heinemann, 1998).
Y. Chen, M. Yamaguchi, M. Wang, and X.-C. Zhang, “Terahertz pulse generation from noble gases,” Appl. Phys. Lett. 91, 251116 (2007). [CrossRef]
V. S. Popov, “Tunnel and multiphoton ionization of atoms and ions in a strong laser field (Keldysh theory),” Phys. Usp. 47, 855-885 (2004). [CrossRef]
V. S. Popov, “Multiphoton ionization of atoms by an ultrashort laser pulse,” JETP Lett. 73, 1-5 (2001). [CrossRef]
This condition can be rewritten as cEω2/8π⪢(α8mc3/32πre3)(J/JH)3/(J/ℏω)2, where α is the fine structure constant, re is the classical electron radius, J/JH is the ratio of the potential of ionization to that of the atom of hydrogen, and J/ℏω is the number of photons required for the ionization. Numerically, this condition means that the laser power flux density should exceed 8.8×1015(J/JH)3/(J/ℏω)2 W/cm2.
N. Karpowicz and X.-C. Zhang, “Coherent terahertz echo of tunnel ionization in gases,” Phys. Rev. Lett. 102, 093001 (2009). [CrossRef] [PubMed]
C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of krypton and xenon and their relation to those of helium and argon,” Proc. R. Soc. London, Ser. A 84, 2805-2807 (1910).
C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of neon,” Proc. R. Soc. London, Ser. A 83, 149-151 (1910). [CrossRef]
S. A. Korff and G. Breit, “Optical dispersion,” Rev. Mod. Phys. 4, 471-503 (1932). [CrossRef]
E. R. Peck and D. J. Fisher, “Dispersion of argon,” J. Opt. Soc. Am. A 54, 1362-1364 (1964). [CrossRef]

Appl. Phys. Lett.

Y. Chen, M. Yamaguchi, M. Wang, and X.-C. Zhang, “Terahertz pulse generation from noble gases,” Appl. Phys. Lett. 91, 251116 (2007). [CrossRef]

J. Exp. Theor. Phys.

L. M. Gorbunov and A. A. Frolov, “Emission of low-frequency electromagnetic waves by a short laser pulse in stratified rarefied plasma,” J. Exp. Theor. Phys. 83, 967-973 (1996).

J. Opt. Soc. Am. A

E. R. Peck and D. J. Fisher, “Dispersion of argon,” J. Opt. Soc. Am. A 54, 1362-1364 (1964). [CrossRef]

J. Opt. Soc. Am. B

X. Lu, N. Karpowicz, and X.-C. Zhang, “Broadband terahertz detection with selected gases,” J. Opt. Soc. Am. B 26, A66-A73 (2009). [CrossRef]

JETP Lett.

V. S. Popov, “Multiphoton ionization of atoms by an ultrashort laser pulse,” JETP Lett. 73, 1-5 (2001). [CrossRef]

Nat. Phys.

M. Kress, T. Löffler, M. D. Thomson, R. Dörner, H. Gimpel, K. Zrost, T. Ergler, R. Moshammer, U. Morgner, J. Ullrich, and H. G. Roskos, “Determination of the carrier-envelope phase of few-cycle laser pulses with terahertz-emission spectroscopy,” Nat. Phys. 2, 327-331 (2006). [CrossRef]

Opt. Commun.

A. Proulx, A. Talebpour, S. Petit, and S. L. Chin, “Fast pulsed electric field created from the self-generated filament of a femtosecond Ti:sapphire laser pulse in air,” Opt. Commun. 174, 305-309 (2000). [CrossRef]

Opt. Express

K. Y. Kim, J. H. Glownia, A. J. Taylor, and G. Rodriguez, “Terahertz emission from ultrafast ionizing air in symmetry-broken laser fields,” Opt. Express 15, 4577-4584 (2007). [CrossRef] [PubMed]

Opt. Lett.

Phys. Rev.

N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, “Optical second-harmonic generation in reflection from media with inversion symmetry,” Phys. Rev. 174, 813-822 (1968). [CrossRef]

Phys. Rev. E

H. Hamster, A. Sullivan, S. Gordon, and R. W. Falcone, “Short-pulse terahertz radiation from high-intensity-laser-produced plasmas,” Phys. Rev. E 49, 671-677 (1994). [CrossRef]
P. Sprangle, J. R. Peñano, B. Hafizi, and C. A. Kapetanakos, “Ultrashort laser pulses and electromagnetic pulse generation in air and on dielectric surfaces,” Phys. Rev. E 69, 066415 (2004). [CrossRef]

Phys. Rev. Lett.

P. Sprangle, E. Esarey, and A. Ting, “Nonlinear theory of intense laser-plasma interactions,” Phys. Rev. Lett. 64, 2011-2014 (1990). [CrossRef] [PubMed]
P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71, 1994-1997 (1993). [CrossRef] [PubMed]
H. Hamster, A. Sullivan, S. Gordon, W. White, and R. W. Falcone, “Subpicosecond, electromagnetic pulses from intense laser-plasma interaction,” Phys. Rev. Lett. 71, 2725-2728 (1993). [CrossRef] [PubMed]
J. Dai, X. Xie, and X.-C. Zhang, “Detection of broadband terahertz waves with a laser-induced plasma in gases,” Phys. Rev. Lett. 97, 103903 (2006). [CrossRef] [PubMed]
X. Xie, J. Dai, and X.-C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005 (2006). [CrossRef] [PubMed]
C.-C. Cheng, E. M. Wright, and J. V. Moloney, “Generation of electromagnetic pulses from plasma channels induced by femtosecond light strings,” Phys. Rev. Lett. 87, 213001 (2001). [CrossRef] [PubMed]
N. Karpowicz and X.-C. Zhang, “Coherent terahertz echo of tunnel ionization in gases,” Phys. Rev. Lett. 102, 093001 (2009). [CrossRef] [PubMed]
G. A. Askar'yan, “Cherenkov radiation from optical pulses,” Phys. Rev. Lett. 57, 2470 (1986). [CrossRef] [PubMed]

Phys. Usp.

V. S. Popov, “Tunnel and multiphoton ionization of atoms and ions in a strong laser field (Keldysh theory),” Phys. Usp. 47, 855-885 (2004). [CrossRef]

Proc. R. Soc. London, Ser. A

C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of krypton and xenon and their relation to those of helium and argon,” Proc. R. Soc. London, Ser. A 84, 2805-2807 (1910).
C. Cuthbertson and M. Cuthbertson, “On the refraction and dispersion of neon,” Proc. R. Soc. London, Ser. A 83, 149-151 (1910). [CrossRef]

Rev. Mod. Phys.

S. A. Korff and G. Breit, “Optical dispersion,” Rev. Mod. Phys. 4, 471-503 (1932). [CrossRef]

Other

This condition can be rewritten as cEω2/8π⪢(α8mc3/32πre3)(J/JH)3/(J/ℏω)2, where α is the fine structure constant, re is the classical electron radius, J/JH is the ratio of the potential of ionization to that of the atom of hydrogen, and J/ℏω is the number of photons required for the ionization. Numerically, this condition means that the laser power flux density should exceed 8.8×1015(J/JH)3/(J/ℏω)2 W/cm2.
J. F. Ready, Effects of High Power Laser Radiation (Academic, 1971).
B. B. Kadomtsev, Collective Phenomena in Plasmas (Pergamon, 1982).
L. A. Artsimovich and R. Z. Sagdeev, Plasma Physics for Physicists (Atomizdat, 1979) (in Russian).
R. Fitzpatrick, The Physics of Plasmas (Lulu, 2008).
L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields, Vol. II of Course of Theoretical Physics, 4th ed. (Butterworth & Heinemann, 1998).
J. Reintjes, Nonlinear Optical Parametric Processes in Liquids and Gases, 2nd ed. (Academic, 1984).
N. Bloembergen, “Recent progress in four-wave mixing spectroscopy,” in Laser Spectroscopy, H.Walther and K.W.Rothe, eds. (Springer, 1979), Vol. IV, pp. 340-348.

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

Figures


Fig. 1	Fig. 2

« Previous Article | Next Article »

OSA is a member of CrossRef.

OpticsInfoBase

Journal of the Optical Society of America B

| OPTICAL PHYSICS

Terahertz emission from a femtosecond laser focus in a two-color scheme

Article Tools

Article Navigation

Abstract

1. INTRODUCTION

2. RADIATION PRESSURE FORCE

2A. Case $E_{2 ω} ∥ E_{ω}$

2B. Case $E_{2 ω} ⊥ E_{ω}$

3. PONDEROMOTIVE FORCE

4. SLOW MOTION

5. ELECTROMAGNETIC EMISSION FROM THE PLASMA FILAMENT

6. EFFECT OF THE TRANSIENT PHOTOCURRENT

7. EXPERIMENTAL RESULTS AND DISCUSSION

8. CONCLUSION

ACKNOWLEDGMENTS

References and links

References

Appl. Phys. Lett.

J. Exp. Theor. Phys.

J. Opt. Soc. Am. A

J. Opt. Soc. Am. B

JETP Lett.

Nat. Phys.

Opt. Commun.

Opt. Express

Opt. Lett.

Phys. Rev.

Phys. Rev. E

Phys. Rev. Lett.

Phys. Usp.

Proc. R. Soc. London, Ser. A

Rev. Mod. Phys.

Other

Cited By

Figures

Select a Conference
Session or Paper #
Year

OpticsInfoBase

Journal of the Optical Society of America B

| OPTICAL PHYSICS

Terahertz emission from a femtosecond laser focus in a two-color scheme

Article Tools

Article Navigation

Abstract

1. INTRODUCTION

2. RADIATION PRESSURE FORCE

2A. Case E 2ω∥ Eω

2B. Case E 2ω⊥ Eω

3. PONDEROMOTIVE FORCE

4. SLOW MOTION

5. ELECTROMAGNETIC EMISSION FROM THE PLASMA FILAMENT

6. EFFECT OF THE TRANSIENT PHOTOCURRENT

7. EXPERIMENTAL RESULTS AND DISCUSSION

8. CONCLUSION

ACKNOWLEDGMENTS

References and links

References

Appl. Phys. Lett.

J. Exp. Theor. Phys.

J. Opt. Soc. Am. A

J. Opt. Soc. Am. B

JETP Lett.

Nat. Phys.

Opt. Commun.

Opt. Express

Opt. Lett.

Phys. Rev.

Phys. Rev. E

Phys. Rev. Lett.

Phys. Usp.

Proc. R. Soc. London, Ser. A

Rev. Mod. Phys.

Other

Cited By

Figures

2A. Case $E_{2 ω} ∥ E_{ω}$

2B. Case $E_{2 ω} ⊥ E_{ω}$