Propagation characteristics of oblique incidence terahertz wave through non-uniform plasma
Chen Antao, Sun Haoyu, Han Yiping, Wang Jiajie, Cui Zhiwei
School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China

 

† Corresponding author. E-mail: yphan@xidian.edu.cn wangjiajie@xidian.edu.cn

Project supported by the National Basic Research Program of China (Grant No. 2014CB340203), the National Natural Science Foundation of China (Grant Nos. 61431010 and 61501350), and the Natural Science Foundation of Shaanxi Province, China (Grant Nos. 2018JM6016 and 2016JM1001).

Abstract

The propagation characteristics of oblique incidence terahertz (THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain (SO-FDTD) method combined with the phase matching condition. The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab. Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions, which makes it a potential candidate to overcome the ionization blackout problem.

1. Introduction

When vehicles fly at hypersonic velocities within the atmosphere, they become enveloped in a plasma sheath that prevents radio communication between the aerospace and the ground control center, which is the well-known blackout problem.[14] In order to mitigate or eliminate the reentry blackout problem, a number of techniques have been proposed, such as high power, high frequencies, low frequencies, lasers, magnetic fields, and so on.[3,57] The terahertz wave (THz), commonly referred to as T-rays, sub-millimeter, or far infrared, refers to the electromagnetic radiation at the frequency from 0.1 THz to 10 THz, has attracted much attention recently partially because it is one of the potential candidates in the high frequencies communication methods. The techniques for the generation of THz wave were developed rapidly in recent years[811] and promising applications of THz waves can be found in various fields.[12] As the basis of communication between the hypersonic vehicles and the ground control center, the interaction between plasma and THz wave (absorption, reflection, and transmission) is a key issue to be analyzed. Over the past few years, the propagation characteristics of the THz wave in plasma have been investigated by different methodologies, such as the analytical method,[1317] experimental method,[18] scattering matrix method (SMM),[19,20] propagation matrix method (PMM),[21] finite-difference time-domain (FDTD) method,[2228] Wentzel–Kramer–Brillouin (WKB) method,[29,30] and so on.

The FDTD method has been widely used to solve the interaction problems between the THz wave and plasma because of its capability in simulating the propagation characteristics of a wideband THz wave through inhomogeneous dispersive plasma. Wang et al.[22,27,28] investigated the propagation properties of the THz wave in a time-varying dusty plasma slab using the auxiliary differential equation finite-difference time-domain (ADE-FDTD) method. Zheng et al.[23] studied the propagation characteristics of the THz wave in unmagnetized plasma by the Z transform finite-difference time-domain (ZT-FDTD) method. Theoretical results were compared with experimental measurements in Ref. [20], where good agreements were achieved. The broadband propagation characteristics of the THz wave in an anisotropic magnetized plasma were investigated by Tian et al.[24] using the JE convolution FDTD (JEC-FDTD) method. Gao et al.[25] investigated the characteristics of reflection, transmission, and absorption of THz waves in a sandwich type microplasma structure by the ZT-FDTD method.

Although studies of propagation characteristics of the THz wave in plasma media can be found in several works, most of the previous research focused on the normally incident wave cases, i.e., the incident plane wave propagates at an angle of zero with respect to the normal of the media interfaces, limited attention was paid to the oblique incidence cases. In practice, the angle of incidence changes over a wide range during flight, which indicates that the oblique incidence wave cases require more consideration. Recently, using the impedance matching method, Rahmani et al.[31] investigated the transmission, reflection, and absorption of s-polarized wave obliquely incident on the plasma slab which has a bell-like electron density distribution. The SMM was adopted by Yuan et al.[32] to analyze the transmission of obliquely incident THz waves in the plasma sheaths covering a blunt coned vehicle, in which the plasma sheaths were generated based on a hypersonic fluid model. Cao et al.[33] studied the propagation characteristics of obliquely incident THz waves though dusty plasma using the PMM, where the electron density distribution of the dusty plasma was assumed to be parabolic. Based on the twice reflection model, Tian et al.[34] investigated the obliquely incident wave through a nonuniform magnetized plasma slab, in which both the electron density and the collision frequency across the plasma were assumed to be in a Gaussian profile. To solve the problem under study, the non-uniform plasma was divided into a series of layers with uniform electron density.[34] Nevertheless, this dividing method fails when one encounters plasma with a high gradient electron density distribution. To overcome this difficulty, the FDTD is a good choice.

The THz wave obliquely incident on a plasma sheath, which is non-uniform in half space, is a two-dimensional problem for the FDTD method. Thus it will cost a lot of memory and computation time in the simulation. To save computational memory and time, in this paper, the FDTD method combined with the phase matching condition is utilized to solve the propagation problem of obliquely incident THz waves through an inhomogeneous plasma. With the help of phase matching condition, the problem is simplified to a pseudo one-dimensional problem, which can greatly reduce the memory required and improve the efficiency. To describe the distribution of the electron density in an inhomogeneous plasma, the Gaussian profile is introduced which corresponds to the flow field around the hypersonic aircraft at a high altitude of space.[35]

The rest of this paper is structured as follows. In Section 2, the physical model of the propagation problem is given, with a brief description of the method used in solving the problem under study. Numerical results and discussions concerning the transmission features of the THz wave through an inhomogeneous plasma are presented in Section 3. Section 4 concludes this paper.

2. Physical model and methods

A geometric sketch of the problem under study is shown in Fig. 1, where a plane THz wave propagates through a plasma sheath around the hypersonic vehicle at an arbitrary incident angle. The plasma sheath is non-uniform in half space, which is inhomogeneous in y direction and uniform in x and z directions. The thickness of the plasma sheath is d in y direction and infinite in the other two directions. The incident THz wave is assumed to propagate through the plasma sheath in the plane of xoy at an incident angle θ with respect to the normal of the plasma sheath interfaces.

Fig. 1. Schematic diagram of a THz wave propagating through nonuniform plasma layer.

As shown in Fig. 1, when the THz wave propagates through the plasma slab in the xoy plane, the electric field E and magnetic field H satisfy E/∂z = 0 and H/∂z = 0, in which /∂z means derivative to z. The fields in the plasma sheath can be divided to two modes: (i) the transverse magnetic (TMz) mode, which is composed of Ez, Hx, Hy, and (ii) the transverse electric (TEz) mode, which is composed of Hz, Ex, Ey. In the following, we will only consider the TEz mode case since the TMz case can be obtained in the same manner.

The fields in plasma can be described by Maxwell’s equations in the frequency domain

where D is the electric flux density, ε0 is the vacuum permittivity, εr is the relative permittivity of the medium, μ0 is the magnetic permeability of vacuum, and ω is the angular frequency of the incident wave. For the medium considered in this paper, i.e., nonmagnetic plasma, the relative magnetic permeability is μr = 1.

For the TEz mode, the propagation equation in the xoy plane can be expressed as

where H0 is the amplitude of the magnetic field, k is the wave number in plasma, and r = x · x + y · y. The time factor exp(jωt) is omitted and the magnetic field component of the incident plane wave in plasma is expressed as
Based on the transformations /∂x → jkx and /∂t → jω, the derivative of Eq. (2) to x in plasma is
Based on the phase matching condition, the phase velocity in the direction parallel to the media interfaces is constant that kx = kair sin θ. Substituting Eq. (6) into Eq. (3) yields

For simplicity, two new auxiliary variables and are defined as

Substituting Eqs. (8) and (9) into Eq. (7) yields
Equations (1) and (10) are then transformed into time domain and written as
where * represents convolution.

The permittivity of unmagnetized plasma is given as[36]

where ε is the infinite-frequency limit of the relative permittivity; ωp = 2πfp, fp is the plasma frequency, and υc is the electron collision frequency.

Substituting Eq. (13) into Eq. (9) yields

The electric flux density Dx and the auxiliary variable can be written as
Comparing Eqs. (13) and (14) with Eqs. (15) and (16), we can obtain that M = N = 2 and

With the transition relationship from the frequency domain to time domain, jω/∂t, we can obtain the expressive forms of Eqs. (15) and (16) in the time domain

If we assume that the function in the time domain has the form

The central difference of Eq. (20) at (n + 1/2)Δt can be approximately expressed as
Now, we introduce the shift arithmetic operators zt,
Combining Eqs. (21) and (22), we can obtain
Comparing Eq. (23) with Eq. (20), we can have
By substituting (24) into Eqs. (18) and (19), and multiplying (zt + 1)N on either side of the equal sign, we can obtain the constitutive relationship of dispersion in the time domain as
Based on the relationship
the solution of can be obtained as

The expression of can be obtained in the same way as

where
For the TE case, the incident magnetic field Hz, i, reflected magnetic field Hz, r, and transmitted magnetic field Hz, t are recorded during the time iteration. By using the Fourier transform, the transmission coefficient T and reflection coefficient R can be obtained by
where f is the frequency of the THz wave. Then the absorption coefficient is calculated by

For the TM case, the transmission coefficient and absorption coefficient can be calculated in the same way by substituting electric field Ez for magnetic field Hz.

3. Numerical results and discussion

Numerical simulations are implemented based on the theoretical treatments presented in Section 2. The incident THz wave of cosine-modulated Gaussian pulse in the SO-FDTD is given by

where fmid = 5.5 × 1011 Hz is the central frequency of the pulse, τ = 0.966/Δf stands for the width of the pulse, and t0 = 4.5τ is the time of the peak. Here Δf is the width of the frequency needed in the simulations and Δf = 9 × 1011 Hz. In what follows, this cosine-modulated Gaussian pulse incident wave is used in all the simulations.

3.1. Validity of the proposed method

Before we start to analyze the propagation characteristics of the THz wave through an inhomogeneous plasma, the validity of the method adopted in this paper is verified by comparing the results obtained by our method with those obtained by an analytical method.[35] The propagation characteristics of THz waves with different incident angles in a homogeneous plasma slab are taken into consideration.

The FDTD problem space consists of 474 cells, and the homogeneous plasma occupies 71–404 cells. The cell size is Δy = 1.5 × 10−5 m, and the thickness of the plasma layer is 5~mm. Both ends of the cell space are treated by uniaxial perfectly matched layer absorbing boundary condition to eliminate unwanted reflections. The plasma parameters are

Comparisons of the transmission coefficients of the THz wave in the plasma slab calculated by our method and those obtained using the analytical solution are displayed in Fig. 2. As shown in Fig. 2, very good agreements are achieved for both TE and TM waves, which partially indicates the correctness of our method.

Fig. 2. Comparisons of the transmission coefficients of a plasma slab versus the incident angles between analytical solutions and SO-FDTD method: (a) TE wave, (b) TM wave.

The plasma density stemmed from the reentry hypersonic flight in the stagnation region could reach 1021 m−3 and the plasma frequency corresponding to this electron density can reach 0.28 THz.[24] In the following simulations, the Gaussian electron density profile is chosen to stand for the variation of electron density, which corresponds to the non-uniform plasma flow field around the hypersonic aircraft at a high altitude of space[35]

where d = y2y1 = 8 cm is the thickness of the plasma, y1 = −d/5, y2 = 4d/5, , , in which α1 and α2 represent the distribution form of electron density. The maximum electron density is assumed to be nemax = 1021 m−3. The distribution of Gaussian profile is plotted in Fig. 3, in which the collision frequency is υc = 1.5 × 1010 Hz, and the plasma frequency ωp is calculated as

Fig. 3. Plasma electron density distribution in Gaussian profile.
3.2. Influences of incident angle

In order to investigate the influence of the incident angle on the propagation characteristics of a THz wave in plasma, the incident angle of the THz wave is varied in steps of 15° from 0° to 90° in the simulations. The transmission and absorption coefficients for TE wave are calculated and shown in Fig. 4 for θ = 0°, 15°, 30°, and 45°. As shown in Fig. 4(a), the transmission coefficient increases as the incident angle decreases or the frequency of the THz wave increases. When the incident wave frequency approaches 1 THz, the transmission coefficient tends to 0 dB, which indicates that the THz wave can be used to release the blackout problem. For the normal incidence case, when the frequency of the incident wave is f = 0.3 THz, the transmission coefficient T ≈ −20 dB. As the incident wave frequency increases, the transmission coefficient increases and tends to 0 dB. However, for θ > 15° cases, when f = 0.3 THz, the transmission coefficient is far less than −20 dB, which indicates that the THz wave cannot penetrate the plasma. As shown in Fig. 4(b), the absorption coefficient increases at first and then decreases as the incident frequency increases. For each incident angle case, there is a maximum absorption coefficient, and the corresponding frequency moves to a larger value as the incident angle increases. It indicates that when a THz wave is used in communication through plasma, the frequency corresponds to the maximum absorption coefficient should be avoided. Furthermore, it is interesting to find that the larger the incident angle is, the smoother the absorption coefficient becomes.

Fig. 4. (a) Transmission and (b) absorption coefficients versus incident angles.
3.3. Influences of thickness

Comparisons of transmission and absorption coefficients for plasma slab of different thickness are displayed in Fig. 5. As shown in Fig. 5, as the thickness of the plasma increases, the transmission coefficient decreases and the absorption coefficient increases. As displayed in Fig. 5(a), in the relatively low frequency range (approximately f < 0.35 THz), the plasma thickness has little effect on the transmission coefficient, and the incident EM wave with frequency less than 0.3 THz cannot penetrate the plasma. Similar to the cases shown in Fig. 4(b), there is an absorption peak for each thickness, as shown in Fig. 5(b). Nevertheless, the frequencies corresponding to the absorption peaks are the same for different thickness of plasma slab. It is interesting to find that the increase of the plasma slab thickness broadens the absorption peak.

Fig. 5. (a) Transmission and (b) absorption coefficients versus thickness of plasma slab.
3.4. Influences of polarization

Figure 6 illustrates the effects of polarization on the transmission and absorption coefficients. As shown in Fig. 6(a), the polarization of the incident wave has little influence on the transmission coefficient for the cases under study. The transmission coefficient increases as the incident angle decreases or the frequency increases, which is the same as we obtained in Fig. 4(a). For the cases of incident angles 30° and 45°, the transmission coefficients for TM wave and TE wave almost overlapped with each other. Figure 6(b) shows that there is a maximum absorption coefficient for each incident angle and the corresponding frequency is the same for different polarizations. The absorption coefficient increases at first and then decreases as the frequency increases. In addition, as shown in Fig. 6(b), the curve of absorption coefficient for TE wave is much smoother than that for the TM wave case.

Fig. 6. (a) Transmission and (b) absorption coefficients versus polarization.
4. Conclusion

The variations of power transmission and absorption coefficients for THz wave obliquely incident on a non-uniform plasma slab are analyzed by the SO-FDTD method combined with the phase matching condition. The factors that may affect the propagation characteristics of THz wave through plasma, including the incident angle of the incident wave, the thickness of the plasma, and the polarization of the wave, are studied. The results reveal that the incident angle can affect the propagation properties greatly, while the polarization of the incident wave has little influence on the propagation properties for the cases under study. When the THz wave is employed to the communication in the ionization blackout situation, the incident angle should be considered, and the frequency corresponding to the maximum absorption should be avoided. This study will be useful to provide a reference for the exploration of potential techniques to overcome ionization blackout problem as well as to provide insights to the interactions between THz wave and plasma.

Reference
[1] Huber P W Evans J S J C J C 1971 AIAA J. 9 1154
[2] Keidar M Kim M Boyd I 2008 J. Spacecr. Rockets 45 445
[3] Rybak J P Churchill R J 1971 IEEE Trans. Aerosp. Electron. Syst. AES-7 879
[4] Morabito D D 2002 The Spacecraft Communications Blackout Problem Encountered During Passage or Entry of Planetary Atmospheres
[5] Zang Q Bai X Ma P Huang J Ma J Yu S Shi H Sun X Liu Y Lu Y 2017 Opt. Lett. 42 687
[6] Hartunian R Stewart G Curtiss T Fergason S Seibold R Shome P 2007 AIAA Atmospheric Flight Mechanics Conference and Exhibit Hilton Head, South Carolina American Institute of Aeronautics and Astronautics 1 10
[7] Gu L Tan Z Y Cao J C 2013 Physics 42 695 in Chinese
[8] Wu M Lang T Shi G Han Z 2018 Opt. Commun. 410 520
[9] Panyaev I Dadoenkova Y Zolotovskii I Sannikov D 2018 Opt. Commun. 426 395
[10] Xie Z He J Wang X Feng S Zhang Y 2015 Opt. Lett. 40 359
[11] Rao Z Wang X Lu Y 2011 Opt. Commun. 284 5472
[12] Li J Pi Y Yang X 2014 Wirel. Commun. Mob. Com. 14 1252
[13] Li S Li J Zhu Z Cui W 2015 Int. Conf. Signal Process. 2398 2402
[14] Yuan C Zhou Z Zhang J W Xiang X 2011 J. Appl. Phys. 109 063305
[15] Yuan C Zhou Z Xiang X Sun H Wang H Xing M Luo Z 2011 Nucl. Instrum. Methods Phys. Res. B 269 23
[16] Yuan C Zhou Z Xiang X Sun H Pu S 2010 Phys. Plasmas 17 113304
[17] Ling Z Qing Z Zhang L S Jun X X 2012 Acta Phys. Sin. 61 245202 in Chinese
[18] Tosun Z Akbar D Altan H 2009 Int. Conf. Infrared Millimeter Terahertz Waves 1 2
[19] Guo L J Guo L X Li J T 2017 Phys. Plasmas 24 022108
[20] Chen J Yuan K Shen L Deng X Hong L Yao M 2016 Prog. Electromagn. Res. 157 21
[21] Wang M Li H Dong Y Li G Jiang B Zhao Q Xu J 2016 IEEE Trans. Antennas Propag. 64 286
[22] Wang M Yu M Xu Z Li G Jiang B Xu J 2015 IEEE Trans. Plasma Sci. 43 4182
[23] Zheng L Zhao Q Liu S Xing X Chen Y 2014 J. Infrared Millim. Te. 35 187
[24] Tian Y Ai X Han Y Liu X 2014 Phys. Plasmas 21 123303
[25] Gao R Yuan C Wang Y Zhou Z Gong D Fang Y Rong X 2013 J. Appl. Phys. 114 123302
[26] Yuan C X Zhou Z X Zhang J W Xiang X L Yue F Sun H G 2011 IEEE Trans. Plasma Sci. 39 1577
[27] Chen W Li J T Dan L 2017 Acta Phys. Sin. 66 84102 in Chinese
[28] Chen W Deng X J Feng J Huang G Y 2014 Acta Phys. Sin. 63 194101 in Chinese
[29] Liu S Zhou T He X Zhou Y Hong W 2006 Int. Symp. Antennas Propag. EM Theory 1 4
[30] Liu S Zhou T Liu M Hong W 2008 J. Syst. Eng. Electron. 19 15
[31] Rahmani Z Moradi H 2018 Optik 155 81
[32] Yuan K Chen J Shen L Deng X Yao M Hong L 2018 IEEE Trans. Plasma Sci. 46 373
[33] Cao Y Li H Wang Z Wu Z 2016 Int. J. Antenn. Propag. 2016 1
[34] Tian Y Yan W Gu X Jin X Li J Li B 2017 AIP Adv. 7 125325
[35] Ling Z 2013 Study ofElectromagnetic Wave Propagation in Spacecraft Plasma Sheath Ph.D. Thesis Chengdu School of Physical Electronics in Chinese
[36] Sullivan D 2013 Electromagnetic Simulation Using the FDTD Method New York IEEE Press