{"id":"97d4cf74-aa09-4a1e-9a0e-a6705669539c","arxiv_id":"2601.05989","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For two-level atoms in a lossy Lorentzian cavity, exact non-Markovian dynamics shows a critical spectral width that grows with atom number and a peak-intensity exponent falling from 2 toward 1.5.","lead":"The paper solves the non-Markovian collective emission of two atoms exactly and simulates up to 1000 atoms in a lossy cavity, finding three regimes: a Markovian burst, critical pulsed emission, and a reabsorbing 'superabsorption' phase. A generalist should care because it predicts that cavity memory is strengthened by atomic cooperativity and that the textbook N² superradiance peak degrades to N^1.5 in the perfect-cavity limit.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite-size extrapolation underpins the 1.5 exponent; the TC limit is exactly solvable and should be used to validate it.","rationale":"The paper is serious: the pseudomode mapping for a Lorentzian reservoir is standard, the weak-symmetry reduction is valid, and the N = 2 analytical solution is a useful contribution. The central concern is not internal inconsistency but the weight placed on finite-size extrapolation. I agree with the reader's weakest assumption. However, the TC limit being exactly solvable means the missing convergence check is readily available. For finite λ, a convergence scan with larger N and tighter temporal tolerances would settle whether ν_m → 1.5. Because the authors' own text already qualifies the λ_crit monotonicity, the abstract overstates one secondary claim, but the main physical picture—three regimes, non-Markovian scaling degradation—remains plausible. Thus no verdict change; the paper should remain CONDITIONAL pending the exact TC check and numerical convergence benchmarks.","tokens_in":21918,"tokens_out":11262,"duration_ms":128453,"concrete_test":"Diagonalize the Tavis-Cummings Hamiltonian in the M = N excitation sector (tridiagonal (N+1)×(N+1) matrix; Eq. 25 with λ = 0), for N = 10^2, 10^3, 10^4, 10^5. From the exact eigenpairs, evaluate I(t) = -ω0 d/dt⟨n_atom⟩(t), locate I_max, and compute the local exponent ν_m via Eq. (28). If ν_m at N = 10^5 is not within 0.02 of 1.5 (or shows a log-type drift), the abstract's 'approaching a subquadratic law' is not supported; if it is, the finite-size extrapolation is confirmed for the perfect-cavity endpoint.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing claim is the announced approach to a subquadratic law, ν ≈ 1.5, in the perfect-cavity limit (abstract; Sec. 4.2, Fig. 7). The evidence is the local log-log exponent ν_m (Eq. 28) computed from peak intensities up to N = 1000, with no error estimate on the N → ∞ extrapolation and no analytic remainder bound. The short-time 'looks like Tavis–Cummings' argument (Sec. 4.2) explains a trend but cannot exclude N^{1.5±ε} or logarithmic corrections, nor can it certify that ν_m at N = 1000 is within any tolerance of the asymptote. This is not merely a technical gap: the TC limit is exactly solvable. For fixed total excitation N, the Tavis-Cummings Hamiltonian is tridiagonal in the Dicke basis, so I_max^TC(N) can be computed to machine precision by diagonalizing an (N+1)×(N+1) matrix for N up to 10^5 or more. The paper does not do this, instead leaving the exponent as a numerical extrapolation. A secondary over-reach is the abstract's unconditional 'increases monotonically' for λ_crit, which Sec. 4.1 explicitly qualifies as 'at least as far as the numerical data shows' and is scanned only to N = 100; if monotonicity fails at larger N, the memory-enhancement-by-cooperativity conclusion is weakened. Both share the same root: asymptotic claims without a convergence certificate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the spontaneous emission of N two-level atoms coupled to a common Lorentzian reservoir, going beyond Markovian and mean-field treatments. It presents an analytical solution for N=2 from the microscopic Schrödinger equation, yielding an exact non-Markovian master equation with jump operators between Dicke states and a regime of eternal non-Markovianity. For larger N, the authors map the Lorentzian reservoir to a damped pseudomode and exploit a weak symmetry to integrate the resulting Tavis–Cummings-type master equation exactly up to N=10^3. They identify three dynamical regimes (Markovian superradiant burst, critical pulsed emission, non-Markovian reabsorption/spontaneous superabsorption), report that the critical spectral width λ_crit separating these regimes increases with N, and claim that the peak-intensity scaling degrades from quadratic toward a subquadratic law, approaching the Tavis–Cummings exponent 1.5 in the perfect-cavity limit. The maximum reabsorbed intensity is also found to scale superlinearly, which the authors call spontaneous superabsorption.","tokens_in":22301,"tokens_out":3455,"duration_ms":36989,"significance":"If the main claims hold, this is a valuable contribution to non-Markovian collective radiation. The N=2 exact solution is a rare analytical result beyond single-excitation problems, and the pseudomode approach with weak-symmetry reduction enables exact treatment of systems far larger than is typical for non-Markovian dynamics. The predicted degradation of superradiant scaling and the emergence of spontaneous superabsorption are falsifiable and could stimulate experimental work. The paper's strengths include a self-contained derivation for N=2, a standard and well-founded pseudomode mapping, and a clear numerical method with explicit complexity statements. However, two asymptotic claims—the 1.5 power law and the monotonicity of λ_crit—are supported only by finite-N extrapolations without convergence certificates; these claims are central to the abstract and conclusions.","major_comments":[{"comment":"The claim that the peak intensity approaches a subquadratic law with exponent ν ≈ 1.5 in the Tavis–Cummings limit rests entirely on the local logarithmic slopes ν_m computed for N up to 1000. No error bar, convergence certificate, or analytic bound is given for the N→∞ limit. Since the λ=0 Tavis–Cummings model is exactly solvable by diagonalizing an (N+1)×(N+1) matrix in the Dicke basis, the authors should verify the 1.5 exponent to much larger N (e.g., N=10^4–10^5) or provide an analytic estimate. Without this, the abstract's 'approaching a subquadratic law' is an extrapolation, not a demonstrated result.","section":"Sec. 4.2, Eq. (28), Fig. 7"},{"comment":"The abstract states unconditionally that λ_crit 'increases monotonically with the number of emitters', but Sec. 4.1 explicitly qualifies this as 'at least as far as the numerical data shows' and the scan only reaches N=100. A monotonicity claim that is load-bearing for the memory-enhancement-by-cooperativity conclusion should either be proven for all N (e.g., by an analytic argument or by extending the numerical scan) or be softened in the abstract to 'increases over the range studied'.","section":"Abstract and Sec. 4.1, Fig. 5 (right)"},{"comment":"The factorization Ξ_TL(t,t') = e^{-λt}ξ(t') is derived after 'assuming that this limit can be interchanged with time derivatives and integrals'. This interchange is a technical assumption that underpins the entire N=2 analytical solution. The manuscript should justify when this is valid for a Lorentzian reservoir (e.g., by dominated convergence or by directly verifying the resulting equations), or at least state explicitly the mathematical conditions. As written, the derivation is self-consistent but relies on an unproven regularity assumption.","section":"Appendix A, Eq. (A.32)"}],"minor_comments":[{"comment":"The definition of ν_m uses N_m and N_{m+1}, but the caption says each point is placed above N_m with N_{m+1} the next data point. For the final point N_m=1000, N_{m+1}=1001, so the quoted local exponent uses a very close pair; this should be stated explicitly to avoid overinterpreting the last data point.","section":"Eq. (28) / Fig. 7 caption"},{"comment":"The derivation of τ_R = 2/(Nγ_M) from the superoperator trace is terse. A brief explanation of how Tr(D) is computed (e.g., in the Dicke basis) would help readers verify the effective √N coupling scaling.","section":"Sec. 4.2, Eq. (29)"},{"comment":"There is a minor inconsistency in sign conventions between Eq. (A.44) and the subsequent Laplace transform expressions; checking and harmonizing the notation would improve readability.","section":"Appendix A, Eq. (A.44)"},{"comment":"The text uses '|mintI|' in the figure caption while the body uses '|min_t I(t)|'. Please unify notation.","section":"Sec. 4.3, Fig. 8"},{"comment":"The paper cites related work on non-Markovian superradiance and pseudomode methods, but does not mention available exact results for the Tavis–Cummings model's photon-emission statistics. A short comparison with known Tavis–Cummings exact solutions would contextualize the 1.5 exponent claim.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a quantum-optics journal and the core methodology is sound. The main reservation is that the headline asymptotic claims outrun the numerical evidence. The authors can likely address this with a modest addition: compute the exact Tavis–Cummings peak intensity for large N to certify the 1.5 exponent, and either prove or appropriately qualify the monotonicity of λ_crit. No concerns about novelty or citation behavior."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look, and worth a serious referee, but the front-page claims are a bit ahead of the evidence. The genuinely new things are the exact N=2 master equation (the derivation in Appendix A is careful and self-consistent) and the pseudomode-plus-weak-symmetry route that gets numerically exact dynamics up to 10^3 emitters. That is real progress: previous multi-emitter non-Markovian treatments were stuck at single-excitation or approximate levels. The three-regime picture (Markovian burst, critical pulsed emission, non-Markovian reabsorption) is plausible and well illustrated, and the \"eternal non-Markovianity\" finding for two atoms follows from the analytic rates. Credit where due: the math is serious, the method is reproducible in principle, and the paper earns its readership. Soft spots, in order. First, the claimed subquadratic scaling, max I ~ N^1.5 in the Tavis-Cummings limit, rests on a finite-N local exponent computed up to N=1000 with no error estimate. That may be right, but the TC limit is an (N+1)-dimensional tridiagonal problem; you can diagonalize it to machine precision at N=10^5 and just check. The short-time argument explains the trend but does not certify the asymptote. Second, the abstract says the critical width \"increases monotonically\" while the text, fairly, says \"at least as far as the numerical data shows\" and the scan stops at N=100. That overreach is small but real. Third, the closed forms for I1 and I2 are omitted \"for brevity\"; fine if a supplement or code comes with the arXiv posting, but as is, the N=2 result cannot be independently spot-checked without redoing the appendix. Also worth noting: despite the \"numerically exact\" claim, no convergence benchmarks or code are shipped, so the N=1000 results are currently assertions. None of these are load-bearing enough to sink the paper. The N=2 analytic solution is likely correct and useful; the pseudomode reduction is clean; the physics is interesting. Who this is for: quantum optics and open-systems people, and anyone working on collective effects in structured reservoirs. They will learn something, and they should be given the chance to referee it. My own verdict: accept conditional on (a) the TC diagonalization check for the 1.5 law, (b) a numerical convergence benchmark, and (c) a less absolute phrasing of the lambda_crit monotonicity. The work is honest and the core is solid; it just needs to be as careful in its asymptotic claims as it is in its algebra.","headline":"Solid exact two-atom result and a clean numerical method; the headline 1.5 exponent and monotonic lambda_crit are extrapolations that need tightening before the central scaling claim is banked.","tokens_in":22759,"tokens_out":666,"would_cite":true,"duration_ms":10094,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V80","81S22"],"pacs":["42.50.Nn","03.65.Yz"],"model":"deepseek-v4-flash","headline":"Exact non-Markovian dynamics of N atoms in a lossy cavity shows that cooperative emission switches from a superradiant N^2 burst to subquadratic scaling (~N^1.5) as cavity memory grows, with spontaneous reabsorption of the emitted field.","keywords":["superradiance","superabsorption","non-Markovian dynamics","Tavis-Cummings model","Lorentzian cavity","pseudomode method","collective emission","open quantum systems"],"falsifier":"Compute the local scaling exponent ν_m for N beyond 1000 (e.g., 2000, 5000) using the same exact method: if the exponent does not approach 1.5 as N grows, or rises back toward 2, the perfect-cavity subquadratic claim fails. Separately, compute λcrit for N > 100; if it does not increase monotonically, the claimed cooperative enhancement of memory is incorrect.","tokens_in":21801,"feed_emoji":"⚛️","tokens_out":3646,"duration_ms":42166,"temperature":0.7,"pith_summary":"This paper tracks exactly how N initially excited two-level atoms in a lossy cavity radiate when the standard Markovian and mean-field approximations are dropped. It claims that the dynamics splits into three regimes depending on the cavity spectral width: a Markovian superradiant burst, a critical pulsed regime, and a non-Markovian regime in which atoms reabsorb part of the emitted field—spontaneous superabsorption. The paper further claims that the critical width separating these regimes grows monotonically with N, so that adding more atoms can push a system into the memory-dominated regime. Finally, it claims that the iconic quadratic scaling of the superradiant peak degrades with system size, approaching a subquadratic law (about N^1.5) in the perfect-cavity limit. This matters because it suggests that superradiance's cooperative enhancement is self-limiting, and that non-Markovian memory can be harnessed for collective absorption without external driving.","feed_headline":"Peak cooperative emission drops from N^2 to N^1.5 as memory grows","feed_subtitle":"Exact dynamics of up to 1,000 emitters shows memory-driven reabsorption replaces the classic superradiant burst.","key_machinery":"The central object is the radiated intensity I(t) and its local scaling exponent ν_m, defined through log-ratios of peak intensities at successive atom numbers. The computational machinery is the pseudomode method: the Lorentzian reservoir is replaced by a single damped bosonic mode, reducing the problem to a Tavis-Cummings-type master equation. A weak symmetry—conservation of total excitation number—block-diagonalizes the density matrix, cutting the cost to O(N^3) and enabling numerically exact results up to 10^3 emitters. For N = 2, a closed-form solution of the coupled integro-differential equations in the thermodynamic limit yields the exact master equation with three jump operators; its","core_discovery":"For N two-level atoms coupled to a Lorentzian cavity, the exact radiated intensity I(t) = -ω0 d⟨n⟩/dt reveals three spectral-width regimes separated by a critical width λcrit: for λ > λcrit the emission is a single delayed superradiant burst; at λ = λcrit the emission is pulsed, halting and resuming at finite times; for λ < λcrit the intensity becomes negative at intervals, signaling reabsorption of previously emitted photons. The authors derive a complete analytical solution for N = 2, giving the exact non-Markovian master equation and showing that at least one canonical decay rate is negative at all times—eternal non-Markovianity. For larger N they use a numerically exact method up to N =","pith_inferences":["If the N^1.5 asymptotic holds, it places a practical ceiling on cavity-based superradiant sources that rely on quadratic enhancement, suggesting that high-finesse cavities may be better suited for energy storage or coherent reabsorption than for maximum radiated power.","The monotonic λcrit(N) trend implies a cooperative enhancement of environmental memory that could be tested experimentally by fixing λ and measuring the onset of intensity revivals as the atom number is increased.","The same pseudomode-plus-symmetry approach could be extended to other structured reservoirs—such as photonic band gaps or multi-mode cavities—where one might observe a similar critical width and a possibly different asymptotic exponent.","The local-exponent method used here could be applied directly to finite-N experimental intensity data to estimate the asymptotic scaling without needing extremely large ensembles."],"forward_implications":["For fixed cavity parameters, increasing the number of atoms can move a system from the Markovian superradiant regime into the non-Markovian reabsorption regime, since the critical spectral width grows with N.","The standard I_max ∝ N^2 superradiant scaling is not an asymptotic law for large N in a lossy cavity; it degrades toward roughly N^1.5 in the perfect-cavity limit.","Spontaneous superabsorption—reabsorption without external driving—emerges naturally in the non-Markovian regime and its peak magnitude scales superlinearly with N.","The exact N = 2 solution provides an analytic benchmark for non-Markovian multipartite open quantum systems, including the first exact demonstration of eternal non-Markovianity in a two-qubit model.","Because the critical spectral width depends on N, engineering Markovian emission by tuning cavity parameters alone is insufficient; the atom number must be accounted for."],"fun_headline_variants":["Superabsorption replaces superradiance as cavity memory strengthens","Exact dynamics shows collective emission shifts from burst to reabsorption","Non-Markovian memory turns superradiant bursts into light reabsorption","Three regimes in cooperative emission: burst, pulse, and reabsorption","Cavity memory degrades superradiant scaling from N^2 to subquadratic"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The claimed subquadratic scaling law rests on the assumption that local exponents computed at finite N (up to 1000) have already converged to their N→∞ values, and the monotonic growth of the critical width is inferred from numerical data only up to N = 100, with no rigorous error bound or proof of the infinite-N limit.","fun_headline_variants_meta":{"raw":{"variants":["Superabsorption replaces superradiance as cavity memory strengthens","Exact dynamics shows collective emission shifts from burst to reabsorption","Non-Markovian memory turns superradiant bursts into light reabsorption","Three regimes in cooperative emission: burst, pulse, and reabsorption","Cavity memory degrades superradiant scaling from N^2 to subquadratic"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000173,"raw_usage":{"total_tokens":1118,"prompt_tokens":746,"completion_tokens":372,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":490,"completion_tokens_details":{"reasoning_tokens":278}},"tokens_in":490,"tokens_out":372,"duration_ms":4276,"temperature":1.0,"reasoning_tokens":278,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T11:29:32.050863+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the local scaling exponent ν_m for N beyond 1000 (e.g., 2000, 5000) using the same exact method: if the exponent does not approach 1.5 as N grows, or rises back toward 2, the perfect-cavity subquadratic claim fails. Separately, compute λcrit for N > 100; if it does not increase monotonically, the claimed cooperative enhancement of memory is incorrect.","supporting_citations":[],"review_version":1}