• <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
    免费看a级黄色片| 免费观看性生交大片5| 国产高清有码在线观看视频| 天堂俺去俺来也www色官网| 在线观看人妻少妇| 久久久久国产网址| 男人舔奶头视频| 日韩三级伦理在线观看| 色婷婷久久久亚洲欧美| 99热这里只有是精品50| 一级毛片黄色毛片免费观看视频| 99热网站在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲精品456在线播放app| 80岁老熟妇乱子伦牲交| 日韩欧美 国产精品| 丰满人妻一区二区三区视频av| 最近中文字幕高清免费大全6| 狂野欧美激情性bbbbbb| 久久久欧美国产精品| 身体一侧抽搐| 精品午夜福利在线看| av卡一久久| 美女视频免费永久观看网站| 国产探花极品一区二区| 精品国产乱码久久久久久小说| 久久久久久久午夜电影| 国产成人aa在线观看| 国产日韩欧美在线精品| av网站免费在线观看视频| 五月玫瑰六月丁香| 男人和女人高潮做爰伦理| 国产 一区 欧美 日韩| .国产精品久久| 亚洲欧美成人综合另类久久久| 中国国产av一级| 国产精品99久久99久久久不卡 | 亚洲精品视频女| 一级片'在线观看视频| 最近手机中文字幕大全| 久久久久久久久久久丰满| 国产成人精品福利久久| 亚洲国产精品专区欧美| 2021少妇久久久久久久久久久| 中文在线观看免费www的网站| 人人妻人人澡人人爽人人夜夜| 精品久久久久久久久亚洲| 视频中文字幕在线观看| 九九爱精品视频在线观看| 成人高潮视频无遮挡免费网站| 亚洲精品日韩av片在线观看| 午夜福利视频精品| 观看免费一级毛片| 天堂中文最新版在线下载 | 免费看日本二区| 麻豆成人av视频| 久久精品国产鲁丝片午夜精品| 久久久久国产网址| 亚洲av男天堂| 亚洲国产精品成人久久小说| 秋霞伦理黄片| 男人舔奶头视频| 激情五月婷婷亚洲| 国产中年淑女户外野战色| 小蜜桃在线观看免费完整版高清| 欧美日韩在线观看h| 欧美一区二区亚洲| 嘟嘟电影网在线观看| 99re6热这里在线精品视频| av播播在线观看一区| 深爱激情五月婷婷| 高清毛片免费看| 国产白丝娇喘喷水9色精品| 久久精品国产自在天天线| 亚洲欧美日韩卡通动漫| 又爽又黄无遮挡网站| 高清视频免费观看一区二区| 自拍偷自拍亚洲精品老妇| 亚洲av电影在线观看一区二区三区 | 97在线人人人人妻| 永久网站在线| 一级片'在线观看视频| 色婷婷久久久亚洲欧美| 亚洲av中文字字幕乱码综合| 在线a可以看的网站| 国产精品成人在线| 欧美少妇被猛烈插入视频| 97热精品久久久久久| 丝瓜视频免费看黄片| 国产白丝娇喘喷水9色精品| 一区二区av电影网| 在现免费观看毛片| 免费av观看视频| 国产精品偷伦视频观看了| 在线观看美女被高潮喷水网站| 色视频www国产| 欧美+日韩+精品| 一区二区三区精品91| 久久人人爽av亚洲精品天堂 | 中国美白少妇内射xxxbb| 久久久亚洲精品成人影院| 亚洲欧洲日产国产| 一级av片app| 特级一级黄色大片| 亚洲国产av新网站| 国产亚洲av嫩草精品影院| 在线观看一区二区三区激情| 日产精品乱码卡一卡2卡三| 嫩草影院精品99| 久久久久性生活片| 免费av不卡在线播放| 一级毛片 在线播放| 亚洲av免费高清在线观看| 国产免费又黄又爽又色| 国产亚洲5aaaaa淫片| 青青草视频在线视频观看| 一级av片app| 久久久久国产网址| 久久人人爽人人爽人人片va| 天堂俺去俺来也www色官网| 日韩国内少妇激情av| 99热国产这里只有精品6| 欧美xxxx黑人xx丫x性爽| 免费黄频网站在线观看国产| 交换朋友夫妻互换小说| 天美传媒精品一区二区| 香蕉精品网在线| 少妇高潮的动态图| xxx大片免费视频| 国产精品99久久99久久久不卡 | 国产精品久久久久久av不卡| 国产淫语在线视频| 亚洲精品乱码久久久久久按摩| 综合色av麻豆| 色视频www国产| 成人高潮视频无遮挡免费网站| 免费播放大片免费观看视频在线观看| 亚洲av免费高清在线观看| 久久这里有精品视频免费| 久久精品人妻少妇| 91久久精品国产一区二区成人| 欧美高清性xxxxhd video| 国产老妇女一区| 亚洲自偷自拍三级| 日韩一区二区三区影片| 色网站视频免费| 亚洲av福利一区| 国产成人aa在线观看| 欧美日本视频| 七月丁香在线播放| 少妇 在线观看| 精品少妇黑人巨大在线播放| 免费不卡的大黄色大毛片视频在线观看| 国产精品熟女久久久久浪| 久久这里有精品视频免费| 国产极品天堂在线| 嫩草影院新地址| 久久精品国产鲁丝片午夜精品| 永久网站在线| 舔av片在线| 亚洲精品国产av蜜桃| 人体艺术视频欧美日本| 综合色av麻豆| 亚洲天堂国产精品一区在线| 中文天堂在线官网| 在线观看免费高清a一片| 亚洲人成网站高清观看| 免费高清在线观看视频在线观看| 男男h啪啪无遮挡| 国产免费又黄又爽又色| 国产黄频视频在线观看| 成人无遮挡网站| 国产淫语在线视频| 久久精品国产自在天天线| 禁无遮挡网站| 国产精品人妻久久久久久| 国产精品无大码| 六月丁香七月| 国产在线男女| 99热网站在线观看| 免费黄色在线免费观看| 熟女电影av网| 国产精品女同一区二区软件| 国产乱人偷精品视频| 人妻夜夜爽99麻豆av| 欧美一区二区亚洲| 尤物成人国产欧美一区二区三区| 日韩欧美 国产精品| 综合色丁香网| 青春草视频在线免费观看| 一区二区三区精品91| 六月丁香七月| 91久久精品电影网| 亚洲精品久久午夜乱码| 午夜精品一区二区三区免费看| 亚洲精品日韩在线中文字幕| 日韩欧美精品免费久久| 欧美老熟妇乱子伦牲交| 国产成年人精品一区二区| 久久久久性生活片| 夫妻性生交免费视频一级片| 欧美一区二区亚洲| 亚洲人成网站高清观看| 99热国产这里只有精品6| 亚洲最大成人手机在线| 国产毛片在线视频| 国产精品国产三级国产av玫瑰| 亚洲成人久久爱视频| 日日摸夜夜添夜夜添av毛片| 国产成人免费观看mmmm| 国产乱来视频区| 久久女婷五月综合色啪小说 | 亚洲天堂国产精品一区在线| 国产极品天堂在线| 性色avwww在线观看| 成人欧美大片| 久久99热这里只有精品18| 国产老妇伦熟女老妇高清| 亚洲精华国产精华液的使用体验| 国产精品秋霞免费鲁丝片| 好男人在线观看高清免费视频| 免费看不卡的av| 亚洲天堂av无毛| 日本一二三区视频观看| 国产成人精品福利久久| 久久人人爽av亚洲精品天堂 | 成年女人看的毛片在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精华国产精华液的使用体验| 蜜臀久久99精品久久宅男| av在线天堂中文字幕| 日本三级黄在线观看| 国产免费福利视频在线观看| 少妇猛男粗大的猛烈进出视频 | 亚洲,一卡二卡三卡| 丝袜喷水一区| 视频区图区小说| 国产人妻一区二区三区在| 十八禁网站网址无遮挡 | 亚洲av成人精品一二三区| 国产精品熟女久久久久浪| 免费看光身美女| 一本色道久久久久久精品综合| a级毛色黄片| 波野结衣二区三区在线| 国产亚洲午夜精品一区二区久久 | 丰满人妻一区二区三区视频av| 国产69精品久久久久777片| 黄色日韩在线| 视频中文字幕在线观看| 秋霞伦理黄片| 国产精品国产三级国产av玫瑰| 麻豆精品久久久久久蜜桃| 大陆偷拍与自拍| 我要看日韩黄色一级片| 久久精品熟女亚洲av麻豆精品| 欧美日韩精品成人综合77777| 久久久久网色| 女人久久www免费人成看片| 一级毛片电影观看| 欧美一级a爱片免费观看看| 免费高清在线观看视频在线观看| 男人和女人高潮做爰伦理| 黄色一级大片看看| 亚洲最大成人av| 亚洲精品一区蜜桃| 亚洲自拍偷在线| 午夜福利视频1000在线观看| 看十八女毛片水多多多| 国产高清国产精品国产三级 | 全区人妻精品视频| 久久99热这里只频精品6学生| 日韩视频在线欧美| 亚洲美女视频黄频| 香蕉精品网在线| 婷婷色麻豆天堂久久| 亚洲综合精品二区| 一个人看视频在线观看www免费| av又黄又爽大尺度在线免费看| 美女视频免费永久观看网站| 久久99精品国语久久久| 亚洲人与动物交配视频| 午夜福利高清视频| 久久久色成人| 一级毛片黄色毛片免费观看视频| 亚洲精品影视一区二区三区av| 80岁老熟妇乱子伦牲交| 久久人人爽人人爽人人片va| 老女人水多毛片| 亚洲精品日韩av片在线观看| 国产精品一二三区在线看| 精品国产露脸久久av麻豆| 免费观看性生交大片5| 久久久欧美国产精品| 日韩大片免费观看网站| 亚洲精品国产av蜜桃| 国产极品天堂在线| 99热全是精品| 天天躁夜夜躁狠狠久久av| 国产免费视频播放在线视频| 亚洲人成网站在线观看播放| 国产黄a三级三级三级人| 国产在线一区二区三区精| 日本与韩国留学比较| 国产精品嫩草影院av在线观看| 欧美激情久久久久久爽电影| 欧美一区二区亚洲| 国产又色又爽无遮挡免| 麻豆国产97在线/欧美| 男女那种视频在线观看| 国产成人91sexporn| 欧美丝袜亚洲另类| 国内少妇人妻偷人精品xxx网站| 六月丁香七月| 成人亚洲精品一区在线观看 | 好男人视频免费观看在线| 人妻一区二区av| 777米奇影视久久| 国产91av在线免费观看| 少妇人妻一区二区三区视频| 青春草国产在线视频| 男女啪啪激烈高潮av片| 丝袜脚勾引网站| 最近中文字幕2019免费版| 午夜亚洲福利在线播放| 尾随美女入室| a级毛片免费高清观看在线播放| tube8黄色片| 国产成人91sexporn| 欧美 日韩 精品 国产| 久久精品人妻少妇| 国产成人aa在线观看| 亚洲欧美精品自产自拍| 成人无遮挡网站| 伦理电影大哥的女人| 日本av手机在线免费观看| 人人妻人人看人人澡| 国产亚洲91精品色在线| 中文在线观看免费www的网站| 美女主播在线视频| 热re99久久精品国产66热6| av网站免费在线观看视频| 国语对白做爰xxxⅹ性视频网站| 人妻一区二区av| 国产精品爽爽va在线观看网站| 夜夜爽夜夜爽视频| 成人国产av品久久久| 亚洲欧美一区二区三区黑人 | av在线天堂中文字幕| 久久精品熟女亚洲av麻豆精品| 免费电影在线观看免费观看| 国产乱人偷精品视频| 视频区图区小说| 亚洲av不卡在线观看| 少妇人妻一区二区三区视频| 国产一级毛片在线| 国产黄频视频在线观看| 少妇熟女欧美另类| 男女边摸边吃奶| 男插女下体视频免费在线播放| 国产精品一区二区三区四区免费观看| 一级毛片久久久久久久久女| 久久这里有精品视频免费| 日日撸夜夜添| 久久精品夜色国产| 日韩国内少妇激情av| 亚洲av国产av综合av卡| 一级毛片久久久久久久久女| 久久精品国产鲁丝片午夜精品| 成人免费观看视频高清| 色吧在线观看| 亚洲欧美中文字幕日韩二区| 亚洲人成网站在线观看播放| 久久国产乱子免费精品| 嫩草影院新地址| 国产成人精品一,二区| freevideosex欧美| 成人无遮挡网站| 免费看不卡的av| 少妇人妻久久综合中文| 免费在线观看成人毛片| 成人欧美大片| 大话2 男鬼变身卡| 久久久精品欧美日韩精品| 日本一二三区视频观看| 国产一区有黄有色的免费视频| 亚洲成人精品中文字幕电影| 性色avwww在线观看| 欧美高清成人免费视频www| 日本色播在线视频| 99热全是精品| 国产真实伦视频高清在线观看| 国产69精品久久久久777片| 永久免费av网站大全| 插阴视频在线观看视频| 插逼视频在线观看| 一本久久精品| 午夜福利在线观看免费完整高清在| 热99国产精品久久久久久7| 欧美三级亚洲精品| 热99国产精品久久久久久7| 寂寞人妻少妇视频99o| 一本久久精品| 亚洲欧美精品自产自拍| 日韩亚洲欧美综合| 人妻夜夜爽99麻豆av| 亚洲国产精品成人久久小说| 免费高清在线观看视频在线观看| 成年免费大片在线观看| 2021天堂中文幕一二区在线观| 午夜激情福利司机影院| 2021天堂中文幕一二区在线观| 欧美丝袜亚洲另类| 国产一区二区在线观看日韩| 亚洲国产精品国产精品| 日日撸夜夜添| 亚洲在线观看片| 自拍欧美九色日韩亚洲蝌蚪91 | 女人被狂操c到高潮| av播播在线观看一区| 国产一级毛片在线| 草草在线视频免费看| 伊人久久国产一区二区| 日日撸夜夜添| 国产日韩欧美在线精品| 久久久国产一区二区| 国产色婷婷99| 日韩av免费高清视频| 波多野结衣巨乳人妻| 一级毛片久久久久久久久女| 国产精品.久久久| 91精品伊人久久大香线蕉| 人人妻人人澡人人爽人人夜夜| 美女内射精品一级片tv| 国产av码专区亚洲av| 日韩亚洲欧美综合| 国产黄色免费在线视频| 国产中年淑女户外野战色| 久久精品人妻少妇| 美女视频免费永久观看网站| 日韩制服骚丝袜av| 国产色婷婷99| 偷拍熟女少妇极品色| 日韩国内少妇激情av| 国产成人a∨麻豆精品| h日本视频在线播放| 久久久久九九精品影院| 永久免费av网站大全| 国产精品国产三级国产专区5o| 国产精品福利在线免费观看| 在线亚洲精品国产二区图片欧美 | 青春草视频在线免费观看| 又粗又硬又长又爽又黄的视频| 国产综合懂色| 日本猛色少妇xxxxx猛交久久| 麻豆国产97在线/欧美| 大又大粗又爽又黄少妇毛片口| 国产黄a三级三级三级人| 欧美变态另类bdsm刘玥| 国产视频首页在线观看| 国产毛片在线视频| 国产免费视频播放在线视频| 日本-黄色视频高清免费观看| 最近2019中文字幕mv第一页| 国产午夜精品久久久久久一区二区三区| 亚洲国产成人一精品久久久| 日韩大片免费观看网站| 男女无遮挡免费网站观看| 一边亲一边摸免费视频| 亚洲精华国产精华液的使用体验| 国产黄片视频在线免费观看| av.在线天堂| 国产精品一区www在线观看| 久久久午夜欧美精品| 成年免费大片在线观看| 日韩中字成人| 国产一区有黄有色的免费视频| 国产亚洲最大av| 内射极品少妇av片p| 国产黄片视频在线免费观看| 色网站视频免费| 久久久色成人| 国产精品无大码| 亚洲精品一区蜜桃| 亚洲欧洲日产国产| 99re6热这里在线精品视频| 好男人视频免费观看在线| 欧美精品一区二区大全| 99热全是精品| 少妇丰满av| 丝袜喷水一区| 精品一区在线观看国产| 亚洲最大成人中文| 亚洲综合精品二区| 婷婷色麻豆天堂久久| 在线亚洲精品国产二区图片欧美 | a级一级毛片免费在线观看| 在线观看美女被高潮喷水网站| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品成人久久小说| 精品人妻视频免费看| 丝瓜视频免费看黄片| av又黄又爽大尺度在线免费看| 午夜视频国产福利| 色视频在线一区二区三区| 国产精品.久久久| 日韩亚洲欧美综合| 97人妻精品一区二区三区麻豆| 亚洲天堂国产精品一区在线| 亚洲精品国产色婷婷电影| 天天一区二区日本电影三级| 午夜福利视频1000在线观看| 中文欧美无线码| 尤物成人国产欧美一区二区三区| 精品久久久久久久久av| 亚洲国产精品999| 毛片一级片免费看久久久久| 大话2 男鬼变身卡| 国产精品蜜桃在线观看| 免费电影在线观看免费观看| 97超碰精品成人国产| 性插视频无遮挡在线免费观看| 我的老师免费观看完整版| 热re99久久精品国产66热6| 99热网站在线观看| 久久久精品免费免费高清| 色视频在线一区二区三区| 成人亚洲欧美一区二区av| 精品一区二区三区视频在线| 寂寞人妻少妇视频99o| 久久精品久久久久久久性| 国产伦精品一区二区三区视频9| 欧美97在线视频| 国产综合精华液| 成人亚洲欧美一区二区av| 乱系列少妇在线播放| 黄片无遮挡物在线观看| 又黄又爽又刺激的免费视频.| 久久人人爽av亚洲精品天堂 | 777米奇影视久久| 伦精品一区二区三区| 亚洲精品国产av成人精品| 亚洲精品自拍成人| 中文字幕制服av| av又黄又爽大尺度在线免费看| 欧美老熟妇乱子伦牲交| 精品久久久精品久久久| 黄片无遮挡物在线观看| 亚洲人成网站高清观看| 中文精品一卡2卡3卡4更新| 人妻一区二区av| 婷婷色综合www| 国产午夜精品久久久久久一区二区三区| tube8黄色片| 噜噜噜噜噜久久久久久91| 少妇高潮的动态图| 国产精品99久久99久久久不卡 | 国产色婷婷99| av在线播放精品| 精品人妻偷拍中文字幕| 久久久久久久大尺度免费视频| 97超视频在线观看视频| 另类亚洲欧美激情| 97热精品久久久久久| 汤姆久久久久久久影院中文字幕| 久久精品国产鲁丝片午夜精品| 97在线视频观看| 国产精品嫩草影院av在线观看| av.在线天堂| 久久精品久久久久久久性| 18禁裸乳无遮挡动漫免费视频 | 亚洲国产最新在线播放| 亚洲真实伦在线观看| 在线观看人妻少妇| 最近中文字幕高清免费大全6| 成人毛片a级毛片在线播放| 日本三级黄在线观看| 国产精品精品国产色婷婷| 一级毛片黄色毛片免费观看视频| 国产精品久久久久久久久免| 一边亲一边摸免费视频| av福利片在线观看| 秋霞在线观看毛片| 中国三级夫妇交换| 伊人久久精品亚洲午夜| 好男人在线观看高清免费视频| 中文欧美无线码| 国产免费又黄又爽又色| 天天躁日日操中文字幕| 神马国产精品三级电影在线观看| 精品亚洲乱码少妇综合久久| 日韩强制内射视频| av福利片在线观看| 性色av一级| 国产亚洲91精品色在线| 成人免费观看视频高清| 国产精品一二三区在线看| 大码成人一级视频| 精品少妇黑人巨大在线播放| 水蜜桃什么品种好| 久久久久久久久大av| 亚洲精华国产精华液的使用体验| 午夜老司机福利剧场| 国产精品无大码| 日韩一区二区视频免费看| 毛片一级片免费看久久久久| 91午夜精品亚洲一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 在线观看美女被高潮喷水网站| 久久久久久九九精品二区国产| 高清毛片免费看| 噜噜噜噜噜久久久久久91| 99精国产麻豆久久婷婷| 免费看日本二区| 精品一区在线观看国产|