{"id":"47ea1ebd-3dc1-4d94-b8df-15910c22e774","arxiv_id":"2506.18224","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A decaying terahertz pulse excites graphene plasmon polaritons whose spatial decay is suppressed, directly demonstrating space-time duality and virtual gain.","lead":"This paper shows that a specially shaped terahertz pulse that fades in time can make graphene surface waves travel without fading in space. It demonstrates a trick called virtual gain in real time, which could help keep optical and electronic signals strong in nanoscale devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The time-dependent q_i that evidences spatial-decay compensation is likely inflated by finite-window fitting of causally truncated, self-normalized x-t slices; Eq. (1) alone cannot produce the observed q_i(t), so the transient fit must be validated against a horizon-aware model.","rationale":"The paper's central observation, a spatially sustained polariton profile inside the causal horizon for a temporally decaying THz pulse, is credible and consistent with known virtual-gain physics; the FDTD model and the direct x-t maps provide real supporting evidence. My concern is narrower and quantitative: the headline number q_i approaching zero is obtained by fitting finite, self-normalized profiles with a steady-state decay model that has no term for the causal horizon. Because the horizon expands with time, a fit that includes the unilluminated region beyond the horizon will return a q_i that falls from roughly Gamma/v_g to zero on the observed timescale, even if the physical profile inside the horizon is already flat. The reader's weakest assumption identified the related issue of broadband spectral content and constant Gamma, and also noted bias from finite self-normalized slices; I agree with that concern but would put the finite-window/horizon artifact first, since it affects the interpretation of the quantitative q_i(t) curve independently of the dispersion model. Eq. (1) predicts a time-independent q_i for fixed omega_i and Gamma, so the observed relaxation in Fig. 3d cannot be a direct steady-state test of Eq. (1); it must be explained by transient/causal effects or by an instantaneous-frequency generalization that the paper does not state. The proposed concrete test, applying the same extraction pipeline to synthetic maps with and without temporal decay and with and without horizon masking, would settle whether the reported q_i(t) reflects virtual gain or is a fitting artifact. The paper would also benefit from releasing the FDTD code and raw x-t data, but that is a reproducibility issue rather than a correctness flaw. Overall, the central claim likely survives, but the quantitative evidence as presented is less secure than the text implies; this reinforces the reader's CONDITIONAL verdict without moving it.","tokens_in":10344,"tokens_out":17134,"duration_ms":185798,"concrete_test":"Run the authors' FDTD solver on three synthetic x-t maps: (i) the measured multi-cycle THz waveform, (ii) a pure exp(-Gamma t)cos(omega_r t) source with Gamma=v_g*0.04 um^-1, (iii) a monochromatic source with no temporal decay. Apply the exact Supp. Sec. 3 q_i-extraction pipeline (same window, self-normalization, horizon masking) to each, and also to lossless versions. If the q_i(t) curves for case (iii) or the lossless sources fall from ~Gamma/v_g to zero as the horizon expands, the reported q_i(t) is a finite-window artifact. Then rerun case (i) with the full frequency-dependent conductivity versus a constant-Gamma model; if q_i(t) is unchanged, the broadband assumption is innocuous.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative support for the claim is the time-dependent imaginary wavevector q_i(t) extracted by fitting each x-t slice to E(x)=cos(q_r x) exp(-q_i x) (Fig. 3, Supp. Sec. 3), culminating in q_i falling below Gamma/v_g=0.04 um^-1 and approaching zero after about 2 ps. This extraction is not a clean test of Eq. (1). Eq. (1) is a steady-state relation for a single complex frequency with constant Gamma, yet the data are transient, causal slices of a broadband pulse. For a source switched on at t=0 the field is exactly zero outside the polaritonic horizon x>v_g t; the paper itself marks this horizon (cyan line, Fig. 2d) and stresses that the sustained region expands with time. If the fit window at early delays includes the unilluminated region beyond the horizon, the fit will absorb the sharp cutoff into a spuriously large q_i; as the horizon moves outward, q_i must relax toward zero even in a perfectly compensated or lossless medium. The self-normalization of each spatial profile to [-1,1] removes the global temporal envelope, so the fit cannot distinguish a spatially flat profile inside the horizon from a profile that is flat only because the overall amplitude is normalized away. The good agreement with the FDTD model does not resolve this if the same horizon and fitting pipeline are baked into the comparison. A second, related gap is that the scalar model and Eq. (1) assume a single constant Gamma across the 0.5-1.5 THz pulse bandwidth; graphene conductivity is frequency-dependent, so spectral reshaping during propagation can mimic reduced spatial decay without complex-frequency compensation. The paper's stated Eq. (1) cannot produce the observed time-dependence of q_i(t) by itself, so the time-dependence is doing the evidential work and needs a dedicated check.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports experiments and modeling on plasmon polaritons in bilayer graphene, claiming that excitation with a temporally decaying THz pulse (a complex-frequency excitation with omega_i/2pi ≈ -0.3 THz) suppresses the spatial decay of propagating polaritons, allowing q_i to fall below the intrinsic loss ratio Gamma/v_g = 0.04 μm^-1 and approach zero after roughly 2 ps. The central theoretical statement is the dispersion relation omega_i - v_g q_i = -Gamma (Eq. 1), which encodes a trade-off between temporal and spatial decay. Supporting evidence includes space-time THz near-field maps, spatial profile fits at successive time delays, a custom FDTD model, and comparison with an analytical model in a parameter-space diagram. The authors frame the work as a first direct visualization of polaritonic space-time duality and 'spatio-temporal virtual gain' without active gain or post-processing.","tokens_in":10614,"tokens_out":2667,"duration_ms":30415,"significance":"If the central claim holds, the result is significant: it would demonstrate in a passive system that temporal shaping of an excitation can compensate spatial losses, with implications for polaritonics, nanophotonics, and wave engineering more broadly. The manuscript is also valuable for combining time-domain THz near-field nanoscopy with a concrete damped-wave model, and for making a falsifiable prediction encoded in Eq. (1). The analytical relation is self-contained and the reported data show the expected qualitative trend of decreasing extracted q_i with delay. However, the quantitative evidence for full spatial-decay compensation rests on a fitting pipeline that is not, as presented, a clean test of Eq. (1). The main strengths are the direct space-time metrology and the explicit model-data comparison; the main weakness is that the extracted time-dependent q_i may be dominated by causal truncation and normalization artifacts rather than by complex-frequency compensation.","major_comments":[{"comment":"The extraction of q_i(t) by fitting each x-t slice to the steady-state monochromatic form E(x)=cos(q_r x)exp(-q_i x) is not a valid test of Eq. (1) on transient causal data: for a source switched on at finite time the field is exactly zero outside the polaritonic horizon x > v_g t, and fitting a decaying exponential over a window that includes this unilluminated region will bias q_i upward at early delays and force q_i to relax toward zero as the horizon advances, even for a constant-loss medium with no complex-frequency compensation. Please demonstrate, using the same FDTD model with a fixed Gamma and a real-frequency source switched on at t=0, that the extracted q_i(t) is not reproduced by the finite-window fitting procedure alone, or adopt a horizon-aware fitting model that accounts for the causal cutoff.","section":"Fig. 3(c,d) and Supplementary Section 3"},{"comment":"The spatial profiles in Fig. 2d are self-normalized to the range [-1,1], so the fits cannot distinguish a truly flat spatial envelope from one that is artificially flattened by normalization; the claim of complete suppression of spatial decay at 3 ps should be backed by an analysis of the unnormalized field amplitudes along constant-delay slices, with the horizon and the overall temporal envelope accounted for explicitly.","section":"Fig. 2d caption and 'Visualizing sustained polaritonic states'"},{"comment":"The model assumes a single complex frequency (0.8-0.3i) THz and a single constant loss rate Gamma, but the measured pulse spans roughly 0.5-1.5 THz and graphene's conductivity, and hence Gamma/v_g, is frequency-dependent; as a result Eq. (1) is not exact for the actual broadband pulse, and the observed q_i suppression could in part reflect transient spectral interference rather than a steady complex-frequency eigenmode. Please quantify this effect by repeating the FDTD simulation with the measured pulse waveform (or a band-limited complex-frequency source) and comparing the extracted q_i(t) with the prediction of Eq. (1) using the actual frequency-dependent Gamma.","section":"Eq. (1) and Fig. 2b"},{"comment":"The experimental points in Fig. 5 are compared with the analytical line for omega_i/2pi=-0.3 ps^-1, but the plotted experimental q_i values are themselves obtained from the same fitting procedure criticized above; if the fitting bias is present, the agreement with Eq. (1) in Fig. 5 is not independent evidence for the model and should be re-evaluated once the extraction is validated on synthetic horizon-aware data.","section":"Eq. (1) and Fig. 5"}],"minor_comments":[{"comment":"The caption references red and purple curves but the figure description is not fully self-contained; please clarify which panel corresponds to which color and ensure all color labels are defined.","section":"Fig. 1 caption"},{"comment":"The sentence 'temporally decaying excitation (omega_i<0) enables spatial amplification (q_i<0)' is only true when omega_i < -Gamma in a lossy medium; please state this condition explicitly to avoid ambiguity.","section":"Outlook section"},{"comment":"The black curves in Fig. 3(c,d) are said to be extracted from measured profiles, while the blue curves are FDTD predictions; the figure would benefit from stating in the caption whether the FDTD curves were processed through the same fitting and self-normalization pipeline as the experimental data.","section":"Fig. 3(d)"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely topic and the experimental dataset is impressive, but the quantitative central claim rests on the q_i(t) extraction, which is vulnerable to a causal-truncation bias that the current manuscript does not rule out. The requested synthetic-data validation and unnormalized-amplitude analysis are feasible within the scope of a revision and would substantially strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know: this is a real experiment, not a simulation paper, and the central effect—suppressed spatial decay of SPPs in bilayer graphene launched by a temporally decaying THz pulse—is directly visible in their x-t maps. The paper deserves a careful referee, but the quantitative q_i(t) extraction has a soft spot that needs addressing.\n\nThe genuinely new thing is the in-operando time-domain THz near-field implementation with direct space-time maps, not the mechanism itself. Virtual gain and complex-frequency excitation are established ideas, and Ref. 10 already showed polariton loss compensation. What this paper adds is a clean, all-optical demonstration in a tunable vdW platform, with data spanning carrier densities and temperatures. Eq. (1) is simple and correctly captures the trade-off between temporal and spatial decay. The inputs—pulse decay rate and intrinsic loss—come from separate measurements, so there is no obvious circularity. Credit also for comparing against an FDTD model, though the model's role is limited by the fitting pipeline, as I note below.\n\nThe main concern is the extraction of q_i(t). Each spatial slice is self-normalized and fit to a monochromatic decaying sinusoid, but the field is causal: it is exactly zero beyond the polaritonic horizon x > v_g t. If the fit window at early delays extends past the horizon, the fit absorbs the sharp cutoff into a spuriously large q_i, and q_i necessarily relaxes toward zero as the horizon expands. That could produce the observed time-dependence even in a perfectly compensated or lossless medium. The paper does not explicitly say whether the fits are restricted to the illuminated region; Supplementary Section 3 presumably has this detail, but it is not in front of me. This needs to be stated and justified. The self-normalization itself is not the issue—it does not flatten a decaying profile—but the finite-window truncation is.\n\nA secondary issue is the use of a single complex frequency (0.8 − 0.3i THz) for a pulse spanning 0.5–1.5 THz. Graphene's conductivity, and hence Gamma, varies across that band, so spectral reshaping could mimic some of the apparent decay suppression. The FDTD agreement helps, but only if the same fitting and horizon treatment are not baked into the comparison. Also, no code or raw data are deposited; the availability statement is the usual boilerplate.\n\nNeither issue sinks the paper. The late-time flat profile over ~20 µm at 3 ps is direct evidence that the effect is real; the artifacts mainly affect the transient, not the steady state. But the quantitative claim that q_i falls below Gamma/v_g and approaches zero rests on the fitting, so that part should be put on firmer ground. This is a solid experimental contribution for the graphene-plasmonics and THz-nanoscopy crowd, and for people working on non-Hermitian wave physics. Send it to peer review, but ask the authors to show fits confined to the causal region and to release the data. If the horizon artifact is ruled out, the paper is publishable essentially as is.","headline":"A credible time-domain realization of virtual gain in graphene plasmon polaritons, with a fitting caveat that should be checked before publication.","tokens_in":11352,"tokens_out":6437,"would_cite":true,"duration_ms":64120,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper demonstrates that a terahertz pulse with an exponentially decaying envelope can cancel the spatial decay of plasmon polaritons in bilayer graphene, achieving sustained propagation in a lossy passive medium.","keywords":["space-time duality","complex-frequency excitation","plasmon polaritons","bilayer graphene","terahertz near-field nanoscopy","virtual gain","polariton loss compensation","spatio-temporal wave engineering"],"falsifier":"Measure the polariton x-t map for the same graphene sample while scanning only the temporal decay rate of the excitation (for example $\\omega_i/2\\pi = -0.1,\\ -0.3,\\ -0.6$ THz) at fixed carrier density and temperature. If Eq. (1) holds, the extracted $q_i$ must fall on the single line $q_i = (\\omega_i+\\Gamma)/v_g$ with one fitted $\\Gamma$; a deviation that grows with pulse bandwidth, or a nonzero residual $q_i$ at $\\omega_i=-\\Gamma$, would show that the constant-$\\Gamma$ assumption fails.","tokens_in":10044,"feed_emoji":"⚡","tokens_out":8796,"duration_ms":84083,"temperature":0.7,"pith_summary":"The paper sets out to prove that the spatial decay of a polariton wave in a lossy material can be suppressed by giving its excitation pulse a matching temporal decay. Using time-domain terahertz near-field nanoscopy, the authors map plasmon polaritons in bilayer graphene with picosecond time resolution and show that the extracted imaginary wavenumber $q_i$ falls below the intrinsic loss rate $\\Gamma/v_g \\approx 0.04\\ \\mu m^{-1}$ and approaches zero within roughly two picoseconds. The governing identity is the complex-frequency dispersion relation $\\omega_i - v_g q_i = -\\Gamma$, which turns spatial decay and temporal decay into exchangeable quantities. If the claim is right, passive lossy media can sustain propagating waves without gain, simply by shaping the excitation envelope, and the recipe transfers to acoustic, photonic, electronic, and other polaritonic platforms.","feed_headline":"Decaying THz pulse cancels spatial decay of graphene polaritons","feed_subtitle":"A pulse that fades in time keeps graphene polaritons from fading in space, without any gain medium.","key_machinery":"The load-bearing object is the complex eigenfrequency of the excitation, encoded in Eq. (1): for a wave obeying a scalar damped wave equation with linear dispersion $\\omega_r = v_g q_r$ and constant amplitude loss rate $\\Gamma$, the imaginary parts of frequency and wavenumber obey $\\omega_i - v_g q_i = -\\Gamma$. This single identity converts the temporal envelope of the driving pulse into the spatial envelope of the propagating wave. In the experiment, the excitation is the multi-cycle THz pulse with complex frequency $(0.8-0.3i)$ THz; its imaginary part supplies the $\\omega_i$ that offsets the material loss, so the observed spatial decay rate $q_i$ is renormalized from $\\Gamma/v_g$ toward zero. The same relation organizes the paper's universal parameter map of normalized loss $\\Gamma/\\omega_r$ versus normalized spatial decay $q_i/q_r$, and it also fixes the boundary of the causally allowed region, the polaritonic horizon with slope $1/v_g$.","core_discovery":"On the paper's own terms, the discovery is the experimental realization of polaritonic space-time duality: a multi-cycle terahertz pulse whose envelope decays exponentially, with $(\\omega_r+i\\omega_i)/2\\pi \\approx (0.8 - 0.3i)$ THz, launches acoustic plasmon polaritons in bilayer graphene whose spatial attenuation is progressively compensated as time evolves. At delays beyond about 1.2 ps the measured $q_i$ falls below the ohmic decay rate $\\Gamma/v_g = 0.04\\ \\mu m^{-1}$, and by 2 ps it approaches zero, leaving a wave of constant oscillation amplitude over a roughly 20-micron span. The mechanism is captured by Eq. (1), which the paper derives from a damped-wave model with linear dispersion: $\\omega_i - v_g q_i = -\\Gamma$. Full compensation, $q_i=0$, occurs at $\\omega_i = -\\Gamma$, so a pulse that fades in time at the same rate the medium loses energy produces a propagation that does not fade in space. The sustained pattern is bounded by a causal 'polaritonic horizon' of slope $1/v_g$; beyond it, the excitation has not yet arrived. The authors call this spatio-temporal virtual gain and stress that it is implemented in-operando, with the pulse itself carrying the complex frequency, rather than synthesized afterwards from multiple real-frequency images.","pith_inferences":["Inference: because Eq. (1) is independent of the microscopic loss mechanism, the same compensation should be observable in phonon and exciton polaritons and in deliberately lossy acoustic or electronic transmission lines; running the same x-t measurement on such a platform would be a direct transfer test the paper does not perform.","Inference: if $\\Gamma$ is constant, a pulse with $\\omega_i < -\\Gamma$ should drive $q_i$ negative, i.e., apparent spatial amplification of a wave that still decays in time; measuring this over-compensation branch would test whether the model extrapolates beyond the compensation point.","Inference: the model's assumption of a single complex frequency could be probed by comparing pulses with identical bandwidth but different decay rates; a systematic scan of $\\omega_i$ at fixed carrier density and temperature would separate the complex-frequency effect from spectral-broadening artifacts.","Inference: the paper's framework suggests temporal analogs, such as using a spatially apodized source envelope to control the temporal decay of a response, which could act as a time lens or temporal grating in the same system."],"forward_implications":["Loss compensation in passive media: any wave system obeying Eq. (1) can sustain propagation without an active gain medium, by matching the temporal decay rate of the source to the intrinsic loss rate.","Universal parameter map: the normalized relation between $\\Gamma/\\omega_r$ and $q_i/q_r$ provides a design rule for choosing pulse envelopes in photonic, plasmonic, acoustic, and electronic waveguides.","Causal bound: the spatially sustained region grows at the group velocity from the source, so the achievable propagation length at a fixed observation time is bounded by $v_g t$; this sets a quantitative limit for virtual-gain schemes.","Real-time operation: because the complex frequency is carried by the physical pulse, the compensation works in-operando and avoids post-processing combinations of multiple real-frequency measurements.","Clean signature: successive spatial profiles within the horizon have constant amplitude while the whole pattern decays uniformly at rate $\\Gamma$ in time, distinguishing virtual gain from genuine material gain."],"supporting_citations":[{"why":"Introduced the virtual-gain concept in which complex-frequency excitation emulates gain in passive systems.","marker":"[1]"},{"why":"Established the theoretical basis that temporal decay of an excitation can emulate gain and circumvent loss bounds in scattering.","marker":"[3]"},{"why":"Demonstrated loss compensation in superlenses using synthesized complex-frequency waveforms, a direct precursor to the present pulse-based approach.","marker":"[5]"},{"why":"Provided the space-time metamaterial framework linking temporal decay to renormalized spatial momentum.","marker":"[9]"},{"why":"Showed that synthesized complex-frequency excitation can compensate losses in polariton propagation; this paper turns that idea into an in-operando demonstration.","marker":"[10]"},{"why":"Established the nano-terahertz space-time mapping methodology used to image polariton worldlines.","marker":"[18]"},{"why":"Supplies the terahertz scanning near-field optical microscopy technique underlying the measurements.","marker":"[23]"},{"why":"Quantified fundamental graphene plasmon losses, providing the intrinsic decay rate used for normalization and comparison.","marker":"[37]"}],"fun_headline_variants":["Time-fading THz pulse suppresses spatial decay of polaritons","Fading pulse yields non-fading graphene polaritons via duality","Polariton space-time duality: pulse decay cancels spatial loss","Graphene polaritons resist decay when excited by fading THz pulse","Spatio-temporal virtual gain via shaped THz pulses for polaritons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the broadband multi-cycle pulse behaves as a single monochromatic complex frequency, $(\\omega_r+i\\omega_i)/2\\pi\\approx(0.8-0.3i)$ THz, and that graphene's loss rate $\\Gamma$ is a constant across the pulse bandwidth; if $\\Gamma$ varies noticeably from 0.5 to 1.5 THz, Eq. (1) is not exact and the measured $q_i$ suppression could be a transient spectral-interference effect rather than a steady complex-frequency eigenmode.","fun_headline_variants_meta":{"raw":{"variants":["Time-fading THz pulse suppresses spatial decay of polaritons","Fading pulse yields non-fading graphene polaritons via duality","Polariton space-time duality: pulse decay cancels spatial loss","Graphene polaritons resist decay when excited by fading THz pulse","Spatio-temporal virtual gain via shaped THz pulses for polaritons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00075,"raw_usage":{"total_tokens":3374,"prompt_tokens":1015,"completion_tokens":2359,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":2264}},"tokens_in":631,"tokens_out":2359,"duration_ms":16040,"temperature":1.0,"reasoning_tokens":2264,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:55:40.814578+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the polariton x-t map for the same graphene sample while scanning only the temporal decay rate of the excitation (for example $\\omega_i/2\\pi = -0.1,\\ -0.3,\\ -0.6$ THz) at fixed carrier density and temperature. If Eq. (1) holds, the extracted $q_i$ must fall on the single line $q_i = (\\omega_i+\\Gamma)/v_g$ with one fitted $\\Gamma$; a deviation that grows with pulse bandwidth, or a nonzero residual $q_i$ at $\\omega_i=-\\Gamma$, would show that the constant-$\\Gamma$ assumption fails.","supporting_citations":[],"review_version":1}