{"id":"15307648-b712-48c8-adad-f7b48f8d8d01","arxiv_id":"1908.04571","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Closed-form scattering coefficients and Faraday rotation are derived for subwavelength arrays of magnetically-biased graphene ribbons, matching full-wave simulations.","lead":"This paper derives closed-form formulas for the reflection and transmission of light through a periodic array of graphene ribbons placed in a magnetic field. The formulas are fast to evaluate and match full-wave simulations, which could speed up the design of tiny terahertz devices such as Faraday rotators and isolators.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The perturbation step from Eq. (18) to Eq. (24) makes the Rxx resonance formally D-independent; a D-sweep at fixed w is needed to test whether this neglected periodic-image term actually matters.","rationale":"The central claim that Eqs. (40)-(43) give accurate closed-form scattering for subwavelength magnetically-biased graphene ribbon arrays is supported by the CST comparison, and the authors are unusually explicit about the subwavelength, low-fill-factor, and low-Fermi-level limitations. The reader's weakest-assumption identification is essentially correct. My concern sharpens it: the specific consequence of assuming the array does not change the eigenvalue equation is that the Rxx pole is formally independent of D, and this is not directly validated by Fig. 4, which changes w rather than D. This is a real soft spot but not a demonstrated failure: the quoted errors at the tested parameters are small, the method is explicitly scoped to the subwavelength regime, and the perturbation assumption is acknowledged in the paper. Thus I would not overturn the ACCEPT verdict; I keep it UNCHANGED while recommending the D-sweep as a cheap, decisive check before relying on the formulas outside the exact geometry tested.","tokens_in":11449,"tokens_out":21939,"duration_ms":233937,"concrete_test":"Fix w = 2 micrometers with Ef = 0.5 eV, tau = 1 ps, B0 = 10 T, and use CST or an independent MoM discretization of Eq. (18) with the exact periodic kernel to find the first Rxx resonance for D = 3, 4, 5, and 8 micrometers. Compare these with Eq. (40). If the full-wave or full-operator resonance shifts by more than a few percent while Eq. (40) predicts no shift, the perturbation ansatz in Eqs. (21)-(24) is the failure point. Repeating the sweep at w/D approximately 0.3, 0.5, and 0.7 would delimit where the central claim actually holds.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing approximation is the reduction of Eq. (18) to the single-ribbon eigenproblem (24). In the quasi-static limit the integral operator in (18) contains the full periodic kernel sum over l of 1/(x-x'-lD), not just 1/(x-x'); the l != 0 terms form a smooth image kernel approximately equal to -(pi^2/3)(x-x')/D^2 plus higher orders. The paper adopts from [31] a perturbation treatment in which these image terms are assumed not to alter the eigenfunctions psi_n or eigenvalues q_n, with all array effects entering only through F and the average current. This is exactly the assumption the paper itself flags on page 5: solution (21) is based on perturbation theory in which the interaction between neighboring ribbons is assumed not to affect the eigenvalue equation (24), and increasing fill factor w/D is expected to degrade accuracy. Because q_n in (24) has no D dependence, the poles of Y_n, and hence the Rxx resonance frequencies in Eq. (40), are D-independent for fixed ribbon width. That is a concrete, falsifiable prediction of the closed forms, but the paper's fill-factor study in Fig. 4 varies w for fixed D and therefore does not isolate the D dependence. If the image kernel shifts resonance positions by more than the reported roughly 0.5-2.7 percent at the tested parameters, the central closed forms inherit a hidden error not captured by the stated k0 D much less than 1 validity condition.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytical method for computing the reflection and transmission matrices and the Faraday rotation angle of a periodic array of magnetically-biased graphene ribbons under normal incidence. Starting from integral equations for the induced surface currents, the authors approximate the periodic Green's function by its constant leading term for the y-component and reduce the x-component equation to the single-ribbon eigenproblem of Ref. [31] with a modified forcing term. This yields the closed-form expressions (40)-(42) for the scattering matrices and (43) for the Faraday rotation. The results are validated against CST full-wave simulations for several parameter sets, with resonance-frequency errors from 0.5% at the first resonance to 2.71% at the fourth. The method is compared with effective-medium theory and a semi-analytical approach, demonstrating that the new closed forms capture resonant features that the earlier methods miss.","tokens_in":11699,"tokens_out":7771,"duration_ms":75899,"significance":"If the results hold, the paper provides a fast, parameter-free analytical tool for designing terahertz devices based on magnetically-biased graphene ribbon arrays, avoiding the computational cost of full-wave simulations and the inaccuracies of effective-medium models. The paper is honest about its limitations, explicitly stating that the perturbation solution (21) assumes that inter-ribbon interactions do not affect the eigenvalue equation (24) and that accuracy degrades with increasing fill factor. The explicit error quantification against CST is a strength, as is the clear statement of the subwavelength validity regime. The main gap is the lack of a D-sweep test, which the major comment addresses.","major_comments":[{"comment":"The eigenvalue equation (24) used in the perturbation solution (21) contains no dependence on the array period D. Since the poles of Yn in (22) are set by qn, the Rxx resonance frequencies predicted by (40) are formally independent of D for a fixed ribbon width w. This is a falsifiable prediction of the closed-form solution, but the paper's validation does not isolate it: Fig. 4 keeps D = 10 µm and varies w, so it tests the effect of the gap D - w but not the D-dependence. To support the central claim that (40)-(42) are accurate for subwavelength arrays, the authors should add a D-sweep at fixed w (e.g., w = 2 µm, D = 3-12 µm, with the same Ef, tau, and B0) and compare the resonance frequencies from (40) with CST, reporting the error as in Fig. 4. This is needed because the neglected image terms in the kernel of (18) could shift the resonances by more than the 0.5-2.71% errors seen in Fig. 2, without violating the stated k0D << 1 condition.","section":"II.C, Eq. (24)"}],"minor_comments":[{"comment":"The perturbation solution (21) is imported from Ref. [31] without an explicit derivation or a statement of the small parameter that controls the perturbation. Please provide a brief justification or a precise reference to the relevant equations in [31] so the reader can assess the conditions under which the interaction between ribbons does not affect the eigenproblem (24).","section":"II.C, text after Eq. (20)"},{"comment":"The caption and the text should specify how the relative error in the resonance frequency is defined (e.g., |f_theory - f_FIT|/f_FIT) and for which resonance (first, n = 1) it is computed.","section":"III.A, Fig. 4 caption"},{"comment":"The derivation of the expansion in k0D and k0|x| is compressed; please show the intermediate steps for the first integral in (A.5) to demonstrate explicitly that it is of order k0D, as claimed.","section":"Appendix A, Eqs. (A.5)-(A.6)"},{"comment":"There are several minor typographical issues: 'in case of normal incidence' should be 'in the case of normal incidence' (Section II.C), 'a analytical' should be 'an analytical' (Section IV), and 'in Fig.5b' should be 'in Fig. 5(b)' (Section III.B).","section":"Throughout"}],"recommendation":"minor_revision","confidential_remarks":"The paper relies on the authors' prior work [31] for the central perturbation solution, and the novelty relative to that work is the extension to magnetic bias and the derivation of the closed-form scattering matrices. The D-sweep validation requested in the major comment would significantly strengthen the acceptance case; the current validation is suggestive but not definitive regarding the range of validity in D."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a real analytical result: closed-form reflection/transmission matrices and Faraday rotation for a periodic array of magnetically biased graphene ribbons, derived from integral equations under quasi-static and perturbation assumptions, and validated against CST with no fitted parameters. Second, the weakest piece is the perturbation step borrowed from the authors' earlier paper [31], and the paper never runs the one numerical test that would directly exercise it—a D-sweep at fixed ribbon width.\n\nThe paper does well what it sets out to do. It extends the single-ribbon integral-equation approach to the anisotropic conductivity tensor with off-diagonal terms, keeps the algebra transparent, and produces equations (40)–(43) that a designer can evaluate in milliseconds. The comparison against the effective-medium approach [26] is fair, and the point that the semi-analytical method in [20] misses the resonances is backed by the figures. The error growth from 0.5% at the first resonance to 2.71% at the fourth is consistent with the stated subwavelength assumption, and the paper says plainly that higher fill factors will hurt accuracy.\n\nNow the soft spots, in proportion. The perturbation solution (21) is imported from [31] without re-derivation. That would be fine if the underlying eigenvalue equation (24) were merely a standard reference, but here that single-ribbon equation carries the entire array interaction, and it has no D dependence. Consequently the closed forms predict Rxx resonance frequencies that are formally independent of D for fixed w—a concrete, falsifiable prediction. The paper's Fig. 4 sweeps fill factor by varying w with D fixed, so it never isolates D. The periodic Green's function's image terms, expanded in the appendix, could shift those resonances by more than the quoted percent-level errors in some parameter regimes. That is not a fatal flaw—the method clearly works in the studied subwavelength cases—but it is the difference between a design tool with a known accuracy domain and one whose domain boundary is untested.\n\nWho is this for? People designing THz Faraday rotators, isolators, or circulators with patterned graphene, and anyone who wants quick design equations instead of running CST every time. It deserves a serious referee. I would recommend a conditional accept: require a D-sweep at fixed w, plus one sentence acknowledging what the O(k0D) corrections in (13) would do to the resonance positions.","headline":"Solid analytical extension of graphene-ribbon scattering to the magnetically biased case, with honest limits and a D-sweep gap worth checking before relying on the closed forms.","tokens_in":12281,"tokens_out":4355,"would_cite":true,"duration_ms":44655,"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":"Closed-form expressions give the full reflection, transmission, and Faraday rotation of an array of magnetically biased graphene ribbons under normal incidence.","keywords":["graphene ribbons","magnetically biased graphene","Faraday rotation","terahertz metasurface","subwavelength array","surface conductivity tensor","reflection and transmission coefficients","magnetoplasmon resonance"],"falsifier":"Run a rigorous full-wave simulation of the same structure at a fill factor near 0.9 (for example $D=10\\,\\mu\\mathrm{m}$, $w=9\\,\\mu\\mathrm{m}$, $E_f=0.3$ eV, $\\tau=1$ ps, $B_0=7.5$ T) and compare the first reflection-resonance frequency and the Faraday rotation angle with equations (40)-(43); if the predicted resonance shifts by more than the few percent the paper reports at lower fill factors, or if extra resonances appear, the constant-term truncation of the periodic Green's function is the point of failure.","tokens_in":11212,"feed_emoji":"🧲","tokens_out":7882,"duration_ms":72824,"temperature":0.7,"pith_summary":"This paper aims to reduce the electromagnetic response of a subwavelength periodic array of magnetically biased graphene ribbons to a small set of closed-form expressions: the reflection and transmission matrices (40)-(42) and the Faraday rotation angle (43). If correct, a designer can predict how such an array transmits, reflects, and rotates the polarization of a normally incident terahertz wave without running a full-wave simulation. The derivation starts from integral equations for the surface currents, replaces the periodic Green's function by its dominant subwavelength term, and treats the array as a perturbation of a single ribbon's current eigenmodes. The resulting formulas reproduce the resonant spectral features that simpler effective-medium models miss, with errors that the paper measures as small in the subwavelength regime.","feed_headline":"Closed-form formulas for Faraday rotation in graphene arrays","feed_subtitle":"Analytic reflection and transmission matrices capture magnetoplasmon resonances that effective-medium models miss.","key_machinery":"The load-bearing object is the quasi-static periodic Green's function, whose leading term for a subwavelength array is $G(x)\\simeq 1/(2jk_0D)$, a constant independent of position. This constant converts the integral equations for the surface currents into algebraic equations for the average currents $I_x$ and $I_y$. The current profile itself is carried by the single-ribbon eigenfunctions $\\psi_n(x)$ and eigenvalues $q_n$ defined by equation (24), with the magnetic field entering through the conductivity tensor components $\\sigma_{xx}$, $\\sigma_{yy}$, and $\\sigma_{xy}$. Each eigenmode contributes an admittance $Y_n=\\sigma_{xx}(2j\\omega\\varepsilon_0/q_n)/(\\sigma_{xx}+2j\\omega\\varepsilon_0/q_n)$, so a resonance appears whenever $q_n\\sigma_{xx}+2j\\omega\\varepsilon_0$ is small; only odd modes couple to the normally incident field. The self-consistent average current then yields the closed-form reflection and transmission coefficients.","core_discovery":"The central claim is that for a subwavelength periodic array of magnetically biased graphene ribbons under normal incidence, the scattering problem has an explicit algebraic answer. With the surface conductivity tensor of biased graphene, the induced $x$-directed current is expanded in the eigenmodes $\\psi_n$ of the single-ribbon integral operator, while the $y$-directed current follows from a simple algebraic relation once the average $x$-current is known. The array enters only through a constant term in the periodic Green's function, $1/(2jk_0D)$, which makes the interaction between ribbons a self-consistent constant field rather than a mode-mixing interaction. This leads to equations (40)-(42) for the four reflection and transmission coefficients and equation (43) for the Faraday rotation angle. The paper argues these expressions agree well with full-wave simulations and, unlike the effective-medium approximation, capture the magnetoplasmon resonances of the array.","pith_inferences":["Because the final formulas are algebraic, they could be inverted: given a target Faraday rotation at a chosen terahertz frequency, one could solve for the bias field, Fermi level, or ribbon width directly rather than by parameter sweeps.","The same perturbation structure, a constant leading term in the periodic Green's function plus single-ribbon eigenmodes, may extend to oblique incidence or to arrays on dielectric substrates, where the constant term would be modified by the substrate's response.","The explicit resonance condition suggests an experimental test: a terahertz transmission measurement through a biased ribbon array should show a rotation peak at the predicted width-dependent frequency, which would also reveal how far the perturbation assumption holds."],"forward_implications":["The reflection and transmission of the array can be evaluated directly from (40)-(42) in terms of material parameters $E_f$, $\\tau$, $B_0$, and geometry $w$, $D$, with no numerical discretization.","The method predicts resonant peaks in the reflection coefficients at frequencies set by the ribbon width, arising from odd current eigenmodes $n=1,3,\\ldots$, features that effective-medium theory cannot reproduce.","The Faraday rotation angle follows from (43) and shows large values near resonance; increasing the magnetic bias increases the rotation and shifts the resonance to higher frequencies.","The method's accuracy is bounded by the subwavelength condition $k_0D\\ll 1$ and by the perturbation assumption that inter-ribbon interaction does not alter the single-ribbon eigenmodes, so accuracy degrades as frequency or fill factor rises."],"supporting_citations":[{"why":"Provides the anisotropic surface-conductivity tensor for magnetically biased graphene that forms the material model throughout the derivation.","marker":"[13]"},{"why":"Supplies the Drude-like approximation for the conductivity tensor components used to simplify the equations.","marker":"[27]"},{"why":"Provides the dyadic Green-function integral-equation formulation connecting surface currents to scattered fields.","marker":"[30]"},{"why":"Gives the single-ribbon integral equation and eigenmode expansion that the array solution perturbs.","marker":"[31]"},{"why":"Supplies the integral identities used in the appendix to reduce the periodic Green's function to its leading constant term.","marker":"[33]"},{"why":"Provides the Faraday-rotation formula (43) and a semi-analytical benchmark whose predictions the new method improves upon.","marker":"[20]"},{"why":"Supplies the effective-medium model used as the comparison that fails to capture the resonant response.","marker":"[26]"}],"fun_headline_variants":["Exact scattering formulas for biased graphene ribbon arrays","Analytic route to Faraday rotation in graphene arrays","Graphene array analysis made exact and fast","Closed-form currents for magnetically biased graphene ribbons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the array is deeply subwavelength ($k_0D\\ll 1$) and that neighboring ribbons interact only through a constant term in the periodic Green's function, leaving the single-ribbon current eigenmodes unchanged; if the period approaches the wavelength or the fill factor becomes large, this premise breaks down.","fun_headline_variants_meta":{"raw":{"variants":["Exact scattering formulas for biased graphene ribbon arrays","Analytic route to Faraday rotation in graphene arrays","Graphene array analysis made exact and fast","Closed-form currents for magnetically biased graphene ribbons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000389,"raw_usage":{"total_tokens":2027,"prompt_tokens":899,"completion_tokens":1128,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":1069}},"tokens_in":515,"tokens_out":1128,"duration_ms":9250,"temperature":1.0,"reasoning_tokens":1069,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:38:47.239616+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a rigorous full-wave simulation of the same structure at a fill factor near 0.9 (for example $D=10\\,\\mu\\mathrm{m}$, $w=9\\,\\mu\\mathrm{m}$, $E_f=0.3$ eV, $\\tau=1$ ps, $B_0=7.5$ T) and compare the first reflection-resonance frequency and the Faraday rotation angle with equations (40)-(43); if the predicted resonance shifts by more than the few percent the paper reports at lower fill factors, or if extra resonances appear, the constant-term truncation of the periodic Green's function is the point of failure.","supporting_citations":[{"cited_title":"Dyadic green’s functions for an anisotropic, non- local model of biased graphene,","cited_arxiv_id":null,"evidence_quote":"Provides the anisotropic surface-conductivity tensor for magnetically biased graphene that forms the material model throughout the derivation."},{"cited_title":"Faraday effect in graphene enclosed in an optical cavity and the equation of motion method for the study of magneto-optical transport in solids,","cited_arxiv_id":null,"evidence_quote":"Supplies the Drude-like approximation for the conductivity tensor components used to simplify the equations."},{"cited_title":"Tai, Dyadic Green functions in electromagnetic theory","cited_arxiv_id":null,"evidence_quote":"Provides the dyadic Green-function integral-equation formulation connecting surface currents to scattered fields."},{"cited_title":"Analytical modeling of graphene ribbons as optical circuit elements,","cited_arxiv_id":null,"evidence_quote":"Gives the single-ribbon integral equation and eigenmode expansion that the array solution perturbs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the integral identities used in the appendix to reduce the periodic Green's function to its leading constant term."},{"cited_title":"Faraday rotation due to excitation of magnetoplasmons in graphene microribbons,","cited_arxiv_id":null,"evidence_quote":"Provides the Faraday-rotation formula (43) and a semi-analytical benchmark whose predictions the new method improves upon."},{"cited_title":"Magnetically-biased graphene-based hyperbolic metasurfaces,","cited_arxiv_id":null,"evidence_quote":"Supplies the effective-medium model used as the comparison that fails to capture the resonant response."}],"review_version":1}