• <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
    国产高清不卡午夜福利| 国产视频一区二区在线看| 久久欧美精品欧美久久欧美| 亚洲成人久久性| 国产探花极品一区二区| 日韩精品有码人妻一区| 免费观看精品视频网站| 99久久成人亚洲精品观看| 亚洲精品国产av成人精品 | 别揉我奶头~嗯~啊~动态视频| 午夜老司机福利剧场| 久久韩国三级中文字幕| 婷婷六月久久综合丁香| 一级黄色大片毛片| 两性午夜刺激爽爽歪歪视频在线观看| 国产激情偷乱视频一区二区| 少妇高潮的动态图| 九九久久精品国产亚洲av麻豆| 亚洲七黄色美女视频| 国产日本99.免费观看| 精品久久久久久久久久免费视频| 亚洲经典国产精华液单| 菩萨蛮人人尽说江南好唐韦庄 | 少妇人妻一区二区三区视频| а√天堂www在线а√下载| 欧美色欧美亚洲另类二区| 免费人成在线观看视频色| 亚洲电影在线观看av| 精品国产三级普通话版| 岛国在线免费视频观看| 国产亚洲欧美98| 亚洲中文日韩欧美视频| 大型黄色视频在线免费观看| 有码 亚洲区| 午夜精品在线福利| 99精品在免费线老司机午夜| 91久久精品电影网| 亚洲av美国av| 亚洲国产色片| 亚洲人成网站在线观看播放| 国产av在哪里看| 麻豆国产97在线/欧美| 国内精品宾馆在线| 精品久久久久久久人妻蜜臀av| 国内少妇人妻偷人精品xxx网站| 波多野结衣高清作品| 99久久精品国产国产毛片| 国产在视频线在精品| 免费av不卡在线播放| 一区福利在线观看| 最后的刺客免费高清国语| 国内久久婷婷六月综合欲色啪| 黄色一级大片看看| 国产蜜桃级精品一区二区三区| 国产高清有码在线观看视频| 国产乱人偷精品视频| 黄色日韩在线| 精品午夜福利在线看| 国产精华一区二区三区| 国产综合懂色| 一区二区三区高清视频在线| 国产精品精品国产色婷婷| 少妇熟女欧美另类| 日韩,欧美,国产一区二区三区 | 男女啪啪激烈高潮av片| 欧美一区二区亚洲| 欧美国产日韩亚洲一区| 日本-黄色视频高清免费观看| 在线免费观看的www视频| 国产成人影院久久av| 级片在线观看| 全区人妻精品视频| 有码 亚洲区| 午夜激情欧美在线| 免费看光身美女| www日本黄色视频网| 日日干狠狠操夜夜爽| 人妻久久中文字幕网| 人妻制服诱惑在线中文字幕| 国产成人91sexporn| 欧美成人一区二区免费高清观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲专区国产一区二区| 国产亚洲91精品色在线| 色视频www国产| 欧美+日韩+精品| 美女大奶头视频| 91在线精品国自产拍蜜月| 欧美成人免费av一区二区三区| 国产精品国产高清国产av| 成人国产麻豆网| 亚洲高清免费不卡视频| 精品一区二区三区av网在线观看| 国产高清视频在线播放一区| 免费看美女性在线毛片视频| 国产久久久一区二区三区| 三级男女做爰猛烈吃奶摸视频| 国产蜜桃级精品一区二区三区| 性色avwww在线观看| 国产色爽女视频免费观看| 身体一侧抽搐| 欧美性感艳星| 久久久久精品国产欧美久久久| 亚洲中文字幕日韩| 九九热线精品视视频播放| aaaaa片日本免费| 国产精品亚洲美女久久久| 在线观看一区二区三区| 好男人在线观看高清免费视频| 精品国内亚洲2022精品成人| 神马国产精品三级电影在线观看| 18禁在线播放成人免费| 成年免费大片在线观看| 欧美最黄视频在线播放免费| 国产欧美日韩一区二区精品| 亚洲欧美日韩卡通动漫| 成年女人毛片免费观看观看9| 亚洲性久久影院| 99久久久亚洲精品蜜臀av| 久久精品国产亚洲av香蕉五月| 69av精品久久久久久| 桃色一区二区三区在线观看| 亚洲精华国产精华液的使用体验 | 男女啪啪激烈高潮av片| 国产精品一区二区三区四区久久| 亚洲精品久久国产高清桃花| 成熟少妇高潮喷水视频| 国产成人a区在线观看| 亚洲无线观看免费| 黄色配什么色好看| 亚洲av中文av极速乱| 日韩欧美一区二区三区在线观看| 日本欧美国产在线视频| 国产探花极品一区二区| 国产欧美日韩精品亚洲av| 日本撒尿小便嘘嘘汇集6| 人妻制服诱惑在线中文字幕| 床上黄色一级片| 久久久久性生活片| 人人妻人人澡人人爽人人夜夜 | 久久精品国产亚洲av天美| 久久久久九九精品影院| 精品久久久久久久久久免费视频| 亚洲欧美日韩无卡精品| 给我免费播放毛片高清在线观看| 久久精品国产亚洲av涩爱 | 午夜福利在线观看吧| 亚洲欧美日韩高清在线视频| 国产真实伦视频高清在线观看| 中国美女看黄片| 亚洲欧美日韩高清在线视频| 欧美日韩精品成人综合77777| 高清毛片免费观看视频网站| 麻豆一二三区av精品| 国产精品无大码| 在线观看午夜福利视频| 91在线观看av| 热99在线观看视频| 少妇人妻精品综合一区二区 | 精品人妻一区二区三区麻豆 | 97人妻精品一区二区三区麻豆| 国产成人a∨麻豆精品| 日本一本二区三区精品| 午夜福利在线在线| 俺也久久电影网| 国产av麻豆久久久久久久| 国产三级在线视频| 99久久中文字幕三级久久日本| 99热精品在线国产| 国产精品电影一区二区三区| 成人毛片a级毛片在线播放| 黑人高潮一二区| 亚洲国产精品成人久久小说 | 亚洲国产色片| 看片在线看免费视频| 亚洲精品亚洲一区二区| 国产精品精品国产色婷婷| 少妇熟女欧美另类| 午夜免费男女啪啪视频观看 | 综合色av麻豆| 欧美性感艳星| 欧美最新免费一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 日韩 亚洲 欧美在线| 久久久久久久久久久丰满| 欧美日韩国产亚洲二区| 人妻制服诱惑在线中文字幕| 亚洲激情五月婷婷啪啪| 国产高清激情床上av| 日韩强制内射视频| 久久久国产成人免费| 午夜久久久久精精品| 少妇高潮的动态图| 2021天堂中文幕一二区在线观| 99在线人妻在线中文字幕| 女人十人毛片免费观看3o分钟| 干丝袜人妻中文字幕| 亚洲国产精品国产精品| 亚洲国产日韩欧美精品在线观看| 日韩欧美 国产精品| 国产精品一区二区三区四区久久| 女人被狂操c到高潮| 成人特级黄色片久久久久久久| 久久久久久国产a免费观看| 中出人妻视频一区二区| 九九在线视频观看精品| 91麻豆精品激情在线观看国产| 毛片一级片免费看久久久久| 久久人人爽人人片av| 美女被艹到高潮喷水动态| 男女啪啪激烈高潮av片| 波多野结衣高清作品| 午夜爱爱视频在线播放| 亚洲色图av天堂| 女生性感内裤真人,穿戴方法视频| 非洲黑人性xxxx精品又粗又长| 日韩欧美在线乱码| 一区二区三区高清视频在线| 国产国拍精品亚洲av在线观看| 美女 人体艺术 gogo| 成人亚洲欧美一区二区av| 色综合站精品国产| 精品久久久噜噜| 亚洲美女黄片视频| 老司机影院成人| 免费看光身美女| 床上黄色一级片| 亚洲国产精品久久男人天堂| 自拍偷自拍亚洲精品老妇| 美女cb高潮喷水在线观看| 搞女人的毛片| 97在线视频观看| 国产久久久一区二区三区| 神马国产精品三级电影在线观看| 麻豆国产av国片精品| 人妻制服诱惑在线中文字幕| 男女边吃奶边做爰视频| 亚洲在线观看片| 97人妻精品一区二区三区麻豆| 麻豆成人午夜福利视频| 国产视频一区二区在线看| 欧美性感艳星| 青春草视频在线免费观看| 国产一区二区在线观看日韩| a级毛色黄片| 亚洲欧美精品自产自拍| 六月丁香七月| 91狼人影院| 深夜精品福利| 国产乱人偷精品视频| 2021天堂中文幕一二区在线观| 久久九九热精品免费| 黄色一级大片看看| 国产黄色小视频在线观看| 黄色欧美视频在线观看| 国产成人91sexporn| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 最后的刺客免费高清国语| 内射极品少妇av片p| 有码 亚洲区| 精品人妻一区二区三区麻豆 | 男女之事视频高清在线观看| 日韩成人伦理影院| 国产一区二区激情短视频| 日本三级黄在线观看| 国产精品久久久久久亚洲av鲁大| 给我免费播放毛片高清在线观看| 99视频精品全部免费 在线| 精品免费久久久久久久清纯| 18+在线观看网站| 黑人高潮一二区| 免费黄网站久久成人精品| 免费看a级黄色片| 国产淫片久久久久久久久| h日本视频在线播放| 狂野欧美激情性xxxx在线观看| 在线天堂最新版资源| 欧美一区二区亚洲| 在线国产一区二区在线| 长腿黑丝高跟| 精品久久久噜噜| 久久婷婷人人爽人人干人人爱| 亚洲最大成人手机在线| av专区在线播放| 色综合亚洲欧美另类图片| 久久久色成人| 在线播放无遮挡| 变态另类丝袜制服| 久久久久久久亚洲中文字幕| 在线免费观看的www视频| 欧美高清性xxxxhd video| 人人妻,人人澡人人爽秒播| 免费观看的影片在线观看| 网址你懂的国产日韩在线| av在线天堂中文字幕| 免费观看精品视频网站| 十八禁网站免费在线| 精品久久久久久久久久久久久| 少妇人妻一区二区三区视频| 女人被狂操c到高潮| 国产三级中文精品| 精品人妻熟女av久视频| 精品一区二区三区人妻视频| 亚洲精品国产av成人精品 | 亚洲欧美中文字幕日韩二区| 人妻夜夜爽99麻豆av| 99riav亚洲国产免费| 少妇熟女aⅴ在线视频| 国产精品爽爽va在线观看网站| 99热精品在线国产| 亚洲内射少妇av| 老师上课跳d突然被开到最大视频| 国产不卡一卡二| 秋霞在线观看毛片| 国产三级中文精品| 一级黄片播放器| 久久热精品热| 久99久视频精品免费| 亚洲国产高清在线一区二区三| 国产精品久久视频播放| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 精品少妇黑人巨大在线播放 | 十八禁网站免费在线| 国产又黄又爽又无遮挡在线| 免费av毛片视频| 亚洲av电影不卡..在线观看| 看非洲黑人一级黄片| 夜夜夜夜夜久久久久| 男人舔奶头视频| 国产欧美日韩一区二区精品| 99久久久亚洲精品蜜臀av| 亚洲四区av| 国产熟女欧美一区二区| 禁无遮挡网站| 麻豆久久精品国产亚洲av| 99久久精品国产国产毛片| 九九在线视频观看精品| 免费观看的影片在线观看| 九九在线视频观看精品| 99九九线精品视频在线观看视频| 日本一本二区三区精品| 国产高清激情床上av| 精品99又大又爽又粗少妇毛片| 亚洲自拍偷在线| 日本黄色片子视频| 日本与韩国留学比较| 韩国av在线不卡| 国产三级中文精品| 国产亚洲欧美98| 亚洲真实伦在线观看| 久久精品国产亚洲网站| 午夜福利18| 日本黄色片子视频| 亚洲精品成人久久久久久| 成年免费大片在线观看| 自拍偷自拍亚洲精品老妇| 欧美色欧美亚洲另类二区| 亚洲精品在线观看二区| 日韩亚洲欧美综合| 黄色一级大片看看| 丰满乱子伦码专区| 日韩欧美国产在线观看| 精品人妻视频免费看| 久久午夜亚洲精品久久| 亚洲性久久影院| av免费在线看不卡| 国产亚洲欧美98| 成年女人看的毛片在线观看| 日本a在线网址| 亚洲丝袜综合中文字幕| 男人舔女人下体高潮全视频| 国产美女午夜福利| 欧美不卡视频在线免费观看| 成年女人永久免费观看视频| 亚洲第一电影网av| 日韩一本色道免费dvd| 亚洲精品影视一区二区三区av| 国产欧美日韩精品一区二区| 亚洲美女视频黄频| 非洲黑人性xxxx精品又粗又长| 在线观看免费视频日本深夜| 午夜亚洲福利在线播放| 国产免费男女视频| 亚洲国产精品成人综合色| 国产毛片a区久久久久| 欧美激情在线99| 精品福利观看| 春色校园在线视频观看| 精品一区二区三区视频在线观看免费| 午夜精品在线福利| 亚洲激情五月婷婷啪啪| 欧美中文日本在线观看视频| 日韩中字成人| 午夜久久久久精精品| 99热只有精品国产| 国产精品野战在线观看| 欧美一区二区亚洲| 白带黄色成豆腐渣| 最后的刺客免费高清国语| 久久精品91蜜桃| 麻豆国产av国片精品| 22中文网久久字幕| 亚洲一区高清亚洲精品| 黄片wwwwww| 亚洲人成网站在线播| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美清纯卡通| 麻豆久久精品国产亚洲av| 熟女人妻精品中文字幕| 国产av在哪里看| 免费看日本二区| 成人国产麻豆网| 国产在视频线在精品| 日韩三级伦理在线观看| 亚洲人成网站在线观看播放| 日本爱情动作片www.在线观看 | 亚洲欧美中文字幕日韩二区| 亚洲真实伦在线观看| 久久草成人影院| 热99在线观看视频| 亚洲在线自拍视频| or卡值多少钱| 黄色一级大片看看| 国产精华一区二区三区| 欧美+亚洲+日韩+国产| 国产亚洲欧美98| 黄片wwwwww| 精品人妻一区二区三区麻豆 | 99久久成人亚洲精品观看| 国产一区亚洲一区在线观看| 久久久久久久久中文| 少妇裸体淫交视频免费看高清| 不卡视频在线观看欧美| 一级毛片我不卡| 1000部很黄的大片| 日本黄色片子视频| 九九久久精品国产亚洲av麻豆| 国产色爽女视频免费观看| 国产成人福利小说| 久久精品国产亚洲av涩爱 | 国产又黄又爽又无遮挡在线| 男女边吃奶边做爰视频| 不卡视频在线观看欧美| 久久久久九九精品影院| 亚洲精品亚洲一区二区| 亚洲av不卡在线观看| 亚洲av.av天堂| 淫妇啪啪啪对白视频| 麻豆国产97在线/欧美| 色5月婷婷丁香| 天美传媒精品一区二区| 99视频精品全部免费 在线| 麻豆成人午夜福利视频| 国产一区二区三区av在线 | 99视频精品全部免费 在线| 国产片特级美女逼逼视频| 日本五十路高清| 中文字幕免费在线视频6| 色综合亚洲欧美另类图片| 中文字幕久久专区| 俄罗斯特黄特色一大片| av天堂中文字幕网| а√天堂www在线а√下载| 在线免费十八禁| 久久精品91蜜桃| 午夜老司机福利剧场| 婷婷精品国产亚洲av| 免费人成视频x8x8入口观看| 日韩国内少妇激情av| 免费大片18禁| 在线播放无遮挡| 亚洲精品国产成人久久av| 在线国产一区二区在线| 久久久精品94久久精品| 黄色欧美视频在线观看| 久久久久久伊人网av| 国产色婷婷99| 欧美日韩乱码在线| 中文亚洲av片在线观看爽| 最近2019中文字幕mv第一页| 黄片wwwwww| 精品一区二区三区av网在线观看| 偷拍熟女少妇极品色| av福利片在线观看| 国产精品一区www在线观看| 国产男靠女视频免费网站| 亚洲精品乱码久久久v下载方式| 波多野结衣高清作品| av在线亚洲专区| 亚洲五月天丁香| 国产精品一区二区性色av| 国产69精品久久久久777片| 国产成人91sexporn| 亚洲在线自拍视频| 国产精品一区二区三区四区久久| 亚洲丝袜综合中文字幕| 日日干狠狠操夜夜爽| 国产极品精品免费视频能看的| 午夜激情欧美在线| 国产伦精品一区二区三区视频9| 99热这里只有是精品50| 黄色视频,在线免费观看| 欧美最新免费一区二区三区| 久久精品夜色国产| 亚洲美女视频黄频| 国产在线男女| 亚洲高清免费不卡视频| 亚洲精品亚洲一区二区| 直男gayav资源| 久久精品综合一区二区三区| 精品福利观看| 狂野欧美白嫩少妇大欣赏| 午夜a级毛片| 欧美潮喷喷水| 色播亚洲综合网| ponron亚洲| 亚洲欧美精品综合久久99| 亚洲国产欧洲综合997久久,| 天天一区二区日本电影三级| 亚洲欧美日韩高清专用| 五月伊人婷婷丁香| 精品国内亚洲2022精品成人| 成人午夜高清在线视频| 国产精品综合久久久久久久免费| 午夜福利成人在线免费观看| 免费高清视频大片| 久久久久国产网址| 国产人妻一区二区三区在| 一级a爱片免费观看的视频| 十八禁国产超污无遮挡网站| 观看美女的网站| 欧美日韩综合久久久久久| 久久精品国产亚洲网站| 国产私拍福利视频在线观看| 欧美日韩国产亚洲二区| 亚洲四区av| 日本色播在线视频| 国产探花在线观看一区二区| 丝袜美腿在线中文| 久久精品人妻少妇| 级片在线观看| 超碰av人人做人人爽久久| 日本-黄色视频高清免费观看| 性欧美人与动物交配| 亚洲18禁久久av| 欧美激情国产日韩精品一区| 老司机影院成人| 欧美bdsm另类| 波多野结衣巨乳人妻| 亚洲精品影视一区二区三区av| 国产三级中文精品| 一本精品99久久精品77| 综合色丁香网| 亚洲成人久久性| 国产精品久久电影中文字幕| 欧美性猛交黑人性爽| h日本视频在线播放| 一级黄片播放器| 国产精品,欧美在线| 搡老熟女国产l中国老女人| 国产亚洲av嫩草精品影院| 午夜福利18| 中文字幕久久专区| 毛片一级片免费看久久久久| 日日啪夜夜撸| 99热这里只有是精品50| av视频在线观看入口| 中文字幕av成人在线电影| 99久久成人亚洲精品观看| 久久综合国产亚洲精品| 麻豆久久精品国产亚洲av| 99热这里只有精品一区| 亚洲最大成人av| 日本黄色视频三级网站网址| 国产精品久久久久久久久免| 国产国拍精品亚洲av在线观看| 国产白丝娇喘喷水9色精品| 免费在线观看成人毛片| 色吧在线观看| 日韩,欧美,国产一区二区三区 | 精品久久久久久久末码| 2021天堂中文幕一二区在线观| 久久久精品大字幕| 夜夜爽天天搞| 日本五十路高清| 天堂网av新在线| 国产精品久久久久久av不卡| 最好的美女福利视频网| 亚洲三级黄色毛片| 美女内射精品一级片tv| 色在线成人网| 日韩一区二区视频免费看| 毛片一级片免费看久久久久| 免费看a级黄色片| 成人综合一区亚洲| 91久久精品国产一区二区成人| 日韩一本色道免费dvd| 国产aⅴ精品一区二区三区波| а√天堂www在线а√下载| 又粗又爽又猛毛片免费看| 丝袜喷水一区| 色播亚洲综合网| 日本与韩国留学比较| 国产精品电影一区二区三区| 日韩欧美精品免费久久| 搡老岳熟女国产| 日韩欧美在线乱码| 久久久色成人| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av电影不卡..在线观看| 狂野欧美激情性xxxx在线观看| 日本黄大片高清| 蜜桃久久精品国产亚洲av|