A fringe jump counting method for the phase measurement in the HCN laser interferometer on EAST and its FPGA-based implementation

Yuan YAO; Yao YANG; Ang TI; Yang SONG; Jiamin ZHANG; Yan WANG; Yao ZHANG; Haiqing LIU; Yinxian JIE

doi:10.1088/2058-6272/ad0dec

Plasma Science and Technology > 2024 > 26(4): 045601. > DOI: 10.1088/2058-6272/ad0dec

Yuan YAO, Yao YANG, Ang TI, Yang SONG, Jiamin ZHANG, Yan WANG, Yao ZHANG, Haiqing LIU, Yinxian JIE. A fringe jump counting method for the phase measurement in the HCN laser interferometer on EAST and its FPGA-based implementation[J]. Plasma Science and Technology, 2024, 26(4): 045601. DOI: 10.1088/2058-6272/ad0dec

Citation:

PDF (1328 KB)

A fringe jump counting method for the phase measurement in the HCN laser interferometer on EAST and its FPGA-based implementation

1.
Institute of Plasma Physics, Chinese Academy of Sciences, Hefei 230031, People’s Republic of China
2.
Anhui Electrical Engineering Professional Technique College, Hefei 230022, People’s Republic of China
3.
Institutes of Physical Science and Information Technology, Anhui University, Hefei 230601, People’s Republic of China

More Information

Author Bio:
Yao YANG: yangyao@ipp.ac.cn
Corresponding author:
Yao YANG, yangyao@ipp.ac.cn
Received Date: August 04, 2023
Revised Date: November 16, 2023
Accepted Date: November 16, 2023
Available Online: April 08, 2024
Published Date: April 10, 2024

Graphical Abstract

Abstract

Abstract

Electron density in fusion plasma is usually diagnosed using laser-aided interferometers. The phase difference signal obtained after phase demodulation is wrapped, which is also called a fringe jump. A method has been developed to unwrap the phase difference signal in real time using FPGA, specifically designed to handle fringe jumps in the hydrogen cyanide (HCN) laser interferometer on the EAST superconducting tokamak. This method is designed for a phase demodulator using the fast Fourier transform (FFT) method at the front end. The method is better adapted for hardware implementation compared to complex mathematical analysis algorithms, such as field programmable gate array (FPGA). It has been applied to process the phase measurement results of the HCN laser interferometer on EAST in real time. Electron density results show good confidence in the fringe jump unwrapping method. Further possible application in other laser interferometers, such as the POlarimeter-INTerferometer (POINT) system on EAST tokamak is also discussed.
- electron density,
- laser-aided interferometer,
- fringe jump unwrapping,
- FPGA,
- EAST

FullText(HTML)

1. Introduction

THz vortex beams with OAM (orbital angular momentum), have potential applications in many fields, including THz wireless communication [1–4], super-resolution imaging [5], manipulating matters [6], spectroscopic probe of magnetism in matter [7], and so have attracted the attention of many researchers. Similar to optical OAM, by using the generation principles of wave-front modulators and THz hologram technology, a large number of methods have been proposed, such as diffractive elements [8], metasurfaces [9, 10], an ultrathin planar THz vortex phase plate [11], a spiral phase plate [5, 12], THz computer-generated holograms [13] and so on. Moreover, based on the nonlinear effects of plasma [14], THz vortex can be generated by beating two vortex lasers in rippled plasma [15, 16], and by transferring the OAM from the plasma vortex to the emitted THz wave [17–19]. Also using plasma, a THz pulse with OAM is obtained by employing two twisted input lasers [20]. What is more, some multiplexing/demultiplexing systems of THz OAM also have been demonstrated successfully [21–23], and the results indicate the potential of THz OAM in improving channel capacity.

Specifically, it is worth noting that the impact of propagation environment on THz OAM communication performance is one of the key points. According to this, Zhao et al discussed the modal coupling and crosstalk of multiple THz OAM beams propagating in turbulence, and the results show that the divergence and turbulence effects in the THz regime are different from those in the mm-wave and optical regime [1]. As we all know, the plasma environment widely exists in nature, laboratory and industry, such as ionosphere, interstellar medium, and laser fusion. In particular, the effects of plasma sheath, which is generated during the re-entering process of spacecraft, on tracking, telemetry and control (TT&C) communication have become a key problem in the application of hypersonic vehicles. In order to overcome the 'blackout' phenomenon caused by plasma sheath, many ways have been proposed, for example, applying a magnetic field [24, 25]. Meanwhile, owing to its good penetrability, THz wave is proved to be a new method to solve 'blackout'. Many studies on the propagation characteristics of THz waves in plasma sheath [26–29] have been reported. Further, considering the existing ablative particles, the sheath can be seen as dusty plasma environments, and the effects of dusty plasma with the ceramic substrate on the propagation of THz waves have been analyzed by the scattering matrix method (SMM) [30]. After discussing the effects of the inhomogeneity of plasma sheath, the incident direction, the flight speed, as well as the applied magnetic field, the conclusion shows that THz waves may raise the propagation performance of electromagnetic waves in plasma sheath. Therefore, in light of the application value of the new dimension-OAM in increasing channel capacity, it is prospective to explore the transmission of THz vortex beams in a plasma environment. In the earlier stage, with the help of the plane wave angular spectrum expansion method, we have carried out the transmission of a THz vortex beam and OAM crosstalk of multiple coaxial THz vortex beams in an unmagnetized plasma slab [31, 32]; however, for magnetized plasma, this method is not applicable. Thus, combining the vector wave function expansion with Fourier transform representation of an electromagnetic wave in magnetized plasma, we have proposed a new method to solve the transmission of a THz vortex beam passing through uniform magnetized plasma [33]. Whereas, for most environments, the distribution of the plasma density is nonuniform, which means that the transmission in non-uniform magnetized plasma is worthy of further study. This is what the work focuses on.

In this work, the typical non-diffraction Bessel vortex beam is chosen as the incident source, and the theoretical model is presented in section 2, including the expansion of electromagnetic fields in plasma regions in terms of cylindrical vector wave functions (CVWFs) and Fourier transform, and derivation of a cascade form of expansion coefficients. Section 3 discusses the effects of the applied magnetized field, the incident angle and the polarization mode upon the characteristics of the transmitted fields. Finally, section 4 summarizes the major conclusions.

2. Theoretical model

Figure 1 depicts a THz Bessel vortex beam impinging upon a multi-layered magnetized plasma slab, where the global coordinate system $O x y z$ is based on the plasma slab. In order to describe the incident, reflected and transmitted beams clearly, three local coordinate systems $O_{ß} x_{ß} y_{ß} z_{ß}$ ( $ß = i, r, or t$ ) are established, in which the positive $z_{ß}$ -axes are along the central axes of the corresponding beams. The incident angle $θ$ is the angle between the $z_{i}$ and $z$ axes. The non-uniform magnetized plasma slab is split into N layers, and the thickness of the jth layer of the magnetized plasma slab is $d_{j} = z_{j} - z_{j - 1} .$ Additionally, the time factor $\exp (- i ω t)$ of the electromagnetic field is omitted.

Figure 1. Schematic diagram of a Bessel vortex beam incident upon a multi-layered magnetized plasma slab.

DownLoad: Full-Size Img PowerPoint

Assuming the direction of the applied magnetic field $B_{0}$ is along the positive z-axis, then, the permittivity tensor of plasma in the jth layer is represented as

${\mathop \varepsilon \limits^ = _j} = \left[ {\begin{array}{*{20}{l}} {{\varepsilon _{j,1}}}&{ - {\varepsilon _{j,2}}}&0\\ {{\varepsilon _{j,2}}}&{{\varepsilon _{j,1}}}&0\\ 0&0&{{\varepsilon _{j,3}}} \end{array}} \right]$

(1)

where $ε_{j, 1} = ε_{0} (1 + \frac{ω_{j, p}^{2} / ω (ω + i v_{j, e})}{[ω_{c}^{2} - {(ω + i v_{j, e})}^{2}]}),$ $ε_{j, 2} = ε_{0} (\frac{i ω_{c} \cdot ω_{j, p}^{2} / ω}{[ω_{c}^{2} - {(ω + i v_{j, e})}^{2}]}),$ $ε_{j, 3} = ε_{0} (1 - \frac{ω_{j, p}^{2} / ω}{(ω + i v_{j, e})}),$ $ω_{c},$ $ν_{j, e},$ $ω_{j, p}$ and $ω$ are the cyclotron frequency, the collision frequency, the angular plasma frequency and the angular frequency of the incident wave, respectively. $ε_{0}$ is the permittivity of vacuum, and $i$ is the imaginary unit.

2.1 Description of the incident THz Bessel vortex beam

The expression of a vector Bessel vortex beam is presented in terms of angular spectrum decomposition method as follows [34]:

$E (r) = \int_{0}^{2 π} E_{0} Q (α, φ) |_{α = α_{0}} e^{i l φ} e^{i k_{0} \cdot r} d φ,$

(2)

where $E_{0}$ and $k_{0}$ are the amplitude and the wave vector of the electric field component of the plane waves propagating over the cone with the half-angle $α_{0}$ and the azimuthal angle $φ,$ respectively. $Q (α, φ)$ is the vector complex polarization function, and $l$ is the topological charge (TC) of the vortex beam.

As a THz Bessel vortex beam is incident obliquely on the plasma slab, making use of CVWFs ( $m_{m λ_{0}}^{(1)},$ $n_{m λ_{0}}^{(1)}$ ) [35] in the global coordinate system, the electromagnetic field vectors of the incident beam can be expanded and written as:

$\begin{matrix} E^{inc} = E_{0} \sum_{m = - \infty}^{\infty} \int_{0}^{π} [A_{m}^{inc} (ζ) m_{m λ_{0}}^{(1)} (h_{0}) \\ + B_{m}^{inc} (ζ) n_{m λ_{0}}^{(1)} (h_{0})] e^{i h_{0} z} d ζ, \end{matrix}$

(3)

$\begin{matrix} H^{inc} = \frac{k_{0}}{i ω μ_{0}} E_{0} \sum_{m = - \infty}^{\infty} \int_{0}^{π} [A_{m}^{inc} (ζ) n_{m λ_{0}}^{(1)} (h_{0}) \\ + B_{m}^{inc} (ζ) m_{m λ_{0}}^{(1)} (h_{0})] e^{i h_{0} z} d ζ, \end{matrix}$

(4)

where $k_{0} = λ_{0} \hat{ρ} + h_{0} \hat{z},$ $|k_{0}| = k_{0},$ $λ_{0} = k_{0} \sin ζ,$ $h_{0} = k_{0} \cos ζ,$ $ζ$ represents the angle between $k_{0}$ and the positive $z$ -axis, $μ_{0}$ is the permeability, and the superscripts $inc$ and (1) stand for the incident beam and the first kind of Bessel function respectively. The coefficients $A_{m}^{inc} (ζ),$ $B_{m}^{inc} (ζ)$ are derived and shown as:

$A_{m}^{inc} (ζ) = \sum_{n = |m|}^{\infty} [{a'}_{m n} c_{m n} (ζ) + {b'}_{m n} d_{m n} (ζ)],$

(5)

$B_{m}^{inc} (ζ) = \sum_{n = |m|}^{\infty} [{a'}_{m n} d_{m n} (ζ) + {b'}_{m n} c_{m n} (ζ)],$

(6)

here $c_{m n} = i^{m - n - 1} τ_{m n} (\cos ζ) / 2 k_{0}, d_{m n} = i^{m - n - 1} π_{m n} (\cos ζ) / 2 k_{0},$ and $({a'}_{m n}, {b'}_{m n})$ are the expansion coefficients of the Bessel vortex beam in terms of spherical vector wave functions in the global coordinate system [36]. $π_{m n} (\cos ζ) = m P_{n}^{m} (\cos ζ) / \sin ζ$ and $τ_{m n} (\cos ζ) = d P_{n}^{m} (\cos ζ) / d ζ$ are the angle functions, where $P_{n}^{m} (\cos ζ)$ is the associated Legendre function. In the following analysis, only the integral interval $[0, π / 2]$ is meaningful.

2.2 CVWFs expansion of electromagnetic fields in anisotropic magnetized plasma

In the jth layer of a multi-layered magnetized plasma slab, the wave equation of the electric field $E_{j} (r)$ in the source-free region is

$\nabla \times \nabla \times E_{j} (r) - ω^{2} μ_{0} {\bar{\bar{ε}}}_{j} \cdot E_{j} (r) = 0 .$

(7)

By using Fourier transform, the electric field in equation (7) can be shown as

$E_{j} (r) = \int_{- \infty}^{+ \infty} \int_{- \infty}^{+ \infty} \int_{- \infty}^{+ \infty} E_{j} (k) e^{i k \cdot r} d k_{x} d k_{y} d k_{z} .$

(8)

After a series of analysis on eigenvalue and eigenvector [37], the electric field in plasma can be expressed as

$E_{j} (r) = \sum_{q = 1}^{2} \int_{0}^{2 π} \int_{0}^{π} E_{jq} (θ_{k}, ϕ_{k}) e^{i k_{jq} \cdot r} k_{jq}^{2} \sin θ_{k} d θ_{k} d ϕ_{k}$

(9)

where $E_{j q} (θ_{k}, ϕ_{k}) = F_{j q} (θ_{k}, ϕ_{k}) f_{j q} (θ_{k}, ϕ_{k})$ is the eigenvector, which is corresponding to the eigenvalue $k_{jq} .$ The subscript $k$ indicates that the integral is solved in the k-domain coordinate system. If we expand $f_{j q} (θ_{k}, ϕ_{k})$ to $f_{j q} (θ_{k}, ϕ_{k}) = \sum_{m = - \infty}^{\infty} {G'}_{j m q} (θ_{k}) e^{i m ϕ_{k}},$ equation (9) can be rewritten as

$\begin{matrix} E_{j} (r) = E_{0} \sum_{q = 1}^{2} \sum_{m = - \infty}^{\infty} \int {G'}_{j m q} (θ_{k}) k_{j q}^{2} \sin θ_{k} d θ_{k} \\ \times \int_{0}^{2 π} F_{j q} (θ_{k}, ϕ_{k}) e^{i k_{j q} \cdot r} e^{i m ϕ_{k}} d ϕ_{k} \end{matrix}$

(10)

Here, the eigenvalues are $k_{j 1}^{2} = (B B_{j} + \sqrt{B B_{j}^{2} - 4 A A_{j} \cdot C C_{j}}) / 2 A A_{j},$ $k_{j 2}^{2} = (B B_{j} - \sqrt{B B_{j}^{2} - 4 A A_{j} \cdot A C_{j}}) / 2 A A_{j},$ where $A A_{j} = a_{j 3},$ $B B_{j} = 2 a_{j 1} a_{j 3} + (a_{j 3} - a_{j 1}) λ_{j q}^{2},$ $C C_{j} = a_{j 3} (a_{j 1}^{2} + a_{j 2}^{2}) - (a_{j 1}^{2} + a_{j 2}^{2} - a_{j 1} a_{j 3}) λ_{j q}^{2},$ $λ_{j q} = k_{j q} \sin θ_{k}$ and $a_{j p} = ω^{2} μ_{0} ε_{j p} (p = 1, 2, 3) .$

By expanding $F_{jq} (θ_{k}, ϕ_{k})$ in terms of CVWFs ( $m_{m λ_{j q}}^{(1)},$ $n_{m λ_{j q}}^{(1)},$ $l_{m λ_{j q}}^{(1)}$ ), and through some reasonable simplification, equation (10) is expressed as

$\begin{matrix} E_{j} = E_{0} \sum_{q = 1}^{2} \sum_{m = - \infty}^{\infty} \int_{0}^{π} G_{j m q} (θ_{k}) [α_{j q} (θ_{k}) m_{m λ_{j q}}^{(1)} \\ + ß_{j q} (θ_{k}) n_{m λ_{j q}}^{(1)} + γ_{j q} (θ_{k}) l_{m λ_{j q}}^{(1)}] e^{i h_{j q} z_{j}} d θ_{k}, \end{matrix}$

(11)

where, $G_{j m q} (θ_{k}) = 2 π i^{m - 1} {G'}_{j m q} (θ_{k}) k_{j q} \sin θ_{k},$ and $h_{j q} = k_{j q} \cos θ_{k} .$ The expansion coefficients are

$\begin{matrix} α_{j q} = - a_{j 2} λ_{j q} h_{j q} / ∆_{j q} \\ ß_{j q} = \frac{- i (k_{j q}^{2} - a_{j 1}) λ_{j q} h_{j q}^{2}}{k_{j q} ∆_{j q}} + \frac{i λ_{j q}}{k_{j q}} \\ γ_{j q} = \frac{(k_{j q}^{2} - a_{j 1}) λ_{j q}^{3} h_{j q}}{k_{j q}^{2} ∆_{j q}} + \frac{λ_{j q} h_{j q}}{k_{j q}^{2}} \\ ∆_{j q} = k_{j q}^{4} - (λ_{j q}^{2} + 2 a_{j 1}) k_{j q}^{2} + a_{j 1} λ_{j q}^{2} + (a_{j 1}^{2} + a_{j 2}^{2}) . \end{matrix}$

(12)

In order to describe the transmission of the electromagnetic waves in each layer of the magnetized plasma slab more vividly, we decompose equation (11) into two parts:

$\begin{matrix} E_{j 1} = E_{0} \sum_{q = 1}^{2} \sum_{m = - \infty}^{\infty} \int_{0}^{\frac{π}{2}} E_{j m q} (ζ) \\ \times [\begin{matrix} α_{j q} (ζ) m_{m λ_{j q}}^{(1)} (h_{j q}) + ß_{j q} (ζ) n_{m λ_{j q}}^{(1)} (h_{j q}) \\ + γ_{j q} (ζ) l_{m λ_{j q}}^{(1)} (h_{j q}) \end{matrix}] e^{i h_{j q} z_{j}} d ζ, \end{matrix}$

(13)

$\begin{matrix} E_{j 2} = E_{0} \sum_{q = 1}^{2} \sum_{m = - \infty}^{\infty} \int_{0}^{\frac{π}{2}} F_{j m q} (ζ) \\ \times [\begin{matrix} - α_{j q} (ζ) m_{m λ_{j q}}^{(1)} (- h_{j q}) + ß_{j q} (ζ) n_{m λ_{j q}}^{(1)} (- h_{j q}) \\ - γ_{j q} (ζ) l_{m λ_{j q}}^{(1)} (- h_{j q}) \end{matrix}] \\ \times e^{- i h_{j q} z_{j}} d ζ, \end{matrix}$

(14)

where equations (13) and (14) denote the downgoing and the upgoing waves in each plasma layer, respectively. Note that the variable $θ_{k}$ is replaced by $ζ .$ $E_{j m q}$ and $F_{j m q}$ are the unknown expanding coefficients.

2.3 Reflected and transmitted beams by a multi-layered magnetized plasma slab

Similar to equations (3) and (4), the electromagnetic fields of the reflected beam by a multi-layered magnetized plasma slab can be expanded as

$\begin{matrix} E^{ref} = E_{0} \sum_{m = - \infty}^{\infty} \int_{0}^{\frac{π}{2}} [C_{m}^{ref} (ζ) m_{m λ_{0}}^{(1)} (- h_{0}) \\ + D_{m}^{ref} (ζ) n_{m λ_{0}}^{(1)} (- h_{0})] e^{- i h_{0} z} d ζ, \end{matrix}$

(15)

$\begin{matrix} H^{ref} = E_{0} \frac{k_{0}}{i ω μ_{0}} \sum_{m = - \infty}^{\infty} \int_{0}^{\frac{π}{2}} [C_{m}^{ref} (ζ) n_{m λ_{0}}^{(1)} (- h_{0}) \\ + D_{m}^{ref} (ζ) m_{m λ_{0}}^{(1)} (- h_{0})] e^{- i h_{0} z} d ζ, \end{matrix}$

(16)

and the electromagnetic fields of the transmitted beam are

$\begin{matrix} E^{tra} = E_{0} \sum_{m = - \infty}^{\infty} \int_{0}^{\frac{π}{2}} [G_{m}^{tra} (ζ) m_{m λ_{0}}^{(1)} (h_{0}) \\ + I_{m}^{tra} (ζ) n_{m λ_{0}}^{(1)} (h_{0})] e^{i h_{0} z} d ζ, \end{matrix}$

(17)

$\begin{matrix} H^{tra} = E_{0} \frac{k_{0}}{i ω μ_{0}} \sum_{m = - \infty}^{\infty} \int_{0}^{\frac{π}{2}} [G_{m}^{tra} (ζ) n_{m λ_{0}}^{(1)} (h_{0}) \\ + I_{m}^{tra} (ζ) m_{m λ_{0}}^{(1)} (h_{0})] e^{i h_{0} z} d ζ, \end{matrix}$

(18)

where $C_{m}^{ref},$ $D_{m}^{ref},$ $G_{m}^{tra}$ and $I_{m}^{tra}$ are the unknown expansion coefficients of the reflected and transmitted beams, respectively.

Based on the boundary conditions of the electromagnetic fields on every interface, we have

$\begin{matrix} \begin{matrix} E_{r}^{inc} + E_{r}^{ref} = E_{11 r} + E_{12 r}, E_{ϕ}^{inc} + E_{ϕ}^{ref} = E_{11 ϕ} + E_{12 ϕ} \\ H_{r}^{inc} + H_{r}^{ref} = H_{11 r} + H_{12 r}, H_{ϕ}^{inc} + H_{ϕ}^{ref} = H_{11 ϕ} + H_{12 ϕ} \end{matrix}\} \\ at z = 0 \end{matrix}$

(19)

$\begin{matrix} \begin{matrix} E_{j 1 ϕ} + E_{j 2 ϕ} = E_{(j + 1) 1 ϕ} + E_{(j + 1) 2 ϕ} \\ H_{j 1 ϕ} + H_{j 2 ϕ} = H_{(j + 1) 1 ϕ} + H_{(j + 1) 2 ϕ} \end{matrix}\} at \\ z = z_{j} (j = 1, 2, \cdot\cdot\cdot, N - 1) \end{matrix}$

(20)

$\begin{matrix} \begin{matrix} E_{N 1 r} + E_{N 2 r} = E_{r}^{tra}, E_{N 1 ϕ} + E_{N 2 ϕ} = E_{ϕ}^{tra} \\ H_{N 1 r} + H_{N 2 r} = H_{r}^{tra}, H_{N 1 ϕ} + H_{N 2 ϕ} = H_{ϕ}^{tra} \end{matrix}\} at \\ z = z_{N} \end{matrix}$

(21)

Substituting equations (3), (4) and (13)–(18) into equations (19)–(21), the relationships of expansion coefficients are obtained as follows:

$[\begin{matrix} κ_{11}^{0} & κ_{12}^{0} & κ_{13}^{0} & κ_{14}^{0} \\ κ_{21}^{0} & κ_{22}^{0} & κ_{23}^{0} & κ_{24}^{0} \\ κ_{31}^{0} & κ_{32}^{0} & κ_{33}^{0} & κ_{34}^{0} \\ κ_{41}^{0} & κ_{42}^{0} & κ_{43}^{0} & κ_{44}^{0} \end{matrix}] [\begin{matrix} E_{1 m 1} \\ E_{1 m 2} \\ F_{1 m 1} \\ F_{1 m 2} \end{matrix}] = [\begin{matrix} A_{m}^{inc} \\ B_{m}^{inc} \\ C_{m}^{ref} \\ D_{m}^{ref} \end{matrix}] at z = 0$

(22)

$\begin{matrix} [\begin{matrix} E_{j m 1} \\ E_{j m 2} \\ F_{j m 1} \\ F_{j m 2} \end{matrix}] = [\begin{matrix} κ_{11}^{j} & κ_{12}^{j} & κ_{13}^{j} & κ_{14}^{j} \\ κ_{21}^{j} & κ_{22}^{j} & κ_{23}^{j} & κ_{24}^{j} \\ κ_{31}^{j} & κ_{32}^{j} & κ_{33}^{j} & κ_{34}^{j} \\ κ_{41}^{j} & κ_{42}^{j} & κ_{43}^{j} & κ_{44}^{j} \end{matrix}] [\begin{matrix} E_{(j + 1) m 1} \\ E_{(j + 1) m 2} \\ F_{(j + 1) m 1} \\ F_{(j + 1) m 2} \end{matrix}] at \\ z = z_{j} (j = 1, 2, \cdot\cdot\cdot, N - 1) \end{matrix}$

(23)

$\begin{matrix} [\begin{matrix} E_{N m 1} \\ E_{N m 2} \\ F_{N m 1} \\ F_{N m 2} \end{matrix}] = [\begin{matrix} κ_{11}^{N} & κ_{12}^{N} & κ_{13}^{N} & κ_{14}^{N} \\ κ_{21}^{N} & κ_{22}^{N} & κ_{23}^{N} & κ_{24}^{N} \\ κ_{31}^{N} & κ_{32}^{N} & κ_{33}^{N} & κ_{34}^{N} \\ κ_{41}^{N} & κ_{42}^{N} & κ_{43}^{N} & κ_{44}^{N} \end{matrix}] [\begin{matrix} G_{m}^{tra} e^{i h_{0} z_{N}} \\ I_{m}^{tra} \frac{h_{0}}{k_{0}} e^{i h_{0} z_{N}} \\ G_{m}^{tra} h_{0} e^{i h_{0} z_{N}} \\ I_{m}^{tra} e^{i h_{0} z_{N}} \end{matrix}] at \\ z = z_{N} \end{matrix}$

(24)

where the expression of the coefficient matrices is shown in the appendix. After a series of simplifications, the cascade form of the unknown expansion coefficients can be obtained and devoted by

$\prod_{j = 0}^{N} [\begin{matrix} κ_{11}^{j} & κ_{12}^{j} & κ_{13}^{j} & κ_{14}^{j} \\ κ_{21}^{j} & κ_{22}^{j} & κ_{23}^{j} & κ_{24}^{j} \\ κ_{31}^{j} & κ_{32}^{j} & κ_{33}^{j} & κ_{34}^{j} \\ κ_{41}^{j} & κ_{42}^{j} & κ_{43}^{j} & κ_{44}^{j} \end{matrix}] [\begin{matrix} G_{m}^{tra} e^{i h_{0} z_{N}} \\ I_{m}^{tra} \frac{h_{0}}{k_{0}} e^{i h_{0} z_{N}} \\ G_{m}^{tra} h_{0} e^{i h_{0} z_{N}} \\ I_{m}^{tra} e^{i h_{0} z_{N}} \end{matrix}] = [\begin{matrix} A_{m}^{inc} \\ B_{m}^{inc} \\ C_{m}^{ref} \\ D_{m}^{ref} \end{matrix}] .$

(25)

Let us make

$\prod_{j = 0}^{n} [\begin{matrix} κ_{11}^{j} & κ_{12}^{j} & κ_{13}^{j} & κ_{14}^{j} \\ κ_{21}^{j} & κ_{22}^{j} & κ_{23}^{j} & κ_{24}^{j} \\ κ_{31}^{j} & κ_{32}^{j} & κ_{33}^{j} & κ_{34}^{j} \\ κ_{41}^{j} & κ_{42}^{j} & κ_{43}^{j} & κ_{44}^{j} \end{matrix}] = [\begin{matrix} κ_{11} & κ_{12} & κ_{13} & κ_{14} \\ κ_{21} & κ_{22} & κ_{23} & κ_{24} \\ κ_{31} & κ_{32} & κ_{33} & κ_{34} \\ κ_{41} & κ_{42} & κ_{43} & κ_{44} \end{matrix}] .$

(26)

Then the expansion coefficients of the reflected and transmitted beams are

$\begin{matrix} G_{m}^{tra} = \frac{A_{m}^{inc} (κ_{22} \frac{h_{0}}{k_{0}} + κ_{24}) - B_{m}^{inc} (κ_{12} \frac{h_{0}}{k_{0}} + κ_{14})}{(κ_{11} + κ_{13} h_{0}) (κ_{22} \frac{h_{0}}{k_{0}} + κ_{24}) e^{i h_{0} z_{N}} - (κ_{21} + κ_{23} h_{0}) (κ_{12} \frac{h_{0}}{k_{0}} + κ_{14}) e^{i h_{0} z_{N}}} \\ I_{m}^{tra} = \frac{A_{m}^{inc} (κ_{21} + κ_{23} h_{0}) - B_{m}^{inc} (κ_{11} + κ_{13} h_{0})}{(κ_{12} \frac{h_{0}}{k_{0}} + κ_{14}) (κ_{21} + κ_{23} h_{0}) e^{i h_{0} z_{N}} - (κ_{22} \frac{h_{0}}{k_{0}} + κ_{24}) (κ_{11} + κ_{13} h_{0}) e^{i h z_{N}}} \\ C_{m}^{ref} = (κ_{31} + κ_{33} h_{0}) e^{i h_{0} z_{N}} G_{m}^{tra} + (κ_{32} \frac{h_{0}}{k_{0}} + κ_{34}) e^{i h_{0} z_{N}} I_{m}^{tra} \\ D_{m}^{ref} = (κ_{41} + κ_{43} h_{0}) e^{i h_{0} z_{N}} G_{m}^{tra} + (κ_{42} \frac{h_{0}}{k_{0}} + κ_{44}) e^{i h_{0} z_{N}} I_{m}^{tra}, \end{matrix}$

(27)

when the inhomogeneous magnetized plasma slab degenerates into a homogeneous one, the results of equation (27) are in conformity with that in [36].

2.4 OAM spectrum of the distorted vortex beam

Owing to the effects of the magnetized plasma slab, the vorticity of the incident Bessel vortex beam will be distorted. For example, except the dominant OAM state (whose TC is the same as that of the incident beam), some other OAM states are often derived in the transmitted beam. Normally, the weight distribution of the OAM states of the distorted beam is used to describe the vorticity of the transmitted beam, and the weight of each OAM state is defined as [33]:

$P_{n} = \frac{W_{n}}{\sum_{n = - \infty}^{\infty} W_{n}}$

(28)

where $W_{n} = \int w_{n} ρ d ρ$ is the electromagnetic field energy of the OAM state $n .$

3. Results and discussion

Considering that the plasma density of the slab is double Gaussian distribution, that is

$N_{e} (z) = \{\begin{matrix} N_{ep} \exp [- a_{1} {(z - z_{p})}^{2}], (0 \leq z \leq z_{p}) \\ N_{ep} \exp [- a_{2} {(z - z_{p})}^{2}], (z > z_{p}) \end{matrix},$

(29)

where, $a_{1}$ and $a_{2}$ represent the density gradient in ascending and descending processes, respectively. $N_{ep}$ is the maximum value of the plasma density, while $z_{p}$ is the corresponding coordinate value. Figure 2 presents the distribution of the plasma density, where $N_{ep} = 1.7 \times 10^{19} m^{- 3},$ $a_{1} = 1000 m^{- 2},$ $a_{2} = 500 m^{- 2},$ $z_{p} = 0.03 m$ and the thickness of the plasma slab is 10 cm.

Figure 2. Double Gaussian distribution of plasma density.

DownLoad: Full-Size Img PowerPoint

Taking the inhomogeneous plasma slab given in figure 2 as an example, the correctness of the method provided in this paper is numerically verified. The plasma slab is split into 100 sublayers. The collision frequency and the magnetic induction intensity are set to 20 GHz and $B_{0} = |B_{0}| = 1 T,$ respectively. As we know, by setting the TC $l$ and the half-cone angle $α_{0}$ to zero, the Bessel vortex beam may degenerate into a plane wave. As a plane wave normally incidence, the variation of power transmission coefficients of the right-handed circularly polarized (RCP) wave with frequency is calculated by both SMM and our method and shown in figure 3. The results indicate that they are in good agreement with each other.

Figure 3. Variation of the transmission coefficients of the RCP wave with the incident frequency.

DownLoad: Full-Size Img PowerPoint

In the following subsection, by using the double Gaussian distribution of the plasma density (shown in equation (29)), the effects of parameters of both the multi-layered magnetized plasma slab and the incident beam on the transmitted beams are analyzed numerically in detail. Without specification, the parameters used in equation (29) are the same as those in figure 2.

3.1 Magnitude and OAM mode of the transmitted beams for a Bessel vortex beam incidence

A comparison of the amplitude profiles of the transmitted electric field components $E_{x}^{tra}$ and $E_{y}^{tra},$ with three different magnetic induction intensities is displayed in figure 4, where Ext and Eyt are corresponding to $E_{x}^{tra}$ and $E_{y}^{tra},$ respectively. The frequency and the incident angle of the x polarized Bessel vortex beam are 0.1 THz and 30°, respectively. The half cone angle and the TC are set to $α_{0} = 5 °$ and $l = 2 .$ The thickness, collision frequency and the number of sublayers of the magnetized plasma slab are the same as those of figure 3.

Figure 4. Effects of magnetic induction intensity

$B_{0}$ upon the transmitted electric field components

$E_{x}^{tra}$ and

$E_{y}^{tra}$ as an x-polarized Bessel vortex beam incidence.

DownLoad: Full-Size Img PowerPoint

The contours shown in figure 4 indicate that as the magnetic induction intensity $B_{0}$ increases, so does the magnitude of the component $E_{y}^{tra} .$ This is because the Faraday rotation angle is proportional to $B_{0},$ resulting in the gradual transformation of the polarization mode of the beam from x-linear polarization to elliptical polarization. Nevertheless, due to the small impact of the increase of $B_{0}$ on the permittivity elements $ε_{j, 1}$ and $ε_{j, 2},$ the differences of the eigenvalues $k_{j q}$ among the three cases are small, that is, the deviation of the propagation direction between two characteristic waves in plasma is small, leading to the tiny effect of $B_{0}$ on the profiles of both $E_{x}^{tra}$ and $E_{y}^{tra} .$

The effects of the incident angle upon the profiles of the total transmitted electric field amplitude ( $E_{T}$ ) are presented in figure 5. The parameters are the same as those used in figure 4, except $B_{0} = 1 T .$ The results reveal that with an increase of $θ,$ the magnitude of $E_{T}$ decreases due to the plasma attenuation, and the center of the transmitted beam is gradually away from the origin of the local coordinate system $O_{t} x_{t} y_{t} z_{t} .$ Moreover, the deviation of the profile of the transmitted electric field from that of the incident field is becoming more and more serious.

Figure 5. Profiles of the transmitted electric field amplitude with different incident angles.

DownLoad: Full-Size Img PowerPoint

For different applied magnetic fields and incident angle parameters, figure 6 shows the distribution of the OAM spectrum of the transmitted beam. Parameters of both the beam and plasma in figures 6(a) and (b) are the same as those of figures 4 and 5, respectively. The results reveal that B ₀ does not affect the OAM state distribution, which is in line with that of a homogeneous plasma slab [33]. However, as the incident angle increases, the phase distribution on the wave front of the transmitted beam gradually deviates from that of the incident vortex beam, which reduces the weight of the dominant OAM state.

Figure 6. Effects of the magnetic induction intensity and the incident angle on OAM spectrum of the transmitted beam, (a) effects of

$B_{0}$ on

$P_{t},$ (b) effects of

$θ$ on

$P_{t} .$

DownLoad: Full-Size Img PowerPoint

3.2 Effects of the magnetic induction intensity upon the polarization mode of the transmitted beam

In this section, for the incident beams with different polarization states, the effects of $B_{0}$ upon the polarization mode of the transmitted beam are simulated numerically. In order to describe the evolution process vividly, figure 7 illustrates the variation of the trajectories of the transmitted electric field vector with $B_{0} .$ The labels on the left side of the diagram x-pol, circular-pol, and radial-pol denote the x-polarization, the circular polarization, and the radial polarization, respectively. The title ${Pol}_{i}$ represents the trajectories of the incident beam, and ${Pol}_{t}$ is corresponding to that of the transmitted beam. The magnitude of magnetic induction intensity $B_{0} = |B_{0}|$ follows ${Pol}_{t} .$ Except for the incident angle $θ = 0 °,$ other parameters of both the incident beam and the magnetic plasma slab are the same as those in figure 4.

Figure 7. Variation of the transmitted electric field vector trajectories with the polarization of the incident beam and the magnetic induction intensity

$B_{0} .$

DownLoad: Full-Size Img PowerPoint

For an x-polarized incident beam incidence, the vector trajectories are parallel to the x-axis. As the increase of $B_{0},$ the angle between the trajectory lines and the x-axis increases and gradually evolves to ellipses and further circles. In other words, with an increase of $B_{0},$ the polarization mode of the transmitted beam changes from linear polarization to elliptical polarization and further to circular polarization. The results are in accord with those shown in figure 3, and can be explained by the Faraday rotation effects of magnetized plasma.

The figure shows that as the incident beam is circular polarization mode, the trajectory lines of the transmitted electric fields are not affected by $B_{0}$ and still circular. It means that similar to plane waves, no matter what large $B_{0}$ is, propagation of circularly polarized vortex beam in magnetized plasma will not produce Faraday rotation effect.

When the incident beam turns into radial polarization, its electric field trajectories are either linear or elliptical. After passing through the magnetized plasma slab, with the increase of $B_{0},$ the electric field trajectory lines of the transmitted field will rotate and change in shape. Furthermore, as $B_{0} = 1 T,$ the trajectory lines become circular shape, that is, the transmitted beam evolves into circular polarization.

In conclusion, like the plane waves, the linear vortex beam propagating in the magnetized plasma will have the Faraday rotation effect, but the circular polarization will not. With the increase of $B_{0},$ the polarization of the transmitted beam varies for both x polarization and radial polarization.

3.3 Coaxial incidence of multiple Bessel vortex beams with different TCs

OAM multiplexing in free space communication is widely concerned in wireless optical communication technology, whereas the mode crosstalk of OAM is one of the key points that must be considered and overcome. In this section, for coaxial incidence of multiple Bessel vortex beams with different TCs, using the method mentioned in the above sections, by modifying the expression of the incident beams to

$\begin{matrix} E^{inc} = E_{0} \sum_{q_{0} = 1}^{N_{0}} \sum_{m = - \infty}^{\infty} \int_{0}^{π} [A_{q_{0} m}^{inc} (ζ) m_{m λ_{0}}^{(1)} (h_{0}) \\ + B_{q_{0} m}^{inc} (ζ) n_{m λ_{0}}^{(1)} (h_{0})] e^{i h_{0} z} d ζ \end{matrix}$

(30)

the magnitude profiles and OAM spectrum characteristics of the transmitted beam can be analyzed. $N_{0}$ is the number of the Bessel vortex beams, and the subscript $q_{0}$ denotes the $q_{0}$ th beam.

Taking coaxial incidence of several Bessel vortex beams with $l = - 2, 2$ and $l = - 3, - 1, 1, 3$ as examples, figure 8 presents the variation of magnitude profiles of the transmitted beam with the maximum value $N_{ep}$ of plasma density. The beam parameters are set to $α_{0} = 8 °,$ $θ = 30 °,$ $f = 0.1 THz,$ and the parameters of the plasma slab are the same as those of figure 2.

Figure 8. Effects of

$N_{ep}$ on magnitude of the transmitted beam for coaxial incident beams, (a)–(c)

$l = - 2, 2,$ (d)–(f)

$l = - 3, - 1, 1, 3 .$

DownLoad: Full-Size Img PowerPoint

The results shown in figure 8 indicate that the magnitude profile of the total transmitted electric field $E_{T}$ is closely related to the TC of the superimposed beams and decreases with the increase of $N_{ep} .$ Comparing with other regions, distortion in the region ( $- 15 λ \leq x \leq 0,$ $- 15 λ \leq y \leq 15 λ$ ) seems more serious. As $N_{ep}$ increases to $8 \times 10^{19} m^{- 3},$ the amplitude of the total transmitted electric field reduces to $10^{- 4} .$ This is because the increase of electron density significantly increases the angular plasma frequency, resulting in the serious attenuation of the electromagnetic wave.

Figure 9 shows the weight distribution of the OAM states of the transmitted beams. Obviously, except the dominant OAM states ( $- 2, 2$ in figure 9(a), $- 3, - 1, 1, 3$ in figure 9(b)), some adjacent states are derived, and with the growth of $N_{ep},$ the weight of the dominant OAM state diminishes (which is in accord with the results in [38]). It means that the mutual OAM coupling and crosstalk are very significant. Furthermore, for the cases of $N_{ep} = 5 \times 10^{19} m^{- 3}$ and $N_{ep} = 8 \times 10^{19} m^{- 3},$ the difference of adjacent TCs of the incident beams plays an important role in the mutual OAM coupling and crosstalk of the transmitted beam, such as the growth of the weight of the OAM state $n = 0$ in figures 9(a) and (b). The mechanism of $N_{ep}$ affecting the OAM spectrum distribution is as follows: the increase of the electron density increases the angular plasma density, reduces the real part and increases the imaginary part of the dielectric tensor element of magnetized plasma, changes the wave number of the characteristic waves in the magnetized plasma, and finally affects the spatial distribution of the phase. The phase disturbance can be characterized by the distribution of OAM states.

Figure 9. OAM state distribution of the transmitted beam for coaxial vortex beams with (a)

$l = - 2, 2,$ (b)

$l = - 3, - 1, 1, 3$ in different plasma slabs.

DownLoad: Full-Size Img PowerPoint

4. Conclusion

A new method to solve the transmission of a THz vortex beam passing through a multi-layered magnetized plasma slab has been provided and extended to the case of multiple coaxial vortex beams incidence. By analyzing the transmission of single and multiple x-polarized Bessel vortex beams through an inhomogeneous magnetized plasma slab with double Gaussian distribution, the major conclusions are obtained as follows. (1) The propagation of a linear polarized vortex beam in magnetized plasma had similar Faraday rotation effects to that of plane waves. With the increase of $B_{0},$ the magnitude of $E_{y}^{tra}$ increased, and the trajectories of the transmitted electric field vector gradually evolved to circular. (2) Owing to the tiny impact of the increase of $B_{0}$ on the permittivity tensor, the phase disturbance caused by the variation of $B_{0}$ was weak, so that the distribution of OAM state had not changed. (3) The increase of the incident angle made the plasma attenuation increase, and the phase distribution of the wavefront was seriously deviated from the incident beam, therefore, the contours and the OAM state of the transmitted beam were distorted dramatically. (4) The increase of the maximum value $N_{ep}$ of the plasma density distorted the profiles of the transmitted field and the distribution of the OAM state severely. On average, the conclusions of the above-mentioned analysis are helpful to qualitatively evaluate the influence of plasma parameters on communication quality, and have theoretical reference value for selecting appropriate beam parameters in THz OAM communication. Meanwhile, it is noteworthy that the method provided in this paper can spread to study the propagation characteristics of other shaped beams in magnetized plasma.

Appendix

The coefficient matrices in equations (22)–(24) are

$\begin{matrix} [\begin{matrix} κ_{11}^{0} & κ_{12}^{0} & κ_{13}^{0} & κ_{14}^{0} \\ κ_{21}^{0} & κ_{22}^{0} & κ_{23}^{0} & κ_{24}^{0} \\ κ_{31}^{0} & κ_{32}^{0} & κ_{33}^{0} & κ_{34}^{0} \\ κ_{41}^{0} & κ_{42}^{0} & κ_{43}^{0} & κ_{44}^{0} \end{matrix}] = {[\begin{matrix} 1 & 0 & 1 & 0 \\ 0 & \frac{h}{k_{0}} & 0 & - \frac{h}{k_{0}} \\ h & 0 & - h & 0 \\ 0 & 1 & 0 & 1 \end{matrix}]}^{- 1} & \times [\begin{matrix} α_{11} & α_{12} & - α_{11} & - α_{12} \\ X^{1} & Y^{1} & - X^{1} & - Y^{1} \\ h_{1} α_{11} & h_{2} α_{12} & h_{1} α_{11} & h_{2} α_{12} \\ \frac{k_{1}}{k_{0}} ß_{11} & \frac{k_{2}}{k_{0}} ß_{12} & \frac{k_{1}}{k_{0}} ß_{11} & \frac{k_{2}}{k_{0}} ß_{12} \end{matrix}] \end{matrix}$

(A1)

$\begin{matrix} [\begin{matrix} κ_{11}^{j} & κ_{12}^{j} & κ_{13}^{j} & κ_{14}^{j} \\ κ_{21}^{j} & κ_{22}^{j} & κ_{23}^{j} & κ_{24}^{j} \\ κ_{31}^{j} & κ_{32}^{j} & κ_{33}^{j} & κ_{34}^{j} \\ κ_{41}^{j} & κ_{42}^{j} & κ_{43}^{j} & κ_{44}^{j} \end{matrix}] = {[\begin{matrix} α_{j 1} e^{i h_{j 1} z_{j}} & α_{j 2} e^{i h_{j 2} z_{j}} & - α_{j 1} e^{- i h_{j 1} z_{j}} & - α_{j 2} e^{- i h_{j 2} z_{j}} \\ X_{j} e^{i h_{j 1} z_{j}} & Y_{j} e^{i h_{j 2} z_{j}} & - X_{j} e^{- i h_{j 1} z_{j}} & - Y_{j} e^{- i h_{j 2} z_{j}} \\ α_{j 1} h_{1}^{j} e^{i h_{j 1} z_{j}} & α_{j 2} h_{j 2} e^{i h_{j 2} z_{j}} & α_{j 1} h_{1}^{j} e^{- i h_{j 1} z_{j}} & α_{j 2} h_{j 2} e^{- i h_{j 2} z_{j}} \\ k_{j 1} ß_{j 1} e^{i h_{j 1} z_{j}} & k_{j 2} ß_{j 2} e^{i h_{j 2} z_{j}} & k_{j 1} ß_{j 1} e^{- i h_{j 1} z_{j}} & k_{j 2} ß_{j 2} e^{- i h_{j 2} z_{j}} \end{matrix}]}^{- 1} \\ \times [\begin{matrix} α_{(j + 1) 1} e^{i h_{(j + 1) 1} z_{j}} & α_{(j + 1) 2} e^{i h_{(j + 1) 2} z_{j}} & - α_{(j + 1) 1} e^{- i h_{(j + 1) 1} z_{j}} & - α_{(j + 1) 2} e^{- i h_{(j + 1) 2} z_{j}} \\ X_{j + 1} e^{i h_{(j + 1) 1} z_{j}} & Y_{j + 1} e^{i h_{(j + 1) 2} z_{j}} & - X_{j + 1} e^{- i h_{(j + 1) 1} z_{j}} & - Y_{j + 1} e^{- i h_{(j + 1) 2} z_{j}} \\ α_{(j + 1) 1} h_{(j + 1) 1} e^{i h_{(j + 1) 1} z_{j}} & α_{(j + 1) 2} h_{(j + 1) 2} e^{i h_{(j + 1) 2} z_{j}} & α_{(j + 1) 1} h_{(j + 1) 1} e^{- i h_{(j + 1) 1} z_{j}} & α_{(j + 1) 2} h_{2}^{j + 1} e^{- i h_{(j + 1) 2} z_{j}} \\ k_{(j + 1) 1} ß_{(j + 1) 1} e^{i h_{(j + 1) 1} z_{j}} & k_{(j + 1) 2} ß_{(j + 1) 2} e^{i h_{(j + 1) 2} z_{j}} & k_{(j + 1) 1} ß_{(j + 1) 1} e^{- i h_{(j + 1) 1} z_{j}} & k_{(j + 1) 2} ß_{(j + 1) 2} e^{- i h_{(j + 1) 2} z_{j}} \end{matrix}] \end{matrix}$

(A2)

$\begin{matrix} [\begin{matrix} κ_{11}^{n} & κ_{12}^{n} & κ_{13}^{n} & κ_{14}^{n} \\ κ_{21}^{n} & κ_{22}^{n} & κ_{23}^{n} & κ_{24}^{n} \\ κ_{31}^{n} & κ_{32}^{n} & κ_{33}^{n} & κ_{34}^{n} \\ κ_{41}^{n} & κ_{42}^{n} & κ_{43}^{n} & κ_{44}^{n} \end{matrix}] & = {[\begin{matrix} α_{N 1} e^{i h_{N 1} z_{N}} & α_{N 2} e^{i h_{N 2} z_{N}} & - α_{N 1} e^{- i h_{N 1} z_{N}} & - α_{N 2} e^{- i h_{N 2} z_{N}} \\ X_{N} e^{i h_{N 1} z_{N}} & Y_{N} e^{i h_{N 2} z_{N}} & - X_{N} e^{- i h_{N 1} z_{N}} & - Y_{N} e^{- i h_{N 2} z_{N}} \\ α_{N 1} h_{N 1} e^{i h_{N 1} z_{N}} & α_{N 2} h_{N 2} e^{i h_{N 2} z_{N}} & α_{N 1} h_{1}^{n} e^{- i h_{N 1} z_{N}} & α_{N 2} h_{2}^{n} e^{- i h_{N 2} z_{N}} \\ \frac{k_{N 1}}{k_{0}} ß_{N 1} e^{i h_{N 1} z_{N}} & \frac{k_{N 2}}{k_{0}} ß_{N 2} e^{i h_{N 2} z_{N}} & \frac{k_{N 1}}{k_{0}} ß_{N 1} e^{- i h_{N 1} z_{N}} & \frac{k_{N 2}}{k_{0}} ß_{N 2} e^{- i h_{N 2} z_{N}} \end{matrix}]}^{- 1} \end{matrix}$

(A3)

where $X_{j} = ß_{j 1} h_{j 1} / k_{j 1} - i γ_{j 1}, Y_{j} = ß_{j 2} h_{j 2} / k_{j 2} - i γ_{j 2},$ and $α_{j q}, ß_{j q}, γ_{j q}$ are presented in equation (12).

Acknowledgments

This study was funded and supported by the Comprehensive Research Facility for Fusion Technology Program of China (No. 2018-000052-73-01-001228), the HFIPS Director’s Fund (No. YZJJKX202301), Anhui Provincial Major Science and Technology Project (No. 2023z020004) and Task JB22001 from the Anhui Provincial Department of Economic and Information Technology.

References (16)

References

[1]	Boboc A et al 2010 Rev. Sci. Instrum. 81 10D538 doi: 10.1063/1.3478146
[2]	Fischer R et al 2010 Fusion Sci. Technol. 58 675 doi: 10.13182/FST10-110
[3]	Shibata Y et al 2011 Plasma Fusion Res. 6 1302136 doi: 10.1585/pfr.6.1302136
[4]	Van Zeeland M A et al 2018 Rev. Sci. Instrum. 89 10B102 doi: 10.1063/1.5037461
[5]	Zhuang G et al 2016 J. Instrum. 11 C02026 doi: 10.1088/1748-0221/11/02/C02026
[6]	Jie Y X et al 2000 Int. J. Infrared Millimeter Waves 21 1375 doi: 10.1023/A:1026496705525
[7]	Shi N et al 2011 Plasma Sci. Technol. 13 347 doi: 10.1088/1009-0630/13/3/14
[8]	Derafshi Z H, Frounchi J and Taghipour H 2010 A high speed FPGA implementation of a 1024-point complex FFT processor In: Proceedings of the Second International Conference on Computer and Network Technology Bangkok: IEEE 2010: 312
[9]	Huang X H, Wang Z H and Chou G Q 2007 New method of estimation of phase, amplitude, and frequency based on all phase FFT spectrum analysis In: Proceedings of the International Symposium on Intelligent Signal Processing and Communication Systems Xiamen: IEEE 2007: 284
[10]	Gil C et al 2008 Rev. Sci. Instrum. 79 10E710 doi: 10.1063/1.2956824
[11]	Ito Y et al 2005 Fusion Eng. Des. 74 847 doi: 10.1016/j.fusengdes.2005.06.241
[12]	Juhn J W et al 2010 Rev. Sci. Instrum. 81 10D540 doi: 10.1063/1.3492381
[13]	Mlynek A et al 2014 Rev. Sci. Instrum. 85 11D408 doi: 10.1063/1.4890574
[14]	Yao Y et al 2017 J. Instrum. 12 C12043 doi: 10.1088/1748-0221/12/12/C12043
[15]	Ding B G et al 2015 Plasma Sci. Technol. 17 837 doi: 10.1088/1009-0630/17/10/05
[16]	Lv S S et al 2021 Sensors 21 4463 doi: 10.3390/s21134463

[1]	Tianyuan HUANG (黄天源), Peiyu JI (季佩宇), Jianjun HUANG (黄建军), Bin YU (于斌), Xuemei WU (吴雪梅). Modification of exposure conditions by the magnetic field configuration in helicon antenna-excited helium plasma[J]. Plasma Science and Technology, 2021, 23(1): 15403-015403. DOI: 10.1088/2058-6272/abcb55
[2]	Yashika GHAI, Nimardeep KAUR, Kuldeep SINGH, N S SAINI. Dust acoustic shock waves in magnetized dusty plasma[J]. Plasma Science and Technology, 2018, 20(7): 74005-074005. DOI: 10.1088/2058-6272/aab491
[3]	Ruggero BARNI, Stefano CALDIROLA, Luca FATTORINI, Claudia RICCARDI. Tomography of a simply magnetized toroidal plasma[J]. Plasma Science and Technology, 2018, 20(2): 25102-025102. DOI: 10.1088/2058-6272/aa9028
[4]	Neda SHAMSIAN, Babak SHIRANI BIDABADI, Hosein PIRJAMADI. Development of a radiographic method for measuring the discrete spectrum of the electron beam from a plasma focus device[J]. Plasma Science and Technology, 2017, 19(7): 75101-075101. DOI: 10.1088/2058-6272/aa632e
[5]	HU Guanghai (胡广海), JIN Xiaoli (金晓丽), YUAN Lin (袁林), ZHANG Qiaofeng (张乔枫), XIE Jinlin (谢锦林), LI Hong (李弘), LIU Wandong (刘万东). Oxide Coated Cathode Plasma Source of Linear Magnetized Plasma Device[J]. Plasma Science and Technology, 2016, 18(9): 918-923. DOI: 10.1088/1009-0630/18/9/08
[6]	YUAN Yuan(袁媛), JIANG Zhonghe(江中和), GUO Weixin(郭伟欣), SUN Xinfeng(孙新锋), HU Xiwei(胡希伟). Mode-Coupling Analysis of Parametric Decay Instability in Magnetized Plasmas[J]. Plasma Science and Technology, 2014, 16(9): 809-814. DOI: 10.1088/1009-0630/16/9/01
[7]	ZHAO Qing(赵青), XING Xiaojun(邢晓俊), XUAN Yinliang(宣银良), LIU Shuzhang(刘述章). The Influence of Magnetic Field on Antenna Performance in Plasma[J]. Plasma Science and Technology, 2014, 16(6): 614-619. DOI: 10.1088/1009-0630/16/6/14
[8]	LI Jiajia, HU Zhanghu, SONG Yuanhong, WANG Younian. Effects of Fast-Ion Injection on a Magnetized Sheath near a Floating Wall[J]. Plasma Science and Technology, 2013, 15(1): 1-6. DOI: 10.1088/1009-0630/15/1/01
[9]	HU Zhidan(胡志丹), SHENG Zhengming (盛政明), Ding Wenjun (丁文君), WANG Weimin (王伟民), DONG Quanli (董全力), ZHANG Jie(张杰), et al. Electromagnetic Emission from Laser Wakefields in Magnetized Underdense Plasmas[J]. Plasma Science and Technology, 2012, 14(10): 874-879. DOI: 10.1088/1009-0630/14/10/04
[10]	ZHANG Shoubiao, XIE Jinlin, HU Guanghai, WANG Fusheng, LI Hong, LIU Wandong, GAO Xiang. An Experimental Platform for the Study on Magnetic Reconnection in Laboratory Plasmas[J]. Plasma Science and Technology, 2011, 13(3): 286-289.

Cited By

Cited by

Periodical cited type(6)

1.	Guo, S., Cen, H., Ouyang, W. et al. Rapid prediction model of terahertz transmission in hypersonic plasma sheath under different flight speeds for different vehicle types. Journal of Physics D: Applied Physics, 2025, 58(8): 085208. DOI:10.1088/1361-6463/ad9dfb
2.	Ma, X., Zeng, Y., Wu, G. et al. Detection of Radar Target Wrapped in Plasma Sheath Based on Intrapulse Energy Focusing. IEEE Transactions on Aerospace and Electronic Systems, 2025. DOI:10.1109/TAES.2025.3528912
3.	Peng, K., Yuan, K., Mao, M. et al. The cause and the evolution of “lens” structure for terahertz signals in hypersonic plasma sheaths. AIP Advances, 2024, 14(7): 075018. DOI:10.1063/5.0219092
4.	Bai, B., Ding, Y., Li, X. et al. Research on the Spatial Distribution Mechanism of the EM Reflection Intensity of a Plasma-Sheath-Covered Target. IEEE Transactions on Plasma Science, 2023, 51(4): 1052-1058. DOI:10.1109/TPS.2022.3224651
5.	Ding, Y., Bai, B., Li, X. et al. Research on EM Shielding Mechanism of the Plasma-Sheath-Covered Target. IEEE Transactions on Plasma Science, 2023, 51(3): 632-640. DOI:10.1109/TPS.2022.3217096
6.	Ouyang, W., Liu, Q., Mao, W. et al. Plasma and Terahertz-Wave Propagation Characteristics in Atmospheric-Pressure Plasma Jet Under Different Operating Conditions. IEEE Transactions on Antennas and Propagation, 2023, 71(3): 2710-2718. DOI:10.1109/TAP.2023.3237284