Flatness and boundness of photonic drumhead surface state in a metallic lattice

Nodal chain (NC) semi-metals have the degeneracy of interlacing rings in their band structure in momentum space. With the projection of degenerate rings towards crystal boundaries, there is a special type of surface dispersion appearing at surface Brillouin zone and termed drumhead surface state (DSS). Previously, experimental investigations on photonic NC and DSS have been done on metallic photonic crystals at microwave frequencies. However, far-field detection of DSS and its coupling to radiative modes in free space have not been studied. In the work, we analyze the photonic DSS in a metallic lattice by angle-resolved far-field reflection measurement and numerical simulation at terahertz (THz) frequencies, and reveal its flatness and boundness in band structure, even in the radiation continuum. Particularly, the DSS band can be tuned being from negatively dispersive via flat to positively dispersive by a single surface parameter, and the DSS at Γ point in surface Brillouin zone is in fact a symmetry-protected bound state in the continuum. Our results might have some potential applications towards THz photonics.


Model and experiment
Our design is a 3D metallic PhC with sub-millimeter feature size, which is fabricated through 3D printing of polymer framework and subsequent chemical plating metallization. Shown in Fig. 1a, the simple cubic lattice structure has the lattice constant, a = 580 µm , and the size of the square rod, d = 130 µm , seeing the unit cell in Fig. 1d. The band structure of the metallic PhC is calculated in a simulation software (CST Microwave Studio), using periodic boundary condition along three dimensions and PEC approximation, and the results are plotted in Fig. 1b,c. In the band structure along high symmetric path in the Brillouin zone, seeing Fig. 1b, we identify three degeneracy points formed by the crossing between band 3 and 4 and being located at 0.341THz. In fact, they belong to a global geometry of NC in the 3D Brillouin zone, seeing Fig. 1c, where three big rings are centered around X points and three small rings around M points considering the periodicity of Brillouin zone 21 . Because of slight frequency dispersion and limited numerical calculating precision, the NC presents a tiny variation in chain diameter, and is not as smooth as mathematical chain model. It is known that the crossing between the 3rd and 4th bands is enforced by mirror symmetry M z : {x, y, z} → {x, y, -z} 21 . Therefore, two eigenmodes of the degeneracy along XM path, illustrated in Fig. 1d, show even and odd parities of electric field (E-field) with respect to the z = 0 plane. The sample is made by two steps of processes. First, a simple cubic lattice of photopolymer framework is produced via 3D printing with printing precision ± 2.5 μm. Then, the polymer structure is coated with a thin layer of nickel through chemical nickel plating, i.e., electroless nickel plating. The resultant nickel layer is ~ 3 μm in thickness, far exceeding the skin depth at experimental frequencies, and effectively converts the framework into a metallic PhC. The total size of the whole metallic lattice is about 30 mm × 30 mm × 5 mm, and its top view along 〈001〉 direction is shown in Fig. 2a.
Due to bulk-edge correspondence in topological semimetals, the DSS is formed in the boundary of the lattice 21 . In Fig. 2b, we schematically plot the 〈001〉 projection of the NC and the resultant DSS at the [001] surface. In this work, we are concentrated on the DSS centered at Γ point in surface Brillouin zone, since it is  www.nature.com/scientificreports/ most readily accessible with far-field detection and is free of diffraction beam under small and moderate incident angles of illumination. Two angles are crucial to implement far-field measurement, being the incident angle θ and the rotating angle φ, as depicted in Fig. 2b. The latter defines the orientation of the incident plane with respect to the 〈001〉 crystal direction, while the former gives rise to the DSS-matching momentum k // = k 0 sinθ, where the wave vector in vacuum k 0 = ω/c with ω being angular frequency and c being the speed of light in vacuum. In experiment, we used a THz time-domain spectral system (Advantest TAS7400) to measure the reflection spectrum of the [001] surface of the sample. Being a pump-probe detecting system, the spectrometer generates and detects the pulsed THz waveform in time domain, and reflection is obtained from the Fourier-transformed signal recorded from the sample surface. The measuring setup is shown in Fig. 2c where the THz beam from the emitting module is obliquely incident upon the sample and the reflected beam is received by the receiving module. The sample is located at the sample stage, and during measurement we rotate it around the normal direction of the sample surface by the angle φ. When ϕ = 0 • , the incident plane is along the ГX direction. Because of the factory-preset p-wave polarization and the modules' bulky volume in the spectrometer, only R pp with θ > 15° can be measured, and the measuring results are plotted in Fig. 2d1-d4, where a flat (θ-independent and φ-independent) reflection dip appears at 0.34THz for four values of φ. It is noted that the dip is located in the region between two bulk band projections. The dip is caused by the DSS, and its angle-independence reveals non-dispersion or flatness of the DSS.
In order to verify the measuring results, we conduct an angle-resolved reflection simulation, where the finite number of unit cells in z direction are modeled and the periodic boundary condition is imposed only in x and y directions. In order to account for the material loss in the metallic PhC and meanwhile minimize the computing time, we set the finite conductivity for nickel through surface impedance parameter. The simulated R pp results are shown in Fig. 2d1-d4, and it is seen that both simulation and experiment are in good agreement after considering a fabrication tolerance of ~ 5 μm. We also simulate R ss where s-wave with E-field perpendicular to the incident plane is incoming and s-wave is outgoing. In sharp contrast to R pp , there is no reflecting dip appearing, as shown in Fig. 3b. Furthermore, R ps and R sp are calculated, and both show nearly null feature in the measuring frequencies as both polarizations are orthogonal to each other. Therefore, the sole sensitivity of the DSS dip to  Fig. 1c) and the drumhead surface state (pink area) in k x -k y plane, where yellow arrows denote far-field polarization of the drumhead surface state and only one quadrant of polarization texture is illustrated due to C 4 symmetry. The incident wave vector (orange arrow) and two incident planes (blue and yellow planes) are also depicted for angle-resolved far-field reflection measurement. θ is the incident angle, and φ is the rotating angle of incident plane with respect to 〈001〉 direction. (c) Photo of the reflection measurement. (d1-d4) Measured and simulated far-field reflection, R pp , with E-field polarized in the incident plane for both incoming and outgoing waves (p-waves). In each panel, left diagrams with θ ranging from 0° to 30° are simulated results, and right ones with θ ranging from 15° to 30° are measured results. Both comparisons are labeled by dashed regions and blue arrows. For simulation to match well with measurement, the size, d, of the metallic rods has been fine-tuned being 125 μm in numerical calculation.  www.nature.com/scientificreports/ p-wave polarization points to a monopole configuration of polarization vectors with monopolar center being Γ point, as illustrated schematically by double arrows in Fig. 2b.

Analysis and discussion
Flatness of photonic drumhead surface state. It has been revealed that the dispersion of surface state in electronic systems is sensitive to the surface decoration, even leading to Lifshitz transition of surface states in some topological semimetals 29,30 . Here, the dispersion of the photonic DSS can also be tuned significantly by a simple surface characteristic, as depicted in Fig. 3d1, that is the extruding length, h, of the sectioned metallic rod on the surface. It is seen from both measured and simulated results in Figs. 2 and 3a that the DSS is a flat (zero-dispersive) band being located at ~ 0.34THz. The flatness takes place for the case of h = 228 μm. Figure 3c shows the effect of the length on the DSS frequency, where the DSS dip shifts towards lower frequencies when increasing h. When h is extended above 315 μm, the dip begins to fuse into the first passband lying in ~ 0.22-0.3THz, representing DSS red-shifts into the bulk band. In contrast, when h is 185 μm, the dip has blue-shifted into the second passband starting from ~ 0.36THz. Figure 3e, Such surface property can be qualitatively explained by a LC circuit model, seeing the inset in Fig. 3c. The capacitance, C, comes essentially from the charges accumulated in the end of the rod, and the inductance, L, is mainly induced by the current flowing along the sectioned rod, sensitive to the length h, seeing the illustration of EM field in Fig. 3d1 31,32 . Therefore, the DSS frequency near Γ point is ω DSS = 1 √ LC . On the other hand, ω DSS far away from Γ point is pinned to the NC frequency ω NC which is determined by bulk geometry of the lattice and is robust to h. Therefore, h can tune ω DSS being either lower than or equal to or higher than ω NC , giving rise to either positively or zero-or negatively dispersive band. The simple physical picture is supported by the simulation, and Fig. 3d2 shows the simulated contour plot of E-field at 0.34THz for R pp with θ = 15°, φ = 0°, and h = 228 μm. The DSS appears clearly on the front [001] surface, upon which the incident wave is impinging, and the field configuration is consistent with the schematic drawing in Fig. 3d1.

Boundness of photonic drumhead surface state.
Another feature of the DSS in simulated and measured results is its vanishing behavior as approaching to Γ point. In other words, it is decoupling from the far field. In fact, the DSS at Γ point is the BIC, where even parity of the mode is incompatible with odd parity of free-space radiative modes, and it belongs to the symmetry-protected BICs. It has been revealed that appearance of BIC is accompanied by a singular point in polarization map 26,27 . In the DSS band, such singularity is seen at Γ point, as illustrated in Fig. 2b, and the monopole configuration of the polarization texture also implies a topological charge of 1.
One technical merit of BIC is extraordinarily high quality (Q-) factor which goes to infinity in ideal case. In order to evaluate the Q-factor, we employ the standard temporal coupled mode theory 4 , where the DSS lying in the bulk bandgap is regarded as single-mode resonator and couples its energy to the radiative mode (p-polarized far field) through reflection, Here, r denotes reflection coefficient, i 2 = −1 , τ r and τ nr characterize energy loss of the resonator due to radiative and non-radiative channel, respectively. τ r , τ nr , and ω DSS are numerically obtained through fitting the simulated angular reflection spectra by Eq. (1), and then the radiative Q-factor Q r = ω DSS τ r /2 and the nonradiative Q-factor Q nr = ω DSS τ nr /2 are calculated. Figure 4a,b shows, respectively, the Q r and Q nr map obtained from the reflection around Γ point where the radial direction away from the center of the map represents the angle θ ranging from 2° to the angular position of the NC projection. Note in Fig. 4a,b that the Q r and Q nr data is absent within the small green circles where θ is less than 2° and the DSS dip in R pp becomes too weak to be fitted accurately. It is seen that compared with the low Q nr , Q r is relatively high, especially near the center region. In particular, Q r rises rapidly when θ decreases, as plotted in Fig. 4c. Such characteristic of Q r manifests the BIC at Γ point, and can also be understood from aforementioned LC picture. The open LC circuit, i.e., sectioned metallic rod, radiates like a dipole antenna at frequency ω DSS , and thus the [001] surface is just an antenna array with a square lattice of z-orientated dipoles. According to the antenna theory, the far-field radiation pattern of the array becomes that of the individual element when the polar angle θ is zero and near zero, because the array factor is equal to constant 33 . According to the power distribution factor, sin 2 θ , of an individual dipole antenna, the far field is null at θ = 0°, and thus the DSS becomes the BIC. A simple comparison between the numerical Q r and the function plot 1/ sin 2 θ displays a qualitative agreement, as shown in Fig. 4c.

Conclusion
In conclusion, we have studies the DSS of the metallic simple cubic lattice at THz frequencies by angle-resolved far-field reflection measurement and numerical simulation. In particular, the DSS persists in the radiation continuum, displays the flatness property in band dispersion, and manifests a symmetry-protected BIC at the center of surface Brillouin zone, which are not identified in previous study. Such properties are of great significance to manipulation of THz wave in planar surfaces and technical applications like THz sensing.