www.nature.

com/scientificreports

OPEN Twisting phase and intensity of

light with plasmonic metasurfaces
Yuchao Zhang, Xiaodong Yang    & Jie Gao

Twisting light in both phase and intensity has recently drawn great interests in various fields related
Received: 11 January 2018
to light-matter interactions such as optical manipulation of particles and quantum entanglement of
Accepted: 12 March 2018 photons. Conventionally, bulky optical components are required to produce such twisted optical beams,
Published: xx xx xxxx which significantly limits their applications in integrated photonics and optical chips. Here, we design
and demonstrate aluminum plasmonic metasurfaces consisting of nanoslit antennas as ultracompact
beam converters to generate the focused twisted beams in both phase and intensity across the visible
wavelength range. The metasurface is encoded with the combined phase profile containing the helico-
conical phase function together with a Fourier transform lens based on the Pancharatnam-Berry
(PB) geometric phase. It is demonstrated that the created twisted beams simultaneously possess
three-dimensional (3D) spiral intensity distribution around the propagation axis and complex phase
structure containing both the central vortex and the peripheral vortex string. Moreover, the twisted
beam exhibits an arithmetic intensity spiral at the focal plane with the maximum photon concentration
located at the leading point of the spiral. Our results show the promising potential for advancing
metasurface-based integrated devices in many applications of light-matter interactions.

Optical vortices having helical phase profiles and phase singularities have drawn considerable attention in many
exciting areas related to light-matter interactions, such as quantum optics1–5, high-resolution imaging6,7, optical
communications8,9, and optical manipulation of particles10–12. The phase profile of optical vortex is expressed by
exp(ilϕ), where ϕ is the azimuthal angle and l represents the topological charge (TC). The intensity distribution of
optical vortex exhibits the doughnut-shaped ring with the ring radius depending on the TC. As the TC is
increased, the vortex ring radius becomes larger and the photon density decreases, which has limitations for
applications requiring both large orbital angular momentum (OAM) and high photon concentration, such as cold
atom rotation13–15 and quantum entanglement of photons16,17. One solution for the above issue is to consider
optical beams with both twisted phase and intensity profiles18 so that TC-independent high photon density can
be realized in the focused areas. One important twisted beam in both phase and intensity is helico-conical (HC)
beam18–22, which can be used to introduce spiral motion on the trapped particles by the transfer of OAM. The
conventional method to generate HC beam requires free-space bulky optical components including spatial light
modulator and Fourier transform lens, which increases the optical system complexity and also limits the photonic
chip integration.
In recent years, plasmonic metasurfaces made of nanoantenna arrays in ultrathin metallic films have provided
a powerful and functional platform for tailoring the phase, intensity and polarization of light at the subwave-
length scale23,24. By introducing the Pancharatnam-Berry geometric phase accompanied with polarization con-
version25–31, metasurfaces have been widely used for building on-chip wavefront shaping devices such as optical
vortex generators32–36, flat optical lenses37–42 compact wave plates43–45, and multiplexed holograms46–51. Although
noble metals like gold and silver have been widely employed in plasmonic metasurfaces in the visible spectrum,
gold has limitations for working below the wavelength of 550 nm due to the interband transition and silver is
susceptible to oxidation and sulphidation under ambient condition. Alternatively, aluminum has been identified
as a promising substitute for gold and silver for plasmonic metasurfaces operating within the visible and near-UV
frequency regions47,52–55 due to its high plasma frequency, chemical and thermal stability thanks to its natively
formed oxidation layer, low cost, and complementary metal oxide semiconductor (CMOS) compatibility.
In this work, we present aluminum plasmonic metasurfaces constructed from nanoslit antenna arrays to gen-
erate HC beams with both twisted phase and intensity profiles in the visible wavelength range from 400 nm to
800 nm. The metasurface with compact area of 50 μm by 50 μm is encoded with the PB phase profile created by

Department of Mechanical and Aerospace Engineering, Missouri University of Science and Technology, Rolla, MO,
65409, USA. Correspondence and requests for materials should be addressed to X.Y. (email: [email protected]) or
J.G. (email: [email protected])

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 1

www.nature.com/scientificreports/

Figure 1. (a) Schematic of the unit cell design of the nanoslit antenna at rotation angle of θ. (b) Simulated
geometric phase shifts for the transmitted converted spin component from the nanoslit antennas at different
rotation angles of 0, 45°, 90°, 135°, and 180°. (c) Simulated electric field |E| distributions of nanoslit antenna
under circular polarizations at 532 nm. (d) A SEM image of the fabricated homogeneous nanoslit antenna array.
(e) Measured and simulated transmission spectra under circular polarization basis.

combining the HC phase function and the phase of Fourier transform lens. It is shown that the produced HC
beams display 3D spiral intensity trajectory around the propagation axis and complex phase structure containing
both the central vortex and the peripheral vortex string. We also demonstrate that the HC beam with varying TC
always exhibits an arithmetic intensity spiral at the focal plane where the maximum photon density is concen-
trated at the leading point of the spiral. Our demonstrated results will provide new opportunities for realizing
metasurface-based photonic devices used in various applications of light-matter interactions such as quantum
information processing, optical trapping, and optical communications.

Results
Design of aluminum metasurface for HC beam generation.  The metasurface with compact size of 50
μm by 50 μm contains 208 by 208 subwavelength nanoslit antennas with different orientation angles. As shown
in Fig. 1(a), the nanoslit antennas are etched in a thin aluminum film with thickness of 35 nm on glass substrate
using focused ion beam (FIB) method. The width and length of each nanoslit antenna is 80 nm and 160 nm
respectively, and the unit cell period is 240 nm. The introduced PB phase from nanoslit antenna is only deter-
mined by the orientation angle θ of the nanoslit. When the incident circularly polarized beam normally transmits
through the anisotropic nanoslit antenna, the transmitted beam will contain both the original spin component
with no phase shift and the converted spin component with the induced PB phase shift of 2θ which is twice as the
nanoslit rotation angle. The overall phase profile encoded in the metasurface is obtained by arranging the nanoslit
antennas with designed rotation angles to form spatially inhomogeneous array. Figure 1(b) plots the simulated
geometric phase shifts for the transmitted converted spin component from the nanoslit antennas at different rota-
tion angles, showing that the phase shift variations of the converted spin component from the nanoslit antennas
are twice as their rotation angles in the broad wavelength range from 400 to 800 nm. The simulated electric field
|E| distributions of the nanoslit antenna under circular polarization basis (left-handed and right-handed circular
polarizations, LCP and RCP) at 532 nm are shown in Fig. 1(c), indicating strong polarization anisotropy of the
nanoslit. The SEM image of the fabricated homogeneous array of nanoslit antenna is shown in Fig. 1(d). Next,
the transmission spectra under circular polarization basis are measured and simulated in Fig. 1(e), where the
original spin component has LCP (blue solid line) and the converted spin component has RCP (red solid line).
It is observed that the converted spin component reaches to the maximum transmission of around 15% near the
wavelength of 550 nm, where the plasmonic resonance occurs. The simulation results (dot lines) match well with
the experimental data, proving the effectiveness of the designed nanoslit antennas.
The metasurface is formed by rotating the nanoslit antenna within each unit cell with the orientation angle
θ(x, y) determined by half of the transmitted PB phase shift ϕPB(x, y) of the converted spin component. The

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 2

www.nature.com/scientificreports/

Figure 2. (a) Phase profile of metasurface to generate HC beam with K = 1 and l = 3. (b) The HC phase
functions with different parameters (upper row) and the corresponding phase contour lines (lower row).
(c) SEM images of the fabricated metasurface sample.

generated PB phase profile ϕPB(x, y) is a superposition of the phase distributions of HC phase and Fourier trans-
form lens, as shown in Fig. 2(a). The PB phase profile is defined by the following functions,
ϕPB(x , y ) = ϕHC(x , y ) + ϕlens(x , y ) (1)

 x 2 + y 2 

ϕHC(x , y ) = l ⋅ arg(x + yi) ⋅ K − 
 r0 
 (2)

π 2 2 2 
ϕlens(x , y ) =  (x + y ) + f − f 
λ  (3)
where x, y represent the position coordinates of the unit cell in the specific column and row in the antenna array
with respect to the center of the array. In the polar coordinates, Eq. (2) is expressed as ϕHC(r, θ) = lθ(K − r/r0),
which defines a HC phase function with the product of a helical phase and a conical phase. The parameters l is the
TC of the HC beam, K takes the value of either 0 or 1, and r0 determines the period of the conical phase. Compare
to common Laguerre–Gaussian optical vortex beams, the HC beam contains nonseparable radial and azimuthal

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 3

www.nature.com/scientificreports/

Figure 3.  Schematic of the experimental setup.

terms, and the generated far-field intensity has spiral structure. Eq. (3) introduces the phase profile of spherical
wave for a Fourier transform lens with the focal length of f , and λ is the operation wavelength.
The PB phase distribution shown in Eq. (1) will directly create a focused HC beam. In our design, the focal
length of the Fourier transform lens is chosen to be f = 200 μm at the wavelength of 532 nm. The HC phase func-
tion contains three parameters including TC parameter l = 1, 3 and 5, the parameter K = 0 or 1, and the conical
phase parameter r0 = 5 μm. Therefore, six metasurface samples are fabricated and characterized for generating
HC beams with the above selected parameters. In Fig. 2(b), the phase distributions (upper row images) and phase
contours (lower row images) of the six different HC phase functions are plotted, where the phase contour lines
are located at ϕHC = 0. It is observed that the HC phase profiles are featured with spiral structures and horizontal
discontinuous cut lines. As the TC increases, the interval between two adjacent contour lines become smaller,
hence the phase gradients along both radial and azimuthal directions increase, which is significantly different
from the common optical vortex beams. Moreover, the HC phase profiles with K = 0 and K = 1 are distinguished
by the phase structure in the central area although both of them have spiral phase structures outside. For K = 0,
the contour lines tend to repel from each other and there is no central vortex structure, while for K = 1, the con-
tour lines join together at the center and there is a vortex structure with TC of l presenting at the center. Figure 2a
shows the designed phase profile with the parameters of K = 1 and l = 3, and the SEM images of the fabricated
metasurface sample are shown in Fig. 2(c).

Measurement of 3D twisted intensity profile and complex phase structure.  The 3D twisted struc-
tures of the generated HC beams by the metasurfaces are then demonstrated. Figure 3 shows the schematic of
experimental setup. The collimated optical beam from a laser diode at the wavelength of 405 nm, 532 nm or
633 nm is first transmitted through a linear polarizer and a quarter-wave plate to create a LCP beam, which is
focused onto the metasurface by a 10X objective lens. The generated focused HC beam is obtained by passing
through the RCP filter including another set of a quarter-wave plate and a linear polarizer. The 3D twisted beam
structure is then observed by a microscope imaging system including a 20X objective lens, a 0.5X tube lens and a
CCD camera (CCD 1), which is placed on a translation stage behind the metasurface. The beam intensity profiles
at different propagation distances are recorded and assembled for constructing the 3D twisted beam tomography.
The interferometry pattern of HC beam is obtained by using the Mach–Zehnder interferometer, where the first
beam splitter (BS1) splits the incident beam into the signal beam to generate the HC beam and the reference
spherical beam. The second beam splitter (BS2) reflects the generated HC beam to the third beam splitter (BS3)
for obtaining the interferometry pattern recorded by another CCD camera (CCD 2).
The measured 3D twisted tomography for the HC beams with four different parameter sets of (K = 0, l = 1),
(K = 1, l = 1), (K = 0, l = 3), and (K = 1, l = 3) at 532 nm is shown in Fig. 4(a), covering the beam propagation
distance from z = 60 μm to z = 200 μm. It is shown that the location of the maximum beam intensity contin-
uously twists around the central propagation axis, and the twisted intensity trajectory is nearly a conical helix
which is controlled by the parameters of HC beam. As the TC of HC beam increases from 1 to 3, the size of the
twisted intensity structure expands almost three times larger. The parameter K only makes the 3D twisted beam
shape slightly changed but almost does not affect the size of the twisted structure. Figure 4(b) displays the meas-
ured cross-sectional intensity profiles for the HC beam with (K = 1, l = 3) at different propagation distances z
at 532 nm, while Fig. 4(c) shows the corresponding calculated results by using the beam diffraction integration
method. The dashed white line represents the projection of twisted trajectory onto transverse plane, and it shows
a spiral shape. As the beam propagates in the free space, its intensity profile also changes with respect to the prop-
agation distance. At the starting position the intensity profile has many side lobes, then the side lobes are reduced
and finally only the main spiral lobe is observed close to the focal plane of z = 200 μm.
Next, the generated HC beam is interfered with the reference spherical beam through a Mach–Zehnder inter-
ferometer at 532 nm. Figure 5 presents the recorded interferometry patterns at z = 50 μm, showing the results
of K = 0 and K = 1 with TC changes as l = 1, 3 and 5. For K = 0, there is no vortex structure presenting in the
beam center, and the peripheral vortex string containing −1 charged single vortices with the total number of l is
observed along the horizontal phase discontinuous cut line, showing the total OAM of −lħ. On the other hand,
for K = 1, besides the −1 charged peripheral vortex chain, there is one center vortex with charge of +l showing
up, so that the total number of vortices is l + 1 and the overall OAM vanishes. For K = 0 the peripheral vortices
with the same charge repel to each other, while for K = 1 the peripheral vortices having the opposite charge as the
central vortex will be attracted to the beam center. It is noted that the observed single-charged vortex chain with
the same TC located at the phase discontinuous cut line in HC beam is different from that of the fractional vortex

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 4

www.nature.com/scientificreports/

Figure 4. (a) From left to right, the 3D twisted tomography for the focused HC beam with parameter set of
(K = 0, l = 1), (K = 1, l = 1), (K = 0, l = 3) and (K = 1, l = 3) at 532 nm. (b) Measured and (c) simulated cross-
sectional intensity profiles for the HC beam with (K = 1, l = 3) at different propagation distances z.

beam, where a vortex chain with alternating TC of +1 and −1 is presented in the phase cut line. Another differ-
ence is that the number of the vortices in the vortex chain is infinite for the fractional vortex beam56, but it is equal
to the TC for the HC beam. And the reason for these differences lies in the distinct phase profiles between the HC
beam with radially varying phase and the fractional vortex beam with radius-independent phase, although both
beams have the phase edge dislocation along radial direction. In addition, the obtained complex phase structures
of HC beams are very different from the common optical vortex beam where only single +l charged vortex is
presented.
Finally, the broadband response of the metasurface is demonstrated at three different visible wavelengths
of 405 nm, 532 nm and 633 nm. Figure 6(a) shows the measured cross-sectional intensity distributions of HC
beams at the focal plane of z = 200 μm at three different wavelengths for all six different metasurface samples
with parameter sets of (K = 0, l = 1), (K = 1, l = 1), (K = 0, l = 3), (K = 1, l = 3), (K = 0, l = 5) and (K = 1, l = 5). It is
observed that the HC beam with varying TC at each wavelength always shows an arithmetic intensity spiral at the
focal plane, indicating the broadband operation capabilities of the metasurface. The size of the arithmetic spiral

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 5

www.nature.com/scientificreports/

Figure 5.  The recorded interferometry patterns of HC beams at z = 50 μm for (a) K = 0 and (b) K = 1. From left
column to right column, l changes from 1, 3 to 5.

gets expanded as either the TC or the wavelength increases. It is found that with the same TC and wavelength, the
parameter K does not affect the size of the arithmetic spiral much, and it only changes the intensity profile near
the beam center. It is also shown that for the HC beam with any TC, the maximum intensity remains concentrated
at the leading point of the arithmetic spiral, which indicates that high photon density keeps localized at a tiny
interaction region and will not spread out as the TC increases. Such unique feature of the HC beam is useful for
enhanced light-matter interactions where both high photon density and large OAM are necessary. Figure 6(b)
plots the simulated intensity and phase profiles at the focal plane of z = 200 μm at 532 nm. The phase profiles for
both K = 0 and K = 1 cases indicate that the peripheral vortex strings containing single-charged vortices are not
observed in the beam area due to the diffraction at the far field. However, for K = 1 case, the original center vortex
with charge of +l is now decomposed into single-charged vortices with the total number of l distributed along the
arithmetic intensity spiral, while the vortex singularity centers are located inside of the spiral curve.

Discussion
In summary, we have demonstrated the direct generation of 3D focused HC beams across the visible wavelength
range by using aluminum plasmonic metasurfaces consisting of nanoslit antennas. The designed metasurface is
encoded with the combined phase profile by superimposing the HC phase function and a Fourier transform lens.
The 3D spiral intensity distributions and complex phase structures of HC beams have been analyzed in detail.
The demonstrated HC beam generation with aluminum metasurfaces will find many potential applications in
broadband OAM-based light-matter interactions such as particle manipulation and transport, quantum optics
and optical communications.

Methods
Simulations.  The simulations shown in Fig. 1 are conducted by the CST Microwave Studio, where periodic
boundary conditions are employed along both x and y directions in the unit cell. The permittivity of aluminum
is taken from experimental data and the refractive index of glass substrate is 1.45. The beam propagation in free
space shown in Fig. 4c is calculated by using the Fresnel-Kirchhoff diffraction integral:
 → →  ikr
U (x , y , z ) =
1
iλ
∬S U(x0, y0) cos( n , r) −2 cos( n , r′)  e r dS
  (4)
where U (x 0, y0 ) is the complex amplitude distribution located at the z = 0 plane with surface area S and normal
direction →n , r′ is the vector between the source point and a point in the z = 0 plane, r is the vector between the
point in z = 0 plane and a point in the plane at the propagation distance z, and k = 2π/λ is the wavevector.

Sample fabrication.  The metasurface is made in a 35 nm-thick aluminum film deposited on a glass sub-
strate using electron-beam evaporation. The nanoslit arrays are milled in the aluminum film using focused ion
beam (FIB) system (FEI Helios Nanolab 600, 30 kV, 93 pA). The metasurface contains 208 × 208 unit cells, and
each unit cell contains a milled nanoslit with size of 160 nm × 8 nm at a specified orientation angle.

Optical characterization.  The transmission spectra through the metasurface under circular polarization
basis in Fig. 1 are measured with a collimated broadband Tungsten-Halogen source, where a combination of a
linear polarizer and an achromatic quarter-wave plate are used to convert the incident light to circularly polarized
wave. The light beam is focused normally onto the sample using a 50× objective lens and the transmitted light

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 6

www.nature.com/scientificreports/

Figure 6. (a) Measured intensity profiles at the focal plane of z = 200 μm. Form the first row to the last row, the
corresponding HC metasurface parameters change with the order: (K = 0, l = 1), (K = 1, l = 1), (K = 0, l = 3),
(K = 1, l = 3), (K = 0, l = 5) and (K = 1, l = 5). And from left column to right column, the laser wavelength
changes with the order 405 nm, 532 nm and 633 nm. (b) Simulated intensity and phase profiles at the focal plane
of z = 200 μm at 532 nm. The yellow spiral curves in the phase profiles represent the intensity trajectories.

is collected by another 10× objective lens to a spectrometer (Horiba, iHR 550). Another set of a quarter-wave
plate and a linear polarizer is used to distinguish the original and converted circular polarization components. A
transparent glass substrate is utilized to normalize the transmission spectra. Since the metasurface operates in a
broad wavelength range from 400 nm to 800 nm, three lasers operating at different wavelengths of 405 nm, 532 nm
and 633 nm are employed in the HC beam experiments, as shown in Fig. 3. The 3D twisted beam structures are
captured by a microscope imaging system with a 20X objective lens, a 0.5X tube lens and CCD camera placed on a
translation stage. The interferometry pattern of HC beam is obtained by using the Mach–Zehnder interferometer.

References
1. Allen, L., Beijersbergen, M. W., Spreeuw, R. & Woerdman, J. Orbital angular momentum of light and the transformation of Laguerre-
Gaussian laser modes. Phys. Rev. A 45, 8185–8189 (1992).
2. Fickler, R. et al. Quantum Entanglement of High Angular Momenta. Science 338, 640–643 (2012).
3. Hiesmayr, B.  C., De Dood, M.  J.  A. & Löffler, W. Observation of Four-Photon Orbital Angular Momentum Entanglement. Phys.
Rev. Lett. 116, 073601 (2016).
4. Marzlin, K. P., Zhang, W. & Wright, E. Vortex coupler for atomic Bose-Einstein condensates. Phys. Rev. Lett. 79, 4728–4731 (1997).
5. Holland, M. J. & Williams, J. E. Preparing topological states of a Bose-Einstein condensate. Nature 401, 568–572 (1999).
6. Fürhapter, S., Jesacher, A., Bernet, S. & Ritsch-MarteSpiral, M. Phase contrast imaging in microscopy. Opt. Express 13, 689–694
(2005).

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 7

www.nature.com/scientificreports/

7. Wang, J., Zhang, W., Qi, Q., Zheng, S. & Chen, L. Gradual edge enhancement in spiral phase contrast imaging with fractional vortex
filters. Sci. Rep. 5, 15826 (2015).
8. Wang, J. et al. Terabit free-space data transmission employing orbital angular momentum multiplexing. Nat. Photon. 6, 488–496
(2012).
9. Richardson, D., Fini, J. & Nelson, L. Space-division multiplexing in optical fibres. Nat. Photon. 7, 354–362 (2013).
10. Grier, D. A revolution in optical manipulation. Nature 424, 810–812 (2003).
11. Padgett, M. & Bowman, R. Tweezers with a twist. Nat. Photon. 5, 343–348 (2011).
12. Albaladejo, S., Marqués, M. I., Laroche, M. & Sáenz, J. J. Scattering Forces from the Curl of the Spin Angular Momentum of a Light
Field. Phys. Rev. Lett. 102, 113602 (2009).
13. Clifford, M. A., Arlt, J., Courtial, J. & Dholakia, K. High-order Laguerre–Gaussian laser modes for studies of cold atoms. Opt.
Commun. 156, 300–306 (1998).
14. Kohel, J. M. et al. Generation of an intense cold-atom beam from a pyramidal magneto-optical trap: experiment and simulation. J.
Opt. Soc. Am. B 20, 1161–1168 (2003).
15. Andersen, M. F. et al. Quantized rotation of atoms from photons with orbital angular momentum. Phys. Rev. Lett. 97, 170406 (2006).
16. Dorfman, K. E., Schlawin, F. & Mukamel, S. Nonlinear optical signals and spectroscopy with quantum light. Rev. Mod. Phys. 88,
045008 (2016).
17. Lo Gullo, N., McEndoo, S., Busch, T. & Paternostro, M. Vortex entanglement in Bose-Einstein condensates coupled to Laguerre-
Gauss beams. Phys. Rev. A 81, 053625 (2010).
18. Daria, V. R., Palima, D. Z. & Glückstad, J. Optical twists in phase and amplitude. Opt. Express 19, 476–481 (2011).
19. Alonzo, C. A., Rodrigo, P. J. & Glückstad, J. Helico-conical optical beams: a product of helical and conical phase fronts. Opt. Express
13, 1749–1760 (2005).
20. Hermosa, N. P. & Manaois, C. O. Phase structure of helico-conical optical beams. Opt. Commun. 271, 178–183 (2007).
21. Hermosa, N., Rosales-Guzmán, C. & Torres, J. P. Helico-conical optical beams self-heal. Opt. Lett. 38, 383–385 (2013).
22. Bareza, N. & Hermosa, N. Propagation dynamics of vortices in helico-conical optical beams. Opt. Commun. 356, 236–242 (2015).
23. Kildishev, A. V., Boltasseva, A. & Shalaev, V. M. Planar Photonics with Metasurfaces. Science 339, 1232009 (2013).
24. Yu, N. & Capasso, F. Flat optics with designer metasurfaces. Nature Materials 13, 139–150 (2014).
25. Berry, M. V. Quantal phase factors accompanying adiabatic changes. R. Soc. London, Ser. A 392, 45–57 (1984).
26. Berry, M. V. The Adiabatic Phase and Pancharatnam’s Phase for Polarized Light. Proc. J. Mod. Opt. 34, 1401–1407 (1987).
27. Niv, A., Biener, G., Kleiner, V. & Hasman, E. Manipulation of the Pancharatnam phase in vectorial vortices. Opt. Express 14,
4208–4220 (2006).
28. Hasman, E., Kleiner, V., Biener, G. & Niv, A. Polarization dependent focusing lens by use of quantized Pancharatnam-Berry phase
diffractive optics. Appl. Phys. Lett. 82, 328–330 (2003).
29. Bomzon, Ze, Kleiner, V. & Hasman, E. Pancharatnam-Berry phase in space-variant polarization-state manipulations with
subwavelength gratings. Opt. Lett. 26, 1424–1426 (2001).
30. Gori, F. Measuring Stokes parameters by means of a polarization grating. Opt. Lett. 24, 584–586 (1999).
31. Zhang, D., Feng, X., Cui, K., Liu, F. & Huang, Y. Identifying Orbital Angular Momentum of Vectorial Vortices with Pancharatnam
Phase and Stokes Parameters. Sci. Rep. 5, 11982 (2015).
32. Karimi, E. et al. Generating Optical Orbital Angular Momentum at Visible Wavelengths Using a Plasmonic Metasurface. Light: Sci.
Appl. 3, e167 (2014).
33. Zeng, J., Li, L., Yang, X. & Gao, J. Generating and Separating Twisted Light by gradient−rotation Split-Ring Antenna Metasurfaces.
Nano Lett. 16, 3101–3108 (2016).
34. Zeng, J., Gao, J., Luk, T. S., Litchinitser, N. M. & Yang, X. Structuring Light by Concentric-Ring Patterned Magnetic Metamaterial
Cavities. Nano Lett. 15, 5363–5368 (2015).
35. Zeng, J., Luk, T. S., Gao, J. & Yang, X. Spiraling light with magnetic metamaterial quarter-wave turbines. Sci. Rep. 7, 11824 (2017).
36. Liu, Y. et al. Generation of perfect vortex and vector beams based on Pancharatnam-Berry phase elements. Sci. Rep. 7, 44096 (2017).
37. Aieta, F. et al. Aberration-Free Ultrathin Flat Lenses and Axicons at Telecom Wavelengths Based on Plasmonic Metasurfaces. Nano
Lett. 12, 4932–4936 (2012).
38. Ni, X., Ishii, S., Kildishev, A. V. & Shalaev, V. M. Ultra-Thin, Planar. Babinet-Inverted Plasmonic Metalenses. Light: Sci. Appl. 2, e72
(2013).
39. Chen, X. et al. Dual-Polarity Plasmonic Metalens for Visible Light. Nat. Commun. 3, 1198 (2012).
40. Khorasaninejad, M. et al. Multispectral Chiral Imaging with a Metalens. Nano Lett. 16, 4595–4600 (2016).
41. Zhang, K., Ding, X., Zhang, L. & Wu, Q. Anomalous three-dimensional refraction in the microwave region by ultra-thin high
efficiency metalens with phase discontinuities in orthogonal directions. New J. Phys. 16, 103020 (2014).
42. Zhang, K. et al. Phase-engineered metalenses to generate converging and non-diffractive vortex beam carrying orbital angular
momentum in microwave region. Opt. Express 26, 1351–1360 (2018).
43. Yu, N. et al. A Broadband, Background-Free Quarter-Wave Plate Based on Plasmonic Metasurfaces. Nano Lett. 12, 6328–6333
(2012).
44. Zhao, Y. & Alu, A. Tailoring the Dispersion of Plasmonic Nanorods to Realize Broadband Optical Meta-Waveplates. Nano Lett. 13,
1086–1091 (2013).
45. Cao, W., Yang, X. & Gao, J. Broadband polarization conversion with anisotropic plasmonic metasurfaces. Sci. Rep. 7, 8841 (2017).
46. Wan, W., Gao, J. & Yang, X. Metasurface holograms for holographic imaging. Adv. Opt. Mater. 5, 1700541 (2017).
47. Wan, W., Gao, J. & Yang, X. Full-Color Plasmonic Metasurface Holograms. ACS Nano 10, 10671–10680 (2016).
48. Ni, X., Kildishev, A. V. & Shalaev, V. M. Metasurface Holograms for Visible Light. Nat. Commun. 4, 2807 (2013).
49. Huang, L. et al. Three-Dimensional Optical Holography Using a Plasmonic Metasurface. Nat. Commun. 4, 2808 (2013).
50. Zheng, G. et al. Metasurface Holograms Reaching 80% Efficiency. Nat. Nanotechnol. 10, 308–312 (2015).
51. Wang, Q. et al. Broadband metasurface holograms: toward complete phase and amplitude engineering. Sci. Rep. 6, 32867 (2016).
52. Cheng, F. et al. Aluminum plasmonic metamaterials for structural color printing. Opt. Express 23, 14552–14560 (2015).
53. Cheng, F. et al. Enhanced structural color generation in aluminum metamaterials coated with a thin polymer layer. Opt. Express 23,
25329–25339 (2015).
54. Li, Z. et al. All-metal structural color printing based on aluminum plasmonic metasurfaces. Opt. Express 24, 20472–20480 (2016).
55. Wang, W., Rosenmann, D., Czaplewski, D. A., Yang, X. & Gao, J. Realizing structural color generation with aluminum plasmonic
V-groove metasurfaces. Opt. Express 25, 20454–20465 (2017).
56. Berry, M. V. Optical vortices evolving from helicoidal integer and fractional phase steps. J. Opt. A: Pure Appl. Opt. 6, 259–268 (2004).

Acknowledgements
The authors acknowledge support from the Office of Naval Research under Grant No. N00014-16-1-2408, and the
National Science Foundation under Grant No. ECCS-1653032 and DMR-1552871. The authors thank the facility
support from the Materials Research Center at Missouri S&T.

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 8

www.nature.com/scientificreports/

Author Contributions
X.Y., J.G. and Y.Z. conceived the idea of the research. Y.Z. performed all the numerical simulations and
experiments. X.Y. and J.G. directed the research. All authors discussed the results and contributed to the
manuscript.

Additional Information
Competing Interests: The authors declare no competing interests.
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 Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article’s Creative Commons license and your intended use is not per-
mitted 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018

SciEnTific REPOrTS | (2018) 8:4884 | DOI:10.1038/s41598-018-23382-7 9