Hostname: page-component-78c5997874-t5tsf Total loading time: 0 Render date: 2024-10-31T23:17:36.321Z Has data issue: false hasContentIssue false

Terahertz Radiation from a Plasma Cylinder with External Radial Electric and Axial Magnetic Fields

Published online by Cambridge University Press:  01 January 2024

L. H. Cao*
Affiliation:
Institute of Applied Physics and Computational Mathematics, Beijing 100088, China HEDPS, Center for Applied Physics and Technology, Peking University, Beijing 100871, China
Wei Yu
Affiliation:
Shanghai Institute of Optics and Fine Mechanics, Shanghai 201800, China College of Engineering Physics, Shenzhen Technology University, Shenzhen 518118, China
M. Y. Yu
Affiliation:
College of Engineering Physics, Shenzhen Technology University, Shenzhen 518118, China Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027, China Institute for Theoretical Physics I, Ruhr University, D-44780 Bochum, Germany
C. Y. Yu
Affiliation:
State Key Laboratory of Drug Research, Shanghai Institute of Materials Media, Chinese Academy of Sciences, Shanghai 201203, China
*
Correspondence should be addressed to L. H. Cao; cao_lihua@iapcm.ac.cn
Rights & Permissions [Opens in a new window]

Abstract

Terahertz (THz) radiation from a plasma cylinder with embedded radial electric and axial magnetic fields is investigated. The plasma density and the electric and magnetic fields are such that the electron plasma frequency is near the electron cyclotron frequency and in the THz regime. Two-dimensional particle-in-cell simulations show that the plasma electrons oscillate not only in the azimuthal direction but also in the radial direction. Spectral analysis shows that the resulting oscillating current pattern has a clearly defined characteristic frequency near the electron cyclotron frequency, suggesting resonance between the cyclotron and plasma oscillations. The resulting far-field THz radiation in the axial direction is also discussed.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © 2021 L. H. Cao et al.

1. Introduction

The frequency regime of Terahertz (THz) radiation lies between that of microwaves and optical waves. Because of shortage of effective sources, the THz regime was often considered as a gap in the usable electromagnetic-wave spectrum. Generation and manipulation of THz radiation (or T-rays) have been extensively investigated because of its relevance in many areas of basic science and application. In particular, it has potential use in medical and security-detection imaging and communication, as well as in investigating molecular dynamics in matter in the condensed phase.

Many approaches have been proposed for generating T-rays, including synchrotron radiation [Reference Nakazato, Oyamada and Niimura1], optical transition radiation [Reference Happek, Sievers and Blum2, Reference Leemans, Geddes and Faure3], and nonlinear frequency upconversion [Reference Savage, Joshi and Mori4]. Optical transition radiation is produced by relativistic electrons as they cross the interface of different dielectric media. Frequency up- or downconversion of electromagnetic radiation can be realized through nonlinear wave-matter and wave-wave interaction. THz radiation can also be produced in intense laser interaction with matter [Reference Leemans, Geddes and Faure3, Reference Hamster, Sullivan, Gordon, White and Falcone5, Reference Hamster, Sullivan, Gordon and Falcone6Reference Hu and Yi14]. T-rays can also be produced by the oscillating electrons in laser-induced wakefields [Reference Cao, Yu and Xu15Reference Wang, Gibbon, Sheng and Li17], Smith–Purcell effect of the electron beams in laser-matter interactions [Reference Gopal, Woldegeorgis and Herzer18], two-color laser gas-plasma interaction [Reference Magesh Kumar, Kumar, Yuan, Sheng and Chen19, Reference Ahmad, Mahmoud, Purohit and Khan20], transition-Cherenkov effect from laser filaments [Reference Hu, Shen, Lei, Li and Xu21], laser-driven wire-guided helical undulators [Reference Tian, Liu and Bai22], and linear and nonlinear mode conversion [Reference Winner, Schomburg and Brand23], as well as laser interaction with strongly magnetized plasmas [Reference Bakhtiari, Esmaeilzadeh and Ghafary24Reference Purohit, Rawat and Rawat27]. However, most such sources of T-rays are rather complex and inconvenient to realize in practice.

We consider in this paper an unusual but simple source of T-rays, namely, a plasma cylinder with embedded axial magnetic electric and radial fields. Two-dimensional (2D) particle-in-cell (PIC) simulations show that electrons in the plasma cylinder oscillate in the azimuthal as well as radial directions. The evolution of the oscillating current distribution is investigated. It is found that resonance-like behavior between the electron oscillations near the plasma and cyclotron frequencies can take place, and the resulting current pattern can produce T-rays in the axial direction.

This paper is organized as follows. In Section 2, the proposed scheme and simulation parameters are described. Section 3 presents the simulation results at the initial stage, and Section 4 considers the long-time behavior. The far-field power of the T-ray emitted from the plasma cylinder is given in Section 5. Section 6 gives the summary.

2. Proposed Scheme and Simulation Parameters

We consider here an unusual source of THz radiation, namely, from a plasma cylinder with external radial electric field and axial magnetic field, as shown in Figure 1. The electric field producing the plasma could be from a radial capacitor discharge and the magnetic field from a solenoid coil outside the cylinder. The electron motion is governed by

(1) d t γ u = e E + 1 c u × B ,
(2) d t γ = e m c 2 u × E ,

where u is the electron velocity, γ is the relativistic factor, and e, m, and c are the electron charge and mass, and the speed of light, respectively. The dynamics of both electrons and (hydrogen) ions are fully included in the PIC simulations [Reference Xu, Chang, Zhuo, Cao and Yue28]. However, as expected, within the picosecond time of interest here, the electron motion is dominant. The constant embedded electric and magnetic fields are E 0 = E 0 r ^ and B 0 = B 0 z ^ , where E 0 = 3.4 × 1010 V/m and B 0 = 11.4 T, which we note are readily realizable. The cylindrical plasma cylinder is of radius R = 3286 μm and initially of uniform density n 0 = 1.26 × 1020 m−3. The temperature of the electrons and ions is 300 eV. The parameters have been chosen such that the electron cyclotron frequency ω c = e B 0 / m c and plasma frequency ω p e = 4 π e 2 n 0 / m 1 / 2 are both in the THz regime, so that resonance of the electron oscillations can be expected.

Figure 1: The plasma cylinder with external radial electric field E 0 and axial magnetic field B 0, where E 0 = 3.4 × 1010 V/m and B 0 = 11.4 T. The (hydrogen) plasma density is n 0 = 1.26 × 1020 m−3.

To investigate the evolution of such a plasma cylinder, we use the 2D3V (two-dimensional in space and three-dimensional in velocity) PIC simulation code LAPINE [Reference Xu, Chang, Zhuo, Cao and Yue28]. The size of the simulation box is 9418 μm × 9418 μm. The spatial mesh has 1024 × 1024 cells, containing 4 × 107 each of electrons and ions, and the simulation time step is 13 fs. As mentioned, the ion dynamics is fully included, as they might play a role in the long-time behavior of the system.

3. Initial Stage of the Evolution

Figure 2 shows the self-consistently induced electric field ΔEr at t = 13.18 ps and t = 16.11 ps. In this beginning stage of the evolution, besides the cyclotron motion, the plasma electrons are driven radially by external electric field E 0 and acquire a radial velocity υ r. The induced electric field ΔEr, of order hundreds V/m, is much weaker than E 0. The Lorentz force from B 0 and υ r also cause the electrons to move in the azimuthal direction, and the resulting azimuthal velocity in turn causes the electrons to move in the radial direction, leading to the oscillations. The much heavier ions have hardly moved (not shown).

Figure 2: Snapshots of the induced electric field ΔEr at (a) t = 13.18 ps and (b) t = 16.11 ps. Note that the induced electric field is much weaker than the embedded electric field.

The evolution of the current density distribution in the early stage, involving about two oscillation periods, is shown in Figure 3. Figure 3(a) shows that, at t = 11.72 ps, the current density near the cylinder center is larger than that near the boundary. Furthermore, the local current vectors (black arrows) indicate that the stronger current is in the anticlockwise direction (note that the electrons near the center rotate in the clockwise direction). For t = 13.18 ps, Figure 3(b) shows that the outward-moving electrons approach the boundary of the circular column. They are then reflected by the sheath there and move inwards toward the center, as can be seen in Figure 3(c). Roughly, Figures 3(a) and 3(c) correspond to near-minimum amplitudes and Figures 3(b) and 3(d) to near-maximum amplitudes, of the current oscillations. The evolution of the current density in Figure 3 shows that besides rotating in the anticlockwise direction, the electrons oscillate in a rather complicated manner.

Figure 3: Snapshots of the current distribution at (a) t =11.72, (b)t =13.18, (c) t =14.65, and (d) t =16.11 ps. The black arrows represent local current vectors. Close inspection shows that the currents oscillate both radially and azimuthally. The white dashed circles mark the boundary of the initial plasma cylinder.

4. Long-Time Behavior

Figure 3 shows the electron dynamics in the initial two oscillation cycles. To see the behavior at longer times, we now focus on a fixed location inside the plasma cylinder. Figures 4(a) and 4(c) show the evolution of the radial and azimuthal currents at (x, y) = (7346.04 μm, 7346.04 μm). One can see that, during the 80 ps long interval, the electron currents J r and J φ have undergone many oscillation cycles. Figures 4(b) and 4(d) show the corresponding frequency spectra. We can see in both the radial and azimuthal spectra a dominant peak at f∼0.35 THz (or ω ∼ 2.2 rad/sec). In fact, their overall spectral profiles are almost identical, suggesting that a 2D normal-mode-like structure has formed.

Figure 4: Evolution of the (a) radial current density J r and (c) azimuthal current density J φ at (x, y) = (7346.04 μm, 7346.04 μm). (b, d) The corresponding frequency spectra.

Figure 5 shows the evolution of the total current density J x in the simulation box. Figure 5(a) shows that the initial oscillations of the total current density are of large amplitude (to be expected since the initial plasma and fields do not form a self-consistent system), and the corresponding rate of decrease of the envelope is − 4.9 A/m2/ps. At later times, the rate of decrease becomes − 0.0012 A/m2/ps. Figure 5(b) shows the corresponding frequency spectrum, showing that the characteristic frequency of the current oscillations at this location is also ∼0.35 THz.

Figure 5: Evolution of the (a) total current density J x and (b) its frequency spectrum.

5. Terahertz Radiation

Using the temporal and spatial current density from the simulation, one can obtain the far-field radiation. The total electromagnetic energy per unit solid angle Ω per unit frequency is given by the following equation [Reference Cao, Yu and Xu15, Reference Jackson29]:

(3) d 2 I d ω d Ω ω 2 4 π 2 c 3 d t d 3 x n ^ × n ^ × J x , t e i ω t n ^ x / c 2 ,

where n ^ is the direction of the observer with respect to the radiating body. The intensity distribution of the radiation in the z-direction is then [Reference Cao, Yu and Xu15]

(4) d 2 I d ω d Ω ω 2 4 π 2 c 3 d t J t e i ω t 2 ,

where J t = d 3 x J x , t . The corresponding result is shown in Figure 6. We can see that the radiation spectrum has a broad peak centered at f 0∼0.35 THz, as to be expected in view of the current density spectra discussed above. The full-width-at-half-maximum bandwidth Δf/f 0 is about 0.057.

Figure 6: The frequency spectrum of the THz radiation. The peak is at ∼0.35 THz, and the bandwidth is ∼5.7%.

6. Summary

In this paper, we have shown that a bounded plasma cylinder with embedded radial electric and axial magnetic fields can emit T-rays if the parameters of the plasma and the electric and magnetic fields are such that the electron plasma and cyclotron frequencies nearly match and are in the THz regime. As a result, resonance-like behavior of the electrons can occur. The electrons not only rotate in the azimuthal direction but also oscillate in both the radial and azimuthal directions. The resulting currents can emit THz radiation. Finally, it may be of interest to note that, in a nonperturbative analysis of the cold fluid and Maxwell equations, where all the electric and magnetic fields are self-consistent, similar coupled oscillation behavior of the charged particles can exist in an unbounded and self-consistently expanding plasma [Reference Karimov, Stenflo and Yu30].

Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 11875091), the National K&D Program of China (Grant no. 2016YFA0401100), and the Science and Technology on Plasma Physics Laboratory at CAEP.

References

Nakazato, T., Oyamada, M., Niimura, N. et al., “Observation of coherent synchrotron radiation,” Physical Review Letters, vol. 63, p. 124, 1989.10.1103/PhysRevLett.63.1245CrossRefGoogle ScholarPubMed
Happek, U., Sievers, A. J., and Blum, E. B., “Observation of coherent transition radiation,” Physical Review Letters, vol. 67, no. 21, p. 2962, 1991.10.1103/PhysRevLett.67.2962CrossRefGoogle ScholarPubMed
Leemans, W. P., Geddes, C. G. R., Faure, J. et al., “Observation of terahertz emission from a laser-plasma accelerated electron bunch crossing a plasma-vacuum boundary,” Physical Review Letters, vol. 91, Article ID 074802, 2003.10.1103/PhysRevLett.91.074802CrossRefGoogle Scholar
Savage, R. L., Joshi, C. Jr., and Mori, W. B., “Frequency upconversion of electromagnetic radiation upon transmission into an ionization front,” Physical Review Letters, vol. 68, no. 7, p. 946, 1992.10.1103/PhysRevLett.68.946CrossRefGoogle ScholarPubMed
Hamster, H., Sullivan, A., Gordon, S., White, W., and Falcone, R. W., “Subpicosecond, electromagnetic pulses from intense laser-plasma interaction,” Physical Review Letters, vol. 71, no. 17, p. 2725, 1993.10.1103/PhysRevLett.71.2725CrossRefGoogle ScholarPubMed
Hamster, H., Sullivan, A., Gordon, S., and Falcone, R. W., “Short-pulse terahertz radiation from high-intensity-laser-produced plasmas,” Physical Review E, vol. 49, no. 1, p. 671, 1994.10.1103/PhysRevE.49.671CrossRefGoogle ScholarPubMed
Carr, G. L., Martin, M. C., Mckinney, W. R., Jordan, K., Neil, G. R., and Williams, G. P., “High-power terahertz radiation from relativistic electrons,” Nature, vol. 420, no. 6912, p. 153, 2002.10.1038/nature01175CrossRefGoogle ScholarPubMed
Sprangle, P., Peñano, J. R., Hafizid, B., and Kapetanakos, C. A., “Ultrashort laser pulses and electromagnetic pulse generation in air and on dielectric surfaces,” Physical Review E, vol. 69, Article ID 066415, 2004.10.1103/PhysRevE.69.066415CrossRefGoogle ScholarPubMed
Yu, T. P., Ma, Y. Y., Chang, W. W., Yin, Y., Tian, C. L., and Shao, F. Q., “Numerical simulation on effect of laser parameters on terahertz radiation,” High Power Laser and Particle Beams, vol. 20, p. 943, 2008, in Chinese.Google Scholar
Li, Y. T., Li, C., Zhou, M. L. et al., “Strong terahertz radiation from relativistic laser interaction with solid density plasmas,” Applied Physics Letters, vol. 100, Article ID 254101, 2012.10.1063/1.4729874CrossRefGoogle Scholar
Varshney, P., Sajal, V., Chauhan, P., Kumar, R., and Sharma, N. K., “Effects of transverse static electric field on terahertz radiation generation by beating of two transversely modulated Gaussian laser beams in a plasma,” Laser and Particle Beams, vol. 32, p. 375, 2014.10.1017/S026303461400024XCrossRefGoogle Scholar
Liao, G. Q., Li, Y. T., Zhang, Y. H. et al., “Demonstration of coherent terahertz transition radiation from relativistic laser-solid interactions,” Physical Review Letters, vol. 116, Article ID 205003, 2016.10.1103/PhysRevLett.116.205003CrossRefGoogle ScholarPubMed
Liao, G.-Q. and Li, Y.-T., “Review of intense terahertz radiation from relativistic laser-produced plasmas,” IEEE Transactions on Plasma Science, vol. 47, no. 6, p. 3002, 2019.10.1109/TPS.2019.2915624CrossRefGoogle Scholar
Hu, K. and Yi, L., “Relativistic terahertz radiation generated by direct-laser-accelerated electrons from laser-foil interactions,” Physical Review A, vol. 102, Article ID 023530,2020.10.1103/PhysRevA.102.023530CrossRefGoogle Scholar
Cao, L., Yu, W., Xu, H. et al., “Terahertz radiation from oscillating electrons in laser-induced wake fields,” Physical Review E, vol. 70, Article ID 046408, 2004.10.1103/PhysRevE.70.046408CrossRefGoogle ScholarPubMed
Sheng, Z. M., Mima, K., Zhang, J., and Sanuki, H., “Emission of electromagnetic pulses from laser wakefields through linear mode conversion,” Physical Review Letters, vol. 94, Article ID 095003, 2005.10.1103/PhysRevLett.94.095003CrossRefGoogle ScholarPubMed
Wang, W. M., Gibbon, P., Sheng, Z. M., and Li, Y. T., “Emission of electromagnetic pulses from laser wakefields through linear mode conversion,” Physical Review Letters, vol. 114, p. 253901, 2015.10.1103/PhysRevLett.114.253901CrossRefGoogle Scholar
Gopal, A., Woldegeorgis, A. H., Herzer, S. et al., “Smith-Purcell radiation in the terahertz regime using charged particle beams from laser-matter interactions,” Laser and Particle Beams, vol. 34, no. 1, p. 187, 2016.10.1017/S0263034615001093CrossRefGoogle Scholar
Magesh Kumar, K. K., Kumar, M., Yuan, T., Sheng, Z. M., and Chen, M., “Terahertz radiation from plasma filament generated by two-color laser gas-plasma interaction,” Laser and Particle Beams, vol. 33, no. 3, p. 473, 2015.10.1017/S0263034615000518CrossRefGoogle Scholar
Ahmad, N., Mahmoud, S. T., Purohit, G., and Khan, F., “Two color laser self focusing and terahertz generation in multi-ion species plasma,” Optik, vol. 158, p. 1533, 2018.10.1016/j.ijleo.2017.12.152CrossRefGoogle Scholar
Hu, G.-Y., Shen, B., Lei, A.-L., Li, R.-X., and Xu, Z.-Z., “Transition-Cherenkov radiation of terahertz generated by super-luminous ionization front in femtosecond laser filament,” Laser and Particle Beams, vol. 28, no. 3, p. 399, 2010.10.1017/S0263034610000388CrossRefGoogle Scholar
Tian, Y., Liu, J., Bai, Y. et al., “Femtosecond-laser-driven wire-guided helical undulator for intense terahertz radiation,” Nature Photonics, vol. 11, no. 4, p. 242, 2017.10.1038/nphoton.2017.16CrossRefGoogle Scholar
Winner, S., Schomburg, E., Brand, S. et al., “Frequency doubling and tripling of terahertz radiation in a gaas/alas superlattice due to frequency modulation of bloch oscillations,” Applied Physics Letters, vol. 77, p. 1259, 2000.10.1063/1.1290141CrossRefGoogle Scholar
Bakhtiari, F., Esmaeilzadeh, M., and Ghafary, B., “Terahertz radiation with high power and high efficiency in a magnetized plasma,” Physics of Plasmas, vol. 24, Article ID 073112, 2017.10.1063/1.4991395CrossRefGoogle Scholar
Rawat, P., Rawat, V., Gaur, B., and Purohit, G., “Generation of terahertz radiation by intense hollow Gaussian laser beam in magnetised plasma under relativistic-ponderomotive regime,” Physics of Plasmas, vol. 24, Article ID 073113, 2017.10.1063/1.4993055CrossRefGoogle Scholar
Singh, R. K., Singh, M., Rajouria, S. K., and Sharma, R. P., “Strong terahertz emission by optical rectification of shaped laser pulse in transversely magnetized plasma,” Physics of Plasmas, vol. 24, Article ID 073112, 2017.Google Scholar
Purohit, G., Rawat, V., and Rawat, P., “Generation of terahertz radiation from beating of two intense cosh-Gaussian laser beams in magnetized plasma,” Laser and Particle Beams, vol. 37, no. 4, p. 415, 2019.10.1017/S0263034619000685CrossRefGoogle Scholar
Xu, H., Chang, W. W., Zhuo, H. B., Cao, L. H., and Yue, Z. W., “Parallel programming of 2(1/2)-dimensional pic under distributed-memory parallel environments,” Chinese Journal of Physics, vol. 19, p. 305, 2002, in Chinese.Google Scholar
Jackson, J. D., Classical Electrodynamics, Chap. 14, Wiley, New York, NY, USA, 3rd edition, 1999.Google Scholar
Karimov, A. R., Stenflo, L., and Yu, M., “Coupled azimuthal and radial flows and oscillations in a rotating plasma,” Physics of Plasmas, vol. 16, Article ID 062313, 2009.Google Scholar
Figure 0

Figure 1: The plasma cylinder with external radial electric field E0 and axial magnetic field B0, where E0 = 3.4 × 1010 V/m and B0 = 11.4 T. The (hydrogen) plasma density is n0 = 1.26 × 1020 m−3.

Figure 1

Figure 2: Snapshots of the induced electric field ΔEr at (a) t = 13.18 ps and (b) t = 16.11 ps. Note that the induced electric field is much weaker than the embedded electric field.

Figure 2

Figure 3: Snapshots of the current distribution at (a) t =11.72, (b)t =13.18, (c) t =14.65, and (d) t =16.11 ps. The black arrows represent local current vectors. Close inspection shows that the currents oscillate both radially and azimuthally. The white dashed circles mark the boundary of the initial plasma cylinder.

Figure 3

Figure 4: Evolution of the (a) radial current density Jr and (c) azimuthal current density Jφ at (x, y) = (7346.04 μm, 7346.04 μm). (b, d) The corresponding frequency spectra.

Figure 4

Figure 5: Evolution of the (a) total current density ∑Jx and (b) its frequency spectrum.

Figure 5

Figure 6: The frequency spectrum of the THz radiation. The peak is at ∼0.35 THz, and the bandwidth is ∼5.7%.