We develop a one-dimensional model of THz emissions induced by laser-driven, time-asymmetric ionization and current oscillations in a hydrogen gas. Our model highlights complex scalings of the THz fields with respect to the laser and gas parameters, in particular, a non-monotonic behavior against the laser parameters. Analytical expressions of the transmitted and reflected fields are presented, explaining the THz spectra observed in particle-in-cell and forward-pulse propagation codes. The backward-propagating THz wave is mainly driven by the electron current oscillations at the plasma frequency, and its resulting spectrum operates below the plasma frequency. The transmitted THz wave is emitted from both plasma current oscillations and photo-ionization. Their respective signal presents a contribution below and around the plasma frequency, plus a contribution at higher frequencies associated to the photo-induced current. The interplay between these two mechanisms relies on the ratio between the propagation length and the plasma skin depth.
© 2014 Optical Society of America
Laser-gas interaction is an elegant way of energy conversion to higher frequencies through up-conversion mechanisms , as well as to lower frequencies through down-conversion processes . The reasonable size and high-repetition rates of the latest laser-based frequency converters make them highly promising for next-generation THz pulse devices [3, 4]. Spanning from the infrared to microwaves, the THz radiation spectrum is attractive for many applications in various fields of physics , biology and medicine , industry , security, communication, remote sensing  and basic science with, for instance, molecular dynamic spectroscopy [9, 10].
Initially attributed to optical rectification and four-wave mixing via a third-order nonlinearity [11, 12], THz emission has been widely discussed in the optics and plasma literature [13–16]. Recently, a consensus has been reached on the so-called two-color scheme, which allows to noticeably increase, by around two orders of magnitude, the THz conversion efficiency in laser-gas interaction above the ionization threshold . The principle of this scheme is to combine a fundamental frequency carried by an ultrashort (femtosecond) pulse with its second harmonic within a highly nonlinear plasma spot. As a result of the subtraction of broad optical bandwidths, new electromagnetic components are generated with spectra centered around THz frequencies. Remarkably, similar THz emissions over large distances have also been reported from long-range filamentation of femtosecond pulses in air or in noble gases (e.g., argon) [18–22]. Recent works also pinpoint local plasma-induced photocurrents as key players for THz generation in the filamentation regime [23, 24]. THz pulses are then emitted during and after the gas ionization. Firstly, bursts of attosecond currents are produced around each ionization event. For a specific dephasing between the two laser beams, these bursts interfere constructively to generate THz pulses of the order of the ionization duration . Moreover, transverse current oscillations are driven at the plasma frequency, which prolongs the THz emission after the laser interaction . The resulting THz field reaches amplitudes of the order of several MV/cm, much higher than those produced by four-wave mixing and comparable to those usually attained in tightly-focused geometries. In two-color filamentation experiments, the pulse intensity is currently limited to 50 – 100TW/cm2. Under these conditions, the THz yield appears to increase with respect to the pump energy, up to some saturation [18, 19]. THz yields up to 100 μJ were recently reported for pump energies at the Joule level . The enhanced THz conversion efficiency (∼ 10−4) was attributed to meter-range filaments that maintain long plasma channels.
Out of a filamentation scenario, Particle-In-Cell (PIC) simulations [26, 28], solving Vlasov-Maxwell equations by means of macro-particles, have shown that rising the laser intensity beyond the previous values does not necessarily increase the THz amplitude. While essentially driven by ionization currents, the THz signal also depends on the noble gas serving as the interaction medium and on the pulse duration and carrier-envelope phase, which impose a saturation value for the ionization current. For instance, the large ionization potential in helium was exploited to obtain THz field amplitudes as high as 1.28 GV/m with moderate laser energies. A careful choice of the plasma and laser parameters is thus required to achieve efficient THz production. When computed inside a plasma of constant density, the radiated field is dominated by ionization currents . By contrast, the contribution of the plasma current oscillations seems to prevail outside the plasma region and/or at high pump intensities [29, 30].
In this regard, the THz spectral signature attributed to photocurrents may be different from that produced by plasma current oscillations along the propagation direction. While the former operate over a broad range of low frequencies [16, 31], the latter more specifically operate around and below the plasma frequency , that spans the THz region for electron densities < 1019 cm−3. Thus, a clear distinction between the potential sources of THz radiation involved in the dynamics of the free electrons still remains to be drawn. Another open issue concerns the THz energy partition between the transmitted and reflected waves. Several studies have assumed that the largest fraction of the THz energy in the photocurrent model is carried by the forward wave component [23–25,31]. Straightforward comparisons between backward and forward spectra computed from a Unidirectional Pulse Propagation Equation (UPPE) model and two-dimensional (2D) Finite-Difference Time-Domain (FD-TD) Maxwell code have validated the unidirectional approach for an argon gas . For pump intensities less than 100 TW/cm2, 2D FD-TD Maxwell simulations predict mostly forward THz emission. This forward radiation is generated predominantly at the ionization front and is not affected by plasma opacity, unlike the backward spectrum. However, when a single laser pulse of much higher intensity (up to 1017 W/cm2) interacts with a 100 μm-thick hydrogen plasma, a noticeable fraction of the THz energy, produced by the plasma current oscillations, can be conveyed by the backward wave . In this regime, alternative key agents for THz generation may be plasma wakefields and the laser wavelength. In this respect, several studies predict strongly enhanced THz yields in the case of single- or two-color pulses operating at long central wavelengths (e.g., 2μm or 4μm) [23, 26]. Here, several effects superimpose, such as the respective influences of the pump laser wavelength and the effective pulse duration, as well as the interference between the low-frequency tail of the pump spectrum with the inner THz spectrum at long enough pump wavelengths.
All of the previously published models of laser-driven THz sources have failed to clarify the connection between the two main “plasma” mechanisms, i.e., ionization-induced photocurrents and current oscillations, particularly the link between forward and backward waves, and their dependencies upon the laser and plasma parameters. The goal of the present work is thus to propose a unified framework accounting both for photocurrents and current oscillations, in order to extend the existing models to laser intensities ≲ 1017 W/cm2. To this purpose, we shall revisit the well-known photocurrent mechanism  by separating the THz component from the laser pump field and by including the plasma current oscillation in a self-consistent way. For the sake of simplicity, we shall introduce two semi-analytical 1D models, either discarding or retaining propagation effects, and restrict our analysis to a hydrogen gas. Section 2 introduces our analytical model, which assumes decoupled THz and laser fields. In section 3, we derive scalings of the THz emission against the laser and plasma parameters within the assumption of negligible propapagation effects. Despite this approximation, our scalings perfectly reproduce the tendencies observed in PIC simulations . The choice of optimal parameters for THz emissions is demonstrated in the case of one and two laser colors. In section 4, propagation effects are included. Analytical formulas for the reflected and transmitted THz waves are derived. We demonstrate that the forward wave carries most of the THz energy, provided that the plasma length exceeds the plasma skin depth. Besides, the dependence of the THz amplitude upon the photocurrent amplitude is clearly evidenced and its limitations are discussed. The competition between photocurrents and plasma current oscillations is analyzed, and shown to depend on the ratio between the propagation length inside the gas and the plasma skin depth. In section 5, the backward and forward spectra inferred from our general model are found in remarkable agreement with both PIC and UPPE simulation results.
2. 1D model for laser-driven THz pulses
In this section, we present a simple semi-analytical model predicting the THz emission in gases for a large set of parameters. Drawing upon Refs. [30, 31], it involves two mechanisms, both attributed to the transverse electron current oscillations induced during and after gas ionization by one- or two-color laser pulses. For the laser intensity range scanned in the present article (∼ 1014 − 1017 W/cm2), both the nonlinear polarization (bound electron response) of the medium and the relativistic effects are neglected. Instead, we focus on the interplay of laser-driven photoionization [4, 17, 31] and transverse current oscillations [26, 30].
In the following, the time and space coordinates, (t, x), are normalized to 1/ω0 and λ0/2π, respectively, where ω0 is the pump laser frequency and λ0 is the pump central wavelength. Let us consider a semi-infinite gas of neutral density na, filling the positive x-region and a transverse electric field, , composed of one- or two-color laser pulses propagating toward the gas. The gas density and the laser electric field are normalized to the critical density and the Compton field meω0c/e, respectively. For low gas densities, na ≪ 1, we can neglect the laser dispersion and energy loss during the ionization processes. We hence assume the laser electric field, generally embedding two colors, to maintain a constant shape as2, 12, 27]. Therefore, for technical convenience, we shall use equal pulse durations as in Refs. [28, 30] in our theoretical analysis. Second harmonic pulse lengths being twice smaller than the pump duration will instead be employed in section 5, when confronting theoretical predictions to direct numerical simulations. Initially, the finite laser pulse is located in front of the gas, between x = −2πτ and 0. Writing the full electric field as 33–35]. For hydrogen, the latter is defined as Equations (4) and (5) recast as Equations (3) and (7) are solved numerically and approximated analytically in the next sections.
3. Non-propagating THz pulses
In this section, we wish to obtain scaling laws for the THz emissions in terms of the plasma and laser parameters. For simplicity, we first neglect the wave propagation ( ), which reduces the above equation to an ordinary differential equation:Fig. 1(a) for a single-color laser pulse (a0 = 0.02, τ = 10, λ0 = 1μm) and negligible collisions (ν = 0). During the ionization process (t < 2πτ), the newborn electrons oscillate in the laser pulse, producing the field at frequency ω0. In contrast to Refs. [4, 17, 31], the remaining current oscillates at the plasma frequency (normalized to ω0), as described in  (note the absence of damping mainly due to the absence of propagation). This is illustrated in the field spectrum (black line) displayed in Fig. 1(b). For the selected gas density, na = 0.0011, complete ionization takes place, so that the main mode emerges at the plasma frequency ωpe ≃ 0.033. We have applied a low-frequency filter, ω < 0.3, to calculate the THz field amplitude. As an example, the filtered field and its spectrum are plotted in red dotted line in Fig. 1(a) and (b), respectively. This methodology has been used to explore the dependencies of the THz field maximum upon the laser and/or gas parameters.
The maximum of the THz field normalized to , versus the laser intensity and for different set of parameters is presented in Fig. 2. On panel (a), the THz field maximum is plotted for different laser wavelengths, keeping constant the number of laser cycles in the pulse. For small laser intensities, the gas is weakly ionized, yielding negligible THz emission. When the laser field is strong enough to fully ionize the gas (I > 2×1014 Wcm2), the THz field increases. At very high intensities, ionization occurs at the very beginning of the pulse and thus, due to weakened photocurrents, the THz field amplitude drops. At leading order, the normalized THz amplitude appears weakly sensitive to the laser wavelength, obeying the scaling28,30]. In physical units the THz field energy is proportional to , where ngas is the gas density. As already reported in Ref. , the THz field amplitude, proportional to a0, “oscillates” when increasing the laser intensity. A similar plot is presented in Fig. 2(b) for other parameters. The black solid line corresponds to a reference case (λ0 = 1μm, τ = 10, na = 0.0011). For a fixed number of laser cycles, the normalized THz field is weakly changed when increasing the gas density by a factor 4 (red dashed line) or when doubling the laser frequency (dashed-dotted blue line). This trend confirms Eq. (10).
In addition, the THz field amplitude is substantially reduced when increasing the number of laser cycles in the pulse, as observed for τ = 40 (green solid line). This suggests (at least) a linear decrease in the THz amplitude with the number of optical cycles τ. To further illustrate this dependency, the THz field amplitude is plotted for the case λ0 = 2μm, τ = 5, and na = 0.0044. This situation corresponds to doubling the laser wavelength, while keeping constant the gas density and the laser pulse duration. The resulting pattern (cyan solid curve) does differ from the reference case (black solid line). The number of oscillations noticeably decreases and their amplitude grows from smaller laser intensities, whereas the peak values of is multiplied by about a factor 2. This result suggests a strong sensitivity of the THz amplitude to the effective number of optical cycles, more pronounced in the case of few-cycle pulses owing to strengthened time-asymmetry and net ionization currents, in agreement with Ref. . Therefore, when considering a single color, we can propose from the previous behaviors a scaling of the THz amplitude as
Apart from this leading-order scaling, the oscillations seen in Fig. 2 introduce a margin of uncertainty of 2 to 4 when doubling the central laser wavelength. This additional refinement completes the above scaling by a factor > 2, which then renders it comparable to the variations reported in Ref. . Here, the THz field amplitudes were seen to drop by a factor ∼10 when passing from 2 μm to 1 μm wavelength using a pump pulse of constant (50 fs) duration and intensity. Similar variations are also compatible with the 14-fold growth in the THz yield reported in Ref.  when comparing two-color filaments with 0.8μm and 2μm pumps of same duration but for different atomic densities (local pressure = 6.44 bars resulting in a net scaling rate of about 15.9).
We have also examined the THz field behavior with the phase shifts ϕ and ϕ2. For one laser color (Fig. 3), the THz field amplitude (black solid line) is obviously π-periodic as . However, a more complex behavior arises when using two colors (red dashed line). The THz field amplitude is now 2π-periodic, with a maximum value reached for ϕ2 ≈ 4.1 rad in hydrogen with r2 = 0.2. Importantly, the THz field amplitude is strongly enhanced when using two colors, as illustrated in Fig. 4(a), where, typically, an order of magnitude increase in intensity is found compared to Fig. 2(a). The optimum THz intensity is obtained for the laser intensity I0 ∼ 3 × 1014 W/cm−2 for τ = 10 and ϕ2 = 0. By comparison with Fig. 2(b), the THz intensity attained with one color at similar pump intensity is at least 100 times weaker, thereby confirming the two orders of magnitude often mentioned in gases [18–20]. Note that the two-color configuration also gives rise to an oscillatory intensity dependence of the THz field. This dependence is sensitive to the phase ϕ2, as illustrated by the red dashed line (ϕ2 = 0.5) and the cyan solid line (ϕ2 = 1.5). By increasing the number of optical cycles to 40, the oscillations relax (see the dashed-dotted green line), while the optimum laser intensity is shifted. Although a scaling in 1/τ may not apply, the THz field still remains sensitive to the number of laser cycles. The interplay between the harmonic amplitude ratio r and the phase ϕ2 is shown in Fig. 4(b). Here again, the THz field amplitude exhibits non-monotonic variations with r. The maximum THz field is obtained for r ≃ 0.85 and ϕ2 ≃ 4.7.
The non-monotonic behavior of the THz field and the previous scalings can be explained with a qualitative analysis of Eqs. (3) and (8). In the short-pulse limit, , the solution of Eq. (8) is approximated as
The oscillatory dependence of the THz field upon the laser intensity, and its sensitivity to the laser duration then follows from estimating the source term31], the electron density evolution can be calculated knowing that the ionization events are localized around the electric field extrema. The electron density is then approximated as Eq. (3), we obtain the recurrence relation
As the laser electric field derives from a potential with a mean value equal to 0 [ ], the THz source term becomes
For one laser color, Eq. (1) recasts as . As τ ≫ 1, the laser envelope derivative is small g′/g ∼ 1/τ ≪ 1, and the field extrema are reached at , k ∈ N. Since g(0) = 0, the vector potential is , and the THz source term is thus proportional to G ∝ naa0/τ, thereby yielding Eq. (11). For one laser color, the higher the number of laser cycles, the weaker the THz field amplitude. The THz amplitude behavior against laser intensity can be explained as follows. For small laser intensities, such as νE [Em,k] ≪ 1, ionization occurs around the laser pulse maximum. As the laser field is anti-symmetric around the maximum, the source term is exactly equal to 0. At higher intensities, there exists a time, tk < πτ, such that νE [Em,k] > 1. In this case, ionization occurs stepwise at several field extrema.
For simplicity, let us assume a two-step ionization: . We look for the periodic zeroes of the THz source. The THz field vanishes if G = 0. As δne,k−1 + δne,k = na and , one finds the relation36] and (18) is weakly sensitive to the laser wavelength, consistently with Fig. 2(a), but depends on g−1 ( ) (we remind that Ea.u. ∝ λ0). This variation with the inverse function of the laser envelope for, e.g., a Gaussian or sin2 profile, accounts for the increase of the oscillation wavelength with the laser intensity seen in Fig. 2. For stronger laser intensities, the ionization events shift from the extremum k to the extremum k − 1 until reaching the first extremum (k = 1). Hence, the number of oscillations is proportional to the number of extrema included in half of the laser duration, as observed in Fig. 2(b) (green solid line) where the number of oscillations has increased by a factor ∼ 4.
When the ionization process is asymmetric, as in the two-color configuration, the scaling of the source term G changes. Assuming ϕ = 0, one hasFig. 4(a), this event corresponds to I ∼ 5 × 1014 W/cm2 for τ = 10. Maximum THz emission varies with the number of laser cycles: for τ = 40, it is shifted down to the value I ∼ 2 × 1014 W/cm2. Reversely, at fixed intensity (e.g., I = 1014 W/cm2) and wavelength, a 4-fold increase in the pulse duration raises the THz intensity by a minimum factor of ∼36. Up to propagation and multidimensional aspects discarded here, such impressive growths are comparable with those reported in Refs. [18, 23].
In conclusion to this section, our simplified 1D model predicts that the ionization-induced THz emission is utterly defined by the electron current amplitude left after the ionization process, Jy ∼ G. As already known, the stronger the ionization asymmetry, the higher the THz field amplitude. However, the origin of the asymmetry and its influence change between the one- and two-color configurations. Besides, as we neglect the THz field propagation, the THz spectrum is mainly located around the plasma frequency. This anomalous result is corrected in the next section, where propagation aspects are taken into account.
4. Propagating THz pulses
This section is devoted to the forward and backward THz emissions, and related modifications of the spectrum. For that purpose, we analytically and numerically solve the interaction of a laser pulse with a semi-infinite gas assuming that the laser pulse remains unperturbed while propagating.
The laser pulse defined by Eq. (1), initially located in the vacuum region, x ∈ [−2πτ, 0], interacts with a semi-infinite gas of density profile na(x) = naH[x]. The THz field equation is given in Eq. (7), with the electron density ne described by the continuity equation, Eq. (3). The latter equation is numerically integrated using a 4th-order Runge Kutta method, while the THz wave equation is solved using a 2nd-order explicit centered scheme. A solution is presented in Fig. 5 in black solid line in the panels (a) and (b) for a single color. The THz field presents a small reflected wave, visible in panel (b) at the distance x ∼ −900, with an amplitude ∼ 10−4 for the chosen parameters, and an intense forward wave of amplitude ∼ 10−2 centered around x = 900 [panel (a)]. The reflected wave is similar to that observed in PIC simulations [26, 30]. Basically, the laser field ionizes the gas, hence producing a net transverse current. As discussed in section 3, this current oscillates at the plasma frequency, which maintains the THz emission during several plasma periods. Inside the gas, the laser propagates and pursues the ionization process. The resulting field δEy propagates also forward and increases quasi-linearly with the distance.
For a more quantitative description, we now derive an approximate formula for the propagating field δEy. Without collisions (ν → 0), the wave equation, Eq. (7), recasts as:Eq. (3) expresses as 28,37]. In the opposite case, to include the ionization effect, we choose to solve Eq. (22) in two distinct domains. In the beam head, β+ ∈ [−2πτ, 0], we neglect the THz field compared to the laser field . This approximation is valid for pulse duration and time such that , i.e., such that the THz dispersion is negligible during the ionization process and the propagation distance. In the second region, behind the laser pulse, the electron plasma density is constant and Eq. (22) is solved using the Laplace transform.
In the beam head, the wave equation is approximated asAppendix A, we find the solution Eq. (24), F is the current integral while G is related to the current density. As F ∼ 2πτG, the second quantity, affected by a linear growth factor in x, prevails in the forward component for x ≫ 2πτ.
Behind the laser pulse ( ), Eq. (22) reduces toEq. (25) becomes (24) becomes . The solution is thus 36], the Bromwitch inversion gives
Expression (28) reproduces the main functional dependencies derived earlier in Ref. , but without the F-function. Here, we go one step beyond. Knowing the solution δEy(x, t) inside the plasma, the reflected wave , satisfying the equation , is deduced following the characteristic x + t = t0, i.e., , where δEy(x, t) is given by Eqs. (24) and (28). For 0 < x + t < 2πτ and x < 0, the reflected wave is
Comparing this expression with Ref. , we can infer that Eq. (30) describes the backward propagation of a distortionless wave oscillating at the plasma frequency with an electron density fixed by nef. Owing to the bi-directional character of the 1D propagator, the THz backward emission solely follows from the current oscillations over the plasma skin depth. At large time t, the backscattered solution behaves as ; its spectrum should convey a dominant mode at plasma frequency accompanied by damped, thus smaller frequencies. In the limit , the first right-hand side term of Eq. (30) is dominant, hence confirming the scaling suggested in section 3. Moreover, at late times t ≫ x + 2πτ and near the plasma border for x > 0, the field simplifies toEq. (24), in the ionization region, the amplitude of the forward-propagating wave is proportional to the source term G and increases linearly with the distance. As we neglect the field dispersion in the beam head, this linear increase is valid for . Behind the laser pulse, where the dispersion is included, the electric field scaling, at a given distance β+ from the beam head, turns into Appendix B). As a result, THz scalings can vary with respect to the location of the propagated THz wave.
To validate our analytical formulae, we plot the solution in Fig. 5 in red solid line. The agreement in the plasma region is remarkable, especially in the beam head [Fig. 5(a)], where the small difference is ascribed to the THz wave dispersion, i.e., the neδEy term in Eq. (20). To illustrate this effect, we have numerically solved Eq. (7) neglecting the term ∂tδayne in the laser region x − t ∈ [−2πτ, 0]. The resulting curve (dashed black line) fairly agrees with the analytical solution both in the beam head and in the plasma region far from x = 0. Although capturing the right oscillation frequency, the analytical expression of the reflected field yields small discrepancies in the amplitude values, which we attribute to the fact that Eq. (28) is valid for an infinite plasma only.
As for the reflected wave, the forward wave oscillates at the plasma frequency behind the beam head. Yet, far from the laser pulse, the wavelength of the emitted field increases. This effect is illustrated in Fig. 6 where the numerical solution is plotted against x and t. In vacuum, the reflected wave propagates at the velocity −1, while in the plasma, the emitted wave propagates at the velocity +1 right behind the laser pulse. Far from the pulse, the longitudinal wave number relaxes to kx = 0. Indeed, according to Eq. (31) and since x ≫ t, the space derivative is ∂xδEy ∼ O(1/t), hence supporting the approximation ∂x ≪ ∂t used in section 2.
To visualize the propagation effect on the THz spectra compared with the “ideal” spectrum of Fig. 1(b), we have run a similar case to Fig. 5 but with a finite rescaled plasma length of 40π. The spectra of the reflected and transmitted waves are plotted in Fig. 7 in red and black solid lines, respectively. The intensity spectrum of the reflected wave is approximately a step function with a cut located at the plasma frequency . Below this frequency, photocurrents weakly contribute to THz emission. Given the plasma thickness, the source term, Eq. (24), mainly depends on xG(x − t) ≫ F(x − t) inside the plasma. As the reflected wave is driven by F(x + t) outside the plasma, the energy carried by the transmitted wave is higher and mainly determined by the ampitude ratio 2xG/F. Besides, its spectrum now extends above the plasma frequency. This THz field is attributed to both the constructive interferences during the ionization process and the plasma current oscillation. The electron current during ionization presents frequencies around 2/τ, which broadens the THz spectrum above the plasma frequency. The spectra deduced from the analytical expressions Eqs. (24), (28) and (30) are plotted in red and black dashed lines for the transmitted and reflected waves, respectively. A good agreement is obtained with the numerical solution, especially concerning the transmitted wave. The reflected wave presents a plateau between 0 and , with a peak pronounced around the plasma frequency. Above this frequency, the spectrum dramatically decreases by many orders of magnitude. The intensity spectrum of the transmitted wave, in black dashed line, presents similar features to the numerical solution. In particular, the contribution of the frequencies higher than ωpe, typical of ionization-induced THz emission, has increased. However, for the gas thickness of 40π, the electron current oscillation at the plasma frequency is the dominant process. Indeed, the spectrum above , mainly depends on the source term xG(x − t), while the low-frequency spectrum depends on the plasma current oscillations, driven by the term .
Hence, for propagation length much larger than , the photocurrent mechanism becomes dominant, whereas for propagation length comparable to, or smaller than the plasma wavelength, as in Fig. 7, the main mechanism is the plasma current oscillation. To confirm this behavior, we have calculated the transmitted spectrum for a thicker gas of L = 100π. The numerical and semi-analytical solutions are plotted in green dashed line and green solid line, respectively. In this regime, the propagation length is such that , which leads to a spectrum contribution higher for , as a result of dominant photocurrents. This tendency is reproduced by the semi-analytical model.
In summary, we have demonstrated within a unified formalism that laser-driven THz radiation can originate from different sources, namely, photoionization [23, 31] and residual current oscillations [26,30]. Our model establishes the missing link between the extreme nonlinear optics and plasma physics communities. The nontrivial interplay between these two mechanisms appears to be highly sensitive to the laser and medium parameters. For a propagation length shorter than the plasma wavelength, THz emission mainly results from plasma current oscillations. For longer gases, THz emission mostly originates from the constructive interferences during gas ionization.
To further validate our model, we have run 1D PIC simulations using the code calder . The code resolves the wave propagation inside the plasma, the particle trajectories including relativistic effects and field ionization . To reproduce the conditions used in Ref. , collisions are neglected. The initial density profile of the hydrogen gas is trapezoidal, with a 90 μm (700c/ω0)-long plateau bordered by 5 μm-long linear ramps. The maximum atomic density is na = 0.0044, yielding an electron density of ne = 1.2 × 1018 cm−3. The initial ion temperature is set to 1 eV. The laser beam has a Gaussian intensity profile with a 35 fs FWHM duration. The laser wavelength is λ0 = 2μm and the maximum intensity is 1017W/cm2. The simulation box is discretized in 9000 cells, each one containing 1000 electrons and 1000 atoms of hydrogen. The numerical resolution is Δx = 16 nm and Δt = 0.048 fs. Calculations have been performed using 3rd-order weight factors, with absorbing conditions for the fields.
The spectra of the transmitted and reflected waves are shown in Fig. 8(a) and compared to those obtained by the simplified wave equation, Eq. (7), assuming an unperturbed laser pulse. At the intensity of 1017W/cm2, both reflected spectra present a plateau at frequencies < ωpe (f ≡ ω/2π < ωpe/2π ≃ 10THz), due to the dominant contribution of the plasma current oscillations. In contrast, the transmitted wave presents a broadband spectrum in the range [ωpe, ω0] (10 ≤ f ≤ 150THz). Despite the single-color setup considered, the ionization occurring at the beginning of the pulse is asymmetric and produces a net transverse current. As discussed previously, this photoinduced current contributes to the spectrum well above the plasma frequency. Since the gas length (∼ 180π, corresponding to 90μm) exceeds the plasma skin depth (15, corresponding to 5μm), this radiation mechanism prevails, so that the THz energy is mainly contained at frequencies . Figure 8(b) compares the resulting THz fields computed below ω = 0.2 (f = 30 THz) from the transmitted and reflected THz waves obtained from the 1D calder code and from our simplified wave equation. This figure displays an excellent agreement between PIC simulations and our theoretical approach.
We have also compared our reduced numerical model to PIC and uppe1d simulations  in a two-color configuration for much weaker intensities. The numerical code solving the Unidirectional Pulse Propagation Equation model in one-dimensional (1D) geometry, uppe1d, is generally used to simulate the forward propagation of ultrashort laser pulses through gas-filled capillaries, accounting for Kerr optical effect and gas ionization. This code neglects relativistic and diffraction effects, but it includes the transverse current oscillation at the plasma frequency. In our simulations the laser intensity is set to 1014 W/cm2 or 1015 W/cm2 with r2 = 0.1, ϕ2 = 0 and λ0 = 1μm. The fundamental and frequency-doubled pulses have Gaussian intensity profiles with ~30 fs and ~15 fs FWHM durations, respectively. A 100 μm-long gas is considered. The other parameters are those used in the one-color case. Here, the Kerr polarisation, having a small impact in 1D geometry, is switched off in uppe1d. The transmitted wave spectra are presented in Figs. 9(a) and 9(b) from weak to strong pump intensity. Albeit overestimated by the simple model, which neglects laser dispersion, the agreement is satisfying. The contributions of both the current oscillation ( ) and the non-symmetric ionization ( ) are correctly reproduced. A striking feature is the increase of the THz spectrum at high intensities reproducing a spectral pattern comparable with that of Fig. 8(a), i.e., THz emission induced through photocurrents assures the merging with the pump field spectrum.
To end with, we have tested the robustness of our theoretical model by performing three-dimensional (3D) simulations that now include transverse diffraction and Kerr self-focusing. The blue dashed-dotted curves in Figs. 9(a) and 9(b) show the spectra computed at 1014 W/cm2 or 1015 W/cm2 pump intensities, respectively, using the uppe3d model of Ref. . Under similar initial conditions, the two-color field propagates over 100μm with a Gaussian transverse amplitude profile of 20μm radius (1/e2). We can observe that the 3D spectra remain close to their 1D counterparts and, thereby, they still remain strongly connected to the physics captured by our simplified model. With a 1014 W/cm2 pump, the laser field slowly diffracts and, in turn, the THz spectrum decreases in intensity. With a 1015 W/cm2 pump, the laser field starts to experience Kerr self-focusing, which contributes to enhance the THz spectrum to some extent.
6. Conclusions and perspectives
In conclusion, we have developed two models of THz emissions induced by laser-gas interaction. The simplest model describes the plasma current oscillation induced during gas ionization. Some scaling laws for the THz emission in terms of the main laser and gas parameters have been derived and found in good agreement with previous numerical studies [26,31]. A more sophisticated model has been developed, which demonstrates that the previous scalings still mostly hold when accounting for propagation effects. This model is able to reproduce the reflected and transmitted wave spectra obtained by self-consistent numerical simulations. Whereas the reflected spectrum is mainly determined by the current oscillation at plasma frequency, the transmitted spectrum is a subtle mix between plasma current oscillation- and ionization-induced THz emissions. For plasma lengths shorter than the plasma skin depth, the former contribution is dominant, while the latter becomes significant in the opposite case. At leading order, the THz field scaling is for one laser color and for two colors, for which the number of optical cycles in the overall pulse matters in a seemingly more complex dependency. Therefore, because the actual THz emission amplitude is oscillating and therefore not monotonic against laser intensity and the number of laser cycles, a simple, systematic dependence of the THz emission upon the laser wavelength may not be achievable. The present analysis is valid in the non-relativistic regime, upon 1D symmetry assumption. A more detailed model is, however, required to prospect the effect of the J × B force and the induced plasma wave on the THz emissions. Furthermore, one should explore different gases from hydrogen, as the THz field strength and related spectra are expected to vary with the ionization (binding) energy and occurence of multiple ionization.
A. Computing Eq. (24)
We look for a solution δEy, in the region β+ ∈ [−2πτ, 0]. Introducing B = ∂sδEy, Eq. (23) becomesEq. (24).
B. Asymptotic behavior Eq. (32)
Behind the beam head, the scaling for the THz amplitude with the crossed thickness can be obtained in the limit t(= s) ≫ β+ ≫ 2πτ, with β+ = cste. As , Eq. (28) recasts as36]. Hence, the THz field amplitude at given distance β+ from the laser pulse scales as
This work was granted access to the HPC resources of TGCC and CINES under the allocation 2013-x2013052707 made by GENCI (Grand Equipement National de Calcul Intensif).
1. F. Brunel, “Harmonic generation due to plasma effects in a gas undergoing multiphoton ionization in the high-intensity limit,” J. Opt. Soc.Am. B 7, 512–526 (1990).
2. K. Y. Kim, A. J. Taylor, J. H. Glownia, and G. Rodriguez, “Coherent control of terahertz supercontinuum generation in ultrafast laser-gas interactions,” Nat. Phys. 2, 605–609 (2008).
3. K. Reimann, “Table-top sources of ultrashort THz pulses,” Rep. Prog. Phys. 70, 1597–1632 (2007). [CrossRef]
4. H. Roskos, M. Thomson, M. Kreß, and T. Löffler, “Broadband THz emission from gas plasmas induced by femtosecond optical pulses: From fundamentals to applications,” Las. Phot. Rev. 1, 349–368 (2007). [CrossRef]
5. M. Tonouchi, “Cutting-edge terahertz technology,” Nat. Photonics 1, 97–105 (2007). [CrossRef]
6. R. M. Woodward, B. E. Cole, V. P. Wallace, R. J. Pye, D. D. Arnone, E. H. Linfield, and M. Pepper, “Terahertz pulse imaging in reflection geometry of human skin cancer and skin tissue,” Phys. Med. Bio. 47, 3853–3863 (2002). [CrossRef]
7. N. Nagai, M. Sumitomo, M. Imaizumi, and R. Fukasawa, “Characterization of electron- or proton-irradiated Si space solar cells by THz spectroscopy,” Semicond. Sci. Technol. 21, 201–209 (2006). [CrossRef]
8. J. Liu, J. Dai, S. L. Chin, and X.-C. Zhang, “Broadband terahertz wave remote sensing using coherent manipulation of fluorescence from asymmetrically ionized gases,” Nat. Photonics 4, 627–631 (2010). [CrossRef]
10. B. M. Fischer, M. Walther, and P. Uhd Jepsen, “Far-infrared vibrational modes of DNA components studied by terahertz time-domain spectroscopy,” Phys. Med. Bio. 47, 3807–3814 (2002). [CrossRef]
11. D. J. Cook and R. M. Hochstrasser, “Intense terahertz pulses by four-wave rectification in air,” Opt. Lett 25, 1210–1212 (2000). [CrossRef]
14. M. Kress, T. Löffler, 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). [CrossRef] [PubMed]
15. 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 (2007). [CrossRef]
16. J. Peñano, P. Sprangle, B. Hafizi, D. Gordon, and P. Serafim, “Terahertz generation in plasmas using two-color laser pulses,” Phys. Rev. E 81, 026407 (2010). [CrossRef]
18. T.-J. Wang, Y. Chen, C. Marceau, F. Théberge, M. Châteauneuf, J. Dubois, and S. L. Chin, “High energy terahertz emission from two-color laser-induced filamentation in air with pump pulse duration control,” Appl. Phys. Lett. 95, 131108 (2009). [CrossRef]
19. T.-J. Wang, S. Yuan, Y. Chen, J.-F. Daigle, C. Marceau, F. Théberge, M. Châteauneuf, J. Dubois, and S. L. Chin, “Toward remote high energy terahertz generation,” Appl. Phys. Lett. 97, 111108 (2010). [CrossRef]
20. T.-J. Wang, J.-F. Daigle, S. Yuan, F. Théberge, M. Châteauneuf, J. Dubois, G. Roy, H. Zeng, and S. L. Chin, “Remote generation of high-energy terahertz pulses from two-color femtosecond laser filamentation in air,” Phys. Rev. A 83, 053801 (2011). [CrossRef]
21. J.-F. Daigle, F. Théberge, M. Henriksson, T.-J. Wang, S. Yuan, M. Châteauneuf, J. Dubois, M. Piché, and S. L. Chin, “Remote THz generation from two-color filamentation: long distance dependence,” Opt. Express 20, 6825–6834 (2012). [CrossRef] [PubMed]
22. O. G. Kosareva, N. A. Panov, R. V. Volkov, V. A. Andreeva, A. V. Borodin, M. N. Esaulkov, Y. Chen, C. Marceau, V. A. Makarov, A. P. Shkurinov, A. B. Savel’ev, and S. L. Chin, “Analysis of Dual Frequency Interaction in the Filament with the Purpose of Efficiency Control of THz Pulse Generation,” J. Infrared Milli Terahz Waves 32, 1157–1167 (2011). [CrossRef]
23. L. Bergé, S. Skupin, C. Köhler, I. Babushkin, and J. Herrmann, “3D Numerical Simulations of THz Generation by Two-Color Laser Filaments,” Phys. Rev. Lett. 110, 073901 (2013). [CrossRef]
24. A. V. Borodin, N. A. Panov, O. G. Kosareva, V. A. Andreeva, M. N. Esaulkov, V. A. Makarov, A. P. Shkurinov, S. L. Chin, and X.-C. Zhang, “Transformation of terahertz spectra emitted from dual-frequency femtosecond pulse interaction in gases,” Opt. Lett 38, 1906–1908 (2013). [CrossRef] [PubMed]
25. I. Babushkin, S. Skupin, A. Husakou, C. Köhler, E. Cabrera-Granado, L. Bergé, and J. Herrmann, “Tailoring terahertz radiation by controlling tunnel photoionization events in gases,” New J. Phys. 13, 123029 (2011). [CrossRef]
26. W.-M. Wang, S. Kawata, Z.-M. Sheng, Y.-T. Li, L.-M. Chen, L.-J. Qian, and J. Zhang, “Efficient terahertz emission by mid-infrared laser pulses from gas targets,” Opt. Lett 36, 2608–2610 (2011). [CrossRef] [PubMed]
27. T. I. Oh, Y. S. You, N. Jhajj, E. W. Rosenthal, H. M. Milchberg, and K. Y. Kim, “Intense terahertz generation in two-color laser filamentation: energy scaling with terawatt laser systems,” New J. Phys. 15, 075002 (2013). [CrossRef]
28. M. Chen, A. Pukhov, X.-Y. Peng, and O. Willi, “Theoretical analysis and simulations of strong terahertz radiation from the interaction of ultrashort laser pulses with gases,” Phys. Rev. E 78, 046406 (2008). [CrossRef]
29. V. B. Gildenburg and N. V. Vvedenskii, “Optical-to-THz Wave Conversion via Excitation of Plasma Oscillations in the Tunneling-Ionization Process,” Phys. Rev. Lett. 98, 245002 (2007). [CrossRef] [PubMed]
30. W.-M. Wang, Z.-M. Sheng, H.-C. Wu, M. Chen, C. Li, J. Zhang, and K. Mima, “Strong terahertz pulse generation by chirped laser pulses in tenuous gases,” Opt. Express 16, 16999–17006 (2008). [CrossRef] [PubMed]
31. I. Babushkin, W. Kuehn, C. Köhler, S. Skupin, L. Bergé, K. Reimann, M. Woerner, J. Herrmann, and T. Elsaesser, “Ultrafast Spatiotemporal Dynamics of Terahertz Generation by Ionizing Two-Color Femtosecond Pulses in Gases,” Phys. Rev. Lett. 105, 053903 (2010). [CrossRef] [PubMed]
32. C. Köhler, E. Cabrera-Granado, I. Babushkin, L. Bergé, J. Herrmann, and S. Skupin, “Directionality of terahertz emission from photoinduced gas plasmas,” Opt. Lett 36, 3166–3168 (2011). [CrossRef] [PubMed]
33. A. Perelomov, V. Popov, and M. Terent’ev, “Ionization of atoms in an alternating electric field,” Sov. Phys. JETP 23, 924–934 (1966).
34. M. Ammosov, N. Delone, and V. Krainov, “Tunnel ionization of complex atoms and of atomic ions in an alternating electromagnetic field,” Sov. Phys. JETP 64, 1191–1194 (1986).
35. G. L. Yudin and M. Y. Ivanov, “Nonadiabatic tunnel ionization: Looking inside a laser cycle,” Phys. Rev. A 64, 013409 (2001). [CrossRef]
36. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover, New York, 1964).
37. B. Nodland and C. J. McKinstrie, “Propagation of a short laser pulse in a plasma,” Phys. Rev. E 56, 7174 (1997). [CrossRef]
38. E. Lefebvre, N. Cochet, S. Fritzler, V. Malka, M.-M. Aléonard, J.-F. Chemin, S. Darbon, L. Disdier, J. Faure, A. Fedotoff, O. Landoas, G. Malka, V. Méot, P. Morel, M. Rabec LeGloahec, A. Rouyer, C. Rubbelynck, V. Tikhonchuk, R. Wrobel, P. Audebert, and C. Rousseaux, “Electron and photon production from relativistic laser plasma interactions,” Nucl. Fusion 43, 629–633 (2003). [CrossRef]
39. R. Nuter, L. Gremillet, E. Lefebvre, A. Lévy, T. Ceccotti, and P. Martin, “Field ionization model implemented in Particle In Cell code and applied to laser-accelerated carbon ions,” Phys. Plasmas 18, 033107 (2011). [CrossRef]
40. I. Babushkin, S. Skupin, and J. Herrmann, “Generation of terahertz radiation from ionizing two-color laser pulses in Ar filled metallic hollow waveguides,” Opt. Express 18, 9658–9663 (2010). [CrossRef] [PubMed]