Multifocal terahertz radiation by intense lasers in rippled plasma

This paper presents a theoretical model for the generation of terahertz radiation by cosh-Gaussian laser beams of high intensity, which are capable of creating relativistic–ponderomotive nonlinearity. We find the components of the terahertz radiation for the relativistic laser plasma interaction, i.e. beating of the two lasers of same amplitude and different frequency in under dense plasma. We plot the electric field profile of the emitted radiation under the effect of lasers index. By creating a dip in peak of the incident lasers’ fields, we can achieve multifocal terahertz radiation.


Introduction
Recently, the high-field THz pulses have been generated from large-scale accelerators using ultrashort electron bunches. For example, peak field of 44 MV/cm has been obtained via coherent transition radiation in the Linac Coherent Light Source [1]. In the laser-based schemes, THz pulses with field strength [8 MV/cm at 1 kHz repetition rate have been generated via two-color laser filamentation [2]. Single-cycle THz pulse with field strength up to 36 MV/cm can be generated by optical rectification of a mid infrared laser in a large-size nonlinear organic crystals assembly [3]. Electromagnetic radiations have been generated by the interaction of ultrashort lasers with gas-jet plasmas [4,5] and water vapor [6]. In addition to these, other methods that have produced high-field THz pulses (but only with components \20 THz) include difference frequency mixing process of two near-infrared lasers in second-order nonlinear crystals [7], and optical rectification (with component limited to \1.5 THz) in lithium niobate (LiNbO 3 ) crystals with tilted laser front [8]. Besides, high-field THz radiations can be generated from relativistic laser irradiated plasmas via various mechanisms [9][10][11][12].
Several experiments employ plasma as a nonlinear medium for the generation of THz radiation by using sub picosecond laser pulses and energetic electron beams [13], since the plasma has the advantage of handling very high field with a remarkable property of being not damaged and having strong nonlinearity [14]. Malik et al. [15] have analytically investigated the THz generation by tunnel ionization of a gas jet with superposed femtosecond laser pulses shone onto it after passing through an axicon. In this mechanism, the temporal evolution of plasma density and hence, oscillatory current density via oscillating dipoles gives rise to the emission of THz radiation. Electron positron plasma has been suggested for high-efficient THz generation by using laser pulse [16]. All the methods suggest that the intense THz radiation can be achieved when very high-intensity lasers are employed in the plasma. However, under that situation, the electron motion will become relativistic and hence, their mass will be varied. This situation is considered in the present article along with the use of coshyperbolic-Gaussian lasers.

Coshyperbolic-Gaussian lasers
We consider linearly polarized cosh-Gaussian electromagnetic beams propagating along the z-axis. Their electric fields are along the y-axis, given by together with j = 1, 2 and p as the index of the lasers. Here E 0 is the amplitude of the fields; b is the beam width of the lasers. Figure 1 shows the profile of the incident laser beam when p C 2 and for different values of the skewness parameter s.

Relativistic ponderomotive force
We use the cold fluid equations for electrons with their number density n 0 and velocity t under the impact of lasers fields. Hence, the electrons oscillatory velocity is obtained as Here we neglect the effect of collisions in the plasma because collision effect is not significant for high-intensity laser beam. The density n obeys the continuity equation on ot þ r Á nt ð Þ ¼ 0. In view of the high-intensity lasers, we consider relativistic motion of the electrons. The relativistic ponderomotive force in the presence of an intense electromagnetic beam can be represented [17,18] as Here m is the rest mass of the electron, c is the velocity of the light, and t 1 and t 2 are the oscillatory drift velocities of the electrons. The use of Eq. (1) in Eq. (2) yields the ponderomotive force as We obtain the following expression for the ponderomotive force: This nonlinear force causes perturbations in the density of the plasma. These perturbations generate space charge potential that leads to the linear density perturbations. The nonlinear density perturbations are obtained as Finally, the nonlinear velocity of the electrons is obtained from equation of motion as Here, x p is the plasma frequency given by q .
The nonlinear current density can be obtained from Here n 2 is the ripple density in the plasma given as n 2 ¼ n a e iaz together with a as the wave number of the ripples produced in the plasma and n a as the amplitude of the density ripples. These density ripples are generally used to achieve the resonance condition. These ripples in density may be produced using various techniques, which involve transmissive ring grating and a patterned mask where the control of ripple parameters might be possible by changing the groove structure, groove period, and duty cycle in such a grating and by adjusting the period and size of the mask [19][20][21][22][23]. Finally, the nonlinear current density is obtained as Calculation of electric field of terahertz radiation The wave equation that governs the emitted THz radiation is given by Here Ẽ THz is the electric field associated with the THz radiation and e is dielectric constant of the plasma amounting to 1 À Here C and D are frequency-dependent terms defined as In order to find the amplitude of the THz radiation, we analyze the expression (9). This relation can be bifurcated into two different expressions as shown in Eq. (10), indicating that it is superposition of two different waves. This is the consequence of the different components of the ponderomotive force. Finally, we achieve the following expression for E THz : Here J and R are the amplitudes of the different wave components that are generated due to the relativistic ponderomotive force nonlinearity in the plasma. The symbols J, R and G are defined as J ¼ ÀGCD; For the high-intensity laser, where a part of its energy is used in heating of the electrons in the plasma, the resonance condition is expected to be modified as where t th is the thermal velocity of the electrons. However, the second term would play a significant role only when k is very large, i.e. for the short wavelength laser pulses. So it is clear that the resonance condition is departed from x = x p when one considers the heating in the plasma and under this situation reduced THz field would be realized (clear from Fig. 3, later).

Results and discussion
This section deals with the results obtained from Eq. (10) for the electric field of THz radiation. The following set of parameters has been used in the numerical analysis: E 0 ¼ 2:7 Â 10 12 V=m; b ¼ 0:01 Â 10 À2 m; s ¼ 0:6; x ¼ 1:15x p ; x 1 ¼ 1:65 Â 10 15 rad/s and x p = 2.0910 13 rad/s. This is seen that out of two components of fields, given by Eq. (10), only one component (second term) dominates and is solely responsible for the THz amplitude. Figure 1 shows the electric field distribution of the incident lasers. This is clear that a dip occurs in the peak of the fields when p [ 2. As the skewness parameter s is increased this dip takes a prominent form in the incident laser profile. Figure 2 shows the output profile (normalized amplitude) of the THz radiation with the normalized distance, when n a = 0.3n 0 . The maximum amplitude occurs at particular values of y/b. Peak positions are symmetric in nature because of skewness parameter s. Maximum amplitudes correspond to the higher SG index lasers. A comparison of Figs. 1 and 2 shows that in general two peaks of the emitted THz radiation are obtained for the lasers having one peak. However, when there is a dip in the peak value of the lasers field, small peaks start generating in the produced radiations. Also the higher peak is further enhanced for the larger dip in the incident laser profile. This is true for the lasers having higher skewness in their profiles (s [ 0.6).
The dependence of THz field on the ripple amplitude n a and beating frequency x is shown in Fig. 3. It is evident that the amplitude E THz is high when n a carries higher values, but it decreases sharply when the resonance condition x * x p is departed. Ripples in plasma density play a positive role to enhance the THz amplitude because large Ac 2 x 1 x 2 À 1 2 3 4 numbers of electrons participate for the generation of nonlinear current so amplitude of THz field increases accordingly.
In order to find the effect of beam width of the lasers on the THz radiation, we have plotted the intensity (normalized by E 0 j j 2 , i.e. the intensity of lasers) with b and p in Fig. 4. Here, we clearly observe that the lasers with lower beam width produce intense radiation and the same is the case with higher index lasers (p). These results are in agreement with the observations made by other researchers [24][25][26].

Comparative study
Kumar et al. [27] have studied the Gaussian lasers beam beating for THz generation. Our results are comparable to their observations when we consider the relativistic case Variation of THz radiation normalized field with beating frequency x and density ripple amplitude n a , for the same parameters as in Fig. 2 and use coshyperbolic-Gaussian lasers. Flat top or super-Gaussian like lasers destabilize as these propagate in the medium leading to two separate peaks in THz profile, whereas Gaussian-like beam retains its shape as it propagates through the medium. Therefore, the THz field in Fig. 2 shows small intensity peak at the center. In nonrelativistic cases the SG index plays important role to enhance the amplitude of the THz [26], but in case of relativistic ponderomotive force the density ripples are found to enhance the amplitude significantly, whereas SG index does not play that much significant role. Main difference between the results obtained in relativistic and nonrelativistic cases is that multiple components appear in Eq. [10] as obtained by Bituk and Fedorov [28] also.

Conclusions
When we consider the relativistic effect in plasma due to the higher intensity lasers, the emitted THz radiation is not found to be focused at a single place, rather it is focused at least at two places. These amplitudes carry different magnitudes. We can create multiple focal THz radiation by creating a dip in the electric field peak of the beating lasers. Other effects of the beam width and density ripples are found to be the similar as in nonrelativistic case.