{"id":"64bc4df8-1a94-4896-a0fe-61e42bea3238","arxiv_id":"2607.13843","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new metric finds operating points where gyrotropic rods scatter maximally nonreciprocally, and arrays of such rods show strongly direction-dependent transmission near diffraction-order onsets.","lead":"This paper uses a cylindrical-wave expansion to quantify and maximize nonreciprocal scattering from magnetically biased plasmonic rods, then arranges them into a grating that transmits differently from opposite sides. A generalist should read it as a design recipe for compact, tunable optical isolators and directional routers that does not yet include an independent numerical validation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The peak collective nonreciprocity is computed at Rayleigh-Wood anomalies with a dipole-truncated grating model that is never checked against an independent full-wave solver; without that check the headline asymmetry may be an artifact of the truncation.","rationale":"The reader's conditional verdict already identifies the dipole truncation, the lossless assumption, and the missing full-wave validation as the key risks. My stress test sharpens this into a single load-bearing concern: the collective effect is claimed specifically at Rayleigh-Wood anomalies, precisely the regime where the truncated dipole model is most fragile and where no independent numerical check exists. The analytic framework is internally coherent: the metric vanishes when ωc=0, energy-conservation sums are close to unity, and the Wood-anomaly mechanism is physically plausible. However, this does not lower the risk. The final field plots in Figs. 9 and 11 appear to be computed from the same approximate Eq. (14), not from a different solver, so they cannot serve as validation. The promised full-wave simulations are absent from the manuscript. The proposed concrete test—a finite-element or RCWA simulation of the two representative designs, with and without damping—directly settles whether the near-anomaly asymmetry survives beyond the dipole approximation. Because the concern is testable and does not by itself prove the result wrong, the appropriate verdict remains conditional: the paper should not be accepted as an unqualified claim of maximum nonreciprocity until that check is performed. Hence I recommend keeping the reader's verdict unchanged.","tokens_in":18569,"tokens_out":7964,"duration_ms":80518,"concrete_test":"Run an independent full-wave simulation (e.g., finite-element or rigorous coupled-wave analysis) of the two selected structures: a=0.07λp, b=8.8a, θ=17°, λ=1.47λp and a=0.07λp, b=10.5a, θ=56°, λ=1.47λp, using the Drude tensor of Eqs. (2)-(3) with ωc=0.05ωp. Extract τ0, τ−1, ρ0, ρ−1 for both illumination sides and compute ΔP from Eq. (17). Compare to Figs. 8 and 10. Repeat with a realistic collision frequency (e.g., ν=0.01ωp). If ΔP drops below ~0.3 or the peaks shift away from the Wood anomaly, the central claim needs qualification. Also rerun the analytical model with multipoles |u|≤2 to check convergence of ΔP near the dashed lines.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that optimized gyrotropic rod arrays transmit light in a totally different way under opposite illumination, with the strongest collective nonreciprocity occurring near the onset of new diffraction orders (Section III.B, Figs. 7-10). The load-bearing assumption is that the dipole-truncated grating model (Eqs. 5, 11-14) remains quantitatively accurate exactly in the near-anomaly regime where ΔP peaks. This is the least secure point for three reasons.\n\nFirst, the high-ΔP regions in Figs. 7(c,d) hug the dashed lines where m=±1 orders become grazing. Equation (14) contains a factor 1/κ_m; as κ_m→0 at the Rayleigh-Wood anomaly, the analytic continuation is delicate. The authors explicitly avoid the exact boundaries because \"our overall approximate approach... faces numerical issues\" (Section III.B), so the maxima they report are inferred from a model they admit is unreliable in the immediately relevant limit. Second, the dipole truncation (u=0,±1) is used both to solve the grating boundary problem and to produce the \"simulated\" field maps in Figs. 9 and 11; those are evaluations of Eq. (14), not independent full-wave simulations. The multipolar check in Fig. 6 uses the same approximate internal field (Eq. 5), so it does not provide independent verification. The abstract and introduction promise \"full-wave numerical simulations,\" but no such comparison appears anywhere in the text. Third, the Drude model is lossless (no collision frequency). At a Wood anomaly, the transition from evanescent to propagating is singular in the lossless limit; even small material damping can smooth that transition and reduce the claimed contrast between opposite sides. Since the effect relies on one order being just barely propagating and the other just barely evanescent, damping is not a minor quantitative correction. Together, these issues make the headline asymmetry a prediction of the approximate model in a regime where that model is least trustworthy and completely un","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an analytical cylindrical-wave framework for TE scattering by gyrotropic (magnetically biased) Drude rods. A single-particle nonreciprocity metric ΔS′ = |S′₁ + S′₋₁| is introduced, maximized over frequency, cyclotron frequency, and rod radius, and linked to a resonance condition in Eq. (7). The optimized rods are then arranged into an infinite periodic grating, and a dipole-truncated multiple-scattering model (Eqs. 11–14) is used to compute direction-dependent transmission/reflection, with emphasis on peaks of the collective nonreciprocity metric ΔP near the onset of new diffraction orders. Field maps and spectral curves for two representative designs (a = 0.07λp, b = 8.8a and b = 10.5a) are presented as evidence of strongly asymmetric response under opposite illumination.","tokens_in":19012,"tokens_out":12385,"duration_ms":126202,"significance":"If the results hold, the work offers a potentially useful design route for all-passive, magnetically tunable nonreciprocal metasurfaces. The analytical treatment of a single gyrotropic cylinder, the explicit optimization over physical parameters, and the mapping of collective nonreciprocity to Wood-anomaly conditions are valuable contributions. The paper also provides transparent definitions of metrics and a detailed multipolar decomposition, which are helpful for physical interpretation. However, the central claims rest on an approximate dipole-truncated grating model in a regime where the authors admit the model is unreliable, and the promised full-wave validation is absent. These issues currently prevent the results from being considered quantitatively established.","major_comments":[{"comment":"Equation (7) is numerically inconsistent as printed. For the operating point ω/ωp = 0.7 and ωc/ωp = 0.05 used throughout the paper, the left-hand side equals approximately 1.53, while the right-hand side is 2.5×10⁻³. The equation therefore cannot describe the blue lines in Fig. 2, which are claimed to follow from (7). Since the single-rod optimization selects operating points based on these resonance loci, this is a load-bearing error. The authors must correct the equation, provide its derivation, and confirm that the blue lines in Fig. 2 correspond to the corrected condition.","section":"Section II.B, Eq. (7)"},{"comment":"The abstract and Section I promise 'full-wave numerical simulations' validating the nonreciprocal response. No such comparison appears in the text. The field maps in Figs. 9 and 11 are evaluations of Eq. (14), not independent full-wave solutions. Moreover, in Section III.B the authors state that the 'overall approximate approach ... faces numerical issues' near the actual peaks of ΔP, yet the headline claim is that peak collective nonreciprocity coincides with the emergence of new diffraction orders. This is precisely the Rayleigh–Wood-anomaly regime where the dipole-truncated model is least secure. An independent full-wave validation, or at least a convergence study with higher azimuthal harmonics, is required to support the central claim.","section":"Section III.B and III.C"},{"comment":"The multipolar decomposition used to argue that electric and magnetic dipole modes dominate is computed from the same approximate internal field (Eq. 5) that underlies the scattering coefficients. Figure 6 therefore does not independently verify the dipole truncation; it restates it. The conclusion that quadrupole contributions are weak is built into the ansatz rather than demonstrated. A comparison with a full-wave solver or a higher-order multipole expansion (e.g., including u = ±2) is needed to confirm that the truncation is quantitatively accurate in the parameter ranges used.","section":"Section II.C, Fig. 6"},{"comment":"The optimization and demonstration are not fully independent. The rods are selected to maximize |S′₁ + S′₋₁|, which directly creates asymmetric scattering patterns; the grating built from these rods then inherits a degree of directional asymmetry. The collective ΔP metric in Eq. (17) is a different observable, but the connection would be more convincing if the authors compared the optimized rods with non-optimized, less asymmetric rods under the same grating parameters, or showed that the collective enhancement exceeds what the single-particle asymmetry alone would produce. This would strengthen the claim that the near-anomaly behavior is a collective effect rather than simply a restatement of the selection criterion.","section":"Section II.B and III.B"}],"minor_comments":[{"comment":"Typos: 'truely' in the Introduction; 'quadrapole' in Section II.C; 'th other' and 'nonlinearity' in Section III.B (should be 'the other' and 'nonreciprocity'). The phrase 'blows up' for Eq. (7) is misleading since ΔS′ is bounded by unity.","section":"Various"},{"comment":"The four polar curves are not labeled with their corresponding rod radii or frequencies. Adding a legend or explicit labels would improve readability.","section":"Fig. 3"},{"comment":"Captions for Fig. 4 and Fig. 7 use a/λ0 while the text uses a/λp; the notation should be made consistent.","section":"Captions"},{"comment":"The statement that 'nonreciprocity gets, on average, boosted' and the phrase 'increasing the maximum nonlinearity' in the discussion of Fig. 7(c) are imprecise; the intended term is likely 'nonreciprocity'.","section":"Section III.B"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely topic and the authors have relevant prior work in the area. However, the numerical inconsistency in Eq. (7) and the absence of the promised full-wave validation are serious. The apology in Section III.B about numerical issues near the very peaks that constitute the main result further weakens confidence. I believe the work can be made publishable if the equation is corrected, full-wave or higher-order convergence evidence is supplied, and the claims are scaled back where the approximation is known to fail."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a competent analytical design study. What's genuinely new: the delta-S' metric, the systematic maps over frequency, cyclotron frequency, and rod radius, and the observation that peak collective nonreciprocity sits near Wood-anomaly boundaries where new diffraction orders emerge. The grating derivation via addition theorem and Poisson summation is standard but not trivial, and the multipolar decomposition supports the dipole picture. Those parts I trust far more than the headline claim.\n\nNow the soft spots, in proportion. First, the abstract and introduction promise full-wave numerical simulations, but none exist in the text. Figures 9 and 11 are field maps evaluated from the same dipole-truncated grating model, Eq. (14), and the multipolar check in Fig. 6 uses the same internal-field approximation, Eq. (5). So the central claim—opposite illumination gives a totally different transmission—is a prediction of one approximate model with no independent verification. This matters because the claimed maxima hug the diffraction-order cutoffs where kappa_m goes to zero and the model itself is admitted to face numerical issues. A separate full-wave solver could easily confirm or refute the contrast. Second, Eq. (7) as printed is numerically inconsistent: at omega/omega_p = 0.7 and omega_c/omega_p = 0.05, the left side is about 1.5 while the right side is 0.0025. Probably a typo, but a reader cannot reproduce the blue curves without a correction or derivation. Third, the Drude model is lossless; near a Wood anomaly the transition from evanescent to propagating is singular in the lossless limit, so material damping is not a minor perturbation to the claimed asymmetry.\n\nI don't think these flaws sink the paper. The metric and optimization maps are reusable, the grating model is clearly stated, and the authors do check power conservation. But the load-bearing \"totally different\" claim should not be repeated as a validated result until an independent full-wave comparison appears. The right referee is someone who runs periodic-solver simulations, and they should be asked to test Fig. 7(d), Fig. 9, and Fig. 11 with a commercial or in-house code, once for lossless and once for lossy Drude.\n\nRecommendation: send to peer review, but with a required full-wave validation section and a correction or derivation of Eq. (7). It is a solid analytical contribution to the nonreciprocal metasurface subfield. I'd bring it to a relevant reading group and would cite it for the metric and the design rule, not for the quantitative contrast predictions.","headline":"Useful analytical design study for gyrotropic rod metasurfaces, but it promises full-wave validation it never delivers, and the central near-anomaly claims rest on the same approximate dipole-truncated model; worth refereeing, needs a real numerical check plus a fix for Eq. (7).","tokens_in":19545,"tokens_out":2444,"would_cite":true,"duration_ms":38382,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.25.Fx","78.20.Ls"],"model":"deepseek-v4-flash","headline":"This paper claims that optimized gyrotropic rods in a periodic array can make an all-passive metasurface transmit light in a totally different way when excited from opposite sides, with the strongest nonreciprocity occurring at the onset of","keywords":["nonreciprocity","gyrotropic rods","metasurface","magneto-optical","time-reversal symmetry breaking","diffraction orders","Wood's anomaly","magnetic bias"],"falsifier":"A full-wave numerical simulation that includes Drude losses and does not truncate the internal field at u = ±1, checking whether the metasurface's transmission from opposite sides still shows the predicted near-anomaly reversal; if the asymmetry largely disappears or the peak moves away from the Wood's anomaly thresholds, the central claim would be falsified.","tokens_in":18509,"feed_emoji":"🔄","tokens_out":4303,"duration_ms":45289,"temperature":0.7,"pith_summary":"The paper develops an analytical framework to design subwavelength all-passive metasurfaces that are strongly nonreciprocal, meaning they transmit light differently when illuminated from opposite sides. For individual gyrotropic rods, it introduces a metric based on the asymmetry between azimuthal dipole modes and shows that maximum nonreciprocity occurs when one dipole coefficient dominates. When these optimized rods are arranged into a periodic grating, the collective nonreciprocity peaks at the thresholds where new diffraction orders emerge, because one illumination direction sees the order as propagating while the other sees it as bound. If correct, this provides a path to tunable flat-optics isolators and directional routers that require no active or nonlinear components, only an external magnetic bias.","feed_headline":"Magnetic bias makes rod grating transmit directionally","feed_subtitle":"Peak nonreciprocity appears when a new diffraction order emerges, opening the path to compact isolators.","key_machinery":"The key machinery is the nonreciprocity metric ΔS' for individual rods, which quantifies the asymmetry between the +1 and −1 azimuthal scattering coefficients, and the collective metric ΔP for the metasurface, which compares transmitted and reflected power across all diffraction orders for opposite illumination sides. The individual-rod analysis uses a dipolar approximation of the internal field, keeping only azimuthal orders u = 0 and u = ±1, and yields a closed-form condition for maximum nonreciprocity. For the array, the grating diffraction orders are computed using the addition theorem and Poisson summation, and the mechanism for enhanced nonreciprocity is the evanescent-to-propagating t","core_discovery":"The central discovery is that magnetically biased plasmonic rods, each small relative to the wavelength, can exhibit a strong directional response quantified by the metric ΔS' = |S'_1 + S'_−1|, which vanishes without magnetic bias and reaches its maximum value of one when only one of the two azimuthal dipole coefficients survives. The authors derive a resonance condition that locates these maxima and use a multipolar decomposition to show the effect is rooted in the overlapping electric and magnetic dipole modes of opposite handedness. In the periodic metasurface, the paper shows that the transmission and reflection spectra change drastically depending on which side the wave arrives from, an","pith_inferences":["The predicted extreme contrasts rely on the lossless Drude plasma assumption; with realistic plasmonic damping, the sharp near-anomaly features would likely broaden and the transmission asymmetry would be reduced, though the qualitative directionality may persist.","The alignment of peak nonreciprocity with Wood's anomaly suggests a general design strategy: any periodic nonreciprocal scatterer, not just gyrotropic rods, may exhibit maximal directional contrast at diffraction-order thresholds, and this could be tested with other meta-atom geometries.","A direct experimental test would compare transmitted power from opposite sides of a prototype array at wavelengths near the predicted threshold; the asymmetry should appear as a sharp function of wavelength and incidence angle, and could be used to validate the design.","The same individual-rod metric could be extended to more complex meta-atoms, such as core-shell or elliptical gyrotropic rods, where geometric asymmetry might further enhance the nonreciprocal response."],"forward_implications":["All-passive subwavelength metasurfaces can act as direction-dependent transmitters or isolators without nonlinear or active components, controlled by an external magnetic field.","A clear design rule emerges: operate near the threshold for the m = −1 or m = −2 diffraction order to maximize collective nonreciprocity.","The individual-rod metric provides a fast screening tool for choosing rod radius, frequency, and magnetic bias to achieve maximal directional scattering.","Reversing the magnetic bias sign swaps the favored diffraction channels, effectively acting as a dynamic switch.","The findings point to compact isolators, directional transceivers, and wavefront routers in flat-optics configurations."],"fun_headline_variants":["Gyrotropic rods unlock max nonreciprocity in metasurfaces","Magnetic bias flips rod grating's transmission direction","Subwavelength rods achieve peak one-way light flow","Directional metasurface from magnetically biased rods","Rod metasurface hits max nonreciprocity with magnetic bias"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The analysis assumes a lossless Drude plasma and truncates the internal field to dipolar orders u = 0 and ±1; if real losses or higher-order multipoles are significant, the predicted maximum nonreciprocity and the sharp contrast near diffraction-order thresholds would change, and the paper does not include the full-wave validation promised in its abstract.","fun_headline_variants_meta":{"raw":{"variants":["Gyrotropic rods unlock max nonreciprocity in metasurfaces","Magnetic bias flips rod grating's transmission direction","Subwavelength rods achieve peak one-way light flow","Directional metasurface from magnetically biased rods","Rod metasurface hits max nonreciprocity with magnetic bias"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000151,"raw_usage":{"total_tokens":1033,"prompt_tokens":735,"completion_tokens":298,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":217}},"tokens_in":479,"tokens_out":298,"duration_ms":26496,"temperature":1.0,"reasoning_tokens":217,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T03:33:11.740878+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A full-wave numerical simulation that includes Drude losses and does not truncate the internal field at u = ±1, checking whether the metasurface's transmission from opposite sides still shows the predicted near-anomaly reversal; if the asymmetry largely disappears or the peak moves away from the Wood's anomaly thresholds, the central claim would be falsified.","supporting_citations":[],"review_version":1}