{"id":"e6f2ed23-0321-4122-87ed-2e360022d9cc","arxiv_id":"2507.12421","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Surface plasmon pulses in a double graphene sheet can propagate stably at fractional pulse areas that are not integer multiples of pi, enabling nanoscale coherent quantum interactions.","lead":"This paper proposes a two-layer graphene waveguide with an embedded atomic ensemble as a platform for nanoscale coherent light-matter interactions. It derives a modified area theorem for graphene surface plasmons and predicts stable pulses with fractional pulse areas, which could enable compact on-chip quantum memory.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Area fixed points are not proof of shape-stable pulses; the paper's own Sec. IV.C caveat concedes the predicted fractional pulses may decay.","rationale":"The reader's weakest-assumption correctly identifies the gap between area stability and temporal-shape stability: the paper derives an area equation but does not show that any actual pulse envelope is a stable solitary solution of the full Maxwell-Bloch equations. I agree with this reading. The paper is honest about the gap in Sec. IV.C, which supports a CONDITIONAL rather than REJECT verdict. The area-theorem derivation itself is standard and the fractional fixed points follow from the model, so the mathematical core has merit; however, the abstract's 'stable propagation' claim overreaches what the area equation alone can establish. The concrete Maxwell-Bloch simulation would settle whether the fractional-area fixed point corresponds to a true shape-preserving pulse or merely to area conservation accompanied by progressive pulse broadening and decay. Since the reader already set the verdict to CONDITIONAL and my analysis does not move it, the appropriate verdict remains unchanged.","tokens_in":20789,"tokens_out":8188,"duration_ms":95493,"concrete_test":"Run a one-dimensional time-domain simulation of the full Maxwell-Bloch system (36) for the sinh-mode small-separation configuration of Sec. IV.C (k_s L=0.5, α1 from Eq. (48), ρ0=1e25 m^-3, d21=5e-32 C.m, ω≈0.5 eV, including the stated relaxation and loss terms), launching a finite-duration sech pulse with initial θ1(0)=3π. Track the peak envelope amplitude, rms pulse duration, and pulse area for propagation distances x≥5/α1. If the area converges to Θ2≈2.46π while the pulse broadens monotonically and the peak decays, the claim of stable fractional-area propagation is not supported; if a shape-preserving envelope emerges, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that graphene SP pulses propagate stably at fractional pulse areas is not established. The modified area theorem of Eqs. (38), (46), and (50) is obtained by integrating the Maxwell-Bloch equations over time, so it constrains only the pulse area θ(x)=R0∫a(x,t)dt. A fixed point of the area ODE is a necessary condition for a stable solitary pulse, but not sufficient; the full Maxwell-Bloch dynamics must also preserve the envelope shape. No such solution is constructed analytically, and no numerical simulation of the full system (36) is reported. The paper itself states in Sec. IV.C that a stable value of the pulse area may not lead to the appearance of graphene plasmon pulses with a stable temporal shape, that the pulse duration will increase to conserve the area, and that this will increase relaxation and eventually cause decay. This is a direct admission that the abstract's 'stable propagation of isolated SP pulses at fractional pulse area values' remains unsubstantiated. The quoted fractional values (Θ2≈2.46π for the sinh mode and Θ2=8π/3 for e^{-ksL}=1/2 in Sec. V) are fixed points of the pulse-area equation only.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a two-graphene-sheet platform in which resonant two-level atoms interact with confined graphene surface-plasmon modes, and it derives a modified pulse-area theorem for this setting. The authors quantize the SP modes, compute weak-field absorption coefficients for three different spatial coupling profiles (sinh, cosh, and exponential), and then solve area-evolution equations to identify fixed points at fractional values of π, including Θ2≈2.46π for the small-separation sinh mode and Θ2=8π/3+4nπ for the exponential mode with e^{-ksL}=1/2. They further propose using these dynamics for a nanoscale CRIB-type quantum memory. The paper's headline claim is that graphene SP pulses can propagate stably at fractional pulse areas, in contrast to the free-space integer-π McCall-Hahn result.","tokens_in":20979,"tokens_out":17785,"duration_ms":170913,"significance":"If the headline claim were fully established, the fractional-area fixed points would be an interesting new nonlinear propagation regime for graphene plasmonics, and the proposed nanoscale quantum-memory cell would be a useful conceptual contribution. The derivation is largely self-contained and analytical: the mode dispersion, field quantization, and area theorem follow from the stated model rather than from fits to external data. The paper also gives closed-form weak-field absorption coefficients and makes falsifiable predictions about which sheet separations and dipole orientations support stable area evolution. However, the current evidence supports only fixed points of the pulse-area ODE, not shape-stable pulse propagation, and the paper itself concedes this limitation. The reported critical values also rely on a printed formula that is inconsistent with the preceding derivation, and the Sec. V generalization contains an internal contradiction.","major_comments":[{"comment":"The printed function F1[y] in Eq. (47), F1[y]=(3 sin y - y cos y)/y², does not follow from Eq. (45). Equation (45) gives a function proportional to (sin y - y cos y)/y², whose zeros satisfy tan y = y, i.e., y≈1.43π and 2.46π, exactly the values Θ1 and Θ2 quoted in the text and in Fig. 6. The zeros of the printed numerator 3 sin y - y cos y instead satisfy tan y = y/3, which gives different roots near 1.29π and 2.38π. In addition, the small-ksL limit of Eq. (40) contains a factor 2 that is not reflected in Eq. (45); this factor does not change the roots but does affect the coefficient α1. The displayed equations should be corrected, and Fig. 6 and the quoted critical values should be tied to the actual expression used.","section":"Sec. IV C, Eqs. (45)-(48)"},{"comment":"The abstract claims stable propagation of isolated SP pulses at fractional pulse areas, but the paper's own caveat in Sec. IV C states that a stable value of the pulse area may not lead to a stable temporal shape, that the pulse duration will increase to preserve the area, and that this will increase relaxation and eventually cause decay. The derivation establishes fixed points of the area ODE (38)/(46) only; it does not solve the full Maxwell-Bloch system (36) for the pulse envelope, nor does it construct a shape-preserving solitary solution. A fixed point of the area equation is necessary but not sufficient for a propagating soliton-like pulse. To support the central claim, the authors should either present numerical solutions of Eq. (36) showing shape-stable propagation at the predicted fractional areas, or explicitly revise the abstract and conclusions to claim area stability rather than stable propagation of pulses.","section":"Sec. IV C and Abstract"},{"comment":"The generalization to e^{-ksL}=1/(2m) is internally inconsistent with the e^{-ksL}=1/2 case. The text states that for e^{-ksL}=1/2 the stable points are Θ2=2π+2π/3+4nπ=8π/3+4nπ, which is consistent with Eq. (54). It then states that for e^{-ksL}=1/(2m) one finds stable fractional values Θ2=(2m+1)π/3. For m=1 this gives Θ2=π, not 8π/3. The inconsistency is not a harmless typo: for m=2, the claimed value Θ2=5π/3 is not a zero of F3 in Eq. (51). The stable roots of sin((1+q)θ/2)sin((1-q)θ/2) with q=1/(2m) form a more complicated family that depends on m through the sign of the derivative at each root. Please provide the correct stability analysis for general m or remove the incorrect generalization.","section":"Sec. V, Eqs. (50)-(54)"}],"minor_comments":[{"comment":"The caption says L=2 nm, while the text just before Fig. 4 says L=5 nm; please make the value consistent.","section":"Fig. 4 caption"},{"comment":"In the middle-layer expression, the coefficient C term is written with (ˆx - ˆz iK||/k2)e^{ks(z-L)}, but the exponential and layer geometry suggest the polarization denominator should be k_s, not k_2.","section":"Eq. (2)"},{"comment":"The notation Θ2=2π+2π/3+4nπ is written with n=1,2,..., but the base value 8π/3 corresponds to n=0; please include n=0 or reindex.","section":"Sec. V"},{"comment":"The conclusion states that stable fractional-area SP pulses occur at 'arbitrary graphene sheet separations', which contradicts Sec. V and Fig. 9, where ksL=9 is reported to have no stable solutions. Please qualify the statement.","section":"Conclusion"},{"comment":"The notation k_sL appears variously as 'ksL', 'k_sL', and 'ksL'; please standardize the typesetting of scaled separation and of the functions F1, F2, F3.","section":"Throughout"},{"comment":"Several references contain 'and at al.' or 'et al' typos (e.g., Refs. [15], [18], [22], [23], [28]); a careful reference cleanup is needed.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The core difficulty is that the headline claim of stable fractional-area pulse propagation goes beyond what the area-theorem analysis establishes. The manuscript would become publishable if the authors correct Eq. (47), repair the Sec. V generalization, and either provide full Maxwell-Bloch pulse simulations or clearly restrict the claim to area stability. The paper is otherwise within the scope of physics.optics and the analytical framework is sound enough to warrant a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is an honest extension of the area theorem to two graphene sheets, and the fractional fixed points are new. But the abstract's 'stable propagation' outruns what the math supports; the authors themselves say in Sec. IV.C that a stable pulse area may not give a stable temporal shape and that the pulse will eventually decay. The stress-test note is on target.\n\nWhat is actually new: it generalizes your prior single-interface surface-plasmon area theorem [51] to the double-sheet geometry, computes the mode structure and quantization, and from the Maxwell-Bloch equations derives modified area equations for sinh, cosh, and exponential coupling. The fractional fixed points (Θ2≈2.46π for the sinh mode, Θ2=(2m+1)π/3 for exponential coupling at e^{-ksL}=1/2, etc.) are genuine and follow from the model, not fitted. The weak-field absorption coefficients of order 1 nm^{-1} with nanoscale cell sizes are plausible and useful. The paper is also unusually candid about limitations—the Sec. IV.C caveat is explicit.\n\nSoft spots, in proportion. The biggest is the one the authors flag: area fixed points are necessary but not sufficient for shape-stable solitary pulses. No full Maxwell-Bloch simulation is shown, so the claimed fractional-pulse propagation is not established. That is a load-bearing gap for the abstract, though not for the area-theorem result itself. Second, there is an apparent inconsistency between Eq. (45) and the printed F1[y] in Eq. (47); the reported Θ1, Θ2 values depend on the expression, and a careful check is needed. This may be a typo or a missing factor, but it is currently unresolved. Third, the quantum memory section is qualitative; it sketches CRIB reversibility and gives cell sizes, but no efficiency or fidelity estimates. Minor typos and notation slips are scattered.\n\nBottom line: the area-theorem derivation and the new fixed-point structure are worth refereeing. The paper should go to peer review, but the authors should be pushed to either simulate the full pulse dynamics or soften the stability claim in the abstract and conclusion. For a reader working on quantum plasmonics or area theorems, this is a useful contribution; for a reader looking for a demonstrated stable pulse, it is conditional.","headline":"Honest extension of the area theorem to double-sheet graphene yields new fractional fixed points, but 'stable propagation' is not established—the authors' own Sec. IV.C caveat says so.","tokens_in":21555,"tokens_out":2955,"would_cite":false,"duration_ms":31163,"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":"A two-graphene-sheet platform can hold surface-plasmon pulses at fractional pulse areas, and the paper derives the modified area theorem behind the prediction.","keywords":["graphene surface plasmons","modified area theorem","Maxwell-Bloch equations","fractional pulse area","resonant atomic ensemble","quantum memory","self-induced transparency","nanoscale light-matter interaction"],"falsifier":"A transmission experiment on a graphene double-sheet sample with a nanoscale gap would settle the claim: launch pulses of varying input area and measure the output area and temporal profile over several absorption lengths. If the fractional fixed points are real, output areas should lock to the predicted $\\Theta_2$ values, whereas growth of the pulse duration without area locking would falsify the propagation claim.","tokens_in":20526,"feed_emoji":"⚛️","tokens_out":9518,"duration_ms":106018,"temperature":0.7,"pith_summary":"This paper proposes a nanoscale platform in which an atomic ensemble sits between two graphene sheets and interacts coherently with graphene surface plasmons. The central result is a modified area theorem for these plasmon pulses, derived from the Maxwell-Bloch equations, whose stable fixed points occur at fractional values of $\\pi$ rather than at the integer multiples required in free space. The authors find stable pulse areas of about $2.46\\pi$ in the small-separation $\\sinh$ coupling regime and $(2m+1)\\pi/3$ for exponential coupling when $e^{-k_s L}=1/2$. If this holds, compact quantum memory cells only a few nanometers long could be integrated with other graphene-based quantum devices. The paper itself cautions that a stable pulse area may not guarantee a stable temporal pulse shape, since preserving the area can lengthen the pulse and increase relaxation.","feed_headline":"Graphene plasmons hold stable pulses at fractional area values","feed_subtitle":"A double-sheet area theorem predicts pulse locking away from integer multiples of pi, enabling nanoscale quantum memory.","key_machinery":"The central machinery is the modified area theorem for graphene surface plasmons, obtained by solving the Maxwell-Bloch equations for the plasmon envelope in the double-sheet waveguide. It reduces to the propagation equation $\\left(\\partial_x + \\gamma_w/2v_g\\right)\\theta(x) = -(\\alpha/2)F[\\theta(x), k_s L]$, where the dimensionless function $F$ encodes the transverse coupling profile: $\\sinh[k_s(z-L/2)]$, $\\cosh[k_s(z-L/2)]$, or $e^{-k_s z}$ depending on dipole orientation and mode symmetry. Stable fixed points are zeros of $F$ with positive slope, and the fractional values of $\\theta$ at those points are the paper's main prediction. Graphene conductivity and sheet separation enter through $k_s L$, the absorption coefficient $\\alpha$, and the group velocity, so the pulse-area behavior is tunable by gating and geometry.","core_discovery":"On the paper's own terms, the discovery is that the interaction of graphene surface plasmons with an inhomogeneously broadened two-level atomic ensemble is governed by an area theorem with spatially inhomogeneous coupling, so stable propagation is no longer tied to integer multiples of $\\pi$. Solving the Maxwell-Bloch equations for three coupling geometries ($\\sinh$, $\\cosh$, and exponential in the coordinate across the graphene gap) yields stable fixed points at fractional pulse areas, for example $\\Theta_2 \\simeq 2.46\\pi$ for the $\\sinh$ mode at small sheet separation and $\\Theta_2 = (2m+1)\\pi/3$ for exponential coupling with $e^{-k_s L}=1/2$. The fractional value reflects the continuous distribution of coupling strengths across the atomic sample: atoms near the sheets see a larger effective pulse area, while atoms in the middle see a smaller one. The authors also show that the double-sheet structure can implement a controlled reversible inhomogeneous broadening (CRIB) quantum memory whose cell size is estimated at $5$-$25$ nm, with the storage lifetime set by atomic coherence rather than by plasmon loss.","pith_inferences":["If area-locking survives full pulse-shape dynamics, input pulses launched above the unstable threshold should converge to the fractional value after propagation, which could be tested by measuring transmitted pulse area versus input area.","The fractional fixed points are effectively a map of the transverse field profile; measuring them for different dipole orientations could serve as a nanoscale probe of the plasmon mode structure.","The same inhomogeneous-coupling mechanism should appear in other double-interface polariton waveguides, so the fractional-area phenomenon may be a general feature of guided polaritons rather than something specific to graphene."],"forward_implications":["If the area theorem holds, graphene surface plasmon pulses can propagate over many absorption lengths with their area locked to a fractional value such as $\\Theta_2 \\simeq 2.46\\pi$, in contrast to the free-space $2\\pi$ soliton condition.","The predicted absorption coefficients are of order $1\\ \\mathrm{nm}^{-1}$, two to three orders of magnitude above graphene losses, making coherent interaction and quantum memory feasible on nanoscale lengths.","A CRIB-based quantum memory cell for surface plasmons would be only about $5$-$25$ nm long, with storage lifetime governed by the atomic ensemble rather than by plasmon decay.","For large sheet separations ($k_s L \\ge 9$) the area equation has no stable fixed points, so stable fractional pulses and the associated memory operation require a tightly confined geometry."],"supporting_citations":[{"why":"Supplies the graphene plasmonics platform properties, including tunable conductivity, strong field confinement, and low losses, on which the proposed structure is built.","marker":"[20]"},{"why":"Establishes low graphene plasmon losses below the Fermi energy, the premise that lets the authors neglect losses and quantize the SP field.","marker":"[25]"},{"why":"Supplies the Maxwell-Bloch and pulse-area derivation for two-level atoms used to obtain the modified area theorem.","marker":"[32]"},{"why":"Supplies the free-space area theorem and the integer-$\\pi$ stability condition that the graphene modification generalizes to fractional values.","marker":"[37]"},{"why":"Supplies the single-interface surface-plasmon area theorem and its unstable $2\\pi$ solution, which the two-sheet geometry extends and stabilizes.","marker":"[51]"},{"why":"Supplies the double-interface surface-polariton dispersion relation used to define the graphene SP modes.","marker":"[53]"},{"why":"Supplies the coupled-mode transcendental equation for double charge-sheet structures that fixes the SP wave numbers.","marker":"[54]"},{"why":"Supplies the validity condition $\\delta t_s \\gg \\lambda/c$ for applying the area-theorem derivation to the short plasmon pulses.","marker":"[61]"},{"why":"Supports keeping only surface-plasmon modes in the Hamiltonian by showing atomic ensembles radiate into SP modes far more strongly than into free-space modes.","marker":"[62]"}],"fun_headline_variants":["Fractional pulse areas stabilize graphene plasmon pulses","Graphene-area theorem breaks the pi rule for plasmon pulses","Nanoscale quantum memory with graphene plasmons at fractional pulse areas","Stable plasmon pulses at fractional areas enable quantum memory","Crossing the gap: graphene plasmons defy the pi-pulse rule"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a stable fixed point of the pulse-area equation means the pulse actually propagates stably in shape; the paper itself cautions (Section IV C) that a stable area can coexist with a growing pulse duration and eventual decay from relaxation.","fun_headline_variants_meta":{"raw":{"variants":["Fractional pulse areas stabilize graphene plasmon pulses","Graphene-area theorem breaks the pi rule for plasmon pulses","Nanoscale quantum memory with graphene plasmons at fractional pulse areas","Stable plasmon pulses at fractional areas enable quantum memory","Crossing the gap: graphene plasmons defy the pi-pulse rule"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000188,"raw_usage":{"total_tokens":1347,"prompt_tokens":974,"completion_tokens":373,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":291}},"tokens_in":590,"tokens_out":373,"duration_ms":3895,"temperature":1.0,"reasoning_tokens":291,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:46:39.426293+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A transmission experiment on a graphene double-sheet sample with a nanoscale gap would settle the claim: launch pulses of varying input area and measure the output area and temporal profile over several absorption lengths. If the fractional fixed points are real, output areas should lock to the predicted $\\Theta_2$ values, whereas growth of the pulse duration without area locking would falsify the propagation claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the graphene plasmonics platform properties, including tunable conductivity, strong field confinement, and low losses, on which the proposed structure is built."},{"cited_title":"Liu and et al., A review of graphene plasmons and its combination with metasurface, J","cited_arxiv_id":null,"evidence_quote":"Establishes low graphene plasmon losses below the Fermi energy, the premise that lets the authors neglect losses and quantize the SP field."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the free-space area theorem and the integer-$\\pi$ stability condition that the graphene modification generalizes to fractional values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the single-interface surface-plasmon area theorem and its unstable $2\\pi$ solution, which the two-sheet geometry extends and stabilizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the double-interface surface-polariton dispersion relation used to define the graphene SP modes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the coupled-mode transcendental equation for double charge-sheet structures that fixes the SP wave numbers."},{"cited_title":"Huttner and S","cited_arxiv_id":null,"evidence_quote":"Supplies the validity condition $\\delta t_s \\gg \\lambda/c$ for applying the area-theorem derivation to the short plasmon pulses."},{"cited_title":"Matloob, R","cited_arxiv_id":null,"evidence_quote":"Supports keeping only surface-plasmon modes in the Hamiltonian by showing atomic ensembles radiate into SP modes far more strongly than into free-space modes."}],"review_version":1}