Strong terahertz pulse generation by chirped laser pulses in tenuous gases

Mechanism of terahertz (THz) pulse generation in gases irradiated by ultrashort laser pulses is investigated theoretically. Quasi-static transverse currents produced by laser field ionization of gases and the longitudinal modulation in formed plasmas are responsible for the THz emission at the electron plasma frequency, as demonstrated by particlein-cell simulations including field ionization. The THz field amplitude scales linearly with the laser amplitude, which, however, holds only when the latter is at a moderate level. To overcome this limitation, we propose a scheme using chirped laser pulses irradiating on tenuous gas or plasma targets, which can generate THz pulses with amplitude 10-100 times larger than that from the well-known two-color laser scheme, enabling one to obtain THz field up to 10MV/cm with incident laser at ∼ 1016W/cm. © 2008 Optical Society of America OCIS codes: (300.6270) Spectroscopy, far infrared; (260.5210) Photoionization; (350.5400) Plasmas; (300.2140) Emission; References and links 1. M. S. Sherwin, C. A. Schmuttenmaer, and P. H. Bucksbaum, DOE-NSF-NIH Workshop on Opportunities in THz Science (http://www.er.doe.gov/bes/reports/files/THz rpt.pdf). 2. B. E. Cole, J. B. Williams, B. T. King, M. S. Sherwin, and C. R. Stanley, “Coherent manipulation of semiconductor quantum bits with terahertz radiation,” Nature (London) 410, 60-63 (2001). 3. Z.-M. Sheng, H.-C. Wu, K. Li, and J. Zhang, “Terahertz radiation from the vacuum-plasma interface driven by ultrashort intense laser pulses,” Phys. Rev. E 69, 025401(R)/1-4 (2004). 4. Z.-M. Sheng, K. Mima, J. Zhang, and H. Sanuki, “Emission of Electromagnetic Pulses from LaserWakefields through Linear Mode Conversion,” Phys. Rev. Lett. 94, 095003/1-4 (2005). 5. W. P. Leemans, C. G. R. Geddes, J. Faure, Cs. Tth, J. van Tilborg, C. B. Schroeder, E. Esarey, G. Fubiani, D. Auerbach, B. Marcelis, M. A. Carnahan, R. A. Kaindl, J. Byrd, and M. C. Martin, “Observation of terahertz emission from a laser-plasma accelerated electron bunch crossing a plasma-vacuum boundary,” Phys. Rev. Lett. 91, 074802/1-4 (2003). 6. D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett. 25, 1210-1212 (2000). (C) 2008 OSA 13 October 2008 / Vol. 16, No. 21 / OPTICS EXPRESS 16999 #99431 $15.00 USD Received 29 Jul 2008; revised 31 Aug 2008; accepted 1 Sep 2008; published 10 Oct 2008 7. M. Kress, T. Loffler, S. Eden, M. Thomson, and H. G. Roskos, “Terahertz-pulse generation by photoionization of air with laser pulses composed of both fundamental and second-harmonic waves,” Opt. Lett. 29, 1120-1122 (2004). 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). 9. X. Xie, J. Dai, and X.-C. Zhang, “Coherent control of THz wave generation in ambient air,” Phys. Rev. Lett. 96, 075005/1-4 (2006). 10. C. D’Amico, A. Houard, M. Franco, B. Prade, A. Mysyrowicz, A. Couairon, and V. T. Tikhonchuk, “Conical Forward THz Emission from Femtosecond-Laser-Beam Filamentation in Air,” Phys. Rev. Lett. 98, 235002/1-4 (2007). 11. Y. Liu, A. Houard, B. Prade, S. Akturk, A. Mysyrowicz, and V. T. Tikhonchuk, “Terahertz Radiation Source in Air Based on Bifilamentation of Femtosecond Laser Pulses,” Phys. Rev. Lett. 99, 135002/1-4 (2007). 12. 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). 13. 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). 14. H. C. Wu, J. Meyer-ter-Vehn, and Z. M. Sheng, “Phase-sensitive terahertz emission from gas targets irradiated by few-cycle laser pulses,” New J. Phys. 10, 043001/1-10 (2008). 15. 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). 16. B. M. Penetrante and J. N. Bardsley, “Residual energy in plasmas produced by intense subpicosecond lasers,” Phys. Rev. A 43, 3100-3113 (1991). 17. M. Chen, Z.-M. Sheng, J. Zheng, Y.-Y. Ma, and J. Zhang, “Developments and applications of multi-dimensional particle-in-cell codes in the investigation of laser plasma interactions,” Chin. J. Comput. Phys. 25, 43-50 (2008). 18. R. Nuter, S. Skupin, and L. Berge, “Chirp-induced dynamics of femtosecond filaments in air,” Opt. Lett. 30, 917-919 (2005). 19. G.-Z. Sun, E. Ott, Y. C. Lee, and P. Guzdar, “Self-focusing of short intense pulses in plasmas,” Phys. Fluids 30, 526-532 (1987). 20. W.-M. Wang and C.-Y. Zheng, “The evolution of the transverse centroid of asymmetric laser field in plasmas with various density distributions,” Phys. Plasmas 13, 053112/1-7 (2006).


Introduction
Intense THz pulses with the field strength up to MV/cm or beyond are expected to be important for nonperturbative THz electro-optics, nonlinear THz spectroscopies, THz extreme nonlinear physics in condensed matters, semiconductors, and nanostructured materials, etc. [1,2], which are mostly unexplored.Such THz pulses are usually obtained from accelerator-based sources, which are still limited by the bandwidth, peak intensity, waveform, as well as the availability to most users.Therefore table-top intense THz sources are attracting significant attention recently.For example, it is found that strong THz radiation can be produced from the laser wakefield in inhomogeneous plasmas by linear mode conversion [3,4] or from the transition radiation at plasma-vacuum boundaries using ultrashort electron bunches [5].THz emission can also be produced from tenuous gas targets irradiated with ultrashort laser pulses, which can be realized even using moderately intense laser pulses [6,7,8,9,10,11].Actually Hamster et al. observed electromagnetic (EM) emission with the same frequency as the electron plasma frequency ω p (around the THz frequency for a plasma electron density n e = m e ω 2 p /4πe 2 = 1.2 × 10 16 cm −3 ) when irradiating helium gas with ultrashort laser pulses in 1993 [12].In 2000, Cook et al. proposed and demonstrated that THz pulses can be generated efficiently when two laser pulses (with one at fundamental and another at the second-harmonic) co-propagate in air [6].Several authors attributed this THz pulse generation to four-wave rectification [6,7,8,9], and among of them Kress et al. found plasma formation is necessary for this THz emission [7].D'Amico et al. proposed the transition-Cherenkov emission from the plasma space charge moving behind the ionization front at light velocity [10,11].Recently, Kim et al. [13] and Wu et al. [14] have suggested that transverse drift currents excited by symmetry-broken laser field ionization of gases are responsible for this THz emission.However, their photocurrent model cannot explain the THz emission frequencies and holds only at moderate laser intensity.Whereas to get multi-MV/cm THz pulses, one has to invoke strong laser pulses.Under high laser intensity, the ionization can be completed within 1-2 laser cycles and in this case the dependence of THz pulse intensity on laser pulse intensity is still unclear.
In this paper we present an analytical model describing the mechanism of the THz emission in gases irradiated by high intensity lasers and propose a new scheme generating strong THz pulses by use of chirped laser pulses.It is found that the quasi-static transverse drift currents formed during the field ionization process and sequently modulated by formed plasmas in the longitudinal direction are responsible for the THz-wave generation with the frequency ω p .The THz wave amplitude is proportional to the drift currents dominated by uncontrollable field ionization of gases for either one-color or two-color lasers.Therefore, the THz wave amplitude is not strengthened monotonically with growing laser amplitude for intense lasers.If one takes chirped laser pulses, however, drift currents can be excited efficiently due to the asymmetric field in each laser cycle and the obtained THz wave amplitude can scale linearly with the incident laser amplitude.

Theoretical model
Interaction of ultrashort laser pulses with tenuous gases includes laser field ionization of gases and the interaction of the laser pulses with formed plasmas.The laser field ionization rates can be calculated by the ADK formula [15,16].For free electrons, they follow the equation of motion given by where v ⊥ is transverse electron velocity normalized by c, a L is laser vector potential normalized by m e c 2 /e.Assume the initial velocity of all newly born free electrons is zero and the jth electron is born at place of x j0 and time of t j0 .Then the transverse velocity of the electron is v ⊥, j (x,t) = a L (x,t) − a L (x j0 ,t j0 ) in a plane laser field.The average transverse velocity at x and t is given by where a L (x 0 ,t 0 ) = ∑ N j=1 a L (x j0 ,t j0 )/N, N is the total electron number at x and t, and a L (x 0 ,t 0 ) should be a function of x − ct.Assume the ionization stops at x f where laser pulses have left or gas atoms have been ionized completely.Then all a L (x 0 ,t 0 ) at x < x f should be statistically the same and therefore quasi-static transverse currents are formed in the wake of the laser pulses.
When laser pulses just leave x at t = x/c, the transverse velocity of all free electrons at x is assumed to be − a L (x 0 ,t 0 ) .At t > x/c, the electron velocity v ⊥ = δ a − a L (x 0 ,t 0 ) , where δ a is the vector potential of newly generated EM waves (THz waves).Here we have assumed δ a = 0 at t < x/c since the newly generated EM wave wavelength is much larger than the laser duration, which will be seen below.Thus, the wave equation of newly generated EM waves in plasmas is given by where ω 2 p = 4πe 2 n e /m e , the formed plasma electron density n e is assumed to be uniform since we only consider the case with plasma wavelength λ p = 2πc/ω p λ 0 and |δ a| 1. Solving Eq. ( 3) by the Laplace transform, one can obtain where } cannot be performed analytically, but it is known that the first term in the exponent represents the propagation effect along the x-direction and the second term means the decaying effect.Equation (4) shows that the THz field E ∝ ∂ δ a/∂t ∝ ω p , i.e., the square root of the plasma density.It also shows that the newly generated EM waves include both ω p and 0−frequency components and have the same polarization as the incident pulses.When the target has a gas-vacuum boundary, the ω p -frequency part within about λ p near the boundary can penetrate into the vacuum effectively.
The above analysis is applicable for gas targets irradiated by either one-color or two-color laser pulses with normalized amplitude a L < 1.When one-color waves are incident, the THz pulses will be produced very inefficiently and have some characteristics of noise.While twocolor waves are incident, the symmetry of positive and negative half-cycle of the fundamental wave is broken [13].The two-color-wave ionization causes drift currents in the same direction in every cycle and thus can emit THz pulses more efficiently than with one-color waves.When the target is a fully ionized plasma slab, either one-color or two-color waves are adopted, the net drift current is null, i.e. a L (x 0 ,t 0 ) = 0, and no THz pulses can be generated.

Simulations
To check the above results, we use one-dimensional (1D) particle-in-cell (PIC) simulations [3,17].In our PIC code the field ionization of gases is included with tunnelling ionization rates calculated by the ADK formula [15,16].For simplicity, only hydrogen gas (H 2 ) is used which generates one kind of plasma density only.A slab target with the gas density n gas = 1.25 × 10 −5 n c and the length L gas = 200λ 0 or 400λ 0 is used starting from x = 500λ 0 , where n c = m e ω 2 0 /4πe 2 is the critical density and ω 0 = 2π/T 0 is the laser frequency.A laser pulse is incident along the +x direction with the vector potential a L = êz a L , where a L = sin 2 (πξ /L t ) • [a 1 cos(k 0 ξ ) + a 2 cos(2k 0 ξ + θ )], k 0 = ω 0 /c, L t is pulse duration, a 1 and a 2 are normalized amplitudes for the fundamental and the second-harmonic waves, respectively, and θ is their relative phase.In the following simulations we take L t = 20λ 0 , a 2 = a 1 /2, and θ = π/2.
Figure 1 is the spacial distribution of transverse currents.In a preformed plasma, no transverse current is left after the laser.However, in a H 2 target a quasi-static negative transverse current is formed.At a later time the current will have a low-frequency oscillating profile [see Fig. 1(c)] because of the plasma modulation.It is indicated that ionization of gases is necessary for the THz emission in this case, which agrees with the above analysis.
The temporal and spacial distribution of THz emission is shown in Fig. 2(a).There are two THz pulses, one going along the +x direction and the other along the −x direction.These THz pulses penetrate into the vacuum and have a waveform with the first cycle strongest, for their current sources within the H 2 plasma decay during radiation.This can be seen more clearly in Fig. 2(b).Figure 2(b) is the temporal evolution of THz pulses observed at 20λ 0 in front of the left gas-vacuum boundary.Their corresponding spectra are plotted by Fig. 2(c).One can see that the THz pulses have frequencies of 0.003ω 0 (1THz), 0.005ω 0 , and 0.01ω 0 , equal to their individual ω p .Note that there is the 0−frequency component because of the asymmetry of the positive and negative peaks of the THz pulse.These are consistent with the previous analysis.For moderately intense incident pulses, experiments in Refs.[6,7,8,9,13] demonstrated that the THz pulse intensity increases almost linearly with incident pulse intensity.For even intense incident pulses, however, it is not clear whether this still holds.By PIC simulations, we are able to study the THz pulse amplitude as a function of incident laser amplitude within a large range.As shown in Fig. 3(a), at a 1 ≤ 0.06 the THz pulse amplitude increases with growing a 1 .However, when one increases a 1 further, the THz pulse amplitude does not increase with a 1 monotonically.The reason is that when a 1 is very small, the ionization process lasts for many laser cycles [see Fig. 3(b)] and the produced electron number and the drift current increase with growing a 1 .When a 1 is not small, the ionization process only occurs in the first few cycles and the leading edge of the incident pulse can ionize the gas completely [see Fig. 3(c)].The peak part of the pulse or the most intense part has no effect on the birth of free electrons or drift current excitation.This can result in the low efficiency in THz emission.When a 1 is increased further, the actual intensity at the part of the pulse ionizing gas may not grow (the part of the pulse ionizing gas is shifted forwards as shown in Figs.3(d) and (e)).As a result, the drift current or the THz emission may not strengthen with the increasing laser intensity.Therefore, for intense incident pulses it is difficult to control the born places of electrons and produced drift currents, and consequently the THz amplitude does not scale with a 1 monotonically.
It is noted that our above simulations are given with H 2 .When a complex gas such as a rare gas, molecular nitrogen, or air is taken, the features of the THz emission should be qualitatively the same as the H 2 gas.For the same incident laser pulse and gas density, the THz frequency should change due to the different ionization thresholds existing.It is possible that the drift currents produced by field ionization can distribute in a larger region throughout the laser pulse, and the resulting THz pulses have a broader frequency spectrum as compared with the simple H 2 gas.

Strong THz emission generation by chirped laser pulses
The unfavorable scaling shown above, which results from the ionization of gases in the twocolor scheme, should be overcome in order to achieve high-intensity THz emission.One can adopt another mechanism rather than the field ionization to generate large drift currents which Here we propose to use a single chirped laser pulse [18] irradiating on a gas or plasma target, so that electrons can get net acceleration in the transverse direction during each laser period.If fixing the chirping value sign of the laser, the net acceleration of electrons will be found in the same direction in all laser periods, which can cause very strong drift currents.After the laser passes, an electron obtains the transverse drift velocity of a L (L t ) − a L j0 ,t j0 ) according to Eq. ( 1), where a L (L t ) is the vector potential of the laser pulse at its trailing edge and a L (x j0 ,t j0 ) is relevant with the ionization process.When a large enough chirping value is taken, |a L (L t )| |a L (x j0 ,t j0 )| for a gas target, which closes to the case in a preformed plasma with a L (x j0 ,t j0 ) ≡ 0. For simplicity, we shall only take plasma targets with a L (L t ) − a L (x j0 ,t j0 ) ≡ a L (L t ) and Eq. ( 4) holds with substituting a L (L t ) for − a L (x 0 ,t 0 ) in this case.Therefore, obtained THz pulse amplitude will be in proportion with a L (L t ) or the incident laser amplitude.
In order to demonstrate the scheme mentioned above in a simple way, we take the chirped pulse electric field E L = êz ε L normalized by m e ω 0 c/e, with where the pulse duration L C t = ( 1 + 4CL 0 t − 1)/2C for the chirping parameter C and L 0 t is the pulse duration at C = 0.For positively-chirped pulses C > 0, L C t < L 0 t and for negatively-chirped pulses C < 0, L C t > L 0 t , when the total period number of laser pulses stays the same at any given C. Taking such chirped pulses, one can easily calculate the net transverse drift velocity dξ for any C.For example, when C = 1.0 and -0.024, there are a L (L C t ) = 0.0115a 0 and −0.0413a 0 , respectively.Then obtained THz pulse amplitude δ a ∝ 0.0115a 0 and −0.0413a 0 for C = 1.0 and -0.024, respectively, according to Eq. ( 4).To test the validity of this analysis, we conduct PIC simulations with C = 1.0 and C = −0.024,taking L 0 t = 10λ 0 .We fix the C or keep the chirped pulse profile and change the pulse amplitude to study the obtained THz pulse amplitude as a function of a 0 , which shown in Fig. 4(a).One see that the THz pulse amplitude increases linearly with a 0 for both positive and negative chirped pulses, which is in good agreement with the previous analysis.Further, the ratio of slopes  4(a).Moreover, compared with the scaling shown in Fig. 3(a), the THz pulse amplitude found for the chirped pulse with C = −0.024 is as 10-100 times large as that in the two-color scheme for a wide range of the incident laser amplitude.When a 0 = 0.067, the THz amplitude is up to 1MV /cm at the given plasma density of 2.5 × 10 −5 n c , corresponding to the frequency of 1.5 THz.
A two-dimensional (2D) PIC simulation is conducted to illustrate the multi-dimensional properties of the THz emission.We take the initial electron density 4 × 10 −4 n c (corresponding to the frequency of 6 THz), a chirped laser pulse with a 0 = 0.1, the transverse radius of 8λ 0 , and C = −0.024,as well as the simulation box of 800λ 0 × 80λ 0 (in x × y). Figure 4(b) is spacial distribution of THz pulses.One can see there are two THz pulses generated with one propagating along the +x direction and another along the -x direction.They have a transverse profile like that of the incident laser, which accords with our previous results that the THz pulse amplitude is proportional to the laser amplitude.The emission has a wavelength around λ p = 50λ 0 and it has the highest intensity in the first cycle with amplitude exceeding 7MV /cm.This amplitude is about 4 times of the one in Fig. 4(a) for the same a 0 , which agrees with Eq. ( 4) that the THz field is proportional to ω p .For the convenience of experimentalists, we rewrite here the 2D simulation parameters with common units that the plasma electron density 4.5 × 10 17 cm −3 , the laser wavelength 800nm, the peak intensity 2.1 × 10 16 W /cm 2 , the transverse radius 6.4μm, the duration 44 f s, and the chirped parameter C = −0.024.
It should be pointed that our 1D and 2D PIC simulation results presented previously have not  completely considered propagation effects of light, especially for large distance propagation.
In cases, the produced THz pulse and the laser pulse can evolve completely.The laser pulse may experience self-focusing provided the critical power P cr [19,20] is exceeded, while the THz pulse may spread very quickly due to small sport size comparable to the THz wavelength.As a result, the linear scaling of the THz emission amplitude with the initial incident laser amplitude could be deviated.

Conclusion
In conclusion, we have studied the mechanism of THz pulse generation from gas targets irradiated by ultrashort intense laser pulses.The quasi-static transverse currents created by laser field ionization and sequent modulation in formed plasmas are responsible for the THz emission.In the well-studied two-color laser scheme, it is found that the THz pulse amplitude scales linearly with the incident laser amplitude only when the latter is moderately large.This presents a significant limitation for the generation high power THz emission over the field strength of MV /cm, necessary for high field THz physics.To overcome this, we propose to use chirped laser pulses, with which even strong THz pulses can be generated with amplitudes scaling linearly with the laser amplitude.

Fig. 1 .
Fig. 1.Spacial distribution of transverse currents in a preformed plasma (a) at t = 60T 0 and in a H 2 target at t = 60T 0 (b) and t = 400T 0 (c), respectively.A two-color laser with a 1 = 0.06 and L gas = 400λ 0 are taken.

Fig. 2 .
Fig. 2. (a)Temporal and spacial distribution of THz pulses from a H 2 gas target with the density n gas = 1.25 × 10 −5 n c .Waveforms (b) and spectra (c) of THz pulses from H 2 targets with different initial gas densities.A two-color laser with a 1 = 0.06 is taken.

Fig. 3 .
Fig. 3. (a)THz pulse amplitude as a function of a 1 .Frames (b)-(e) show the number of new born electrons produced via laser field ionization against the laser propagation coordinate x − ct for different a 1 and a 2 = a 1 /2.Two-color lasers and a H 2 target are taken.

Fig. 4 .
Fig. 4. (a)THz pulse amplitude as a function of the chirped laser amplitude for two different chirping parameters, where the electron density is 2.5 × 10 −5 n c ; (b)Spacial distribution of THz pulses simulated by 2D PIC, where the electron density is 4 × 10 −4 n c and a chirped laser pulse with a 0 = 0.1, the transverse radius of 8λ 0 , and C = −0.024 is used.Plasma targets are taken.