{"id":"06baf79a-9c5d-4ea7-a5b6-0f7e74b6eb8f","arxiv_id":"2607.06888","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"In a waveguide-coupled atomic array under Bragg conditions, finite photon propagation delay causes atom-atom correlations to localize near the initially excited atom with a correlation length scaling as (γτ)^{-1/2}.","lead":"This paper shows that when atoms in a waveguide interact with a time delay (because photons take time to travel between them), correlations between atoms get trapped near the initially excited atom instead of spreading out. The correlation length shrinks as a power law with the delay time, offering a new way to control where quantum information goes in light-matter systems.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Power-law scaling derived from a Bloch ansatz on an infinite translationally-invariant lattice, but numerics use N=21 with open boundaries; finite-size and boundary effects on the exponent are not controlled.","rationale":"The reader correctly identifies the Bragg/perfect-coupling idealization as a weakness, and this is a legitimate concern for experimental relevance. However, the authors explicitly acknowledge this limitation in their Discussion section, and the central analytical claim (the power-law scaling) is a mathematical result derived under stated assumptions. The more load-bearing concern for the correctness of the central claim is whether the numerical verification actually tests the same regime that the analytical derivation addresses. The derivation assumes an infinite, translationally invariant lattice with a pure Bloch mode; the numerics use a finite, open-boundary system with a localized initial condition. The smallest-delay data points — which most strongly constrain the power-law exponent — are precisely where finite-size effects are most dangerous (ξ ≈ 10 sites in an N=21 array). If those points are contaminated, the agreement between the analytical -0.5 exponent and the numerical -0.48 could be partially coincidental. That said, the analytical derivation itself is clean and parameter-free, and the large-delay regime (where ξ is small and finite-size effects are negligible) should still follow the power law. So even if the small-delay points are contaminated, the core result likely holds in the regime where it is most physically interesting (large delays, strong localization). This is why the verdict remains CONDITIONAL rather than being downgraded — the concern is about the rigor of the numerical verification, not about the analytical result itself. The reader's focus on the Bragg condition is about physical applicability, which is a separate axis from correctness of the stated mathematical claim.","tokens_in":13694,"tokens_out":934,"duration_ms":454447,"concrete_test":"Recompute the correlation length ξ(γτ) for N=21, 41, 61, and 81 atoms at γτ = 0.01, 0.02, 0.05 (the small-delay regime where ξ is largest). If the extracted ξ changes by more than 10% between N=41 and N=81 at any of these points, the finite-size contamination of the power-law fit is confirmed. Additionally, verify that g(j) at γτ=0.01 has converged to zero before the array boundary (j=10 for N=21); if g(j) is still non-negligible at the boundary, the exponential fit is extracting a correlation length from a truncated distribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytical derivation of ξ ∝ (γτ)^{-1/2} (Eq. S22) proceeds by substituting a Bloch-wave ansatz ϕ_j = e^{ikjd} into Eq. S17, which assumes a translationally invariant system. The sum over m runs from -∞ to +∞ (Eq. S18), corresponding to an infinite lattice. However, the numerical simulations that confirm the scaling (Fig. 4, showing ξ ∼ (γτ)^{-0.48}) use N=21 atoms with open boundaries and a central initially excited atom. The paper states that 'the simulation ensures that the solution converges and that boundary effects are eliminated,' but provides no quantitative evidence for this claim. The concern is twofold: (1) For large delays, ξ becomes small (approaching 1 lattice site), so N=21 should be sufficient — but for small delays, ξ grows, and the correlation function g(j) may be truncated by the finite array before reaching the asymptotic regime. The log-log fit in Fig. 4 spans γτ from roughly 0.01 to 1.0; at γτ=0.01, the predicted ξ ≈ 10, which is comparable to the half-array size of 10 sites. This means the smallest-delay data points — which anchor the power-law fit — could be contaminated by finite-size effects, artificially flattening the slope. (2) The Bloch ansatz assumes translational invariance, but the initial condition (single central excited atom) and open boundaries break this symmetry. The steady state is not a single Bloch mode but a superposition; the evanescent component's localization length could differ from the single-mode prediction. The paper does not verify that the steady-state spatial profile matches the assumed Bloch form.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This manuscript studies atom-atom correlations in a finite array of atoms coupled to a waveguide under the Bragg condition, focusing on the regime where the photon propagation delay between adjacent atoms is non-negligible. Starting from a single centrally excited atom, the authors solve the delay-differential equations of motion (Eq. 3) numerically and find that the steady-state excitation localizes near the initially excited atom. The central quantitative result is that the spin-spin correlation length $g(j) = e^{-|j|/ξ}$ scales as $ξ ∝ (γτ)^{-1/2}$, which is derived analytically in the Supplemental Material (Eqs. S14–S22) via a Bloch-wave ansatz on an infinite lattice and confirmed numerically (Fig. 4, giving exponent $-0.48$). The paper also examines the interplay between delay-induced localization and positional disorder (Fig. 5).","tokens_in":14118,"tokens_out":1417,"duration_ms":141860,"significance":"The paper addresses a timely question in waveguide QED: how retardation effects modify correlation transport in collective atomic arrays. The analytical derivation of the power-law scaling $ξ ∝ (γτ)^{-1/2}$ is a genuine strength—it is parameter-free, follows from a standard Bloch ansatz into the exact delay equation, and the geometric series evaluation is exact. The numerical exponent $-0.48$ agrees well with the analytical $-1/2$. The study of disorder-induced localization and its interplay with delay is a natural extension. The results are falsifiable and the model is clearly specified. The main conceptual novelty is the identification of a steady-state (rather than transient) correlation localization mechanism in non-chiral wQED, which connects to broader localization phenomena.","major_comments":[{"comment":"Finite-size control for the power-law fit (Fig. 4 and Supplemental Eq. S22): The analytical derivation assumes an infinite, translationally invariant lattice (the sum in Eq. S18 runs from $-∞$ to $+∞$), but the numerics use $N=21$ with open boundaries and a central initially excited atom. The manuscript states that 'the simulation ensures that the solution converges and that boundary effects are eliminated,' but provides no quantitative evidence. At the smallest delays in the fit range ($γτ ≈ 0.01$), the predicted $ξ ≈ 10$ is comparable to the half-array size of 10 sites, meaning the correlation function $g(j)$ may be truncated before reaching the asymptotic regime. This could artificially flatten the fitted slope. The authors should either (a) demonstrate convergence by showing that the fitted exponent is stable when $N$ is increased (e.g., $N=31, 41$), or (b) restrict the fit range toγ","section":null},{"comment":"Steady-state existence and robustness: The true steady state with non-zero excitation retention depends critically on the exact Bragg condition ($ω_0 τ = 2π$) and perfect emitter-waveguide coupling (no loss to non-guided modes). The Discussion acknowledges that deviations lead to eventual excitation loss, but the main text and figures present the localization as a steady-state phenomenon without quantifying how sensitive the power-law scaling is to small deviations from the Bragg condition. Since the Bragg condition $ω_0 τ = 2π$ ties the delay to the atomic frequency, it would strengthen the paper to show (even briefly) how the correlation length behaves for a small detuning from the Bragg condition, or to state explicitly that the power-law applies only in the quasi-stationary window.","section":null},{"comment":"Disorder section (Fig. 5): The disorder study states that 'arbitrarily close to the Bragg condition, the atomic excitation always decays at infinitely long times; however, this decay is slow enough to allow us to compute the correlation lengths at long times as a quasi-stationary process.' However, no timescale or criterion for 'long times' is specified, and it is unclear whether the quasi-stationary $ξ$ depends on the measurement time. The authors should specify the time at which $ξ$ is extracted and demonstrate that it is in a well-defined quasi-stationary plateau, or at least state the sensitivity of the extracted $ξ$ to this choice.","section":null}],"minor_comments":[{"comment":"The disorder scaling exponents in Fig. 5 are described as $σ^{-0.75}, σ^{-0.5}, σ^{-0.3}, σ^{-0.2}$ 'with increasing τ,' but the specificτ values corresponding to each curve are not labeled in the figure or caption. Please add labels.","section":null},{"comment":"In the caption of Fig. 4, the inset is described as showing 'the exponential fit to the spin-spin correlation function for the point indicated by an arrow,' but no arrow is visible in the main panel. Please add the arrow or clarify.","section":null},{"comment":"The Markovian steady-state formula $|β_c(t→∞)|^2 = ((N-1)/N)^2$ is given in the main text; the Supplemental Material derives it (Eq. S12), but the main text does not clearly state that this is for the initially excited atom specifically. A brief clarification would help.","section":null},{"comment":"The phrase 'the localization effects stem from the interplay between the interaction delay timeτand the temporal coherence length of the spontaneously emitted photons, given by $γ^{-1}$' could be made more precise: the scaling $ξ ∝ (γτ)^{-1/2}$ is dimensionally $√(γτ)$ in the denominator, but the physical interpretation of the geometric mean is not elaborated.","section":null},{"comment":"Reference [16] is cited as 'Phys. Rev. A113, 013701 (2026)' — please verify the year and volume.","section":null},{"comment":"The notation $g^{(6)}$ in the inset of Fig. 4 is unclear; the main text defines $g(j)$ in Eq. (5) without a superscript. Please reconcile.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The reader's stress-test concern about finite-size effects on the power-law exponent is the most substantive issue. The concern is valid: at $γτ=0.01$, $ξ≈10$ is indeed comparable to the half-array size, and the paper's claim of converged boundaries is unsupported. However, this is addressable with additional numerics or a restricted fit range, so it does not rise to the level of a major revision. The analytical derivation itself is clean and the core result is likely robust. The paper is a solid contribution to the non-Markovian wQED literature."},"author_rebuttal":null,"desk_editor":{"model":"glm-5.2","letter":"The main thing to know: this paper derives a parameter-free power law, ξ ∝ (γτ)^{-1/2}, for the steady-state atom-atom correlation length in a Bragg-configured waveguide QED array with delayed interactions, and confirms it numerically (exponent -0.48). The derivation in the Supplemental Material is clean — Bloch-wave ansatz into the delay equation, exact geometric series, short-delay expansion, no fitted constants. This is a genuine new result for bidirectional (non-chiral) waveguide QED steady states, extending the transient correlation localization that Windt et al. showed for chiral systems. The authors earn credit for deriving the scaling analytically rather than just fitting it, and for noting that concurrence and coherence functions inherit the same scaling in the single-excitation regime. The disorder analysis, showing that delay and disorder independently localize correlations, is a nice addition even though it is purely numerical. The reader's assessment of soundness (7) and novelty (6) is about right. The conditional verdict is slightly too cautious — the analytical result stands on its own merits. Now the soft spots. The stress-test concern about finite-size effects is the real issue. The analytical derivation assumes an infinite translationally invariant lattice, but the numerics use N=21 with open boundaries. At the smallest delays (γτ ≈ 0.01), the predicted ξ ≈ 10, which is right at the half-array size. The paper asserts boundary effects are eliminated but shows no evidence — no convergence study, no comparison with larger N. The small-delay data points anchor the power-law fit, so if they are contaminated by finite-size effects, the exponent could shift. This does not sink the paper because the analytical derivation is independent of the numerics, but it weakens the claimed numerical confirmation. The second concern — that the steady state requires the Bragg condition plus perfect waveguide coupling — is real but explicitly acknowledged by the authors. They are honest that deviations introduce decay channels and the localization becomes transient. This limits experimental relevance but does not undermine the theoretical contribution. The disorder exponents lack analytical backing, which is a minor gap. This paper is for waveguide QED theorists and people working on non-Markovian quantum optics. It deserves a serious referee who should ask for a finite-size convergence check and a brief discussion of how robust the scaling is when the Bragg condition is imperfectly met. The core analytical result is solid and worth publishing.","headline":"Clean power-law scaling for correlation localization in delayed waveguide QED, but finite-size control of the numerics is thin and the steady state requires idealized conditions","tokens_in":14709,"tokens_out":577,"would_cite":true,"duration_ms":72701,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.50.Nn","42.50.Ct","03.67.-a"],"model":"glm-5.2","headline":"Photon delay traps correlations: ξ scales as (γτ)^{-1/2}","keywords":[],"falsifier":"If the Bragg condition is not exactly satisfied or if atoms couple to non-guided modes, the retained excitation would eventually decay and no true steady-state correlation length would exist, making the power-law scaling inapplicable to the asymptotic regime.","tokens_in":13905,"feed_emoji":"📏","tokens_out":1082,"duration_ms":162906,"temperature":0.7,"pith_summary":"This paper studies a one-dimensional array of atoms coupled to a waveguide, where adjacent atoms are spaced so that the photon propagation phase between them satisfies the Bragg condition (a mirror-like resonance). When one atom is initially excited, the excitation partially spreads through the array via photon-mediated interactions. In the standard Markovian limit—where photon travel time is treated as negligible—the excitation distributes evenly across all atoms. The authors show that when the finite photon propagation time τ between neighbors is retained (a non-Markovian, delay-differential regime), a fraction of the excitation instead remains trapped near the initially excited atom. The atom-atom correlation function decays exponentially with distance from the central atom, defining a correlation length ξ. The central analytical result is that this correlation length obeys a power law: ξ is proportional to (γτ)^{-1/2}, where γ is the single-atom decay rate into the waveguide. This scaling is derived by substituting an evanescent Bloch-wave ansatz into the delay equations of motion and expanding in the short-delay limit, yielding a localization length set by the geometric mean of the photon coherence time γ^{-1} and the delay time τ. The power law is confirmed numerically (ξ ~ (γτ)^{-0.48}). Positional disorder further shortens the correlation length, with delay and disorder acting as independent but combinable localization mechanisms.","feed_headline":"","feed_subtitle":"In a Bragg-spaced waveguide-QED array, finite photon travel time causes atom-atom correlations to localize near the source with a power-law-","key_machinery":"The machinery consists of delay-differential equations of motion for the atomic excitation amplitudes β_j(t), which include retarded arguments β_{j'}(t − τ_{j,j'}) weighted by Bragg-condition phase factors e^{-iω₀τ}. The correlation length is extracted from the spin-spin correlation function g(j) = ⟨σ₀^(z) σ_j^(z)⟩ − ⟨σ₀^(z)⟩⟨σ_j^(z)⟩, which decays exponentially as e^{-|j|/ξ}. The power-law scaling is derived by inserting a Bloch-wave ansatz with imaginary wavevector (evanescent profile) into the delay equations, evaluating the resulting geometric series in the dispersion relation, and expanding in the short-delay limit |η| ≪ 1 to obtain κd ≈ √(2γτ), giving ξ ∝ (κd)^{-1} ∝ (γτ)^{-1/2}.","core_discovery":"The central discovery is that delayed photon-mediated interactions in a Bragg-configured waveguide-QED array cause steady-state atom-atom correlations to localize exponentially around the initially excited atom, with a correlation length that scales as the inverse square root of the dimensionless delay parameter γτ. This is distinct from the Markovian regime, where correlations spread uniformly across the entire array. The scaling ξ ∝ (γτ)^{-1/2} emerges analytically from the structure of the delay-differential equations under the Bragg condition, connecting the photon coherence time to the spatial extent of correlation confinement.","pith_inferences":[],"forward_implications":["The power-law scaling ξ ∝ (γτ)^{-1/2} provides a tunable knob for controlling the spatial range of correlations and entanglement in waveguide-QED arrays by adjusting interatomic spacing (and thus τ) or by engineering the waveguide group velocity v.","Because the Markovian limit is never perfectly realized in extended arrays, the intrinsic delay-induced localization identified here sets a practical lower bound on correlation transport even in nominally Markovian systems, affecting quantum network design.","The interplay between delay-induced localization and disorder-induced Anderson-like localization suggests a two-parameter landscape for engineering localized quantum states, potentially relevant for decoherence-free subspaces and quantum memory protocols.","The resemblance to bound states in the continuum (BIC) in extended arrays opens a path toward photonic bound states in non-Markovian many-emitter systems, which remain experimentally and theoretically unexplored beyond few-atom configurations."],"fun_headline_variants":["Delayed photon interactions localize atom-atom correlations in waveguide QED","Photon delay confines correlations to source atom in Bragg waveguide array","Non-Markovian delay shrinks correlation length in waveguide-coupled atoms","Correlation length scales as inverse square root of photon delay in QED array","Finite photon travel time localizes excitations near source atom in waveguide"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The existence of a true steady state with non-zero retained excitation depends on the Bragg condition (the photon propagation phase between adjacent atoms being an exact multiple of 2π) combined with perfect coupling of atoms to the guided modes and no loss to unguided radiation. The authors acknowledge that deviations from either condition introduce additional decay channels that prevent a true steady state, reducing the localization to a transient quasi-stationary effect to","fun_headline_variants_meta":{"raw":{"variants":["Delayed photon interactions localize atom-atom correlations in waveguide QED","Photon delay confines correlations to source atom in Bragg waveguide array","Non-Markovian delay shrinks correlation length in waveguide-coupled atoms","Correlation length scales as inverse square root of photon delay in QED array","Finite photon travel time localizes excitations near source atom in waveguide"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":519,"prompt_tokens":422,"completion_tokens":97,"prompt_tokens_details":null},"tokens_in":422,"tokens_out":97,"duration_ms":25876,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T23:29:43.024645+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the Bragg condition is not exactly satisfied or if atoms couple to non-guided modes, the retained excitation would eventually decay and no true steady-state correlation length would exist, making the power-law scaling inapplicable to the asymptotic regime.","supporting_citations":[],"review_version":1}