• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Numerical study on the modulation of THz wave propagation by collisional microplasma photonic crystal

    2020-11-10 03:01:46ShuqunWU吳淑群YuxiuCHEN陳玉秀MingeLIU劉敏格LuYANG楊璐ChaohaiZHANG張潮海andShaobinLIU劉少斌
    Plasma Science and Technology 2020年11期
    關(guān)鍵詞:張潮

    Shuqun WU(吳淑群),Yuxiu CHEN(陳玉秀),Minge LIU(劉敏格),Lu YANG(楊璐),Chaohai ZHANG(張潮海) and Shaobin LIU(劉少斌)

    1 College of Automation Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,People’s Republic of China

    2 College of Electronic and Information Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,People’s Republic of China

    Abstract

    Keywords:microplasma,plasma photonic crystal,terahertz wave,electron density

    1.Introduction

    Terahertz(THz)waves with a wavelength range of 30 μm to 3 mm have attracted considerable attention,due to their promising advanced applications,such as chemical spectroscopy[1],biomedical diagnostics[2],6 G communications[3],and astronomical observations[4].Photonic crystals(PCs)for THz wave modulation have been widely studied for over a decade[5].PCs possess photonic bandgaps(PBGs),within which the propagation of THz waves is not permitted,and which are usually used to confine and guide the desired propagating THz wave with a low degree of loss[6].Once a PC is developed,PBGs dependent on the structure of the PCs are barely capable of being tuned.To enhance the tunability of PBGs for THz wave modulation,plasma is introduced into the PCs,forming an artificial periodic structured composition of plasma and dielectrics,known as a plasma photonic crystal(PPC)[7-9].By electronically controlling the ignition and the physical properties of the plasma,the PBGs of the PPC become reconfigurable and tunable.Compared to other methods based on the external electric field[10],temperature[11],magnetic field[12],and pressure[13],the advantages of the PPC are the wide tunability of plasma permittivity from positive to negative values,and its rapid reconfiguration capability,facilitated by electronically switching the plasma on/off.These advantages of the PPC have been effectively demonstrated in relation to the modulation of the microwave propagation[14-18].

    Microplasma is one promising candidate for the modulation of THz wave propagation in PPCs.In 2015,Askari et al[19]computationally analyzed the THz bandgap of a one-dimensional PPC with collision-free microplasma.They found that the width and the central frequency of the THz stopband gap are strongly dependent on the electron density and the plasma layer thickness.In 2016,PaneerChelvam et al[20]numerically modeled the interactions between a collisionless microdischarge cell and an electromagnetic wave at a frequency of 0.1 THz.It was revealed that due to the epsilon-zero critical density resonance,significant scattering,and the power deposition in the over-dense plasma is observed.In 2017,Yang et al[21]experimentally observed a narrowband attenuation at 157 GHz in a PPC composed of a microplasma jet array at atmospheric pressure.The narrow stopband width was 1 GHz and the maximum time-averaged attenuation of the electromagnetic wave was 5%,representing a surprising deviation from the simulated results.Subsequently,this group further demonstrated a possible application of a 3D PPC in relation to dynamic band-stop filters in the 120-170 GHz region[22].In these previous modeling results for a microplasma photonic crystal(MPPC),the collision frequency for momentum transfer between electrons and atoms or molecules in the PPC is considered to be negligible as compared to the incident wave frequency.However,this is not true for atmospheric pressure microplasmas,where the collision frequency attained is as high as several THz.Tian et al[23]reported that the uniformity of the collision frequency in high-temperature plasma for satellite reentry has a great impact on the THz absorption spectra.Therefore,collision frequency should be taken into consideration in PPC modeling at high pressure for the modulation of THz wave propagation.

    In this work,the numerical modeling of a one-dimensional collisional MPPC is implemented by means of the transfer matrix method(TMM).In addition,the role of the collision frequency in plasma relative permittivity is theoretically analyzed.The effects of electron density,pressure,and plasma width on the THz band structure of MPPCs are also investigated.

    2.Model and theory

    2.1.Numerical model

    Figure 1 shows a schematic of the one-dimensional collisional MPPC.Here,five quartz layers and four uniform microplasma layers are alternated and periodically located.The width(dq)and the dielectric constant(εq)of the quartz layer are 50 μm and 4.2,respectively.The width(dp)of the microplasma layer is 100 μm,which is variable in section 3.3.The lattice constant(d)of the MPPC is the sum of the widths of the quartz layer and the microplasma layer.The number of cycles of the MPPC is 4.A THz wave in transverse electromagnetic(TEM)mode is vertically incident to the MPPC.

    Figure 1.Schematic of a one-dimensional collisional MPPC.Microplasma layer(red),quartz layer(grey).

    2.2.Plasma relative permittivity

    For non-magnetized plasma in weakly ionized gas,the expression for the plasma relative permittivity(εr)can be derived from the well-known Drude model[24]:

    where ω,ωpe,and νmdenote the frequency of the incident electromagnetic wave,the plasma frequency,and the collision frequency for momentum transfer between electrons and atoms or molecules,respectively.In this case,ωpecan be obtained from the plasma equation in[25]:

    where ne,e,me,and ε0represent the electron density in plasma,the electron charge,the electron mass,and the vacuum permittivity,respectively.Here,νmis proportional to the density of the gas,and is dependent on the electron temperature,which is given by[26]

    where nmis the density of background gas(proportional to the gas pressure),p is the gas pressure,and Teis the electron temperature.Here,veis the thermal electron velocity,proportional to the square of the electron temperature,andσtris the momentum transfer collision cross section,which is also dependent on the electron temperature.For a gas pressure of 101 kPa(atmospheric pressure)and an electron temperature of 1 eV,the collision frequency is estimated to be 6.07×1012rad s?1.

    Figure 2.Map of Re(εr)as a function of electromagnetic wave frequency(ω)and electron density.(a)Collision-free plasma,νm=0,(b)collisional plasma,νm=6.07×1012 rad s?1.

    Based on equation(1),the value of the plasma relative permittivity is controllable on the complex plane by changing the plasma frequency and the collision frequency independently.Consequently,the attenuation and the phase shift of the THz waves can be varied.To clarify the role of collision frequency in this process,maps of the real part of the plasma relative permittivity,dependent on the electromagnetic wave frequency and the electron density for the collision-free plasma and the collisional plasma,are plotted in figure 2.For ω in the range of 0.1-0.5 THz,when the electron density increases from 1013cm?3to 1017cm?3,the Re(εr)decreases from 1 to a negative value.The minimums of Re(εr)are?804 for the collision-free plasma in figure 2(a)and?7.51 for the collisional plasma in figure 2(b),respectively.In other words,the addition of the collision frequency compensates for the reduction in Re(εr)for high electron density.

    Figure 3 shows the values of the plasma relative permittivity,controlled by changing the gas pressure and the electron density independently.It is interesting to observe that both Im(εr)and Re(εr)decrease monotonically with the increase in electron density.When the gas pressure increases from 30 kPa to 202 kPa at the fixed electron density of 5×1016cm-3,the Im(εr)first decreases,reaching?8.04 at a gas pressure of 50.5 kPa,increasing as the gas pressure further increases.On the other hand,Re(εr)increases with an increase in the gas pressure.Therefore,both Im(εr)and Re(εr)can be varied by changing the electron density and gas pressure on demand.

    To elucidate plasma-wave coupling based on plasma relative permittivity,a map of the loss tangent,defined by tanδ=Im(εr)/Re(εr),is plotted in figure 4.For an incident wave frequency of 0.1-0.5 THz,as the electron density increases from 1013cm-3to 1017cm-3,the tanδ decreases from 0 to a negative value,then abruptly reverses to a positive value.The transition region(region 1 in figure 4)corresponds to the value of the Re(εr)approaching to zero,which yieldsThis is useful for determining the electron density in plasma.

    2.3.Microplasma photonic crystal

    The MPPC consists of two dielectrics with different permittivities.One of these is microplasma.These two dielectrics are arranged alternately,as shown in figure 1.Due to the Bragger reflection,a photonic bandgap occurs in the MPPC,where electromagnetic wave transmission is forbidden.The relationship between the central wavelength of the bandgap,the refractive index,and the dimension parameters is given by:

    Figure 3.Real and imaginary parts of the plasma relative permittivity for different gas pressures and electron densities.The frequency of the incident electromagnetic wave is 0.5 THz.

    where λ0is the central wavelength; m is an integer,which represents the formation of the periodic bandgap structure in the MPPC.Here,nqand npdenote the refractive indexes of quartz and plasma,respectively,and dqand dpare the widths of the quartz and the plasma,respectively.These refractive indexes are given by

    By substituting the corresponding central frequency of the bandgapf0=c/λ0,an expression related to the central frequency and the relative dielectric permittivity is obtained:

    Figure 4.(a)Map showing the loss tangent as a function of incident wave frequency(ω)and electron density in collisional plasma.(b)The relationship between the loss tangent and the electron density at an incident wave frequency of 0.2 THz,corresponding to the dashed line in figure 4(a).The collision frequency is 6.07×1012 rad s?1.The grey strip(region 1)indicates that the value of Re(εr)approaches to zero.

    Therefore,by tuning the microplasma permittivity and the dimensional parameters,the bandgap characteristics of the MPPC can be changed accordingly.

    2.4.Transfer matrix method

    The TMM is used to calculate the dispersion relation and the transmission of a vertically incident THz wave through the PPC,as shown in figure 1.This well-known method is based on the analytic solution of the 2D Maxwell’s wave equation,described in detail in[27].Therefore,only a brief introduction to this method is presented here.The electric field and the magnetic field at both sides of the Nth layer can be linked by the transfer matrix MNas follows:

    where N equates to 4 in figure 1,and the transfer matrix MNis given by

    In this way,a single matrix can be defined for each dielectric layer;the relationship between the incident THz wave and the output THz wave can be therefore be calculated by

    where M is the final transfer matrix for the one-dimensional collisional PPC with N dielectric layers.Here,A,B,C,and D are dependent on the permittivity of the dielectric layer and the dimensional parameters of the PPC.The power transmission coefficient of the THz wave and the dispersion relation are obtained from the final transfer matrix,given by

    where η0and ηN+1denote the impedance of the THz wave at the left side and the right side of the PPC in figure 1,respectively,k denotes the wave number and δq,δp,ηq,and ηpare derived from

    where c and μ0are the speed of light and the vacuum permeability,respectively.

    3.Numerical results

    3.1.Effects of collision frequency

    To investigate the effects of the collision frequency,the THz transmission coefficient and the dispersion relation,both with and without consideration of the collision frequency,are calculated by the TMM,as shown in figure 5.Three THz stopbands at central frequencies of 0.85 THz,2.29 THz,and 3.75 THz are observed in both MPPCs,labeled as 1st BG,2nd BG,and 3rd BG.The bandwidths of these stopbands,defined by the full width at half maximum of the transmission coefficient,decrease as the central frequency increases.Moreover,a significant attenuation of the transmission coefficient is observed in the case of collisional microplasma.The lower the incident frequency,the larger the attenuation of the transmission coefficient will be.Similarly,when the collision frequency of 6.07×1012rad s?1(corresponding to the atmospheric pressure)is taken into account in the plasma relative permittivity,the dispersion relation is also significantly changed,indicating that the band structure of the PPC should be affected.

    Given that the collision frequency is proportional to the gas pressure,we now investigate the effects of collision frequency on the THz transmission coefficient and the central frequency of the 1st BG via variations in the gas pressure.In figure 6(a),when the gas pressure increases,the transmission coefficient is shown to increase inside the 1st BG and on the left side of the 1st BG,but decreases at the right side of the 1st BG.On the other hand,the transmission coefficient around the 2nd BG and the 3rd BG decreases monotonically with the increase in gas pressure,as shown in figures 6(b)and(c).In addition,as the gas pressure increases from 50.5 kPa and 202 kPa,the central frequency shifts of the 1st BG,the 2nd BG,and the 3rd BG are 81 GHz,30.7 GHz,and 9.8 GHz,respectively.In figure 7,when the gas pressure increases from 50.5 kPa to 202 kPa,the central frequency of the 1st BG decreases monotonically from 0.871 THz to 0.79 THz.Therefore,the collision frequency has a great impact on the attenuation of the transmission coefficient and the band structure of the 1st BG for atmospheric pressure microplasmas.It should be noted that due to the strong attenuation of the THz wave propagation,the band structure of the 1st BG is severely distorted,resulting in the bandwidths for the 1st BG at different gas pressure being barely identifiable.

    Figure 5.The THz transmission coefficient,and the dispersion relation of the MPPC for both collision-free and collisional microplasma.(a),(b)Collision-free microplasma,with an electron density of 1016 cm-3.(c),(d)Collisional microplasma,νm=6.07×1012 rad s?1,with an electron density of 1016 cm-3.

    Figure 6.Transmission coefficients for(a)1st BG,(b)2nd BG,and(c)3rd BG as a function of the incident wave frequency at different gas pressures.The electron density is 1016 cm-3.

    Figure 7.The central frequency of the 1st BG as a function of gas pressure.The electron density is 1016 cm-3.

    Figure 8.Transmission coefficient of the 1st BG as a function of incident wave frequency,with different electron densities.The collision frequency is 6.07×1012 rad s?1.An electron density of 0 represents vacuum.

    3.2.Effects of electron density

    Electron density is one of the key parameters in the study of plasma,and can be used to control the plasma’s relative permittivity by means of plasma frequency,as shown in equation(2).Figure 8 shows the transmission coefficient in the 1st BG as a function of electron density.When the electron density increases from 0 to 1016cm-3,the shift of the central frequency of the 1st BG is 110 GHz,and the transmission coefficient around the 1st BG decreases by at least one order of magnitude over all bandwidths.Figure 9 shows the relationship between the electron density,the central frequency,and the bandwidth of the 1st BG.When the electron density increases from 1013cm-3to 1015cm-3,both the central frequency and the bandwidth of the 1st BG remain almost unchanged at around 0.74 THz and 0.42 THz,respectively.However,as the electron density further increases from 1015cm-3to 1016cm-3,both the central frequency and the bandwidth increase significantly,up to 0.85 THz and 0.62 THz,respectively.In figure 10,as the electron density increases from 0 to 1016cm?3,the shifts of the central frequencies of the 2nd BG and the 3rd BG are 67 GHz and 53 GHz,respectively.The transmission coefficients at the central frequencies in the 2nd BG and the 3rd BG are attenuated by less than three times.In summation,the higher the central frequency of the THz stopband,the less significant the impact caused by the electron density will be.In addition,it should be pointed out that electron density has a minor impact on the THz band structure of the MPPC if the electron density is less than 1015cm-3.

    Figure 9.Central frequency and bandwidth of the 1st BG as a function of electron density.Here,the collision frequency is 6.07×1012 rad s?1.

    3.3.Effects of plasma width

    As illustrated by the results above,the electron density required for the THz wave modulation is greater than 1015cm?3.It is not easy to develop such a low-temperature,high electron-density plasma at high pressure.The recent development of microplasmas generated inside a capillary or a cavity is anticipated to meet these strict requirements[28-30].However,the dimensional parameters of microplasmas are strictly limited,and have a significant influence on the electron density of the microplasmas[28,31].It is therefore essential to investigate the effects of plasma width on the THz transmission coefficients of MPPCs.In figure 11,when the plasma width increases from 50 μm to 70 μm(i.e.larger lattice constant),the central frequencies for the 1st,2nd,and 3rd BGs are shifted toward a low frequency,and the transmission coefficients duly decrease.

    Figure 10.Transmission coefficients for(a)2nd BG and(b)3rd BG as a function of incident wave frequency,with different electron densities.The collision frequency for both is 6.07×1012 rad s?1.

    4.Discussion

    By introducing microplasma into the photonic crystal,a widely tunable MPPC is obtained.For example,as the electron density of the microplasma increases from 1015cm-3to 1016cm-3,the central frequency of the 1st BG shifts toward high frequency by 110 GHz,and the bandwidth of the 1st BG increases by about 200 GHz.The tunable range of the band structure caused by microplasma is much greater than that by controlling the temperature or the electric field to change the conductivity of the metamaterial[10,11].This method of THz wave modulation is probably useful in chemical spectroscopy[1],biomedical diagnostics[2],and astronomical observations[4].

    Figure 11.The THz transmission coefficients of the MPPCs with different plasma widths.(a)50 μm,(b)60 μm,and(c)70 μm.The width of the quartz remains constant.Here,the collision frequency is 6.07×1012 rad s?1,and the electron density is 1016 cm-3.

    4.1.The mechanism of THz wave modulation

    With regard to the modulation of the central frequency of the THz stopband,based on equation(6),the central frequency of the stopbands in the MPPC is related to the relative permittivity of the microplasma.As shown in equations(1)and(2),the relative permittivity of the microplasma is dependent on electron density and collision frequency,which are proportional to the gas pressure.The relative permittivity of the microplasma is a complex,in contrast to the permittivity of quartz.The imaginary component of plasma permittivity represents the energy dissipation of the electromagnetic wave,and the real part of the plasma permittivity is related to the band structure of the MPPC[26].In order to clarify the relationship between gas pressure,electron density,and the central frequency of the THz stopband,the imaginary component of the plasma permittivity is neglected.By substituting equations(1)and(2)into equation(6),the central frequency can be given by

    Figure 12.Theoretical central frequency of the 1st BG as a function of(a)gas pressure,and(b)electron density.

    With the knowledge of the values of εq,dq,dp,m,e,me,ε0,and c,the relationship between the central frequency of the 1st BG,its electron density,and the gas pressure can be calculated theoretically,as shown in figure 12.As the gas pressure increases from 50.5 kPa to 202 kPa,the central frequency of the first THz stopband decreases from 0.95 THz to 0.78 THz,which is slightly higher than that for TMM in figure 7.When the electron density increases from 1013cm-3to 1016cm-3,the central frequency of the first THz stopband increases from 0.74 THz to 0.86 THz,which is in good agreement with the results calculated via TMM in figure 9(a).

    With regard to the modulation of the bandwidth of the THz stopband,the bandwidth is usually dependent on the ratio of the refractive index of these two dielectrics,and their dimensional parameters.In photonic crystals,the higher the ratio of the refractive index(nq/np),the larger the bandwidth will be[32].The refractive index of the plasma is determined by the plasma permittivity,expressed as[26]:

    Similarly,with a fixed collision frequency of 6.07×1012rad s?1,the relationship between the refractive index of the plasma and the electron density in relation to different incident wave frequencies is calculated theoretically.Figure 13 indicates the way in which,as the electron density increases from 1013cm-3to 1015cm?3,the refractive index of the plasma decreases slowly from 1.0 to 0.97 for a wave frequency of 0.8 THz,resulting in an almost constant ratio of the refractive index of these two dielectrics.However,when the electron density further increases up to 1016cm?3,the refractive index of the plasma decreases significantly,from 0.97 to 0.8,leading to a rapid increase in nq/npvalues,and therefore an increase in the bandwidth of the stopband.These results are consistent with the numerical results calculated by TMM in figure 9(b).Therefore,by tuning the electron density and the gas pressure,the permittivity of the microplasma can be varied,resulting in the modulation of the structure of the THz stopbands in MPPCs.

    Figure 13.Theoretical refractive index of plasma as a function of electron density,with different incident wave frequencies.

    4.2.The mechanism of THz wave attenuation

    From the perspective of electromagnetic wave propagation,collisional plasma is a dissipative medium,which leads to the attenuation of THz wave propagation.This is because the electrons gain energy from the electric field of the THz wave,and then transfer it to the molecules or atoms through frequent elastic collisions[33].The absorption index is given by[26]:

    The absorption index is therefore shown to be greatly dependent on the third term,,corresponding to the imaginary part of the plasma permittivity.For electron densities of 1015cm-3and 1016cm-3,the plasma frequencies are 0.28 THz and 0.9 THz,respectively.The collision frequency at atmospheric pressure is 0.96 THz.On one hand,when ω?νmand ω?ωpe,equation(15)can be simplified as follows:

    Therefore,the higher the electron density and the collision frequency are,the larger the absorption index is,and the stronger the attenuation of THz wave propagation will be.This exact phenomenon can be observed for both the 2nd and 3rd BGs in figures 6(b)and(c),and figure 10.On the other hand,the frequency of the incident THz wave in the 1st BG is close to the plasma frequency,but less than the collision frequency at high pressure.By assuming ω?νm,the equation(15)can be simplified as follows:

    It is therefore clear that an increase in collision frequency leads to a decrease in the absorption index,resulting in a weak attenuation of THz wave propagation.This is why a higher transmission coefficient in the 1st BG is obtained as the gas pressure increases.From the perspective view of plasma physics,the higher the gas pressure is,the shorter the mean free path will be.For ω?νm,several elastic collisions between electrons and atoms or molecules occur during one period of THz wave.When these collisions become more frequent,it becomes increasingly difficult for the electrons to gain energy from the electric field,due to the shorter mean free path.

    For microplasma generated at atmospheric pressure,the collision frequency for momentum transfer between electrons and atoms or molecules is in the THz range,which cannot be neglected in terms of the plasma’s relative permittivity.In this work,taking the collision frequency into account,the characteristics of the THz band structure of the MPPC are investigated for the first time.In figure 5,taking into account the collision frequency at 101 kPa,the THz transmission coefficient and the dispersion relation of the MPPC make a difference.Further investigation shows that collision frequency and electron density have a great impact on both the transmission coefficient and the band structure of the first THz stopband for atmospheric pressure microplasmas.This differs from the results in previous research on the subject of microwave attenuation due to low-pressure volume plasma[34,35],where the collision frequency has a negligible effect on the band structure of the PPC in the microwave range,but a great impact on the attenuation of microwave propagation in plasma.In low-pressure volume plasma[34,35],since the collision frequency is negligible as compared to the incident wave frequency,the real part of the plasma relative permittivity is almost independent of the collision frequency,and the imaginary part of the plasma relative permittivity is proportional to the collision frequency,based on equation(1).

    4.3.Microplasma for the modulation of THz wave

    With regard to the characteristics of plasma employed for the modulation of a THz wave,it is interesting to note that in order to have a significant impact on the band structure of the MPPC,the electron density of the plasma must be at least 1015cm?3.It is almost impossible to realize such low-temperature plasma with high electron density at low pressure.However,several reported methods,based on microdischarge,nanosecond pulsed discharge,RF discharge,or microwave discharge represent several possible pathways to generating high-ionized plasma at atmospheric pressure.Our recent research has demonstrated atmospheric pressure microplasma sources generated inside capillaries,with an electron density of 1015-1017cm-3[28,31,36].These results suggest the possibility of achieving wide tunability of the THz band structure in MPPCs.Moreover,in relation to plasma width,simulation results show that the stopband structure of the plasma photonic crystal exists for plasma widths of tens of micrometers,disappearing as the plasma width further increases up to 800 micrometers.This indicates that the characteristic size of the plasma should ideally be tens to hundreds of micrometers in order to achieve successful modulation of a THz wave.In other words,microplasma at atmospheric pressure is well-suited to the modulation of THz waves.

    In addition,although the high gas pressure and high electron density of the MPPC can significantly alter the bandwidth and the central frequency of THz stopbands,there is a concomitant sizeable reduction in the transmission coefficient of the THz wave,as shown in figures 6(a)and 8.For example,the transmission coefficient decreases by at least one order of magnitude for an electron density increase of up to 1016cm?3,which may be not appropriate for the modulation of THz waves in real applications.Therefore,taking full account of the attenuation and modulation of THz wave propagation,the ideal electron density of microplasma at atmospheric pressure is suggested to be in the range from 1×1015-6×1015cm-3,where the attenuation of the THz wave propagation by the MPPC is then less than 10 dB.

    5.Conclusion

    Taking collision frequency into consideration,the THz band structure and the transmission characteristics of a onedimensional MPPC are simulated via TMM.At high pressure,both the electron density and the collision frequency for momentum transfer between electrons and atoms or molecules have a significant impact on the attenuation of the THz wave and the THz band structure,by altering the plasma permittivity.On increasing the electron density from 1015cm-3to 1016cm-3,the central frequency of the 1st BG shifts toward high frequency by 110 GHz,and the bandwidth of the 1st BG increases by about 200 GHz,far higher than the results achieved by previously recorded methods[10-13].Moreover,as the gas pressure increases,the transmission coefficient increases at low frequencies in the incident THz wave,but decreases at higher frequencies of incident THz wave,where ω?νmand ω?ωpe.On the other hand,the transmission coefficient decreases monotonically with an increase in the electron density.In addition,to achieve wide tunability of the THz stopbands of MPPCs,the electron density must be as high as 1015cm?3,and the ideal feature size of the plasma should be submillimeter.This study indicates that high electron-density microplasma at atmospheric pressure is likely to be suitable for the modulation of THz wave propagation in MPPCs.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China(No.51977110).

    猜你喜歡
    張潮
    張潮與黃周星交游略論
    Plasma propagation in single-particle packed dielectric barrier discharges: joint effects of particle shape and discharge gap
    幽夢影
    消失的趙三
    Characteristics of the propagation of partial discharge ultrasonic signals on a transformer wall based on Sagnac interference
    雅致的生活
    兇手叫門露馬腳
    張潮和池米(短篇小說)
    雅致的生活
    夜來臨
    文學(xué)港(2014年7期)2015-01-04 22:51:27
    别揉我奶头 嗯啊视频| 精品久久久噜噜| 青春草国产在线视频| 大香蕉97超碰在线| 日日摸夜夜添夜夜爱| 国产精品日韩av在线免费观看| 午夜激情久久久久久久| 亚洲av.av天堂| 亚洲国产欧美在线一区| 亚洲av国产av综合av卡| 成人毛片a级毛片在线播放| av网站免费在线观看视频 | 三级经典国产精品| 网址你懂的国产日韩在线| 不卡视频在线观看欧美| 亚洲精华国产精华液的使用体验| 国产又色又爽无遮挡免| 久久亚洲国产成人精品v| 国产亚洲精品av在线| 赤兔流量卡办理| 人人妻人人看人人澡| 18+在线观看网站| 在线播放无遮挡| 亚洲无线观看免费| 欧美人与善性xxx| 美女被艹到高潮喷水动态| 秋霞伦理黄片| 赤兔流量卡办理| 18+在线观看网站| 99re6热这里在线精品视频| 免费看光身美女| 九草在线视频观看| 三级男女做爰猛烈吃奶摸视频| 免费看光身美女| 国产精品一区二区三区四区免费观看| 亚洲在久久综合| 一个人看的www免费观看视频| 久久精品国产亚洲av涩爱| 成人高潮视频无遮挡免费网站| 伊人久久精品亚洲午夜| 伦理电影大哥的女人| 国产精品一区二区在线观看99 | 国产精品无大码| 26uuu在线亚洲综合色| 3wmmmm亚洲av在线观看| 日本色播在线视频| 亚洲伊人久久精品综合| 国产成人91sexporn| 久久久久久久久久黄片| 一级片'在线观看视频| 我要看日韩黄色一级片| 久久97久久精品| 午夜爱爱视频在线播放| 日日摸夜夜添夜夜爱| 日韩一区二区三区影片| 久久精品国产鲁丝片午夜精品| 中文字幕免费在线视频6| 国产亚洲精品久久久com| 天天躁日日操中文字幕| 中文在线观看免费www的网站| 水蜜桃什么品种好| 51国产日韩欧美| 亚洲精品第二区| 国产亚洲一区二区精品| 国产69精品久久久久777片| 国产精品国产三级国产专区5o| a级毛色黄片| 国产精品蜜桃在线观看| 国产91av在线免费观看| 色尼玛亚洲综合影院| av黄色大香蕉| 十八禁国产超污无遮挡网站| a级一级毛片免费在线观看| 天堂av国产一区二区熟女人妻| 欧美97在线视频| 赤兔流量卡办理| 伦精品一区二区三区| 人妻制服诱惑在线中文字幕| 国产探花极品一区二区| 久久久久国产网址| 少妇高潮的动态图| 国产成人精品福利久久| 国产精品1区2区在线观看.| 高清午夜精品一区二区三区| 国产免费福利视频在线观看| 2022亚洲国产成人精品| 国产精品无大码| 国产伦精品一区二区三区视频9| 高清日韩中文字幕在线| 久久人人爽人人爽人人片va| 最近视频中文字幕2019在线8| 日本爱情动作片www.在线观看| 日本爱情动作片www.在线观看| 嫩草影院入口| 一本久久精品| 亚洲人与动物交配视频| 一级av片app| 天堂影院成人在线观看| 国产成人午夜福利电影在线观看| 亚洲欧美成人精品一区二区| 国产探花极品一区二区| 插逼视频在线观看| 亚洲精品自拍成人| 看非洲黑人一级黄片| 亚州av有码| 噜噜噜噜噜久久久久久91| 一级毛片久久久久久久久女| 国产在视频线精品| 国产中年淑女户外野战色| 欧美极品一区二区三区四区| 免费看av在线观看网站| 午夜免费激情av| 久久久久精品性色| 嘟嘟电影网在线观看| 亚洲av成人精品一区久久| 久久亚洲国产成人精品v| 欧美不卡视频在线免费观看| 国产有黄有色有爽视频| 蜜臀久久99精品久久宅男| 日本wwww免费看| 国产亚洲午夜精品一区二区久久 | 色哟哟·www| 最近中文字幕高清免费大全6| 美女内射精品一级片tv| 成人性生交大片免费视频hd| 晚上一个人看的免费电影| 国产日韩欧美在线精品| 国产成人午夜福利电影在线观看| 波野结衣二区三区在线| 国产免费一级a男人的天堂| 2021天堂中文幕一二区在线观| 日韩中字成人| 一级片'在线观看视频| 久久6这里有精品| 国产黄色视频一区二区在线观看| 精品久久久久久成人av| 国产探花在线观看一区二区| 国产伦一二天堂av在线观看| av福利片在线观看| 免费高清在线观看视频在线观看| 99热网站在线观看| 18禁在线播放成人免费| 日本欧美国产在线视频| 亚洲av不卡在线观看| 天堂网av新在线| 十八禁网站网址无遮挡 | 91久久精品电影网| 国产真实伦视频高清在线观看| 七月丁香在线播放| 国产高清国产精品国产三级 | 永久网站在线| 97人妻精品一区二区三区麻豆| 欧美精品国产亚洲| 女人十人毛片免费观看3o分钟| 亚洲av免费高清在线观看| 成人亚洲精品av一区二区| 欧美一级a爱片免费观看看| 欧美人与善性xxx| 精品一区二区三区人妻视频| 亚洲熟女精品中文字幕| 精品午夜福利在线看| 偷拍熟女少妇极品色| 好男人在线观看高清免费视频| 中文字幕亚洲精品专区| 久久鲁丝午夜福利片| 寂寞人妻少妇视频99o| 久久久久精品久久久久真实原创| 久久久国产一区二区| 日本av手机在线免费观看| 在线免费十八禁| 狠狠精品人妻久久久久久综合| 日韩欧美一区视频在线观看 | 国产有黄有色有爽视频| 国产精品无大码| 一边亲一边摸免费视频| 97热精品久久久久久| 久久草成人影院| 一级毛片aaaaaa免费看小| h日本视频在线播放| 纵有疾风起免费观看全集完整版 | 欧美zozozo另类| 久久精品夜色国产| 国产综合懂色| 内射极品少妇av片p| 日韩av不卡免费在线播放| 热99在线观看视频| 国产视频首页在线观看| 欧美成人一区二区免费高清观看| 亚洲精品第二区| 中文字幕人妻熟人妻熟丝袜美| 国产探花极品一区二区| 中文天堂在线官网| 高清日韩中文字幕在线| 一级毛片久久久久久久久女| 国产淫片久久久久久久久| av线在线观看网站| 尾随美女入室| 国产一级毛片七仙女欲春2| 久久精品国产自在天天线| 国产亚洲最大av| 最后的刺客免费高清国语| 欧美成人精品欧美一级黄| 国产精品人妻久久久久久| 永久网站在线| 美女国产视频在线观看| 一级毛片久久久久久久久女| 精品一区二区三区视频在线| 亚洲精品一区蜜桃| 大香蕉久久网| 国产三级在线视频| 高清av免费在线| 国产精品久久视频播放| 人妻少妇偷人精品九色| 国产在视频线在精品| 伦理电影大哥的女人| 深夜a级毛片| 国产高清国产精品国产三级 | 国产色爽女视频免费观看| 日韩一区二区视频免费看| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 乱人视频在线观看| 亚洲av免费高清在线观看| videos熟女内射| 国产精品.久久久| 3wmmmm亚洲av在线观看| 国产精品精品国产色婷婷| 狂野欧美白嫩少妇大欣赏| 日韩电影二区| 亚洲电影在线观看av| 国产精品久久久久久精品电影| 一夜夜www| 晚上一个人看的免费电影| 熟妇人妻久久中文字幕3abv| 亚洲精品日韩在线中文字幕| 97在线视频观看| 久久久色成人| 国产综合懂色| 夫妻性生交免费视频一级片| 热99在线观看视频| 欧美zozozo另类| 久久久久久久国产电影| 日韩在线高清观看一区二区三区| 亚洲精品自拍成人| 色哟哟·www| 成人av在线播放网站| 联通29元200g的流量卡| 寂寞人妻少妇视频99o| 我的女老师完整版在线观看| 99久国产av精品国产电影| 久久精品熟女亚洲av麻豆精品 | 国产亚洲午夜精品一区二区久久 | 青春草国产在线视频| 人人妻人人看人人澡| 国产精品一区二区性色av| 国产一区亚洲一区在线观看| 午夜福利视频1000在线观看| 在线观看人妻少妇| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品专区欧美| 能在线免费看毛片的网站| 丝瓜视频免费看黄片| 伊人久久精品亚洲午夜| 日本午夜av视频| 视频中文字幕在线观看| 91精品国产九色| 欧美不卡视频在线免费观看| 美女主播在线视频| 舔av片在线| 人妻系列 视频| 欧美人与善性xxx| 91aial.com中文字幕在线观看| 青春草亚洲视频在线观看| 国产单亲对白刺激| 久久久欧美国产精品| 午夜日本视频在线| 亚洲精品国产成人久久av| 麻豆精品久久久久久蜜桃| 日本熟妇午夜| 少妇被粗大猛烈的视频| 天美传媒精品一区二区| 精华霜和精华液先用哪个| 夫妻午夜视频| 日韩亚洲欧美综合| 亚洲人成网站高清观看| 成人毛片60女人毛片免费| xxx大片免费视频| 午夜爱爱视频在线播放| 18+在线观看网站| av免费观看日本| 日日干狠狠操夜夜爽| 午夜精品国产一区二区电影 | 亚洲av电影不卡..在线观看| 女人久久www免费人成看片| 日韩制服骚丝袜av| 亚洲婷婷狠狠爱综合网| 亚洲欧美中文字幕日韩二区| 国产乱人偷精品视频| 亚洲无线观看免费| 国产三级在线视频| 99久久精品热视频| 69人妻影院| 亚洲欧美中文字幕日韩二区| 国产在线一区二区三区精| 亚洲经典国产精华液单| 搡女人真爽免费视频火全软件| 最后的刺客免费高清国语| 十八禁国产超污无遮挡网站| 你懂的网址亚洲精品在线观看| 我的女老师完整版在线观看| 亚洲av免费高清在线观看| 大又大粗又爽又黄少妇毛片口| av在线蜜桃| 国国产精品蜜臀av免费| 人妻系列 视频| 国产成人freesex在线| 搡女人真爽免费视频火全软件| av国产久精品久网站免费入址| 中文在线观看免费www的网站| 精品一区二区三区人妻视频| 丰满乱子伦码专区| 亚洲欧美中文字幕日韩二区| 波多野结衣巨乳人妻| av在线观看视频网站免费| 中国国产av一级| a级毛色黄片| 内射极品少妇av片p| 欧美日韩精品成人综合77777| 欧美日韩国产mv在线观看视频 | 国产精品日韩av在线免费观看| 国产精品国产三级国产专区5o| 伊人久久精品亚洲午夜| 国产一区二区在线观看日韩| 69av精品久久久久久| 老师上课跳d突然被开到最大视频| 国产精品一区二区性色av| 狂野欧美激情性xxxx在线观看| 国产精品伦人一区二区| 国产精品一区二区在线观看99 | 国产成人a∨麻豆精品| 日本一本二区三区精品| freevideosex欧美| 久久6这里有精品| 别揉我奶头 嗯啊视频| 精品国产露脸久久av麻豆 | 日本免费a在线| 中文字幕免费在线视频6| 精品国产露脸久久av麻豆 | 国产精品美女特级片免费视频播放器| 亚洲国产日韩欧美精品在线观看| 成人国产麻豆网| 久久精品国产亚洲网站| 成人亚洲欧美一区二区av| 我要看日韩黄色一级片| 中文字幕人妻熟人妻熟丝袜美| 亚洲真实伦在线观看| 久久午夜福利片| 欧美潮喷喷水| 亚洲欧美精品专区久久| 一级毛片我不卡| 国产伦精品一区二区三区四那| 少妇裸体淫交视频免费看高清| 久久久久久久久久久免费av| 日本爱情动作片www.在线观看| 九九久久精品国产亚洲av麻豆| 日韩欧美精品免费久久| 久久久久久久亚洲中文字幕| 午夜日本视频在线| 卡戴珊不雅视频在线播放| 99热这里只有是精品50| 亚洲av中文av极速乱| 国产高清有码在线观看视频| 国产综合懂色| 肉色欧美久久久久久久蜜桃 | 午夜日本视频在线| 丝袜美腿在线中文| 一级毛片aaaaaa免费看小| 国产一区二区三区av在线| 性插视频无遮挡在线免费观看| 日韩欧美国产在线观看| 午夜福利网站1000一区二区三区| 美女xxoo啪啪120秒动态图| 又爽又黄a免费视频| 国内少妇人妻偷人精品xxx网站| 欧美日韩视频高清一区二区三区二| 伊人久久国产一区二区| 亚洲精品日韩在线中文字幕| 久久99热这里只频精品6学生| 国产淫语在线视频| 亚洲欧美日韩无卡精品| 国产亚洲av嫩草精品影院| 人人妻人人看人人澡| 九九久久精品国产亚洲av麻豆| 三级毛片av免费| 亚洲成人av在线免费| 听说在线观看完整版免费高清| 亚洲精品一区蜜桃| 天天躁夜夜躁狠狠久久av| 国产成人a∨麻豆精品| 国产淫语在线视频| 又大又黄又爽视频免费| 熟女人妻精品中文字幕| 99久久九九国产精品国产免费| 精品久久久久久久久av| 亚洲不卡免费看| 一区二区三区高清视频在线| 嫩草影院新地址| 亚洲人与动物交配视频| 亚洲,欧美,日韩| 国产国拍精品亚洲av在线观看| 18禁动态无遮挡网站| 国产精品一区二区性色av| 最近中文字幕高清免费大全6| 国产激情偷乱视频一区二区| 日韩精品有码人妻一区| 99热全是精品| 免费观看av网站的网址| 一夜夜www| 秋霞伦理黄片| 99热网站在线观看| 日韩一区二区三区影片| 熟女电影av网| 亚洲av一区综合| 精品久久久噜噜| 又大又黄又爽视频免费| 亚洲人与动物交配视频| 熟女电影av网| 国产精品久久视频播放| 国内精品一区二区在线观看| 一边亲一边摸免费视频| 99热这里只有精品一区| 久久久久久九九精品二区国产| 日韩av在线大香蕉| 亚洲精华国产精华液的使用体验| www.色视频.com| 韩国高清视频一区二区三区| 亚洲欧美一区二区三区黑人 | 99久久精品热视频| a级毛色黄片| 国产精品1区2区在线观看.| 国产亚洲最大av| 免费黄频网站在线观看国产| 免费观看a级毛片全部| 日韩电影二区| 久久99热这里只有精品18| 欧美丝袜亚洲另类| 国产成人精品福利久久| 大又大粗又爽又黄少妇毛片口| 麻豆成人av视频| 黑人高潮一二区| 极品教师在线视频| 麻豆久久精品国产亚洲av| 嫩草影院新地址| 夜夜爽夜夜爽视频| 久久精品国产自在天天线| 狂野欧美白嫩少妇大欣赏| 国产高清不卡午夜福利| 国产精品女同一区二区软件| 精品人妻一区二区三区麻豆| 免费观看av网站的网址| 一级片'在线观看视频| 国产91av在线免费观看| 日韩欧美 国产精品| av在线播放精品| 熟妇人妻不卡中文字幕| 大陆偷拍与自拍| 久久久久久国产a免费观看| 亚洲最大成人手机在线| videos熟女内射| 亚洲精品日韩av片在线观看| 男人爽女人下面视频在线观看| 99久久九九国产精品国产免费| 精品久久国产蜜桃| 别揉我奶头 嗯啊视频| 精品人妻偷拍中文字幕| 国产成人a∨麻豆精品| 国产成人一区二区在线| 日韩大片免费观看网站| 麻豆精品久久久久久蜜桃| 国产av码专区亚洲av| 一级毛片久久久久久久久女| 国产大屁股一区二区在线视频| 2021天堂中文幕一二区在线观| 国产精品嫩草影院av在线观看| 国产毛片a区久久久久| 久久久精品94久久精品| 观看免费一级毛片| 久久久久精品久久久久真实原创| 久久久久网色| 人人妻人人澡人人爽人人夜夜 | 国产成人aa在线观看| 肉色欧美久久久久久久蜜桃 | 毛片女人毛片| 网址你懂的国产日韩在线| 欧美日韩视频高清一区二区三区二| 小蜜桃在线观看免费完整版高清| 欧美精品国产亚洲| 国产色爽女视频免费观看| 少妇丰满av| 成人二区视频| 2022亚洲国产成人精品| 2021少妇久久久久久久久久久| 久久久久国产网址| 久久久精品欧美日韩精品| 最近最新中文字幕大全电影3| av播播在线观看一区| 国产精品久久久久久精品电影小说 | 国产大屁股一区二区在线视频| 五月伊人婷婷丁香| 亚洲精品成人av观看孕妇| 永久网站在线| 午夜福利在线观看吧| 十八禁网站网址无遮挡 | 99热网站在线观看| 中文字幕免费在线视频6| 免费高清在线观看视频在线观看| 99久国产av精品国产电影| 亚洲成色77777| 少妇猛男粗大的猛烈进出视频 | 免费观看精品视频网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色综合亚洲欧美另类图片| 男人狂女人下面高潮的视频| 免费高清在线观看视频在线观看| 能在线免费观看的黄片| 肉色欧美久久久久久久蜜桃 | 亚洲av一区综合| 午夜免费男女啪啪视频观看| 国产精品一区二区三区四区久久| 午夜久久久久精精品| 成人国产麻豆网| 18禁动态无遮挡网站| 久久99蜜桃精品久久| 久久精品久久久久久久性| 久久99精品国语久久久| 三级国产精品片| 亚洲精品亚洲一区二区| 国产精品国产三级国产专区5o| 日日啪夜夜撸| 国产精品一二三区在线看| 国产av在哪里看| 麻豆精品久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 麻豆成人午夜福利视频| 午夜福利高清视频| 少妇人妻精品综合一区二区| 午夜日本视频在线| 三级国产精品片| 街头女战士在线观看网站| 久久精品久久久久久久性| 美女大奶头视频| 又爽又黄无遮挡网站| 久久国产乱子免费精品| 晚上一个人看的免费电影| 精品国产三级普通话版| 亚洲av在线观看美女高潮| 久久99热这里只有精品18| 最近最新中文字幕大全电影3| 久久久久久久久久黄片| 成人综合一区亚洲| 亚洲av成人av| 久久久久久久久中文| 91久久精品国产一区二区三区| 久久久久免费精品人妻一区二区| 69人妻影院| 亚洲成色77777| 插阴视频在线观看视频| 欧美三级亚洲精品| 国产精品熟女久久久久浪| 国内少妇人妻偷人精品xxx网站| 欧美xxxx性猛交bbbb| 汤姆久久久久久久影院中文字幕 | 国产午夜精品一二区理论片| 乱人视频在线观看| 欧美激情国产日韩精品一区| 日本爱情动作片www.在线观看| 伦理电影大哥的女人| 亚洲高清免费不卡视频| 在线天堂最新版资源| 国产精品熟女久久久久浪| av网站免费在线观看视频 | 亚洲人成网站在线观看播放| 国产精品爽爽va在线观看网站| 中国美白少妇内射xxxbb| 欧美日韩视频高清一区二区三区二| 搞女人的毛片| 在线观看av片永久免费下载| 国产av国产精品国产| 99视频精品全部免费 在线| 老司机影院成人| 久久久久久九九精品二区国产| 欧美xxxx性猛交bbbb| 亚洲熟女精品中文字幕| 国产精品蜜桃在线观看| 日韩精品有码人妻一区| 老师上课跳d突然被开到最大视频| 97超视频在线观看视频| 成人欧美大片| 青春草亚洲视频在线观看| 亚洲无线观看免费| 91精品伊人久久大香线蕉| 人体艺术视频欧美日本| av女优亚洲男人天堂| 汤姆久久久久久久影院中文字幕 | 国产探花极品一区二区| 日本av手机在线免费观看| 三级经典国产精品| 婷婷色综合www| 国内少妇人妻偷人精品xxx网站| 久久99热6这里只有精品| 亚洲精品日韩av片在线观看| 色尼玛亚洲综合影院| 五月玫瑰六月丁香|