Citation:

High-efficiency, large-area, topology-optimized metasurfaces


  • Light: Science & Applications  8, Article number: (2019)
More Information
  • Corresponding author:
    Jonathan A. Fan (jonfan@stanford.edu)
  • Received: 09 December 2018
    Revised: 30 April 2019
    Accepted: 05 May 2019
    Published online: 29 May 2019

doi: https://doi.org/10.1038/s41377-019-0159-5

  • Metasurfaces are ultrathin optical elements that are highly promising for constructing lightweight and compact optical systems. For their practical implementation, it is imperative to maximize the metasurface efficiency. Topology optimization provides a pathway for pushing the limits of metasurface efficiency; however, topology optimization methods have been limited to the design of microscale devices due to the extensive computational resources that are required. We introduce a new strategy for optimizing large-area metasurfaces in a computationally efficient manner. By stitching together individually optimized sections of the metasurface, we can reduce the computational complexity of the optimization from high-polynomial to linear. As a proof of concept, we design and experimentally demonstrate large-area, high-numerical-aperture silicon metasurface lenses with focusing efficiencies exceeding 90%. These concepts can be generalized to the design of multifunctional, broadband diffractive optical devices and will enable the implementation of large-area, high-performance metasurfaces in practical optical systems.
  • 加载中
  • [1] Yu, N. F. & Capasso, F. Flat optics with designer metasurfaces. Nat. Mater. 13, 139–150 (2014). doi: 10.1038/nmat3839
    [2] Pahlevaninezhad, H. et al. Nano-optic endoscope for high-resolution optical coherence tomography in vivo. Nat. Photonics 12, 540–547 (2018). doi: 10.1038/s41566-018-0224-2
    [3] Chen, K. et al. Moiré nanosphere lithography. ACS Nano 9, 6031–6040 (2015). doi: 10.1021/acsnano.5b00978
    [4] Tittl, A. et al. Imaging-based molecular barcoding with pixelated dielectric metasurfaces. Science 360, 1105–1109 (2018). doi: 10.1126/science.aas9768
    [5] Silva, A. et al. Performing mathematical operations with metamaterials. Science 343, 160–163 (2014). doi: 10.1126/science.1242818
    [6] Mueller, J. P. B. et al. Metasurface polarization optics: independent phase control of arbitrary orthogonal states of polarization. Phys. Rev. Lett. 118, 113901 (2017). doi: 10.1103/PhysRevLett.118.113901
    [7] Iyer, P. P., Pendharkar, M. & Schuller, J. A. Electrically reconfigurable metasurfaces using heterojunction resonators. Adv. Opt. Mater. 4, 1582–1588 (2016). doi: 10.1002/adom.201600297
    [8] Ding, F. et al. Broadband high-efficiency half-wave plate: a supercell-based plasmonic metasurface approach. ACS Nano 9, 4111–4119 (2015). doi: 10.1021/acsnano.5b00218
    [9] Arbabi, A. et al. Dielectric metasurfaces for complete control of phase and polarization with subwavelength spatial resolution and high transmission. Nat. Nanotechnol. 10, 937–943 (2015). doi: 10.1038/nnano.2015.186
    [10] Yang, J. & Fan, J. A. Topology-optimized metasurfaces: impact of initial geometric layout. Opt. Lett. 42, 3161–3164 (2017). doi: 10.1364/OL.42.003161
    [11] Lalau-Keraly, C. M. et al. Adjoint shape optimization applied to electromagnetic design. Opt. Express 21, 21693–21701 (2013). doi: 10.1364/OE.21.021693
    [12] Molesky, S. et al. Inverse design in nanophotonics. Nat. Photonics 12, 659–670 (2018). doi: 10.1038/s41566-018-0246-9
    [13] Frandsen, L. H. et al. Broadband photonic crystal waveguide 60° bend obtained utilizing topology optimization. Opt. Express 12, 5916–5921 (2004). doi: 10.1364/OPEX.12.005916
    [14] Borel, P. I. et al. Topology optimization and fabrication of photonic crystal structures. Opt. Express 12, 1996–2001 (2004). doi: 10.1364/OPEX.12.001996
    [15] Piggot, A. Y. et al. Inverse design and demonstration of a compact and broadband on-chip wavelength demultiplexer. Nat. Photonics 9, 374–377 (2015). doi: 10.1038/nphoton.2015.69
    [16] Xiao, T. P. et al. Diffractive spectral-splitting optical element designed by adjoint-based electromagnetic optimization and fabricated by femtosecond 3D direct laser writing. ACS Photonics 3, 886–894 (2016). doi: 10.1021/acsphotonics.6b00066
    [17] Sell, D. et al. Ultra-high-efficiency anomalous refraction with dielectric metasurfaces. ACS Photonics 5, 2402–2407 (2018). doi: 10.1021/acsphotonics.8b00183
    [18] Sell, D. et al. Periodic dielectric metasurfaces with high-efficiency, multiwavelength functionalities. Adv. Opt. Mater. 5, 1700645 (2017). doi: 10.1002/adom.201700645
    [19] Lin, Z. et al. Topology-optimized multilayered metaoptics. Phys. Rev. Appl. 9, 044030 (2018). doi: 10.1103/PhysRevApplied.9.044030
    [20] Lalanne, P. & Morris, G. M. Highly improved convergence of the coupled-wave method for TM polarization. J. Opt. Soc. Am. A 13, 779–784 (1996).
    [21] Hugonin, J. P., & Lalanne, P. Reticolo Software for Grating Analysis. (Institut d'Optique, 2005).
    [22] Knop, K. Rigorous diffraction theory for transmission phase gratings with deep rectangular grooves. J. Opt. Soc. Am. 68, 1206–1210 (1978). doi: 10.1364/JOSA.68.001206
    [23] Li, L. F. Use of Fourier series in the analysis of discontinuous periodic structures. J. Opt. Soc. Am. A 13, 1870–1876 (1996). doi: 10.1364/JOSAA.13.001870
    [24] Coppersmith, D. & Winograd, S. Matrix multiplication via arithmetic progressions. J. Symb. Comput. 9, 251–280 (1990). doi: 10.1016/S0747-7171(08)80013-2
    [25] Patterson, D. A. & Hennessy, J. L. Computer Organization and Design. 5th edn. (Morgan Kaufmann, Waltham: MK, 2014).
    [26] Strehl, K. Aplanatische und fehlerhafte abbildung im fernrohr. Z. für Instrum. 15, 362–370 (1895).
    [27] Arbabi, A. et al. Subwavelength-thick lenses with high numerical apertures and large efficiency based on high-contrast transmitarrays. Nat. Commun. 6, 7069 (2015). doi: 10.1038/ncomms8069
    [28] Sell, D. et al. Large-angle, multifunctional metagratings based on freeform multimode geometries. Nano Lett. 17, 3752–3757 (2017). doi: 10.1021/acs.nanolett.7b01082
    [29] Yang, J. J., Hugonin, J. P. & Lalanne, P. Near-to-far field transformations for radiative and guided waves. ACS Photonics 3, 395–402 (2016). doi: 10.1021/acsphotonics.5b00559
    [30] Chew, W. C. & Weedon, W. H. A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates. Microw. Opt. Technol. Lett. 7, 599–604 (1994). doi: 10.1002/mop.4650071304
    [31] Chew, W. C., Jin, J. M. & Michielssen, E. Complex coordinate stretching as a generalized absorbing boundary condition. Microw. Opt. Technol. Lett. 15, 363–369 (1997). doi: 10.1002/(SICI)1098-2760(19970820)15:6<363::AID-MOP8>3.0.CO;2-C
    [32] Teixeira, F. L. & Chew, W. C. Differential forms, metrics, and the reflectionless absorption of electromagnetic waves. J. Electromagn. Waves Appl. 13, 665–686 (1999). doi: 10.1163/156939399X01104
    [33] Edee, K., Granet, G. & Plumey, J. P. Complex coordinate implementation in the curvilinear coordinate method: application to plane-wave diffraction by nonperiodic rough surfaces. J. Opt. Soc. Am. A 24, 1097–1102 (2007). doi: 10.1364/JOSAA.24.001097
    [34] Stratton, J. A. & Chu, L. J. Diffraction theory of electromagnetic waves. Phys. Rev. 56, 99–107 (1939). doi: 10.1103/PhysRev.56.99
    [35] Yang, J. J. & Fan, J. A. Analysis of material selection on dielectric metasurface performance. Opt. Express 25, 23899–23909 (2017). doi: 10.1364/OE.25.023899
    [36] Sell, D. et al. Visible light metasurfaces based on single-crystal silicon. ACS Photonics 3, 1919–1925 (2016). doi: 10.1021/acsphotonics.6b00436
    [37] Palmer, C. & Loewen, E. Diffraction Grating Handbook. 6th edn. (New York: Newport Corporation, Rochester, 2005).
    [38] Yang, J. J., Sell, D. & Fan, J. A. Freeform metagratings based on complex light scattering dynamics for extreme, high efficiency beam steering. Ann. der Phys. 530, 1700302 (2018). doi: 10.1002/andp.201700302
    [39] Chung, J. W. et al. Seamless on-wafer integration of Si(100) MOSFETs and GaN HEMTs. IEEE Electron Device Lett. 30, 1015–1017 (2009). doi: 10.1109/LED.2009.2027914
    [40] Shrestha, S. et al. Broadband achromatic dielectric metalenses. Light.: Sci. Appl. 7, 85 (2018). doi: 10.1038/s41377-018-0078-x
    [41] Chen, W. T. et al. A broadband achromatic metalens for focusing and imaging in the visible. Nat. Nanotechnol. 13, 220–226 (2018). doi: 10.1038/s41565-017-0034-6
    [42] Paniagua-Dominguez, R. et al. A metalens with a near-unity numerical aperture. Nano Lett. 18, 2124–2132 (2018). doi: 10.1021/acs.nanolett.8b00368
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(5)

Article Metrics

Article views(377) PDF downloads(74) Citation(0) Citation counts are provided from Web of Science. The counts may vary by service, and are reliant on the availability of their data.

High-efficiency, large-area, topology-optimized metasurfaces

  • 1. Department of Electrical Engineering, Stanford University, Stanford, CA 94305, USA
  • 2. Department of Applied Physics, Stanford University, Stanford, CA 94305, USA
  • 3. Université Clermont Auvergne, Institut Pascal, BP 10448, F-63000 Clermont-Ferrand, France
  • 4. CNRS, UMR 6602, Institut Pascal, F-63177 Aubière, France
  • Corresponding author:

    Jonathan A. Fan, jonfan@stanford.edu

doi: https://doi.org/10.1038/s41377-019-0159-5

Abstract: Metasurfaces are ultrathin optical elements that are highly promising for constructing lightweight and compact optical systems. For their practical implementation, it is imperative to maximize the metasurface efficiency. Topology optimization provides a pathway for pushing the limits of metasurface efficiency; however, topology optimization methods have been limited to the design of microscale devices due to the extensive computational resources that are required. We introduce a new strategy for optimizing large-area metasurfaces in a computationally efficient manner. By stitching together individually optimized sections of the metasurface, we can reduce the computational complexity of the optimization from high-polynomial to linear. As a proof of concept, we design and experimentally demonstrate large-area, high-numerical-aperture silicon metasurface lenses with focusing efficiencies exceeding 90%. These concepts can be generalized to the design of multifunctional, broadband diffractive optical devices and will enable the implementation of large-area, high-performance metasurfaces in practical optical systems.

Introduction
  • Metasurfaces are optical devices that utilize subwavelength-scale structuring to shape and manipulate electromagnetic waves1. They are powerful complements to bulk refractive and scalar diffractive optics and have a wide range of potential applications in imaging2, lithography3, sensing4, and computing platforms5. To date, a broad range of metasurface design concepts have emerged. The most widely used methods, which we will term "conventional methods, " sample the desired phase profile using discrete phase-shifter elements to form a nanoscale phased array (Fig. 1a). These methods utilize a library of simple, physically intuitive building blocks, including anisotropic waveguides6, Mie resonators7, plasmonic resonators8, and dielectric transmit arrays9, and can quickly produce macroscale device designs. However, these approaches lack the necessary degrees of freedom for realizing high-efficiency in devices that are designed for large-angle deflections, multiple functions, and broadband responses10, thereby preventing metasurfaces from being practically applied in many contexts.

    Fig. 1  Strategies for metasurface design.

    a Conventional approaches sample the desired phase profile at discrete points and specify phase shifting elements to form a nanoscale phased array. b Our approach is to divide the desired phase profile into wavelength-scale, linear sections and use topology optimization to design each section individually. c Computation time versus device size for topology-optimized metasurfaces that are designed using two approaches: direct optimization of the entire metasurface (orange) and optimization of the metasurface after division into 3λ-wide sections (green)

    Topology optimization is an alternative design method that can overcome the efficiency limitations of conventional approaches11. Inverse design methods, including objective-first and adjoint-based topology optimization12, have produced photonic crystals13, 14, optical demultiplexers15, and spectral splitters16 with outstanding performance resulting from the devices' unusual geometries and nonintuitive optical dynamics. More recently, topology optimization has led to high-performance metasurfaces that have a broad range of capabilities, such as high-efficiency light deflection with nearly arbitrary input and output angles17, spectral sorting of plane waves to distinct diffraction orders18, and light focusing with field-curvature correction19. While devices that are designed using topology optimization support enhanced efficiencies, they require extensive computational resources that dramatically increase with the size of the device. Hence, current topology-optimized devices are either microscale in dimension or limited to periodic structures with microscale unit cells.

    We report a conceptually new approach to metasurface design that extends the high performance of topology-optimized devices to macroscopic areas in a computationally efficient manner. The concept is schematically illustrated in Fig. 1b and is comprised of three parts. First, we discretize the desired phase profile into a series of wavelength-scale, linear sections. Next, we use topology optimization to design metasurface elements that optimally scatter light as required by each section of the phase profile. Finally, we stitch all the elements together to form a full metasurface. Our approach produces devices that operate more efficiently than conventional designs by accounting for and optimizing the near-field optical coupling between neighboring nanostructures. Our design method is also substantially more computationally efficient than existing topology optimization methods.

Results
  • To benchmark the improvements in computational efficiency that are afforded by our approach, we perform adjoint-based topology optimization on metagratings that are made of silicon ridges. These periodic metasurfaces are designed to deflect light of a specified wavelength λ to the +1 diffraction order. We perform electromagnetic simulations for topology optimization using rigorous coupled wave analysis (RCWA), which is also known as the Fourier modal method20-23, on a personal computer and plot the time that is required for optimizing a full device as a function of its width L in Fig. 1c. The simulation time scales approximately as ${\cal{O}}\left({L^{2.4}} \right)$, which corresponds to the general scaling trend for electromagnetic solvers that utilize standard matrix multiplication and inversion algorithms24. These trends demonstrate that prohibitively large computational resources are required for directly optimizing devices that are many times larger than the wavelength.

    If we instead divide the metagrating into sections of width d, we have L/d sections to optimize. The total computation time now scales as ${\cal{O}}\left({d^{2.4} \cdot L/d} \right) = {\cal{O}}\left({d^{1.4} \cdot L} \right)$, which is a linear function of the total size. This expression also shows that the computation time decreases with decreasing section size. In practice, we found that using 3λ-wide sections minimizes the total computation time, as the benefits from using smaller sections are outweighed by the computational overhead in the electromagnetic solver. We plot the computation time for a metasurface that has been optimized using 3λ-wide sections in Fig. 1c. The observed reduction in the computation time enables us to optimize devices with dimensions that are many times larger than the wavelength using realistic computational resources. For example, via our method, we can produce millimeter-scale topology-optimized metasurfaces that operate at near-infrared wavelengths in less than 1 day using a personal computer. In contrast, optimizing the full device at once would take nearly 1 year and require intractably large amounts of memory.

    We realize additional savings in computation time via our method when using multiple computing cores because individual segments can be optimized on different cores. Our design concept allows the segments to be treated independently, thereby enabling the optimizations to be parallelized without issues concerning race conditions and synchronization25. With N computing nodes, metasurfaces that are subdivided into N sections can be optimized in the same amount of time as it takes to optimize a single segment.

    These reductions in computational complexity also apply to fully three-dimensional topology-optimized metasurfaces, which exhibit even more severe scaling trends: the time it takes to optimize a metasurface of size L × L all at once scales approximately as ${\cal{O}}\left({L^{4.8}} \right)$. This trend indicates the necessity of linear sectioning for these more intricate design problems, which would reduce the runtime to a more reasonable ${\cal{O}}\left({L^2} \right)$.

  • The process of approximating a curvilinear phase profile with a series of linear sections introduces wavefront error, which is denoted as ϵ (Fig. 1b, inset). However, this error has a negligible impact on the overall metasurface performance if the sections are sufficiently small. To analyze the effect of linearizing a general curvilinear phase profile ϕ(x), we locally describe each section of ϕ(x) at location x0 using a 2nd-order Taylor series expansion as follows:

    $$ \phi \left( x \right) \approx \phi \left( {x_0} \right) + \phi ' \left( {x_0} \right)\left( {x - x_0} \right) + \frac{1}{2}\phi '' \left( {x_0} \right)\left( {x - x_0} \right)^2 $$ (1)

    This section can be approximated as a line of slope $\phi ' \left({x_0} \right)$ with a phase offset of ϕ(x0) + Δϕ, which incurs an error of ${\it{\epsilon }}\left(x \right) = \phi '' \left({x_0} \right)\left({x - x_0} \right)^2/2 - {\mathrm{\Delta }}\phi$. Given a section of length $d$ at position x0, the root-mean-square (RMS) wavefront error is minimized when ${\mathrm{\Delta }}\phi = \phi '' \left({x_0} \right)d^2/24$ and is expressed as follows:

    $$ {\it{\epsilon }}_{{{rms}}}\left( {{\mathrm{\Delta }}x} \right) = \frac{1}{{12\sqrt 5 }}\phi '' \left( {x_0} \right)d^2 $$ (2)

    We use this result to analyze the impact of linearization on a focusing cylindrical lens, which enables us to benchmark the device performance using well-established metrics in lens design. To quantify the performance of lenses that are constructed with linear sections, we use the Strehl ratio26, which is a metric that compares the diffraction efficiency of our lenses to that of an ideal lens. For a lens that focuses normally incident light, the ideal phase profile is $\phi \left(x \right) = \left({2\pi /\lambda } \right)({f - \sqrt {f^2 + x^2} })$27, where f is the focal length and λ is the wavelength. A lens with a Strehl ratio of 0.98, which corresponds to an RMS wavefront error of λ/50, can be realized if we use linear sections that are no larger than the following:

    $$ d\, < \, 0.73\sqrt {f\lambda } $$ (3)

    This equation provides a practical and quantitative guide for linearizing phase profiles in a manner that minimizes the phase error. As an example, consider a cylindrical lens with a focal length of f = 36λ and a numerical aperture (NA) of 0.7 (Fig. 2a). Using Eq. 3, we expect that linearizing the phase profile with segments that are smaller than 4.4λ will have a negligible impact on the performance. To verify, we simulate lenses that are linearized with various segment lengths and calculate the field intensities at the focal plane. These field intensity profiles are plotted in Fig. 2b, according to which the lenses that are linearized with section lengths of 2λ and 4λ are nearly indistinguishable from the ideal lens.

    Fig. 2  Impact of phase profile linearization on the cylindrical lens performance.

    a Cylindrical lenses of focal length 36λ and NA of 0.7 are constructed using linear sections of length d for various values of d. b Line scans of the field intensity at the focal planes of the lenses. Lenses that are linearized with sections that are smaller than 4λ have peak intensities that are within 1% of that of the ideal lens

    Our approach to sectioning readily extends to three-dimensional phase profiles, which can be approximated as series of planar tiles. This approximation is discussed in detail in the Supplementary Section. For a hyperboloid that is discretized into tiles with dimensions d × d, an RMS wavefront error of λ/50 can be realized if:

    $$ d\, < \, 0.61\sqrt {f\lambda } $$ (4)
  • To design metasurface elements that have the desired linear phase profiles, we utilize adjoint-based topology optimization28. Adjoint-based optimization is an iterative algorithm that modifies the device's dielectric constant distribution, namely, ε(x), to maximize a figure of merit (FoM). Our objective is to optimize a device that scatters normally incident electromagnetic waves in a desired direction with electric field amplitude Etgt and phase ϕtgt. To compute the FoM, we run a forward simulation in which waves that are incident onto the metasurface element scatter in the desired direction with field amplitude Efwd and phase ϕfwd. Near-to-far-field transformations from the forward simulations are used to evaluate Efwd and ϕfwd29. The FoM describes the difference between the current and desired responses and has the following form:

    $$ FoM = - A_1\left[ {\left| {\boldsymbol{E}}_{\text{tgt}} \right|^2 - \left| {\boldsymbol{E}}_{\text{fwd}} \right|^2} \right]^2 - A_2\left[ {\arg \left( {e^{i\left( {\phi _{\text{tgt}} - \phi _{\text{fwd}}} \right)}} \right)} \right]^2 $$ (5)

    The terms A1 and A2 are weights that balance how strongly the FoM is biased toward optimizing the amplitude and phase, respectively. To determine how ε(x) should be modified to improve the FoM each iteration, we perform a pair of forward and adjoint simulations and record the electric fields in the device for each excitation condition. These fields are used to calculate δFoM, which is the gradient of the FoM with respect to the dielectric constant at each position x:

    $$ \begin{array}{l}\delta FoM = {2}A_{1}\left( {\left| {\boldsymbol{E}}_{\text{tgt}} \right|^{2} - \left| {\boldsymbol{E}}_{\text{fwd}} \right|^2} \right){\cal{R}}e\left\{ {{\boldsymbol{E}}_{\text{fwd}} \cdot \delta {\boldsymbol{E}}^\ast } \right\}\\ - {2}A_{2}\left( \phi _{\text{tgt}} - \phi _{\text{fwd}} \right)\frac{1}{{\left| {\boldsymbol{E}}_{\text{fwd}} \right|^2}}{\cal{I}}m\left\{ {\boldsymbol{E}}_{\text{fwd}} \cdot \delta {\boldsymbol{E}^ \ast } \right\}\end{array} $$ (6)

    where $\delta {\boldsymbol{E}}$ is a function of the adjoint field and represents the variations of the field in the target direction in response to variations of the refractive index within the device. A more detailed discussion of the adjoint optimization method that is applied here is provided in the Supplementary Section.

  • To apply these concepts to the design of isolated, finite-sized device elements, we have developed an aperiodic Fourier modal method (AFMM), which is a hybrid method that combines a solver for periodic systems with perfectly matched layers (PMLs). The key challenge of implementing PMLs involves describing both the periodic incident plane wave (the input field) and the aperiodic scattered field (the output field) of the isolated device within the same formalism. To address this challenge, we introduce a hybrid method that combines a Fourier basis, Maxwell's equations in complex coordinates, and the Stratton-Chu integral formalism30-34.

    A metasurface is typically composed of a single layer of patterned material. The patterned material can be expressed as a distribution of the relative permeability, namely, ϵ(x, y), and the permeability, namely, μ(x, y), on the xy-plane. Along the thickness of the device in the z-direction, the device cross-section is constant. In this case, it can be shown from Maxwell's equations that the transverse electric fields satisfy the following eigenvalue equation:

    $$ - \gamma ^2\left[ {\begin{array}{*{20}{c}} {{\boldsymbol{E}}_x} \\ {{\boldsymbol{E}}_y} \end{array}} \right] = {\cal{L}}_{EH}{\cal{L}}_{HE}\left[ {\begin{array}{*{20}{c}} {{\boldsymbol{E}}_x} \\ {{\boldsymbol{E}}_y} \end{array}} \right] $$ (7)

    where ${\cal{L}}_{EH}$ and ${\cal{L}}_{HE}$ are differential operators that are based on ε and μ and the electric field solution can be expressed as ${{\Phi }}\left({x, y} \right)e^{ - ik\gamma z}$. The Fourier modal method can easily be used to expand the input fields, which are assumed to be periodic, into a Fourier basis:

    $$ {{\Phi }}_{in}\left( {x, y, z} \right) = \mathop {\sum }\limits_p A_pe^{ - ik\gamma _pz}\mathop {\sum }\limits_{nm} \delta _{nm, p}e^{ - ik\alpha _nx}e^{ - ik\beta _my} $$ (8)

    Next, the PML boundary conditions are described similarly. The PMLs can be introduced via a change of coordinates $\left({x, y, z} \right) \to (\tilde x, \tilde y, \tilde z)$30, 31:

    $$ \left\{ {\begin{array}{*{20}{c}} {\tilde x\left( x \right) = \left( {\chi _x - i\eta _x} \right)\left( {x - x_0} \right) + x_0} \\ {\tilde y\left( y \right) = \left( {\chi _y - i\eta _y} \right)\left( {y - y_0} \right) + y_0} \end{array}} \right\} $$ (9)

    Here, the parameter χ controls the scaling of the PML layers, while the parameter η controls the PML attenuation. This transformation is also useful for expressing the output scattered fields, which are computed as eigenfunctions of Maxwell's equations in complex coordinates:

    $$ {{\Phi }}_s\left( {x, y, z} \right) = \mathop {\sum }\limits_p B_pe^{ - ik\gamma _pz}\mathop {\sum }\limits_{nm} {{\Phi }}_{nm, p}e^{ - ik\alpha _nx}e^{ - ik\beta _my} $$ (10)

    Finally, the Stratton-Chu integral equation computes the radiated field in all of space34. In this study, we design devices that are comprised of nanoridges, which can be described in this analysis by a one-dimensional Fourier basis.

  • To demonstrate adjoint-based optimization in the design of metasurface elements, we design a 2.5λ-wide element that scatters incident TM-polarized light at a 20° angle with a phase response of π/2. The dielectric distribution and the amplitude and phase scattering profiles of the element at successive iterations in the optimization process are shown in Fig. 3a. The optimization begins with a random dielectric continuum with values between air and silicon. After a few iterations, the continuum begins to strongly scatter light at the desired angle. The final metasurface section is a binary structure of silicon in air that possesses a peak scattering amplitude and phase response matching with the targeted values. The full width at half maximum of the scattering peak is consistent with that expected from light that is diffracting from a 2.5λ-wide aperture; hence, the element is performing directional scattering near its physical limits.

    Fig. 3  Metasurface element optimization.

    a At the beginning of the optimization process, the initial dielectric distribution is a random dielectric continuum. The dashed lines indicate the desired scattering angle of 20° and phase of π/2. After 10 iterations, the scattering profile is already highly directional. After 100 iterations, the optimization is complete and the metasurface section is a binary structure of silicon and air that supports the desired scattering metrics. b An intensity plot that shows the scattered fields of the optimized metasurface element

    An intensity plot of the scattered fields from the fully optimized section (Fig. 3b) shows strong near-field coupling between neighboring nanostructures; hence, optimal near-field coupling is responsible for mediating strong scattering in the desired direction. As a method of gradient descent, adjoint-based topology optimization is a local optimizer and is sensitive to the initial dielectric distribution10. To obtain high-performance elements for a desired scattering angle and phase target, we perform ten optimizations with various initial dielectric distributions and select the best result.

  • Combining everything into a proof-of-concept demonstration, we stitch together optimized metasurface elements to construct cylindrical metalenses. We design these metalenses to focus TM-polarized light at a wavelength of 640 nm. To enable device operation at visible wavelengths, we use 250-nm-thick crystalline silicon, which has relatively low absorption compared to polycrystalline and amorphous silicon, but much higher index contrast than materials such as titanium dioxide35, 36.

    First, we design and simulate 64-μm-wide metalenses with NAs that range from 0.2 to 0.9. We divide the metalenses into sections that are 2 μm wide, which is below the phase error limit of Eq. 3 and near the optimal size for efficient computation (Fig. 1c). Further reductions in section size lead to degradation of the device performance. The reason can be traced to our design of each metasurface element, which is optimized in isolation with PML boundary conditions. When the elements are stitched together to produce a device, the optical fields that are guided by a single element have evanescent tails and can couple to a neighboring element in a parasitic manner. Smaller section sizes require more elements to be stitched together to produce a desired metasurface, thereby resulting in more boundaries and more parasitic coupling. Below a section size of 2λ, the device performance begins to degrade and below a section size of 1λ, the aperiodic boundary conditions in our optimizer are no longer valid.

    There are a few approaches for addressing the issue of stitching error. One is to keep the section size relatively large compared to the wavelength. Another is to perform boundary optimization on the stitched regions to eliminate the stitching error. A third approach, which we use here, is to separate silicon structures from other sections by a gap of at least 0.2λ, thereby reducing the near-field coupling between sections. To ensure a reduction in the stitching error with this scheme, we simulate stitched sections to check for spurious diffraction and redesign the sections in the event of excess error.

    The efficiencies of our simulated metalenses are summarized in Fig. 4a. The absolute efficiency is defined as the amount of power that is contained in the principal lobe of the focus compared to that of an ideal lens with 100% transmission. The relative efficiency, or focusing efficiency, compares the power in the principal lobe to that of an ideal lens that transmits the same amount of power as the device37. This efficiency corresponds to the efficiency of the diffraction process, as it removes the effects of absorption from the material and reflection at the metalens interface.

    Fig. 4  Simulation results of optimized metalenses.

    a Relative and absolute efficiencies for metalenses that were designed with various numerical apertures. b A full-field simulation of a metalens with an NA of 0.9

    The efficiency plots show that the relative efficiencies are consistently high, namely, above 93%, with minimal drop-off in performance as the NA increases. This trend is unlike that of conventional metalenses, where the efficiency decreases with increasing NA because conventional designs cannot efficiently deflect light at large angles27, 38. The absolute efficiencies of the metalenses all exceed 75%, with approximately 10% of the light reflected from the metalens and 10% absorbed by the silicon. Reflection losses can be reduced via the use of more intricate three-dimensional silicon nanostructures, while absorption losses can be minimized by designing silicon-based devices for longer wavelength operation28. A simulated field profile of a metalens with an NA of 0.9 is shown in Fig. 4b, which demonstrates that the lens focuses strongly with minimal spurious diffraction.

  • We design, fabricate, and characterize 200-μm-wide metalenses with NAs of 0.2, 0.5, and 0.8. To prepare crystalline silicon thin films on glass, we use hydrogen silsesquioxane to bond silicon-on-insulator wafers onto Pyrex wafers under high temperature and pressure36, 39. After removing the silicon handle wafer and the buried oxide layer, we pattern and etch the devices via standard electron beam lithography and dry etching techniques. We characterize the metalenses by collimating polarized, monochromatic light from a tunable white-light laser onto the devices and imaging the light at the focal plane with a ×100 objective (NA = 0.9) and a CCD sensor.

    Scanning electron microscope images of the center of a representative device are shown in Fig. 5a and show silicon nanostructures that exhibit smooth and vertical sidewalls. The metalenses all have relative efficiencies that exceed 89% and absolute efficiencies that exceed 67% for all NAs, which are within 10% of the simulated values (Fig. 5b). All the metalenses exhibit diffraction-limited performance, as shown by the theoretical and experimental intensity plots in Fig. 5c-e. The device with an NA of 0.8 can focus light at a wavelength of 640 nm to a spot with a beam waist of 340 nm. The central lobes of the foci are all much stronger than the side lobes; hence, the focusing efficiency is high.

    Fig. 5  Experimental characterization of fabricated metalenses.

    a Scanning electron micrographs of the metalenses with tilted and top-down views. b Relative and absolute efficiencies of the fabricated metalenses, along with their simulated values. c-e Intensity line scans at the focal planes of the metalenses with NAs of 0.2, 0.5, and 0.8, respectively, along with comparisons with the simulated lenses. f-h Efficiencies of the fabricated metalenses as a function of the wavelength, along with their simulated values. These metalenses have the same design but were fabricated on a different sample than those that were measured in (a-e)

    The metalenses maintain reasonably high efficiencies for wavelengths that range from 580 to 700 nm, as shown in Fig. 5f-h. As these lenses are not designed to be achromatic, the focal length shifts with the wavelength. The simulated shifts are shown in the Supplementary section in Fig. S3. Future work will focus on generalizing our design approach to include achromatic functionality, which can be addressed by modifying the optimizer's figure of merit to include multiple wavelengths18. The figure of merit can be specified so that each metasurface section deflects all wavelengths in the same direction and realizes the correct dispersion for ensuring constructive interference at the focus40, 41.

Discussion
  • In summary, we present a computationally efficient method for designing large-area, topology-optimized metasurfaces. Our approach of approximating a desired phase profile with a series of linear segments renders the design problem computationally tractable while minimally impacting the device performance. Wavelength-scale scattering elements that are designed via topology optimization can strongly scatter light at a desired angle and phase, thereby serving as high-performance building blocks for larger metasurfaces.

    The metalenses that are demonstrated here have limited functionality and other heuristic design methods are available that can produce high-performance, high-numerical-aperture metalenses42. What differentiates our method is that our design platform can be generalized systematically to multifunctional, broadband, and even multilayer metasurfaces with large areas and high efficiencies. To practically extend topology optimization in these directions, continued improvements in electromagnetic solver speed will be required, particularly for fully three-dimensional devices. For multilayer devices, new innovations will be needed to address the stitching error, which can compound due to interactions between neighboring elements. We envision that these design methods will enable the realization of compact and monolithic electromagnetic devices that exhibit high efficiencies and practical functionality across the electromagnetic spectrum.

Materials and methods
  • Crystalline silicon on glass wafers was prepared via an approach that is similar to that done by Sell, et al.36. A Pyrex and silicon-on-insulator (SOI) wafer were bonded under high temperature and pressure using hydrogen silsesquioxane. The silicon handle wafer was etched away using SF6 plasma and the buried oxide layer was etched away using HF solution. The wafer was diced into pieces for later use. The metasurfaces were patterned onto a piece in the AR-P-6200 resist via electron beam lithography. Aluminum oxide was evaporated onto a piece and removed using liftoff to form a hard mask. The silicon was etched using a Cl2 and HBr plasma. The aluminum oxide was removed using a solution of HCl at 50 ℃. Additional details can be found in the supplementary information.

  • A collimated beam of light from a tunable white light laser is filtered with a longpass filter and polarized with a linear polarizer. Then, the beam passes through a weak cylindrical lens to focus it onto the metasurface aperture without significantly changing the angle of incidence. The metasurface focuses the light. Next, the focal plane is imaged onto a CMOS detector using a 0.9 NA, 100x objective and a tube lens.

    To evaluate the lens efficiencies, a line scan of the focal plane image is obtained. The lens relative efficiency is calculated by integrating over a small aperture around the main focal lobe and comparing it to what would be expected for an ideal lens of the same numerical aperture. The lens absolute efficiency is calculated by multiplying the relative efficiency by the measured transmission. A schematic diagram of the optical setup and additional characterization details can be found in the supplementary information.

  • We used Reticolo, a rigorous coupled-wave analysis solver, as a basis to perform optimizations21. The final metasurface lens designs were simulated using Lumerical FDTD. The efficiencies were calculated by determining the optical power that is contained in the main lobe of the focus and comparing it to that expected from an ideal lens.

Acknowledgements
  • This work was supported by the U.S. Air Force under Award Number FA9550- 18-1-0070, the Office of Naval Research under Award Number N00014-16-1- 2630, and the David and Lucile Packard Foundation. DS was supported by the National Science Foundation (NSF) through an NSF Graduate Research Fellowship. SD was supported by the Department of Defense (DoD) through an NDSEG Research Fellowship. This work was performed in part at the Stanford Nanofabrication Facility (SNF) and the Stanford Nano Shared Facilities (SNSF), which are supported by the National Science Foundation as part of the National Nanotechnology Coordinated Infrastructure under award ECCS- 1542152.

Author' contributions
  • T.P., D.S., K.E. and J.Y. formulated the design algorithms. T.P., D.S., E.W. and S.D. fabricated the devices. D.S. performed the measurements and analyzed the data. T.P. and J.F. wrote the manuscript with input from all authors.

Conflict of interest
  • The authors declare that they have no conflict of interest.

Supplementary information
Reference (42)

Catalog

    /

    DownLoad:  Full-Size Img PowerPoint
    Return
    Return