{"id":"c3340959-24df-4df3-b36a-35d48bc01cb6","arxiv_id":"1908.06916","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A new resonant-state expansion method calculates the complete set of optical resonances of a periodic photonic-crystal slab from a homogeneous-waveguide basis.","lead":"Photonic-crystal slabs are usually simulated with scattering-matrix or finite-difference codes; this paper instead expands their optical resonances in a complete basis of homogeneous-slab modes, turning Maxwell's equations into a matrix diagonalization. The method promises a complete set of resonances with no missing or spurious modes, and is checked against the scattering matrix method for a periodically modulated slab.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix F shows PC-RSE fails at b=a, the full-thickness-modulation case most relevant to real slabs; the claimed completeness and generality of the method is not established over its stated parameter range.","rationale":"The reader's weakest assumption correctly identifies the Mittag-Leffler basis completeness as the load-bearing element and notes the admitted b=a failure. My stress-test pass confirms that this is the single most consequential limitation of the paper's central claim. The claim of asymptotic exactness and guaranteed completeness is stated without the b<a qualifier, yet the only fully independent verification is performed away from the boundary (b=a/2), and the one exploratory step toward the boundary (b=0.95a) already shows degraded convergence. The b=a case is not an exotic corner: many practical PC slabs are modulated through their full thickness. The paper deserves credit for stating the limitation explicitly in Appendix F and for providing a nontrivial analytic benchmark in Appendix D, which is genuine independent support for the method in the homogeneous-perturbation limit. However, that support does not validate the boundary-touching regime. The M=5 SMM baseline is a secondary weakness, but the b=a failure is more load-bearing because it attacks the generality of the completeness claim itself. Since the paper already received a CONDITIONAL verdict and my concern is the same one the reader identified, no change in verdict is needed: the paper should remain conditional until the b=a regime is either fixed, convincingly tested, or explicitly excluded from the scope of the claims.","tokens_in":26282,"tokens_out":7047,"duration_ms":80208,"concrete_test":"Recompute the PC-RSE for the same slab parameters as Fig. 3 but with b=a (full-thickness modulation) and compare against an independent SMM calculation at M=7 and M=9, varying the cut-mode density F=0.5, 1, and 2. If the relative errors in the RS frequencies are below ~1e-3 and decrease with basis size N, the b=a restriction can be lifted; if the errors remain elevated or depend strongly on F, the central claim must be explicitly restricted to b<a.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Sec. IV) is that PC-RSE is asymptotically exact and returns a complete set of resonant states controlled by a single truncation parameter, with no missing or spurious modes, for planar photonic-crystal structures. That claim rests on the Mittag-Leffler expansion of the homogeneous-slab Green's function remaining a convergent basis once the periodic perturbation is added, leading to the eigenvalue problem in Eq. (23). Appendix F explicitly reports the opposite for the case b=a: \"Presently, it prevents the PC-RSE from being used with exactly b=a,\" and even at b=0.95a the relative errors are up to an order of magnitude larger than at b=a/2 (Fig. 16). Since Eq. (32) defines the modulation for any b<=a, and full-thickness etched gratings are a standard photonic-crystal-slab geometry, this is not a remote corner of parameter space. The demonstrated cases (b=a/2, d=2pi/5) therefore do not support the unqualified claim of a complete set of RSs for a photonic-crystal slab; the basis-completeness assumption underlying Eq. (23) is unverified precisely when the perturbation reaches the slab boundary, and the paper itself leaves this as requiring further study.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a photonic-crystal resonant-state expansion (PC-RSE) for planar photonic-crystal slabs. The method uses as a basis the resonant states and discretized branch-cut modes of a homogeneous planar slab, taken over all Bragg channels, and maps Maxwell's equations for a periodic permittivity perturbation onto a linear matrix eigenvalue problem, Eq. (23)/(28). The central results are the derivation of this eigenvalue problem, numerical verification against an analytic core-shell solution for homogeneous perturbations (Appendix D) and against the scattering-matrix method for a one-dimensionally periodic dielectric slab (Section III C), and an application tracing the formation of bound states in the continuum and quasi-guided modes as the periodic modulation amplitude increases.","tokens_in":26453,"tokens_out":7524,"duration_ms":76255,"significance":"If the stated claims hold, the PC-RSE would be a significant methodological contribution: it replaces a nonlinear eigenproblem for open photonic-crystal resonances by a linear matrix diagonalization, uses a basis that is constructed analytically from the homogeneous slab, involves no fitted parameters, and exhibits fast convergence (roughly 1/N^3) in the tested cases. The formal derivation from Maxwell's equations to Eqs. (21)-(23) is clean and is a genuine strength of the paper, as is the use of an analytic core-shell solution for independent verification. The numerical demonstrations of BIC formation and basis-mode contributions are also informative. The central claim of completeness and asymptotic exactness is, however, currently overstated relative to what is demonstrated, because the method is shown in Appendix F to fail in an important parameter regime (b = a).","major_comments":[{"comment":"The unqualified conclusion in Sec. IV that the PC-RSE is 'asymptotically exact' and 'guarantees completeness, i.e. has no missing or spurious modes' is not supported over the parameter range defined by Eq. (32). Appendix F explicitly states that the method cannot be used with exactly b = a, and that at b = 0.95a the relative errors are up to an order of magnitude larger than at b = a/2 (Fig. 16). Since b = a (modulation reaching the slab boundaries) is a standard geometry for etched photonic-crystal slabs and is included in the model via Eq. (32), this is not a peripheral corner of parameter space. The basis-completeness assumption underlying Eq. (23) is therefore unverified precisely when the perturbation extends to the slab boundary, and the paper itself notes that 'the ML series changing its convergence properties on the borders of the system ... requires a further study.' The authors should either extend the method to cover b = a, or explicitly qualify the abstract and conclusions so that the completeness claim is restricted to the demonstrated regime b < a. A concrete test for the b = a case would be a comparison with SMM at larger M, together with a study of convergence of the PC-RSE eigenvalues as N increases.","section":"Appendix F; Sec. IV; Eq. (32)"},{"comment":"The quantitative verification of the PC-RSE against the SMM uses M = 5 Bragg channels as the 'exact' reference. The SMM at finite M is itself an approximation, and the observed discrepancies near the cuts (e.g., region 3 in Fig. 3) may be due to the limited SMM accuracy rather than to PC-RSE error. To make the accuracy claim quantitative, the paper should show convergence of the SMM reference with increasing M, report errors relative to a converged SMM result, or explicitly state that the comparison is at fixed, low M and that the quoted errors are upper bounds. Without this, the statement that the PC-RSE is 'unprecedentedly accurate' is not fully quantified.","section":"Section III C; Fig. 4"}],"minor_comments":[{"comment":"Equation (C8) appears to contain a sign typo: the expression B_n(e^{-iq_n z} + (-1)^n e^{-iq_n z}) has two identical exponential factors; presumably one of them should be e^{+iq_n z}, consistent with Eq. (29).","section":"Eq. (C8)"},{"comment":"The conclusion that the method 'depends on a single parameter, the truncation frequency omega_max' is not quite accurate, because the basis also depends on the cut-mode fraction F, which is set to F = 1 based on a numerical optimization study in Appendix D2. The role of F should be acknowledged in the conclusions.","section":"Sec. IV; Appendix D2"},{"comment":"The phrase 'unprecedented accuracy' in the abstract is a strong claim that is not directly compared with other numerical methods in the paper; it would be safer to describe the accuracy as 'high' and quantify it in the specific examples.","section":"Abstract"},{"comment":"In Fig. 4(b), the text states that N_tot ≈ 12000 is used as the reference, but the figure labels are not entirely clear; please ensure that the reported reference basis size is consistent between text, caption, and legend.","section":"Fig. 4"},{"comment":"The notation in Eq. (23) uses an integral symbol to denote summation over the continuum of cut modes, but this is not defined explicitly at that point. Please add a sentence clarifying that the integral represents the discretized/continuum cut-mode contribution, as is done later in Appendix D2.","section":"Sec. II; Eq. (23)"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a technically sound method in the tested regime and is likely publishable after major revision. The principal barrier is the mismatch between the broad claims in the abstract and conclusions and the explicit limitation reported in Appendix F. I see no citation or novelty concerns beyond the authors' own references to prior RSE work. If the authors re-scope the claims or provide a satisfactory treatment of b = a, the paper would be acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real content here is better than the packaging. Neale and Muljarov have done something not present in the cited RSE literature: a working resonant-state expansion for a periodically modulated slab that includes multiple Bragg channels and discretized cut modes, and they test it honestly. The derivation from Maxwell's equations to the linear eigenvalue problem (23) is clean, and the convergence tests against the analytic core-shell solution are real evidence. The observation that all diagonal perturbation matrix elements vanish for a purely periodic modulation, and that this accelerates convergence, is a nice non-obvious point. The mode-decomposition analysis of BICs and quasi-guided modes is physically useful.\n\nThe soft spots are in proportion. The main one is not hidden: Appendix F admits the PC-RSE cannot currently be used at b = a, and by b = 0.95a the errors are up to an order of magnitude larger. Full-thickness modulation is the standard etched-grating geometry, so this is not a corner case. That means the conclusion's unqualified statement that the approach is asymptotically exact and yields a complete set of RSs for a photonic-crystal slab overstates what has been shown. The claim should be qualified to the demonstrated b < a range, or the boundary problem should be fixed. This is a real limitation, not a fatal one: the paper already shows the method works well for b = a/2 and a range of periods.\n\nSecond, the SMM comparison uses only M = 5 Bragg channels, so the \"exact\" reference is itself approximate. That weakens the quantitative error claims, but the convergence-to-self behavior in Fig. 4(b) and the analytic core-shell verification in Appendix D partly compensate. An independent reference with more channels would settle it.\n\nMinor: no code is provided, so replication requires rebuilding the basis from scratch. And the abstract's \"new paradigm\" and \"unprecedented accuracy\" are more than the demonstrated scope supports.\n\nCitation pattern looks fine. The distinction from Refs. 19 and 20 is correct.\n\nI would send this to peer review. It is a legitimate method paper with a real derivation and meaningful numerical verification, and the main caveat is disclosed in the paper itself. The right referee will ask for the b = a limitation to be either fixed or prominently qualified, and for a stronger external benchmark. Worth reading for anyone doing QNM/RSE calculations on patterned slabs or BIC engineering.","headline":"A genuine resonant-state expansion for periodic slabs with clean derivation and honest convergence tests, but the paper's own Appendix F shows it fails at full-thickness modulation, so the completeness claim is overbroad as stated.","tokens_in":27042,"tokens_out":2832,"would_cite":true,"duration_ms":32184,"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":"Expanding a photonic-crystal slab's resonances in the exact resonant states of a homogeneous slab yields a complete, asymptotically exact set with no missing or spurious modes.","keywords":["resonant-state expansion","photonic crystal slab","bound states in the continuum","quasi-guided modes","Mittag-Leffler expansion","cut modes","Bragg channels","scattering-matrix method"],"falsifier":"For a slab with $\\epsilon = 6$, $b = 0.95a$, $d = 2\\pi/5$, $\\beta = 1$, and $p = 0$, compute the resonant frequencies with a scattering-matrix code using more Bragg channels than $M = 5$; if the PC-RSE results near $\\mathrm{Re}(\\omega a) = 5$ do not keep moving toward those values as $N$ grows, or if the scattering-matrix calculation finds modes the PC-RSE's fixed-$N$ diagonalization misses, the claim of guaranteed completeness fails in practice.","tokens_in":26008,"feed_emoji":"💡","tokens_out":6082,"duration_ms":62931,"temperature":0.7,"pith_summary":"Photonic-crystal slabs are open optical systems: light can leak out, so their natural eigenmodes are resonant states with complex frequencies, and a computation that misses one mode or adds a fake one corrupts the spectrum. This paper claims that all such states of a periodically modulated dielectric slab can be obtained by expanding them in the analytically known resonant states of an unmodulated slab, provided the branch-cut continuum is added as discretized cut modes. The resulting photonic-crystal resonant-state expansion (PC-RSE) is a linear matrix eigenvalue problem with one convergence parameter, the truncation frequency, and the paper demonstrates relative errors falling roughly as $1/N^3$. If the claim holds, the approach turns a previously delicate mode-searching task into a routine diagonalization and supplies the complete mode set needed to interpret spectra.","feed_headline":"Plain slab modes yield complete photonic-crystal resonances","feed_subtitle":"A single frequency cutoff controls accuracy, with errors shrinking by an order of magnitude per basis doubling.","key_machinery":"The load-bearing object is the Mittag-Leffler expansion of the homogeneous-slab dyadic Green's function, written as a sum over resonant states and cut modes. From that expansion the wave function of any photonic-crystal slab state is built from basis functions $F_n(z;p+g)e^{i(p+g)x}$, and the perturbation matrix $V^{gg'}_{nn'}$ is formed from the Fourier coefficients of the periodic permittivity change. Substitution turns Maxwell's equations into the linear eigenvalue problem $\\omega\\sum(\\delta+V)c = \\omega^g_n c^g_n$. The same completeness that makes the Green's function series converge is what guarantees that the diagonalized matrix returns every perturbed state and no spurious ones; the vanishing diagonal elements for a purely periodic modulation with zero mean make first-order perturbation vanish and accelerate convergence.","core_discovery":"The central claim is that the resonant-state expansion, previously developed for finite resonators and inhomogeneous waveguides, can be extended to infinite planar photonic-crystal slabs by using the homogeneous slab as the basis and treating the periodic permittivity modulation as a perturbation. Because the slab is periodic in one direction, the basis must include resonant states from many Bragg channels, labelled by reciprocal-lattice vectors $g$, and the Green's function of the homogeneous slab acquires branch cuts in the complex-frequency plane; the paper handles these cuts by discretizing them into artificial cut modes included on the same footing as ordinary resonant states. The result is the linear eigenvalue problem (28), whose diagonalization yields a complete set of perturbed resonant states with no missing and no spurious modes, controlled by a single truncation frequency $\\omega_{\\mathrm{max}}$. Numerical verification against the scattering-matrix method for a dielectric slab with a harmonic periodic modulation shows agreement that improves as roughly $N^{-3}$ with basis size, and mode-expansion analysis shows that bound states in the continuum emerge from pairs of degenerate waveguide modes, while the companion quasi-guided modes differ by their coupling to leaky modes of the zeroth Bragg channel.","pith_inferences":["If the completeness guarantee survives in two-dimensionally periodic structures, the method could become a workhorse for inverse design of metasurfaces, where current methods must painstakingly verify mode sets.","The exact $b=a$ failure boundary suggests a natural independent test: profiles with the modulation reaching the slab surface should be checked with other methods to determine whether the slowdown is true incompleteness or only slower convergence fixable by enlarging the cut-mode set.","Because vanishing diagonal elements speed convergence, designs that keep the periodic modulation at zero mean should converge faster than those with a uniform component; this is an implicit design rule one could test in cavity optimization.","The BIC-QGM pair analysis implies that coupling to the zeroth Bragg channel's leaky modes is what turns a would-be BIC into a high-Q quasi-guided mode; engineering that coupling by symmetry breaking or by moving the modulation toward the boundary gives a tunable Q-factor knob."],"forward_implications":["The complete, spurious-free mode set makes transmission, reflection, scattering, and extinction computable as superpositions of resonant states, without background fit parameters.","Only one parameter ($\\omega_{\\mathrm{max}}$) controls accuracy, so parameter sweeps over modulation strength, period, and layer thickness become automated matrix diagonalizations, useful for optimizing photonic-crystal cavities.","Bound states in the continuum are shown to be symmetry-protected states formed from waveguide modes, with quasi-guided partners distinguished by the presence of zeroth-channel leaky-mode components; calculations can now target BICs deliberately.","Cut modes, representing Rayleigh-Wood anomalies, are not optional decoration: including about one cut mode per resonant state ($F\\approx 1$) restores $1/N^3$ convergence for modes near the cuts.","The same formalism carries over to TM polarization and, the paper argues, to oblique incidence and two-dimensional periodicity without changing its structure."],"supporting_citations":[{"why":"Introduced the resonant-state expansion mapping Maxwell's equations to a linear eigenvalue problem, which this paper extends to periodic slabs.","marker":"[26]"},{"why":"Showed how branch cuts of the Green's function are included via discretized cut modes in inhomogeneous waveguides.","marker":"[14]"},{"why":"Provided the RSE for a homogeneous planar slab in the k-representation and the normalization used for basis states.","marker":"[17]"},{"why":"Established the convergence scaling and matrix formulation for open 2D systems used here.","marker":"[27]"},{"why":"Treated cut discretization for 2D open systems, adapted here for the slab cuts.","marker":"[31]"},{"why":"Gave the generalized normalization and operator form for magnetic and bi-anisotropic systems used in Eq. (20).","marker":"[34]"},{"why":"Is the scattering-matrix method used as the numerical reference for verification.","marker":"[7]"},{"why":"Represents the incomplete guided-mode expansion that motivates adding leaky modes and cut modes to the basis.","marker":"[36]"},{"why":"Is the Mittag-Leffler theorem that justifies the Green's function expansion.","marker":"[51]"}],"fun_headline_variants":["Photonic-crystal resonances from plain slab modes","Waveguide basis nails photonic-crystal resonances","Slab modes unlock complete photonic-crystal resonances","New paradigm: RSE for photonic crystals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire method rests on the premise that the resonant states of the unmodulated slab, supplemented by the discretized cut modes, form a complete basis for the modulated slab; the paper reports that convergence changes and deteriorates when the modulation layer reaches the slab boundary ($b = a$), so exactly that limiting case is not covered.","fun_headline_variants_meta":{"raw":{"variants":["Photonic-crystal resonances from plain slab modes","Waveguide basis nails photonic-crystal resonances","Slab modes unlock complete photonic-crystal resonances","New paradigm: RSE for photonic crystals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000825,"raw_usage":{"total_tokens":3603,"prompt_tokens":934,"completion_tokens":2669,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":2606}},"tokens_in":550,"tokens_out":2669,"duration_ms":22483,"temperature":1.0,"reasoning_tokens":2606,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:31:20.865714+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a slab with $\\epsilon = 6$, $b = 0.95a$, $d = 2\\pi/5$, $\\beta = 1$, and $p = 0$, compute the resonant frequencies with a scattering-matrix code using more Bragg channels than $M = 5$; if the PC-RSE results near $\\mathrm{Re}(\\omega a) = 5$ do not keep moving toward those values as $N$ grows, or if the scattering-matrix calculation finds modes the PC-RSE's fixed-$N$ diagonalization misses, the claim of guaranteed completeness fails in practice.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is the Mittag-Leffler theorem that justifies the Green's function expansion."},{"cited_title":"Vollmer et al., Appl","cited_arxiv_id":null,"evidence_quote":"Introduced the resonant-state expansion mapping Maxwell's equations to a linear eigenvalue problem, which this paper extends to periodic slabs."},{"cited_title":"Gamow, Z","cited_arxiv_id":null,"evidence_quote":"Showed how branch cuts of the Green's function are included via discretized cut modes in inhomogeneous waveguides."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provided the RSE for a homogeneous planar slab in the k-representation and the normalization used for basis states."},{"cited_title":"Rosenblit, P","cited_arxiv_id":null,"evidence_quote":"Established the convergence scaling and matrix formulation for open 2D systems used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Treated cut discretization for 2D open systems, adapted here for the slab cuts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gave the generalized normalization and operator form for magnetic and bi-anisotropic systems used in Eq. (20)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is the scattering-matrix method used as the numerical reference for verification."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Represents the incomplete guided-mode expansion that motivates adding leaky modes and cut modes to the basis."}],"review_version":1}