Photonic helicoid-like surface states in chiral metamaterials

We investigate the photonic topological phases in chiral metamaterials characterized by the magnetoelectric tensors with diagonal chirality components. The underlying medium is considered a photonic analogue of the topological semimetal featured with a Weyl cone and a cylindrical surface in the frequency-wave vector space. As the ’spin’-degenerate condition is satisfied, the photonic system can be rearranged as two hybrid modes that are completely decoupled. By introducing the pseudospin states as the basis for the hybrid modes, the photonic system is described by two subsystems in the form of spin-orbit Hamiltonians of spin 1, which result in nonzero spin Chern numbers that determine the topological properties. Surface modes at the interface between vacuum and the chiral metamaterial exist in their common gap in the wave vector space, which are analytically formulated by algebraic equations. In particular, the surface modes form a pair of spiral surface sheets wrapping around the Weyl cone, resembling the helicoid surface states that occur in topological semimetals. At the Weyl frequency, the surface modes contain two Fermi arc-like states that concatenate to yield a straight line segment.


Ruey-Lin Chern
We investigate the photonic topological phases in chiral metamaterials characterized by the magnetoelectric tensors with diagonal chirality components.The underlying medium is considered a photonic analogue of the topological semimetal featured with a Weyl cone and a cylindrical surface in the frequency-wave vector space.As the 'spin'-degenerate condition is satisfied, the photonic system can be rearranged as two hybrid modes that are completely decoupled.By introducing the pseudospin states as the basis for the hybrid modes, the photonic system is described by two subsystems in the form of spin-orbit Hamiltonians of spin 1, which result in nonzero spin Chern numbers that determine the topological properties.Surface modes at the interface between vacuum and the chiral metamaterial exist in their common gap in the wave vector space, which are analytically formulated by algebraic equations.In particular, the surface modes form a pair of spiral surface sheets wrapping around the Weyl cone, resembling the helicoid surface states that occur in topological semimetals.At the Weyl frequency, the surface modes contain two Fermi arc-like states that concatenate to yield a straight line segment.
Topological phases are new phases of matter characterized by integer quantities known as topological invariants, which remain constant under arbitrary continuous deformations of the system.Quantum Hall (QH) state 1 is the very first example of two dimensional (2D) topological phase, belonging to the class with broken time-reversal (TR) symmetry because of the presence of a static magnetic field.Quantum spin Hall (QSH) state [2][3][4] is a different 2D topological phase without the magnetic field and preserves the TR symmetry, where the spin-orbit coupling is responsible for the topological characters.Topological properties of the QH states are characterized by TKNN invariants or Chern numbers 5 , while those of the QSH states are characterized by Z 2 invariants 2 or spin Chern numbers 6 .Theoretical concepts developed in the QSH states are generalized to three dimensions (3D), leading to the more general class of 3D topological insulators 7,8 .
One remarkable feature of the QSH state is the emergence of gapless edge states inside the bulk band gap.The propagation direction of the edge states is locked by the spin 9 , which enables topologically protected edge states that propagate unidirectionally without back scattering 10 .As the edge states are protected by the bulk topology, they are insensitive to small perturbations that do not change the topology.Similar to the case of 2D topological phases, gapless surface states appear inside the band gap between two topologically distinct bands in 3D topological insulators 11,12 , which can be realized in both TR broken 13,14 and TR invariant systems [15][16][17] .In contrast to 3D topological insulators that are gapped topological phases, 3D gapless topological phases are a new type of phases known as topological semimetals [18][19][20][21][22] .
A majority of the topological semimetals are characterized by Weyl degeneracies, which are degeneracies between topologically inequivalent bands.The main signature of 3D gapless topological phases is the appearance of Weyl points existing in the systems that lack TR symmetry, inversion symmetry, or both.The Weyl points are understood as the monopoles of Berry curvature in the momentum space that carry quantized topological charges, which are equal to the topological invariants of the system.A useful perspective on the Weyl semimetals is to view them as the transitional state between a topological insulator and a trivial insulator 22 .An important feature of the Weyl points is the existence of Fermi arcs that connect the Weyl points, corresponding to the topologically protected surface states that are robust against disorder.In particular, surface states may form a spiral surface sheet that connects the upper and lower bulk cones, which are protected from being gapped by nonsymmorphic symmetries and termed as helicoid surface states 23 .
A spin-like quantity called pseudospin can be formed by the linear combination of electric and magnetic fields when a specific degenerate condition between the electric and magnetic parameters is satisfied 32 .As a result, the photonic system can be described by an effective Hamiltonian consisting of two subsystems for the pseudospin states [32][33][34] , and the Kramers pair can be formed in the photonic system.In the presence of chirality or bianisotropy that emulates the spin-orbit coupling, a topological phase can be constructed in the photonic system [52][53][54][55][56] .
In the present study, we investigate the photonic topological phases in chiral metamaterials characterized by the magnetoelectric tensors with diagonal chirality components.Bulk modes of the underlying medium are represented by two decoupled quadratic equations when a certain symmetry of the material parameters is included.As the 'spin'-degenerate condition 32,34,38 is satisfied, the bulk modes are featured with a Weyl cone and a cylindrical surface in the frequency-wave vector space.The electromagnetic duality allows for Maxwell's equations to be decoupled into two subsystems for the hybrid modes, which are defined as the linear combinations of electric and magnetic fields.By introducing the pseudospin states as the basis for the hybrid modes, the photonic system can be described by a pair of spin-orbit Hamiltonians of spin 1 [53][54][55][56][57][58][59] that respect the fermioniclike pseudo TR symmetry.The topological properties of the photonic system are determined by nonzero spin Chern numbers calculated from the eigenfields of the Hamiltonians.Surface modes at the interface between vacuum and the chiral metamaterial exist in their common gap in the wave vector space, which are analytically formulated by algebraic equations.In particular, the surface modes are represented by a pair of spiral surface sheets wrapping around the Weyl cone, resembling the helicoid-like surface states that occur in the topological semimetals 23 .At the Weyl frequency, the surface modes contain two Fermi arc-like states that concatenate to yield a straight line segment.

Results
Bulk modes.Consider a general bianisotropic medium characterized by the constitutive relations: where ε , µ , ξ and ζ are the frequency-dependent relative permittivity, permeability, and magnetoelectric tensors, respectively.Treating the combined electric field E = E x , E y , E z T and magnetic field H = H x , H y , H z T as a six-component vector, where T denotes the transpose, Maxwell's equations for the time-harmonic electromagnetic fields (with the time convention e −iωt ) are written in matrix form as where I is the 3 × 3 identity matrix and H ′ = η 0 H , with η 0 = √ µ 0 /ε 0 .Let the medium be lossless ( ε = ε † , µ = µ † , and ξ = ζ † , where † denotes the Hermitian conjugate) and reciprocal ( ε = ε T , µ = µ T , and ξ = −ζ T , where T denotes the transpose) 60 , which implies that ε = ε * , µ = µ * , ξ = −ξ * , and ζ = −ζ * , where * denotes the complex conjugate.In the present study, we further assume that the permittivity, permeability, and magnetoelectric tensors are uniaxial: ε = diag(ε t , ε t , ε z ) , µ = diag(µ t , µ t , µ z ) , and ξ = −ζ = diag(iγ t , iγ t , iγ z ) , where ε n , µ n , and γ n (n = t, z) are real-valued quantities, and diag(•, •, •) denotes a 3 × 3 diagonal matrix.The medium with the purely imaginary magnetoelectric tensors is known as the chiral medium.Here, the chirality parameter γ n ( n = t, z ) appears in the diagonal elements of the magnetoelectric tensors ξ and ζ , which means that the magnetoelectric couplings of the chiral medium occur in parallel directions.The underlying medium can be synthesized by metallic helices oriented along three perpendicular directions 61 .In the chiral medium, the inversion symmetry is broken because of the chirality 44,62 , whereas the TR symmetry is preserved 24 .
The existence of a nontrivial solution of E and H requires that the determinant of the 6 × 6 matrix in Eq. (3) be zero, giving rise to the characteristic equation of bulk modes as where k 0 = ω/c .This is a bi-quadratic equation that incorporates the coupling between transverse electric and transverse magnetic modes.If 4) can be decoupled as a product of two quadratic equations as In the isotropic case, where ε t = ε z ≡ ε , µ t = µ z ≡ µ , and γ t = γ z ≡ γ , Eq. ( 5) can be simplified to (1) where n ± = √ εµ ± |γ | .There exists a critical condition: |γ | = √ εµ , where n − = 0 and the corresponding bulk mode is reduced to a point.In case γ = 0 , Eq. ( 6) is further simplified to Note that the features of bulk modes may change with the frequency for a dispersive medium (usually the case of metamaterials), depending on the choice of frequency range.In the neighborhood of a reference frequency where εn is posi- tive definite 57 .A similar relation is valid for µ n ( n = t, z ).We further assume that the chirality parameter γ varies smoothly around ω ref and can be treated as a constant in the analysis 32,53,58 .
Spin-orbit Hamiltonians The electromagnetic duality of Maxwell's equations dictates that the matrix in Eq. ( 3) holds a symmetric pattern when the 'spin'-degenerate condition ε = µ 32,34,38 is satisfied.This allows us to rewrite Eq. ( 3) as where H ± 0 = ∓ωε + i ck × I + ωξ and F ± = E ± iH ′ are the hybrid modes that linearly combine the electric and magnetic fields.Note that F + and F − are completely decoupled and determined by two subsystems ( 3 × 3 matrices) with a similar form.By introducing the pseudospin states ψ ± = U −1 ψ± as the basis for the hybrid modes, where ψ± = − T and U = diag εz /ε t , 1, εz /ε t , Eq. ( 8) can be formulated as a pair of eigensystems when the frequency dispersion of the medium near the reference frequency ω ref is taken into account.In the isotropic case, where ε t0 = ε z0 ≡ ε and εt = εz ≡ ε , the eigensystems for Eq. ( 8) are given by (see Methods A) where and ) being the spin matrices of spin 1.Note that Eq. ( 9) is expressed as an eigensystem with ±δω being the eigenvalue.The Hamiltonian H ± in Eq. ( 10) represents the spin-orbit coupling k • S of spin 1, which is mathematically equivalent to the Hamiltonian of a magnetic dipole moment in the magnetic field 57 .

Topological invariants
The topological properties of the spin-orbit Hamiltonians H ± can be characterized by the topological invariants based on the eigenfields.For this purpose, we calculate the Berry flux over a closed surface in the wave vector space.The eigensystem for the Hamiltonian H ± in Eq. (10): is solved to give the eigenvalues σ ± and eigenvectors ψ σ ± ( σ = ±1, 0 ).Here, the eigenvalue σ ± is related to δω in Eq. ( 9) as σ ± = D ± ± δω .Based on Eq. ( 11), the Chern numbers are calculated to give (see Methods B) The nonzero C σ ( σ = ±1 ) characterizes the topological properties of the system, where σ refers to the helicity (or handedness) of the pseudospin states.In particular, the edge or surface states at the interface between two distinct topological phases are topologically protected, which means that their existence is guaranteed by the difference in band topology on two sides of the interface.In this system, the total Chern number C = σ C σ = 0 and the spin Chern number C spin = σ σ C σ = 4 , which are consistent with the quantum spin Hall effect of light 51 .The topological invariants remain constant under arbitrary continuous deformations of the system.The topological properties in the isotropic case will be retained when a certain anisotropy is included in the system.For a more general anisotropic case, the exact calculation of topological invariants can be obtained by numerical integration of the Berry curvatures 63 .
Pseudo time-reversal symmetry.The Hamiltonian for Maxwell's equations [cf.Eq. ( 3)] in the chiral medium, which is lossless and reciprocal, is TR invariant under T b , that is, where (6) ) is the bosonic TR operator for photons 24 , with K being the complex conjugation, and ⊗ denotes the tensor product.The Hamiltonian H m , however, is not TR invariant under T f , that is, , where T f = iσ y K (with T 2 f = −1 ) is the fermionic TR operator for electrons 24 .Nevertheless, the combined Hamiltonian formed by two spin-orbit Hamiltonians H ± [cf.Eq. (10)] is TR invariant under T p , that is, where and T p is the fermionic-like pseudo TR operator having the same form of T f .The pseudo TR operator T p is inspired by noticing that E ↔ H during the TR operation, which is defined as Here, σ x = (0, 1; 1, 0) , σ y = (0, −i; i, 0) , and σ z = diag(1, −1) are the Pauli matrices.The pseudo TR symmetry of the combined Hamiltonian H c is crucial in constructing the Kramers pair and determining the topological phases in photonic systems of spin 1, allowing for the existence of spin-polarized counterpropagating helical edge states as in electronic systems.
Surface modes.Topological phase transition occurs between two distinct topological phases.In the present study, the topological phase transition is manifest on the existence of surface modes at the interface between vacuum and the chiral metamaterial.Here, we begin with the linear combination of eigenfields based on Maxwell's equations, leaving the weighting coefficients as unknown variables.Applying the boundary condition at the interface, a condition can be arrived by requiring these coefficients to be nontrivial, which gives rise to the characteristic equation of surface modes.
Let the xz plane be an interface between vacuum ( y > 0 ) and the chiral medium ( y < 0 ) characterized by ε n = ε , µ n = µ , and γ n = γ ( n = t, z ), where the surface modes may exist.According to Maxwell's boundary conditions: the continuity of tangential electric and magnetic field components at the interface, the characteristic equation of surface modes can be analytically formulated by using the eigenfields of bulk modes on two sides of the interface, given by (see Methods C) z is the normal wave vector component (to the interface) in vacuum, k (1) z are the normal wave vector components in the chiral medium, and the superscripts (1) and (2) refer to two independent polarizations.

Discussion
Bulk modes.Figure 1a shows the equifrequency surfaces of the bulk modes for the chiral metamaterial in the wave vector space based on Eq. ( 5).In the present study, we assume that n ± t n ± z > 0 so that the bulk modes are described by two elliptic equations [cf.Eq. ( 5)].This condition is crucial to form the photonic Weyl system in the chiral metamaterial, which will be discussed later (cf.Results: Photonic Weyl system.).As a result, the bulk modes are represented by two concentric ellipsoids in the wave vector space.Note that the bulk modes for opposite sign of the chirality parameter are identical because of the symmetry about γ [cf.Eq. ( 5)].Here, the material parameters are arranged such that n The bulk modes are therefore either completely inside or completely outside the vacuum dispersion spheroid: , as shown in Fig. 1(b) for the bulk modes on the half space ( k y > 0 ).Note also that the bulk modes in Fig. 1a, b are represented by the same ellipsoids, although n − s > 0 ( s = t, z ) for the former and n − s < 0 for the latter.The wave propagations in the two cases, however, are different in the issue of negative refraction and backward wave 64,65 .In the isotropic case, the inner bulk mode at the critical condition: |γ | = √ εµ (cf.Results: Bulk modes.) is reduced to a point at the origin.
Recall that the effective Hamiltonian in the present problem consists of two subsystems of the hybrid modes.Each subsystem is described by the spin-orbit Hamiltonian with spin 1 (cf.Results: Spin-orbit Hamiltonians.)and characterized by nonzero topological invariants (cf.Results: Topological invariants.).In this regard, the chiral metamaterial is considered a photonic analogue of the topological phase.

Surface modes.
Figure 2 shows the surface modes at the interface between vacuum ( y > 0 ) and the chiral metamaterial ( y < 0 ) in the k x -k z plane based on Eq. ( 17).The bulk modes at k y = 0 are overlaid in the same (15) www.nature.com/scientificreports/plots.For clarity, we discuss the surface modes in the isotropic case, where ε n = ε and µ n = µ ( n = t, z ), and the analytical formulation for the surface modes are available.In particular, the surface modes are represented by a pair of curve segments symmetric about the center, which are located in the second and fourth quadrants for γ < √ εµ [cf.Fig. 2a], and the first and third quadrants for γ > √ εµ [cf.Fig. 2b].Note that the surface modes and bulk modes 'merge' at the points: (k x , k z ) = ±| √ εµ − γ |k 0 , 0 for the chiral metamaterial and (0, ±k 0 ) for vacuum.
The bulk modes for opposite sign of the chirality parameter are identical because of the symmetry about γ [cf.Eq. ( 4) or (5)].The surface modes are located in the common gap of bulk modes in the wave vector space and tangent to the bulk modes 18,22 , including the chiral metamaterial (black solid contour) and vacuum dispersion circle: k 2 x + k 2 z = k 2 0 (gray dashed contour).This feature follows from the fact that the surface modes must convert seamlessly into the bulk modes as they approach their termination points 66 .The evanescent depth of the surface mode grows until at the point where the surface mode merges with the bulk mode 22 .The bulk modes on the vacuum side are topologically trivial, while on the chiral medium side they are topologically nontrivial with nonzero topological invariants (cf.Results: Topological invariants.).The surface modes correspond to the topological phase transition between two distinct topological phases in the momentum space 53,67 , their existence being guaranteed by the bulk-edge correspondence.In particular, the Hamiltonian of the photonic system respects the pseudo TR symmetry (cf.Results: Pseudo time-reversal symmetry.),leading to the topological protection of photonic edge or surface states.www.nature.com/scientificreports/Photonic Weyl system.Let the frequency dependence of the chiral medium be characterized by the Lorentz dispersion models: , which are usually employed in the study of metamaterials 68 .Here, ω p is the plasma frequency and ω 0 is the resonance frequency.The chirality parameter is given by γ = � γ ωω p / ω 2 − ω 2 0 , where . This model guarantees that the energy density in the underlying medium is positive definite (see Methods D).
Figure 3a shows the dispersion of bulk modes for the chiral metamaterial in the frequency-wave vector space.The bulk modes consist of a conic surface at the center of a cylindrical surface.In the present configuration, the mate r i a l p ar ame te rs are ar r ange d su ch t hat ε , where the inner bulk modes are reduced to a point at (k t , k z , ω) = (0, 0, ω 1 ) .This is the condition that fulfills both the 'spin'-degenerate condition: ε = µ (cf.Results: Spin-orbit Hamiltoni- ans.) and the critical condition: |γ | = √ εµ (cf.Results: Bulk modes.) in the present medium, which also forms the point-like degeneracy in the bulk modes.Here, ω 1 is the transition frequency between the bulk modes with n − > 0 and n − < 0 (cf.Results: Bulk modes.), and the bulk modes for either ω > ω 1 or ω < ω 1 are represented by similar elliptic curves.The former and the latter touch at a degenerate point, forming the conic surface for the inner bulk mode with n − , while the outer bulk mode with n + (always positive) is a cylindrical surface.In this situation, the dispersion branch of the inner bulk modes resembles the linear crossing of valence and conduction bands in the Weyl semimetal 71 , with the crossing point known as the Weyl point and the associated frequency ω 1 as the Weyl frequency.The topological charge associated with the Weyl point is consistent with the nonzero topological invariants of the system (cf.Results: Topological invariants.).Note that the inner bulk mode is reduced to a single point at the Weyl frequency.In this regard, the underling medium is considered a photonic analogue of the type-I Weyl semimetal 22 .
Note that in the absence of chirality ( γ = 0 ), the bulk modes are featured with the Dirac cone with fourfold degeneracy at the Dirac point: (k x , k y , ω) = (0, 0, ω 1 ) in the wave vector space [cf.Eq. ( 7)].In the presence of chi- rality ( γ = 0 ), the inversion symmetry is broken (cf.Results: Bulk modes.) and the fourfold degeneracy is lifted.As a result, the bulk modes are featured with a Weyl cone with twofold degeneracy at the Weyl point, along with a cylindrical surface [cf.Eq. ( 6)].The topological charge carried by the Weyl point is consistent with the nonzero topological invariants C ± = ±2 of the present system (cf.Results: Topological invariants.).Here, the topological charge ±2 is associated with the unconventional spin-1 Weyl point with threefold linear degeneracy 46,[72][73][74] .The net chirality vanishes in the Weyl semimetal, agreeing with the fact that the total Chern number is zero (cf.Results: Topological invariants.).Note also that there is only one Weyl point in the system, which is similar to the case of chiral metamaterial with C 3 rotational symmetry 46 .The feature of single Weyl point can be observed in the recent experimental study, where the Weyl point is surrounded by charged nodal walls 75 .There, however, may exist two 42,53 or four 44,55 Weyl points, depending on the constitutive relations of the metamaterials.A number of different Weyl points can be found in gyroid photonic crystals 40 and non-symmorphic phononic crystals 74 .Very recently, higher-dimensional topological systems are proposed and expected to possess properties that their lower-dimensional counterparts do not support, such as Weyl surfaces and Weyl arcs 76,77 .
Figure 3b shows the dispersion of surface modes at the interface between vacuum and the chiral metamaterial in the frequency-wave vector space.For comparison, the bulk modes at constant frequencies for k y = 0 are outlined in gray mesh and overlaid in the same plot.Different from the surface modes in topological insulators that exist in the frequency (energy) band gap, the surface modes in gapless topological semimetals are defined in the region free of bulk modes at the same frequency (energy) 22 .Because of the frequency dependence of material parameters, the dispersion of surface modes is shown to be bended.In particular, the surface modes form a pair of spiral surface sheets wrapping around the Weyl cone.This feature is similar, though not identical, to the  www.nature.com/scientificreports/helicoid edge states that occur in the topological semimetals 23,44,72,78 .Similar features of helicoid-like edge states have also been observed in chiral metamaterials 46 .Notice that the surface modes move to the other side of the k z axis as the frequency goes across the Weyl frequency ω 1 , where the dispersion of surface sheets experiences a smooth transition.At the Weyl frequency, the edge state that connects the Weyl point forms the so-called Fermi arc.The surface mode at the Weyl frequency is considered the Fermi arc-like edge state.In the present configuration, the two spiral surface sheets contact at the Weyl point and two associated Fermi arc-like edge states concatenate to yield a straight line segment.
In conclusion, we have investigated the photonic topological phases in chiral metamaterials characterized by the magnetoelectric tensors with diagonal chirality components.The underlying medium is considered a photonic analogue of the topological semimetal featured with a Weyl cone and a cylindrical surface in the frequencywave vector space.Surface modes at the interface between vacuum and the chiral metamaterial exist in their common gap in the wave vector space, which are analytically formulated by algebraic equations.In particular, the surface modes form a pair of spiral surface sheets wrapping around the Weyl cone, resembling the helicoid surface states that occur in the topological semimetals.

Methods
Spin-orbit Hamiltonians.The wave equation for the hybrid modes F ± = E ± iH ′ in Eq. ( 3) can be rewritten as where and ψ± = − T are pseudospin states that include ±π/2 phase difference between the transverse hybrid field components (with respect to the optical axis of the medium) 57 .In the neighborhood of a reference frequency ω ref , ε n ( n = t, z ) can be approximated as where εn is positive definite 57 .Taking into account of the frequency dispersion of the medium near the reference frequency, Eq. ( 18) is rearranged as a pair of eigensystems: where and ψ ± = U −1 ψ± with U = diag εz /ε t , 1, εz /ε t .In the isotropic case, where ε t0 = ε z0 ≡ ε and εt = εz ≡ ε , Eq. ( 21) is simplified to where and D ± = ±ω ref (ε ± γ /ε) .Here, v = c/ε , k = k x x + k y ŷ + k z ẑ , S = S x x + S y ŷ + S z ẑ , with being the spin matrices of spin 1, and * denoting the complex conjugate.where k x = |k| sin θ cos φ , k y = |k| sin θ sin φ , and k z = |k| cos θ have been used.Here, θ and φ are the polar and azimuthal angles, respectively, on the surface S: The eigensystem for the Hamiltonian H ± : is solved to give the eigenvalues σ ± = σ v|k| ( σ = ±1, 0 ) and the normalized eigenvectors as Note here that the eigenvalue σ ± is related to δω in Eq. ( 24) as σ ± = D ± ± δω .Based on Eqs. ( 29) and ( 30), the Berry connections A σ ± = −i ψ σ ± ∇ψ σ ± are obtained as The Berry curvatures F σ = ∇ × A σ ± are then given by Integrating over the sphere S, the Chern numbers C σ = 1 2π S F σ • ds are calculated to give Surface wave equation.According to Maxwell's equations, the eigenfields on either side of the interface ( y = 0 ) are given by the nontrivial of E and H [cf. Eq. ( 3)] or the null space of H m [cf.Eq. ( 14)].On vacuum side (say, y > 0 ), we have where k z is the normal (to the interface) wave vector component in vacuum.On the chiral medium side ( y < 0 ), the eigenfields are given by where n ± = √ εµ ± γ , and k (1) z are the normal wave vector components in the chiral medium.Note that the eigenfields in Eqs. ( 35)- (38) share the common tangential wave vector components k x and k z across the interface, as a direct consequence of the phase matching of electromag- netic fields.For the surface waves to exist on the vacuum side ( y > 0) , k y should be purely imaginary with a positive value, so that the waves decay exponentially away from the interface.On the chiral medium side ( y < 0 ), k , H (4) = i η 0 ε µ E (4) , where n = x, z and C 1 , C 2 , C 3 , C 4 are constants.The existence of a nontrivial solution of these constants requires that the determinant of the 4 × 4 matrix obtained from Eqs. ( 39) and (40) be zero, which gives the characteristic equation of surface modes as Electromagnetic energy density.The time-averaged energy density in a lossless medium is given by 79 where and V = ε 0 E x , ε 0 E y , ε 0 E z , µ 0 H x , µ 0 H y , µ 0 H z T , with V † being the Hermitian conjugate of V.The energy density must be positive definite, which implies that both the trace and determinant of M are positive: Based on the Lorentz dispersion model for the present medium (cf.Results: Photonic Weyl system.), these quantities become and both of which are positive in the present study. ( https://doi.org/10.1038/s41598-023-40926-8

yθ e ∓iφ sin θ 0 e ±iφ sin θ 0 e ∓iφ sin θ 0 e ±iφ sin θ − √ 2
should be purely imaginary with a negative value for a similar reason.The tangential electric and magnetic field components are continuous at the interface: