{"id":"a106c7b1-0ca5-452a-865d-e812d8ee4522","arxiv_id":"2412.13826","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A rational-approximation framework with a centrally placed line source computes selected waveguide resonance modes without computing the full mode spectrum.","lead":"This paper computes selected guided modes of a hollow-core photonic crystal fiber by fitting a rational function to simulated scattered fields from a carefully placed line source. The approach avoids computing the many cladding and higher-order modes that standard eigensolvers produce, removing a costly mode-filtering step.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The framework's pole-selection step equates 'significant pole influence' with 'strong source coupling,' but a high-Q mode with weak coupling can dominate the rational approximation, so the method does not guarantee computing only relevant modes without a formal selection rule.","rationale":"The reader's weakest assumption targets the source-coupling premise: a significant pole implies a significantly coupled mode, and the chosen line source isolates the fundamental mode. My concern sharpens this: pole influence is not equivalent to coupling because the influence also depends on the pole's distance to the sampling contour. This makes the selection step load-bearing in a way the reader did not fully separate out. The paper's own conclusion already limits the framework to systems where the wanted mode is 'much more relevant' than others, but no criterion is given to test that condition or to pick peaks automatically. The accuracy and convergence results for the single HC-PCF example are strong and machine-checkable through the provided data, so I do not dispute the numerical demonstration; however, the generality of the 'efficient and accurate computation of the fundamental resonance mode' claim is conditional on formalizing pole selection and demonstrating that source coupling, not spectral proximity, drives the outcome. The efficiency claim is also unquantified, but the pole-selection issue is more fundamental because it affects correctness, not just speed. Therefore the existing CONDITIONAL verdict is appropriate, and my analysis does not move it.","tokens_in":7060,"tokens_out":8400,"duration_ms":87716,"concrete_test":"Use the published dataset [21] to list all poles from the AAA approximation with 40 sampling points. For each pole, compute the vector residue a_n from Eq. (4), the imaginary part |Im(neff_n)|, and the true source-coupling coefficient c_n = |∫ J·E_n dV| from the Arnoldi reference modes. If any pole with |a_n| above the visually selected threshold has c_n orders of magnitude smaller than the fundamental mode's c_n, then pole selection is controlled by Q rather than source coupling. Alternatively, rerun the framework with the line source displaced into a cladding air hole; if neff_1 still emerges as a dominant peak, the method's selectivity is not due to source placement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section II, the paper states: 'When a pole z_pole has a significant influence on the rational approximation r(z), then we assume that a_n and z_pole are a good approximation to an eigenpair... This means that the resonance mode a_n ... has a significant coupling with the source term s(z).' This equates pole influence on r(z) over the real sampling interval [0.995,1] with physical source coupling. But the contribution of a pole to r(z) scales with its residue divided by the distance to the sampling contour, not with coupling alone. A resonance with an extremely small imaginary part will produce a sharp, dominant peak even with weak excitation, while a low-Q mode with strong coupling may appear only as broad background. In the HC-PCF example, Im(neff1)=3.376e-9, so the fundamental mode naturally dominates; the example therefore cannot validate the general claim that the source 'has a significant coupling with the mode' while other modes have negligible coupling. The paper selects 'two significant peaks' by inspection, without an objective threshold. Thus the central claim that the framework computes only relevant modes is not established for arbitrary microstructured waveguides; it depends on spectral proximity and on an unformalized peak-picking step.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a framework for computing selected resonance modes of microstructured photonic waveguides by combining a specially placed line source with AAA rational approximation of the resulting scalar scattering data. For the hollow-core photonic crystal fiber example, the AAA-derived eigenvalues converge to Arnoldi reference values with real-part relative errors below 10^-14, imaginary-part errors below 10^-5 (and 10^-8 for the second mode), and mode-field errors below 10^-5. The authors argue that the source coupling avoids the computation of cladding and higher-order modes and thus eliminates post-processing mode filtering.","tokens_in":7325,"tokens_out":3223,"duration_ms":34701,"significance":"If the central claim is sustained, the approach would be practically useful because it targets relevant modes directly while using only black-box scattering solves, and its reproducibility assets are strong: source code and simulation data are deposited in an open data publication, and the numerical convergence studies are clearly reported. The main limitation is that the example itself cannot establish the general coupling assumption, since the demonstrated fundamental mode has an extremely small imaginary part and therefore dominates the rational approximation for spectral reasons, not necessarily because of superior source coupling.","major_comments":[{"comment":"The term \"efficient\" is load-bearing in the title and abstract, but the manuscript reports no runtime, no number of scattering solves, and no comparison with the Arnoldi computation. Reporting 40 sampling points is not sufficient; a wall-clock time or solver-cost comparison with the Arnoldi reference computation would be needed to support the efficiency claim.","section":"Abstract and Section I"},{"comment":"The assumption that a pole with significant influence on r(z) implies an eigenpair with significant source coupling is explicitly stated but never validated. The HC-PCF example has Im(neff1)=3.376e-9, so the fundamental mode's pole dominates the rational approximation because of its very small imaginary part, not necessarily because of exceptional source coupling. To support the general claim, the authors should either prove or numerically test the selection rule for cases with weak coupling and high quality factor, or with strong coupling and low quality factor.","section":"Section II, paragraph after Eq. (4)"},{"comment":"The selection of \"two significant peaks\" is made by visual inspection, and no objective criterion is given for distinguishing relevant poles from spurious poles, background continuum, or poles outside the sampled interval. Since the central claim is that only relevant modes are computed, a threshold on residue magnitude, pole proximity, or another well-defined measure is required, together with a demonstration that the selected set is stable with respect to that threshold.","section":"Section III B, Fig. 3"},{"comment":"The reference Arnoldi solutions are obtained with the same FEM solver, material model, and computational domain as the scattering data used in the AAA approximation. Therefore the reported convergence demonstrates consistency between two methods sharing a common discretization, not absolute accuracy. An independent check, such as FEM mesh refinement or comparison with an analytic waveguide benchmark, is needed to substantiate the claim of accurate computation.","section":"Section III A and III B, Figs. 3 and 4"}],"minor_comments":[{"comment":"The caption says the mode E2 corresponds to neff1; this should read neff2.","section":"Section III C, Fig. 4 caption"},{"comment":"The residue formula as typeset appears to list two factors separated by a comma in the denominator; the intended expression is likely a single quotient, and the notation should be corrected for clarity.","section":"Section II, Eq. (4)"},{"comment":"The random projection vector y is drawn from a uniform distribution, but the manuscript does not report the random seed or the variation of the results over different draws of y; this information should be included so that the numerical experiment is reproducible.","section":"Section III B"},{"comment":"The term \"special light sources\" is used informally in the introduction and conclusion; since the actual source is a singular line source, the terminology should be defined consistently at first use.","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a short application note built on the authors' previous work in Ref. [19], and its reproducibility practices are good. The numerical demonstration is convincing as a consistency test, but the central generality claim depends on an unproven source-coupling assumption and an unformalized peak-selection step. I do not see integrity problems, but the efficiency claim needs quantitative support and the selection rule needs an objective formulation before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nShort version: this is a credible incremental numerical paper, and I would send it out. It takes the AAA-with-special-source idea from the authors' own Ref [19] and demonstrates it on a hollow-core photonic crystal fiber, getting neff1 = 0.9993596784939 + 3.376e-9 i with real-part errors below 1e-14 and mode errors below 1e-5 against Arnoldi. That is a concrete, reproducible result, and the Zenodo deposit backs it up.\n\nWhat the paper does well: the accuracy validation is honest. The AAA eigenvalues are outputs of a rational fit to scattering data, not fitted to the Arnoldi values; the reference is used only for comparison. The convergence study over number of sampling points is a real check. The source-placement argument is intuitive and works in the example. The authors are also candid that the core computation uses the approach of Ref [19]; self-citation here is legitimate because that is where the algorithm is defined.\n\nWhere it is soft: the headline word \"efficient\" has no numbers behind it. We never learn runtime or memory for AAA versus computing 512 Arnoldi modes, so the practical selling point is unquantified. The pole-selection step is inspection: \"two significant peaks\" is not a defined criterion, and details on the other eigenvalues are pushed to the data publication. That matters because the framework's generality rests on the assumption that significant poles correspond to strongly coupled modes. The stress-test concern about high-Q weak-coupling modes is correct in principle: the example has the fundamental mode with Im(neff) = 3.4e-9, so it naturally dominates the rational approximation. That does not invalidate the demonstrated result, but it does mean the general claim that only relevant modes are computed should be stated as a heuristic that holds when the source is well chosen. Using the same FEM solver for both scattering data and Arnoldi reference is a mild shared-discretization circularity, not a fatal one.\n\nWho should read it: people in computational photonics who want targeted resonance-mode calculations and are deciding between eigensolvers and scattering-based rational approximation. They will find a useful data point, but should not expect a general-purpose algorithm with formal guarantees.\n\nRecommendation: send to peer review. Ask the authors to benchmark against a standard eigensolver, to specify an automatic peak-selection rule, and ideally to add a second waveguide where the spectrum is denser or the desired mode is not the lowest-loss one. The paper is solid enough to warrant that revision.","headline":"A credible, well-documented numerical demonstration of AAA-based resonance-mode computation for a hollow-core fiber; the method is incremental, and the two real soft spots are the unquantified efficiency claim and the informal pole-selection step.","tokens_in":7842,"tokens_out":2601,"would_cite":true,"duration_ms":25805,"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":"This paper shows that applying AAA rational approximation to the scattered field of a single centrally placed line source recovers the fundamental resonance mode of a microstructured waveguide while skipping the cladding and higher-order…","keywords":["AAA algorithm","rational approximation","resonance modes","photonic waveguide","hollow-core photonic crystal fiber","eigenvalue computation","source coupling","Maxwell equations"],"falsifier":"Run the same AAA workflow with a line source moved into the cladding or shaped like a higher-order mode profile, and compare the returned poles against the full reference spectrum: a significant pole that matches a cladding mode but misses the fundamental mode would confirm the output is source-selected, while a significant pole with no nearby reference eigenvalue would refute the claim that significant poles are eigenpairs.","tokens_in":6832,"feed_emoji":"💡","tokens_out":8165,"duration_ms":67319,"temperature":0.7,"pith_summary":"This paper claims that the fundamental resonance mode of a microstructured photonic waveguide can be computed directly, without computing the cladding and higher-order modes that crowd the spectrum. The method applies AAA rational approximation to the scattered field produced by a line source placed at the waveguide center, so that the poles of the rational approximation correspond to modes the source actually excites. On a hollow-core photonic crystal fiber from the literature, the approach returns the fundamental mode's effective index with real-part errors below $10^{-14}$ and imaginary-part errors below $10^{-5}$ relative to reference eigensolver computations, and mode-field errors below $10^{-5}$. If this holds, it removes the need to solve a large eigenproblem and filter irrelevant modes, and it makes sensitivity information available at negligible extra cost.","feed_headline":"One line source isolates the fundamental waveguide mode","feed_subtitle":"Rational AAA fit of the scattered field skips cladding and higher-order modes entirely.","key_machinery":"The load-bearing object is the AAA (adaptive Antoulas-Anderson) rational approximation in barycentric form, $r(z) = n(z)/d(z)$, whose poles $z_{\\mathrm{pole},n}$ and residues $a_n$ are read off directly. The paper reuses the scalar weights and poles in a vector-valued residue formula to assemble the mode field from the finite-element coefficient vectors of the scattered field. The second mechanism is source selection: a line source at the center of the hollow core couples strongly to the fundamental mode and negligibly to cladding and higher-order modes, so the rational approximation is dominated by the wanted eigenpair.","core_discovery":"The central claim is that one scalar projection of scattering data is enough to isolate a wanted resonance mode, provided the illuminating source couples strongly to that mode. Concretely, the paper solves the time-harmonic Maxwell scattering problem for an $x$-polarized line source at the center of the fiber at 40 sampling values of the effective index, forms the scalar function $y^T E_x$ with a random vector $y$, and fits it with AAA rational approximation. The dominant poles of that fit are interpreted as eigenvalues, and the vector-valued residue formula reconstructs the mode field. For the hollow-core photonic crystal fiber, this yields $n_{\\mathrm{eff}1} = 0.9993596784939 + 0.000000003376i$ and $n_{\\mathrm{eff}2} = 0.996754264645 + 0.00000190093i$, matching reference eigenvalues to $10^{-14}$ in the real part and to $10^{-5}$--$10^{-8}$ in the imaginary part. The identifying property of the fundamental mode is its central-core localization, which gives the on-axis source a dominant coupling to it.","pith_inferences":["The same source-coupling principle should let one target any specific mode, not just the fundamental one, by engineering an incident field whose spatial profile overlaps that mode and is mostly orthogonal to the others.","A direct test of the method's generality is to replace the line source with a Gaussian beam or a measured higher-order-mode profile and check whether the returned pole tracks the intended mode's eigenvalue.","In inverse design, the low-cost sensitivities could make this procedure a practical objective-function evaluator, optimizing the waveguide geometry to place a selected mode's eigenvalue at a target value without re-solving the full spectrum."],"forward_implications":["Only the relevant modes are computed: in the demonstration, two core-localized modes emerge from 40 scattering solves instead of the 512 eigenvalues the reference eigensolver produces.","No mode-filtering post-processing is needed, because undesired cladding and higher-order modes never enter the rational approximation.","Because the approach relies on scattering solves, sensitivities of the eigenvalues with respect to geometry or material parameters come at negligible extra cost via algorithmic differentiation.","Accuracy can be pushed further by choosing complex sampling points near the physical eigenvalues, as the paper notes.","Other source types, such as multiple line sources, a fundamental-mode field of a single-mode fiber, or a Gaussian beam, can be used in the same framework."],"supporting_citations":[{"why":"Supplies the AAA rational approximation algorithm that produces poles and residues from the sampled scattering data.","marker":"[14]"},{"why":"Establishes the source-coupling assumption linking significant poles to eigenpairs and provides the vector-valued residue formula used to reconstruct modes.","marker":"[19]"},{"why":"Introduces the hollow-core photonic crystal fiber geometry used as the demonstration system.","marker":"[11]"},{"why":"Provides the random-vector projection used to reduce the vector-valued field to the scalar function fed into AAA.","marker":"[15]"},{"why":"Supplies the iterative eigenvalue algorithm used to produce reference eigenvalues and modes.","marker":"[22]"},{"why":"Supplies the resonance-mode computation context and reference methods that motivate the comparison.","marker":"[6]"},{"why":"Supplies a finite-element eigensolver implementation used as part of the reference solution route.","marker":"[7]"}],"fun_headline_variants":["Scalar fit picks out a single waveguide mode","AAA rational fit targets one resonance mode","One projection isolates the core mode","Line source plus AAA rational fit selects mode","Compute a single mode via scalar AAA fit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes that whenever a pole has a significant influence on the rational approximation, its vector residue is an accurate eigenvector, because the source couples strongly to that mode; modes the source does not excite are assumed to stay out of the approximation.","fun_headline_variants_meta":{"raw":{"variants":["Scalar fit picks out a single waveguide mode","AAA rational fit targets one resonance mode","One projection isolates the core mode","Line source plus AAA rational fit selects mode","Compute a single mode via scalar AAA fit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000301,"raw_usage":{"total_tokens":1698,"prompt_tokens":867,"completion_tokens":831,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":483,"completion_tokens_details":{"reasoning_tokens":767}},"tokens_in":483,"tokens_out":831,"duration_ms":6378,"temperature":1.0,"reasoning_tokens":767,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:44:59.956415+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same AAA workflow with a line source moved into the cladding or shaped like a higher-order mode profile, and compare the returned poles against the full reference spectrum: a significant pole that matches a cladding mode but misses the fundamental mode would confirm the output is source-selected, while a significant pole with no nearby reference eigenvalue would refute the claim that significant poles are eigenpairs.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the AAA rational approximation algorithm that produces poles and residues from the sampled scattering data."},{"cited_title":"Evaluation of resonances: adaptivity and AAA rational approximation of randomly scalarized boundary integral resolvents","cited_arxiv_id":"2405.19582","evidence_quote":"Establishes the source-coupling assumption linking significant poles to eigenpairs and provides the vector-valued residue formula used to reconstruct modes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the hollow-core photonic crystal fiber geometry used as the demonstration system."},{"cited_title":"Nakatsukasa, O","cited_arxiv_id":null,"evidence_quote":"Provides the random-vector projection used to reduce the vector-valued field to the scalar function fed into AAA."},{"cited_title":"Binkowski, F","cited_arxiv_id":null,"evidence_quote":"Supplies the iterative eigenvalue algorithm used to produce reference eigenvalues and modes."},{"cited_title":"Lalanne, W","cited_arxiv_id":null,"evidence_quote":"Supplies the resonance-mode computation context and reference methods that motivate the comparison."},{"cited_title":"Lalanne, W","cited_arxiv_id":null,"evidence_quote":"Supplies a finite-element eigensolver implementation used as part of the reference solution route."}],"review_version":1}