{"id":"2622fbc4-2a05-425e-b09b-58618b5a5f21","arxiv_id":"2501.07979","paper_version":6,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new symmetry-based domain decomposition, with rotational symmetry order increasing outward in radial zones, enables full-wave topology optimization of metalenses up to 3000 wavelengths in diameter.","lead":"The authors introduce a computational design method, zoned discrete axisymmetry, that cuts the cost of full-wave metalens optimization from scaling with lens area to scaling roughly with diameter. They use it to topology-optimize millimeter- and centimeter-scale 3D freeform metalenses, including color-corrected visible lenses, with simulated efficiencies above published state of the art.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zone-independence approximation is not validated for the high-NA, large-aperture designs whose efficiencies anchor the state-of-the-art claim.","rationale":"I read the paper in good faith and find the scaling argument for ZDA (n growing with radius, constant-per-zone wedge area, linear total cost) to be internally coherent and a genuine contribution. The central risk is not the scalability math but whether the independently simulated zones, when assembled, reproduce the field of the real monolithic lens. Every reported efficiency is a zone-composited simulation result, so all state-of-the-art comparisons inherit this approximation. The prior work cited for the <1% error bound does not obviously cover the high-NA visible regime or the very large number of zones used here, and no full-lens validation is presented. This matches the reader's weakest assumption exactly. I recommend keeping the CONDITIONAL verdict: the method and scaling claim are plausible, but the performance claims should be treated as conditional pending a monolithic full-wave check or experiment. The concrete test proposed is feasible for a downscaled aperture and would directly settle whether the zone-independence assumption introduces a material error.","tokens_in":12830,"tokens_out":11233,"duration_ms":121858,"concrete_test":"Downscale the RGB design to a diameter where a full-lens 3D FDTD simulation is feasible (e.g., 100λ–200λ aperture while keeping zone widths at 25λ), re-optimize or simply re-use the corresponding inner zones, and compare the full-lens simulated PSF and absolute focusing efficiencies at 488/532/658 nm against the zone-composited values. If the efficiencies differ by more than a few percent (or the PSF shape changes materially), the zone-independence assumption is not adequate for high-NA visible designs and the reported 33.1% average must be treated as an upper bound rather than a validated state-of-the-art result. A cheaper adjoint test is to simulate two adjacent zones together without the intervening PML and quantify the field difference in each zone.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The ZDA method's reported efficiencies (58.4% LWIR, 33.1% RGB, 12% six-wavelength) all come from simulations in which each radial zone is solved independently with PML boundaries, then assembled by far-field summation. The paper imports the validity of this decomposition from Ref. 13 ('less than 1% error ... when utilizing >10λ-wide zones'), but Ref. 13's test was not made at NA=0.8 in the visible, and the present designs use zones 25λ–47λ wide with thicknesses of only 0.6λ–1λ, where the physical argument for locality (Refs. 24–25) applies, but the quantitative error bound is not re-established here. There is no full-lens simulation or experiment anywhere in the paper showing that the assembled PSF or efficiency matches a monolithic simulation. Since the headline 'outperforms the state of the art' is based on these zone-composited efficiencies, an overestimate of even a few percent in the assembled focusing efficiency—or a PSF change from inter-zone coupling—would directly weaken the central claim. The scaling claim is internally sound; the weak point is the performance validation path.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces \"zoned discrete axisymmetry\" (ZDA) as a symmetry-reduction and domain-decomposition strategy for full-wave inverse design of large-area metalenses. The lens is divided into concentric supra-wavelength annular zones, and the azimuthal symmetry order n is increased roughly linearly with radius so that the simulated wedge area per zone is nearly independent of radius and the total computational cost scales almost linearly with lens diameter. Each zone is simulated independently in a GPU-accelerated cylindrical-coordinate FDTD solver with Bloch boundary conditions and PML terminations, and a differentiable adjoint/automatic-differentiation framework is used for freeform topology optimization. The authors report three designs: a 0.8-cm-diameter LWIR metalens (NA=0.3, 58.4% absolute focusing efficiency), a 1.05-mm RGB-achromatic metalens (NA=0.8, 33.1% average efficiency), and a 1.96-mm six-wavelength metalens (NA=0.3, 12% average efficiency), and they argue that these results outperform state-of-the-art metalenses. Permittivity profiles for all designs are provided as supplementary files.","tokens_in":13051,"tokens_out":14662,"duration_ms":159152,"significance":"If the ZDA approximation is valid, this is a substantial methodological contribution: it offers a route to full-wave, non-LPA topology optimization of millimeter- and centimeter-scale metalenses with near-linear cost scaling, and it provides concrete, reproducible designs with stated efficiencies. The GPU cylindrical-FDTD implementation and the hybrid time/frequency-domain adjoint-AD integration are valuable building blocks for the community. The main caveat is that the reported performance numbers are produced by the same zone-decomposed PML-terminated solver used for optimization, and the accuracy of this decomposition is not independently verified for the specific high-NA, wide-aperture parameter ranges. The significance of the efficiency claims therefore depends on closing that validation gap.","major_comments":[{"comment":"The headline efficiencies (58.4% LWIR, 33.1% RGB, 12% six-wavelength) are all computed from the same zone-decomposed, PML-terminated solver that generated the designs. The accuracy of the decomposition is imported from Ref. 13, whose <1% error bound is cited for >10λ-wide zones, but the conditions of that test are not reproduced here. The present designs use zone widths of 25λ and 47λ, thicknesses of roughly 0.6λ–1λ, and NAs up to 0.8, and the assembled lenses contain abrupt changes in the symmetry order n at zone boundaries that are not represented when each zone is independently terminated by PMLs. No full-lens simulation or experimental measurement is reported for any assembled design. Because the \"outperform the state of the art\" claim rests on these efficiency values, an independent validation—for example, a full-wave simulation of a smaller lens with the same NA and zone widths, or a fabricated and measured sub-aperture—is needed before the central claim is established.","section":"Zoned discrete axisymmetry (ZDA); Results and Discussion"},{"comment":"The state-of-the-art comparison in Table 1 is not yet on a common footing. The absolute focusing efficiency is defined in the LWIR section as the fraction of total incident power within 3 FWHMs, but the RGB and six-wavelength efficiencies are reported without restating this definition, and the table does not define the column header R or the normalization (incident vs. transmitted power). Literature entries in Table 1 use a variety of efficiency definitions, so the assertion that the present designs outperform the state of the art needs a side-by-side comparison with the same integration aperture and normalization for each cited work. At minimum, please specify the integration radius in wavelengths for every design and state whether each value is absolute or relative efficiency.","section":"Table 1; Millimeter-scale poly-achromatic metalenses in the visible"},{"comment":"The claim that discrete axisymmetry is \"critical\" and that the performance gain over continuous axisymmetry is due to the enlarged design space is supported by a single comparison (58.4% vs. 25.1%). The text states that the two designs share diameter, NA, and materials, but it does not report the optimization schedules (number of iterations, filter/binarization ramps, starting points) for the two runs. Without evidence that both designs were optimized to comparable convergence, the attribution of the 2.3x efficiency gap specifically to the additional azimuthal degrees of freedom is not fully established.","section":"Centimeter-scale metalens at long-wave infrared"}],"minor_comments":[{"comment":"The column header R is not defined; from the text it appears to denote lens radius rather than diameter, and the caption should state this explicitly.","section":"Table 1"},{"comment":"The implementation of the continuously axisymmetric first zone (n approaching infinity) is not described; please clarify how this limit is represented in the Bloch-mode FDTD solver and how the m=±1 excitation is handled near r=0.","section":"Differentiable FDTD in cylindrical coordinates"},{"comment":"The term \"full-wave\" is used for a method that still relies on the zone-independence approximation; please qualify it (e.g., \"full-wave within each zone\") so that readers do not infer that the assembled lens is simulated without any domain-decomposition approximation.","section":"Introduction and Abstract"},{"comment":"There are two occurrences of \"for for\" in this subsection; please correct the typos.","section":"Millimeter-scale poly-achromatic metalenses in the visible"},{"comment":"The paper does not state how the FDTD throughput of one billion voxels per second was measured; a sentence describing the benchmark procedure would improve reproducibility.","section":"Centimeter-scale metalens at long-wave infrared"},{"comment":"The text refers to the optimized variables as \"pixels\" and to the designs as \"3D freeform\" and \"volumetric\"; please clarify explicitly whether the permittivity varies along z within the metalens layer or whether the design is a 2D patterning in (r, phi) with fixed thickness.","section":"Results and Discussion"}],"recommendation":"major_revision","confidential_remarks":"This is a promising manuscript with a genuinely interesting algorithmic idea and concrete, reproducible designs. The main obstacle to acceptance is the absence of an independent accuracy check for the zone-decomposition approximation at the specific parameter ranges used for the headline efficiencies. I recommend requesting a validation study—either a full-wave simulation of a downscaled lens with the same NA and zone widths, or an experimental PSF/efficiency measurement of a fabricated sub-aperture—and a standardized efficiency comparison with the cited state-of-the-art works. The self-citation to Ref. 13 is appropriate, but the quantitative transfer of its error bound to the present settings needs to be demonstrated in this paper rather than assumed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nThe thing to know: the ZDA construction is a real algorithmic idea. By letting the rotational symmetry order n grow with radius, each radial zone keeps a nearly constant simulation area, so full-wave FDTD-based topology optimization scales roughly linearly with diameter instead of quadratically. The paper works out the Bloch-mode cylindrical FDTD, the r=0 regularization, and the choice of n(ℓ) to suppress higher azimuthal diffraction orders, then demonstrates three large designs: 0.8 cm LWIR at 58.4% simulated efficiency, 1.05 mm RGB at 33.1%, and 1.96 mm six-wavelength at 12%. That is a legitimate step beyond continuous axisymmetry and beyond locally periodic approximation for freeform design.\n\nCredit where due: the comparison against a continuous-axisymmetry design for the LWIR case (25.1% vs 58.4%) is a nice control showing the extra azimuthal DOFs matter. The paper also ships the permittivity data for the designed lenses as supporting information, which is more than most papers in this area do. The area-counting scaling argument is internally consistent.\n\nThe soft spot is the validation path. Every headline efficiency is computed with the same zoned/PML solver that produced the design, and the zone-independence assumption is carried over from the authors' earlier overlapping-domains paper (ref 13), which reported <1% error only for >10λ-wide zones in a different regime. Here the visible RGB lens uses 25λ-wide zones at NA=0.8 and thickness around one wavelength; the locality argument from Miller and Li-Hsu bounds nonlocal spread by thickness, but the quantitative error bound is not re-established for this regime. No full-lens simulation, no experiment, and no independent reconstruction of the assembled PSF appear anywhere. Also, the adjoint-AD implementation is deferred to \"elsewhere,\" so the method is not fully reproducible from the text—only the final designs are. The \"outperform the state of the art\" claim is fair only if read as \"state-of-the-art simulated efficiencies under the ZDA model,\" not as demonstrated performance.\n\nWho benefits: anyone working on large-area metalens inverse design. The ZDA idea is likely to be picked up and tested independently. The paper deserves a serious referee, with the key request being either a full-lens simulation for at least the smaller RGB design or a careful error-bounded justification of the zone decomposition at high NA. My own verdict is conditional until that validation appears.\n\nRecommendation: send to peer review.","headline":"Zoned discrete axisymmetry is a genuine linear-scaling full-wave design idea, but the headline efficiencies rest on an unvalidated zone-independence assumption and need independent confirmation before the state-of-the-art claim is taken at face value.","tokens_in":13593,"tokens_out":2630,"would_cite":true,"duration_ms":27377,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Zoned discrete axisymmetry makes full-wave inverse design of centimeter-scale metalenses scale nearly linearly with diameter, without locally periodic approximations.","keywords":["zoned discrete axisymmetry","metalens inverse design","topology optimization","full-wave FDTD","freeform meta-optics","achromatic metalens","GPU-accelerated simulation","computational imaging"],"falsifier":"Take one of the published designs (for example the 0.8-cm LWIR lens), run a full-lens simulation or experimental measurement with all zones present, and compare the focal-spot efficiency and point-spread function with the zone-decomposed prediction; a discrepancy beyond a few percent in focusing efficiency, or a visible change in focal position, would show the independent-zone assumption fails at the relevant numerical aperture.","tokens_in":12642,"feed_emoji":"🔬","tokens_out":6983,"duration_ms":66278,"temperature":0.7,"pith_summary":"This paper introduces \"zoned discrete axisymmetry\" (ZDA), a way to split a wide-aperture metalens into a few concentric radial zones and give each zone its own degree of rotational symmetry, increasing outward. The central claim is that this makes full-wave, freeform inverse design scale nearly linearly with lens diameter, instead of with area, while keeping the full 3D design freedom that continuous axisymmetry sacrifices. On that basis the authors report topology-optimized designs of a 0.8-cm long-wave-infrared lens with 58.4% absolute focusing efficiency, a 1.05-mm RGB-achromatic lens with 33.1% average efficiency, and a 1.96-mm six-wavelength lens with 12% average efficiency, all exceeding their chosen state-of-the-art baselines. A sympathetic reader would care because full-wave Maxwell optimization at millimeter and centimeter scales was previously thought impractical without approximations such as the locally periodic approximation.","feed_headline":"Full-wave metalens design now scales nearly linearly with diameter","feed_subtitle":"Zone-by-zone full-wave solver designs cm-scale lenses that beat locally periodic approximations.","key_machinery":"The central mechanism is \"zoned discrete axisymmetry\" (ZDA): the lens is divided into concentric radial zones of width much larger than the wavelength, and the ℓ-th zone is made periodic in azimuth with period $2\\pi/n_{\\ell}$, where $n_{\\ell}$ grows roughly linearly with radius. Each zone's simulation is a cylindrical-coordinate wedge terminated by absorbing layers, with Bloch mode expansions in the azimuthal coordinate; a differentiable FDTD solver with adjoint-based gradients then performs topology optimization inside each zone independently. This object carries the argument because it reduces the total computational domain from area-scaling to diameter-scaling while retaining genuinely three-dimensional freeform patterns.","core_discovery":"Using a GPU-accelerated finite-difference time-domain solver in cylindrical coordinates, the paper states that each supra-wavelength radial zone can be simulated independently with absorbing boundaries and azimuthal Bloch periodicity, so the simulated volume (and hence computational cost) grows roughly linearly with radius. The n-fold symmetry order n is chosen to grow approximately linearly with radius, keeping the azimuthal period subwavelength and suppressing spurious diffraction orders while adding about ten times more degrees of freedom than continuous axisymmetry. The paper's demonstration designs—a 0.8-cm LWIR lens at NA 0.3, a 1.05-mm RGB lens at NA 0.8, and a 1.96-mm six-wavelength lens at NA 0.3—are claimed to outperform prior large-area metalenses, and the LWIR design beats its continuously axisymmetric counterpart (58.4% vs 25.1% focusing efficiency).","pith_inferences":["An implication the authors leave implicit: if inter-zone independence is as clean as claimed, the method could be parallelized across zones almost without limit, so wall-clock time should drop further with more GPUs.","A testable extension would apply ZDA to an end-to-end imaging objective rather than focal-spot intensity; the freeform, non-Fresnel patterns typical of computational imaging are precisely where LPA is least reliable, so ZDA's advantage should be largest there.","The same zone-with-growing-n construction could be adapted to discrete rotational symmetries for polarization or orbital-angular-momentum sensing, though the paper does not demonstrate those devices."],"forward_implications":["Full-wave inverse design becomes practical for metalenses thousands of wavelengths across, without invoking the locally periodic approximation.","Adding azimuthal degrees of freedom raises achievable focusing efficiency: in the paper's comparison, the discrete-axisymmetry LWIR lens reaches 58.4% versus 25.1% for the continuous-axisymmetry lens under the same parameters.","Achromatic and poly-achromatic visible metalenses with diameters around 1600–3000 wavelengths can be designed at average focusing efficiencies of 33% and 12%, outperforming LPA-based designs at similar diameters.","All target wavelengths are extracted from one forward and one adjoint time-domain simulation, so adding colors costs little extra compute.","The optimized lenses can produce super-oscillatory focal spots (FWHM 1.34λ versus 1.6λ for the ideal Airy disk), which could aid computational super-resolution imaging."],"supporting_citations":[{"why":"Establishes the supra-wavelength zone-decomposition premise: PML-terminated zones wider than 10λ give far-field error below 1% versus full simulation.","marker":"[13]"},{"why":"Supplies the fullwave Maxwell inverse-design method and axisymmetric metalens baseline that ZDA generalizes.","marker":"[14]"},{"why":"Provides the continuous-axisymmetry fullwave metalens framework whose design space ZDA enlarges.","marker":"[17]"},{"why":"Introduces the ad hoc periodic azimuthal variation that ZDA replaces with rigorous fullwave optimization.","marker":"[8]"},{"why":"State-of-the-art LPA-based inverse-designed large-scale meta-optics baseline to which the visible RGB and poly-achromatic designs are compared.","marker":"[5]"},{"why":"RGB-achromatic meta-optics baseline for visible metalenses, used as a performance comparison.","marker":"[4]"},{"why":"Provides the FDTD formulation (staggered Yee grid, leapfrog update, stretched-coordinate PMLs) used by the solver.","marker":"[26]"},{"why":"Supplies the hybrid time/frequency-domain adjoint sensitivity analysis that makes freeform gradient-based optimization possible.","marker":"[28]"},{"why":"Underpins the GPU-accelerated simulation throughput used to estimate speedup potential.","marker":"[23]"}],"fun_headline_variants":["Metalens design goes linear: full-wave, zone by zone","Zoned axisymmetry makes full-wave metalens design scale","Full-wave metalens design beats LPA, scales linearly","Zone-by-zone full-wave solver designs cm-scale lenses","Full-wave metalenses: near-linear scaling with diameter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes that adjacent radial zones interact only weakly through the lens, so each zone can be optimized on its own with artificial absorbing edges; if near-field coupling across zone boundaries is actually significant, the assembled lens will not perform as the zone-by-zone simulations predict.","fun_headline_variants_meta":{"raw":{"variants":["Metalens design goes linear: full-wave, zone by zone","Zoned axisymmetry makes full-wave metalens design scale","Full-wave metalens design beats LPA, scales linearly","Zone-by-zone full-wave solver designs cm-scale lenses","Full-wave metalenses: near-linear scaling with diameter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001343,"raw_usage":{"total_tokens":5441,"prompt_tokens":915,"completion_tokens":4526,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":531,"completion_tokens_details":{"reasoning_tokens":4454}},"tokens_in":531,"tokens_out":4526,"duration_ms":35571,"temperature":1.0,"reasoning_tokens":4454,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:31:19.098708+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take one of the published designs (for example the 0.8-cm LWIR lens), run a full-lens simulation or experimental measurement with all zones present, and compare the focal-spot efficiency and point-spread function with the zone-decomposed prediction; a discrepancy beyond a few percent in focusing efficiency, or a visible change in focal position, would show the independent-zone assumption fails at the relevant numerical aperture.","supporting_citations":[],"review_version":1}