{"id":"0a6c0c65-0488-483e-adb1-bd3d7890f9c3","arxiv_id":"2607.06517","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"Cylindrical-harmonic expansions plus Graf translations, solved by a four-level adaptive high-precision solver, yield verified high-accuracy scattering solutions for arbitrary non-overlapping parallel circular cylinders including dense subwavelength clusters.","lead":"A semi-analytical multipole method with adaptive high-precision linear algebra solves 2D EM scattering by arbitrary clusters of parallel circular cylinders to controlled accuracy, even for dense subwavelength gaps. It supplies benchmark near- and far-field quantities for plasmonic nanowire clusters that grid methods struggle to resolve.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The multipole–Graf construction is classical; the paper’s contribution is the carefully engineered adaptive linear-algebra layer and the multi-test verification suite that together make the method a reliable benchmark for denser numerical schemes. The reader’s identification of the g\to0 limitation is accurate but already disclosed and does not undermine the claim for finite gaps. Condition-number growth is monitored and countered by extended-precision / exact arithmetic, and forward-error control is demonstrated by independent physical tests rather than by residual alone. No load-bearing derivation gap, circularity, or untested regime that would reverse the ACCEPT verdict was found. Independent re-implementation remains desirable for further confidence but is not required to sustain the present assessment.","tokens_in":28689,"tokens_out":436,"duration_ms":5790,"concrete_test":"Independently re-implement the assembly of the symmetry-reduced block matrix A for the aluminum trimer (Table I parameters, N=18) and solve at solver level 3 (mp.dps=120); recompute the optical-theorem residual and the three independent C_sca values of Sec. VII E. If all relative discrepancies remain ≲10^{-6}, the strongest claim is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's weakest_assumption correctly flags the excluded g=0 singularity and the need for N to grow as gaps shrink, but this is already stated and bounded by the paper (Sec. II D, V A, Table V). The central claim—that finite multipole truncation plus adaptive high-precision solution of the truncated system yields controlled, multistage-verified accuracy for non-touching configurations, including densely packed subwavelength ones—is supported by the residual, optical-theorem, boundary-condition, energy-balance, and symmetry tests of Sec. VII E. No internal inconsistency or unacknowledged gap in the derivation of the block system (5.1)–(5.2) or the four-level solver appears. The geometric restriction to circular cylinders is an explicit scope limit, not a hidden flaw.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript presents a semi-analytical method for 2D electromagnetic scattering by an arbitrary finite ensemble of parallel, non-overlapping homogeneous circular cylinders of arbitrary radii, complex permittivities, and transverse positions. The field is expanded in cylindrical harmonics about each cylinder axis; multiple scattering is closed via Graf’s addition theorem, yielding a block linear system for the modal coefficients (Eqs. 3.18, 5.1–5.2). The system is truncated at multipole order N and solved by a four-level adaptive procedure (float64 LU, equilibration with iterative refinement, mpmath high precision, and exact elimination over Gaussian rationals), with condition-number monitoring. The method is demonstrated on a subwavelength aluminum nanotrimer (Hz polarization, λ = 116 nm, a = 10 nm, g = 5 nm): cross sections, scattering indicatrix, and Poynting-vector streamlines are computed and subjected to a multi-test verification suite (truncation convergence, boundary residuals, optical theorem, energy-flux balance, full-vs-reduced symmetry).","tokens_in":28851,"tokens_out":1277,"duration_ms":30491,"significance":"If the accuracy claims hold, the work supplies a carefully engineered, high-precision multipole solver for a classical but still application-relevant class of problems (plasmonic oligomers, nanowire clusters, metamaterial lattices of circular rods). The analytic skeleton is standard, but the multistage adaptive linear-algebra strategy and the unusually thorough internal verification suite (Sec. VII E, Tables V–VI) are genuine strengths: residual control, optical-theorem and energy-balance checks at the 10^{-6}–10^{-7} level, and explicit handling of extreme ill-conditioning (κ up to ~10^{44}) make the scheme useful both for parametric studies and as a benchmark generator for FEM/FDTD/BEM codes. Scope limits (circular cross-sections, non-touching cylinders, moderate P) are stated honestly. The aluminum-trimer Poynting maps and screening analysis are of independent nano-optics interest.","major_comments":[{"comment":"The abstract and Sec. V claim controlled accuracy for densely packed subwavelength configurations (including g/a down to ~0.01 in Table V). The full physical verification suite of Sec. VII E (boundary-condition residual, three-route indicatrix, optical theorem ~10^{-6}, energy-flux balance, symmetry) is reported in detail only for the working case g/a = 0.5 (Table VI). For g/a ≲ 0.1 the paper mainly documents κ(A) and the solver level needed for ρ ≲ 10^{-10}, while noting that larger N may be required. Please add at least one denser-gap case (e.g. g/a = 0.1 or 0.01) with the same physical diagnostics (optical theorem, energy balance, boundary residual at the surfaces) so that the dense-packing claim rests on the same multistage evidence as the main example.","section":null},{"comment":"Sec. V A and the discussion after Table V correctly note that N must grow as gaps shrink and that κ grows rapidly with N. The forward-error control argument then relies on high-precision solution of the truncated system plus a posteriori physical tests. A short, explicit statement of the practical protocol—how N is increased until the physical observables (not only ρ) stabilize for a target g/a—would make the “controlled accuracy” claim fully operational for readers who wish to reproduce dense-pack runs.","section":null}],"minor_comments":[{"comment":"The comparison with FEM/FDTD in Sec. VIII is qualitative. A single quantitative cross-check (e.g. Q_ext or near-field |S|max for the same trimer against a commercial or open multipole/FEM code) would strengthen the positioning without changing the paper’s scope.","section":null},{"comment":"Footnotes on the sign of s_n relative to the Bohren–Huffman convention are easy to miss. A one-sentence remark in the main text near Eqs. (3.13) and (7.3) would reduce the risk of mis-comparison with the literature.","section":null},{"comment":"Figure 1: the color scale for |S| and the streamline density are informative, but a brief note in the caption on how many seed points survive thinning and whether streamlines are integrated through the interfaces would aid reproducibility of the topology discussion in Sec. VII B.","section":null},{"comment":"Eq. (5.4) for mp.dps is given as a heuristic; the trimer runs use a fixed mp.dps = 120. Stating which choice was used for each row of Table V would clarify the numerical protocol.","section":null},{"comment":"Minor typographical/consistency items: “kissing” cylinders citation [21] is fine but the academic-interest remark in Sec. II D could be shortened; ensure consistent use of k_0 vs k0 and of Gothic S_p for the scattering-coefficient matrix versus the Poynting vector S.","section":null},{"comment":"References [4–10] establish prior multipole/Graf work; a one-sentence contrast (what those solvers do not do regarding conditioning/precision) in the Introduction would sharpen the novelty claim without overstating it.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid methods contribution with unusually careful numerics; the main risk is over-claiming dense-pack accuracy relative to the verification actually shown. Scope is narrow (circular cylinders only), which the authors acknowledge—fit is good for an optics/computational-electromagnetics journal that values benchmark-quality solvers. No integrity or citation-pattern concerns. I would not block on the missing external code comparison if the denser-gap physical tests are added."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The multipole-plus-Graf skeleton is classical (Twersky, Linton–Martin, Schäfer, Mackowski, all cited). What is new is the four-level adaptive solver—equilibration + iterative refinement + mpmath + exact Gaussian-rational elimination—together with systematic condition-number monitoring and a multi-test verification suite. That package lets them keep controlled accuracy down to g/a ~ 0.01 where κ reaches 10^44, which is the practical contribution.\n\nThe analytic derivation (Secs. II–IV) is clean and standard. The numerical claims are unusually well backed: truncation convergence, boundary residuals ~10^{-5}, three independent routes to the indicatrix, optical-theorem error ~10^{-6}, energy-flux balance, and full-vs-reduced symmetry agreement at 10^{-13}. The aluminum-trimer example (λ=116 nm, a=10 nm, g=5 nm) is a sensible demonstration of near-field energy redistribution and screening, not a claim of record enhancement. Free parameters are only N, working precision, and residual thresholds; no fitted physics enters the central claims.\n\nSoft spots are real but proportional. Scope is limited to non-touching circular cylinders; the g=0 singularity is explicitly excluded. No code is released, so independent re-implementation would raise confidence further. Cost still scales as O((PM)^3), so the method is aimed at moderate clusters. None of these undercut the stated claims for the class of problems they treat.\n\nThis is for people who need benchmark near-field solutions for plasmonic nanowire oligomers or who want a high-accuracy reference against which to test FEM/FDTD. It does not open a new physical regime, but it is a reliable tool paper. I would send it to peer review; the verification work is thorough enough to deserve referee time. Worth citing if you need a controlled 2-D cylinder-cluster solver or a validation case.","headline":"Solid, carefully engineered multipole solver for circular-cylinder clusters; the real advance is the adaptive high-precision linear algebra that tames extreme ill-conditioning for dense subwavelength gaps.","tokens_in":29468,"tokens_out":501,"would_cite":true,"duration_ms":6479,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A multipole expansion plus high-precision linear algebra yields controlled-accuracy solutions for scattering by arbitrary clusters of parallel circular cylinders, including dense subwavelength packs.","keywords":["multiple scattering","Mie theory","cylindrical harmonics","Graf addition theorem","plasmonics","Poynting vector","optical theorem","nanowires"],"falsifier":"Compute the aluminum-trimer fields and optical theorem residual at N=18 and g/a=0.5 with an independent high-order method (or exact rational elimination) and check whether the boundary residual stays below ~10^{-5} and the optical-theorem discrepancy below ~10^{-6}; if either error exceeds those levels by an order of magnitude, the controlled-accuracy claim fails.","tokens_in":29583,"feed_emoji":"⚡","tokens_out":732,"duration_ms":8461,"temperature":0.7,"pith_summary":"The paper sets out a semi-analytical method for two-dimensional electromagnetic scattering by any number of parallel, non-overlapping circular cylinders of arbitrary radii and complex permittivities. Fields are expanded in cylindrical harmonics about each cylinder axis; Graf's addition theorem converts multiple scattering into a linear system for the modal coefficients. That system is solved with condition-number monitoring and, when needed, extended-precision or exact arithmetic, then checked by a multistage convergence suite. The claim is that the procedure delivers numerically verified, controlled-accuracy solutions across a wide parameter range, including densely packed subwavelength geometries where ordinary floating-point solvers fail. A worked example—an aluminum nanotrimer under normal-incidence H_z polarization—computes cross-sections, the scattering pattern, and Poynting-vector streamlines that show energy redistribution and localized surface enhancement. A sympathetic reader cares because the method supplies exact boundary satisfaction and analytic near- and far-field evaluation for a class of geometries that appear repeatedly in plasmonics and metamaterial design, while remaining cheap enough for parametric studies.","feed_headline":"High-precision multipoles tame dense cylinder clusters","feed_subtitle":"Controlled-accuracy scattering solutions hold even for subwavelength aluminum nanowire trimers","key_machinery":"The block system (5.1)–(5.2) assembled from single-cylinder scattering matrices S_p and Graf translation matrices T_pq, solved by the four-level adaptive procedure that escalates precision according to the condition number κ(A).","core_discovery":"After truncation of the multipole series, the multiple-scattering problem for an arbitrary ensemble of parallel circular cylinders reduces to a block linear system whose diagonal blocks are single-cylinder scattering coefficients and whose off-diagonal blocks are products of those coefficients with Graf translation matrices; a four-level adaptive solver that monitors the condition number and escalates from ordinary LU through equilibration, arbitrary-precision arithmetic, and exact elimination over the Gaussian rationals produces residuals small enough that physical diagnostics (boundary conditions, optical theorem, energy balance) remain at the 10^{-5}–10^{-6} level even for subwavelength g","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Adaptive multipoles solve dense cylinder EM scatter to 1e-6","Graf translations yield verified residuals for nanowire packs","Block multipole systems tame arbitrary parallel-cylinder clusters","Precision harmonics map flux redistribution in Al nanowire trimers","Escalating solvers hold optical theorem for subwavelength cylinders"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"That a finite multipole cutoff together with high-precision solution of the truncated system is enough to keep the forward error under control for arbitrarily small but non-zero gaps between cylinders.","fun_headline_variants_meta":{"raw":{"variants":["Adaptive multipoles solve dense cylinder EM scatter to 1e-6","Graf translations yield verified residuals for nanowire packs","Block multipole systems tame arbitrary parallel-cylinder clusters","Precision harmonics map flux redistribution in Al nanowire trimers","Escalating solvers hold optical theorem for subwavelength cylinders"]},"model":"grok-4.5","effort":"low","cost_usd":0.005102,"raw_usage":{"total_tokens":1437,"prompt_tokens":832,"num_sources_used":0,"completion_tokens":83,"cost_in_usd_ticks":51020000,"prompt_tokens_details":{"text_tokens":832,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":522,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":832,"tokens_out":83,"duration_ms":7121,"temperature":1.0,"reasoning_tokens":522,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T00:13:03.943777+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compute the aluminum-trimer fields and optical theorem residual at N=18 and g/a=0.5 with an independent high-order method (or exact rational elimination) and check whether the boundary residual stays below ~10^{-5} and the optical-theorem discrepancy below ~10^{-6}; if either error exceeds those levels by an order of magnitude, the controlled-accuracy claim fails.","supporting_citations":[],"review_version":2}