# Radio wave attenuation by a large-scale photonic crystal sculpture

David Röhlig, Vincent Laude, Ralf Zichner, Felix Thieme, Angela

Thränhardt, Thomas Blaudeck

► To cite this version:

David Röhlig, Vincent Laude, Ralf Zichner, Felix Thieme, Angela Thränhardt, et al. Radio wave attenuation by a large-scale photonic crystal sculpture. Scientific Reports, 2025, 15, pp.12317. 10.1038/s41598-025-95986-9. hal-05221418

HAL Id: hal-05221418

https://hal.science/hal-05221418v1

Submitted on 25 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Radio wave attenuation by a large-scale photonic crystal sculpture

David Röhlig<fnref n="1" />✉, Vincent Laude<fnref n="2" />, Ralf Zichner³,⁴, Felix Thieme³, Angela Thränhardt<fnref n="1" /> & Thomas Blaudeck⁴,⁵

In this work, we present an artwork that embodies the by far largest photonic crystal ever published, operating at radio frequencies. We demonstrate both theoretically and experimentally the existence of band gaps. While photonic crystals are typically investigated on the nano and micro scale, our aim is to present not only a remarkable example of art with scientific relevance but also to highlight the potential for large-scale applications that have so far been underrated.

**Keywords** Photonic crystal, Electromagnetic crystal, Band gap, FEM, Radio waves, Mobile communication

The splendor of a butterfly's wing, illuminated by the sun, unveils not only aesthetic magnificence but also the conceptual beauty of the underlying mechanism: the intricate chitin nanostructure of its scales interacts with light<fnref n="1" />⁻³, blurring the lines between science and art. These structures, if periodic, are commonly recognized as photonic crystals, first described in 1987 both by Yablonovitch⁴ and John⁵, playing pivotal roles not only in nature but also in diverse technological applications⁶. They consist of different dielectric materials or even metals⁷,⁸ arranged in a periodicity comparable to the wavelength of light.

Remarkably, these crystals are predominantly examined at the meso- or microscopic level, possibly due to human fascination for visible light, yet their potential applications extend beyond this scope. The situation is different in the realm of acoustics: here, the concept of phononic crystals has been coined as analogous to photonic crystals⁹,¹⁰. In particular, a publication exists for a comparable art installation made of a two-dimensional periodic structure of cylinders made from stainless steel, showcasing the formation of an acoustic band gap for the propagation of sound waves and hence the observation of a frequency-dependent sound attenuation¹¹.

For electromagnetic waves, it is uncommon to consider the concept of periodicity at such large scales. Around the turn of the millennium, there was a strong interest in photonic structures operating at radio frequencies, as emphasized by Yablonovitch and others¹²,¹³. Surprisingly, there still appears to be no literature addressing photonic crystals, or in this context sometimes referred to as electromagnetic crystals¹²⁻¹⁴, with lattice constants exceeding a few centimeters. They have been studied in the millimeter to centimeter range for one-dimensional multilayer structures¹⁵⁻¹⁷, as well as for two⁻¹⁸ and three-dimensional structures¹⁹⁻²¹, operating primarily in the radio wave regime.

In this article, we present an investigation of an unusual example of a two-dimensional metallic photonic crystal that far exceeds the scale of previously reported structures. Its periodicity is relevant to radio wave frequencies pertinent to network standards such as 4 or 5G²². The structure under study is the artwork “Denk- und Wahrnehmungsmodell zum Phänomen der Farbe” (translated as “Thinking and Perception Model on the Phenomenon of Colour”) from 1998 created by artist Stefan Nestler and located in front of the central lecture hall of Chemnitz University of Technology. This piece of art consists of 187 periodically arranged hollow steel columns, standing upright and in parallel. Interestingly, configurations involving rectangular cavities within a metallic waveguide are well-known to function as bandpass filters²³,²⁴. More notably, the two-dimensional structure of the stelae formally satisfies the criteria of a finite photonic crystal and therefore joins the ranks of art with scientific relevance. Our objective is to adopt a cross-scale perspective on photonic crystals, further unraveling their potential and application possibilities.

<footnote n="1">¹Institute of Physics, Chemnitz University of Technology, Reichenhainer Straße 70, 09126 Chemnitz, Germany.</footnote>

<footnote n="2">²Université Marie et Louis Pasteur, CNRS, FEMTO-ST Institute, 15B avenue des Montboucons, 25030 Besançon, France.
³Faculty of Electrical Engineering and Information Technology, Chemnitz University of Technology, Reichenhainer Straße 70, 09126 Chemnitz, Germany.
⁴Fraunhofer Institute for Electronic Nano Systems (ENAS), Technologie-Campus 3, 09126 Chemnitz, Germany.
⁵Research Center for Materials, Architectures and Integration of Nanomembranes (MAIN), Chemnitz University of Technology, Rosenbergstraße 6, 09126 Chemnitz, Germany.</footnote>

<footnote>✉email: david.roehlig@physik.tu-chemnitz.de</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Artwork “Denk- und Wahrnehmungsmodell zum Phänomen der Farbe” by Stefan Nestler from 1998, consisting of 187 colored steel stelae.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Visualization of the sculpture (generated with Blender ^25), viewed from above, consisting of a 14 × 14 grid of steel stelae, with nine missing elements along one of the primary diagonals. Two measurement scenarios are indicated by dashed lines, corresponding to $S^1$ and $S^2$.</figcaption></figure>

## Modeling the structure

The stelae are arranged on a 4 × 4 square meter metal plate, forming a 14 × 14 array with a lattice constant of 0.24 meters (see Fig. 1). Along one of the main diagonals, nine columns are missing, as illustrated in Fig. 2. The measured heights of the stelae range from 0.99 to 3.97 meters, each with a base area of 120 × 120 millimeters. They are hollow, and anchored in metal rails. Exposure to local weather conditions has resulted in minor surface inhomogeneities. The overall wall thickness is approximately 5 millimeters.

For describing wave phenomena within a photonic crystal, we begin with an eigenvalue equation, which can be directly derived from Maxwell's equations by assuming time harmonic field solutions $F(r, t) = F(r)e^{i\omega t}$. In two dimensions, the governing equation can be decoupled for two distinct modes, yielding:
$$ -\nabla \cdot \left( a(\mathbf{r}) \nabla F_z(\mathbf{r}) \right) = b(\mathbf{r}) \frac{\omega^2}{c^2} F_z(\mathbf{r}), \qquad (1) $$

where $\omega$ is the angular frequency and $c$ represents the speed of light in vacuum. This formulation applies to both TM modes, where the magnetic field is oriented along the $z$ direction (perpendicular to the plane of periodicity), and TE modes, where the electric field is out of plane. Depending on the mode, the parameters $a(\mathbf{r})$ and $b(\mathbf{r})$ are defined to be:

<table><thead><tr><th>mode</th><th>$F_z(\mathbf{r})$</th><th>$a(\mathbf{r})$</th><th>$b(\mathbf{r})$</th></tr></thead><tbody><tr><td>TM</td><td>$H_z(\mathbf{r})$</td><td>$\epsilon(\mathbf{r})^{-1}$</td><td>1</td></tr><tr><td>TE</td><td>$E_z(\mathbf{r})$</td><td>1</td><td>$\epsilon(\mathbf{r})$</td></tr></tbody></table>

## Band structure calculations

For calculating the band structures, we employed the finite element method (FEM), which necessitates a variational formulation. To derive it, we begin by multiplying the eigenvalue Eq. (1) with a test function $\tilde{F}_z$. Next, we integrate over the domains $\Omega$ within the unit cell and apply the divergence theorem. Given that, in the context of band structures, we are exclusively concerned with the propagation of Bloch waves within an infinite periodic structure, we incorporate Bloch wave solutions $F_z(\mathbf{r}, t) = F_z(\mathbf{r}) \exp(i\omega t - i\mathbf{k} \cdot \mathbf{r})$ for both the field component and the test function, so that we derive the following expression:

$$ \begin{aligned} 0 = & \frac{\omega^2}{c^2} \int_{\Omega} \tilde{F}_z b F_z \, d\mathbf{r} \\ & - \int_{\Omega} (\nabla + i\mathbf{k}) \tilde{F}_z \cdot [a(\nabla - i\mathbf{k})F_z] \, d\mathbf{r}. \end{aligned} \qquad (2) \\ & + \int_{\partial\Omega_i} \tilde{F}_z [a(\nabla - i\mathbf{k})F_z] \cdot \mathbf{n} \, ds \quad \forall \tilde{F}_z. $$

Functions $\tilde{F}_z$ and $F_z$ in this equation are understood as the periodic part of Bloch waves. We could, in principle, decompose the integrals over the distinct subregions $\Omega$ of the air matrix and the inclusion; however, as we shall see shortly in the next section, it suffices to mesh only the matrix. The final term represents integrals along the internal boundaries $\partial\Omega_i$, incorporating the normal vectors $\mathbf{n}$. At the external boundaries $\partial\Omega_e$, periodic boundary conditions are applied and lead to the cancelling of boundary terms along $\partial\Omega_e$.

## Boundary conditions

Although the scatterer material exhibits dispersion, we neglect the explicit frequency dependence by applying metallic boundary conditions, which, as will be demonstrated later, show good agreement with experimental results. Specifically, we assume an infinitely high conductivity ($\sigma \to \infty$) within the metal, such that the electric field vanishes at the material boundary $\partial\Omega_i$ ($E \to 0$). This is heuristically supported by Ohm's law, $j = \sigma E$, for a finite current density $j$. Given that the tangential component of the electric field must remain continuous across the material interface $\partial\Omega_i$, the condition $E_1 \times \mathbf{n} = E_2 \times \mathbf{n}$, where $\mathbf{n}$ denotes the normal vector, must hold. Assuming the second region behaves as a so called perfect electric conductor (PEC), the field $E_2$ is zero, resulting in the following condition at the boundary $\partial\Omega_i$:

$$ E_1 \times \mathbf{n} = 0. \qquad (3) $$

This considerably simplifies the problem. For the calculations, we can now treat the hollow steel inclusions as filled, since the fields can not enter the PEC. Moreover, there is no need to mesh its domain within the unit cell, as it can be entirely described by the boundary integral. Given that the matrix material is air, we set $\epsilon = 1$, which leads to $a = b = 1$ in the weak formulation, as shown in Eq. (2). Consequently, the weak form is identical for both mode types, except for the application of boundary conditions. For TE modes, we simply apply Dirichlet boundary conditions, forcing $F_z = E_z = 0$ along the boundary. For TM modes, we must consider Ampère's law in the time-harmonic regime, $\nabla \times \mathbf{H} = (\sigma - i\omega\epsilon)\mathbf{E}$. By inserting it into the PEC condition Eq. (3), we obtain:

$$ 0 = \mathbf{n} \cdot (\nabla - i\mathbf{k})H_z. \qquad (4) $$

In the weak form Eq. (2), this is a Neumann boundary condition. The boundary integral in Eq. (2) then naturally vanishes.

## Transmission properties

To investigate the transmission properties of the artwork, we return to the variational formulation (2). Since – unlike for band structures – we are only concerned with a finite structure, we adhere to the previous derivation but $F_z$ now represents the electric or magnetic field and not a Bloch wave. By prescribing a specific value on a defined line within the simulation cell, as indicated in Fig. 5, we can introduce a source. In this case, the left-hand side of the equation is no longer zero. Since for a plane wave, $\nabla F_z = -i\mathbf{k}_0 F_z$, we obtain:

$$ \tilde{F}_z (a \nabla F_z) \cdot \mathbf{n} = - \tilde{F}_z a F_z i \mathbf{k}_0 \cdot \mathbf{n} \qquad (5) $$
for the expressions inside the boundary integrals. To enable plane waves to freely exit the simulation cell, we impose radiation boundary conditions at the external cell borders $\partial\Omega_e$ in the x and y directions, putting $k_0 \cdot n = \omega/c$. Consequently, we arrive at:

$$ \int_{\text{source}} \tilde{F}_z F_0 \frac{i a \omega}{c} ds = - \frac{\omega^2}{c^2} \int_{\Omega} \tilde{F}_z b F_z dr + \int_{\Omega} \nabla \tilde{F}_z \cdot (a \nabla F_z) dr + \int_{\partial \Omega_e} \tilde{F}_z F_z \frac{i a \omega}{c} ds \quad \forall \tilde{F}_z. \quad (6) $$

Note that Dirichlet (for TE waves) or Neumann (for TM waves) boundary conditions are still imposed along the boundaries $\partial\Omega_i$ defining the inclusions. Here, $F_0$ denotes the initial field generated by the source. To ultimately determine the frequency response, we integrate along the line of the sensor position, as depicted in Fig. 5:

$$ T = \frac{1}{F_0 l} \int_{\text{sensor}} F_z \, ds. \quad (7) $$

Both the sensor and the source have a length $l$, corresponding to a dipole length in the experiment. For the measurements, we employ logarithmic antennas composed of numerous dipoles of varying sizes and distances (from the antenna base) – each responsible for different frequency ranges. In the simulation, these geometric aspects were taken into account. Consequently, depending on the frequency, both the position and size of the sensor regions vary according to the geometry of our antenna model. Since, for TE modes, the dipoles are oriented out of the plane of periodicity, in a two-dimensional simulation, this can only be considered for TM modes. For a more comprehensive understanding of the theoretical foundations discussed in this section – including radiation boundary conditions – the reader is referred to Ref. 26.

## Experiment

For the exploration of the properties of this finite photonic crystal, we focus on transmission spectra, as they are particularly suitable for comparing theoretical predictions with experimental results. The setup comprises a vector network analyser (model: N5245A from Agilent Technologies) and two logarithmic periodic broadband antennas (model: STLP 9149 from Schwarzbeck), connected via coaxial cables (model: elspec Duratest 150 Series from Teledyne Technologies). The antennas were mounted on tripods at a height of 1.06 meters, each positioned one meter from the outer stele row, with the antennas directed towards each other. Using the network analyzer, we measured the scattering parameters $S$ for two distinct orientations, as indicated by the dashed lines in Fig. 2, both corresponding to the X (or Y) direction in reciprocal space. The horizontal line is denoted by $S^1$, and the vertical line by $S^2$. Data were recorded from sensor one to sensor two ($S_{21}$) and the reverse ($S_{12}$), while each curve is derived by averaging 10 individual measurements. The frequency response function is then defined as the mean of the different scattering parameters, as follows:

$$ T(f) = \frac{1}{4} (S_{12}^1 + S_{21}^1 + S_{12}^2 + S_{21}^2). \quad (8) $$

The same definition is employed for the response spectrum, which was calculated using FEM.

## Results

The transmission characteristics of the periodic stele artwork, analyzed as a function of frequency, for both TM and TE modes, are shown in Fig. 3 and 4. Overall, the experimental measurements reveal pronounced frequency-dependent variations of the response function, aligning well with the simulated FEM spectra. This is consistent with the anticipated activity of a photonic crystal, as represented by the band structure. For both polarizations, we observe the occurrence of band gaps, within which waves cannot propagate due to the exclusive presence of evanescent solutions.

The band structure for TM modes, shows characteristic features such as directional Bragg gaps, opening at the high-symmetry points $\Gamma$ and $X$; they are evident as spectral dips in the response function. The measurements exhibit a strong agreement not only with the band structure but also with the simulated response. Notably, the cell of the latter closely resembles the experimental setup, as illustrated in Fig. 5 for the magnetic field.

For TE modes, the remarkable phenomenon of a cutoff band at 1 GHz, which can be interpreted as an effective plasma frequency of the periodic stelae system 27, leads to a significant omnidirectional response reduction at lower frequencies (as illustrated in Fig. 4). Around this band, a slight increase is recorded, which delineates the low-frequency band gap from the subsequent complete gap between 1 and 1.5 GHz. The increase in response is not much higher than in the gap regions, as coupling modes originating from the antenna with the flat-band modes within the artwork proves particularly challenging in the experiment – primarily due to the low group velocity. The second band is clearly indicated by a response maximum, followed by a dip that precisely corresponds to the third frequency gap.
<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Theoretical band structure predictions for TM polarization, obtained using FEM on a triangular mesh of P2 elements with a resolution of 25 nodes per lattice constant at both the external and internal boundaries, are compared both to the simulated and measured frequency response function of the stele artwork. The simulation cell, depicted in Fig. 5, features the same resolution as that used for the band structure calculation. Measurements were taken along the Γ-X path in reciprocal space under TM polarization, wherein the magnetic field is oriented perpendicular to the plane of periodicity. The spectra fall within a frequency range of 0.5 – 1.5 GHz.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Theoretical band structure predictions are compared with the measured and simulated frequency response function, analogous to Fig. 3, but for TE polarization, where the electric field is oriented perpendicular to the plane of periodicity.</figcaption></figure>

Regarding the overall influence of defects, it was found that the different directions are indistinguishable in the simulated spectra. This finding is corroborated by the experiment, which reveals almost no difference between forward and backward measurements (*S*<sub>12</sub> or *S*<sub>21</sub>) in the observed response.

In general, we can state that there is a consistent agreement between theory and experiment, especially when considering the numerous approximations and potential sources of error. These primarily include the following factors: the finiteness of the structure, the nine defects, the varying stele heights in the third dimension, and the fact that the frequency dependence of the permeability and permittivity function of the metal were not accounted for in the two-dimensional band structure calculations. Additionally, the scatterers are hollow and mounted on four interconnected metal ground plates. Due to the public nature of the measurement location, several external sources of influence were encountered, including: a 2 × 3 meter glass plate positioned adjacent to the setup, transient passers-by, bicycle parking stalls about four meters from the source, a parking area with vehicles approximately 200 meters behind. Furthermore, 5G equipment emitting a pulsed signal is installed on the roof of a neighboring building, at a considerable distance of approximately 150 meters – well outside the main radiation direction. The influence remains rather significant – despite the time-varying nature of the signal, no peaks could be observed on the network analyzer.
<figure><img src="image_7.png" /><figcaption>Fig. 5. The normalized field distribution for TM modes at a frequency of 1.25 GHz was computed using the FEM with P2 elements with a resolution of 25 nodes per lattice constant, by applying radiation boundary conditions in both the x- and y-directions. This simulation cell serves as the basis for calculating the frequency response curves shown in Fig. 3 and 4. The source is defined on the left, while the frequency response is evaluated on the right, corresponding to the measurement scenario $S_{21}^1$.</figcaption></figure>

## Conclusions

With this article, we have demonstrated that art can convey additional meanings beyond the artist's original intent. The sculpture by Stefan Nestler subjects radio waves to the principles of periodicity, which are confirmed by both photonic band structure and frequency response calculations that align well with the experiment. Our findings reveal that within the steele structure, there are frequency regions where electromagnetic wave propagation is confined to evanescent modes. With these results, we aim to draw attention to the potential of applying the concept of photonic crystals on scales much larger than those previously reported, thereby opening new avenues for specific applications in filtering, guiding, or localizing waves that fall within the telecommunications spectrum.

## Data availability

All numerical and experimentally obtained data are provided within the supplementary information files.

Received: 15 November 2024; Accepted: 25 March 2025

Published online: 10 April 2025

## References

1. Biró, L. et al. Living photonic crystals: Butterfly scales — Nanostructure and optical properties Materials Science and Engineering: C **27**, 941 eMRS 2006 Symposium A: Current Trends in Nanoscience - from Materials to Applications. https://doi.org/10.1016/j.msec.2006.09.043 (2007).

2. Corkery, R. W. & Tyrode, E. C. On the colour of wing scales in butterflies: Iridescence and preferred orientation of single gyroid photonic crystals. *Interface Focus*. **7**, 20160154. https://doi.org/10.1098/rsfs.2016.0154 (2017).

3. Thomé, M., Richalot, E. & Berthier, S. Light guidance in photonic structures of Morpho butterfly wing scales. *Applied Physics A*. **126**, 778. https://doi.org/10.1007/s00339-020-03948-x (2020).

4. Yablonovitch, E. Inhibited spontaneous emission in solid-state physics and electronics. *Phys. Rev. Lett.* **58**, 2059. https://doi.org/10.1103/PhysRevLett.58.2059 (1987).

5. John, S. Strong localization of photons in certain disordered dielectric superlattices. *Phys. Rev. Lett.* **58**, 2486. https://doi.org/10.1103/PhysRevLett.58.2486 (1987).

6. Butt, M., Khonina, S. & Kazanskiy, N. Recent advances in photonic crystal optical devices: A review. *Optics & Laser Technology*. **142**, 107265. https://doi.org/10.1016/j.optlastec.2021.107265 (2021).

7. Scalora, M. et al. Transparent, metallo-dielectric, one-dimensional, photonic band-gap structures. *J. Appl. Phys.* **83**, 2377. https://doi.org/10.1063/1.366996 (1998).

8. Zamudio-Lara, A. et al. Characterization of metal-dielectric photonic crystals Optical Materials **29**, 60 https://doi.org/10.1016/j.optmat.2006.03.026 selected papers from the Second Topical Meeting on Nanostructured Materials and Nanotechnology NANOTECH-2005 (2006).

9. Kushwaha, M. S., Halevi, P., Dobrzynski, L. & Djafari-Rouhani, B. Acoustic band structure of periodic elastic composites. *Phys. Rev. Lett.* **71**, 2022. https://doi.org/10.1103/PhysRevLett.71.2022 (1993).

10. Vasileiadis, T. et al. Progress and perspectives on phononic crystals. *J. Applied Phys.* **129**, 160901. https://doi.org/10.1063/5.0042337 (2021).

11. Martínez-Sala, R. et al. Sound attenuation by sculpture. *Nature*. **378**, 241. https://doi.org/10.1038/378241a0 (1995).

12. Scherer, A., Doll, T., Yablonovitch, E., Everitt, H. & Higgins, J. Guest editorial: Electromagnetic crystal structures, design, synthesis, and applications. *J. Lightwave. Technol.* **17**, 1928. https://doi.org/10.1109/JLT.1999.1673018 (1999).

13. Yablonovitch, E. Photonic crystals: semiconductors of light. *Scientific American*. **285**, 47. https://doi.org/10.1038/scientificamerica n1201-46 (2001).
14. Ghanem, F., Delisle, G., Denidni, T. & Ghanem, K. A metallic wire electromagnetic crystal structure for radar applications. in *Ultra-Wideband, Short Pulse Electromagnetics*, Vol. 9. 413–421. https://doi.org/10.1007/978-0-387-77845-7_48 (Springer, 2010).

15. Aghajamali, A., Suthar, B., Nayak, C. & Srivastava, S. K. Bandgaps of microwave photonic crystals: Study of quasi-periodic metamaterial multilayers. *Mater. Sci. Eng.* **295**, 116547. https://doi.org/10.1016/j.mseb.2023.116547 (2023).

16. Darzi, B. & Ghadi, A. Radar absorbing photonic-crystal structure containing magnetic Co0.9Fe2.1O4 ferrite nanocomposite and dielectric materials. *Optik*. **268**, 169814. https://doi.org/10.1016/j.ijleo.2022.169814 (2022).

17. Mohamed, A. G., Elsayed, H. A., Mehaney, A., Aly, A. H. & Sabra, W. Transmittance properties of one-dimensional metamaterial nanocomposite photonic crystal in GHz range. *Sci. Rep.* **12**, 18331. https://doi.org/10.1038/s41598-022-21455-2 (2022).

18. Darthy, R. R., Venkateswaran, C., Subramanian, V., Ouyang, Z. & Yogesh, N. Accessing new avenues of photonic bandgaps using two-dimensional non-Moiré geometries. *Sci. Rep.* **13**, 17077. https://doi.org/10.1038/s41598-023-44385-z (2023).

19. Takahashi, S. et al. Microwave hinge states in a simple cubic lattice photonic crystal insulator. *Phys. Rev. B* **109**, 125304. https://doi.org/10.1103/PhysRevB.109.125304 (2024).

20. Lu, Z., Shi, S., Schuetz, C. A., Murakowski, J. A. & Prather, D. W. Three-dimensional photonic crystal flat lens by full 3D negative refraction. *Opt. Exp.* **13**, 5592. https://doi.org/10.1364/OPEX.13.005592 (2005).

21. Üpping, J., Miclea, P., Wehrspohn, R., Baumgarten, T. & Greulich-Weber, S. Direction-selective optical transmission of 3D fcc photonic crystals in the microwave regime. *Photonics and Nanostructures - Fundamentals and Applications* **8**, 102. special Issue PECS 8. https://doi.org/10.1016/j.photonics.2009.11.002 (2010).

22. Pi, Z. & Khan, F. An introduction to millimeter-wave mobile broadband systems. *IEEE Communications Magazine* **49**, 101. https://doi.org/10.1109/MCOM.2011.5783993 (2011).

23. Dahiya, A. & Deshwal, D. Design Fundamentals: Iris Waveguide Filters Versus Substrate Integrated Waveguide (SIW) Bandpass Filters. *Proceedings of International Conference on Artificial Intelligence and Applications. Advances in Intelligent Systems and Computing*, Vol. 1164. https://doi.org/10.1007/978-981-15-4992-2_19 (Springer, Singapore, 2021).

24. Aghanim, A., Lasri, R., Oulhaj, O. & Chekenbah, H. Design of a waveguide bandpass filter with inductive irises for a satellite application E3S Web of Conferences. **351** https://doi.org/10.1051/e3sconf/202235101059 (2022).

25. Community, B. O. Blender - a 3D modelling and rendering package, Blender Foundation, Stichting Blender Foundation, Amsterdam (2018), http://www.blender.org, Blender 4.0.2.

26. Laude, V. *Phononic Crystals: Artificial Crystals for Sonic, Acoustic, and Elastic Waves*. https://doi.org/10.1515/9783110641189 (De Gruyter, 2020).

27. Brand, S., Abram, R. A. & Kaliteevski, M. A. Complex photonic band structure and effective plasma frequency of a two-dimensional array of metal rods. *Phys. Rev. B* **75**, 035102. https://doi.org/10.1103/PhysRevB.75.035102 (2007).

## Acknowledgements

This work was partially funded by a Ph.D. scholarship awarded by the Free State of Saxony (D. R.). Support by the TUCculture2025 initiative of Chemnitz University of Technology, by the EIPHI Graduate school (contract “ANR-17-EURE-0002”) and by the Bourgogne-Franche-Comté Region is acknowledged. We thank Stefan Nestler for his approval to employ the artwork for both the modelling and the experiments in the scope of this manuscript.

## Author contributions

David Röhlig: Conceptualization, methodology (simulation), data curation, analysis, visualization, interpretation, writing (original draft), writing (editing), funding acquisition. Vincent Laude: Resources, methodology (simulation), validation, analysis, writing (review), funding acquisition. Felix Thieme: Methodology (experiment), analysis, writing (contribution). Ralf Zichner: Resources, methodology (experiment), validation, analysis, writing (review), supervision. Angela Thranhardt: Resources, methodology (simulation), validation, analysis, writing (review), funding acquisition, supervision. Thomas Blaudeck: Conceptualization, project administration, interpretation, writing (contribution), writing (review), funding acquisition.

## Funding

Open Access funding enabled and organized by Projekt DEAL.

## Declarations

### Competing interests

The authors declare no competing interests.

### Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-025-95986-9.

**Correspondence** and requests for materials should be addressed to D.R.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
**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/.

© The Author(s) 2025