Periodic Subwavelength Research Articles

Interference mechanism of plasma self-organization in transparent dielectrics under the intense femtosecond laser pulse exposure

A new mechanism of plasma self-organization in transparent dielectrics with wide bandgap exposed to the intense tightly focused laser radiation was revealed, which causes the generation of 3D periodic ring structures with subwavelength period both along the laser pulse propagation and in the radial direction. The mechanism involves formation of dense plasma burst in the pre-focal region that provides efficient scattering of the incident wave. The interference of a plane incident laser wave in the focal region and a divergent reflected one will form the standing wave pattern with local minima and maxima of laser field both in the direction of the incident wave propagation and perpendicular to it producing the ring patterns of effective ionization regions in the dielectric volume. Analytical and numerical simulations of the process of laser wave scattering on a near-spherical plasma object with dimensions both smaller and larger than the laser radiation wavelength are performed to verify the proposed model.

Ultrahigh reflectivity photonic crystal membranes with optimal geometry

Photonic crystal (PhC) structures with subwavelength periods are widely used for diffractive optics, including high reflectivity membranes with nanoscale thickness. Here, we report on a design procedure for 2D PhC silicon nitride membrane mirrors providing optimal crystal geometry using simulation results obtained with rigorous coupled-wave analysis. The Downhill Simplex algorithm, a robust numerical approach to finding local extrema of a function of multiple variables, is used to optimize the period and hole radius of PhCs with both hexagonal and square lattices, as the membrane thickness is varied. Following these design principles, nanofabricated PhC membranes made from silicon nitride have been used as input couplers for an optical cavity, resulting in a maximum cavity finesse of 33 000, corresponding to a reflectivity of 0.999 82. The role played by the spot size of the cavity mode on the PhC was investigated, demonstrating the existence of an optimal spot size that agrees well with predictions. We find that, compared to the square lattice, the hexagonal lattice exhibits a spectrally wider reflective range, less sensitivity to fabrication tolerances, and higher reflectivity for membranes thinner than 200 nm, which may be advantageous in cavity optomechanical experiments. Finally, we find that all of the cavities that we have constructed exhibit well-resolved polarization mode splitting, which we expect is due primarily to a small amount of anisotropic stress in the silicon nitride and PhC asymmetry arising during fabrication.

Interaction of highly focused intense laser radiation with transparent dielectrics: Subwavelength nanogratings formation

Femtosecond laser writing of birefringent subwavelength nanolattices in dielectrics has been studied for almost two decades since it reveals a number of applications for optical memory devices, optical waveguides, microfluidic channels, etc. In this work, a numerical study of the formation of plasma quasiperiodic nanostructures in a fused silica in the propagation direction of a focused laser pulse is carried out. It is shown that the focused beam creates a dense plasma, which provides an effective reflection of the incident laser pulse, leading to the formation of a standing wave. In the bundles of standing wave, an effective ionization emerges, which forms plasma gratings with a subwavelength period. The conducted modeling allows us to determine the conditions under which the proposed regime of material nanostructuring is possible. It is shown that the energy absorbed in plasma gratings will ensure the substance melting, which leads to the formation of “frozen” quasiperiodic nanopatterns.

Compact unidirectional waveguide grating emitter with enhanced wavelength sensitivity based on the hybrid plasmonic mode

Waveguide grating antennas are widely adopted in beam-steering devices, typically enabling the beam steering in longitudinal direction within a two-dimensional scanning optical array by changing the input wavelength. However, traditional waveguide grating antennas suffer from limited tuning range due to low dispersion of the gratings. In this paper, a compact silicon grating waveguide antenna array is proposed with enhanced wavelength sensitivity by introducing a periodically modulated hybrid plasmonic mode. The hybrid plasmonic mode is supported by the hybrid plasmonic waveguides (HPWs) composed of silicon waveguides and periodic subwavelength silver strips. In order to convert the guided waves to the radiated waves, a series of silicon emitting segments are deposited above the HPWs. Additionally, the horizontally arranged array of HPWs also acts as a reflector of the downward radiation, resulting in an effective unidirectional emission. Through the optimization of physical parameters, the proposed antenna array achieves a wavelength-length tuning efficiency up to 0.3°/nm within the wavelength range of 1500∼1600 nm, exhibiting a significant improvement compared with traditional ones. Moreover, an average upward emissivity exceeding 80% with a maximum value of 89% within the 100 nm bandwidth is demonstrated through the numerical simulations. The proposed compact antenna array provides an alternative solution in realizing large-scale integrated high-tuning-efficiency optical beam-steering devices.

Comprehensive review of the directed design of optomechanical crystal cavities using intelligent algorithms

Optomechanical crystal cavities are devices based on optomechanical interactions to manipulate photons and phonons on periodic subwavelength structures, enabling precise measurement of the force and displacement. The performance of the target structures varies when applied to different applications. Optomechanical crystal cavities now rely on an empirical forward design, which is inefficient. Therefore, a desired shift is toward directed design with a “problem-oriented” strategy. The directed optimization problem’s nonconvex nature and extensive parameter space necessitate substantial computational resources, driving the need for intelligent algorithms in a sub-wavelength structure design. Intelligent algorithms can surpass the constraints of traditional methods and discover novel structures that are effective in different materials, topologies, modes, and wavelengths. This paper provides an extensive overview of intelligent algorithms for guiding the directed design of optomechanical crystal cavities. It presents a systematic classification of 15 algorithmics, including, but not limited to, topology algorithms, particle swarm optimization algorithms, convolutional neural networks, and generative adversarial networks. The article provides a comprehensive review and thorough analysis of the principle and current application state, as well as the advantages and disadvantages of each intelligent algorithm. By using these intelligent algorithms, researchers can enhance the efficiency and accuracy of optimizing optomechanical crystal cavities in a broader design space.

Stacked Polarizing Elements for Controlling Parameters of Surface Relief Gratings Written in Photosensitive Materials.

Photosensitive materials are widely used for the direct fabrication of surface relief gratings (SRGs) without the selective etching of the material. It is known that the interferometric approach makes it possible to fabricate SRGs with submicron and even subwavelength periods. However, to change the period of the written SRGs, it is necessary to change the convergence angle, shift a sample, and readjust the interferometric setup. Recently, it was shown that structured laser beams with predetermined, periodically modulated polarization distributions can also be used to fabricate SRGs. A structured laser beam with the desired polarization distribution can be formed with just one polarizing optical element-for example, the so-called depolarizer, a patterned micro-retarder array. The use of such stacked elements makes it possible to directly control the modulation period of the polarization of the generated laser beam. We show that this approach allows one to fabricate SRGs with submicron periods. Moreover, the addition of q-plates, elements effectively used to generate cylindrical vector beams with polarization singularities, allows the efficient formation of fork polarization gratings (FPGs) and the fabrication of higher-order fork-shaped SRGs. Full control of the parameters of the generated FPGs is possible. We demonstrate the formation of FPGs of higher orders (up to 12) by only adding first- and second-order q-plates and half-wave plates to the depolarizers. In this work, we numerically and experimentally study the parameters of various types of SRGs formed using these stacked polarizing elements and show the significant potential of this method for the laser processing of photosensitive materials, which often also serve as polarization sensors.

Innovative Application of Subwavelength Periodic Polystyrene Microspheres as Saturable Absorbers in Nonlinear Optics

Polystyrene (PS) possesses numerous remarkable properties like high transparency, impressive mechanical strength, and a large specific surface area, making it an excellent mask plate or template for the preparation of anti-opal structures. Moreover, it should be noted that PS also exhibits exceptional nonlinear properties due to the subwavelength periodic configuration. In this paper, a self-assembled PS microsphere photonic crystal saturable absorber (PSM-SA) has been proposed and fabricated. It exhibits impressive properties including high stability, high damage threshold, high refractive index, and large specific surface area. It is suggested that the periodic structure of PS in the film has a significant impact on the photonic band gap, resulting in excellent adjustable optical nonlinear characteristics. By integrating PSM-SA into a self-built ring fiber laser system, a Q-switched laser with a pulse width of approximately 2 μs and a repetition rate of 40 kHz at a wavelength of 1562 nm is obtained. These findings demonstrate its potential for enabling efficient and adjustable nonlinear optical functionalities in various optical devices, contributing to the expanding realm of PS microsphere photonic crystals and their significant impact on advancing nonlinear optics technology.

2D compound structures with deep subwavelength period on silicon fabricated by double time delayed femtosecond laser beams

In this work, two orthogonally polarized femtosecond laser beams are employed to irradiate a p-doped silicon wafer with an electrical resistivity of 0.008 Ω·cm. It is interesting to find that 2D compound structures composed of sub-wavelength periodic ripples and deep sub-wavelength nanodot array can be produced when proper laser fluence and time delay between the dual laser beams are used. The formation of the periodic ripples can be explained by the interference between the preceding incident laser and it induced surface plasmon polaritons (SPPs). The periodic nanodot array has a period down to ∼200 nm and the radius of the nanodot is ∼30 nm, most of which appear at the boundary between the ditch and ridge of the ripple. During the ripples’ formation, the residual melting silicon is most probably located at the boundary between the ditch and ridge of the ripple. Furthermore, the period of the nanodot array is roughly equal to the perimeter of the nanodot. Therefore, it is considered that the dot array may be generated due to the Rayleigh-Taylor instability of the melting silicon. It is also noted that these nanodots are all uniformly arranged along vertical lines, indicating that the subsequent incident laser may break the stochastic characteristic of the Rayleigh-Taylor instability and produce the 2D periodic dot array. The thermo-hydrodynamical process combined with the interference effect between SPPs and the incident laser can benefit the formation of complex surface structures with versatile functions.

A polarization‐insensitive ultra‐wideband absorber based on hybrid structure

AbstractWith the continuous improvement of detection technology, the survivability of targets on the battlefield is facing great challenges. Consequently, microwave absorbers with stealth characteristics are of important research value to modern military. Aim at solving the problems of narrow bandwidth or complex structure, an implementation model of ultra‐wideband (UWB) absorber based on hybrid structure is proposed. In this letter, a polarization‐insensitive ultra‐wideband absorber (PIUWA) is designed, manufactured, and characterized, which has a good absorptivity of more than 90% ranges from 4 to 24.53GHz (the fractional bandwidth [FBW] is 144%) by triggering multiple resonance peaks in the hybrid structure. In addition, the subwavelength periodicity of the PIUWA is theoretically beneficial to the angular stability under full‐wave polarizations. The operating mechanism of the PIUWA is demonstrated by the equivalent circuit model and the correctness of the design is verified by experimental measurements. This letter contributes to the approach of designing and fabricating UWB microwave absorbers with high absorptivity, good angular stability, and simpler assembly procedure.

Design of Alvarez Beam Scanning Reflectarray With Inversely Proportional Focal Length

In this letter, a beam scanning Alvarez reflectarray (RA) metasurface with its focal length ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">F</i> ) inversely proportional to the phase distribution on the Alvarez RA aperture is presented for ±45 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">o</sup> beam scanning coverage and less than 3-dB gain loss. The beam scanning mechanism is based on mechanically tilting the Alvarez RA panel to the desired direction to steer the beam. Compared to a conventional RA where <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">F</i> is directly proportional to the parabolic phase distribution, the proposed Alvarez RA shows a wider 1-dB and 3-dB gain bandwidth, higher aperture efficiency (AE <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">)</i> , and wider scanning range with less than 3-dB gain loss. The proposed RA has adopted an Alvarez phase distribution which is a summation of two shifted cubic phase profiles. Pancharatnam- Berry (PB) meta-atoms of sub-wavelength periodicity were utilized to provide the 360 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">o</sup> phase range and to enhance the bandwidth of the RA. An Alvarez RA consisting of 38 × 38 PB meta-atoms has been designed, fabricated, and tested for ±45 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">o</sup> beam scan coverage with less than 3-dB gain loss. The simulated and measured results show that the gain of the proposed Alvarez RA metasurface is 29.1 dBi at 19.5 GHz, which corresponds to an AE of 44.2%. The 3-dB gain variation of the Alvarez RA is from 15 GHz to 23.4 GHz, which corresponds to a 3-dB gain bandwidth of 43.7%. The 1-dB gain bandwidth was 27.8% from 17.3 GHz to 22.9 GHz.

Realization of Single-Layer Fourier Phased Metasurfaces for Wideband RCS Reduction

An efficient and fast strategy to design and realize single layer Fourier phased metasurfaces for wideband radar cross section (RCS) reduction when illuminated by a circular polarization (CP) plane wave is proposed in this letter. The scattering phase (between 0 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">o</sup> and 360 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">o</sup> ) required at each unit cell of the proposed metasurfaces was computed using the Fourier phase formula in which the focal length ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">F</i> ) is inversely proportional to the phase distribution. Pancharatnam-Berry (PB) phase theory was applied with unit cells of subwavelength periodicity to further enhance the scattering and RCS reduction characteristics. The proposed wideband Fourier phased metasurface has a square shape and contains 30 × 30 PB unit cells with subwavelength periodicity of 5 mm ≈ 0.26λ <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">16GHz</sub> . Both simulation and measured results show that the proposed Fourier phased metasurfaces can achieve more than 10 dB of RCS reduction under normal incidence of CP plane wave regardless of the value of <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">F</i> . Under oblique incidence, more than 10 dB of RCS reduction was maintained for incident angles up to 60°. In addition, the single-layer Fourier phased metasurface features wideband 10 dB RCS reduction bandwidth from 10 GHz to 24 GHz with a thickness of only 2 mm. This resulted in an 82.3% fractional bandwidth (FBW) of RCS reduction which is higher than other designs reported in the literature. The proposed design strategy provides a promising way to design and realize metasurfaces for wideband and stable RCS reduction performance without the need to use a computationally complex and/or time-consuming and slow running optimization algorithm.

Numerical study of feature-distribution effects foranti-reflection structured surfaces on binary gratings.

Suppressing Fresnel reflections from dielectric boundaries using periodic and random antireflection structured surfaces (ARSSs) has been vigorously studied as an alternative to thin film coatings for high-power laser applications. A starting point in the design of ARSS profiles is effective medium theory (EMT), approximating the ARSS layer with a thin film of a specific effective permittivity, which has features with subwavelength transverse-scale dimensions, independent of their relative mutual positions or distributions. Using rigorous coupled-wave analysis, we studied the effects of various pseudo-random deterministic transverse feature distributions of ARSS on diffractive surfaces, analyzing the combined performance of the quarter-wave height nanoscale features, superimposed on a binary 50% duty cycle grating. Various distribution designs were investigated at 633nm wavelength for TE and TM polarization states at normal incidence, comparable to EMT fill fractions for a fused silica substrate in air. The results show differences in performance between ARSS transverse feature distributions, exhibiting better overall performance for subwavelength and near-wavelength scaled unit cell periodicities with short auto-correlation lengths, in comparison to equivalent effective permittivity designs that have less complicated profiles. We conclude that structured layers of quarter-wavelength depth and specific feature distributions can outperform conventional periodic subwavelength gratings as antireflection treatments on diffractive optical components.

Improving transmittance of long-wave infrared guided-mode resonant filters

The long-wave infrared (LWIR) spectral region spanning from 8 to 12 μm is useful for many scientific and industrial applications. Many of these applications require use of either a bandpass or a bandstop filter that can be realized by the guided-mode resonance (GMR) effect with subwavelength periodic features in layered dielectric materials transparent in the LWIR. The GMR filters operating in the LWIR region are fabricated by depositing an amorphous germanium (Ge) film to form a zero-contrast (ZC) waveguide-grating (WGG) on a polished zinc selenide (ZnSe) substrate. In general, the backside of a ZnSe substrate with refractive index 2.41 is uncoated causing a 17% Fresnel-reflection loss in the light transmitted through the filter due to a large impedance mismatch at the ZnSe/air interface. Because we use such filters in the LWIR laser experiments for notch filtering, to improve the filter transmittance we used ZnSe substrates coated on one-side with broadband antireflection coating (ARC) covering the 7 to 12 μm spectral range to fabricate GMRFs with one-dimensional (1D) Ge ZC WGG. We employed high-spatial resolution e-beam lithography and reactive-ion etching nanofabrication techniques to achieve high-performance large-area (12 × 12 mm2) 1D notch filters with subwavelength periods. We characterized polarization dependent spectral performance of the prototype filters with both coherent and incoherent incident light using a tunable quantum cascade laser system that spans the 7 to 12 μm region, and a Fourier transform infrared spectrometer with collimated incident beam to achieve close to 15% improvement in the peak transmittance as well as significant reduction in coherent noise compared to our earlier results with GMRFs without ARC. Here, we present the filter design simulation and measurement results.

Numerical investigation of sound transmission through single and double walls with periodic arrays of resonators

Lightweight partition walls are widely used in buildings, especially nowadays when sustainability is of a major interest. At the same time noise pollution is becoming a growing problem across the globe. All this implies the need to find possible solutions for sound insulation improvement of the lightweight partitions. This paper is focused on exploring the possibilities of exploiting the dynamic effects of periodic subwavelength arrays of resonators for sound insulation enhancement. In order to improve the sound transmission loss in a specific frequency region, periodic resonant units are tuned and introduced to single and double panels. It can be shown that the resonance mode of the resonators couples with the host plate vibration in the way of breaking the mass law and overcoming some phenomena like coincidence effect and mass-air-mass resonance. The investigation of the considered panels under acoustic field excitation is conducted using different methods. Sound transmission loss curves and dispersion diagrams are calculated. Finite element method is employed for this purpose, but also the theory of three-dimensional elasticity and well-known transfer matrix method (TMM). Several designs of resonators are proposed.

Genetic optimization of Brillouin scattering gain in subwavelength-structured silicon membrane waveguides

On-chip Brillouin optomechanics has great potential for applications in communications, sensing, and quantum technologies. Tight confinement of near-infrared photons and gigahertz phonons in integrated waveguides remains a key challenge to achieving strong on-chip Brillouin gain. Here, we propose a new strategy to harness Brillouin gain in silicon waveguides, based on the combination of genetic algorithm optimization and periodic subwavelength structuration to engineer photonic and phononic modes simultaneously. The proposed geometry is composed of a waveguide core and a lattice of anchoring arms with a subwavelength period requiring a single etch step. The waveguide geometry is optimized to maximize the Brillouin gain using a multi-physics genetic algorithm. Our simulation results predict a remarkable Brillouin gain exceeding 3300 W−1m−1, for a mechanical frequency near 15 GHz.

High performance antenna-on-chip inspired by SIW and metasurface technologies for THz band operation

In this paper, a high-performance antenna-on-chip (AoC) is implemented on gallium arsenide (GaAs) wafer based on the substrate integrated waveguide (SIW) and metasurface (MTS) technologies for terahertz band applications. The proposed antenna is constructed using five stacked layers comprising metal-GaAs-metal-GaAs-metal. The conductive electromagnetic radiators are implemented on the upper side of the top GaAs layer, which has a metallic ground-plane at its underside. The metallic feedline is implemented at the underside of the bottom GaAs layer. Dual wrench-shaped radiators are framed by metallic vias connected to the ground-plane to create SIW cavity. This technique mitigates the surface waves and the substrate losses, thereby improving the antenna’s radiation characteristics. The antenna is excited by a T-shaped feedline implemented on the underside of the bottom GaAs substrate layer. Electromagnetic (EM) energy from the feedline is coupled to the radiating elements through the circular and linear slots etched in the middle ground-plane layer. To mitigate the surface-wave interactions and the substrate losses in the bottom GaAs layer, the feedline is contained inside a SIW cavity. To enhance the antenna’s performance, the radiators are transformed into a metamaterial-inspired surface (i.e., metasurface), by engraving periodic arrangement of circular slots of sub-wavelength diameter and periodicity. Essentially, the slots act as resonant scatterers, which control the EM response of the surface. The antenna of dimensions of 400 × 400 × 8 μm3 is demonstrated to operate over a wide frequency range from 0.445 to 0.470 THz having a bandwidth of 25 GHz with an average return-loss of − 27 dB. The measured average gain and radiation efficiency are 4.6 dBi and 74%, respectively. These results make the proposed antenna suitable for AoC terahertz applications.

Topology-optimized ventilation barrier for mid-to-high frequency ultrabroadband sound insulation

An ultrabroadband ventilated acoustic barrier, composed of periodic subwavelength units (with a thickness of 50 mm), is put forward here by a density-based topology optimization method. We propose a unit with a cross channel as the initial unit for the topology optimization process. Compared with the initial one, the optimized barrier reduces the unit volume filling ratio by nearly 50 %. The efficient sound attenuation band is shifted to a lower frequency by nearly 665 Hz, and the broadband is broadened by nearly 870 Hz. Then, a simplification of the optimized unit is carried out, and the influence of the structure parameters on the acoustic attenuation performance is discussed. Numerical simulations and experiments are conducted to verify the sound-blocking performance. The optimized and simplified structures can block 90 % of the incident energy in [1420 Hz, 5140 Hz] and [1571 Hz, 5150 Hz], respectively, and both show excellent omnidirectional sound insulation performance. In addition, the ventilative capabilities are further evaluated by simulating the airflow across the two structures. The proposed barrier may provide new insight into the acoustic applications calling for soundproofing in free surroundings.

Surface Waves on Non-conducting Structured Interfaces

Surface waves are demonstrated theoretically to propagate along periodically corrugated surfaces made of non-conductive lossless materials with positive permittivity. An analytic derivation of the surface wave dispersion relation is presented in the case of the deeply subwavelength period of a corrugated structure using impedance boundary conditions (IBCs) at the interfaces. Thus obtained dispersion relation is verified numerically and limitations as well as the physicality of IBCs when modeling surface waves are discussed. Finally, suitable materials, potential experimental realization, and sensing applications of surface waves in the terahertz spectral range are discussed.

Threshold effect on the femtosecond laser-induced periodic subwavelength structure: An analytical approach

For linearly polarized femtosecond laser-induced periodic nanostructures on a metal surface, the underlying physics of self-formation patterns is yet to be fully explained. We demonstrated here that the stabilized periodic interface electrostatic field (Coulomb field) arising from the surface phase-locking resonant dominates the self-formation of ripple structures. The kinetics of ultrafast laser ablation on solid surface were discussed. The threshold fluence and upper limit fluence for self-formation of ripple structures were calculated. The underlying mechanisms were clarified with both simulation and experimental results. Analytical formula for the interspace between ripples (i.e., wavelength of surface Langmuir wave) was developed. As illustrated by the experimental results, the competition of the two ablation mechanisms, namely surface Coulomb ablation and direct ultrafast laser ablations, was shown to result in distinct structuring patterns for a typical Gaussian pulse laser irradiation on the target surface at the appropriate laser energy fluence.

Sub-wavelength patterned pulse laser lithography for efficient fabrication of large-area metasurfaces

Rigorously designed sub-micrometer structure arrays are widely used in metasurfaces for light modulation. One of the glaring restrictions is the unavailability of easily accessible fabrication methods to efficiently produce large-area and freely designed structure arrays with nanoscale resolution. We develop a patterned pulse laser lithography (PPLL) approach to create structure arrays with sub-wavelength feature resolution and periods from less than 1 μm to over 15 μm on large-area thin films with substrates under ambient conditions. Separated ultrafast laser pulses with patterned wavefront by quasi-binary phase masks rapidly create periodic ablated/modified structures by high-speed scanning. The gradient intensity boundary and circular polarization of the wavefront weaken diffraction and polarization-dependent asymmetricity effects during light propagation for high uniformity. Structural units of metasurfaces are obtained on metal and inorganic photoresist films, such as antennas, catenaries, and nanogratings. We demonstrate a large-area metasurface (10 × 10 mm2) revealing excellent infrared absorption (3–7 μm), which comprises 250,000 concentric rings and takes only 5 minutes to produce.