Abstract
We study the dispersions of the guided modes in the continuous uniaxial crystal slab waveguide and engineer their degeneracies through dielectric anisotropy. By switching the uniaxial positivity and negativity, we can obtain distinctive nodal types, point and line, for the lowest degeneracy in frequency. The mirror symmetry protections, and , are pointed out, and the degeneracy properties are intuitively analyzed through comparing the approximate slopes of the guided modes. Our results reveal a link between the lowest nodal types and the positivity/negativity of the uniaxial crystal, and provide a new approach to regulate the topology of degeneracy in two-dimensional photonic bands.
Introduction
Recently, a research focus is concentrated on nodal degeneracy in band diagrams in periodical structures [1–5], such as quantum materials, photonic crystals, and phononic crystals, because new physics and novel applications are anticipated arising from peculiar band degeneracies including point degeneracy [6–16], line/loop degeneracy [17–28], nodal chain degeneracy [29–31], nodal surface degeneracy [32, 33] and so on.
Band degeneracy is generally enforced by symmetry in the physical system. It is well known that photonic guided waves in confined structures can be classified into transverse electric (TE) and transverse magnetic (TM) modes in terms of mirror symmetry [our meanings of “TE” and “TM” adapted to classical waveguide theory, see [34]. Taking an example of an isotropic dielectric slab waveguide, both modes evolve out from the light cone in free space, and their dispersion curves are rapidly asymptotic to the light cone in dielectrics as propagation constant increasing. Because of the asymptotic parallelism, the TE and TM dispersions cross rarely to form the degeneracy beyond the free space light cone, which is schematically depicted in Figure 1B. On the other hand, the slopes of the TE and TM modes in the dispersion diagram can be tailored in a polarization-distinguishable way through introducing the uniaxial anisotropy to the dielectric slab. The tailoring mechanism is rooted on the refractive index difference which is manifested likewise in the propagation of ordinary and extra-ordinary light in a uniaxial bulk crystal [35], seeing the index ellipsoid in Figure 1.
FIGURE 1
In the study, we start with an isotropic dielectric slab, where TE and TM modes are not degenerate. By changing the dielectric constant into the uniaxial permittivity tensor and tuning the component of the tensor along the propagation direction, we show the slope of TM modes can either increase or decrease significantly while maintaining the slope of TE modes, which corresponds to the positively and negatively uniaxial anisotropy, respectively. Thus, the crossing between TE and TM modes can be engineered, as illustrated in Figure 1. Furthermore, the mirror symmetry along the out-of-slab direction imposes an extra enforcement on the degeneracy and leads to distinctive nodal types, Dirac point (DP) and Dirac line (DL), in the positive and negative uniaxial cases for the lowest degeneracy in frequency.
Guided mode in uniaxial crystal slabs
Here, we consider a two-dimensional (2D) infinite (along y- and z-directions), uniaxial crystal slab (finite thickness d = 2 mm in the x-direction) with non-magnetic permeability (; being the permeability in vacuum). The slab is located in free space where the wave is assumed propagating along the z-direction, and the permittivity tensor has the diagonal form with () being the refractive index of ordinary (extra-ordinary) light and being the permittivity in vacuum. The dielectric principal axis in the uniaxial crystal is spanned by , shown in Figure 2A, where denotes the in-plane rotation of the optical axis. When , the anisotropic permittivity tensor of the slab can be expressed in the coordinate system as:where is the permittivity along the direction with relative value , , and . The dielectric loss of the permittivity is neglected in the study. The time harmonic waves that propagate in the z-direction can be expressed as:where represents the component of wave vector in the direction, has been assumed for homogeneity in the y-direction, and is angular frequency. Since the system shows the mirror symmetry , the guided wave is cataloged into two polarization modes, TE with non-zero electric field perpendicular to the mirror plane (any xz plane due to uniformity in the y-direction) and TM with non-zero electric field parallel to the mirror plane (any xz plane due to uniformity in the y-direction) [Ref. 34]. By expressing the field components and matching the boundary conditions on the two surfaces of the slab, we can get the characteristic equations for TE mode:
FIGURE 2
and for TM mode:
Above, Eqs 4, 6 are the dispersion relation in regime I, α is the imaginary part of perpendicular component of wave vector in theregime and , and with the speed of light in vacuum , which satisfies:
For concreteness, we choose , and switch the value of for negatively uniaxial case () and positively uniaxial case (). It is has known that 2D dielectric slabs have been extensively used as basic waveguides in microwave engineering and devices where a broad horizon of dielectric materials, e.g., high-k printed circuit board (PCB) and ceramics, may offer various permittivity including such values [36, 37]. The calculated results are shown in Figures 2B,C. Within a qualitative physical picture, we approximate the slope of TM modes beyond the light line to be roughly , and the effective index , where the effective permittivity for TM mode can be regarded to some degree as special average of ; because the electric field is oriented along both x- and z-direction. In contrast, the slope, being roughly , of TE modes is related to in terms of the y-orientation of electric field. Therefore, given and when switching only from the negatively uniaxial case () to positively uniaxial case (), we see that the slopes of the TE modes remain almost unchanged whereas the TM ones change expectedly in the dispersion diagram, as comparing Figure 2B with Figure 2C.
It is noted that the first-order TM mode (TM1) crosses with the second-order TE mode (TE2) for the negatively uniaxial case (), as labeled by point in Figure 2B, and that it crosses with the first-order TE mode for positively uniaxial case (), as labeled by point in Figure 2C, where TM curves are generally less steep than TE ones. For instance, the slope of TM1 near the point is estimated as that corresponds to ∼25, whereas the slope of TE1 near the point is estimated as that corresponds to ∼16. Therefore, the decrease in the slope of TM modes with leads to the switch of one of degenerating bands from TE2 (point ) to TE1 (point ). In addition, as increasing , more TM modes appear in Figure 2C, and TM3 is crossing with TE2, which gives rise to more degeneracies at higher frequencies (see Section A in Supplementary Materials).
Type-II Dirac degeneracy in uniaxial crystal slabs
In order to exhibit the complete dispersion structure around the degenerate points, we need calculate the band diagram . In the calculation, we first rotate the in-plane dielectric principal axes around the -direction with the angle , as shown in Figure 2A. The non-diagonalized permittivity tensor after rotation can be written as
in which
Then, we assume that the waves still propagate along the z-direction and express the electric fields and magnetic fields in different regions. Because the mirror symmetry is broken under the rotation, the guided modes are no longer pure TE or TM mode, but are their combination which we call hybrid mode. The characteristic equation for hybrid modes is solved by matching boundary conditions, which gives us the dispersion . Finally, the dispersion is obtained through a standard map from polar coordinate to Cartesian coordinate. Although the calculation is based on the rotated dielectric principal axis, the results are the same as those of rotating the coordinate system while maintaining the dielectric principal axis, because both are the equivalent description of rotation.
We first analyze the degenerate point in negatively uniaxial case with , where it is the lowest degeneracy in frequency, as displayed in Figure 2B. We calculate the dispersion when varying the angle , and the degeneracy is seen to become gapped upon being nonzero, as shown in Figures 3A–C. Combining all -cut plots, we can achieve the three-dimensional (3D) view of band diagram in momentum space, as depicted in Figure 3D. The band structure around exhibits the characteristic of two over-tilted cones, and thus the degeneracy is just the type-II DP. The gapping reason is that the two modes display the like parity in the mirror symmetry (the symmetry classifies all modes as even or odd parity with respect to the mirror plane x = 0, see Section B in Supplementary Materials) and simultaneously is broken upon . Such 2D type-II DP has been observed in the artificially designed metasurfaces with periodic metallic patterns at microwave frequencies [38–40]. In contrast, our system is of no discrete translational symmetry, but continuous in space.
FIGURE 3
Next, let’s focus our attention on the degenerate point in positively uniaxial case with . Interestingly, in the cut plot with different , the degeneracy persists, which forms a section of line in momentum space, as plotted in Figure 4. According to the slope of two crossing bands, such degeneracy is the Type-II DL. The degeneracy is protected by the mirror symmetry , in which the two bands display the opposite parity (see Section B in Supplementary Materials), and would be gapped if the uniaxial crystal slab lies in an asymmetric background, (see Section C in Supplementary Materials).
FIGURE 4
It is also noted from Figure 4C that the two bands, responsible for the Dirac line, almost coincide with each other when comes to 30 . If we increase the rotation angle further, the DL degeneracy will lift. The physical reason can be understood qualitatively from an effective permittivity point of view. Although being hybrid in nature upon , the two bands are dominated, respectively, by the TE1 and TM1 modes and can be considered as their descendants. As increasing the angle , the slope of TE-like dispersion curve will decrease from at to at where the effective permittivity is the -mediated average between and , and we make the approximation (the relation expressed for extra-ordinary light in the bulk crystal [Ref. 35]). Consequently, the slope of TE-like dispersion may become comparable to that of TM-like dispersion upon some angle, and appear less than the latter beyond the angle, and such slope difference does not lead to the degeneracy any more. The threshold angle can be estimated through requiring the slope equality , where we assume after noticing insignificant -dependence of TM-like dispersions in Figure 4. Thus, the threshold angles are solved as and , which are quite close to the numerical values and from COMSOL simulation and correspond to the and directions, respectively. Overall, this kind of nodal line is essentially different from the loop degeneracy [Ref. 40] because it cannot construct a closed loop in momentum space.
Discussions
When investigating the propagation of the guided modes along the direction perpendicular to the optical axis (e.g., the propagation along y-axis while the optical axis being z-axis, and the mirror symmetry giving rise to the TE/TM modes), we find that the relative magnitude between the slopes of TM and TE modes will swap, because TE waves have the electric field polarized along the optical axis or and TM waves have the electric field components along both x- and y-directions. Thus, the DL degeneracy from the lowest two modes (even and odd parity) will occur to the negatively uniaxial case, and the DP degeneracy from the higher-order modes will appear in the positively uniaxial case (see Section D in Supplementary Materials).
In addition, we present a microwave metamaterial design which approximates at lower frequencies a positively uniaxial crystal with dispersive permittivity component, and the band degeneracy shows the similarity to and the difference from the case of the continuous crystal slab (see Section E in Supplementary Materials).
Conclusion
In conclusion, by controlling positivity/negativity of the refractive index ellipsoid, we can obtain the nodal point and nodal line degeneracies for the guided modes on the uniaxial crystal slab waveguide. Furthermore, the point and line characteristics and their connections with the refractive index ellipsoid can be swapped through switching the propagation direction. Our results link the band degeneracy with positivity/negativity of the uniaxial crystal, and provide a new approach to regulate the topology of degeneracy in 2D photonic bands.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
BH conceived and supervised the research; XP performed the research; HL, WD, XZ, KX, CH, and GW assisted in analyzing the data; and BH, XP, and HL wrote the manuscript.
Funding
This work was supported by the Natural Science Foundation of China (NSFC) (Grant No. 12074279), the Major Program of Natural Science Research of Jiangsu Higher Education Institutions (Grant No.18KJA140003), and the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions.
Conflict of interest
Chuandeng Hu was employed by Shenzhen Fantwave Tech Co., Ltd
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2022.1095669/full#supplementary-material
References
1.
ArmitageNPMeleEJVishwanathA. Weyl and Dirac semimetals in three-dimensional solids. Rev Mod Phys (2018) 90:015001. 10.1103/revmodphys.90.015001
2.
LiSYuZ-MYaoYYangSA. Type-II topological metals. Front Phys (2020) 15(4):43201. 10.1007/s11467-020-0963-7
3.
LuLJoannopoulosJDSoljačićM. Topological photonics. Nat Photon (2014) 8:821–9. 10.1038/nphoton.2014.248
4.
KimMJacobZRhoJ. Recent advances in 2D, 3D and higher-order topological photonics. Light Sci Appl (2020) 9:130. 10.1038/s41377-020-0331-y
5.
ParkHDGaoWZhangXOhSS. Nodal lines in momentum space: Topological invariants and recent realizations in photonic and other systems. Nanophotonics (2022) 11(11):2779–801. 10.1515/nanoph-2021-0692
6.
LuLWangZ-YYeD-XRanL-XFuLJoannopoulosJDet alExperimental observation of Weyl points. Science (2015) 349:622–4. 10.1126/science.aaa9273
7.
GaoWYangBLawrenceMFangFBérilBZhangS. Photonic Weyl degeneracies in magnetized plasma. Nat Commun (2016) 7:12435. 10.1038/ncomms12435
8.
YangBGuoQTremainBBarrLEGaoWLiuHet alDirect observation of topological surface-state arcs in photonic metamaterials. Nat Commun (2017) 8:97. 10.1038/s41467-017-00134-1
9.
YangBGuoQTremainBLiuRBarrLEYanQet alIdeal Weyl points and helicoid surface states in artificial photonic crystal structures. Science (2018) 359:1013–6. 10.1126/science.aaq1221
10.
QiuPQiuWLinZChenHRenJWangJ-Xet alDouble Dirac point in a photonic graphene. J Phys D: Appl Phys (2017) 50:335101. 10.1088/1361-6463/aa7bc2
11.
XiaoMLinQFanS. Hyperbolic Weyl point in reciprocal chiral metamaterials. Phys Rev Lett (2016) 117:057401. 10.1103/physrevlett.117.057401
12.
YangYGaoZFengXHuangY-XZhouPYangSAet alIdeal unconventional Weyl point in a chiral photonic metamaterial. Phys Rev Lett (2020) 125:143001. 10.1103/physrevlett.125.143001
13.
LiuG-GZhouPYangYXueHRenXLinXet alObservation of an unpaired photonic Dirac point. Nat Commun (2020) 11:1873. 10.1038/s41467-020-15801-z
14.
LuoLDengWYangYYanMLuJHuangXet alObservation of quadruple Weyl point in hybrid-Weyl phononic crystals. Phys Rev B (2022) 106:134108. 10.1103/physrevb.106.134108
15.
HuangXLaiYHangZHZhengHChanCT. Dirac cones induced by accidental degeneracy in photonic crystals and zero-refractive-index materials. Nat Mater (2011) 10:582–6. 10.1038/nmat3030
16.
XuCMaGChenZ-GLuoJShiJLaiYet alThree-dimensional acoustic double-zero-index medium with a fourfold degenerate Dirac-like point. Phys Rev Lett (2020) 124:074501. 10.1103/physrevlett.124.074501
17.
LinJYHuNCChenYJLeeCHZhangX. Line nodes, Dirac points, and Lifshitz transition in two-dimensional nonsymmorphic photonic crystals. Phys Rev B (2017) 96:075438. 10.1103/physrevb.96.075438
18.
YangS-YYangHDerunovaEParkinSSPYanBAliMN. Symmetry demanded topological nodal-line materials. Adv Phys (2018) 3:1414631. 10.1080/23746149.2017.1414631
19.
FengBFuBKasamatsuSItoSChengPLiuC-Cet alExperimental realization of two-dimensional Dirac nodal line fermions in monolayer Cu2Si. Nat Commun (2017) 8:1007. 10.1038/s41467-017-01108-z
20.
FengBZhangR-WFengYFuBWuSMiyamotoKet alDiscovery of Weyl nodal lines in a single-layer ferromagnet. Phys Rev Lett (2019) 123:116401. 10.1103/physrevlett.123.116401
21.
WangS-SYuZ-MLiuYJiaoYGuanSShengX-Let alTwo-dimensional nodal-loop half-metal in monolayer MnN. Phys Rev M (2019) 3:084201. 10.1103/physrevmaterials.3.084201
22.
YouJ-YChenCZhangZShengX-LYangSASuG. Two-dimensional Weyl half-semimetal and tunable quantum anomalous Hall effect. Phys Rev B (2019) 100:064408. 10.1103/physrevb.100.064408
23.
GaoWYangBTremainBLiuHGuoQXiaLet alExperimental observation of photonic nodal line degeneracies in metacrystals. Nat Commun (2018) 9:950. 10.1038/s41467-018-03407-5
24.
QiuHQiuCYuRXiaoMHeHYeLet alStraight nodal lines and waterslide surface states observed in acoustic metacrystals. Phys Rev B (2019) 100:041303. 10.1103/physrevb.100.041303
25.
DengWLuJLiFHuangXYanMMaJet alNodal rings and drumhead surface states in phononic crystals. Nat Commun (2019) 10:1769. 10.1038/s41467-019-09820-8
26.
XiongZZhangR-YYuRChanCTChenY. Hidden symmetry enforced nexus points of nodal lines in layer-stacked dielectric photonic crystals. Light Sci Appl (2020) 9:176. 10.1038/s41377-020-00382-9
27.
HuMZhangYJiangXQiaoTWangQZhuSet alDouble-bowl state in photonic Dirac nodal line semimetal. Light Sci Appl (2021) 10:170. 10.1038/s41377-021-00614-6
28.
DengW-MChenZ-MLiM-YGuoC-HTianZ-TSunK-Xet alIdeal nodal rings of one-dimensional photonic crystals in the visible region. Light Sci Appl (2022) 11:134. 10.1038/s41377-022-00821-9
29.
YanQLiuRYanZLiuBChenHWangZet alExperimental discovery of nodal chains. Nat Phys (2018) 14(5):461–4. 10.1038/s41567-017-0041-4
30.
WangYZhouXLiSZhangWHuCLuWet alFlatness and boundness of photonic drumhead surface state in a metallic lattice. Sci Rep (2021) 11:8684. 10.1038/s41598-021-88004-1
31.
WuXFanHLiuTGuZZhangR-YZhuJet alTopological phononics arising from fluid-solid interactions. Nat Commun (2022) 13:6120. 10.1038/s41467-022-33896-4
32.
XiaoMYeLQiuCHeHLiuZFanS. Experimental demonstration of acoustic semimetal with topologically charged nodal surfaceSci. Adv. (2020) 6:eaav2360.
33.
YangYXiaJ-PSunH-XGeYJiaDYuanS-Qet alObservation of a topological nodal surface and its surface-state arcs in an artificial acoustic crystal. Nat Commun (2019) 10:5185. 10.1038/s41467-019-13258-3
34.
JoannopoulosJDMeadeRDJohnsonSGJoshuaN. Photonic crystals: Molding the flow of light. 2nd ed.Princeton, NJ: Princeton University Press (2011).
35.
SalehBEATeichMC. Fundamentals of photonics. 2nd ed. New York, NY: John Wiley & Sons (2007).
36.
PozerDMAppendixG. Microwave engineering. 2nd ed. New York, NY: John Wiley & Sons (1998).
37.
ROGERSPCB MAIN (2020). Rogers laminates in stock. Available at: https://www.rogerspcb.com.cn/.
38.
HuCLiZTongRWuXXiaZWangLet alType-II Dirac photons at metasurfaces. Phys Rev Lett (2018) 121:024301. 10.1103/physrevlett.121.024301
39.
LiHHuCZhouXLuWHouBWenW. Conical diffraction from type-II Dirac point at metasurfaces. EPL (2020) 130:17007. 10.1209/0295-5075/130/17007
40.
LiHHuCJiangJ-HWuJWenWHouB. Photonic type-III nodal loop and topological phase transitions at bilayer metasurfaces. Front Mater (2022) 9:909381. 10.3389/fmats.2022.909381
Summary
Keywords
nodal line, nodal points, uniaxial crystal, topological transformation, band degeneracy, symmetry, guided mode
Citation
Pan X, Li H, Dong W, Zhou X, Xing K-A, Hu C, Wang G and Hou B (2023) Nodal degeneracy of guided modes in uniaxial crystal slabs. Front. Phys. 10:1095669. doi: 10.3389/fphy.2022.1095669
Received
11 November 2022
Accepted
15 December 2022
Published
06 January 2023
Volume
10 - 2022
Edited by
Zhiwei Guo, Tongji University, China
Reviewed by
Biao Yang, National University of Defense Technology, China
Xiao-Dong Chen, Sun Yat-sen University, China
Updates
Copyright
© 2023 Pan, Li, Dong, Zhou, Xing, Hu, Wang and Hou.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bo Hou, houbo@suda.edu.cn
This article was submitted to Optics and Photonics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.