Efficiency Analysis of Terahertz Radiation Produced under Different Pumping Conditions of Two Femtosecond Laser Pulses ()
1. Introduction
The terahertz (THz) vortex beam is a light beam with a spiral wave front [1]. Because of its distinctive ability to carry an Orbital Angular Momentum (OAM) that is attributed to the correlated spatial phase distribution [2], the vortex THz beam is running in various real-life applications, ranging from the THz wireless communication [3] to the detection of the astrophysical images [4], in addition to the super-resolution imaging [5], control the chirality of twisted metal nanostructures [6], and the electron bunch acceleration [7]. Basically, the production of THz vortex beams is maintained by two main mechanisms which are: 1) the wavefront modulation device mechanism; in which the vortex phase plates [8] and the THz hologram technology [9] [10] method are applied, and 2) the direct excitation of the THz vortex helicity mechanism; where the optical rectification (OR) [11], the difference frequency generation (DFG) [12] [13], and the spatially periodically and non-periodically modulation plasma [14] [15] method are employed. Among these mentioned methods, peculiar attention has been given to the spatially non-periodically modulation plasma method for its advantage to achieve a high intensity and broadband regular THz field behind the threshold damage of bulk materials [15]-[17].
In the non-periodically modulated plasma method, by the filamentation of a fundamental femtosecond beam (
) and its second harmonic (
) in air plasma an ultra-broad band THz angularly accelerated vortex beam is induced [18]. The THz beam production by this filamentation is govern by two models, namely the Fourth-Wave Mixing (FWM) [19] model; where the mixing process
or
occurs, and the Photocurrent (PC) model [20]; in which by an asymmetric coupling between the two fs beams an oscillating optical current
is generated to emit a terahertz radiation. Based on the FWM and PC model, extensive researches have been conducted to study the production of the THz vortex beam at different pumping conditions of the two input fs beams. The central purpose of these researches is to examine the properties of the produced THz vortex radiation in order to underline the optimum input beams conditions for an efficient THz vortex beam production. At first, by adjusting the relative phase between two-color vortex laser beams [21], a controllable THz necklace-shaped vortex beam is demonstrated. Afterwards, using few cycle input vortex beams, various MIR-infrared vortex beams are generated [22], as at the low-input-frequency a MIR vortex beam with a simple
stepwise phase invariance with the azimuthal angle (
) is induced, where at the high-input frequency an anisotropic MIR vortex beam that carrying OAM with nonlinearly invariance phase profile is produced. Then after, with a first harmonic Gaussian (
) and a second harmonic Laguerre-Gaussian beam (
) that is carrying a vortex topical charge
with an intensity modulation feature along the azimuthal angle [23], a unique THz vortex beam structure is formed.
In advanced investigations, by tunning the chirping parameters of two few-cycle input fs beams in different input intensity regions [24], different vortex beam structures with diverse properties are established, for instance in low-input intensities where the Kerr effect (FWM) dominates over the plasma effect (PC) a THz vortex beam of uniform patterns with a constant angular velocity is formed, in medium-input intensities where the Kerr and Plasma effects are comparable an Angular Accelerated Vortex Beams (AAVBs) with linear phase-azimuthal angle dependence is generated, and in high-input intensities where the plasma effect is dominated an AAVBs with strongly nonlinear azimuthal angle dependence on the relative phase profile is induced. Moreover, by the combination of a fundamental Laguerre-Gaussian fs vortex beam (FH) with a second harmonic Gaussian beam (SH) at variable relative amplitudes and conjugated topological [25] (TC)
, mm-scale length necklace THz beam with stepwise phase profile is generated, where at large scale of this combination both of AAVBs with linear and nonlinear phase profile can be obtained. In this article, we analyze the efficiency of produced THz vortex beams at different pumping conditions of two fs laser pulses, these conditions are the equal and unequal two fs beams amplitudes filamentation and the different input intensity regions for each filamentation case. The essential objective of this analysis is to determine the optimum pumping conditions for an efficient THz vortex beam production. In section 2, we present the physical model employed to conduct this study by listing its basic equations and explaining its physical assumptions. In section 3, we present for each of the equal and unequal beams amplitudes filamentation case, simulation results for the transverse intensity distributions and the relative-azimuthal angle relation of an induced THz vortex beam in three different input intensity regions.
2. The Physical Mode
In our efficiency analysis for the THz vortex beam induced under different pumping conditions of two input fs beams, the fundamental beam (FH) is a Gaussian (nonvortex) beam with wavelength
(
) and the second harmonic beam (SH) is a Laguerre-Gaussian vortex beam with
(
). The initial scalar combined field of the conventional vortex of these beams is given by
(1)
where
is the spatial part,
and
is the temporal part,
where
are the initial real amplitudes,
are the initial pulse durations of the FH and SH beams, respectively,
is the fundamental frequency,
is the relative phase between the FH and SH.
is the topological part,
where
is the azimuthal angle and
is the topological charge (TC).
In our numerical calculations [26], both the FH and SH beam have the same initial pulse duration
fs, and for the nonvortex FH beam
and for the SH vortex beam
. Thus, the initial scalar combined field can be finally written as:
(2)
Our air plasma is the Nitrogen gas
at initial density
that is Tunneling Ionizing (TI) by the input femtosecond beams. Herein, the ionization rate is the standard ADK formula [27].
(3)
is the atomic frequency,
are the ionization potential of Hydrogen and Nitrogen molecules, respectively.
is the optical field strength in atomic unit, where
is the electric field of the applied fs beams and
is the atomic unit of the electric field.
3. Simulation Results
The filamentation of two fs beams in air plasma is the principle mechanism of our terahertz radiation production [18]. As known, this filamentation is demonstrated when the self-focusing by the Kerr nonlinearly and the defocussing by the nonlinear plasma ionization effects is balance. Essentially, this balance is achieved at the namely known the clamping intensity [28], the clamping intensity (
) is the peak intensity that is constricted inside the plasma filament when the length of the filament starts to increase. The clamping intensity value is the propagation conditions and the plasma structures dependence, for example: in case of an applied polarized beams,
is the simulated value [28] at
, while
is the proper value at
under the influences of the freeman resonance [29]. In a similar to our filamentation condition and air plasma structure, the
is input intensity period [30] where the clamping intensity is effective and the THz radiation is generated. In a correlation between this effective clamping intensity period and the Kerr and plasma nonlinearity effects, in the THz vortex beams production researches this period is classified into: 1) the low input intensity region (
) where the input intensity is less than the clamping intensity and the Kerr nonlinearity effects is predominated, 2) the medium-input intensity region wherein the input intensity nearly equal the clamping intensity (
) and the Kerr and plasma nonlinearity effects are comparative and competitive, and 3) the high-input intensity region (
) at which the input intensity larger than the clamping intensity and the plasma nonlinearity is the dominated effects. In this research, we are crucially interested to conduct an elaborated analysis for the spatial and temporal properties of the produced THz vortex beams for the equal and unequal input beams amplitudes filamentation, in each of these listed input intensity regions.
3.1. The THz Vortex Beams Production by Two fs Beams of Equal Amplitudes
In Figure 1(a) we present the transverse intensity distribution of a THz vortex beam induced by two fs beams with equal amplitudes in the low-input intensity region. As shown in this distribution, a doughnut shape intensity patterns (ring-shaped) with zero intensity in the center is demonstrated. As clearly seen in this demonstration, the intensity patterns are spatiality symmetrically distributed and shaped out of a main ring that is largely spreed over the transverse space and preserves the major non-distortion intensity of the induced THz vortex beam. In addition to this, the main ring is surrounding with few numbers of rings with much larger width and lower intensity. In fact, the ring-shaped THz vortex beam is previously observed in a number of studies, in one of this study [24] that has been conducted at the similar pumping conditions of our study, ring patterns with zero central intensity is also formed, while in another study [25] where the FH beam is LG and the SH is Gaussian, ring patterns with non-zero intensity center is produced, the SH (
) Gaussian input beam is the responsible of this non-zero intensity center as it is explained in this study.
![]()
![]()
Figure 1. The transverse intensity distribution (left) and the phase azimuthal angle relation (right) of a THz vortex beam induced by the filamentation of two fs beams with equal amplitudes in the low-input intensity region.
To analyze the temporal properties of the induced ring-shaped THz vortex beam, we display the phase
and the azimuthal angle
relation of this shape in Figure 1(b). As seen in this figure, the
is linearly evolving with
, the linear
relation indicates that
and the phase-azimuthal angle is non-modulated dependence which implies that the demonstrated ring-shaped patterns displayed in Figure 1(a) has no Orbital Angular Momentum (OAM). It is important to bear in mind that, although the Kerr nonlinearity is predominated in the low-input intensity region, the refractive index variation
value due to this nonlinearity; where
is the nonlinear refractive index in air and
is peak intensity, is not adequate to persuade
variation and modulate the
dependence in this region, thus
is linearly (
) evolving with
as displayed in Figure 1(b), accordingly spatiality symmetric non-rotating ring-shaped THz vortex beam with the doughnut intensity patterns is demonstrated in Figure 1(a).
In Figure 2(a), we re-present the transverse intensity distribution of a THz beam induced by two equal fs beams amplitudes in the medium-input intensity region, a distinct development in the induced THz vortex beam structure is exploited in this figure. As illustrated in this distribution and in comparison with Figure 1(a), the ring-shaped is modified into two-petals-shaped intensity patterns, the modified patterns are spatiality asymmetrically distributed over a circular ring that has approximately the same radius of the main ring displayed in Figure 1(a), as well as, each petal preserves a higher intensity that is confined in a narrower transverse space. To interpret this modified THz vortex beam structure, it is necessary to clarify that in the medium-input intensity region the Kerr and plasma nonlinearity effects are competitive and comparative, other than that each
value due to the Kerr nonlinearity and
value due to the plasma nonlinearity; where
is the effective volume of the plasma [31], is high enough [24], on that account the impacts of the total refractive variation
on the phase variation [31] (
) are highly considered. Within this interpretation and because of the direct dependence of the generated THz beam amplitude
on the resulted phase variation
, the preliminary impact of this dependence is the developed spatiality asymmetric two-petals-shaped intensity patterns structure formed in Figure 2(a).
![]()
![]()
Figure 2. The transverse intensity distribution (left) and the phase azimuthal angle relation (right) of a THz vortex beam induced by the filamentation of two fs beams with equal amplitudes in the medium-input intensity region.
The additional impact of the resulted
variation on the temporal properties of the induced two-petals-shaped THz vortex beam is analyzed in Figure 2(b) where the phase-azimuthal angle relation of a THz beam induced by two equal fs beams amplitudes in the medium-input intensity region is display. As noted in this figure, the phase
is nonlinearly evolving with
, i.e.,
, within this context, the phase is azimuthally modulated which implies that the two-petal patterns exposed in Figure 1(a) has an Orbital Angular Momentum (OAM). Thereupon, the direct impacts of the phase variation
resulted by the total refractive variation
due to the Kerr and plasma nonlinearity effects in the medium-input intensity region is the developed rotated anisotropic two-petals-shaped THz vortex beam demonstrated in Figure 2(a).
It has been recognized from the above results that depending on the input intensity region, the phase variation
resulted due to the influential effects of the Kerr and plasma nonlinearity and its relation with the azimuthal angle
is the origin of developing the THz vortex beam structure and modifying its spatial-temporal properties. In the high-input intensity region, because
[24], the kerr and plasma nonlinearities effects are no-more comparative and the plasma nonlinearity is the dominated effects. Be noted that, the plasma nonlinearity has a distinguished characteristic which is the tunneling ionization is achieving at the peak intensity of the input beams, upon that the ultimate energy is delivered to the air plasma at that moment. The direct impact of this feature is clearly observed in Figure 3(a), as noted the intensity patterns are to a great extent spatially asymmetrically distributed over a narrower ring, moreover the intensity of these patterns is even more reinforced where its maximum values are confined in narrower spatial space.
![]()
![]()
Figure 3. The transverse intensity distribution (left) and the phase azimuthal angle relation (right) of a THz vortex beam induced by the filamentation of two fs beams with equal amplitudes in the high-input intensity region.
To explore the impact of the plasma nonlinearly domination in the high-input region and its peculiar feature on the temporal properties of the induced two-petals-shaped, in Figure 3(b) the phase as function of the azimuthal angle of a THz vortex beam induced by two equal amplitudes fs beam in the high-input intensity region is shown. As clear in this figure and in comparison with Figure 2(b), the phase
is strongly nonlinearly evolving with
, this rigorous nonlinearity leads to the highly modulated
dependence which escalates the OAM acceleration. Thus, in the high-input intensity region, more efficient two-petal-shaped THz vortex beam with more enhancement intensity and accelerated rotation is established as shown in Figure 3(a).
3.2. The THz Vortex Beams Production by Two fs Beams of
Unequal Amplitudes
The filamentaion of two unequal fs beams amplitudes in air plasma is a nonlinear dynamics of a crucial importance in the long-range fs beams propagating [32]. Numerous processes, such as spontaneous symmetry breaking [33] and the cross-correlation technique [34] are taking place during this filamentation, even so, the energy exchange between the two unequal fs beams amplitudes is the most vital among these processes for its boarded applications in supercontinum generation [35], lightning and discharge triggering [36], remote sensing [37], and THZ radiation production [38]. The energy exchange is a promising nonlinear process in controlling the multi fs beams propagation, in any case, the exchange rate and direction depend on the relative time delay, the initial chirps, the input intensities, the intersection angle, and the relative polarization of the two fs beams. The plasma mediated forward Simulated Raman Scattering (SRS) [39], the Traveling Plasma Grating (TPG) [40]-[42], and the Two Beams Coupling (TBC) are the primitive mechanisms to govern the energy exchange. In our study, the TBC is the employed mechanism since it is appropriate whenever a frequency difference between the pumping beams is applied. In this section, we analyze the efficiency
of the produced THz vortex beams induced by the filamentaion of two unequal fs beams amplitudes, the role of emerged energy exchange process on the efficiency and the spatial-temporal properties of the produced THz vortex beam in the three input intensity regions will be investigated.
In Figure 4(a), we present the transverse intensity distribution of a THz beam induced by the superposition of two fs beams with unequal amplitudes in the low-input intensity region. As depicted in this figure, a necklace-shaped THz vortex beam with four identical petals is demonstrated. It is noteworthy that in the THz vortex beams researches, the petal-like-shaped and ring-shaped are the predominated and the most excepted structures, even so, these structures have different and separated intensity patterns distribution. In contrast, the petal-like-shaped patterns are centrally located, while the ring-shape patterns that have lower intensity are spreading over a larger area. Due to the overlapping of these two structures [43] in the presence of the energy exchange in the low-input intensity region, spatiality symmetric necklace-shaped THz vortex beam with four identical petals is demonstrated as seen in Figure 4(a).
The necklace-shaped THz vortex beam is an interesting light beam structure that has a growing applications ranging from cold atom manipulation to optically driven micromachines [44] [45]. Overall, the necklace THz vortex beam with and without OAM has been theoretically and experimentally predicated [46] [47], on that account to examine the temporal properties of our induced necklace THz vortex beam, the phase and the azimuthal angle relation for the two fs beam with unequal amplitudes in the low-intensity region is displayed in Figure 4(b). As noted in this figure, a
-stepwise phase profile is observed for each of the two lobes that are separated at
. The observed profile marks that this necklace-shaped beam has no OAM, in addition the abruptly change of this profile at
implies the reverse polarity
of the induced THz beam phase for each lobe. The non-rotating spatially symmetric intensity patterns is anticipated in the low-input intensity region, as in this region, although the Kerr effect is prevailed,
value is small to develop
variation, furthermore the energy exchange magnitude and rate is not sufficient enough to produce a phase variation and then an azimuthal angle modulation. As a result, the non-accelerated spatially symmetric necklace-shaped THz vortex beam with four identical petals that have equal maximum intensity is demonstrated in Figure 4(a).
![]()
![]()
Figure 4. The transverse intensity distribution (left) and the phase azimuthal angle relation (right) of a THz vortex beam induced by the filamentation of two fs beams with unequal amplitudes in the low-input intensity region.
In Figure 5(a), we re-present the transverse intensity distribution of a THz vortex beam generated due to the superposition of two fs beams with unequal amplitudes in the medium-input intensity region. Although a necklace-shaped THz vortex beam with four-petals intensity patterns is remained in this figure, in comparison with Figure 4(a), a notable reformation in respect to the spatial properties of these patterns is emerged. As clear in Figure 5(a), the four-petals are spatially asymmetrically distributed over narrower ring, where each petal is randomly displaced and preserves more refinement maximum intensity that is confined in a narrower space. As explained before, the
resulted by to the total refractive variation
due the kerr and plasma nonlinearity in the medium-input intensity region is origin of this spatiality asymmetric reshaping, meanwhile the presence of the energy exchange, additional nonlinear effects should be regarded which is the Cross Phase Modulation [48] (XPM). The XPM is a nonlinear temporal-spatial modulation that is established by the overlapping between the two unequal fs beams amplitudes, during this overlapping the phase of one beam modulates the counterpart one as these beams propagating through the nonlinear medium. In principle, the XPM induces arbitrarily velocity, peak beam intensity, and relative phase modulation
through the refractive index variation
. As matter of fact, XPM that is the propagating beam intensity dependence
, therefore it is considered as part of kerr nonlinearity. For this reason, in the present study the XPM nonlinearity is summed with the conventional Kerr nonlinearity effect, accordingly, the refractive index variation due to the kerr effects is updated by
.
![]()
![]()
Figure 5. The transverse intensity distribution (left) and the phase azimuthal angle relation (right) of a THz vortex beam induced by the filamentation of two fs beams with unequal amplitudes in the medium-input intensity region.
After we have affirmed the spatial asymmetric properties of the induced necklace-shaped THz vortex beams in the medium-input intensity region, the temporal properties of this shape in this region should be examined as well. In Figure 5(b) we display phase-azimuthal angle relation of a THz beam induced due to the superposition of two fs beams with unequal amplitudes in the low-input intensity region. Although the phase profile is jumping up at the separated lobe point
, in this figure, per each lobe the
is nonlinear evolving with
, on that
is not constant, and thus this induced reshaped has OAM. Thereby, the necklace-shape THz vortex beam formed in Figure 5(a) is spatially anisotropic and angularly momentum accelerating.
As we have previously explained, in the high-input intensity region the plasma nonlinearity is the dominated effect and the phase variation
is mostly resulted due to
. For the two unequal beams amplitudes filamentation case, in the high-input intensity region, due to plasma nonlinearity domination and the presence of the energy exchange, the resulted
has straightforward impact on the oscillation optical current
that consequently influences on the induced patterns intensity (
) through the intensity transport equation [49]
(4)
where
patterns intensity of the energy exchange during the propagation. To look into the consequences of the updated
impact and influence on the spatial properties of the induced necklace-shaped THz vortex beam, we display in Figure 6(a) the transverse intensity distribution of an induced THz vortex beam by the superposition of two fs beams with unequal amplitudes in the high-input intensity region. In comparison with Figure 5(a), the spatial properties of the induced necklace THz vortex beam is further developing, as represented, the four-petals of the remained necklace-shaped are highly spatially asymmetrically distributed over more narrower ring, beside each petal is largely randomly displaced and preserves extra refinement maximum intensity that is confined in more tight space.
![]()
![]()
Figure 6. The transverse intensity distribution (left) and the phase azimuthal angle relation (right) of a THz vortex beam induced by the filamentation of two fs beams with unequal amplitudes in the high-input intensity region.
To follow up the impact and influence of the updated
on the temporal properties on this induced beam, in Figure 6(a) we show phase-azimuthal angle relation of two fs beams with unequal amplitudes in the high-input intensity region. As illustrated in this figure and in a comparison with Figure 5(a), per each lobe, the
is more nonlinearly evolving, hence more efficient and highly accelerated THz vortex necklace-shaped THz vortex beam with largely anisotropic petals distribution that preserve more confined maximum intensity for each petal is demonstrated in Figure 6(a).
4. Conclusion
We analyzed the spatial-temporal properties of THz vortex beams induced under an equal and unequal two fs beams amplitudes filamentation conditions, and in three different input intensity regions for each condition. Numerical simulations have revealed that in the equal beams amplitudes filamentation case, a spatiality symmetric non-rotated ring-shaped THz vortex beam is demonstrated in the low-input intensity region, while in the medium-input intensity region where the Kerr and plasma nonlinearity effects start to be influential, the ring-shaped THz vortex beam is modified into higher intensity anisotropic two-petals-shaped that is associated with OAM, meanwhile in the high-input intensity region where the plasma nonlinearly is the main key player, the properties of this modified anisotropic beam is more enhanced to preserve more higher intensity at more accelerated OAM. In the unequal beams amplitudes filamentation case and in the presence of the energy exchange between these amplitudes, in the low-input intensity region spatiality symmetric non-rotated necklace-shaped THz vortex beam with four identical petals is demonstrated, even though in the medium-input intensity region where the updated Kerr and plasma nonlinearity effects are comparatively influential, the necklace-shaped is remained but with spatiality anisotropic non identical four-petals-shaped at larger intensity and angular momentum acceleration, on the other hand in the high-input intensity region where the plasma nonlinearly is dominated, more refinement necklace-shaped THz vortex beam is formed at more confined and magnified intensity with an increasing angular momentum acceleration. Regardless of the filamentation condition and the demonstrated THz vortex beam shaped, in the high-input intensity region and due to the dominated plasma nonlinearity effects, more efficient THz vortex beam is produced.