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

    Valley-dependent topological edge states in plasma photonic crystals

    2023-03-15 00:53:54JianfeiLI李健飛ChenZHOU周晨JingfengYAO姚靜鋒ChengxunYUAN袁承勛YingWANG王瑩ZhongxiangZHOU周忠祥JingwenZHANG張景文andAnatolyKUDRYAVTSEV
    Plasma Science and Technology 2023年3期
    關(guān)鍵詞:王瑩

    Jianfei LI (李健飛), Chen ZHOU (周晨), Jingfeng YAO (姚靜鋒),2,*,Chengxun YUAN (袁承勛),2,*, Ying WANG (王瑩),2,Zhongxiang ZHOU (周忠祥),2, Jingwen ZHANG (張景文) and Anatoly A KUDRYAVTSEV

    1 School of Physics, Harbin Institute of Technology, Harbin 150000, People’s Republic of China

    2 Heilongjiang Provincial Key Laboratory of Plasma Physics and Application Technology,Harbin 150000,People’s Republic of China

    3 Physics Department, St.Petersburg State University, St.Petersburg 198504, Russia

    Abstract Plasma photonic crystals designed in this paper are composed of gas discharge tubes to control the flow of electromagnetic waves.The band structures calculated by the finite element method are consistent with the experimental results which have two distinct attenuation peaks in the ranges of 1-2.5 GHz and 5-6 GHz.Electromagnetic parameters of the plasma are extracted by the Nicolson-Ross-Weir method and effective medium theory.The measured electron density is between 1 × 1011cm-3 and 1 × 1012 cm -3, which verifeis the correctness of the parameter used in the simulation,and the collision frequency is near1.5 ×1010 Hz.As the band structures are corroborated by the measured scattering parameters, we introduce the concept of photonic topological insulator based on the quantum Valley Hall effect into the plasma photonic crystal.A valley-dependent plasma photonic crystal with hexagonal lattice is constructed, and the phase transition of the valley K(K ′ ) occurs by breaking the spatial inversion symmetry.Valley-spin locked topological edge states are generated and excited by chiral sources.The frequency of the non-bulk state can be dynamically regulated by the electron density.This concept paves the way for novel,tunable topological edge states.More interestingly,the Dirac cone is broken when the electron density increases to 3.1 × 1012 cm-3 ,which distinguishes from the methods of applying a magnetic field and changing the symmetry of the point group.

    Keywords: plasma photonic crystal, valley-dependent, topological state, electron density

    1.Introduction

    Plasma is widely used in stealth applications, broadband absorber, and plasma antenna fields by virtue of its extraordinary physical properties such as complex permittivity,electrical conductivity, and reconfigurability [1-3].The electron density and collision frequency substantially affect the reflection, absorption, and refraction of electromagnetic waves when electromagnetic waves propagate in plasma.The combination of plasma and metamaterials shows promising applications in dynamically controlling the flow of electromagnetic waves.Hojo first proposed the concept of plasma photonic crystals, arranging the plasma and quartz periodically in one-dimension.The frequency gap and cut-off band were calculated by solving the Maxwell equations, and the transmittance profile could be modulated by the electron density [4, 5].This successfully exploited the merits of plasma and dielectric photonic crystals.Because of the tunable and reconfigurable properties of plasma, the study of plasma photonic crystals has been extended from microwave to terahertz bands, and the research dimension has also been extended from one-dimension to three-dimension theoretically [6, 7].In the experiment, Sakai constructed twodimensional plasma photonic crystals in experiments that included plasma columns generated by the dielectric barrier discharge and the background air.The transmittance had a clear attenuation peak at 73.83 GHz when the transverse electric(in-plane electric field)mode was excited,which was consistent with the results of theoretical calculations using the plane-wave expansion method [8, 9].In addition, they designed plasma metamaterials where double metallic helices were inserted in the waveguide periodically to produce negative permeability.At the same time, the metal helices were applied with low-frequency high voltage to produce microplasma (i.e.negative permittivity at specific frequencies).The combination of the two designs produced a double negative refraction property [10].A tunable negativerefractive-index device was designed by Cappelli, which consisted of plasma discharge tubes and double split-ring resonators.Negative permittivity and permeability were achieved in the range of 2-2.47 GHz, and the transmittance increased with increasing current over a specific range [11].Liu et al designed a plasma photonic crystal with annular lattices, a self-organized form generated by dielectric barrier discharge.The transition from annular lattice to core-annular and concentric-annular lattice was regulated by the applied voltage,and the forbidden band between 28 and 30 GHz was verified experimentally [12].Fan et al used dielectric barrier discharge to form honeycomb plasma photonic crystal.The lattice structure can be regulated by the applied voltage,which is a form of self-organization [13].In the previous study of our group, we used dielectric columns and a plasma background to form a plasma photonic crystal with a square lattice.By dynamically modulating the electron density of the plasma, a triple degenerated Dirac cone was obtained at the centre of the Brillouin zone.When electromagnetic waves are incident at the frequency of the Dirac point, the plasma photonic crystal has zero permittivity and permeability [14].

    With the successful development of photonic crystals,the photonic topological insulator analogous to electronic systems has attracted extensive research.Topologically-protected edge modes based on the quantum Hall effect and the quantum-spin Hall effect have promising applications in the field of waveguides, coupled resonators, chip development,and biosensors [15-17].Besides, the topological effect can also be achieved by adjusting the valley degrees of freedom of the valley photonic crystal, which has been proposed in recent years.By breaking the spatial inversion symmetry,topological phase transitions arise at theK(K′ ) valley in the Brillouin zone,and valley-spin locked topological edge states were detected theoretically and experimentally [18-20].However, topological edge states in all-dielectric photonic crystals can only be observed in a specific frequency range with fixed lattice structures and electromagnetic parameters,and the electromagnetic modes are not switchable.Several materials such as liquid crystals and barium titanate have been used to explore the reconfigurable topological states[21,22].Qiu et al used plasmonic to form honeycomb lattices.Topological angular states were verified by changing the radius of the column or the position of the plasmonic [23].Zhao et al filled nematic liquid crystals into dielectric columns to construct valley photonic crystals while topological phase transitions were produced by changing the position of the liquid crystal columns.In theory, the anisotropic permittivity of liquid crystals can be used to modulate and encode topological states [24].Meanwhile, silicon columns and liquid crystal backgrounds were combined into valley photonic crystals designed by Shalaev, which realized the modulation of topological states by changing the anisotropic permittivity of liquid crystals [25].

    Figure 1.(a) The experimental setup of a plasma photonic crystal with square lattices, (b) the measured S21 parameters at different discharge currents.

    In this work,we investigate valley-dependent topological edge states in plasma photonic crystals composed of gaseous plasma, dielectric layers, and air background.Firstly, a twodimensional plasma photonic crystal is constructed with discharge tubes,and the transmittance of electromagnetic waves is measured to verify the simulated band structures, which confirms the feasibility of constructing photonic crystals using gaseous plasma.Next,we test the scattering parameters(S-parameters) of the discharge plasma within the transverse electro-magnetic (TEM) cell and extract the electromagnetic parameters via the Nicolson-Ross-Weir (NRW) method and effective medium theory.Finally, the tunable topological edge states are realized theoretically by constructing plasma valley photonic crystals.

    2.Plasma photonic crystal with square lattice

    Plasma photonic crystals enable the modulation of the electromagnetic wave since the electromagnetic parameters can be dynamically controlled by external conditions such as applied voltages, gas pressure and magnetic fields.At the same time,the band structures of plasma photonic crystals can be calculated precisely as in the case of all-dielectric photonic crystals.In particular,it is necessary to accurately acquire the electromagnetic parameters of plasma.In plasma physics, a common way to generate plasma is capacitively coupled discharge.The fluorescent lamp provides beneficial conditions for discharge plasma as it is filled with an Ar-Hg mixture, and the pressure is 3 Torr.A two-dimensional plasma photonic crystal with square lattice is composed of 28 lamps.The experimental setup is shown in figure 1(a),and the primitive cell is a square lattice with a lattice constant of 30 mm as shown in the inset.Two horn antennas with 1-14 GHz are placed on both sides of the plasma photonic crystal to measure the S-parameters, and the measurement system is surrounded by wave-absorbing sponges.

    All the lamps are connected in series to ensure that each lamp has the same discharge state and produces the same electron density.It is known that the lamp has a low maintained voltage but a high excited voltage.The power source of CTP-2000K is used to support the discharge of the lamps.The pulsed high frequency oscillator is adopted to generate high voltages of 15 kV and very low currents.Then all the lamps can be ignited by the pulsed high frequency oscillator,and the discharge current is regulated by the power source.It must be noted that the power source line needs a diode in series to isolate the high voltage generated by the pulsed high frequency oscillator.The measured results at different discharge currents are shown in figure 1(b), two distinct attenuation peaks appear in the ranges of 1-2.5 GHz and 5-6 GHz, and the magnitude of attenuation increases with increasing discharge current.However,in the high-frequency region, the attenuation peak shifts rather than changes in magnitude at different currents, which verifies the existence of two forbidden bands in the plasma photonic crystal.

    Next, the experiments to precisely measure the plasma electron density and collision frequency are performed in the TEM cell.The NRW method is used to extract the electron density and collision frequency from the measured scattering parameters [26, 27].This approach of electromagnetic wave diagnosis does not require the electron kinetic temperature,the electron energy distribution function,or the collision cross section[28].Figure 2(a)shows the top view of the measuring apparatus, which consists of a TEM cell and eight mercurybased fluorescent lamps.The length of the TEM cell is 430 mm, and the height is 60 mm.At the centre of the TEM cell, the transverse electromagnetic wave is formed with a length of 152 mm [29, 30].Eight fluorescent lamps with an outer radius of 6.33 mm and an inner radius of 5.6 mm are neatly arranged in the central area of the TEM cell, and the energy flow of electromagnetic waves is perpendicular to the sample as shown in figure 2(b).The red line represents the electric field and the blue line represents the magnetic field.So, the effective medium theory base on the Garnett rule can suitably extract the complex permittivity and permeability of plasma [31, 32].The effective sample is shown as the black dashed line in the figure, which consists of plasma, quartz tube, and air.The volume ratio of the plasma to the effective sample isf= 0.5188.NRW method is commonly used in testing electromagnetic parameters of the sample.In this paper, the TEM cell can be viewed as a coaxial transmission line which accurately measures scattering parameters in a closed field compared to free space.According to Nicolson’s algorithm,the permittivity and permeability of the sample can be derived from the measured S-parameters [26, 27]:

    Figure 2.(a)Experimental diagram of the S-parameter measurement in the TEM cell,(b)a schematic of the effective sample in the TEM cell,(c)the measured S21 over a wide frequency ranges,(d)the measured S11 over a wide frequency range,(e)the electron density calculated by the NRW method and effective medium theory, (f) the collisional frequency calculated by the NRW method and effective medium theory.

    where reflection coefficient isΓ =(∣ Γ∣≤ 1),and the relation between reflection coefficient and scattering parameters isX=are permeability, permittivity, length of the sample, radian frequency of the incident wave, and the wavelength in the sample, respectively.In the experiment, eight lamps are connected in series to ensure the same discharge state of each fluorescent lamp and are driven by the CTP-2000K.

    The measured S21is shown in figure 2(c).When the root mean square of the current reaches 38 mA, all lamps are discharged.However, S21is attenuated by 20 dB on average as shown in the green line,which is not only the effect of the plasma but more importantly the Mercury vapor.The extracted electron density and collision frequency from this parameter do not truly reflect the properties of the plasma.Then,the collision frequencies calculated for this case are too large, such as 29-74 GHz, as seen in the [33].After that, the attenuation of S21gradually increases with the increasing current.It is worth noting that the discharge current increases to 128 mA and a distinct cut-off region appears at less than 5 GHz, in which the permittivity of plasma is negative, as shown in the left part of the red line.As seen in figure 2(d),S11is almost the same when theIrms=64 mA andIrms=128 mA,which means that the collisional absorption of electromagnetic waves by the plasma is weakened.Then,the effective medium theory is combined with the NRW algorithm to accurately diagnose the plasma parameters.To exclude the effect of mercury vapor,the complex permittivity of the effective sample is used as the background medium when the current is 38 mA.The classical Maxwell Garnett approximation for extracting the complex permittivity of plasma is adopted as follows [32]

    whereεeffis the total effective permittivity of the sample extracted from the S-parameters using the NRW method,ε1is the permittivity of plasma,ε2is the permittivity of background medium, andfis the volume ratio of the plasma to the background medium (in this casef= 0.5188).The measured S-parameters are used as a background medium to eliminate the effect of mercury vapor whenIrms=38 mA,and the complex permittivity of the plasma is obtained.

    Then, the electron density and collision frequency of plasma can be derived from Drude equation:

    whereωpeis the Langmuir frequency andνcis the electronneutral collision frequency.The Langmuir frequency can be denoted aswherenis the electron density,eis the unit charge, andmeis the mass of the electron.

    In the previous literatures, the analysis of the collision frequency is ambiguous when electromagnetic waves propagate in the plasma generated by fluorescent lamps.Here,we can precisely extract the plasma parameters using effective medium theory.The plasma cut-off frequency is around 5 GHz due to the low transmittance in the region less than that frequency.From figure 2(e), it can be seen that the electron density increases with current.When the current reaches 128 mA, the electron density is around 8 × 1011cm-3.The collision frequency is around 15 GHz as shown in figure 2(f).It is found that the electron density and collision frequency of the plasma fluctuate with the radian frequency of the incident wave.Because the measured S-parameters consist of a real and an imaginary part.For example, the real part of S21iscos (?2-?1) and the imaginary part of S21issin (?2-?1) where?nis the phase that contains the radian frequency.The obtained permittivity and permeability of plasma must fluctuate with frequency.Therefore, the extracted electron density and collision frequency fluctuate in a small range with the frequency of incident wave.

    The calculation of the dispersion relationship is extremely valuable and can accurately predict the transport properties of electromagnetic waves.The critical issue is that experimental results are necessary to validate the simulation.It is reasonable to ignore the collision frequency in substance because the electromagnetic energy will be concentrated at the high permittivity region when the electromagnetic wave propagates in photonic crystals[34].The band structure of the plasma photonic crystal is calculated using the finite element method (FEM) by COMSOL software.Before discharge, the band structure of the photonic crystal composed of fluorescent lamps is shown in the left part of figure 3.There is no forbidden band in the Γ-X direction which is the direction of electromagnetic wave propagation in the experiment.When the plasma is added to the lamps, the calculation of the band structure becomes complicated due to the dispersive properties of plasma.The eigenfrequency interface and global ODE and DAE interfaces must be combined to find the eigenvalues in the primitive cell.

    Here we use the collisionless cold plasma approximation whose permittivity is expressed asεr(ω) = 1-and the plasma electron density of 5 × 1011cm-3is adopted.The calculation results are shown in the right part of figure 3.It is obvious to see two forbidden bands in the ranges of 0-1.9 GHz and 4.9-5.5 GHz, which is highly consistent with the experimental results.An omnidirectional forbidden band appears in the low-frequency band, which is a typical feature of plasma photonic crystals and has been extensively studied in[9,35-37].The collisionless cold plasma approximation is reasonable through theoretical and experimental studies of plasma photonic crystals with a square lattice.These results show that the plasma photonic crystal composed of gaseous plasma has a distinct band gap characteristic, which will provide the basis for tunable devices based on plasma.Meanwhile, the novel effects in optics such as topology,pseudo-diffusive transport, and zero refraction have yet to be investigated in plasma photonic crystals.Next, we will focus on topological edge states based on the valley Hall effect.

    Figure 3.Band structures of plasma photonic crystal calculated by the finite element method at different electron densities.

    Figure 4.(a)The plasma photonic crystal with a honeycomb lattice,(b) the calculated band structure at 8 × 1011 cm -3.

    3.Valley photonic crystal based on gaseous plasma

    In the previous section, we obtained the electron density in the fluorescent lamps by the NRW method and demonstrated that the collisionless cold plasma approximation is reasonable through S21in plasma photonic crystal.Then, the valley plasma photonic crystal is designed to realize the topological edge states, which possess dynamically tunable properties.We design the plasma photonic crystal with the honeycomb lattice as shown in figure 4(a).The grey area is the primitive cell with a lattice constant ofa= 30 mm.Each column consists of a ring-shaped solid dielectric material with a permittivity of 11.7, whose inner radius is0.154aand outer radius of0.2a, and gaseous plasma.

    They are distributed periodically in the x and y directions,and the length is infinite in the z direction.The electron density used in the simulation is 8 × 1011cm-3.Columns of types 1,2,and 3 are connected in series,and columns of types 4, 5, and 6 are connected in series, which ensures that they can be controlled separately.In this work, the TM mode of electromagnetic wave with out-of-plane electric fields is considered,and the band structure of plasma photonic crystal is calculated as shown in figure 4(b).

    TheC3vsymmetry and time-reversal symmetry are preserved when all tubes are filled with gaseous plasma,and the Dirac cone is formed at theKvalley andK′ valley in the reduced Brillouin zone of plasma photonic crystal.The horizontal axis represents the wave vector in the reciprocal lattice space.For the calculation,the wave vector k is scanned at the high symmetry points of the reduced Brillouin zone to obtain the eigenfrequencies.For more intuitive results, the vertical axis uses frequency units rather than the normalized frequencyωa2πcwhereωis the radian frequency andcis the light speed.When the dispersive material is introduced into the calculation of the band structures,an omnidirectional band gap appears in the entire Brillouin zone below 2.08 GHz.

    Valley cleavage occurs when the spatial inversion symmetry is broken, whose lower left column or upper right one in the primitive cell is filled with gaseous plasma,as shown in figures 5(a) and (c).In this case, an omnidirectional bandgap(4-4.3 GHz) arises here because the symmetry of the honeycomb lattice is reduced fromC3vtoC.3The width of the bandgap can be modulated dynamically by the electron density of the plasma.The eigenstates of band 1 and band 2 atKandK′ valleys are shown in figure 5(b) when columns of types 1, 2, and 3 are filled with plasma.The coloured area represents the phase distribution ofEz.The arrow direction represents the Poynting vector direction and the size represents the intensity.

    When the Dirac cone is opened by the changes of spatial inversion symmetry, theK(K′ ) points of the band structures will appear to have phases with either left- or right-handed polarization.At this time, there is circularly polarized orbital angular momentum at the valley, which can be described by topological charge.Topological charge is defined bywhere arg(E z) represents the phase of electric field andLrepresents the closed path around the singularity of electric field distributions.The right-handed circular polarized (RCP,l= 1) angular momentum at theKvalley and the left-handed circular polarized (LCP,l= -1)angular momentum at theK′ valleys are identified in band 1(see figure 5(b)).As for band 2,the LCP angular momentum at theKvalley and RCP angular momentum at theK′ valleys occur.When the filled form of the plasma is changed to its opposite, the phase transition occurs as shown in the eigenstate in figure 5(d).The topological charge describes the vortex properties of the eigenmode.The value represents the number of times the phase changes from0 to2π, and the symbol represents the direction of phase change.For example, the phase variation around the electric field singularity for K point of band 1 is shown in figure 5(e)when lower left column is filled with plasma.The horizontal axis represents the position around the electric field singularity,and the vertical axis represents the wrapped phase ofEz.It can be seen that the phase positively changes by one cycle along the closed path,which corresponds to a topological charge of+1.At the same time, the direction of the Poynting vector is clockwise,which maintains a better consistency.As for the K point of band 2, the phase negatively changes by one cycle along the closed path corresponding to the topological charge of-1, and the direction of the Poynting vector is counter clockwise as shown in figure 5(b).

    Figure 5.(a)The calculated band structure when plasma fills the lower left column,(b)phase distribution and Poynting vector at valley K( K′ )when plasma fills the lower left column,(c)the calculated band structure when plasma fills the upper right column,(d)phase distribution and Poynting vector at valley K(K ′ ) when plasma fills the upper right column, (e) when lower left column is filled with plasma, the phase variation around the singularity for K point of band 1, (f) when lower left column is filled with plasma, the phase variation around the singularity for K point of band 2.

    Figure 6.(a) The project band structure of supercell, (b) the z-component of electric field distribution when the frequency of the electromagnetic wave is 4.16 GHz,(c)the non-bulk state varies with plasma electron density,(d)a point-like chiral source with a clockwise decreasing phase is placed at the interface, (e) a point-like chiral source with a clockwise increasing phase is placed at the interface.

    We define structure 1,which has the lower left column in the primitive cell filled with plasma.The structure 2 denotes that the upper right column in the primitive cell is filled with plasma.Splicing structures 1 and 2 together will support a valley-dependent topological edge state.The supercell is constructed as shown in the inset of figure 6(a), with the upper part being plasma-filled columns of types 1, 2, and 3 and the lower part being plasma-filled columns of types 4, 5,and 6.The Floquet conditions are applied to the boundary in the x-direction of the supercell, and the scattering boundary conditions are used for the boundary in the y-direction.Then,the project band structure is calculated from the supercell(see figure 6(a)).The black line represents the dispersion relationship of the passing bands (bulk states) in the plasma photonic crystal.In addition, the non-bulk state arises in the common bandgap as the blue line.For figure 6(b),we choose 4.16 GHz on the non-bulk state and can see that the electric field is mainly localized at the interface of the two structures,while it decays rapidly in the interior.When the wave vector is fixed tokx=0.22(the red point in figure 6(a)), the frequency of the topological edge state is regulated dynamically by the electron density as shown in figure 6(c).The topological properties are still maintained due to the change in electron density that does not cause the phase transition to occur.The eigenfrequency of the non-bulk state increases with the electron density, which expands the scope of application.To excite the valley-dependent topological edge states,point-like chiral sources are constructed by four antennas with the phases set to 0,π, 3When the phases of the four antennas decrease clockwise with an excitation frequency of 4.16 GHz,electromagnetic energy is transmitted in the left direction along the edge as shown in the Ezdistribution of figure 6(d).The path of electromagnetic wave propagation in the backward direction is suppressed.Also,the electromagnetic energy does not leak into the interior of the plasma photonic crystal, which gives the system a high signal-to-noise ratio.However, the electromagnetic wave propagation in the right direction is excited when the phase of the antenna increases clockwise as shown in figure 6(e).The sources with different chirality excite electromagnetic waves propagate in opposite directions.In application areas, when electromagnetic waves propagate through the device,they are strongly reflected when encountering impurities or defects,thus affecting the transmission efficiency.Topological edge states possess the properties of unidirectional propagation,backscatter suppression, and immunity to impurities, which can solve the problem perfectly.Based on the topological edge state theory, topological lasers, optical waveguides,optical isolators,and optical modulators have been developed one after another.The plasma photonic crystal designed in this paper not only realizes topological edge states but also has tunable features, which will have a wide prospective application in tunable microwave devices.

    Figure 7.(a)When the electron density is raised to 8 × 1012 cm -3, the Dirac cone is broken and an omnidirectional band gap appears,(b)the critical density at which the Dirac point is broken,(c)the calculated band structure when plasma fills the lower left column at 8 × 1012 cm-3,(d) the edge state without topological protection.

    By modulating the electron density in a small range, we obtain tunable topological edge states and the topological properties are preserved.But the results change dramatically when the electron density is extremely high.As is known to us, the photonic crystal with honeycomb lattice has a linear dispersion relation at theK(K′ ) point, which is known as Dirac cone.The Dirac cone has good stability and can maintain its existence even when εxxand εyyare very different.The topological phase transition must be accompanied by the closing and reopening of the Dirac cone.A commonly method to break the Dirac cone is to apply a magnetic field outside a photonic crystal made of magneto-optical material,which breaks the time-reversal symmetry of the system [38].Another mechanism is to change the honeycomb lattice to a triangular lattice, which makes the spatial inverse symmetry of the system break.This also causes the Dirac cone to open[16].In the previous section, the Dirac cone is formed atK(K′ ) point when all tubes are filled with plasma in a primitive cell.When the electron density increases by an order of magnitude, the Dirac cone is opened while the band structure is raised overall (see figure 7(a)).The Dirac point changes when the electron density varies continuously as shown in figure 7(b).

    From the results, the critical electron density that makes the Dirac cone break is 3.1 × 1012cm-3.The permittivity of the plasma at this point is -19.24.Therefore, the negative permittivity can be used as an alternative method to break the Dirac cone independently of applying a magnetic field or changing the lattice structure.However, the band gap obtained by increasing the electron density is not topologically nontrivial.The specific proof is shown in figure 7(c).We increase the electron density of figure 5(a) to 8 × 1012cm-3.It can be seen from the results that the band gap is widened and its topological properties are changed.Since the Poynting vector at point K in the two bands does not have symmetry as in figure 5(b).When the chiral source is placed at the interface, the electromagnetic wave propagates along two directions at an electron density of 8 × 1012cm-3as shown in figure 7(d).The edge states in this case are defective states and not protected by topology.Although the method of breaking the Dirac cone by increasing the electron density does not lead to a topological phase transition of the system, the effect caused by the negative permittivity deserves to be studied in depth.

    4.Conclusions

    The plasma photonic crystal is the intersection of plasma and photonic crystal disciplines, which both have the properties of photonic crystals and plasma.It can be used for tunable filters and waveguides, etc by changing the discharge plasma parameters and the applied magnetic field.In past research,the band structures of plasma photonic crystals have been focused on different primitive cells, defect structures, and spatial dimensions, which need to be expanded in terms of theory and application.This work is devoted to investigating the tunable topological states in plasma photonic crystals that cannot be achieved in all-dielectric photonic crystals.

    We experimentally measure the electromagnetic parameters of the discharge plasma by the NRW method and TEM cell.The electron density extracted by the effective medium theory is between 1 × 1011cm-3and 1 × 1012cm-3,and the collision frequency between electrons and neutral atoms is near1.5 ×1010Hz.The precisely measured collision frequencies are much lower than the results in the literature such as 9× 1010Hzby excluding the absorption effect of mercury vapor.After that, 28 discharge tubes form a twodimensional plasma photonic crystal while the transmittance is measured at different discharge currents.The band structure of the plasma photonic crystal was calculated by COMSOL software, demonstrating the band gaps in the ranges of 0-1.9 GHz and 4.9-5.5 GHz in the direction of wave propagation.The collisionless cold plasma approximationis used in the simulation, and the electron density is 8 × 1011cm-3as measured by the TEM cell.The experimental results agree very well with the simulation results, which proves the rationality of this approximation.Finally, we design a honeycomb lattice of plasma photonic crystals theoretically.The valley-dependent topological phase transition is achieved by breaking the C3vsymmetry whose position of plasma is changed in the primitive cell.The edge states are generated at the interface of two plasma photonic crystals with different valley topological properties,and the direction of wave propagation is controlled by the chiral source.The frequency corresponding to the topological edge states can be dynamically regulated by the plasma parameter.More importantly, there is a critical electron density to break the Dirac cone, which originates from the negative permittivity of the plasma.This is worthy of further study in manipulating the flow of electromagnetic waves.

    Acknowledgments

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

    猜你喜歡
    王瑩
    鋼琴性能對音樂創(chuàng)作風格的影響
    音樂探索(2022年2期)2022-05-30 21:01:37
    巧用比較策略,突破學生的學習難點
    Machine learning combined with Langmuir probe measurements for diagnosis of dusty plasma of a positive column
    王瑩作品
    王瑩作品賞析
    王瑩作品
    王瑩作品賞析
    萊儷青年藝術(shù)獎獲獎者 王瑩:《租賃一平方米》的力量
    Theoretical research on the transport and ionization rate coefficients in glow discharge dusty plasma
    王瑩作品
    丁香六月天网| 美女午夜性视频免费| 母亲3免费完整高清在线观看| 免费观看av网站的网址| 国产在线免费精品| 亚洲人成77777在线视频| 国产亚洲欧美在线一区二区| 日韩熟女老妇一区二区性免费视频| 美女高潮到喷水免费观看| 国产亚洲欧美在线一区二区| 老司机影院毛片| 叶爱在线成人免费视频播放| 亚洲伊人久久精品综合| 日日爽夜夜爽网站| 亚洲欧美日韩另类电影网站| 久久久久久人人人人人| 尾随美女入室| 母亲3免费完整高清在线观看| 中文精品一卡2卡3卡4更新| 亚洲中文日韩欧美视频| 欧美 日韩 精品 国产| 亚洲第一青青草原| 一区二区三区精品91| 午夜久久久在线观看| 免费在线观看视频国产中文字幕亚洲 | 亚洲 欧美一区二区三区| 夜夜骑夜夜射夜夜干| 午夜91福利影院| 亚洲av电影在线进入| 91成人精品电影| 女人被躁到高潮嗷嗷叫费观| 成人国产一区最新在线观看 | 中国美女看黄片| 如日韩欧美国产精品一区二区三区| 制服诱惑二区| 久久精品久久久久久噜噜老黄| 少妇人妻久久综合中文| 少妇人妻久久综合中文| 91麻豆av在线| 成年女人毛片免费观看观看9 | 男人舔女人的私密视频| 在线精品无人区一区二区三| 国产精品成人在线| 久久国产亚洲av麻豆专区| 麻豆乱淫一区二区| 久久精品国产亚洲av涩爱| 后天国语完整版免费观看| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美精品综合一区二区三区| 午夜福利一区二区在线看| cao死你这个sao货| 伦理电影免费视频| 国产精品久久久久久精品古装| 亚洲,欧美精品.| 视频在线观看一区二区三区| 9色porny在线观看| 亚洲精品国产av蜜桃| 国产欧美亚洲国产| 免费少妇av软件| 久久鲁丝午夜福利片| 精品人妻熟女毛片av久久网站| 欧美日韩一级在线毛片| 日本av免费视频播放| 高清av免费在线| 七月丁香在线播放| 中文字幕最新亚洲高清| 国产高清videossex| 91老司机精品| xxx大片免费视频| 久久久国产一区二区| 久久ye,这里只有精品| 91老司机精品| 美国免费a级毛片| 大片免费播放器 马上看| 五月天丁香电影| 免费日韩欧美在线观看| 亚洲欧美激情在线| 久久久久久久大尺度免费视频| 免费看不卡的av| 亚洲七黄色美女视频| 中文字幕人妻熟女乱码| 女人爽到高潮嗷嗷叫在线视频| 黑人巨大精品欧美一区二区蜜桃| 免费女性裸体啪啪无遮挡网站| 一本—道久久a久久精品蜜桃钙片| 99久久99久久久精品蜜桃| 99久久人妻综合| 成在线人永久免费视频| 一级片'在线观看视频| 国产老妇伦熟女老妇高清| 天天躁狠狠躁夜夜躁狠狠躁| 51午夜福利影视在线观看| 中文精品一卡2卡3卡4更新| 丝袜脚勾引网站| 成人午夜精彩视频在线观看| 最黄视频免费看| 啦啦啦视频在线资源免费观看| 狠狠精品人妻久久久久久综合| 久久精品人人爽人人爽视色| 可以免费在线观看a视频的电影网站| 国语对白做爰xxxⅹ性视频网站| 无遮挡黄片免费观看| 免费女性裸体啪啪无遮挡网站| 午夜激情久久久久久久| 秋霞在线观看毛片| 女人久久www免费人成看片| 日本wwww免费看| netflix在线观看网站| av在线app专区| 黄频高清免费视频| 国产男女内射视频| 国产一级毛片在线| 性色av乱码一区二区三区2| 亚洲综合色网址| 亚洲欧美一区二区三区国产| 久久久欧美国产精品| 亚洲人成电影观看| 亚洲中文日韩欧美视频| 两个人免费观看高清视频| 丝袜人妻中文字幕| 真人做人爱边吃奶动态| 久久亚洲精品不卡| 婷婷色综合大香蕉| 久久精品亚洲熟妇少妇任你| 精品一品国产午夜福利视频| 天天操日日干夜夜撸| 丰满人妻熟妇乱又伦精品不卡| 国产日韩欧美视频二区| 90打野战视频偷拍视频| √禁漫天堂资源中文www| 在线观看www视频免费| 亚洲国产日韩一区二区| 精品熟女少妇八av免费久了| 久久ye,这里只有精品| 一级毛片我不卡| 2021少妇久久久久久久久久久| 精品福利观看| 男人添女人高潮全过程视频| 婷婷色麻豆天堂久久| 黄色怎么调成土黄色| 欧美日韩黄片免| av有码第一页| 国产精品免费大片| 中文字幕制服av| 最近中文字幕2019免费版| 777米奇影视久久| 亚洲欧美精品综合一区二区三区| 精品一品国产午夜福利视频| 亚洲第一青青草原| 日韩av免费高清视频| 一边亲一边摸免费视频| 日本av手机在线免费观看| av国产久精品久网站免费入址| 天堂俺去俺来也www色官网| 欧美+亚洲+日韩+国产| 国产日韩欧美在线精品| 精品人妻熟女毛片av久久网站| 纵有疾风起免费观看全集完整版| 操出白浆在线播放| 一级毛片电影观看| 在现免费观看毛片| 丝袜喷水一区| 日本av免费视频播放| 丰满少妇做爰视频| 尾随美女入室| 久久久久久亚洲精品国产蜜桃av| 一级黄色大片毛片| 久久久国产精品麻豆| 久久久久久亚洲精品国产蜜桃av| 欧美亚洲日本最大视频资源| 亚洲一码二码三码区别大吗| 91老司机精品| av欧美777| 国产黄频视频在线观看| 精品卡一卡二卡四卡免费| 久久精品熟女亚洲av麻豆精品| 国产一区二区三区av在线| 精品卡一卡二卡四卡免费| 日韩精品免费视频一区二区三区| 国产国语露脸激情在线看| 国产片特级美女逼逼视频| 男女床上黄色一级片免费看| 久久久久视频综合| 久久人人爽人人片av| 国产三级黄色录像| 亚洲国产中文字幕在线视频| 欧美成人精品欧美一级黄| 久久久国产一区二区| 又紧又爽又黄一区二区| 午夜激情av网站| 青春草视频在线免费观看| 99精国产麻豆久久婷婷| 欧美成狂野欧美在线观看| 男女边吃奶边做爰视频| 日韩免费高清中文字幕av| 亚洲av成人不卡在线观看播放网 | 国产日韩欧美在线精品| 成人免费观看视频高清| 久久99精品国语久久久| 熟女少妇亚洲综合色aaa.| 久久精品亚洲熟妇少妇任你| 日韩 欧美 亚洲 中文字幕| 日本av免费视频播放| 中文字幕精品免费在线观看视频| 黄色视频不卡| 精品欧美一区二区三区在线| 亚洲精品国产av蜜桃| 男的添女的下面高潮视频| 亚洲精品久久午夜乱码| 欧美黑人精品巨大| 亚洲国产av影院在线观看| 久久久久国产精品人妻一区二区| 亚洲av电影在线进入| 国产视频首页在线观看| av线在线观看网站| 日本av免费视频播放| 日本欧美国产在线视频| 免费高清在线观看日韩| 欧美日韩视频高清一区二区三区二| 国产一区有黄有色的免费视频| 在线观看免费日韩欧美大片| 国产片特级美女逼逼视频| 九草在线视频观看| 久久午夜综合久久蜜桃| 久久久国产欧美日韩av| 久久久精品国产亚洲av高清涩受| 秋霞在线观看毛片| 人成视频在线观看免费观看| 夜夜骑夜夜射夜夜干| 在线精品无人区一区二区三| 国产高清videossex| 精品人妻一区二区三区麻豆| 操出白浆在线播放| 一边摸一边做爽爽视频免费| 亚洲av日韩在线播放| 高清不卡的av网站| 久久国产精品男人的天堂亚洲| 国产一区二区三区av在线| 久久久久网色| 捣出白浆h1v1| 欧美另类一区| 美女中出高潮动态图| 中文字幕人妻熟女乱码| 性色av一级| 国产亚洲精品久久久久5区| 青青草视频在线视频观看| 午夜影院在线不卡| 女人精品久久久久毛片| 亚洲 国产 在线| 亚洲激情五月婷婷啪啪| 成年人午夜在线观看视频| 国产欧美亚洲国产| 欧美日韩视频高清一区二区三区二| 美女中出高潮动态图| 国产成人av教育| 一级a爱视频在线免费观看| 久久亚洲国产成人精品v| 美女主播在线视频| 人人妻人人添人人爽欧美一区卜| 亚洲午夜精品一区,二区,三区| 男女无遮挡免费网站观看| 久久久亚洲精品成人影院| 国产高清videossex| 亚洲熟女毛片儿| 国产精品久久久人人做人人爽| 一级毛片黄色毛片免费观看视频| 日韩一卡2卡3卡4卡2021年| 国产免费视频播放在线视频| 男女午夜视频在线观看| 美女脱内裤让男人舔精品视频| 黄色片一级片一级黄色片| 国产免费一区二区三区四区乱码| 亚洲av在线观看美女高潮| 久久人妻福利社区极品人妻图片 | 精品亚洲乱码少妇综合久久| 视频区欧美日本亚洲| 亚洲伊人久久精品综合| 色婷婷久久久亚洲欧美| 亚洲国产精品国产精品| 国产片内射在线| 日韩 亚洲 欧美在线| avwww免费| 丝袜美腿诱惑在线| 色视频在线一区二区三区| 国产真人三级小视频在线观看| 国产极品粉嫩免费观看在线| 久久久久久久大尺度免费视频| 国产成人影院久久av| 男人爽女人下面视频在线观看| 夫妻性生交免费视频一级片| 啦啦啦视频在线资源免费观看| 亚洲国产精品999| 国产精品三级大全| 久久精品aⅴ一区二区三区四区| 日韩一卡2卡3卡4卡2021年| 欧美成人精品欧美一级黄| 国产黄色免费在线视频| 亚洲欧美日韩另类电影网站| 亚洲精品自拍成人| 97人妻天天添夜夜摸| 中文字幕人妻丝袜制服| 一本综合久久免费| 色网站视频免费| 黄色毛片三级朝国网站| 精品少妇黑人巨大在线播放| 欧美大码av| 欧美av亚洲av综合av国产av| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久精品古装| 男女边摸边吃奶| 亚洲精品国产av成人精品| 中文字幕人妻丝袜一区二区| 免费av中文字幕在线| 国产一区有黄有色的免费视频| 精品一区二区三区四区五区乱码 | 在线天堂中文资源库| 高清av免费在线| 国产精品麻豆人妻色哟哟久久| 亚洲人成电影免费在线| 天天躁夜夜躁狠狠躁躁| 亚洲欧美成人综合另类久久久| 久久精品久久久久久久性| 视频在线观看一区二区三区| 一级黄色大片毛片| 丝袜人妻中文字幕| 韩国高清视频一区二区三区| 国产黄色免费在线视频| 麻豆乱淫一区二区| 久久九九热精品免费| 午夜两性在线视频| 嫁个100分男人电影在线观看 | 日本av免费视频播放| 日韩一区二区三区影片| 久久性视频一级片| 天堂俺去俺来也www色官网| 中国美女看黄片| 水蜜桃什么品种好| 国产av一区二区精品久久| 别揉我奶头~嗯~啊~动态视频 | 亚洲男人天堂网一区| 中文字幕最新亚洲高清| 亚洲国产精品999| av天堂久久9| 国产不卡av网站在线观看| √禁漫天堂资源中文www| avwww免费| 大片电影免费在线观看免费| 美女午夜性视频免费| 可以免费在线观看a视频的电影网站| 男女午夜视频在线观看| 最新在线观看一区二区三区 | 成人手机av| 亚洲熟女毛片儿| 久久久亚洲精品成人影院| 欧美精品亚洲一区二区| cao死你这个sao货| 侵犯人妻中文字幕一二三四区| 天天操日日干夜夜撸| 中文乱码字字幕精品一区二区三区| 亚洲欧美精品自产自拍| 亚洲欧美中文字幕日韩二区| 51午夜福利影视在线观看| av电影中文网址| 天天添夜夜摸| 国产精品久久久久成人av| 美女脱内裤让男人舔精品视频| 99久久综合免费| 成人免费观看视频高清| 9191精品国产免费久久| 中文字幕另类日韩欧美亚洲嫩草| 亚洲自偷自拍图片 自拍| 亚洲av电影在线进入| 99久久99久久久精品蜜桃| 国产男女超爽视频在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲九九香蕉| 女性被躁到高潮视频| 精品国产乱码久久久久久男人| 亚洲国产欧美网| 高清视频免费观看一区二区| 91老司机精品| 国产亚洲欧美在线一区二区| 丰满少妇做爰视频| 国产精品二区激情视频| 91麻豆av在线| 女人高潮潮喷娇喘18禁视频| 看十八女毛片水多多多| 夫妻午夜视频| av线在线观看网站| 久久久亚洲精品成人影院| 国产精品一区二区在线不卡| 日本欧美国产在线视频| 女人爽到高潮嗷嗷叫在线视频| av在线老鸭窝| 纵有疾风起免费观看全集完整版| 久久 成人 亚洲| 午夜免费观看性视频| 日本vs欧美在线观看视频| 建设人人有责人人尽责人人享有的| 国产熟女欧美一区二区| 国产成人av激情在线播放| 国产又爽黄色视频| 操出白浆在线播放| 秋霞在线观看毛片| 亚洲欧美精品自产自拍| 美女福利国产在线| 免费看不卡的av| 国产97色在线日韩免费| 悠悠久久av| 久久天躁狠狠躁夜夜2o2o | 男女无遮挡免费网站观看| 亚洲精品国产一区二区精华液| 一区二区日韩欧美中文字幕| 青春草亚洲视频在线观看| 午夜福利,免费看| 校园人妻丝袜中文字幕| 天天添夜夜摸| 曰老女人黄片| 国产精品国产三级专区第一集| 热re99久久精品国产66热6| 精品熟女少妇八av免费久了| 国产亚洲精品第一综合不卡| 午夜久久久在线观看| 蜜桃在线观看..| 久久人人爽av亚洲精品天堂| 亚洲中文av在线| 亚洲av片天天在线观看| 男女之事视频高清在线观看 | 性高湖久久久久久久久免费观看| 高清欧美精品videossex| 亚洲黑人精品在线| 国产日韩欧美视频二区| 日本一区二区免费在线视频| 亚洲成人免费av在线播放| 日韩中文字幕视频在线看片| 一区二区日韩欧美中文字幕| 亚洲欧美日韩高清在线视频 | 在线精品无人区一区二区三| 王馨瑶露胸无遮挡在线观看| 大话2 男鬼变身卡| 欧美精品亚洲一区二区| 最黄视频免费看| 亚洲av成人不卡在线观看播放网 | 啦啦啦在线免费观看视频4| 久久精品国产a三级三级三级| 777久久人妻少妇嫩草av网站| 欧美日本中文国产一区发布| 欧美日韩亚洲国产一区二区在线观看 | 一区在线观看完整版| 在现免费观看毛片| 国产一区二区三区综合在线观看| 首页视频小说图片口味搜索 | 免费观看a级毛片全部| 老司机亚洲免费影院| 亚洲精品美女久久久久99蜜臀 | 乱人伦中国视频| 久久天躁狠狠躁夜夜2o2o | 性色av乱码一区二区三区2| 欧美久久黑人一区二区| 女人高潮潮喷娇喘18禁视频| 国产一区二区三区av在线| 99热全是精品| 国产一区二区激情短视频 | 亚洲av美国av| 最新在线观看一区二区三区 | 亚洲免费av在线视频| 国产免费又黄又爽又色| 久久人人爽av亚洲精品天堂| 天天添夜夜摸| 久热爱精品视频在线9| 成人黄色视频免费在线看| 嫩草影视91久久| 各种免费的搞黄视频| 国产又色又爽无遮挡免| 黄色毛片三级朝国网站| 满18在线观看网站| 韩国高清视频一区二区三区| 日本vs欧美在线观看视频| 一区福利在线观看| 青青草视频在线视频观看| 国产男女超爽视频在线观看| 日韩 亚洲 欧美在线| 别揉我奶头~嗯~啊~动态视频 | 国产女主播在线喷水免费视频网站| 建设人人有责人人尽责人人享有的| 国产成人av教育| 国产淫语在线视频| 亚洲 欧美一区二区三区| 狠狠婷婷综合久久久久久88av| 久久久久国产精品人妻一区二区| 18禁黄网站禁片午夜丰满| 波多野结衣av一区二区av| 日本猛色少妇xxxxx猛交久久| 国产亚洲一区二区精品| 69精品国产乱码久久久| 天天影视国产精品| 91老司机精品| 国产伦人伦偷精品视频| 性色av一级| 国产成人免费观看mmmm| 亚洲精品一二三| 99国产精品99久久久久| 老鸭窝网址在线观看| 国产成人欧美| 国产男人的电影天堂91| 亚洲精品美女久久av网站| 免费在线观看视频国产中文字幕亚洲 | 国产成人影院久久av| 韩国精品一区二区三区| 久久久精品94久久精品| 人人妻人人澡人人爽人人夜夜| 中文乱码字字幕精品一区二区三区| 成人手机av| 日本91视频免费播放| 中国国产av一级| 视频区欧美日本亚洲| 777久久人妻少妇嫩草av网站| 99香蕉大伊视频| 美女午夜性视频免费| xxxhd国产人妻xxx| 久久99一区二区三区| 妹子高潮喷水视频| 一区二区av电影网| 另类亚洲欧美激情| 9191精品国产免费久久| 汤姆久久久久久久影院中文字幕| 女人高潮潮喷娇喘18禁视频| 中国国产av一级| 美女脱内裤让男人舔精品视频| 日韩一区二区三区影片| 最新在线观看一区二区三区 | 男女边吃奶边做爰视频| 久久久国产精品麻豆| 一二三四在线观看免费中文在| 少妇人妻 视频| 午夜福利一区二区在线看| 亚洲成人手机| 日韩大码丰满熟妇| 国产亚洲av片在线观看秒播厂| 两个人看的免费小视频| 免费观看av网站的网址| 老司机靠b影院| 可以免费在线观看a视频的电影网站| 国产成人一区二区在线| 国产在视频线精品| 国产黄频视频在线观看| 国产在线观看jvid| 黄色怎么调成土黄色| 大型av网站在线播放| 中文字幕制服av| 国产精品 欧美亚洲| 亚洲av电影在线进入| 亚洲欧美激情在线| 少妇的丰满在线观看| 国产成人一区二区在线| 真人做人爱边吃奶动态| 夫妻性生交免费视频一级片| 欧美精品啪啪一区二区三区 | 亚洲欧洲精品一区二区精品久久久| 国产亚洲av片在线观看秒播厂| 人人妻人人澡人人看| 99久久人妻综合| 国产男女内射视频| www.999成人在线观看| 99九九在线精品视频| 人人妻人人添人人爽欧美一区卜| 好男人电影高清在线观看| 美女主播在线视频| 精品亚洲乱码少妇综合久久| 国产成人欧美在线观看 | 一级,二级,三级黄色视频| 丝袜在线中文字幕| 18禁裸乳无遮挡动漫免费视频| 久久精品熟女亚洲av麻豆精品| 天堂8中文在线网| 美女高潮到喷水免费观看| 我要看黄色一级片免费的| 日韩一本色道免费dvd| 晚上一个人看的免费电影| 国产国语露脸激情在线看| 丁香六月天网| 男女国产视频网站| 欧美黑人欧美精品刺激| 午夜福利视频精品| 亚洲九九香蕉| 亚洲五月色婷婷综合| 午夜老司机福利片| 美女中出高潮动态图| 成人国语在线视频| kizo精华| 9191精品国产免费久久| 欧美在线一区亚洲| 飞空精品影院首页| 中文字幕色久视频| 久久精品aⅴ一区二区三区四区| 久久久久视频综合| 人妻 亚洲 视频| 亚洲,一卡二卡三卡| 90打野战视频偷拍视频| 女人爽到高潮嗷嗷叫在线视频| 美女大奶头黄色视频| 亚洲精品中文字幕在线视频| 欧美成狂野欧美在线观看| 久久人人97超碰香蕉20202| 99热网站在线观看| 考比视频在线观看| 欧美大码av| www.999成人在线观看| 久久久久久久久久久久大奶| 在线 av 中文字幕| svipshipincom国产片| 亚洲成av片中文字幕在线观看| 9191精品国产免费久久| 精品福利永久在线观看| 久久精品久久久久久久性|