Research ArticleOPTICS

Bandwidth-unlimited polarization-maintaining metasurfaces

See allHide authors and affiliations

Science Advances  29 Jan 2021:
Vol. 7, no. 5, eabe1112
DOI: 10.1126/sciadv.abe1112


Any arbitrary state of polarization of light beam can be decomposed into a linear superposition of two orthogonal oscillations, each of which has a specific amplitude of the electric field. The dispersive nature of diffractive and refractive optical components generally affects these amplitude responses over a small wavelength range, tumbling the light polarization properties. Although recent works suggest the realization of broadband nanophotonic interfaces that can mitigate frequency dispersion, their usage for arbitrary polarization control remains elusively chromatic. Here, we present a general method to address broadband full-polarization properties of diffracted fields using an original superposition of circular polarization beams transmitted through metasurfaces. The polarization-maintaining metasurfaces are applied for complex broadband wavefront shaping, including beam deflectors and white-light holograms. Eliminating chromatic dispersion and dispersive polarization response of conventional diffractive elements lead to broadband polarization-maintaining devices of interest for applications in polarization imaging, broadband-polarimetry, augmented/virtual reality imaging, full color display, etc.


Polarization, one of the key parameters of electromagnetic wave, plays an important role in many areas, such as three-dimensional stereoscopic displays, quantum computation, chemical sensing, etc. Conventional generation method of arbitrary polarization involves bulky optical components combining polarizers and waveplates. In recent years, metasurfaces with optically thin thickness have been proposed by locally designing subwavelength structures to control all of the parameters of the electromagnetic wave, such as phase (1, 2), frequency (3), amplitude (4), propagation direction (5, 6), and polarization (710), leading to various applications, including lenses (11), cloaking (12), sensing (13, 14), information security (15), etc. The manipulation of polarization state of light using metasurfaces usually relies on designing the structural properties of anisotropic (16) and chiral nanostructures (17). The superposition of two orthogonal polarizations has been proposed to generate versatile output polarization, simply by controlling the offset between the two polarization components (18, 19). However, to generate full arbitrary polarization, the amplitude difference of the two polarization components has to be precisely controlled (20). Amplitude control has generally been performed through the nanostructure geometric sizes, at the expense of introducing a huge dispersion of the amplitude response, compromising the device broadband operation, as schematically shown in Fig. 1 (A and B) .

Fig. 1 Linear superposition of CP beams using phase gradient metasurfaces.

(A) A linear superposition of right CP (RCP) and left CP (LCP) with different sizes of metastructures. (B) Schematic representation of the dispersive polarization response of the traditional polarization-controlling metasurfaces. (C) Broadband polarization-maintaining design consisting in combination of RCP, LCP, and LCP (or RCP, RCP, and LCP) realized with a uniform size metastructure. (D) Schematic representation of the broadband polarization-maintaining metasurfaces. (E) Angular nondispersive design with metasurface doublet. (F) Schematic representation of the angular nondispersive metasurfaces.

Here, we propose a general method to maintain the full-polarization state of diffracted signal, essentially over an unlimited wavelength range, using appropriate superposition of orthogonal circular polarization (CP) (see schematic in Fig. 1, C and D). We also prove that our approach is compatible with basic chromatic dispersion compensation using a simple metasurface doublet (see schematic in Fig. 1, E and F), thus demonstrating a broadband white-light holographic image with broadband polarization-maintaining performances. The control of polarization helicity, i.e., left CP (LCP) or right CP (RCP), and the direction of the scattered light are controlled through an arrangement of spatially oriented birefringent geometric phase metastructures. This peculiar method to address wavefront ranging from basic beam deflection to complex holographic field distribution (2128) relies on replication and rotation of a unique nanostructure. In general, the intrinsic dispersive modal response of the nanostructure creates strong dispersive response, but in this contribution, we report an interesting superposition regime for which the dispersive contribution of the antenna does not comprehend the overall polarization response. The method consists in superposing: (i) diffracted fields with the same handiness of CP to control the output light amplitude, essentially relying on simple interference condition, and (ii) an additional phase-retarded diffracted field oscillating on the opposite handiness, to access full polarization of the output light. Previously, vectorial meta-hologram able to address the problem of polarization-dependent image projection has been proposed, increasing the design degrees of freedom to enable polarization multiplexing (2934). A reflective meta-hologram based on diatomic metasurface has been proposed for the full polarized meta-hologram at oblique incidence (35, 36). Full polarized meta-hologram can also be realized by combining geometric phase and propagation phase, which suffers from narrow bandwidth (3739). Here, we applied our concept of polarization-maintaining metasurfaces to encode complex phase profiles with geometric phase gradient metasurfaces and experimentally demonstrated arbitrary polarization holographic images with broadband characteristics.


Design method of polarization-maintaining metasurfaces

Our general concept of the proposed broadband polarization-addressing metasurfaces is illustrated in Fig. 1 (C and D) and compared with respect to the dispersive design in Fig. 1 (A and B). A 1-μm-thin GaN film on the x-y plane is etched to form rectangular pillars suitably disposed and oriented to acquire the desired geometric phase and shape a given incident CP into a deflected beam with opposite helicity. The excitation input light beam propagating along the −z direction is chosen to be linearly polarized in the x direction and can thus be decomposed into two CP beams with equal amplitude and phase. Note that the incident light can be arbitrary polarization except for two pure CP beams (see more details in notes S3 and S4). As considered in the literature (22), a single Pancharatnam−Berry (PB) phase gradient array is able to direct the two incident CP beams into deflected beams with opposite helicity accordingly to their incident polarization spin. In this first example (Fig. 1A), the dispersive case, the first metastructure line in the bottom of the metasurface is arranged counterclockwise with an orientation angle increment of δd to deflect the incident RCP into an LCP beam propagating at an angle of θtd = arcsin (λδdP), where λ is the operating wavelength and P is the period of the metastructures. The complex transmission of LCP beam is a1(λ)eiφ1(λ)eiL (see eq. S3) with starting rotation angle (SOA) of δL, where the conversion efficiency a1(λ) and propagation phase φ1(λ) are usually dispersive and related to the size of the bottom metastructures, while the geometric phase −2δL is only related to the rotation angle of the metastructures. The additional metastructures disposed in the top line of the metasurface are arranged clockwise with the same angle increment of δd to convert the LCP into an RCP beam propagating at the same angle of θtd. The transmitted RCP has the complex amplitude of a2(λ)eiφ2(λ)eiR according to eq. S3 with SOA of δR, where a2(λ) and φ2(λ) are related to the size of the top metastructures. Therefore, the output state of polarization (SoP) can be described asnLR=a1(λ)eiφ1(λ)ei2δL++a2(λ)eiφ2(λ)ei2δR(1)

The ellipticity angle becomes χ=12arcsina2(λ)2a1(λ)2a2(λ)2+a1(λ)2, and the azimuth angle is ψ=φ2(λ)φ1(λ)2+δR+δL. Since a1,2(λ) and φ1,2(λ) are usually dispersive, the output polarization is related to the wavelength as shown in the schematic of Fig. 1B.

To address this problem, we propose an interesting superposition using three lines of phase gradient metasurface as shown in Fig. 1C. The broadband polarization-maintaining design consists of controlling the diffraction of one LCP and two RCP with SOA of δL1, δR1, and δR2. The key idea is to homogenize the dispersion effect using a unique metastructure, i.e., only one nanopillar with a uniform size, for all lines, so that all transmitted LCP and RCP beams share the same conversion efficiency a(λ) and propagation phase φ(λ). According to eq. S3, the output polarization can be written asnLRR=a(λ)eiφ(λ)ei2δL1++a(λ)eiφ(λ)ei2δR1+a(λ)eiφ(λ)ei2δR2=a(λ)eiφ(λ)(ei2δL1++ei2δR1+ei2δR2)(2)

It can be seen that the dispersive term a(λ)eiφ(λ) factorizes, leading to the generated polarization only related to the rotation angle of δL1, δR1, and δR2, yielding to a nondispersive polarization behavior with unlimited bandwidth as shown in the schematic of Fig. 1D. Note that the deflection efficiency, which depends on the polarization conversion efficiency of the building block, changes as a function of the incident wavelength. Further efforts in designing broadband half-waveplate nanostructure could be realized to improve the device performance over much larger spectral range. Since the diffracted angles of the converted light beam θtd (eq. S7) of both the dispersive and nondispersive polarization design are angular dispersive, a metasurface doublet is designed to eliminate the angular dispersion as shown in Fig. 1 (E and F) (see more details in note S5).

Simulated results and dispersion analysis

We performed the numerical simulation of realistic nanostructures using full-wave finite-difference time domain. The metastructure consists of GaN nanopillars on a sapphire substrate as shown in Fig. 2 (A and B). The period of the metastructure is chosen to avoid spurious diffraction effects in the substrate as P = 300 nm, and the thickness of the GaN nanopillar is h = 1 μm (see more simulated results in fig. S5). To better explain the difference between our proposed broadband polarization-addressing method using a unique nanostructure size and the conventional dual-nanostructure reconstruction method, i.e., dispersive superposition obtained with different nanostructure sizes (20), we carried out numerical simulation for both methods as shown in Fig. 2. We first numerically calculated the polarization conversion efficiency a(λ) (red and black lines) and propagation phase φ(λ) (blue and green lines) with two different sizes of GaN nanopillars as shown in Fig. 2C. As expected, one can see that both a(λ) and φ(λ) are highly dependent on the size of the nanostructures. Thereafter, we simulated the azimuth angle and ellipticity angle of the output polarization as a function of the incident wavelength for both cases (Fig. 2, D and F, for dispersive case, calculated from Eq. 1, and polarization-maintaining case, calculated from Eq. 2, respectively) and the corresponding SoP on Poincaré sphere (Fig. 2, E and G, for dispersive case and polarization-maintaining case, respectively). The dispersive superposition relying on two nanostructures with different sizes (Lu1 = 170 nm, Lv1 = 70 nm, and δL = 0° for the bottom line and Lu2 = 220 nm, Lv2 = 120 nm, and δR = 0° for the top line) exhibits strong wavelength-dependent SoP (Fig. 2E), since a1,2(λ) and φ1,2(λ) have different dispersive responses according to Fig. 2C. However, our design with a uniform size of GaN nanopillars (Lu2 = 220 nm and Lv2 = 120 nm) oriented at different starting angles to address the chosen SoP eliminates dispersion even in the presence of strong nanostructure dispersive response as shown in Fig. 2 (F and G) with δL1 = 0°, δR1 = 22.5°, and δR2 = 22.5°.

Fig. 2 Numerical simulation of the metasurfaces.

(A) Top view and (B) perspective view of one metastructure. P = 300 nm and h = 1000 nm. (C) Simulated results of the polarization conversion efficiency a(λ) (red and black lines) and propagation phase φ(λ) (blue and green lines) with two different sizes of GaN nanopillars. Red and blue curves: Lu1 = 170 nm and Lv1 = 70 nm. Black and green curves: Lu2 = 220 nm and Lv2 = 120 nm. (D) Calculated azimuth angle and ellipticity angle of the output polarization and (E) the corresponding SoP on Poincaré sphere as a function of incident wavelength from 475 to 675 nm with dispersive superposition consisting of two different sizes of GaN nanopillars with zero SOA, indicating that the traditional linear superposition of opposite polarization leads to dispersive polarization response. (F) Calculated azimuth and ellipticity angles of the output polarization and (G) the corresponding SoP on Poincaré sphere as a function of incident wavelength from 475 to 675 nm with the proposed nondispersive superposition consisting of a uniform size of GaN nanopillars with Lu2 = 220 nm, Lv2 = 120 nm, δL1 = 0°, δR1 = 22.5°, and δR2 = 22.5°, indicating that an assembly of a generic building block could realize broadband polarization-maintaining behavior with unlimited bandwidth.

Full-polarization generation

Our superposition design not only addresses nondispersive polarization properties, but it can also generate arbitrary and full polarization with a uniform size of the metastructures. To simplify the discussion, we propose to refer it as LRR configuration (the design using one LCP and two RCP) as shown in Fig. 3A, consisting of phase gradient supercells with controllable SOA of δL1, δR1, and δR2. Since the output polarization is related to the relative phase difference between these three components, we fix δL1 as zero. According to Eq. 2, the dispersive term a(λ)eiφ(λ) factorizes, so the output SoP becomes ∣nLRR〉 = ∣+〉 + (eiR1 + eiR2)∣−〉 = ∣+〉 + aeiδ∣−〉, where a and δ are the amplitude and phase of (eiR1 + eiR2), respectively. Therefore, the ellipticity angle becomes χ=12arcsina21a2+1, and the azimuth angle is ψ=δ2. The ellipticity and azimuth angles as a function of δR2 are shown in Fig. 3 (B and C), respectively. Since 0 ≤ a ≤ 2, the ellipticity angle of the SoP has an upper limit of 12arcsin35 (i.e., π4χ12arcsin35). The corresponding Poincaré sphere with δR1 = 0°, 30°, and 60° are shown in Fig. 3 (D to F), respectively. A full azimuth angle is realized as shown in Fig. 3G by controlling both δR1 and δR2 from 0° to 180°. Therefore, the generated SoPs cover the bottom sphere cap with π4χ12arcsin35. Likewise, the RLL configuration (one RCP and two LCP) shown in Fig. 3 (H to N) can generate SoPs that cover the top sphere cap with 12arcsin35χπ4 as shown in Fig. 3N. In this way, depending on the choice of the desired diffracted state of polarization, one could choose in the LRR and RLL configurations, leading to full-polarized and broadband properties of output light. Note that using twofold superposition with a uniform metastructure, the dispersion term in Eq. 1 can also be factorized, i.e., a1(λ)eiφ1(λ) = a2(λ)eiφ2(λ), leading to a nondispersive property but only for linear output polarization.

Fig. 3 Design of full-polarization pixels using geometric phase metasurface.

(A) One supercell of the metasurface consists of one LCP and two RCP beams. (B) Ellipticity angle and (C) azimuth angle of the output light as a function of rotation angle δR2 for fixed δL1= 0°. Black curve: δR1 = 0°; red curve: δR1 = 30°; blue curve: δR1 = 60°. The corresponding Poincaré sphere when (D) δR1 = 0°, (E) δR1 = 30°, and (F) δR1 = 60°. (G) By controlling both δR1 and δR2 from 0° to 180°, the SoP covers the bottom sphere cap of Poincaré sphere with π4χ12arcsin35. (H) One supercell of metasurface consists of one RCP and two LCP beams. (I) Ellipticity angle and (J) orientation angle of the output light as a function of rotation angle δL2 for fixed δR1 = 0°. Black curve: δR1 = 0°; red curve: δR1 = 30°; blue curve: δR1 = 60°. The corresponding Poincaré sphere when (K) δR1 = 0°, (L) δR1 = 30°, and (M) δR1 = 60°. (N) By controlling both δL1 and δL2 from 0° to 180°, the SoP covers the top sphere cap of Poincaré sphere with 12arcsin35χπ4. (G) and (N) indicate that the superposition of CPs diffracted by array of a generic PB building block could produce any arbitrary SoP, without having to carefully engineers the metastructure geometry for the desired application.

Unlimited bandwidth polarization-maintaining meta-hologram

As a proof of principle, we design four elliptical polarizations using the above mentioned LRR and RLL configurations as shown in Fig. 4. The deflector is pixelated into pixel array containing four metastructures in x direction and three lines in y direction, per pixel. By encoding the holographic phase profiles into pixelated metasurfaces, we are able to generate holographic images with arbitrary polarization and broadband characteristics. The details of the rotation angles and the corresponding ellipticities and azimuth angles are shown in the table S1. We encode the holographic phase profiles to exhibit four playing card suits of “spade,” “heart,” “club,” and “diamond.” Scanning electron microscopy (SEM) images of the fabricated metasurfaces with top and tilt view is shown in Fig. 4 (A and B). The measured holographic images illuminated with coherent laser source at λ = 600 nm are shown in Fig. 4C (see more details of fabrication processes and measurement setup in Materials and Methods and figs. S9 and S11). The measured SoP is plotted as blue dots in the Poincaré sphere as shown in Fig. 4D, which agrees well with the designed SoP denoted by the red dots. The broadband response of the polarization is demonstrated in Fig. 5. The holographic images of spade with the incident wavelength ranging from λ = 475 nm to 675 nm are shown in the left panels of Fig. 5 (A to I). To characterize the polarization on a broadband wavelength range, we have selected a Thorlabs mounted Superachromatic quarter waveplate (AQWP05M-600) with rotation angle of θλ/4 = 67.5° and a broadband Thorlabs linear polarizer (WP25M-VIS) with rotation angle of θLP = 49.07° that were positioned between the meta-hologram and the projector to block the holographic images of designed SoP (details of setup are shown in the fig. S11). With this choice of waveplate/polarizer orientations, holographic images are almost totally blocked in the entire visible range as shown in the right panels of Fig. 5 (A to I), thus confirming the exceptional broadband response of our polarization-maintaining design approach. The viewing angle of the holographic image is related to the operating wavelength as Θ = 2 arcsin (λ/2Ph), where Ph is the pixel size of the hologram. Therefore, the blue holographic image has smaller size than the red holographic image. To further quantify experimentally the SoP, we projected the holographic images into a broadband polarimeter, which indicated that the SoP is correctly maintained around the designed value as shown in Fig. 5J. Since the pixel size of the meta-hologram is larger than the operating wavelength, there are higher-order images generated, which decrease the efficiency of the interested order as measured in Fig. 5K. It is possible to increase the efficiency by decreasing the pixel size to subwavelength scale. We have added a simulation proving that higher efficiency with subwavelength pixel size is achievable, referring to fig. S6. Note that the unlimited bandwidth of polarization-maintaining behavior is conserved for all of the input polarization (see details in note S2). The only limitation of our design approach concerns the addressing of arbitrary polarization state that are realized for all of the input polarization except for two pure input CP beams (see details in note S4).

Fig. 4 Full polarized meta-hologram.

(A) Top view and (B) tilt view of the SEM images of the metasurfaces. The red rectangular highlight area in (A) shows one pixel of the meta-holograms. (C) Measured holographic images (illumination at λ = 600 nm) of the meta-holograms. (D) Measured SoP on Poincaré sphere. Red dots, designed value; blue dots, measured value. Scale bars, 1 μm (A and B).

Fig. 5 Broadband polarization-maintaining properties of the nondispersive metasurface.

(A to I) Measured holographic images with the incident wavelength from 475 to 675 nm (left: without analyzer; right: with broadband quarter waveplate at θλ/4 = 67.5° and linear polarizer (LP) at θLP = 49.07° to block the holographic images). The holographic images in the right panels are almost totally blocked in the entire visible range, confirming the exceptional broadband response of out polarization–maintaining design approach. (J) Measured SoP at different wavelengths agrees well to the designed value in Fig. 2G. (K) Measured efficiency is higher than 10% across the entire visible range.

Angular nondispersive wavefront shaping using cascaded metasurface

According to eq. S7, the deflection angle is, in general, subjected to the incident wavelength, i.e., in other words, the deflection function of the device remains chromatically dispersive. However, this angle-dependent dispersive response can be treated using appropriate combination of metasurfaces with refractive dispersive materials (40), which would resolve this issue (see more details in note S5). The idea to address chromatic angular dispersion is straightforward and consists in using cascaded parallel metasurfaces, designed for α = 0° (see definition of α in fig. S7A), as shown in Fig. 6A. Thus, the transmitted angle can be described from eq. S27 asθt=arcsin(λδdπP+λδdπP)(3)

Fig. 6 Angular nondispersive wavefront shaping using cascaded metasurface doublet.

(A) Schematic of angular nondispersive design using two layers of metasurfaces on both sides of a sapphire substrate. The phase gradient of the two-layer metasurfaces are dφdx=2δdP on the top and dφdx=2δdP on the bottom with δd=δd=30°. A linearly broadband light is normally illuminated from the top metasurface and diffracted to the bottom metasurface. Two spots are generated from the top deflector as shown in the green optical path. (B) Transmitted angle with single-layer metasurface (red curve) and double-layer metasurfaces (blue curve). A nondispersive angle response is realized in unlimited bandwidth using double-layer metasurfaces. Measurement results of the holographic image of spade using (C) single-layer metasurface and (D) double-layer metasurfaces. The holograms are imaged at the entrance slit of a spectrometer for spectrum characterization. A slit with 20 μm in width is placed to select a line of the image light at specific position of P1, P2, … P5. The corresponding spectra for (E) dispersive case and (F) nondispersive case at position P1 to P5 show a red shift of the bandwidth for large angle in (E), while (F) spectra are essentially broadband.

It can be seen that when the increment rotation angle of the two metasurfaces are opposite (δd=δd), i.e., when the angular dispersive properties of each metasurfaces compensate, a constant transmission angle θt = 0° is obtained, leading to a perfect angular nondispersive metasurface with unlimited bandwidth and without any additional approximations as shown in Fig. 6B. The fabrication details and results are shown in fig. S10. The same optical setup of fig. S11 is used here to capture the angular nondispersive holographic image, but the transmitted angle of the image is zero. The experimental results of the dispersive and nondispersive holographic images are shown in Fig. 6 (C and D), respectively. Two spots observed at the exit of the doublet, as shown in Fig. 6D, are induced by the first layer of deflector as shown in the green light path in Fig. 6A (explaining the weaker intensity on right spot with respect to left counterpart). The holographic images are then projected at the entrance slit of a spectrometer to perform spectrum characterization at specific position along the image located at P1, P2, … P5 (which is chosen from left edge to the right edge of the image) as shown in Fig. 6 (C and D). As it is expected for the dispersive case, the spectrum of the hologram red shift according to the angular dispersion law for larger deflection angles, as shown in Fig. 6E. Therefore, P1 has more intensity in the blue range, while P5 has more in the red range. Note that the relative intensities in P1 to P5 are changing because of the size of the actual portion of the image considered in the experiment. Comparatively, the nondispersive and polarization-maintaining case, all spectra always remain broadband at all the positions across the image as shown in Fig. 6F.


In conclusion, we have demonstrated a general method of full-polarization generation and bandwidth-unlimited polarization-maintaining diffractive planar optics. It relies on the superposition of orthogonal polarizations generated by PB phase metasurfaces composed of a unique element to eliminate the dispersive response of the structures. By linear superposition of amplitude-controlled polarization handiness, using RLL and LRR configurations, full-polarization generation is established over an exceptionally unlimited bandwidth, experimentally demonstrated across the entire visible range from 475 to 675 nm. The realization of a metasurface doublet, involving an additional beam deflector with opposite deflection properties, is shown to address both the angular and the polarization chromatic dispersion over an extremely large bandwidth. The capability of generating and maintaining any arbitrary state of polarization on the entire visible range would lead to a plethora of promising applications in full-color displays, augmented/virtual reality display, broadband polarization camera, vector beam generation, visible light communication, etc.


Device fabrication

The fabrication processes are shown in fig. S9. A molecular beam epitaxy (MBE) Riber system is used to grow GaN thin film with 1 μm in thickness on a double-side polished c-plan sapphire substrate. The GaN nanopillars are patterned using conventional electron beam lithography processes. We spin-coat ~180 nm of polymethyl methacrylate (PMMA) resist (495A4) on the GaN thin film and bake it on a hot plate at a temperature of 125°C. The PMMA resist is exposed with designed patterns using electron beam lithography at 20 keV (Raith Elphy Plus, Zeiss Supra 40) and developed using 3:1 isopropyl alcohol:methyl isobutyl ketone solution. Subsequently, a 50-nm-thick Ni is deposited on the sample using E-beam evaporation. After lift-off process by immersing the sample into acetone solution for 2 hours, a Ni hard mask is obtained. The GaN nanorods are created by reactive ion etching (Oxford System) with a plasma composed of Cl2CH4Ar gases, followed by chemical etching with 1:1 H2O2:H2SO4 solution to remove the Ni hard mask, revealing a GaN nanopillar array.

Optical setup

The optical setup for meta-hologram characterization is shown in fig. S11. A laser beam propagates through a broadband linear polarizer (WP25M-VIS) with axis of transmission in horizontal direction to generate linear polarized input beam. It is weakly focused by an achromatic lens with a focal length of 50 mm onto the meta-hologram, which is mounted on a three-dimensional translation stage. The holographic image is projected onto a projector placed 10 cm away from the meta-hologram. A broadband quarter waveplate with fast axis at angle θλ/4 and a broadband linear polarizer with axis of transmission at angle θLP are used to block selected images. We assume an incident beam Ein propagates through the selected quarter waveplate and linear polarizer and describe the output electric field Eout as Eout=ALP(θLP)Aλ/4(θλ/4)Ein. Selected rotation angles θλ/4 and θLP are chosen to block the specific output images (i.e., Eout=0).


Supplementary material for this article is available at

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Funding: We acknowledge funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 639109). P.G. acknowledges support from French defense procurement agency under the ANR ASTRID Maturation program, grant agreement number ANR-18-ASMA-0006-01. Author contributions: Q.S. and P.G. conceived the idea and carried out the experiment. Q.S. performed the numerical simulation and optical characterization of meta-hologram. Q.S., S.C., and V.B. contributed to the nanofabrication. S.V., B.D., and P.D.M. performed the GaN MBE growth. Q.S., S.K., and P.G. wrote the manuscript. P.G. supervised and coordinated the project. All the authors discussed the results and approved the paper. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the corresponding author.

Stay Connected to Science Advances

Navigate This Article