{"id":"5a5120d0-61cb-440a-9378-e5b7aa09a1f1","arxiv_id":"2608.10075","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Power-law correlated incoherent pumping of a superradiant laser preserves coherent, ultra-narrow emission down to pump-range exponent alpha≈1 while reducing the required pump rate by up to a factor N^{1-alpha}.","lead":"Superradiant lasers that use a spatially correlated pump, with rates decaying as a power law of exponent alpha, keep their ultra-narrow linewidth for every alpha and remain fully coherent down to alpha around 1. This means fully coherent laser light does not require strictly local pumping, and the pump intensity (and the recoil heating it causes) can be reduced by up to a factor N^{1-alpha}.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The g2→1 at α≈1 crossover rests on unbenchmarked TWA data: SM S2 validates TWA only in the local-pump limit, at N=40 and w̃≲1/4, while the finite-α g2 minima lie outside that range with unknown bias sign.","rationale":"The reader's weakest assumption—that TWA's quantitative error at the relevant pump rates could invalidate the g2→1 extrapolation—is the same point I regard as load-bearing. I add that the SM S2 benchmark is only in the local-pump limit, so for finite α not even the sign of the bias is known; this is precisely the regime where the crossover is claimed. The direction of the local-pump bias (overestimation) would protect the qualitative coherence claim if it persisted, but that persistence is exactly what is unverified. The proposed small-N exact test directly controls the sign and N-dependence of the bias and the location of the minimum. Because the reader already assigned CONDITIONAL with medium confidence and this concern does not point to an internal contradiction or a reason to reject the physics, the verdict should remain CONDITIONAL; I therefore recommend UNCHANGED relative to the reader's verdict.","tokens_in":17757,"tokens_out":14030,"duration_ms":142067,"concrete_test":"For N=6,8,10,12 and α=0.5,0.7,0.9,1.0,1.1, compute the exact steady-state g2_min on the same 10-point w̃∈[0.1,1] grid using full Lindblad exact diagonalization or quantum trajectories (QuTiP/PIQS handles permutation symmetry only in the α=0 and α=∞ limits, so finite α requires the full Hilbert space for these small N), and run the TWA with identical parameters. Compare Δg2=g2_TWA−g2_exact at the minimum and the location w̃* of the minimum. If Δg2 remains positive and non-increasing with N and w̃* matches TWA, the extrapolated g2 limit is an upper bound and the α≈1 conclusion is supported; if Δg2 changes sign, grows with N, or shifts the minimum, the 1/N extrapolation in Fig. 3 is uncontrolled and the crossover boundary should be relabeled as an upper bound pending quantitative TWA diagnostics. Report 95% bootstrap confidence intervals for the N=10^3 and 10^4 g2_min points.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that fully coherent emission, g2→1 in the thermodynamic limit, survives at least down to α≈1. This is obtained by 1/N fits to g2_min values from TWA simulations at finite α. The method's only exact benchmark (SM S2) is in the local-pump limit, where TWA is quantitatively reliable for w̃≲1/4 and overestimates g2 and linewidth at stronger pumping. The relevant g2 minima, however, occur at pump rates comparable to the local-pump optimum w̃≈1/2 (SM S3), and more importantly the finite-α regime itself has no benchmark: permutation symmetry is broken, so the semiclassical noise treatment of the correlated pump is not controlled. If the TWA bias changes sign or grows with N once 0<α<∞, the monotone decrease of g2_min with α and its apparent approach to 1 in Fig. 3 could be a semiclassical artifact rather than a physical crossover. The paper's own caveat that TWA cannot capture intrinsically quantum effects without a well-defined mean-field order parameter applies precisely to the α≈1 coherent-emission regime. No error bars or convergence diagnostics are reported for the extrapolated limits, so the boundary at α≈1 is not quantitatively constrained.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a chain of N two-level atoms in the bad-cavity regime, subject to collective loss and an incoherent pump whose rates decay with interatomic distance as a power law, w_ij = w/(|i-j|+1)^alpha. Using a GPU-accelerated truncated Wigner approximation (TWA), the authors map the steady-state magnetization, intensity, second-order coherence g^(2), and linewidth over the rescaled pump rate w̃ = w_max^nu/(ΓN) and the exponent α, for systems up to N=10^4. Their central claims are: (i) the drive requirement for lasing is set by the largest eigenvalue of the pump matrix, reducing the required pump rate by N^{1−α} for α<1, by log N at α=1, and by ζ(α) for α>1; (ii) ultra-narrow linewidth persists for all α; and (iii) fully coherent emission, g^(2)→1 in the thermodynamic limit, survives at least down to α≈1, giving a continuous interpolation between collective and local pumping. The paper also advertises GPU-accelerated TWA as a practical tool for driven-dissipative spin systems beyond permutation symmetry.","tokens_in":18002,"tokens_out":6926,"duration_ms":68472,"significance":"If the results hold, the eigenvalue-scaling argument for the drive reduction is a clean and practically relevant result, and the identification of α≈1 as the boundary of full coherence would establish a new trade-off axis for superradiant lasers, with free-space emission (α≈1) singled out as a promising operating point. The manuscript has clear strengths: Eq. (3) is a parameter-free mathematical scaling statement; the collapse of the lasing region onto the rescaled pump rate is an empirical observation rather than a fit; the data and code are made available; and the GPU benchmark is detailed and reproducible. The paper is also commendably explicit about the known limitations of TWA, including its overestimate of g^(2) and linewidth at strong pumping. However, the central coherence crossover at α≈1 rests on TWA data in a regime that has not been benchmarked, and the finite-size extrapolations are quoted without uncertainties, so the quantitative boundary is not yet firmly established.","major_comments":[{"comment":"The central coherence result, g^(2)→1 for α≳0.7–0.9 in the thermodynamic limit, is obtained from TWA data for finite α, but the only quantitative benchmark (SM S2) is for the local-pump limit α=∞ at N=40 and for w̃≲1/4, where TWA \"reproduces the exact solution\". The g^(2)_min data in Fig. 3 are sampled over w̃∈[0.1,1], and the finite-α regime is not benchmarked at all, so the known overestimate of g^(2) in the local strong-pump limit does not by itself control the error in the finite-α region. Please benchmark TWA against an exact small-N Lindblad solution for several finite α (e.g., N=6–12), report the value of w̃ at which g^(2)_min occurs and the TWA error at that point, and discuss whether the overestimate direction is expected to persist for finite α.","section":"§Methods; SM S2; Figs. 2–3"},{"comment":"The 1/N extrapolations of g^(2)_min are quoted as numbers (1.15 at α=0.5, 1.04 at α=0.7, \"even closer\" at α≥0.9) without statistical uncertainties, fit residuals, or sensitivity analysis. Since the extrapolated values are close to the coherent limit 1 and the TWA bias is of comparable magnitude, the statement that full coherence survives \"at least down to α≈1\" is not quantitatively constrained. Add error bars from trajectory resampling, report the fit form and residuals, and test sensitivity to the largest-N points and to the range of N included.","section":"Fig. 3(b) and §Results"},{"comment":"The statement that correlated pumping \"parametrically suppress[es] the recoil heating\" is an inference from the reduced required pump rate; no recoil-heating term or temperature observable is modeled. Since the abstract motivates the work with recoil heating, either include a minimal model of the heating (e.g., a momentum-diffusion or temperature observable coupled to the pump rate) or restate the result as a reduction of the required drive intensity, with the heating reduction presented as a prediction rather than a demonstrated outcome.","section":"§Results, paragraph on recoil heating"},{"comment":"The ultra-narrow linewidth claim for all α rests on a spectral estimate whose low-frequency resolution is set by the finite evolution time τ=5/Γ, yet no frequency resolution or convergence check with τ is reported. Without this information, \"ultra-narrow\" is not quantitatively defined and the comparison across α is not controlled. Please report the resolution δω=2π/τ, show Δν as a function of τ for representative α, and state the smallest resolvable linewidth.","section":"Fig. 2(d) and §Methods"}],"minor_comments":[{"comment":"The main text states w_max^ν ∼ wζ(α) for α>1, but the footnote correctly notes that this limit requires α−1 ≫ (log N)^{−1}. Make this crossover condition explicit in the main text, since the α=1 log N regime and the α>1 constant regime meet at this scale.","section":"Eq. (3) and footnote 32"},{"comment":"The caption of Fig. 3 says the data \"place an upper bound near α≈1 on the crossover to coherent emission,\" while the abstract states that full coherence survives down to α≈1. Clarify whether α≈1 is an upper or a lower bound on the coherent-emission region and use consistent language.","section":"Fig. 3 caption and Abstract"},{"comment":"The benchmark against the exact solution is shown only at N=40; provide at least one additional exact comparison at a different N for the local-pump limit, or state explicitly that permutation-invariant solvers are limited in N, so that the N-dependence of the TWA error is documented.","section":"SM S2"},{"comment":"The timing benchmark fixes τ_r=5 and adjusts ΓN so that ΓN·Δt is constant; state whether this protocol affects the relative performance for small N, where the physical value of ΓN would be much smaller than in the large-N runs.","section":"End Matter, GPU benchmark"},{"comment":"The sampling of w̃ is described as ten equally spaced values between 1/10 and 1; specify whether the same ten values are used for all α and N, and whether g^(2)_min is selected from these points or obtained by interpolation.","section":"§Results and Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope, and the eigenvalue-scaling result is likely robust, but the main coherence claim depends on TWA data in an unbenchmarked regime. I see no grounds for rejection, provided the authors add small-N exact benchmarks for finite α, uncertainty estimates for the g^(2) extrapolations, and either model the recoil heating or soften the corresponding claim. If the benchmarks are not feasible, the abstract's claim about α≈1 should be explicitly downgraded to a TWA-based extrapolation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper introduces a one-parameter family of correlated incoherent pumps interpolating between the collective and local limits, and uses GPU-accelerated TWA to map the steady-state phase diagram up to N=10^4. The central new claim is that fully coherent emission (g^(2)→1) survives at least down to α≈1, while the pump rate needed for optimal lasing drops by log N at α=1 (and by N^(1−α) for α<1), matching the free-space far-field envelope. That combination—a practical drive reduction with preserved coherence—is genuinely new and addresses a real trade-off in superradiant lasers.\n\nWhat the paper does well: the model is minimal and motivated, the eigenvalue-scaling argument behind Eq. (3) is clean and parameter-free, and the data collapse of the lasing threshold onto the rescaled pump rate is convincing. The GPU acceleration is a real technical contribution, and the authors are honest about TWA's limitations—they benchmark it against the exact solution in the local-pump limit and state clearly where it overestimates g^(2) and linewidth. The code and data are released, which is a plus.\n\nThe soft spots are exactly where the stress-test lands. The quantitative claim that g^(2)→1 at α≈1 rests on TWA simulations in the finite-α regime, where permutation symmetry is broken and the method is not benchmarked. The only exact benchmark in SM S2 is for N=40 in the local-pump limit, and TWA is quantitatively reliable only for w̃≲1/4 there, while the g^(2) minima are sampled over w̃∈[0.1,1]. No error bars are given for the 1/N extrapolations, so the location of the crossover is not tightly constrained. The recoil-heating suppression is also inferred from the reduced pump rate, not modeled; that is a framing issue more than a flaw, but it should be labeled as an inference. None of this is fatal—the qualitative picture is plausible, and the main drive-reduction result does not depend on TWA at all, since it follows from the eigenvalue scaling. But the abstract and title overstate the certainty of the g^(2) crossover.\n\nThis paper deserves a serious referee. The model and the numerical phase diagram are worth engaging with, and the GPU approach is broadly useful. A referee should press for TWA-accuracy diagnostics in the finite-α regime, error bars on the extrapolated g^(2) limits, and either softer claims about the α≈1 boundary or additional benchmarks. I'd bring it to a reading group and cite it for the model and the drive-reduction scaling, even while remaining cautious about the coherence crossover.","headline":"A clean, well-executed numerical study of a new power-law correlated pump family that makes a plausible case for a coherent-lasing crossover near α≈1, though the key g^(2) extrapolation rests on TWA data outside the benchmarked range.","tokens_in":18574,"tokens_out":2647,"would_cite":true,"duration_ms":25507,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V80","81-08"],"pacs":["42.50.Ct","42.55.-f"],"model":"deepseek-v4-flash","headline":"Power-law correlated pumping relaxes the superradiant-laser drive requirement: the optimal pump rate drops by $N^{1-\\alpha}$ for $\\alpha<1$ and by $\\log N$ at $\\alpha=1$, while ultra-narrow emission persists for all $\\alpha$ and…","keywords":["superradiant laser","correlated pumping","power-law pump kernel","truncated Wigner approximation","GPU acceleration","second-order coherence","linewidth","open quantum systems"],"falsifier":"Compute $g^{(2)}_{\\min}$ at $\\alpha=0.9$ and $\\alpha=1$ for $N$ from 40 to a few hundred using an exact or non-semiclassical method, such as a tensor-network or symmetry-adapted solver, and compare the $1/N$ extrapolation: if the extrapolated value saturates above 1, or the minimum shifts in rescaled pump rate relative to the truncated Wigner scan, the claim that fully coherent emission survives at $\\alpha\\approx1$ is falsified. In experiment, a chain of atoms with a tunable power-law correlated repumping kernel measuring $g^{(2)}$ versus $N$ at $\\alpha=1$ would settle the same question.","tokens_in":17510,"feed_emoji":"⚛️","tokens_out":10448,"duration_ms":86407,"temperature":0.7,"pith_summary":"Superradiant lasers store optical coherence in the atomic medium rather than the cavity, but the incoherent pump that sustains inversion forces a trade-off: local pumping gives coherent light only at a pump rate that grows linearly with atom number $N$ and heats the medium through recoil, while fully collective pumping removes that scaling but caps the coherence at $g^{(2)}=6/5$. The paper argues that a spatially correlated pump with power-law rates $w_{ij}=w/(|i-j|+1)^\\alpha$ interpolates between these extremes and breaks the trade-off. Using GPU-accelerated truncated Wigner dynamics for up to $10^4$ atoms, it finds that the pump rate needed for optimal lasing is reduced by a factor $N^{1-\\alpha}$ for $\\alpha<1$, by $\\log N$ at $\\alpha=1$, and by $\\zeta(\\alpha)$ for $\\alpha>1$ relative to the local-pump requirement, while ultra-narrow superradiant emission persists for all $\\alpha$. Coherence improves as the pump becomes shorter-ranged, with $g^{(2)}\\to1$ surviving at least down to $\\alpha\\approx1$, so fully coherent light does not require local pumping, and $\\alpha=1$ coincides with the far-field envelope of free-space dissipative couplings.","feed_headline":"Power-law pump cuts laser drive by log N, keeps light coherent","feed_subtitle":"The pump rate needed for optimal lasing drops by log N while ultra-narrow emission persists for all α.","key_machinery":"The central object is the pump-rate matrix $w_{ij}=w/(|i-j|+1)^\\alpha$ on a chain of $N$ two-level atoms, whose exponent $\\alpha$ tunes from a rank-one collective pump ($\\alpha=0$) to a full-rank local pump ($\\alpha\\to\\infty$). The argument is carried by the eigenvalue spectrum of this matrix: its largest eigenvalue, proportional to the inverse Kac factor $K^{-1}(N,\\alpha)=\\sum_j w_{ij}$, sets the lasing-region boundary and produces the $N^{1-\\alpha}$, $\\log N$, and $\\zeta(\\alpha)$ reductions in required pump rate, while the breaking of permutation symmetry at finite $\\alpha$ opens the Hilbert-space sectors that allow $g^{(2)}\\to1$. The computational machinery is the truncated Wigner approximation with independent stochastic trajectories evolved in parallel on GPUs, which lets the authors evaluate steady-state magnetization, intensity, second-order coherence, and linewidth for up to $10^4$ atoms.","core_discovery":"The paper's central claim is that correlated pumping of power-law type directly relaxes the drive requirement of the superradiant laser. Relative to the standard local-pump value, the rate needed to reach optimal lasing is reduced by a factor $\\sim N^{1-\\alpha}$ for $\\alpha<1$, by $\\log N$ at $\\alpha=1$, and by $\\zeta(\\alpha)$ for $\\alpha>1$, parametrically suppressing the recoil heating that limits current implementations. This reduction does not compromise the defining feature of the SR laser: superradiant emission with an ultra-narrow linewidth persists for all values of $\\alpha$. What it costs is coherence, which improves as the pump becomes shorter ranged: fully coherent emission, $g^{(2)}\\to1$, survives at least down to $\\alpha\\approx1$, and finite-size extrapolation suggests it may extend below. The three properties decouple: the drive requirement tracks the largest eigenvalue of the pump matrix $w_{ij}$, coherence tracks how strongly the pump breaks permutation symmetry, and the ultra-narrow linewidth is inherited from collective loss with minor sensitivity to pump range.","pith_inferences":["If the $1/N$ extrapolation is trusted, the coherent regime may extend below $\\alpha\\approx1$; finite-size data in the paper suggest but do not prove this, so an exact or experimental measurement of $g^{(2)}$ at smaller $\\alpha$ would test whether a true threshold exists.","The eigenvalue-scaling argument implies a design rule the paper does not state explicitly: any pump kernel whose largest eigenvalue grows slower than $N$ would relax the drive requirement, so shaping the kernel's spectrum rather than its real-space power law could further suppress recoil heating.","The paper's kernel is a purely positive envelope; realistic Green-tensor kernels add oscillatory phases and coherent dipole-dipole exchange, and whether those preserve the coherence-drive trade-off is a testable question for waveguide or free-space implementations.","Because the truncated Wigner approximation overestimates $g^{(2)}$ and linewidth at strong pumping, the quantitative boundary of the coherent region may shift when evaluated with methods that retain operator noncommutativity; the qualitative decoupling of drive, coherence, and linewidth is the part most likely to survive."],"forward_implications":["Free-space dissipative couplings correspond to $\\alpha\\approx1$, so a repumping scheme built from free-space emission would already lower the required pump rate by a factor $\\sim\\log N$ relative to the local-pump SR laser, without further dissipation engineering.","For $\\alpha>1$ the required pump rate is reduced by the constant factor $\\zeta(\\alpha)$, so even short-ranged correlated pumping softens the drive requirement while the lasing region broadens as $\\alpha$ grows.","Fully coherent emission with $g^{(2)}\\to1$ does not require local pumping; a power-law pump with $\\alpha\\gtrsim1$ combines coherence with a sublinear or constant pump-rate scaling.","The three defining properties of the SR laser decouple: drive requirement is set by the largest eigenvalue of the pump matrix, coherence by the degree of permutation-symmetry breaking, and linewidth by the collective loss channel.","GPU-accelerated truncated Wigner dynamics makes full steady-state phase diagrams for thousands of spins practical and applies to arbitrary pump or decay matrices, so the same scan can be repeated for other correlated-dissipation kernels."],"supporting_citations":[{"why":"Supplies the minimal bad-cavity SR-laser model and the local-pump result that the optimal pump rate scales as $\\Gamma N$.","marker":"[11]"},{"why":"Provides the fully collective-pump limit where coherence is capped at $g^{(2)}\\to6/5$, the benchmark the correlated pump must improve on.","marker":"[14]"},{"why":"Supplies the truncated Wigner approximation used for all simulations, including the response-function method for linewidth.","marker":"[20]"},{"why":"Supports the eigenmode decomposition of the pump matrix and the physical relevance of correlated dissipative kernels in atomic chains.","marker":"[25]"},{"why":"Underlies the inverse-Kac-factor estimate for the largest eigenvalue of the power-law pump matrix that yields the $N^{1-\\alpha}$, $\\log N$, and $\\zeta(\\alpha)$ scalings.","marker":"[28]"},{"why":"Provides the standard SR-lasing behavior of $g^{(2)}$ and linewidth used for comparison.","marker":"[35]"},{"why":"Supplies the exact permutation-invariant solver used in the Supplementary benchmark at $N=40$.","marker":"[37]"},{"why":"Supplies the semi-implicit integration scheme used to evolve the GPU trajectories.","marker":"[46]"}],"fun_headline_variants":["Correlated pump slashes laser drive, keeps superradiant light","Power-law pump cuts drive by log N, preserves ultra-narrow emission","Lower lasing drive with correlated pump, coherence intact","Pump correlations relax drive, sustain superradiant linewidth","Short-range pump reduces drive, still coherent at α≈1"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's central quantitative conclusions rest on the semiclassical simulation method being reliable for large systems over the pump-rate range that locates the best coherence; the method is benchmarked against an exact solution only at $N=40$ and the paper states it overestimates $g^{(2)}$ and linewidth at stronger pumping.","fun_headline_variants_meta":{"raw":{"variants":["Correlated pump slashes laser drive, keeps superradiant light","Power-law pump cuts drive by log N, preserves ultra-narrow emission","Lower lasing drive with correlated pump, coherence intact","Pump correlations relax drive, sustain superradiant linewidth","Short-range pump reduces drive, still coherent at α≈1"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000436,"raw_usage":{"total_tokens":2286,"prompt_tokens":1080,"completion_tokens":1206,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":696,"completion_tokens_details":{"reasoning_tokens":1119}},"tokens_in":696,"tokens_out":1206,"duration_ms":9610,"temperature":1.0,"reasoning_tokens":1119,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:15:17.869576+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $g^{(2)}_{\\min}$ at $\\alpha=0.9$ and $\\alpha=1$ for $N$ from 40 to a few hundred using an exact or non-semiclassical method, such as a tensor-network or symmetry-adapted solver, and compare the $1/N$ extrapolation: if the extrapolated value saturates above 1, or the minimum shifts in rescaled pump rate relative to the truncated Wigner scan, the claim that fully coherent emission survives at $\\alpha\\approx1$ is falsified. In experiment, a chain of atoms with a tunable power-law correlated repumping kernel measuring $g^{(2)}$ versus $N$ at $\\alpha=1$ would settle the same question.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the fully collective-pump limit where coherence is capped at $g^{(2)}\\to6/5$, the benchmark the correlated pump must improve on."},{"cited_title":"Hosseinabadi, O","cited_arxiv_id":null,"evidence_quote":"Supplies the truncated Wigner approximation used for all simulations, including the response-function method for linewidth."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the eigenmode decomposition of the pump matrix and the physical relevance of correlated dissipative kernels in atomic chains."},{"cited_title":"Campa, T","cited_arxiv_id":null,"evidence_quote":"Underlies the inverse-Kac-factor estimate for the largest eigenvalue of the power-law pump matrix that yields the $N^{1-\\alpha}$, $\\log N$, and $\\zeta(\\alpha)$ scalings."},{"cited_title":"Meiser and M","cited_arxiv_id":null,"evidence_quote":"Provides the standard SR-lasing behavior of $g^{(2)}$ and linewidth used for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the semi-implicit integration scheme used to evolve the GPU trajectories."}],"review_version":1}