Effect of ion stress on properties of magnetized plasma sheath

1.   Introduction

With the development of modern science and technology, terahertz (THz) waves that have not been fully exploited and utilized have attracted more and more people's attention [1, 2]. Today, THz technology is widely used in medical diagnosis [3, 4], space literature [5], security detection [6–8] and military fields [9]. However, a great leap forward in THz technology can only be achieved by finding efficient THz radiation sources and detectors [10, 11]. In a recent paper by Wang et al, theoretical research and analysis based on the nonlinear optical difference frequency principle show that the THz radiation source with optical parameter effect can effectively generate the THz signal [12]. As early as nearly 30 years ago, Dyakonov and Shur et al had shown that stable DC currents in ballistic field-effect transistors (FETs) could be unstable against the generation of high-frequency plasma waves due to the wave reflection from the device boundaries, which provides a new mechanism of electromagnetic radiation [13]. Aizin et al extended the hydrodynamic model to the case considering stubs, which proved that the instability of plasma in FETs can be controlled by optimizing boundary conditions, the velocity of plasma waves, etc [14].

Plasma is an ionized gas composed of charged particles as well as neutral particles. The plasma has a very high electrical conductivity, and a strong coupling effect with the electromagnetic field. Plasma waves are caused by a wide variety of collective motion modes. There are elastic restoring forces, thermal pressure, electrostatic force and magnetic field force in plasmas. Because of the interaction among different forces, plasma waves include electromagnetic waves (shear waves) caused by electromagnetic force, static waves (longitudinal waves) caused by electrostatic force, and sound waves caused by thermal pressure. The plasma wave traveling in the FET channel is reflected and amplified when it spreads to the drain. Then the plasma wave becomes unstable, which leads to electromagnetic radiation. When the gate length reaches to submicron, the frequencies of electromagnetic radiation belong to the THz range [12].

Spin is an intrinsic motion caused by the intrinsic angular momentum of particles, which is an intrinsic property of particles [15]. The contribution of intrinsic spin to the ponderomotive force was explored in a magnetized plasma by using an expression containing both the classical as well as the spin-induced ponderomotive force and the results show that a strong spin polarization in a plasma can be induced by the spin-ponderomotive force [16]. In a recent paper by Brodin [17], the properties of solitary structures were investigated taking into account the microscopic spin properties of the electrons and it is found that the quantum spin effects will induce the solitary structures. Lin et al analyzed the effect of the spin effect on the propagation of electromagnetic waves in any direction in magnetized plasmas, and found that the spin effect has a significant effect on low-frequency modes [18]. A fully consistent quantum-mechanical description of a spin-polarized plasma was constructed by using a generalized matrix Wigner function [19]. There has been many investigations on single-particle spin effects in external field models, for example, intense laser field [20]. In recent years, the study on the excitation of collective modes in spin systems has also aroused curiosity in the scientific community [21, 22]. It is shown by previous studies that there are many novel features and certain advantages in THz radiation in the presence of the spin effect of electrons in the channel of nano-FET.

Experimental and numerical research has demonstrated that the THz plasma waves become unstable in the presence of the asymmetric boundary conditions that the impedance between the gate and the drain approaches infinity, and the impedance between the gate and the source tends to zero Ohm. In fact, there are some difficulties in the realization of asymmetric boundary conditions, which leads to the attention of other boundary conditions. The instability of THz plasma waves was developed using a compact model with tunable narrow-channel regions with an increased width, and the plasmonic boundary conditions are matched at the interfaces between different regions in the channel [14]. In 2018 [23], by linearization of hydrodynamic equations describing the collective behavior of two-dimensional electron gas in the channel of high electron mobility transistor, the THz radiation was firstly explored with the impact of non-ideal termination impedances at the source and the drain, and numerical results showed that the DS instability would be affected by non-ideal boundary conditions. However, to the authors' best knowledge, very few publications address the issue of the characteristics of THz plasma waves traveling in the channel of FETs with non-ideal boundary conditions, spin effects and quantum effects.

In our work, the instability of THz plasma waves in FETs under the non-ideal boundary conditions is investigated with the joint action of spin effects, quantum effects and the external magnetic field. In order to study the features of the THz plasma waves, we started from the quantum fluid dynamics equation and spin equation, and obtained dispersion relation with non-ideal boundary conditions and quantum effects. The influences of drain capacitance, source capacitance, spin effects, quantum effect, channel width and electron drift velocity on the instability of plasma THz wave in FETs under the joint action of non-ideal boundary conditions are discussed in detail.

2.   Theoretical model

With the rapid development of nanotechnology, the characteristic size of electronic devices is becoming smaller and smaller, that is, from the micron to nanometer. The characteristic size of traditional semiconductor devices is in the order of microns, and the particle property of electrons is dominant. When the characteristic size of the device is equal to or smaller than the de Broglie wavelength of electrons or the mean free path of electrons, and the volatility of electrons is dominant, and then we get all sorts of quantum effects. At Fermi temperature, the electron state changes from the Maxwell–Boltzmann distribution to the Fermi–Dirac distribution [24]. In this case, the motion of the electron must be described using the equations of quantum fluid mechanics with quantum effects. So, in order to investigate the instability of THz plasma waves with spin effects and quantum effects, the collective behavior of electron fluid is governed by the quantum fluid dynamics equations [25, 26] with spin effects [27, 28] and Maxwell's equations:

where $n$ represents the number density of the electrons, $\stackrel{⃑}{v}$ is the fluid velocity of the electrons, $e$ represents the charge of the electron, $m$ represents the mass of the electron; $\stackrel{⃑}{E}$ and $\stackrel{⃑}{B}$ represent the electric field strength and the magnetic field strength, respectively; $\mu$ represents Bohr magnetic moment, the smallest element of the Bohr magnetic moment is the Bohr magneton. $P$ stands for quantum Fermi pressure; ${\mu }_0$ and ${\mathit{ϵ}}_0$ are used to represent the permeability and electrical conductivity in vacuum, respectively; and the current density $J$ is equal to $\stackrel{⃑}{J}=-e{n}_0\stackrel{⃑}{v};$ electron spin-magnetization current density ${\stackrel{⃑}{J}}_{\mathrm{m}}=\mathrm{\nabla }\times \stackrel{⃑}{M}$ (here $\stackrel{⃑}{M}$ is the spin-magnetization vector in the macroscopic state.), and $\stackrel{⃑}{M}=\frac{n_{\mathrm{e}}{u}_{\mathrm{B}}^2}{k_{\mathrm{B}}{T}_{\mathrm{Fe}}}\stackrel{⃑}{B}$ (here the ${\mu }_{\mathrm{B}}=\frac{eħ}{2m}$ represents the Bohr magneton, $T_{\mathrm{Fe}}$ represents Fermi temperature, $k_{\mathrm{B}}$ indicates the Boltzmann constant, $n_{\mathrm{e}}$ represents the electron gas density at the Fermi temperature.).

The evolution equation describing the spin effects in quantum plasma is as follows:

Because the scale of the electron is larger than the Larmor radius of the electron in the magnetohydrodynamic model, we can ignore the quadratic term of the spin $\stackrel{⃑}{S}$ in the spin evolution equation. For the lower order, when it is far below the cyclotron frequency of the electron, the inertia term of the electron can also be ignored and $\stackrel{⃑}{S}+\stackrel{⃑}{B}=0$ can be obtained, so we can get the spin expression:

here, $\eta \left(x\right)=\tanh \left(x\right)\left(x=\frac{{\mu }_{\mathrm{B}}{B}_0}{K_{\mathrm{B}}{T}_{\mathrm{e}}}\right)$ is the Langevin parameter, which is a physical quantity describing the size of the spin. We see from the spin formula that spin is a function of the Langevin parameter. The Langevin parameter depends on the magnetic field, which is a variable when the magnetic field changes. $T_{\mathrm{e}}$ represents the electron temperature, $B_0$ is the strength of the external magnetic field.

To study the stability of THZ plasma waves, we set $n=n_0+n_1,$ $v_x=v_0+v_{x_1},$ $v_y=v_{y_1},$ $E_x=E_{x_1},$ $E_y=E_{y_1},$ $B=B_0+B_{z1},$ $n_0,$ $v_0,$ $v_0$ are constants. In order to linearize equations (1)–(6), we set $n_1,$ $v_{x1},$ $v_{y1}$ ~ $\exp \left(-\mathrm{i}\omega t+\mathrm{i}kx+\mathrm{i}qy\right),$ where $k$ and $q$ are the wave vectors along the X-axis and Y-axis directions, respectively. The linearized equations can be written as:

where ${\omega }_{\mathrm{ce}}=\frac{e{B}_0}{m},$ $v_{\mathrm{Fe}}=\sqrt{\frac{2T_{\mathrm{Fe}}}{m}},$ $B_{z1}={\displaystyle \frac{1}{\omega }}\left(k{E}_{y1}-q{E}_{x1}\right).$

To study the features of THz plasma waves simply, the linearized expressions of equations (9)–(13) are as follows,

where $K=kL,$ $Q=qL,$ $\mathrm{\Omega }=\frac{\omega L}{v_{\mathrm{Fe}}},$ $\beta =\frac{v_0}{v_{\mathrm{Fe}}},$ ${\tilde{v}}_{x1}=\frac{v_{x1}}{v_{\mathrm{Fe}}},$ ${\tilde{v}}_{y1}=\frac{v_{y1}}{v_{\mathrm{Fe}}},$ ${\tilde{\omega }}_{\mathrm{ce}}=\frac{{\omega }_{\mathrm{ce}}L}{v_{\mathrm{Fe}}},$ ${\tilde{E}}_{x1}=E_{x1}\frac{eL}{m{v}_{\mathrm{Fe}}^2},$ ${\tilde{E}}_{y1}=E_{y1}\frac{eL}{m{v}_{\mathrm{Fe}}^2},$ $H=\frac{{ħ}^2}{4m^2}\frac{1}{v_{\mathrm{Fe}}^2{L}^2},$ ${\tilde{n}}_1=\frac{n_1}{n_0},$ $\tilde{c}=\frac{c}{v_{\mathrm{Fe}}},$ $\chi \equiv 1-\frac{n_0{\mu }_0{\mu }_{\mathrm{B}}^2}{k_{\mathrm{B}}{T}_{\mathrm{e}}},$ ${\omega }_{\mathrm{pe}}={\frac{e^2{n}_0}{m{\epsilon }_0}}^{\frac{1}{2}},$ $\widetilde{{\omega }_{\mathrm{pe}}}=\frac{{\omega }_{\mathrm{pe}}L}{v_{\mathrm{Fe}}}.$

According to equations (14)–(18), we obtain the dimensionless dispersion relation:

In equation (19), $\mathrm{\Omega }=\mathrm{\Omega }\text{'}+\mathrm{i}\mathrm{\Omega }\text{'}\text{'},$ $\mathrm{\Omega }\text{'}$ represents the plasma oscillation frequency, and $\mathrm{\Omega }\text{'}\text{'}$ represents the instability gain of the THz plasma waves. If the value of the instability gain is greater than zero, the THz plasma waves are unstable. Instead, if the value of the instability gain is less than zero, the THz plasma waves are stable. Furthermore, $K$ represents the dimensionless wave vector in the X-direction, and $Q$ represents the dimensionless wave vector in the Y-direction, $H$ represents the contribution of the quantum effect, $\widetilde{{\omega }_{\mathrm{pe}}}$ is the dimensionless plasma frequency. $\eta$ is the Langevin parameter, which represents the spin distribution under thermal equilibrium. ${\tilde{\omega }}_{\mathrm{ce}}$ is used to characterize the dimensionless electron cyclotron frequency of the magnetic field.

From equation (19), we can obtain the seven roots of $K.$ Since the waves propagating from the source to the drain will be reflected, the incident wave and the reflected wave are mutually symmetrical, so we should choose two symmetric roots $K_1$ and $K_2.$ The solution of linearized equation (9) takes the form:

here, $I_{\omega }=-eW\left(v_0{n}_1+n_0{v}_{x1}\right)$ is the channel current, the voltage distribution is $V_{\omega }=\delta {\phi }_{\omega },$ $W$ is used to denote the width of the channel.

DS instable theory indicated that the THz plasma waves become unstable in the presence of the asymmetric boundary conditions that the impedance between the gate and the drain approaches infinity, and the impedance between the gate and the source tends to zero Ohm. Then, Mona et al have come up with other boundary conditions that source capacitance and drain capacitance are attached in high electron mobility transistors [23]. These boundary conditions agree with equations $V\left(0,t\right)=-\mathrm{i}{Z}_{\mathrm{S}}I\left(0,t\right)$ and $V\left(L,t\right)=\mathrm{i}{Z}_{\mathrm{D}}I\left(L,t\right),$ which are also called non-ideal boundary conditions. (Here $Z_{\mathrm{S}}$ and $Z_{\mathrm{D}}$ indicate source and drain impedance, respectively.) The non-ideal boundary conditions can be written in the form:

here $C_{\mathrm{S}}$ is the source capacitance, $C_{\mathrm{D}}$ is the drain capacitance, and $L$ is the channel length.

According to equations (20)–(23), we obtain:

We shall write the above expression as:

here, $K_1=k_{1}L,$ $K_2=k_{2}L,$ $\widetilde{C_{\mathrm{S}}}={\displaystyle \frac{C_{\mathrm{S}}}{C{L}^2}},$ $\widetilde{C_{\mathrm{D}}}={\displaystyle \frac{C_{\mathrm{D}}}{C{L}^2}}.$

3.   Numerical simulation

In order to investigate the propagation of THz plasma waves in the channel of FET with the combined contribution of spin effects, the quantum effects and non-ideal boundary conditions, we numerically solve equations (19) and (25). The features of instability of THz plasma waves are discussed in two different cases.

Case 1 Source capacitance is greater than the drain capacitor

Figure 1 shows the trend of the oscillation frequency $\mathrm{\Omega }\text{'}$ and the instability gain $\mathrm{\Omega }\text{'}\text{'}$ of the THz plasma waves changing with the drift velocity $\beta$ for different $C_{\mathrm{S}}$ and $C_{\mathrm{D}},$ when the source capacitance $C_{\mathrm{S}}$ is larger than the gate capacitance $C_{\mathrm{D}}.$ As shown in figures 1(a) and (c), the plasma oscillation frequency decreases with increasing drift velocity $\beta ,$ and when the drift velocity $\beta$ is certain, the oscillation frequency of the THz plasma waves decreases with the increasing source capacitance $C_{\mathrm{S}}$ and drain capacitance $C_{\mathrm{D}}.$ From figure 1(b), it can be seen that the instability gain increases with the drift velocity $\beta ;$ when the drift velocity $\beta$ is a fixed value, the larger the source capacitance $C_{\mathrm{S}},$ the greater the instability gain increases; and the effect of $C_{\mathrm{S}}$ on the instability gain is not obvious when $\beta$ is small, however, as $\beta$ continuously increases, the effect of $C_{\mathrm{S}}$ on the instability gain becomes increasingly apparent. As can be seen from figure 1(d), the instability gain of the THz plasma waves increases as $\beta$ grows. When $\beta$ is a fixed value, the larger the drain capacitance $C_{\mathrm{D}},$ the smaller the instability gain. As we know, plasma waves propagating along the channel of FET are reflected from the drain and amplified, then become unstable. When plasma waves propagate along the channel, there will be a change in energy due to the variation of physical factors. The energy loss will reduce the instability of plasma waves. So, we can conclude that the increase in source capacitance will reduce the energy loss, and the increase in drain capacitance will strengthen the energy loss.

Figure  1.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $\beta$ for different $C_{\mathrm{S}},$ $C_{\mathrm{D}}$ with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,\eta =0.1,\mu$ $=0.4,{\omega }_{\mathrm{pe}}=0.6,H=0.001,W=30,C_{\mathrm{D}}=1,C_{\mathrm{S}}=20.$

DownLoad: Full-Size Img PowerPoint

The characteristics of FET will change in the presence of the quantum effects. Figure 2 shows the effect of quantum effects $H$ on the oscillation frequency and the instability gain of the THz plasma waves when the source capacitance $C_{\mathrm{S}}$ is larger than the drain capacitor $C_{\mathrm{D}}.$ It can be seen from figure 2(a) that the oscillation frequency increases with quantum effect. Figure 2(b) shows that the instability gain increases with increasing drift velocity $\beta ,$ when the quantum effects $H$ is certain; and if the drift velocity $\beta$ is a constant value, the larger the quantum effects $H,$ the greater the instability gain. The increase of instability gain represents the increase of radiant power. It is clear that quantum effects can reduce the energy loss of the plasma during the oscillation. As the size of electronic devices gets smaller and smaller, quantum effects become more pronounced. We can change the oscillation frequency and instability gain of the THz plasma waves by changing the size of the device. This is consistent with the results of Knap et al, that the oscillation frequency increases as the device size decreases [29].

Figure  2.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $\beta$ for different $H$ with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,\eta =0.1,\mu$ $=0.4,{\omega }_{\mathrm{pe}}=0.6,W=30,C_{\mathrm{S}}=20,C_{\mathrm{D}}=1.$

DownLoad: Full-Size Img PowerPoint

Figure 3 shows the effect of the channel width $W$ on the oscillation frequency and the instability gain when the source capacitance $C_{\mathrm{S}}$ is larger than the drain capacitor $C_{\mathrm{D}}.$ From figure 3(a), we can find that the plasma oscillation frequency increases with the channel width. During our discussion, we found that the influence of the channel width $W$ on the instability gain is complicated. Such cases are shown in figure 3(b). Figure 3(b) depicts that for a certain $\beta ,$ the instability gain increases firstly and after crossing a critical value in terms of $\beta$ it stars decreasing gradually. From the above results, it can be obtained that the current will become unstable in the channel of FET with non-ideal boundary conditions and there is the emission of THz waves when the source capacitance is larger than the drain capacitor.

Figure  3.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $W$ and $\beta$ with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,\eta =0.1,\mu$ $=0.4,{\omega }_{\mathrm{pe}}=0.6,H=0.001,C_{\mathrm{S}}=20,C_{\mathrm{D}}=1.$

DownLoad: Full-Size Img PowerPoint

Figure 4 shows the effect of the spin effect on the oscillation frequency and the instability gain when the source capacitance is larger than the drain capacitor. We can see from figure 4 that the oscillation frequency and instability gain with the spin effect are larger than that without the spin effect. That is, the spin can not only enhance the oscillation frequency but also increase the instability gain. Obviously, the presence of the spin reduces the energy loss during the plasma oscillation and amplification.

Figure  4.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $\beta$ with and without spin effects with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,H=0.001,{\omega }_{\mathrm{pe}}$ $=0.6,W=30,C_{\mathrm{S}}=20,C_{\mathrm{D}}=1.$

DownLoad: Full-Size Img PowerPoint

Since $\eta$ affects the spin value, we change $\eta$ to observe the change of oscillation frequency and instability gain. The results show that the size of $\eta$ has a very small effect on the oscillation frequency and the instability. Compared with the results in [30], the oscillation frequency and instability gain of THz plasma wave increase with increasing $\eta .$ In other words, under non-ideal boundary conditions, the oscillation frequency and instability gain of THz plasma waves increase with the appearance of spin, but the size of spin does not have a significant effect on the oscillation frequency and instability gain. However, under asymmetric boundary conditions, the effect of the size of the spin on the oscillation frequency and instability gain is evident.

Case 2 Drain capacitor is greater than the source capacitance.

Figures 5–8 give an example of the instability of THz plasma waves in the channel of FET with spin effects, quantum effects and non-ideal boundary conditions when the source capacitance $C_{\mathrm{S}}$ is smaller than the drain capacitance $C_{\mathrm{D}}.$ As can be seen from figures 5–8, in two cases, that is, the source capacitance is smaller than the drain capacitance, and the source capacitance is bigger than the drain capacitance, the influence of spin effects, quantum effects, the drain capacitance, the source capacitance and the width of the channel on the plasma oscillation frequency is the same. As shown in figures 5(b), 6(b), 7(b) and 8(b) when the source capacitance $C_{\mathrm{S}}$ is less than the drain capacitance $C_{\mathrm{D}},$ the value of the resulting instability gain $\mathrm{\Omega }\text{'}\text{'}$ is less than zero, which indicates that the THz wave oscillation is stable in this case. The larger the absolute value of $\mathrm{\Omega }\text{'}\text{'}$ means the more stable THz plasma waves. From figure 6(b), we can see that the value of the instability gain is negative when the source capacitance is smaller than the drain capacitor, indicating that the oscillation of the THz wave is stable in this case. As the quantum effect increases, the instability gain becomes small and the THz waves become more stable. That is, in this case, THz radiations do not occur. When the source capacitance increases the instability of THz plasma waves slows down.

Figure  5.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $\beta$ for different $C_{\mathrm{S}},$ $C_{\mathrm{D}}$ with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,\eta =0.1,\mu$ $=0.4,{\omega }_{\mathrm{pe}}=0.6,H=0.001,W=30,C_{\mathrm{D}}=20,C_{\mathrm{S}}=1.$

DownLoad: Full-Size Img PowerPoint

Figure  6.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $\beta$ for different $H$ with ( $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,\eta =0.1,\mu$ $=0.4,{\omega }_{\mathrm{pe}}=0.6,W=30,C_{\mathrm{S}}=1,C_{\mathrm{D}}=20$ ).

DownLoad: Full-Size Img PowerPoint

Figure  7.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $W$ and $\beta$ with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,\eta =0.1,\mu$ $=0.4,{\omega }_{\mathrm{pe}}=0.6,H=0.001,C_{\mathrm{S}}=1,C_{\mathrm{D}}=20.$

DownLoad: Full-Size Img PowerPoint

Figure  8.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $\beta$ with and without spin effects with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.3,c=1000,H=0.001,{\omega }_{\mathrm{pe}}$ $=0.6,W=30,C_{\mathrm{S}}=1,C_{\mathrm{D}}=20.$

DownLoad: Full-Size Img PowerPoint

Similarly, as shown in figures 7(b) and 8(b), the effects of the channel width $W$ and spin on the instability gain are indicated, and the instability gain in this case is less than zero, which indicates that the THz wave oscillation is stable in this case.

Under asymmetric boundary conditions, the source impedance $Z_{\mathrm{S}}=0$ and drain impedance $Z_{\mathrm{D}}=\infty .$ However, under non-ideal boundary conditions, the $Z_{\mathrm{S}}$ and $Z_{\mathrm{D}}$ take the finite values. These impedances are the result of the presence of capacitance at the source or drain. The larger induced impedance of the bonding wires is effectively shunted by these impedances, and defines the termination impedance at both ends of the plasma cavity [31]. Figure 9 shows the effect of two different boundary conditions, namely, asymmetric boundary conditions and non-ideal boundary conditions, on the oscillation frequency and the instability gain of the plasma THz waves. From figure 9(a), we can see that the plasma oscillation frequency in the non-ideal boundary conditions is larger than that in the asymmetric boundary conditions. It is clear that the change in the source and drain impedance caused the change in the oscillation frequency. As shown in figure 9(b), the instability gain in asymmetric boundary conditions is higher than that in the non-ideal boundary conditions. Thus, under non-ideal boundary conditions, finite source and drain impedances suppress instability. That is, non-ideal boundary conditions can increase the oscillation frequency of the plasma THz waves, while asymmetric boundary conditions can improve the instability gain of the plasma THz waves.

Figure  9.  Trend diagram of $\mathrm{\Omega }\text{'}$ and $\mathrm{\Omega }\text{'}\text{'}$ versus $\beta$ for asymmetric boundary conditions and non-ideal boundary conditions with $q=0.1,\chi =0.5,{\omega }_{\mathrm{ce}}=0.1,c=1000,\eta =0.1,\mu =0.4,{\omega }_{\mathrm{pe}}$ $=0.6,H=0.001,W=30,C_{\mathrm{D}}=1,C_{\mathrm{S}}=20.$

DownLoad: Full-Size Img PowerPoint

4.   Conclusion

In conclusion, the linear characteristics of THz plasma waves in the channel of FET with non-ideal boundary conditions are considered by linearizing Maxwell's equations and the quantum fluid dynamics equations containing the spin effect and quantum effect. Summing up the results, it can be concluded that in a FET, if the source capacitance is greater than the drain capacitance, the THz plasma waves become unstable, which indicates that there is THz emission in this case. When the drift velocity is fixed, the oscillation frequency increases with the increase in quantum effect, and decreases with the increase of source capacitance and drain capacitance, while the instability gain increases with the increase of source capacitance and quantum effect, and decreases with the increase of drain capacitance. In contrast, if the drain capacitance is greater than the source capacitance, the value of the instability gain is less than zero, that is, the THz plasma waves are stable. Meanwhile, the quantum effect and spin effect enhance the emission of THz waves, and the boundary conditions can affect the oscillation frequency and the instability gain of plasma THz waves. In [30], the effects of spin, quantum effects, drift velocity and other factors on the oscillation frequency and instability gain of THz plasma waves are investigated under asymmetric boundary conditions. Compared with the results in [30], we found that the influence of boundary conditions is very important. The oscillation frequency with asymmetric boundary conditions is smaller than that under non-ideal boundary conditions, but the instability gain of THz plasma waves becomes lower under non-ideal boundary conditions. The findings suggest that this study could also be useful for the development of THz source and detector.