{"id":"e514f230-2298-4fa1-8f30-605c96fd8b80","arxiv_id":"2608.03759","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"AFLOW-EMERALD is a new open-source solver that unifies scattering-matrix, RCWA, and complex-k plane-wave methods for layered metamaterials.","lead":"This paper presents AFLOW-EMERALD, an open-source Python and MATLAB framework for simulating light propagation in layered optical and plasmonic materials. It combines scattering-matrix, rigorous coupled-wave, and plane-wave expansion methods to compute spectra, fields, and band structures for photonic crystals and hyperbolic metamaterials.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (24) is an unverified finite-stack quantization rule applied to lossy media, with an indexing inconsistency; it can mislabel VPP mode orders even if the core solvers are correct.","rationale":"The reader's weakest-assumption analysis is correct: Eq. (24) is the least-supported link in the paper's central mode-identification workflow. The rest of the theoretical framework (SMM, PWE, RCWA) follows standard formulations and the examples are internally plausible, so I do not see a more fundamental correctness flaw that would justify a stronger verdict. The main caveat is that Eq. (24) is exact only in the lossless limit and for real k_z; applying it to a lossy HMM without a derivation or a validation study is a genuine soft spot. I also note the indexing mismatch between ℓ=1,...,N in Eq. (24) and m=0,...,N−1 in the text, which makes the overlay ambiguous. This does not change the reader's CONDITIONAL verdict, because it is an addressable, example-specific issue rather than a reason to reject the software claim outright. A direct numerical test of dip energies versus isoline intersections would settle whether the concern lands.","tokens_in":17360,"tokens_out":39463,"duration_ms":327331,"concrete_test":"Run the released code on the Ag/TiO2 HMM with N=3, 6, and 9 bilayers at the Fig. 4c parameters; for each reflectance dip in rta_energy, count the field nodes in the fields output and compare the dip energy with the intersection of the complex-k band map and the Eq. (24) isoline k_z=mπ/(NΛ_z). Repeat for the lossless TiO2/SiO2 N=5 stack, where the rule should be exact. If the lossy dip energies deviate from the isoline intersections by more than the line spacing π/(NΛ_z), or the node counts disagree with m, then the finite-stack mode-overlay is not reliable for design.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV.B.4 introduces Eq. (24), k_z = ℓπ/(NΛ_z) for ℓ=1,...,N, as the resonance condition for finite stacks, without derivation or citation. It is then used in Section V.D to overlay the VPP isolines in Fig. 7b and, together with the field profile, to label the mode in Fig. 6b as a first-order VPP. The exact lossless finite-stack condition is sin(N k_z Λ_z)=0, which yields real k_z=mπ/(NΛ_z); but the paper applies Eq. (24) to a lossy, dispersive Ag/TiO2 HMM, where k_z is complex and the real-part isolines need not coincide with actual RTA resonances. The indexing is also inconsistent: Eq. (24) states ℓ=1,...,N, while the text says the mode index m ranges from 0 to N−1 and plots 'first five' values for N=6. If the overlay is offset, the advertised 'which VPP order is excited' analysis in Figs. 6b and 7b is mislabeled even though the underlying RTA and field calculations may still be correct.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript describes AFLOW-EMERALD, an open-source Python/MATLAB package for simulating electromagnetic propagation in finite and periodic layered structures. The code combines a scattering-matrix method (SMM) with rigorous coupled-wave analysis (RCWA) for laterally patterned systems and plane-wave expansion (PWE) with a complex-k formulation for photonic bandstructures of dispersive, lossy media. A YAML workflow imports dielectric data from experimental, literature, or first-principles sources. The paper presents three examples: a TiO2/SiO2 photonic crystal (RTA spectra and bands), a Ag/TiO2 hyperbolic metamaterial with a grating coupler (RTA, fields, complex-k bands), and an air/Ag surface plasmon polariton (angle-resolved reflectance compared with experiment). The finite-stack mode identification relies on a proposed resonance quantization rule, Eq. (24), which is used to overlay volume-plasmon-polariton (VPP) mode isolines and to label mode orders in the field profiles.","tokens_in":17622,"tokens_out":4532,"duration_ms":42732,"significance":"If the implementation is correct, AFLOW-EMERALD would be a practically useful open-source tool that connects realistic finite-stack optics with infinite-periodic bandstructure analysis, including complex-k dispersion in lossy media. The paper's strengths include the modular architecture, the use of external dielectric datasets with no parameter fitting, the demonstration of a Kretschmann SPP reflectance minimum that matches experimental data, and the public availability of the code. The main advertised capabilities beyond standard RTA, however, rest on the finite-stack quantization rule of Eq. (24), which is asserted without derivation or validation, and the paper lacks independent benchmarks for the RCWA, PWE, and complex-k solvers. These gaps are significant because the central claim is that the released code reliably supports material-geometry co-design.","major_comments":[{"comment":"The finite-stack quantization rule is introduced without derivation or citation and contains an indexing inconsistency. Eq. (24) states k_z = ℓπ/(NΛ_z) for ℓ = 1,...,N, but the surrounding text and Figure 7b refer to \"first five k_z resonance values (m = 0–4)\" for N = 6, with \"mode index m ranging from 0 to N–1.\" These two prescriptions are different: Eq. (24) excludes m = 0 and includes ℓ = N, whereas the text and figure exclude ℓ = N and include m = 0. Moreover, for the lossy, dispersive Ag/TiO2 system used in Section V.D, k_z is complex; the real-part isolines given by Eq. (24) need not coincide with the actual resonances of a finite stack. In the lossless case, the finite-stack condition is sin(N k_z Λ_z) = 0, which gives real k_z = mπ/(NΛ_z) for integer m, not the set in Eq. (24). Because this rule is used to overlay the VPP isolines in Figure 7b and to label the mode in Figure 6b as a first-order VPP, the advertised \"which mode is excited\" analysis is not reliable as presented. The authors should provide a derivation or correct the condition, fix the indexing, and demonstrate with RTA spectra and field profiles that the isolines track the resonances in the lossy case.","section":"Section IV.B.4, Eq. (24) and Section V.D"},{"comment":"The external validation is limited to a single experimental comparison: the air/Ag SPP reflectance dip at one angle and energy in Figure 5. No analytic benchmark (e.g., Fresnel reflection from a single slab), no independent-solver comparison for the RCWA grating calculation or the complex-k PWE bandstructure, and no convergence study with respect to the plane-wave truncation parameter (halfnpw) are provided. The performance claims in Section IV.C (0.03 s per energy point, 10–100 minutes for 2D PBS maps, up to one order-of-magnitude GPU speedup, memory below 1–2 GB) are stated without hardware details, basis sizes, or convergence criteria. For a software paper whose central claim is that the released code correctly and stably implements these solvers for design use, these gaps should be filled: analytic checks, comparisons with established codes or published grating efficiencies, and convergence tests with respect to basis size.","section":"Section V and Section IV.C"}],"minor_comments":[{"comment":"The word \"Bruillouin\" in the inset label is misspelled; it should be \"Brillouin.\"","section":"Figure 7 caption"},{"comment":"The text says \"For a stack of N layers,\" but Section V.D uses N = 6 bilayers (periods). Clarify whether N is the number of periods or the number of individual layers, since Eq. (24) depends on this distinction through Λ_z.","section":"Section IV.B.4"},{"comment":"For reproducibility, the performance measurements should specify the processor, memory, MATLAB engine version, and the convergence criterion used to select the plane-wave basis size in each reported timing.","section":"Section IV.C"},{"comment":"The interpolation of dielectric-function data is described only as \"polynomial interpolation\"; specify the polynomial order and the treatment of data points outside the supplied energy range.","section":"Section IV.A"},{"comment":"The repository is given as a GitHub URL without a version tag or DOI; a versioned release would improve reproducibility of the specific results shown in the manuscript.","section":"Data and Code Availability"}],"recommendation":"major_revision","confidential_remarks":"The main risk is Eq. (24): it is unverified, has an indexing inconsistency, and is used to label VPP modes. This is fixable but central, so I recommend major revision rather than rejection. The lack of benchmarks is also significant for a software paper, but the underlying SMM/RCWA/PWE formulation appears standard and the SPP comparison is a positive indicator. I would be satisfied if the authors supply a derivation or a clearly labeled heuristic for Eq. (24), correct the mode indexing, and add at least one independent-solver or analytic benchmark for the complex-k and RCWA paths."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, read this if you do photonic/plasmonic simulations and want a convenient open-source package. The paper is a software-tool paper: it wraps standard SMM/RCWA/PWE methods into one Python/MATLAB framework, with a YAML-driven interface, complex-k photonic bandstructures, and finite-stack resonance overlays. The open-source release, documentation, and example scripts are real and will genuinely help people who want quick RTA spectra or field profiles for layered heterostructures. The physics examples are sensible, the SPP angle matches the experimental reference, and the authors use external dielectric data without fitting parameters, so the workflow is self-contained. Credit where earned: this is a useful integration of known techniques, not a new method.\n\nThe main soft spot is Eq. (24), the finite-stack quantization condition k_z = l*pi/(N*Lambda_z). It appears in Section IV.B.4 without derivation or citation and is then used to overlay VPP mode isolines in Fig. 7b and label the mode in Fig. 6b as first-order. For a lossy, dispersive HMM, k_z is complex; using real-part isolines to identify resonances is at best an approximation. The indexing also seems off: Eq. (24) says l=1,...,N, but the text says mode index m=0,...,N-1 and the figure plots \"first five\" for N=6. If this rule is not justified, the advertised capability of assigning VPP orders in finite stacks is unreliable, even if the underlying spectra and field calculations are correct. The stress-test note captures this accurately; it lands. A proper fix would be to derive the quantization for finite lossy stacks, or to compare the isoline crossings against field-node counts from an independent solver.\n\nOther issues are minor in comparison: only one experimental validation, no independent-solver benchmarks, and performance numbers given without hardware or convergence details. The MATLAB dependency will annoy some open-source users, but the package still runs with standard Python dependencies. None of these are fatal; they are standard referee requests for a tool paper.\n\nWho is this for: anyone doing layered photonic crystals, hyperbolic metamaterials, or plasmonic multilayers who wants a ready-made solver rather than writing their own. It deserves a serious referee, with the request to fix or properly frame Eq. (24) and add at least one independent benchmark. I would support publishing it after those revisions.\n\nI would not personally cite it in the next year unless I start working with that exact material stack, but I would bring it to a reading group discussing metamateiral design tools.","headline":"A genuinely useful open-source tool for layered metamaterial simulations, with one load-bearing assumption (Eq. 24) that needs a derivation or a caveat before the mode-labeling claims can be trusted.","tokens_in":18115,"tokens_out":1592,"would_cite":false,"duration_ms":16899,"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":"The paper presents an open-source electromagnetic solver that computes optical spectra, field maps, and complex-k photonic band structures of layered metamaterials within one framework.","keywords":["electromagnetic simulation","scattering-matrix method","plane-wave expansion","rigorous coupled-wave analysis","photonic band structure","complex-k dispersion","hyperbolic metamaterials","plasmon polaritons"],"falsifier":"A concrete check is to take Ag/TiO2 stacks with $N=2,4,6,8$ periods, compute full-field magnetic profiles at the energies where the quantization rule predicts resonance crossings, and count the field nodes inside the stack. If the node counts do not follow the predicted mode order, or if the predicted crossings do not line up with dips in the reflectance spectrum, the quantization rule is the point of failure.","tokens_in":17188,"feed_emoji":"💡","tokens_out":13353,"duration_ms":113296,"temperature":0.7,"pith_summary":"This paper presents an open-source computational framework that simulates light propagation through layered and periodically patterned (meta)materials, from finite stacks of a few films to infinite periodic crystals. It claims that one unified solver, built on a scattering-matrix method, plane-wave expansion, and rigorous coupled-wave analysis, can reliably produce optical spectra, spatial field distributions, and photonic band structures, including complex-k dispersion in lossy and dispersive media. The intended payoff is a direct bridge between material data and device geometry, so that dielectric photonic crystals, plasmonic multilayers, and hyperbolic metamaterials can be co-designed rather than tuned by trial and error. The paper tests the claim on a TiO2/SiO2 photonic crystal, an Ag/TiO2 hyperbolic multilayer with a grating coupler, and an air/silver surface-plasmon interface, where the computed resonance angle matches experiment.","feed_headline":"Open-source solver maps light modes in layered metamaterials","feed_subtitle":"One code bridges finite-stack spectra, field maps, and photonic band structures.","key_machinery":"The machinery is the scattering-matrix method with auxiliary zero-thickness vacuum gaps inserted between physical layers, giving each layer a self-contained scattering matrix; the layer matrices are cascaded through the star product into a global S-matrix. Plane-wave expansion treats the infinite periodic problem, and for absorbing dispersive media the code uses an inverse-dispersion, complex-k formulation that solves for $\\beta_z$ at fixed frequency instead of solving for frequency at fixed wavevector. Laterally patterned layers enter through rigorous coupled-wave analysis, which expands fields and permittivity in Floquet-Bloch harmonics, and the finite-stack mode selector is the quantization rule $k_z = \\ell\\pi/(N\\Lambda_z)$, whose isolines are drawn over the band-structure map to mark excitable resonances.","core_discovery":"The central discovery claimed is that a single modular implementation can treat finite and infinite layered structures on equal footing: the scattering-matrix method gives numerically stable reflectance, transmittance, absorptance, and field maps for realistic finite stacks, while plane-wave expansion, extended to a complex-k inverse-dispersion formulation, gives the Bloch modes of the corresponding infinite periodic medium. The two sides are connected by a finite-stack quantization condition, $k_z = \\ell\\pi/(N\\Lambda_z)$ for $\\ell=1,\\dots,N$, which selects which Bloch modes a stack of $N$ periods can actually support. Overlaying these $k_z$ isolines with the in-plane momenta supplied by a grating labels each resonance, for instance identifying a reflectance dip as a volume plasmon-polariton of a given order, and the field profiles confirm the mode order by the number of nodes. The paper argues that with this link, spectra, near-field maps, and band dispersion form one interpretable description for designing layered metamaterial devices.","pith_inferences":["A natural test beyond the paper's examples is to compare the predicted mode orders from the quantization rule against full-field node counts across different numbers of periods; this would reveal whether the rule is exact or merely a convenient approximation.","The same framework could be extended to anisotropic, magneto-optical, or nonlinear layers by replacing the scalar permittivity with a tensor or field-dependent response, which the modular architecture appears to allow.","The reported speed of about 0.03 seconds per energy point suggests a high-throughput screening use that the paper does not yet demonstrate.","The isoline-plus-grating-harmonic picture offers an inverse-design route: tune grating period, filling factor, and stack thickness so that a chosen grating harmonic crosses a chosen $k_z$ isoline at the target energy."],"forward_implications":["Finite multilayer spectra can be computed without the exponential instability of transfer-matrix methods, even when layers are thick, lossy, or strongly impedance-mismatched.","Complex-k band maps with grating harmonics and $k_z$ isolines overlaid identify which surface and volume plasmon-polariton modes a given finite stack can excite and assign each mode an order.","Because material data can be imported from experiment, literature, or first-principles calculations through a simple YAML workflow, the same geometry can be tested with different dielectric models.","Angle-resolved reflectance reproduces the measured surface-plasmon resonance angle for the air/silver interface, supporting use of the code for experimental design.","Released as open-source, modular software, the solver is positioned for integration into larger design pipelines for photonic crystals, plasmonic multilayers, and hyperbolic metamaterials."],"supporting_citations":[{"why":"Supplies the scattering-matrix construction used for multilayer diffraction.","marker":"[11]"},{"why":"Provides the improved scattering-matrix formulation with auxiliary gaps that yields numerical stability.","marker":"[12]"},{"why":"Defines rigorous coupled-wave analysis, which extends the method to laterally patterned gratings.","marker":"[17]"},{"why":"Introduces the complex-k inverse-dispersion treatment used for absorbing dispersive media.","marker":"[27]"},{"why":"Supplies the inverse-dispersion method for computing complex photonic band diagrams.","marker":"[28]"},{"why":"Defines the star product used to cascade layer scattering matrices into the global S-matrix.","marker":"[30]"},{"why":"Supplies the dispersive TiO2 permittivity used in the photonic-crystal RTA example.","marker":"[32]"},{"why":"Supplies the dispersive SiO2 permittivity used in the photonic-crystal RTA example.","marker":"[33]"},{"why":"Provides the experimental surface-plasmon resonance data used as the accuracy benchmark.","marker":"[36]"},{"why":"Provides the experimental silver dielectric function used as input for the plasmonic simulations.","marker":"[37]"}],"fun_headline_variants":["Solver unifies finite and infinite layered light modes","Open-source tool links spectra, fields, and band structure","AFLOW-EMERALD: one framework for layered metamaterial design","Bridging finite stacks and Bloch modes in one solver","Quantized modes tie finite stacks to infinite bands"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that in a finite stack of $N$ periods the allowed longitudinal wavevectors are exactly the discrete values $k_z = \\ell\\pi/(N\\Lambda_z)$ for $\\ell=1,\\dots,N$; the paper states this rule without derivation and uses it to label which resonances are excited and what order they have. If the rule is approximate, the mode assignments would be wrong even if every computed spectrum were correct.","fun_headline_variants_meta":{"raw":{"variants":["Solver unifies finite and infinite layered light modes","Open-source tool links spectra, fields, and band structure","AFLOW-EMERALD: one framework for layered metamaterial design","Bridging finite stacks and Bloch modes in one solver","Quantized modes tie finite stacks to infinite bands"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000399,"raw_usage":{"total_tokens":2105,"prompt_tokens":981,"completion_tokens":1124,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":1042}},"tokens_in":597,"tokens_out":1124,"duration_ms":10267,"temperature":1.0,"reasoning_tokens":1042,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:46:46.037686+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check is to take Ag/TiO2 stacks with $N=2,4,6,8$ periods, compute full-field magnetic profiles at the energies where the quantization rule predicts resonance crossings, and count the field nodes inside the stack. If the node counts do not follow the predicted mode order, or if the predicted crossings do not line up with dips in the reflectance spectrum, the quantization rule is the point of failure.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the scattering-matrix construction used for multilayer diffraction."},{"cited_title":"Figotin and I","cited_arxiv_id":null,"evidence_quote":"Introduces the complex-k inverse-dispersion treatment used for absorbing dispersive media."},{"cited_title":"Takagi, S","cited_arxiv_id":null,"evidence_quote":"Provides the experimental surface-plasmon resonance data used as the accuracy benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental silver dielectric function used as input for the plasmonic simulations."}],"review_version":3}