pith. sign in

arxiv: 2604.06110 · v1 · submitted 2026-04-07 · ⚛️ physics.optics

Refractive Index Robustness of Metalenses

Dongyoung Lee , Jisoo Kyoung This is my paper

Pith reviewed 2026-05-10 18:51 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords metalensesrefractive index robustnessphase-reset boundarieswavefront engineeringdielectric metasurfacesfocal stabilityenvironmental perturbations
0
0 comments X

The pith

Dielectric metalenses keep their focal spot stable when refractive index varies because the locations of their 2π phase-reset boundaries do not move.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that dielectric metalenses maintain nearly constant best-focus position and spot size across a wide range of refractive index values. This occurs because index changes deform the phase only locally while leaving the physical positions of the 2π phase-reset boundaries fixed, so the overall wavefront gradient stays the same. Readers should care because real optical systems experience index shifts from temperature or other environmental factors, and this built-in stability could reduce the need for active compensation in compact metalens designs. The work also finds that the same boundary invariance produces quasi-scale-invariant focusing even when material dispersion is present under uniform scaling.

Core claim

Parametric sweeps of refractive index reveal that the metalens focal profile shows negligible deviation in position and size. The origin of this robustness is the structural invariance of the zone boundaries: the spatial locations of the 2π phase-reset boundaries remain stationary, preserving the effective wavefront gradient despite local phase deformations induced by the index change. This invariance further supports a quasi-scale-invariant focusing behavior that follows a predictable linear trend under uniform geometric scaling even in the presence of material dispersion.

What carries the argument

The stationary spatial locations of the 2π phase-reset boundaries that keep the effective wavefront gradient unchanged when local phases shift with refractive index.

Load-bearing premise

Local phase changes from index variation do not move the physical positions of the 2π phase-reset boundaries or add unmodeled scattering, absorption, or fabrication effects that would alter the zone structure.

What would settle it

A measured shift in best-focus position or spot size when refractive index is varied in a fabricated metalens while holding geometry fixed would falsify the claimed robustness.

Figures

Figures reproduced from arXiv: 2604.06110 by Dongyoung Lee, Jisoo Kyoung.

Figure 3
Figure 3. Figure 3: Refractive-index-dependent modification of the meta-atom phase response and the resulting metalens phase profile. (a) Transmission phase as a function of pillar diameter 𝒅 for several refractive indices around the nominal value 𝒏𝟎 = 𝟐. 𝟎𝟑𝟗𝟒, showing that the phase response varies systematically with refractive index. (b) For 𝒏 > 𝒏𝟎 , comparison of the phase profile at 𝒏 = 𝒏𝟎 , the expected trend under refr… view at source ↗
Figure 5
Figure 5. Figure 5: (a) Geometric comparison between the Fresnel lens and the metalens. (b) Magnified view illustrating the alignment [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
read the original abstract

Metalenses have emerged as a powerful platform for compact wavefront engineering; however, their performance stability under refractive index fluctuations induced by environmental perturbations, such as temperature shifts, remains a critical concern. Here, we demonstrate the intrinsic refractive index robustness of dielectric metalenses and elucidate its physical origin. By parametrically sweeping the refractive index, we observe that the metalens maintains a stable focal profile with negligible deviations in best-focus position and spot size over a broad range of variations. We identify that this robustness arises from the structural invariance of the zone boundaries: despite index-induced local phase deformations, the spatial locations of the 2{\pi} phase-reset boundaries remain stationary, thereby maintaining the effective wavefront gradient. Furthermore, we reveal that this robustness enables a "quasi-scale-invariant" focusing behavior, where the focusing performance follows a predictable linear trend under uniform geometric scaling even in the presence of material dispersion. Our findings suggest that metalenses can maintain stable focusing behavior against refractive index variations that may arise from unavoidable environmental perturbations in practical optical systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper claims to demonstrate the intrinsic refractive index robustness of dielectric metalenses. Through parametric sweeps of the refractive index, the focal profile remains stable with negligible deviations in best-focus position and spot size. The physical origin is identified as the stationarity of the 2π phase-reset boundaries, which maintains the effective wavefront gradient despite local phase deformations from index changes. It further reveals a quasi-scale-invariant focusing behavior under uniform geometric scaling even with material dispersion.

Significance. If substantiated, this result would be of high significance for the optics community, particularly in designing metalenses resilient to environmental variations such as temperature fluctuations affecting refractive index. The elucidation of the boundary invariance as the mechanism provides a valuable physical insight, and the quasi-scale-invariant property could inform scaling laws in dispersive media. The use of parametric sweeps to observe the effect is a strength.

major comments (1)
  1. [Abstract and physical origin discussion] The assertion that 'the spatial locations of the 2π phase-reset boundaries remain stationary' (as stated in the abstract) is central to the robustness claim. However, this appears to assume a uniform phase response to δn across meta-atoms. In reality, for dielectric pillars, n_eff is diameter-dependent due to mode confinement, so δφ/δn varies spatially. This would distort the phase profile and likely shift the reset boundaries or change local slopes, as per the relation φ(d) = (2π/λ)(n_eff(d)−1)h. The manuscript should clarify if the sweeps use geometry-specific phase maps from full-wave simulations or an approximate model, and show explicit evidence (e.g., phase profiles before/after index shift) that boundaries do not move.
minor comments (2)
  1. The abstract lacks quantitative details such as the range of index variations swept, error bars on focal position and spot size, and specific simulation parameters or software used.
  2. Consider adding references to previous studies on metalens sensitivity to index or temperature variations for context.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and constructive comments on our manuscript. We address the major comment point-by-point below. We have revised the manuscript to include additional clarifications and evidence as requested.

read point-by-point responses

  1. Referee: [Abstract and physical origin discussion] The assertion that 'the spatial locations of the 2π phase-reset boundaries remain stationary' (as stated in the abstract) is central to the robustness claim. However, this appears to assume a uniform phase response to δn across meta-atoms. In reality, for dielectric pillars, n_eff is diameter-dependent due to mode confinement, so δφ/δn varies spatially. This would distort the phase profile and likely shift the reset boundaries or change local slopes, as per the relation φ(d) = (2π/λ)(n_eff(d)−1)h. The manuscript should clarify if the sweeps use geometry-specific phase maps from full-wave simulations or an approximate model, and show explicit evidence (e.g., phase profiles before/after index shift) that boundaries do not move.

    Authors: We appreciate the referee highlighting this potential subtlety in the physical mechanism. The refractive index sweeps were indeed performed by recalculating the transmission phase for each meta-atom using full-wave simulations with the updated refractive index, employing geometry-specific phase maps rather than a uniform approximation. Although δφ/δn does vary with meta-atom diameter due to mode confinement, our simulations reveal that the 2π phase-reset boundaries remain stationary. This is because the meta-atom diameters are selected to match the required wrapped phase profile at the design index, and the resulting phase deformations preserve the locations where the phase crosses multiples of 2π. To provide explicit evidence, we have added phase profile comparisons in the revised manuscript (see new Supplementary Figure S2), which demonstrate the invariance of the boundary positions before and after the index variation. We have also updated the abstract and the physical origin discussion to clarify the use of full-wave simulations and reference this evidence. revision: yes

Circularity Check

0 steps flagged

No circularity: robustness derived from direct parametric observation of geometric invariance

full rationale

The paper derives refractive index robustness by parametrically sweeping the material index in simulations and directly observing that focal profiles remain stable while 2π phase-reset boundary locations stay fixed despite local phase shifts. This stationarity is presented as an emergent geometric property of the zone structure, verified independently of the focal metrics themselves. No equations reduce the result to a fitted parameter renamed as prediction, no self-citations bear the central load, and no ansatz or uniqueness theorem is smuggled in; the claim rests on the simulation outputs rather than self-referential definition.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard electromagnetic assumptions for phase calculation in dielectrics without introducing new free parameters or postulated entities; robustness emerges from geometry rather than fitted constants.

axioms (1)
  • standard math Electromagnetic wave propagation follows standard phase accumulation rules in lossless dielectrics
    Invoked when relating refractive index to local phase shift and zone boundaries.

pith-pipeline@v0.9.0 · 5467 in / 1155 out tokens · 45482 ms · 2026-05-10T18:51:28.093929+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages

  1. [1]

    The future and promise of flat optics: a personal perspective,

    F. Capasso, "The future and promise of flat optics: a personal perspective," Nanophotonics 7(6), 953–957 (2018)

    work page 2018

  2. [2]

    Review for optical metalens based on metasurfaces: fabrication and applications,

    Z. Hu, M. Gu, Y. Tian, C. Li, M. Zhu, H. Zhou, B. Fang, Z. Hong, and X. Jing, "Review for optical metalens based on metasurfaces: fabrication and applications," Microsyst Nanoeng 11(1), 189 (2025)

    work page 2025

  3. [3]

    On-chip integration of achromatic metalens arrays,

    Y. Zhang, X. Jiang, G. Qu, J. Han, C. Li, B. Bo, Q. Ruan, Z. Liu, Q. Song, and S. Xiao, "On-chip integration of achromatic metalens arrays," Nat Commun 16(1), 7485 (2025)

    work page 2025

  4. [4]

    Metalens-integrated compact imaging devices for wide-field microscopy,

    B. Xu, H. Li, S. Gao, X. Hua, C. Yang, C. Chen, F. Yan, S. N. Zhu, and T. Li, "Metalens-integrated compact imaging devices for wide-field microscopy," AP 2(6), 066004 (2020)

    work page 2020

  5. [5]

    Large area metalenses: design, characterization, and mass manufact uring,

    A. She, S. Zhang, S. Shian, D. R. Clarke, and F. Capasso, "Large area metalenses: design, characterization, and mass manufact uring," Opt. Express, OE 26(2), 1573–1585 (2018)

    work page 2018

  6. [6]

    Towards high -throughput large-area metalens fabrication using UV-nanoimprint lithography and Bosch deep reactive ion etching,

    C. A. Dirdal, G. U. Jensen, H. Angelskår, P. C. V. Thrane, J. Gjessing, and D. A. Ordnung, "Towards high -throughput large-area metalens fabrication using UV-nanoimprint lithography and Bosch deep reactive ion etching," Opt. Express, OE 28(10), 15542–15561 (2020)

    work page 2020

  7. [7]

    Wafer-Scale Functional Metasurfaces for Mid-Infrared Photonics and Biosensing,

    A. Leitis, M. L. Tseng, A. John-Herpin, Y. S. Kivshar, and H. Altug, "Wafer-Scale Functional Metasurfaces for Mid-Infrared Photonics and Biosensing," Advanced Materials 33(43), 2102232 (2021)

    work page 2021

  8. [8]

    Large-area metasurface on CMOS-compatible fabrication platform: driving flat optics from lab to fab,

    N. Li, Z. Xu, Y. Dong, T. Hu, Q. Zhong, Y. H. Fu, S. Zhu, and N. Singh, "Large-area metasurface on CMOS-compatible fabrication platform: driving flat optics from lab to fab," Nanophotonics 9(10), 3071–3087 (2020)

    work page 2020

  9. [9]

    Production -Ready Double-Sided Fabrication of Dual -Band Infrared Metaoptics Using Deep-Ultraviolet Lithography,

    K. Sun, X. Yan, J. Scott, J. -Y. Ou, J. N. Monks, and O. L. Muskens, "Production -Ready Double-Sided Fabrication of Dual -Band Infrared Metaoptics Using Deep-Ultraviolet Lithography," ACS Nano 19(42), 37218–37229 (2025)

    work page 2025

  10. [10]

    A Metalens with a Near-Unity Numerical Aperture,

    R. Paniagua-Domí nguez, Y. F. Yu, E. Khaidarov, S. Choi, V. Leong, R. M. Bakker, X. Liang, Y. H. Fu, V. Valuckas, L. A. Krivitsky, and A. I. Kuznetsov, "A Metalens with a Near-Unity Numerical Aperture," Nano Lett. 18(3), 2124–2132 (2018)

    work page 2018

  11. [11]

    Ultrahigh Numerical Aperture Metalens at Visible Wavelengths,

    H. Liang, Q. Lin, X. Xie, Q. Sun, Y. Wang, L. Zhou, L. Liu, X. Yu, J. Zhou, T. F. Krauss, and J. Li, "Ultrahigh Numerical Aperture Metalens at Visible Wavelengths," Nano Lett. 18(7), 4460–4466 (2018)

    work page 2018

  12. [12]

    High -efficiency high-numerical-aperture metalens designed by maximizing the efficiency limit,

    S. Li, H.-C. Lin, and C. W. Hsu, "High -efficiency high-numerical-aperture metalens designed by maximizing the efficiency limit," Optica, OPTICA 11(4), 454–459 (2024)

    work page 2024

  13. [13]

    High-performance achromatic flat lens by multiplexing meta-atoms on a stepwise phase dispersion compensation layer,

    J. Lin, J. Chen, J. Zhang, H. Liang, J. Li, and X.-H. Wang, "High-performance achromatic flat lens by multiplexing meta-atoms on a stepwise phase dispersion compensation layer," Light Sci Appl 14(1), 110 (2025)

    work page 2025

  14. [14]

    Engineering metallic nanostructures for plasmonics and nanophotonics,

    N. C. Lindquist, P. Nagpal, K. M. McPeak, D. J. Norris, and S. -H. Oh, "Engineering metallic nanostructures for plasmonics and nanophotonics," Rep. Prog. Phys. 75(3), 036501 (2012)

    work page 2012

  15. [15]

    Recent advances in planar optics: from plasmonic to diel ectric metasurfaces,

    P. Genevet, F. Capasso, F. Aieta, M. Khorasaninejad, and R. Devlin, "Recent advances in planar optics: from plasmonic to diel ectric metasurfaces," Optica, OPTICA 4(1), 139–152 (2017)

    work page 2017

  16. [16]

    Metalenses at visible wavelengths: Diffractio n-limited focusing and subwavelength resolution imaging,

    M. Khorasaninejad, W. T. Chen, R. C. Devlin, J. Oh, A. Y. Zhu, and F. Capasso, "Metalenses at visible wavelengths: Diffractio n-limited focusing and subwavelength resolution imaging," Science 352(6290), 1190–1194 (2016)

    work page 2016

  17. [17]

    Metasurface eyepiece for augmented reality,

    G.-Y. Lee, J.-Y. Hong, S. Hwang, S. Moon, H. Kang, S. Jeon, H. Kim, J.-H. Jeong, and B. Lee, "Metasurface eyepiece for augmented reality," Nat Commun 9(1), 4562 (2018)

    work page 2018

  18. [18]

    Compact folded metasurface spectrometer,

    M. Faraji-Dana, E. Arbabi, A. Arbabi, S. M. Kamali, H. Kwon, and A. Faraon, "Compact folded metasurface spectrometer," Nat Commun 9(1), 4196 (2018)

    work page 2018

  19. [19]

    Metamaterials and Metasurfaces for Sensor Applications,

    Y. Lee, S.-J. Kim, H. Park, and B. Lee, "Metamaterials and Metasurfaces for Sensor Applications," Sensors 17(8), (2017)

    work page 2017

  20. [20]

    Role of refractive index in metalens performance,

    E. Bayati, A. Zhan, S. Colburn, M. V. Zhelyeznyakov, and A. Majumdar, "Role of refractive index in metalens performance," App l. Opt., AO 58(6), 1460–1466 (2019)

    work page 2019

  21. [21]

    Silica Meta-Optics: When High Performance Does Not Need a High Index,

    L. Sacchi, A. Palmieri, V. Mishra, J.-S. Park, M. Piccardo, and F. Capasso, "Silica Meta-Optics: When High Performance Does Not Need a High Index," Nano Lett. 25(50), 17448–17457 (2025)

    work page 2025

  22. [22]

    Thermo -optic coefficient of silicon at 1550 nm and cryogenic temperatures,

    J. Komma, C. Schwarz, G. Hofmann, D. Heinert, and R. Nawrodt, "Thermo -optic coefficient of silicon at 1550 nm and cryogenic temperatures," Appl. Phys. Lett. 101(4), 041905 (2012)

    work page 2012

  23. [23]

    Measurements of the refractive indices and thermo -optic coefficients of Si 3N4 and SiO x using microring resonances,

    A. Arbabi and L. L. Goddard, "Measurements of the refractive indices and thermo -optic coefficients of Si 3N4 and SiO x using microring resonances," Optics Letters, Vol. 38, Issue 19, pp. 3878-3881 (2013)

    work page 2013

  24. [24]

    Photoelastic Constants of Vitreous Silica and Its Elastic Coefficient of Refractive Index,

    W. Primak and D. Post, "Photoelastic Constants of Vitreous Silica and Its Elastic Coefficient of Refractive Index," J. Appl. Phys. 30(5), 779– 788 (1959)

    work page 1959

  25. [25]

    Stress optic coefficient and stress profile in optical fibers,

    N. Lagakos, R. Mohr, and O. H. El-Bayoumi, "Stress optic coefficient and stress profile in optical fibers," Appl. Opt., AO 20(13), 2309–2313 (1981)

    work page 1981

  26. [26]

    Sensitive measurement of optical nonlinearities using a single beam,

    M. Sheik-Bahae, A. A. Said, T.-H. Wei, D. J. Hagan, and E. W. Van Stryland, "Sensitive measurement of optical nonlinearities using a single beam," IEEE Journal of Quantum Electronics 26(4), 760–769 (1990)

    work page 1990

  27. [27]

    Dispersion of bound electron nonlinear refraction in solids,

    M. Sheik-Bahae, D. C. Hutchings, D. J. Hagan, and E. W. Van Stryland, "Dispersion of bound electron nonlinear refraction in solids," IEEE Journal of Quantum Electronics 27(6), 1296–1309 (1991)

    work page 1991

  28. [28]

    Thermal effects on metalenses,

    D. Lee and J. Kyoung, "Thermal effects on metalenses," J. Appl. Phys. 137(11), 113101 (2025)

    work page 2025

  29. [29]

    Multiwavelength polarization -insensitive lenses based on dielectric metasurfaces with meta-molecules,

    E. Arbabi, A. Arbabi, S. M. Kamali, Y. Horie, and A. Faraon, "Multiwavelength polarization -insensitive lenses based on dielectric metasurfaces with meta-molecules," Optica, OPTICA 3(6), 628–633 (2016)

    work page 2016

  30. [30]

    Broadband mid -infrared frequency comb generation in a Si3N4 microresonator,

    K. Luke, Y. Okawachi, M. R. E. Lamont, A. L. Gaeta, and M. Lipson, "Broadband mid -infrared frequency comb generation in a Si3N4 microresonator," Opt. Lett., OL 40(21), 4823–4826 (2015)

    work page 2015

  31. [31]

    Highly -efficient full-color holographic movie based on silicon nitride metasurface,

    M. Yamaguchi, H. Saito, S. Ikezawa, and K. Iwami, "Highly -efficient full-color holographic movie based on silicon nitride metasurface," Nanophotonics 13(8), 1425–1433 (2024)

    work page 2024

  32. [32]

    Hecht, Optics, Fifth edition (Pearson, 2017)

    E. Hecht, Optics, Fifth edition (Pearson, 2017)

    work page 2017

  33. [33]

    M. Born, E. Wolf, and A. B. Bhatia, Principles of Optics: Electromagnetic Theory of Propagation, Interference, and Diffraction of Light , Seventh (expanded) anniversary edition, 60th anniversary edition (Cambridge University Press, 2019)

    work page 2019

  34. [34]

    Influence of temperature on diffractive lens performance,

    G. P. Behrmann and J. P. Bowen, "Influence of temperature on diffractive lens performance," Appl. Opt., AO 32(14), 2483–2489 (1993)

    work page 1993

  35. [35]

    Photonic Crystals: Molding the Flow of Light - Second Edition,

    J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, "Photonic Crystals: Molding the Flow of Light - Second Edition," in (Princeton University Press, 2008)

    work page 2008

  36. [36]

    J. A. Stratton, Electromagnetic Theory, IEEE Press Series on Electromagnetic Wave Theory (Wiley, 2007)

    work page 2007

  37. [37]

    Chromatic dispersion in photonic crystal fibers: fast and accurate scheme for calculation,

    J. Læ gsgaard, A. Bjarklev, and S. E. B. Libori, "Chromatic dispersion in photonic crystal fibers: fast and accurate scheme for calculation," J. Opt. Soc. Am. B, JOSAB 20(3), 443–448 (2003)

    work page 2003

  38. [38]

    Improving image control by correcting the lens-heating focus drift,

    B.-J. Cheng, H.-C. Liu, Y. Cui, and J. Guo, "Improving image control by correcting the lens-heating focus drift," in Optical Microlithography XIII (SPIE, 2000), 4000, pp. 818–826

    work page 2000

  39. [39]

    4H -SiC Metalens: Mitigating Thermal Drift Effect in High -Power Laser Irradiation,

    B. Chen, X. Sun, X. Li, L. Cai, D. Zhao, K. Du, M. Pan, and M. Qiu, "4H -SiC Metalens: Mitigating Thermal Drift Effect in High -Power Laser Irradiation," Advanced Materials 37(3), 2412414 (2025). Supplementary information

    work page 2025

  40. [40]

    Refractive Index-Dependent Phase Profiles of a Convex-Plano Refractive Lens Fig. S1. Wrapped phase profiles of a convex –plano refractive lens for refractive indices 𝑛 = 1.8, 2.0394, and 2.2. As the refractive index varies, the phase gradient changes, while the overall phase curvature remains similar. To provide a comparative reference, we analyzed a conv...

    work page

  41. [41]

    Additional NA cases Fig. S2. Additional NA cases. (a,b) High -NA metalens (𝐷 = 50 μm, 𝑓 ≈ 20 μm, 𝑁𝐴 ≈ 0.78) and (c,d) Low -NA metalens (𝐷 = 50 μm, 𝑓 ≈ 160 μm, 𝑁𝐴 ≈ 0.15). (a,c ) best-focus position and spot size (FWHM) versus refractive index 𝑛. (b,d) normalized x-z intensity map at 𝑛 = 1.8,2.0394, and 2.2 (red dashed lines: best-focus positions). To asse...

    work page

  42. [42]

    We define the scaling factor 𝑠 relative to the baseline design at 𝑠 = 1

    Scaling protocol for quasi-scale-invariant focusing This section summarizes the definition of uniform geometric scaling used in the main text and clarifies the scope of parameters that were jointly scaled in the FDTD simulations. We define the scaling factor 𝑠 relative to the baseline design at 𝑠 = 1. The wavelength and all length-type geometric parameter...

    work page

  43. [43]

    Geometry (design): For the metalens case, period of unit cell, height of Si3N4 pillars, lens diameter 𝐷, diameter of meta-atoms 𝑑, thickness of substrate; for the refractive lens case, lens diameter 𝐷, center thickness, and radius of curvature

    work page

  44. [44]

    Simulation conditions: FDTD region span 𝐿𝑥, 𝐿𝑦, 𝐿𝑧, size and position of field monitor for near -field data, size and position of source (source-to-metalens distance), z-scan range for best-focus search, z-scan step size, focal-plane monitor size, and mesh cell size (set relative to 𝜆(𝑠))

    work page

  45. [45]

    The table below lists representative parameters jointly adjusted under the scaling applied in this work

    Total simulation time. The table below lists representative parameters jointly adjusted under the scaling applied in this work. Table S3. Summary of representative scaled parameters. Scale 𝑠 Wavelength 𝜆 [nm] Lens diameter [μm] Refractive index of SiO2 substrate Height of Si3N4 pillar [μm] Period of unit cell [μm] 0.5 316.5 25 1.4835 0.75 0.17 1 633 50 1....

    work page