Abstract
Topological photonic crystals, which offer topologically protected and back-scattering-immune transport channels, have recently gained significant attention for both scientific and practical reasons. Although most current studies focus on dielectric materials with weak dispersions, this study focuses on topological phases in dispersive materials and presents a numerical study of Chern insulators in gaseous-phase plasma cylinder cells. We develop a numerical framework to address the complex material dispersion arising from the plasma medium and external magnetic fields and identify Chern insulator phases that are experimentally achievable. Using this numerical tool, we also explain the flat bands commonly observed in periodic plasmonic structures, via local resonances, and how edge states change as the edge termination is periodically modified. This work opens up opportunities for exploring band topology in new materials with non-trivial dispersions and has potential radio frequency (RF) applications, ranging from plasma-based lighting to plasma propulsion engines.
Similar content being viewed by others
Introduction
The discovery that topological phases can exist beyond electronics1,2,3 has garnered significant attention in other wave systems such as photonics4,5,6, plasmonics7,8,9, polaritonics10, acoustics11,12, and even water waves13,14. Of particular interest is the quantum anomalous Hall effect, also known as Chern insulators, which can offer topologically protected transport channels and immunity against back-scattering at the interface with normal insulators. Such transport channels have both scientific inquiry and practical applications in fields such as optical communication, low-loss waveguides, and circulators.
To create photonic Chern insulators, time-reversal symmetry must be broken while retaining near-Hermiticity. This requires a permittivity or permeability tensor that breaks reciprocity, such that \(\varvec{\varepsilon }^{\text {T}}\ne \varvec{\varepsilon }\) or \(\varvec{\mu }^{\text {T}}\ne \varvec{\mu }\). The properties of the underlying materials, along with their geometric aspects, determine the topological invariants (Chern numbers) of the electromagnetic band gaps and their transport properties.
Most current studies of photonic Chern insulators have relied on gyromagnetic materials (e.g. yttrium iron garnet2,15,16,17) and external magnetic fields. These materials are mostly dielectric in nature, meaning their permeability remains positive in the frequency range of interest, resulting in Bloch modes that are delocalized in the photonic crystals. There has been recent interest in gyroelectric materials, which have the potential for large Faraday effects in magnetized plasmas in metals. This focus has centered on continua18 and plasmonic crystals7,8,9, where the Chern insulator phase is mostly composed of coupled plasmonic resonances localized at individual sites with Drude-like material responses and negative permittivity.
Here, we present a numerical study of photonic Chern insulators in plasma crystals with gaseous-phase plasma, which simultaneously exhibit both extended photonic bands and localized plasmonic modes in the RF regime. Our design is based on a 2D crystal of plasma cells placed in an external magnetic field. Without the magnetic field, the responses of the plasma elements are Drude-like, and the associated structure is known to support coexisting de-localized and localized modes in suitable polarizations19,20,21. We explore the time-reversal broken generalization, which exhibits an interesting interplay between the Drude and Lorentzian dispersion due to the applied magnetic field causing cyclotron motions in the plasma. We propose a plasma crystal design that features a Chern insulator gap between de-localized photonic bands, coexisting with nearby dense groups of flat bands associated with localized plasmons of fixed-handedness. On termination of this crystal, we observe a rich interplay between localized plasmonic bands and de-localized chiral edge states inside the Chern insulating gap. Finally, we explore how local and de-localized edge modes evolve under continuous deformations of the interface between the Chern insulator and perfect magnetic conductors. The mode evolution can then be interpreted as a manifestation of the filling anomalies.
Our work is organized as follows: first, we review the plasma dispersion without and with an external magnetic field. Then, we present the band structure of plasma crystals without and with an external magnetic field and their associated topological invariants. Next, we explain the origin of the observed flat bands in calculations as localized surface plasmon polariton resonances. Afterward, we explore how the chiral edge state dispersion evolves when the edge termination changes. Finally, we discuss the limitations existing in our calculations and practical aspects related to the experimental verification of our proposal.
Results
Dispersion of plasma without magnetic field: drude model
We start by presenting a comprehensive review of the dispersion properties of the plasma permittivity (\(\varepsilon\)) by employing the widely accepted Drude model22. In our analysis, we made the fundamental assumption that the positive ions within the plasma are significantly heavier and thus remain immobile. Consequently, the contribution to the volume current density (J) is exclusively attributed to the movement of electrons: \({\textbf {J}} = -n e {\textbf {v}}\). In this equation, n represents the volume density of electrons, and \(-e\) corresponds to the electron charge of an individual electron. The equation of motion for electrons in the plasma crystals reads:
Here, \(m_{\textrm{e}}\) is the electron mass, \(\gamma\) is the damping rate, and \({\textbf {E}}\) is the electric field. For harmonic solutions at a fixed angular frequency \(\omega\), all temporal derivatives can be substituted via \(\partial _t \rightarrow -\textrm{i}\omega\). Accordingly, the frequency-dependent plasma conductivity \(\sigma\) can be written as
Here, \(\omega _{\textrm{p}} = \sqrt{\frac{ne^2}{m_{\textrm{e}}\varepsilon _0}}\) is the plasma frequency. Noting that the volume current density is also related to the electric polarization: \({\textbf {J}} = \partial _t {\textbf {P}} = -\textrm{i}\omega {\textbf {P}}\), the Drude permittivity of plasma \(\varepsilon\) can be defined as:
Dispersion of magnetized plasma
Following these steps, we calculated the permittivity tensor describing gaseous phase plasma placed in an external magnetic field (Fig. 1a). Following Eq. (3), we need to re-write the conductivity tensor (\(\bar{\bar{\sigma }}\)), based on the updated equation of motion for electrons. Considering the Lorentz force, the equation of motion becomes:
Re-writing the equation in circular bases in the xy plane: \((J_+, J_-, J_z) = (\frac{J_x+\textrm{i}J_y}{\sqrt{2}}, \frac{J_x-\textrm{i}J_y}{\sqrt{2}}, J_z)\) and \((E_+, E_-, E_z) = (\frac{E_x+\textrm{i}E_y}{\sqrt{2}}, \frac{E_x-\textrm{i}E_y}{\sqrt{2}}, E_z)\), both matrices \(\bar{\bar{\sigma }}\) and \(\bar{\bar{\varepsilon }}\) become diagonal. Specifically, the conductivity tensor \(\bar{\bar{\sigma }}\) can be expressed as:
Here \(\omega _{\textrm{c}} = \frac{eB}{m_{\textrm{e}}}\) is the cyclotron resonance. Accordingly, the permittivity tensor of the magnetized plasma \(\bar{\bar{\varepsilon }}\) has the following form in a Cartesian basis:
where the column vectors of the matrix U label the directions of the optical principle axes:
Meanwhile, along the principal axes, the material dispersion can be expressed as:
This outcome can be intuitively understood as follows: when subjected to an external magnetic field along the \(\hat{\textbf{z}}\) direction, electrons undergo cyclotron motion, forming orbits at a constant angular frequency of \(\omega _{\textrm{c}}\) within the xy plane. In a co-rotating (counter-rotating) reference frame synchronized with the electrons, the induced electric field preserves the same (opposite) circular polarization but is slightly shifted in frequency by \(\omega _{\textrm{c}}\) (\(-\omega _{\textrm{c}}\)), yielding the distinctive properties \(\varepsilon _+\) and \(\varepsilon _-\), respectively. On the other hand, the electric field applied in the z direction remains unaffected by the cyclotron motion. Consequently, the dispersion of \(\varepsilon _z\) maintains the conventional characteristics of the standard Drude dispersion.
Next, we numerically compute the material dispersion as functions of frequency \(f = \omega /2\pi\) using the typical values in gaseous phase plasma. The results are shown in Fig. 1, where the external magnetic field is set at 0.054 T and the corresponding cyclotron resonance is at \(f_{\textrm{c}} = \omega _{\textrm{c}}/2\pi = 1.5\,\text {GHz}\). The plasma frequency is controlled by the number density of electrons, which is set to be a practical value of \(n=3.1\times 10^{11}\,\text {cm}^{-3}\) throughout the calculations. Accordingly, the plasma frequency is at \(f_{\textrm{p}} = \omega _{\textrm{p}}/2\pi = 5\,\text {GHz}\). The real and imaginary parts of the permittivity are shown in blue and red, respectively. The parameter \(\gamma\) refers to the damping rate, which mostly originates from the electron-ion collisions at room temperature. For helium plasma gas, the damping rate is roughly linearly proportional to the gas pressure: \(\gamma = 0.318\,\text {GHz\,Torr}^{-1}\times p\). Using a practical pressure of \(p= 0.314\,\text {Torr}\), the damping is set to be \(\gamma = 0.1\,\text {GHz}\) throughout our calculations.
As expected from Eq. (8), both \(\varepsilon _+\) and \(\varepsilon _-\) are affected by the cyclotron resonance and deviate from the standard Drude model (\(\varepsilon _z\)). For example, at low frequencies (\(f \approx 0\)), \(\varepsilon _+\) is positive and diverges as 1/f; \(\varepsilon _{-}\) is negative and diverges as 1/f; meanwhile, the Drude model \(\varepsilon _z\) is negative and diverges as \(1/f^2\). These fundamental differences in scaling lead to challenges when fitting the dispersion to standard formalism in commercial software, as described later.
We note that our description of the plasma medium is limited by a few approximations. Overcoming these approximations would lead to modifications of our results and will be discussed elsewhere. First, only electrons are assumed to move under the influence of external electromagnetic fields, while ions are assumed to be always stationary. Second, our permittivity ignores non-local effects, which leads to a frequency gap between surface plasmon polaritons traveling in opposite directions23,24,25. Finally, we neglect the inhomogeneous broadening effect, and thus the bandwidth of our permittivity is defined solely by electron damping.
Band structures of plasma crystals without external magnetic field
Utilizing the dispersion equation outlined in Eq. (8)26, we compute the band structure associated with a square lattice comprising plasma cylinders positioned within an air medium. Our investigation primarily centered around quadratic point degeneracy, a distinctive phenomenon safeguarded by both spatial symmetry and time-reversal symmetry. This quadratic point is important in our study, as it arises from the combined effects of spatial symmetry (\(C_4\)) and time-reversal symmetry (T), which is broken when T is broken.
The plasma photonic crystal unit cell is shown in Fig. 2a, where the lattice constant a is 6 cm, and the radius of the cylinder r is 1.5 cm. The external magnetic field along the z direction preserves the mirror symmetry in z (\(\sigma _z\)) and separates the electromagnetic modes into two mode types: to avoid possible terminology confusion, we follow the definition of Sakoda27 to distinguish the E-polarization case (with the \(\textbf{E}\)-field parallel to the z-axis) from the H-polarization case (with the \(\textbf{H}\)-field parallel to the z-axis). We also use the E-case and H-case for short. Given that the gyroelectric response has response components only in the xy-plane but not along the z-direction [Eq. (5)], we focus solely on the H-case and disregard the associated E-case.
Without breaking time-reversal symmetry (T), i.e., \(B = 0\), the plasma dispersion follows the standard Drude model, and the H-case can be calculated using a standard finite element method (FEM). The results are plotted along high-symmetry lines in the Brillouin zone (Fig. 2b). As shown, a pair of quadratic degeneracies are found at the M point in the Brillouin zone around 2.4 GHz, which is protected by the 90-degree rotation symmetry \(C_{4z}\) and T. Specifically, the two modes, marked as ‘\(+\)’ and ‘−’, have \(C_{4z}\) indices of \(\pm i\),. They are connected to each other by T. The phases of the corresponding mode profiles, \(\mathop {\textrm{arg}}(H_z)\), confirm the \(C_{4z}\) indices of the two modes. We note that a set of flat bands is observed in the calculation (blue ribbon) with an upper-frequency bound of \(f_{\textrm{p}}/\sqrt{2} = 3.5\,\text {GHz}\), which is further discussed in the next section.
Band structure of the magnetized-plasma photonic crystals and Chern insulators
In the presence of an external magnetic field, the behavior of the plasma dispersion (\(\varepsilon _{+}\) and \(\varepsilon _{-}\)) deviates from the conventional Drude model. This alteration introduces complexities that pose challenges for accurately representing the band structure using standard material dispersion models available in commercial software solutions. As a result, we modify a standard Finite Element Method (FEM) technique to effectively capture the intricacies of the band structure, especially concerning band crossings that demand precise resolution. Following conventional methodologies (e.g., in28,29,30,31,32), we applied Floquet periodic boundary conditions in the 2D spatial domain. This technique allowed us to address a quadratic eigenvalue problem within a square unit cell configuration, as illustrated in Fig. 2a. To solve this problem, we discretized the unit cell to formulate a weak FEM approach. Our computational process leveraged COMSOL Multiphysics™, a versatile software platform adept at constructing customized systems of coupled equations. For the numerical eigenvalue computations, we integrated the Portable Large Scale Eigenvalue Package (P_ARPACK)33. P_ARPACK, which is a parallel implementation derived from the ARPACK software34. P_ARPACK employs the Implicitly Restarted Arnoldi Method (IRAM) and is particularly suited for tackling large sparse eigenvalue problems, even when constrained by specific eigenvalue counts. The combination of COMSOL Multiphysics™ and P_ARPACK facilitated the development of a scalable and efficient eigenvalue solver, making it feasible to tackle the customized quadratic eigenvalue problems arising from our research. The adaptability of the software, coupled with its user-friendly interface, allowed us to fine-tune key parameters of the IRAM-based solver according to our requirements. Furthermore, the efficiency of employing P_ARPACK within the COMSOL framework for addressing customized quadratic eigenvalue problems has been demonstrated by multiple studies, exemplified by references31,35,36. This substantiates the reliability and advantages of our chosen approach for solving complex eigenvalue problems in our research context.
In contrast to the known methods largely employing the auxiliary equations for the polarization vector, we couple the weak-form equations for the current density and the \(\textbf{E}\)-field, achieving a stable performance with flexible tracking of a desired number of bands and unambiguous resolution of the band crossings. We have provided a general equation for current density under an external magnetic field. In the H-case, we directly use the External Current Density Interface with an auxiliary algebraic equation (AE) for in-plane vectors,
We then employ a normalized current density \(\textbf{j}_{xy} = \textbf{J}_{xy}/\left( \varepsilon _0\omega _{\textrm{p}}^2\right)\) in a normalized Drude model
The weak form is obtained by integrating the dot product of Eq. (10) with an arbitrary test function \(\textbf{j}_{xy} = \mathop {\textrm{test}}(\textbf{j}_{xy})\) over the Drude material domain,
An auxiliary equation for the current density is introduced to the COMSOL framework through the Weak Contribution interface,
where \(\textbf{E}_{xy} = \mathop {\textrm{test}}(\textbf{E}_{xy})\) is the \(\textbf{E}\)-field test function and \(k_{\textrm{p}}=\omega _{\textrm{p}}/{c}\) is the plasma wave number. The approach exhibits good error convergence with accuracy controlled through the meshing density and the FE order. In contrast with31,35,36,37,38 we completely exclude the polarization vector, reducing the order of the auxiliary equations and improving the numerical accuracy at \(\omega \rightarrow 0\). The full details on the numerical implementation and verification of the IRAM-based eigensolver are not the main focus of the paper and these details will be published elsewhere.
An example of the calculated band structures is shown in Fig. 3a when the external magnetic field is set to be \(B = 0.054\,\text {T}\). As time-reversal symmetry T is broken, the M-point degeneracy is lifted, opening an 8% full energy gap, from 2.18 to 2.37 GHz (green ribbon). As the structure still maintains \(C_{4z}\) symmetry, the Chern number C of the first band is necessarily non-trivial39, since
where \(C_{2z}\) is the phase change of the mode profile after a \(\pi\)-rotation. \(C_{2z} (\Gamma )=1\) since the phase of an electromagnetic wave is locked at zero frequency, and \(C_{2z} (M)=-1\) because it originates from the time-reversal symmetry breaking of the degenerate modes shown in Fig. 2b. As a result, the first gap highlighted in green corresponds to a Chern insulator40 and supports unidirectional transport channels, as shown next. Due to the magnetic field, the flat bands split into two regions (two blue ribbons), as explained in detail in the next section.
Flat bands from localized surface plasmon polariton resonances
In this section, we elucidate the origin of the flat bands observed above, from the viewpoint of localized surface plasmon polaritonic (SPP) resonances. We note that such flat band features are also commonly observed elsewhere, such as in metallic photonic crystals28. The one unusual feature is related to the splitting of flat band regions under an external magnetic field (Fig. 3a).
It is more straightforward to understand the flat bands if we consider the local SPP resonances supported by a single plasma cylinder7,8. As the cylinder has full rotation symmetry, the SPP resonances can be labeled by different azimuthal numbers, m, corresponding to different angular momenta. A few interesting features can be observed in Fig. 3b. First, the resonances with negative m (orange dots, rotating counter-clockwise) are at higher frequencies than the resonances with positive m (blue dots, clockwise). Second, the local SPP resonance frequency generally increases with |m|, approaching different bounds near the two ends: \((\sqrt{f_{\textrm{c}}^2+2f_{\textrm{p}}^2}+f_{\textrm{c}})/2=4.4\,\text {GHz}\) when m approaches \(-\infty\) and \((\sqrt{f_{\textrm{c}}^2+2f_{\textrm{p}}^2}-f_{\textrm{c}})/2=2.9\,\text {GHz}\) when m approaches \(+\infty\). Finally, the resonances with larger |m|s are better localized in space than the resonances with smaller |m|. Such a trend can be well observed in the comparison between the mode profile of \(m=-5\) (more localized) versus the model profile of \(m=1\) (more extended, inset of Fig. 3b). Taken together, at large |m|, the SPP resonances are tightly confined to individual cylinders, resulting in minimal modal overlaps. This limited overlap leads to reduced dispersion, which, in turn, results in the formation of flat bands. Thus, each flat band corresponds to a different m. Furthermore, the frequency of the flat bands (blue ribbons) also splits into two regions, approaching the two frequency bounds mentioned above.
Here we note that the number of calculated flat bands using our numerical method increases with increased mesh density—a common feature also observed in the literature28—although the frequency of the flat bands is always confined to the blue regions. Furthermore, the reliable azimuthal number of calculated modes is always limited by the numerical resolution. In our specific setting, modes with \(|m|>10\) are no longer reliable. See Supplementary Information for more details.
Finally, we verify the local SPP resonance frequencies in each plasma cylinder using the SPP dispersion along the interface between air and magnetized plasma, where a good agreement is found (Fig. 3c). The momentum of each local SPP resonance is determined as m/r, where m is the azimuthal number and r is the radius of the cylinder.
Evolution of the chiral edge state dispersion with changing edge termination
While the existence of chiral edge states (CES) is guaranteed at the interface between a Chern insulator and a trivial insulator, their exact dispersion depends on the details of the interface. In this section, we continuously change the interface configuration and study how the CES dispersion evolves accordingly. Our finding suggests that CES dispersion essentially reflects that a localized plasma resonance emerges at the interface when the plasma is cut through, and the frequency of the antenna state decreases with the shrinking of the plasma region.
The dispersion of Chiral Edge States (CES) is computed through interactions at the interfaces between a Chern insulator super-cell and a pair of perfect magnetic conductors (PMC), as visualized in Fig. 4. Each unit cell within the Chern insulator preserves the identical design and parameters presented in Fig. 3. The arrangement entails a fixed lower interface (blue), while the upper interface (red) is systematically adjusted by incrementing the distance d from 0 to a. Beginning with \(d/a=0\), a significant outcome is observed in the super-cell dispersion-two CES bands materialize. One is located at the upper interface (red), and the other emerges at the lower interface (blue). As the ratio d/a increases, the CES located at the upper interface (red) progresses in frequency, ascending until the PMC interacts with the plasma cell at \(d/a = 0.25\). With further increments in d/a, the high-frequency section of the CES remains relatively unchanged while the low-frequency component continues ascending. This trend persists until \(d/a=0.5\) where the PMC bisects the plasma cell. As d/a continues to rise, the CES dispersion transforms into a progressively flatter profile.
At \(d/a{=}0.7\), the CES dispersion undergoes a pivotal change-dividing into two bands. One of these bands traverses the topological gap and retains CES characteristics, while the other is a trivial band existing outside the topological gap. A second splitting transpires at \(d/a=0.735\), yielding an additional trivial band. As the distance d/a increases further, both trivial bands descend in frequency, crossing through the initial bulk continuum and eventually vanishing at zero frequency \(f=0\) when \(d/a=0.75\). See Supplementary Information for the case when \(d/a=0.74996\). Upon reaching \(d/a{=}1\), the interface configuration regresses to the \(d/a=0\) arrangement, leading the CES dispersion to revert to its initial configuration.
In essence, this detailed description captures the intricate evolution of CES dispersions as the distance d is varied in relation to the unit cell size a. The interplay between interface configurations and distance ratios results in a rich spectrum of phenomena and transitions.
Discussion
Our simulation has a few limitations that should be noted. Firstly, the ions are assumed to be stationary and not moving, i.e., only electrons are allowed to move under the external fields. Introducing ions’ motion would lead to an effective mass and lower plasma frequency. Secondly, our permittivity model ignores non-local effects, which results in a frequency gap between surface plasmon polaritons traveling in opposite directions23,24. Thirdly, we have neglected the inhomogeneous broadening effect in the plasma, and the bandwidth of our permittivity is contributed solely by electron damping. Introducing the inhomogeneous broadening effect would further broaden the energy bands and reduce the effective size of the band gaps. It is also worth noting that the cyclotron frequency of the plasma at 1.5 GHz is much larger than the electron damping rate at 0.1 GHz, and hence the homogeneous broadening effect contributes solely to the broadening of the energy bands without qualitatively changing the band topology.
Despite these limitations, we believe that our proposal is feasible to demonstrate in an experiment. The required magnetic field of 0.314 T can be achieved, even over large areas, via commercial electromagnets or permanent magnets. Meanwhile, the required carrier number density of \(3.1\times 10^{11}\,\text {cm}^{-3}\) and the pressure of 0.314 Torr are also within the typical range in experiments.
Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
References
Haldane, F. D. M. & Raghu, S. Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry. Phys. Rev. Lett. 100, 013904. https://doi.org/10.1103/PhysRevLett.100.013904 (2008).
Wang, Z., Chong, Y., Joannopoulos, J. D. & Soljačić, M. Observation of unidirectional backscattering-immune topological electromagnetic states. Nature 461, 772–775. https://doi.org/10.1038/nature08293 (2009).
Lu, L., Joannopoulos, J. D. & Soljačić, M. Topological photonics. Nat. Photonics 8, 821–829. https://doi.org/10.1038/nphoton.2014.248 (2014).
Sigalas, M. M., Chan, C. T., Ho, K. M. & Soukoulis, C. M. Metallic photonic band-gap materials. Phys. Rev. B 52, 11744–11751. https://doi.org/10.1103/PhysRevB.52.11744 (1995).
He, L., Addison, Z., Mele, E. J. & Zhen, B. Quadrupole topological photonic crystals.. Nat. Commun. 11, 3119. https://doi.org/10.1038/s41467-020-16916-z (2020).
Lu, J., He, L., Addison, Z., Mele, E. J. & Zhen, B. Floquet topological phases in one-dimensional nonlinear photonic crystals. Phys. Rev. Lett. 126, 113901. https://doi.org/10.1103/PhysRevLett.126.113901 (2021).
Jin, D. et al. Topological magnetoplasmon.. Nat. Commun. 7, 13486. https://doi.org/10.1038/ncomms13486 (2016).
Jin, D. et al. Infrared topological plasmons in graphene. Phys. Rev. Lett. 118, 245301. https://doi.org/10.1103/PhysRevLett.118.245301 (2017).
Pan, D., Yu, R., Xu, H. & de Abajo, F. J. G. Topologically protected Dirac plasmons in a graphene superlattice. Nat. Commun. 8, 1243. https://doi.org/10.1038/s41467-017-01205-z (2017).
He, L., Wu, J., Jin, J., Mele, E. J. & Zhen, B. Polaritonic Chern insulators in monolayer semiconductors. Phys. Rev. Lett. 130, 043801. https://doi.org/10.1103/PhysRevLett.130.043801 (2023).
Yang, Z. et al. Topological acoustics. Phys. Rev. Lett. 114, 114301. https://doi.org/10.1103/PhysRevLett.114.114301 (2015).
Xue, H., Yang, Y. & Zhang, B. Topological acoustics. Nat. Rev. Mater. 7, 974–990. https://doi.org/10.1038/s41578-022-00465-6 (2022).
Delplace, P., Marston, J. B. & Venaille, A. Topological origin of equatorial waves. Science 358, 1075–1077. https://doi.org/10.1126/science.aan8819 (2017).
Perrot, M., Delplace, P. & Venaille, A. Topological transition in stratified fluids. Nat. Phys. 15, 781–784. https://doi.org/10.1038/s41567-019-0561-1 (2019).
Wang, Z., Chong, Y. D., Joannopoulos, J. D. & Soljačić, M. Reflection-free one-way edge modes in a gyromagnetic photonic crystal. Phys. Rev. Lett. 100, 013905. https://doi.org/10.1103/PhysRevLett.100.013905 (2008).
Skirlo, S. A., Lu, L. & Soljačić, M. Multimode one-way waveguides of large Chern numbers. Phys. Rev. Lett. 113, 113904. https://doi.org/10.1103/PhysRevLett.113.113904 (2014).
Liu, G.-G. et al. Topological Chern vectors in three-dimensional photonic crystals. Nature 609, 925–930. https://doi.org/10.1038/s41586-022-05077-2 (2022).
Silveirinha, M. G. Bulk-edge correspondence for topological photonic continua. Phys. Rev. B 94, 205105. https://doi.org/10.1103/PhysRevB.94.205105 (2016).
Sakoda, K. et al. Photonic bands of metallic systems. I. Principle of calculation and accuracy. Phys. Rev. B 64, 045116. https://doi.org/10.1103/PhysRevB.64.045116 (2001).
Ito, T. & Sakoda, K. Photonic bands of metallic systems. II. Features of surface plasmon polaritons. Phys. Rev. B 64, 045117. https://doi.org/10.1103/PhysRevB.64.045117 (2001).
Raman, A. & Fan, S. Perturbation theory for plasmonic modulation and sensing. Phys. Rev. B 83, 205131. https://doi.org/10.1103/PhysRevB.83.205131 (2011).
Swanson, D. G. Plasma Waves, chap. 2 2nd edn. (CRC Press, Boca Raton, 2003).
Kuzmiak, V. & Maradudin, A. A. Photonic band structures of one- and two-dimensional periodic systems with metallic components in the presence of dissipation. Phys. Rev. B 55, 7427–7444. https://doi.org/10.1103/PhysRevB.55.7427 (1997).
Gangaraj, S. A. H. & Monticone, F. Do truly unidirectional surface plasmon-polaritons exist?. Optica 6, 1158–1165. https://doi.org/10.1364/OPTICA.6.001158 (2019).
Buddhiraju, S. et al. Absence of unidirectionally propagating surface plasmon-polaritons at nonreciprocal metal-dielectric interfaces. Nat. Commun. 11, 674. https://doi.org/10.1038/s41467-020-14504-9 (2020).
Howard, J. Introduction to Plasma Physics I—Collated Notes, 106 (Flooved.com, 2013), 1 edn. https://archive.org/details/flooved3263.
Sakoda, K. Optical Properties of Photonic Crystals 20th edn. (Springer, Berlin, 2005). https://doi.org/10.1007/b138376.
Kuzmiak, V., Maradudin, A. A. & Pincemin, F. Photonic band structures of two-dimensional systems containing metallic components. Phys. Rev. B 50, 16835–16844. https://doi.org/10.1103/PhysRevB.50.16835 (1994).
Kuzmiak, V. & Maradudin, A. A. Distribution of electromagnetic field and group velocities in two-dimensional periodic systems with dissipative metallic components. Phys. Rev. B 58, 7230–7251. https://doi.org/10.1103/PhysRevB.58.7230 (1998).
Degirmenci, E. & Landais, P. Finite element method analysis of band gap and transmission of two-dimensional metallic photonic crystals at terahertz frequencies. Appl. Opt. 52, 7367–7375. https://doi.org/10.1364/AO.52.007367 (2013).
Yan, W., Faggiani, R. & Lalanne, P. Rigorous modal analysis of plasmonic nanoresonators. Phys. Rev. B 97, 205422. https://doi.org/10.1103/PhysRevB.97.205422 (2018).
Xiao, W. & Sun, J. Band structure calculation of photonic crystals with frequency-dependent permittivities. J. Opt. Soc. Am. A 38, 628–633. https://doi.org/10.1364/JOSAA.412235 (2021).
Maschhoff, K. J. & Sorensen, D. C. P_ARPACK: An efficient portable large scale eigenvalue package for distributed memory parallel architectures. In Applied Parallel Computing Industrial Computation and Optimization (eds Waśniewski, J. et al.) 478–486 (Springer, Berlin, 1996).
Lehoucq, R. B., Sorensen, D. C. & Yang, C. ARPACK Users’ Guide (Society for Industrial and Applied Mathematics, 1998). https://epubs.siam.org/doi/pdf/10.1137/1.9780898719628.
Davanço, M., Urzhumov, Y. & Shvets, G. The complex Bloch bands of a 2D plasmonic crystal displaying isotropic negative refraction. Opt. Exp. 15, 9681–9691. https://doi.org/10.1364/OE.15.009681 (2007).
Fietz, C., Urzhumov, Y. & Shvets, G. Complex k band diagrams of 3D metamaterial/photonic crystals. Opt. Exp. 19, 19027–19041. https://doi.org/10.1364/OE.19.019027 (2011).
Raman, A. & Fan, S. Photonic band structure of dispersive metamaterials formulated as a Hermitian eigenvalue problem. Phys. Rev. Lett. 104, 087401. https://doi.org/10.1103/PhysRevLett.104.087401 (2010).
Parisi, G., Zilio, P. & Romanato, F. Complex Bloch-modes calculation of plasmonic crystal slabs by means of finite elements method. Opt. Exp. 20, 16690–16703. https://doi.org/10.1364/OE.20.016690 (2012).
Fang, C., Gilbert, M. J. & Bernevig, B. A. Bulk topological invariants in noninteracting point group symmetric insulators. Phys. Rev. B 86, 115112. https://doi.org/10.1103/PhysRevB.86.115112 (2012).
Chong, Y. D., Wen, X.-G. & Soljačić, M. Effective theory of quadratic degeneracies. Phys. Rev. B 77, 235125. https://doi.org/10.1103/PhysRevB.77.235125 (2008).
Acknowledgements
A.V.K. acknowledges fruitful discussions with M.G. Silveirinha. This work was partly supported by the Air Force Office of Scientific Research through Grant FA9550-21-1-0299. Work by T.C. is supported by Villum Fonden (42106).
Author information
Authors and Affiliations
Contributions
Y.J. and B.Z. conceived the project. C.Q., Y.J., and A.V.K. performed numerical simulations assisted by J.J. and T.C. C.Q., Y.J., A.V.K., and B.Z. wrote the paper with input from all authors and C.Q. prepared figures 1–4. All authors discussed the results. B.Z. supervised the project.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Qian, C., Jiang, Y., Jin, J. et al. Topological electromagnetic waves in dispersive and lossy plasma crystals. Sci Rep 13, 20445 (2023). https://doi.org/10.1038/s41598-023-47848-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-47848-5
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.