{"id":"0d8cf0f6-f528-44e9-ae03-fe7fe52fecee","arxiv_id":"2505.24052","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a two-dimensional electron gas, collective superradiant decay of spin qubits occurs for emitter spacings up to sqrt(λ_F v_F/Δ), about a micrometer for typical NV centers.","lead":"The authors study a group of spin qubits releasing energy into the electron gas of a nearby metal, rather than into light. They find the qubits can emit collectively, with superradiance and subradiance happening under spacing rules different from those for photons, and that the emitted electron current should form a spiral wave.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Coherent-state superradiance: SM Eq. (29) and main Eq. (8) double the collective decay rate; exact master-equation result for N=2 is R=γ0+γa/2, not γ0+γa.","rationale":"I read the paper as a careful conditional theory of collective emission into a 2D electron gas. The central claims are the two length scales and the spiral current. The two-length-scale structure is robust within the stated clean, non-interacting, d=0 model. However, the derivation of the coherent-state superradiance rate contains a concrete internal inconsistency: SM Eq. (28) and the master-equation calculation give a collective coefficient of 1/4 per ordered pair, while SM Eq. (29) and main Eq. (8) use 1/2. This is not a matter of convention: for N=2 the published formula predicts R=γ0+γa, while the exact solution of Eq. (2) is R=γ0+γa/2. The error does not destroy the qualitative result, since both formulas produce the same parametric scale sqrt(λ_F v_F/Δ), but it changes the prefactor of the superradiant enhancement and the precise threshold. The reader's identified weakest assumption (disorder / finite mean free path) is a genuine limitation and is explicitly acknowledged in the Discussion; I do not treat it as invalidating the paper. The factor-of-two error is more decisive because it is internal and verifiable, and it should be corrected before final acceptance. Hence I recommend conditional acceptance: the scientific message is sound, but the quantitative formula and the abstract's 'requires' wording need revision.","tokens_in":21025,"tokens_out":38529,"duration_ms":381198,"concrete_test":"Derive R for two dipoles in the coherent state directly from Eq. (2) and its adjoint: ∂t⟨σ_z^j⟩ = -γ0 - (1/2)γa for j=1,2, giving R=γ0+γa/2. If this contradicts SM Eq. (29), which gives γ0+γa, the factor error is confirmed. Then propagate the corrected coefficient into Eq. (8): R≈Nγ0/2+(Nπn/2)∫rdrγ(r), and verify whether the superradiance threshold shifts from a<λ'_SR to a<sqrt(π/2)λ'_SR.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The concrete load-bearing issue is an internal factor-of-two error in the coherent-state superradiance rate. From the adjoint of Eq. (2), ∂t⟨σ_z^j⟩ = -∑_n γ_{jn}(⟨σ_+^nσ_-^j⟩+⟨σ_+^jσ_-^n⟩). For the product state (|+⟩+|−⟩)/√2, this yields R = -1/2∑_j∂t⟨σ_z^j⟩ = Nγ0/2 + (1/4)∑_{n≠m}γ_{nm} (ordered sum). SM Eq. (28) is consistent with this, but SM Eq. (29) and main-text Eq. (8) replace the 1/4 coefficient with 1/2, doubling the collective term. For two dipoles, the exact result is R=γ0+γa/2, whereas Eq. (29) gives γ0+γa. Propagating the corrected coefficient into Eq. (8) gives R≈Nγ0/2+(Nπn/2)∫rdrγ(r), rather than the published R≈Nγ0/2+Nπn∫rdrγ(r). The length scale λ'_SR=sqrt(λ_F v_F/Δ) survives, so the qualitative central claim is intact, but the abstract's 'requires' phrasing overstates what is a sufficient condition, and quantitative predictions based on Eq. (8) are off by a factor of 2. The separate and explicitly acknowledged disorder/mean-free-path limitation is real but does not undermine the conditional theoretical claim.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies collective dissipative emission from an ensemble of color-center qubits (modeled as magnetic dipoles) coupled to the magnetic-field noise of a two-dimensional electron gas. The central results are two characteristic length scales: subradiance requires emitter separations below the Fermi wavelength (or below the height d when d > λ_F), while superradiance for coherently initialized emitters requires separations below sqrt(λ_F v_F / Δ), with v_F the Fermi velocity and Δ the emitter frequency. The authors derive a Lindblad master equation with spatially nonlocal decay rates obtained from the equilibrium transverse-current noise spectrum, compute the rates for a magnetostatic coupling model, verify the spectrum numerically (Fig. 2), and analyze the initial-time superradiance indicators for both fully-excited and coherent-superposition initial states. They also compute the transient electron-current response to the decaying macrospin and predict a two-arm spiral current pattern.","tokens_in":21383,"tokens_out":22135,"duration_ms":186690,"significance":"If the quantitative statements are corrected, this is a worthwhile and novel extension of Dicke physics to fermionic baths. The paper is self-contained: the length scales emerge from Fermi-surface properties and the magnetostatic coupling matrix, with no fitted parameters, and the analytic noise approximations are checked numerically. The emitted-current spiral is a distinctive, in-principle falsifiable prediction. The authors also explicitly flag the main physical caveats (disorder, Fermi-liquid effects, zero-sound modes, and temperature), which strengthens the presentation. The main weaknesses are two algebraic errors in the superradiance section that affect displayed equations and quantitative thresholds; the qualitative length-scale claims survive these corrections.","major_comments":[{"comment":"The coherent-state superradiance rate contains a factor-of-two error. For the product state (|+> + |->)/sqrt(2), direct evaluation of the adjoint master equation gives R = N γ0/2 + (1/4) Σ_{n≠m} γ_nm. This is consistent with SM Eq. (28), since |ψ_+ ψ_-|^2 = 1/4, but SM Eq. (29) and the main-text Eq. (8) replace the 1/4 coefficient with 1/2. For N=2 the exact result is R = γ0 + γ_a/2, not γ0 + γ_a. Consequently Eq. (8) and the estimates following it (for example R ≈ N n γ0√A λ_F/2) overestimate the collective contribution by a factor of two. The length scale sqrt(λ_F v_F/Δ) survives, but the quantitative predictions and the abstract's 'requires' phrasing need to be corrected.","section":"SM Eq. (29) and main-text Eq. (8)"},{"comment":"The stated value λ_SR^2 ≈ λ_F^2 ln(4E_F/ℏΔ) does not follow from Eq. (7). Substituting C(Δ,q) = a/q for q ∈ (Δ/v_F, 2k_F) into Eq. (7) yields πλ_SR^2 = π ln(4E_F/ℏΔ)/(2k_F^2), i.e. λ_SR^2 = ln(4E_F/ℏΔ)/(2k_F^2). With the paper's convention λ_F ≈ 2/k_F this is λ_SR^2 ≈ (λ_F^2/8) ln(4E_F/ℏΔ), a factor of 8 smaller than the published expression (and a factor 8π^2 smaller if one uses λ_F = 2π/k_F). The same algebra applied to the top-hat model reproduces the SM result λ_SR^2 = 4q_max^{-2}, so the discrepancy is specific to the 1/q spectrum. This changes the density threshold for fully-excited superradiance by the same factor and weakens the claimed contrast with the subradiance threshold n >~ λ_F^{-2}.","section":"Main text, after Eq. (7)"}],"minor_comments":[{"comment":"There are several typographical errors: the abstract begins 'Correlated emission of light offer' (should be 'offers'), the author line contains 'Tserko vnyak', and the Fig. 3 caption has 'A video showing the evolution be found' (missing 'can').","section":"Abstract and title page"},{"comment":"The symbol γ0 is used in Eq. (8) before it is defined in the main text; define γ0 = γ(r=0) explicitly at first use.","section":"Main text, Eq. (8)"},{"comment":"The sentence 'so long as λ_F < √A < v_F/Δ' uses √A as a length, although A is introduced as an area; clarify that √A is the radius of the region.","section":"Superradiance section"},{"comment":"The word 'requires' overstates the status of the condition a < sqrt(λ_F v_F/Δ), because the paper itself notes that ∂tR(t=0)>0 is sufficient but not known to be necessary for a superradiant burst; rephrase to distinguish the sufficient condition from the open question about intermediate-time dynamics.","section":"Abstract and superradiance discussion"}],"recommendation":"major_revision","confidential_remarks":"The reader's favorable assessment is largely justified: the central qualitative claims are sound and the paper is a solid contribution to the field. However, the two errors in the Superradiance section are in central, displayed equations, and the second one changes the fully-excited superradiance threshold by an order-eight factor. These are local and fixable, so I recommend requiring a corrected derivation and a recheck of all prefactors before publication, rather than rejection. The self-citations to prior derivations of the coupling matrix element are appropriate and do not raise novelty concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First paper I know that takes correlated emission into a fermionic bath seriously, and the qualitative physics survives close reading. The key new idea is that a Fermi gas offers a continuum of wavevectors through the particle-hole spectrum, so cooperative effects are set by two length scales rather than one: subradiance wants spacing below λ_F, superradiance of coherence-initialized dipoles wants spacing below sqrt(λ_F v_F/Δ). For NV parameters that works out to roughly a micron, which is experimentally accessible. The derivation is transparent: a Lindblad master equation, a well-defined magnetostatic coupling, and an analytic noise spectrum checked numerically in Fig. 2. No fitted parameters. The spiral current transient is a genuinely new observable, and the static limit connects nicely to Landau diamagnetism. The paper also flags its own limitations—disorder and mean-free-path, and the open question about intermediate-time superradiance—which is honest.\n\nThe soft spot is real and specific. In the superradiance section, SM Eq. (28) gives for a product of equal superpositions R = 1/2∑_k γ_kk + 1/4∑_{n≠m} γ_nm. But SM Eq. (29) and the main-text Eq. (8) quote 1/2 on the off-diagonal sum, doubling the collective term. For two dipoles, the exact result is R = γ0 + γa/2, not γ0 + γa. With the correct coefficient, Eq. (8) becomes R ≈ Nγ0/2 + (Nπn/2)∫ r dr γ(r), off from the published expression by a factor of two. The length scale λ'_SR = sqrt(λ_F v_F/Δ) survives because it comes from where the integrand ∫ r dr γ(r) saturates, not from the prefactor. So the central claim is intact, but quantitative predictions using Eq. (8) are wrong by a factor of two, and the abstract's 'requires' wording overstates what is actually a sufficient condition for enhancement. This is a fixable error, not a conceptual collapse.\n\nWho it is for: quantum optics and condensed matter theorists working on collective dissipation, and NV experimentalists interested in coupling to two-dimensional conductors. It deserves a serious referee. My recommendation: send it to review, but ask the authors to correct the factor-of-two inconsistency between SM Eq. (28) and Eq. (29), and to soften the abstract's 'requires' to something like 'can occur when'. The paper will be a worthwhile contribution once that is done.","headline":"Novel and mostly sound theory of correlated emission into a Fermi gas, but a factor-of-two error in the coherent-state superradiance rate needs fixing.","tokens_in":21865,"tokens_out":8696,"would_cite":true,"duration_ms":68118,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Correlated emission into a two-dimensional electron gas is governed by two length scales: subradiance needs emitter spacing below the Fermi wavelength, while coherent superradiance works up to roughly a micron for typical nitrogen-vacancy…","keywords":["correlated emission","superradiance","subradiance","nitrogen-vacancy centers","two-dimensional electron gas","magnetic noise","transverse current fluctuations","quantum dissipation"],"falsifier":"Measure the nonlocal decay rate $\\gamma(r)$ between two color centers at known separation and height above a two-dimensional conductor: if $\\gamma(r)$ does not stay near the local rate for $r<\\lambda_F$ and then decay as $1/r$ out to $v_F/\\Delta$, the length-scale picture fails. A more direct check initializes both centers in a coherent superposition with spacing $\\lambda_F < a < \\sqrt{\\lambda_F v_F/\\Delta}$ and looks for a total decay rate exceeding the independent rate; no enhancement would refute the superradiance claim.","tokens_in":20885,"feed_emoji":"🌀","tokens_out":13437,"duration_ms":132668,"temperature":0.7,"pith_summary":"This paper asks when an ensemble of spin qubits, such as nitrogen-vacancy centers sitting above a metal film, emit energy cooperatively into the surrounding conduction-electron gas rather than independently. Electrons differ from photons because a fixed emitter frequency can deposit energy into electron-hole pairs over a continuum of wave vectors, so the noise spectrum has two natural length scales: the Fermi wavelength $\\lambda_F$ and the longer scale $v_F/\\Delta$ set by the emitter splitting $\\Delta$ and the Fermi velocity $v_F$. The authors find that subradiance requires emitters to be spaced closer than $\\lambda_F$, whereas superradiance, for emitters initialized with phase coherence, persists up to a spacing of roughly $\\sqrt{\\lambda_F v_F/\\Delta}$, about a micron for typical nitrogen-vacancy parameters. They further show that the energy leaving the spins takes the form of a two-armed spiral of transverse electron current moving at the Fermi velocity. If the picture is right, a two-dimensional conductor becomes a controllable dissipative environment for collective quantum behavior, with possible use in entangling spin qubits.","feed_headline":"Coherent spins superradiate into an electron gas at micron spacings","feed_subtitle":"Fermi wavelength governs subradiance; a micron-scale de Broglie length governs superradiance.","key_machinery":"The load-bearing object is the decay-rate matrix $\\gamma_{nm} = \\gamma(r_n - r_m)$, whose off-diagonal entries determine whether the ensemble decay is cooperative, together with its Fourier source, the magnetic-field noise spectrum $C^{-+}(\\omega, q)$. For a Fermi gas the noise is supported on the particle-hole continuum $q \\in (\\omega/v_F, 2k_F)$, and the relevant coupling is the magnetostatic matrix element for transverse current fluctuations, $V_{k,k+q}(0) = (\\hat{q}+i\\hat{d})(2\\pi e\\hbar/cm)e^{-qd}(\\hat{q}_\\perp \\cdot k)$. This element makes the noise spectrum approximately $\\propto 1/q$ in the continuum, which produces the constant, then $1/r$, then oscillatory behavior of $\\gamma(r)$. The superradiance length $\\lambda'_{SR} = \\sqrt{\\lambda_F v_F/\\Delta}$ comes from integrating that $1/r$ tail; the spiral current follows from the transverse-current response function, whose paramagnetic part has a traveling-wave factor $e^{i\\omega \\rho/v_F}$ and whose radial and azimuthal components fall off as $\\rho^{-3}$ and $\\rho^{-2}$, giving ballistic propagation at $v_F$.","core_discovery":"The central claim is that collective emission into a metallic electron gas is set by two length scales rather than one. In the Markovian weak-coupling limit, the ensemble is described by a Lindblad equation whose off-diagonal rates $\\gamma_{nm}$ are fixed by the magnetic-field noise spectrum $C^{-+}(\\Delta, q)$, and for a two-dimensional non-interacting Fermi gas this spectrum is nonzero for wave vectors $q \\in (\\Delta/v_F, 2k_F)$ and, for transverse current fluctuations, decays as $1/q$. Transforming that spectrum to real space gives a decay-rate matrix $\\gamma(r)$ that stays nearly constant for $r \\lesssim \\lambda_F$, decays as $1/r$ for $\\lambda_F \\lesssim r \\lesssim v_F/\\Delta$, and oscillates beyond $v_F/\\Delta$. Subradiance follows from the constant region: two emitters closer than $\\lambda_F$ form a dark singlet. Superradiance splits into two cases. For fully excited, incoherent emitters, the early-time acceleration of the decay is positive only when the density exceeds $\\lambda_{SR}^{-2} \\approx \\lambda_F^{-2}/\\ln(4E_F/\\hbar\\Delta)$; for emitters initialized in a coherent superposition, the $1/r$ tail of $\\gamma(r)$ accumulates over the whole region up to $v_F/\\Delta$, so enhancement survives up to a mean spacing $\\lambda'_{SR} = \\sqrt{\\lambda_F v_F/\\Delta}$, evaluated at height $d = 0$. Finally, the paper claims that the emitted energy is carried by a transverse-current wave with a two-armed spiral shape, whose radial periodicity is $v_F/\\Delta$ and whose packet width is set by the collective decay time.","pith_inferences":["If the Fermi surface is anisotropic, the noise cutoff $2k_F$ depends on the direction of $q$, so the thresholds for subradiance and superradiance should become direction-dependent; this could make correlated-decay rates a directional probe of Fermi-surface shape.","In a strongly interacting Fermi liquid, a transverse-current zero-sound mode would add a bosonic peak to the noise spectrum and might create a second superradiant channel with its own, likely longer, length scale; this is a testable Fermi-liquid calculation.","Because the paper shows that $\\partial_t R>0$ is sufficient but not necessary for a superradiant burst, exact small-$N$ Lindblad simulations could decide whether incoherent fully excited ensembles also superradiate at intermediate times at densities below $\\lambda_{SR}^{-2}/\\pi$.","The predicted spiral current wave could be looked for with scanning nitrogen-vacancy magnetometry above a two-dimensional conductor, provided the micron-scale periodicity can be resolved in time."],"forward_implications":["Two nitrogen-vacancy centers separated by less than the Fermi wavelength should form a nearly dark state, with collective decay suppressed below the single-emitter rate.","Coherently initialized ensembles should superradiate at mean spacings up to $\\lambda'_{SR} = \\sqrt{\\lambda_F v_F/\\Delta}$, roughly $1\\,\\mu$m for $\\Delta \\approx 1$ GHz, with an enhanced total rate scaling as $n\\lambda'^2_{SR}$ times the single-emitter rate.","The decay energy is carried away by transverse electron-current waves with a two-armed spiral structure; the spiral's radial period $v_F/\\Delta$ is a direct fingerprint of the emitter frequency.","In dirty conductors with mean free path $l < v_F/\\Delta$, the superradiance condition likely tightens to spacings below $l$ rather than below $\\lambda'_{SR}$.","At densities above $\\min(\\lambda_F^{-2}, d^{-2})$ subradiant states become typical; below that they still exist but are rare and localized on pairs closer than $\\lambda_F$."],"supporting_citations":[{"why":"Defines the cooperative-emission concepts of superradiance and subradiance that the paper generalizes from bosonic to fermionic environments.","marker":"[1]"},{"why":"Shows that dissipative coupling can generate Bell states between spin qubits, the application motivating correlated emission.","marker":"[12]"},{"why":"Reports a single NV center sensing Johnson noise and ballistic transport in a normal metal, establishing the experimental geometry and measured coupling rate.","marker":"[16]"},{"why":"Demonstrates nanoscale conductivity imaging with an NV center above a conductor, another instance of the modeled setup.","marker":"[17]"},{"why":"Supplies nitrogen-vacancy center properties and parameters used for the numerical estimates and length scales.","marker":"[28]"},{"why":"Establishes that an NV above a conductor couples predominantly to transverse current fluctuations, motivating the matrix element in Eq. (5).","marker":"[30]"},{"why":"Contains the derivations of the coupling matrix element, noise spectrum, superradiance conditions, and spiral current response used in the Letter.","marker":"[33]"},{"why":"Provides the early-time decay-acceleration criterion used to define the superradiance length scale.","marker":"[35]"}],"fun_headline_variants":["Two length scales decide emission into electron gas","Electron bath drives correlated spin emission waves","Fermi wavelength governs subradiance in electron seas","Color centers emit spiral currents in electron bath","Superradiance across micron spacings via electron noise"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions assume a clean, non-interacting two-dimensional electron gas whose magnetic noise comes only from transverse current fluctuations, with the superradiance result evaluated at emitters sitting at height $d = 0$; if disorder, longitudinal currents, or electron-electron interactions change the noise spectrum, the quoted length scales shift.","fun_headline_variants_meta":{"raw":{"variants":["Two length scales decide emission into electron gas","Electron bath drives correlated spin emission waves","Fermi wavelength governs subradiance in electron seas","Color centers emit spiral currents in electron bath","Superradiance across micron spacings via electron noise"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000227,"raw_usage":{"total_tokens":1534,"prompt_tokens":1070,"completion_tokens":464,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":686,"completion_tokens_details":{"reasoning_tokens":393}},"tokens_in":686,"tokens_out":464,"duration_ms":5294,"temperature":1.0,"reasoning_tokens":393,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:36:48.209068+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the nonlocal decay rate $\\gamma(r)$ between two color centers at known separation and height above a two-dimensional conductor: if $\\gamma(r)$ does not stay near the local rate for $r<\\lambda_F$ and then decay as $1/r$ out to $v_F/\\Delta$, the length-scale picture fails. A more direct check initializes both centers in a coherent superposition with spacing $\\lambda_F < a < \\sqrt{\\lambda_F v_F/\\Delta}$ and looks for a total decay rate exceeding the independent rate; no enhancement would refute the superradiance claim.","supporting_citations":[{"cited_title":"Correlated emission of electron-current waves","cited_arxiv_id":"2505.24052","evidence_quote":"Defines the cooperative-emission concepts of superradiance and subradiance that the paper generalizes from bosonic to fermionic environments."},{"cited_title":"Correlated decay sensitive to kF might then be able to probe the shape of the Fermi surface","cited_arxiv_id":null,"evidence_quote":"Establishes that an NV above a conductor couples predominantly to transverse current fluctuations, motivating the matrix element in Eq. (5)."},{"cited_title":"Zhang and Y","cited_arxiv_id":null,"evidence_quote":"Contains the derivations of the coupling matrix element, noise spectrum, superradiance conditions, and spiral current response used in the Letter."}],"review_version":1}