{"id":"8c904195-fce9-43a7-a9f7-25298c05cf46","arxiv_id":"2506.19228","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A numerically optimized Rydberg atom chain can transfer single excitations at the quantum speed limit, with extrapolated entanglement distribution over roughly 54 atoms at 300 K and 90 atoms at 4 K.","lead":"This paper uses numerical simulations to show how off-resonant laser driving can turn a chain of Rydberg atoms into a spin channel that transfers a single excitation from one end to the other at the quantum speed limit. The authors identify control parameters that maximize the transfer probability and estimate that entanglement could be distributed over chains of about 54 atoms at room temperature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The scalar decay factor exp(-Gamma t_pi) is an uncontrolled approximation for the dissipative Rydberg dynamics; it neglects population recycling through the ground state and therefore biases p*_pi(L) and the extrapolated L_c.","rationale":"The paper's core methodology is sound: the effective XX mapping, the inversion procedure, and the unitary Rydberg simulations are internally consistent, and the paper honestly presents the extrapolation as a suggestion. The reader's conditional verdict correctly identifies the scalar decay approximation as the weakest link, and my analysis agrees: the factorization exp(-Gamma t_pi) is exact only in the idealized single-excitation, excitation-conserving limit. The native Rydberg Hamiltonian does not conserve excitation number, and the unitary simulations already show non-negligible m=0 and m>=2 populations, so the dissipative dynamics cannot be reduced to a global prefactor without a controlled error estimate. This issue is load-bearing because the headline chain lengths L_c = 54 and 90 are obtained by exponentiating a decay rate derived from this prefactor. A concrete Monte Carlo wavefunction calculation for L = 16 and L = 20 would settle whether the approximation introduces a systematic bias large enough to change L_c. My verdict remains CONDITIONAL, matching the reader's assessment, so no verdict change is needed.","tokens_in":12700,"tokens_out":7904,"duration_ms":90729,"concrete_test":"Run Monte Carlo wavefunction simulations of the dissipative Rydberg chain with jump operators sqrt(Gamma)|g_l><r_l| and the same Gamma(n,theta), at L = 16 and L = 20, using the optimal delta0*, n*, and spacings from Sec. III. Compare the resulting p_pi with the reported values obtained by multiplying unitary results by exp(-Gamma t_pi). If the relative difference exceeds about 10%, re-extract p*_pi(L) and the exponential decay rate; this directly determines whether L_c = 54/90 survives. If full MCWF is too costly at L = 20, run it at L = 16 and check the trend.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. III A multiplies unitary transport probabilities by exp(-Gamma(n,theta) t_pi) to account for radiative decay. This factorization is exact only if the dynamics is confined to the single-excitation subspace and the Hamiltonian conserves excitation number. Neither condition holds: the native Rydberg Hamiltonian contains Omega/2 sigma^x terms, and Fig. 3b shows substantial excursions through the m=0 and m>=2 subspaces even before decay. Under a Lindblad master equation with jump operators sqrt(Gamma)|g><r|, population that decays to the ground state can be re-excited by the same off-resonant drive, and population passing through m=2 decays at rate 2Gamma. The scalar factor applies a uniform decay to the entire final single-excitation amplitude, ignoring how long each amplitude component spent in the ground or doubly-excited sectors. The resulting error is not controlled by an explicit small parameter; for the parameters used (Omega ~ MHz, delta0 ~ tens of MHz, Gamma ~ kHz), the re-excitation probability over t_pi is of order (Omega/delta0)^2 Gamma t_pi, comparable to the correction itself. Because L_c is extracted from the exponential decay of p*_pi(L) in Sec. IV, this systematic error propagates directly into the headline claim of L_c = 54 at 300 K and 90 at 4 K. The extrapolation from L <= 20 to L = 54/90 also assumes the exponential trend persists and that a physical control solution exists at those lengths, but the decay-factor error is the more fundamental modeling gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies single-Rydberg-excitation transport along one-dimensional chains of 87Rb atoms. The authors solve the inverse problem of choosing atom positions and site-selective detunings so that the effective XX model generated by the off-resonantly driven Rydberg Hamiltonian reproduces the linear-spectrum perfect-transport couplings of Christandl et al. They then compare the transport probability p_pi for the Heisenberg XX model, the Rydberg model with nearest-neighbor interactions, and the full long-range Rydberg model, optimizing the global detuning delta_0 and principal quantum number n under fixed experimental constraints (Delta x_min = 3 micrometers, P_er = 4 W). Radiative decay is included by multiplying the no-decay transport probabilities by exp(-Gamma(n,theta) t_pi). They report optimal probabilities above 50% for chains up to L = 20 and extrapolate a maximum useful channel length L_c = 54 at 300 K and L_c = 90 at 4 K, corresponding to entanglement distribution over hundreds of microns.","tokens_in":13006,"tokens_out":6198,"duration_ms":65090,"significance":"The work is useful because it turns the known perfect-transport solution into a concrete, experimentally actionable parameter recipe for Rydberg arrays, and the comparison among XX, NN-Rydberg, and LRI-Rydberg models cleanly exposes the delta_0 trade-off between perturbation-theory breakdown and decay-induced loss. The use of ARC interaction data and explicit experimental constraints (trap spacing, laser power, temperature-dependent lifetimes) makes the predictions falsifiable and provides a concrete baseline for experiments. However, the headline length scale L_c rests on a post-hoc decay factor and on an extrapolation beyond the simulated chain lengths, so the quantitative claims need stronger support before they can be taken at face value.","major_comments":[{"comment":"The treatment of radiative decay by multiplying no-decay transport probabilities by exp(-Gamma(n,theta) t_pi) is not a controlled approximation for the parameters optimized here. The native Hamiltonian contains Omega/2 sigma^x terms, and Fig. 3b shows substantial excursions into the m=0 and m>=2 subspaces at the optimal delta_0; the dynamics is therefore not confined to the single-excitation subspace. Under a Lindblad master equation with jump operators sqrt(Gamma)|g><r|, population that decays to the ground state can be re-excited by the same off-resonant drive, and population passing through the doubly excited sector decays at a rate 2Gamma. A scalar multiplier applied at the final time ignores both effects, and the error is not controlled by an explicitly identified small parameter at the optimized detuning. Because L_c in Sec. IV is extracted from the exponential decay of p*_pi(L), this approximation propagates directly into the central quantitative claim. Please benchmark the scalar factor against full master-equation simulations for representative L, n, and delta_0, or provide a rigorous bound on the resulting error.","section":"Sec. III A"},{"comment":"The extrapolated values L_c = 54 (300 K) and L_c = 90 (4 K) are obtained from fits to data for L <= 20. This extrapolation assumes that the approximately exponential decay of p*_pi(L) continues without new degradation mechanisms and that a physical solution of the inversion problem exists at those chain lengths. The manuscript should either simulate larger chains, report the fit form with confidence intervals, or explicitly present L_c as an indicative extrapolation rather than a quantitative prediction. As written, the claim that entanglement can be distributed over 'more than fifty atoms' is not supported by direct simulation.","section":"Sec. IV, Fig. 6d"},{"comment":"The statement that 'as delta_0 approaches resonance, the effective interaction coefficients increase as J_max ~ delta_0^2' is inconsistent with Eq. (11). For fixed Omega, Eq. (11) gives |J_max| proportional to 1/(|delta_0|(1+|delta_0|/V_max)), which grows only as 1/|delta_0| near resonance and as 1/delta_0^2 far from resonance. The scaling narrative in Sec. III B should be corrected, since it is used to motivate the central trade-off between transport speed and excitation injection.","section":"Sec. III B, Eq. (11)"}],"minor_comments":[{"comment":"In the sentence 'when the state of the chain is initialized in an eigenstate of the N_j sigma_j^z operator', N_j is undefined; the operator is presumably the product over sites of sigma_j^z, and the notation should be corrected.","section":"Sec. II"},{"comment":"The colors used for 'approximation errors', 'radiative decay', and 'long-range interactions' are not identified in the caption or in the main text; the hierarchy diagram would be clearer if each source were explicitly labeled with its color.","section":"Fig. 2"},{"comment":"Please report the functional form and parameters of the exponential fit used to obtain L_c, together with confidence intervals, so that the extrapolation is reproducible and its uncertainty is visible.","section":"Sec. IV, Fig. 6d"},{"comment":"The approximate expression for lambda(n) in Eq. (13) should be checked for consistency: if lambda is J_NNN_max / J_NN_max, the denominator should be evaluated with the NN values and the numerator with the NNN values in a way that is explicit about which optimized detunings are used.","section":"Eqs. (12)-(13)"}],"recommendation":"major_revision","confidential_remarks":"The central modeling weakness is the scalar decay factor; if the authors can benchmark it against a master-equation treatment for representative parameters and temper the extrapolation claims, the paper would make a solid contribution. The delta_0 scaling error is easy to fix but should not be left in the published version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this as a numerical design study, not a new-physics paper. The work that holds up is the short-chain optimization; the headline numbers, L_c = 54 at 300 K and 90 at 4 K, are softer than they look.\n\nThe genuinely new pieces are the inversion that solves for atom positions and site-selective detunings realizing the Christandl linear-spectrum couplings in the native Rydberg model, the principal-quantum-number scan showing an optimum at n* = 71 set by next-nearest-neighbor interactions, and the extrapolated critical lengths. The short-chain results carry real weight: the full Rydberg simulation is a legitimate check on the effective XX mapping, the delta0 trade-off between perturbation-theory breakdown and decay is cleanly demonstrated, and the parameters are experimentally actionable.\n\nThe main soft spot is the decay treatment. Section III A multiplies the unitary transport probability by exp(-Gamma t_pi), which is exact only if the dynamics stays in the single-excitation subspace. Their own Fig. 3b shows substantial m=0 and m>=2 population, and the Hamiltonian contains sigma^x drive terms, so population that decays to the ground state can be re-excited and doubly-excited components decay faster. There is no explicit small parameter controlling the factorization error, and since L_c is read directly from the exponential decay of p_pi(L), the 54-atom and 90-atom claims inherit whatever bias the factor carries. This does not damage the L <= 20 results, where decay is a modest correction, but it makes the headline channel-length numbers conditional on a proper Lindblad simulation or experiment. The extrapolation beyond L = 20 also has no uncertainty bands, although the authors hedge with 'suggest,' which is fair.\n\nSecondary issues: no code or data are shipped, which hurts reproducibility for a purely numerical exercise; choosing delta0 and n on the same reported p_pi is a mild selection effect; and I would have liked one master-equation run at representative parameters just to quantify the error of the decay-factor ansatz.\n\nThe citation pattern is healthy: Yang et al. and Christandl are the right anchors, and the self-citations point to directly relevant work. This paper deserves a serious referee. If I were the editor I would send it out and ask the referee to require that the decay factorization be validated or removed before the L_c numbers go out as predictions.","headline":"Solid numerical design study for Rydberg perfect transport — the short-chain optimization and the n* = 71 trade-off hold up, but the L_c = 54/90 claims rest on a scalar decay factor and an unquantified extrapolation.","tokens_in":13564,"tokens_out":6968,"would_cite":true,"duration_ms":67103,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.67.Hk","03.67.Mn","32.80.Ee"],"model":"deepseek-v4-flash","headline":"A chain of neutral atoms excited to Rydberg states can be tuned to transport a single excitation from end to end at the quantum speed limit, distributing entanglement across chains of more than fifty atoms.","keywords":["Rydberg atom arrays","perfect state transfer","quantum speed limit","spin chains","entanglement distribution","Heisenberg XX model","Rydberg blockade","quantum transport"],"falsifier":"Measure the optimal transport probability $p_\\pi^*$ in a Rydberg chain of $L\\approx 30$ atoms using the reported control parameters; if the result lies significantly below the extrapolated curve $p_\\pi^*(L)\\sim\\exp(-t_\\pi^*(L)/\\tau(n^*))$ derived from $L\\le 20$, then the decay-factorization assumption or the extrapolation to $L_c=54$ and $90$ is wrong. A more direct test would be to build a chain of $L_c$ atoms at 300 K and check whether the two-spin concurrence is actually positive.","tokens_in":12433,"feed_emoji":"⚛️","tokens_out":8373,"duration_ms":77690,"temperature":0.7,"pith_summary":"This paper claims that a chain of neutral atoms excited to Rydberg states can be used as a coherent quantum channel that transfers a single excitation from one end to the other at the quantum speed limit, thereby distributing entanglement across the chain. The authors show that by tuning interatomic distances and site-selective detunings, the native Rydberg spin model maps onto the Heisenberg XX model with the linear-spectrum coupling profile known to give perfect state transfer. Numerical simulations then optimize the transport probability over the global detuning and the principal quantum number, balancing faster transport against the breakdown of perturbation theory and the growth of next-nearest-neighbor couplings. The central quantitative result is that optimal transport probabilities exceed 50%—the threshold for nonzero concurrence—for chains up to an extrapolated 54 atoms at 300 K and 90 atoms at 4 K, spanning hundreds of microns.","feed_headline":"Rydberg chains carry entanglement past 50 atoms","feed_subtitle":"Tuned spacings and detunings turn Rydberg chains into quantum channels up to 54 atoms at room temperature.","key_machinery":"The load-bearing mechanism is the linear-spectrum solution for perfect state transfer in the Heisenberg XX chain: couplings $J_{l,l+1}$ proportional to $\\sqrt{l(L-l)}$ with a constant on-site potential make a single excitation from site $l$ perfectly refocus at the mirror site $L+1-l$ at time $t_\\pi=\\pi\\bar{L}/(4J_{\\max})$, saturating the quantum speed limit. The paper's contribution is to solve the inversion problem that realizes this profile using Rydberg interactions: the effective flip-flop rate $J_{ll'}$ from second-order perturbation theory is matched to the linear-spectrum values by varying interatomic distances and atom-selective detunings, with the central spacing fixed at the minimum and the central detunings set to zero. The effective Hamiltonian is the Heisenberg XX model $\\hat{H}_{\\rm XX}=\\sum_l \\mu_l \\hat{n}_l + \\sum_{l<l'} J_{ll'}(\\hat{\\sigma}_l^+\\hat{\\sigma}_{l'}^- + \\hat{\\sigma}_l^-\\hat{\\sigma}_{l'}^+)$, whose single-excitation sector is equivalent to a tight-binding chain. The numerical optimization then sweeps the global detuning and principal quantum number to maximize $p_\\pi$.","core_discovery":"The paper's central discovery is that the perfect-transport condition of spin-chain quantum communication can be physically realized in Rydberg atom chains, not just assumed as a coupling profile. Solving the inversion problem for the set of control parameters $\\{\\Delta_l, V_{ll'}\\}$ under the mirror-symmetry constraint produces chains whose spacings and local detunings reproduce the linear-spectrum couplings $J_{l,l+1}=2J_{\\max}\\sqrt{l(L-l)}/\\bar{L}$ with uniform on-site potential. Under realistic driving powers and minimum spacing, the transport probability $p_\\pi$ is maximized at an optimal global detuning $\\delta_0^*$ that trades the speedup $t_\\pi\\sim\\delta_0^{-2}$ against excitation injection, and at an optimal principal quantum number $n^*$ that trades lifetime, interaction strength, and driving strength. For $L=16$ the optimum is $n^*=71$; the degradation at larger $n$ is attributed to next-nearest-neighbor couplings whose relative strength $\\lambda(n)$ diverges. With radiative decay included as $\\exp(-\\Gamma(n,\\theta)t_\\pi)$, the probability decays exponentially with chain length, giving the extrapolated critical lengths $L_c=54$ (300 K) and $L_c=90$ (4 K).","pith_inferences":["The same inversion strategy could in principle be applied to other effective exchange mechanisms the paper lists, such as resonant dipole-dipole interactions or Rydberg dressing; the paper does not compute the resulting $p_\\pi$, so this is an open quantitative question.","If the transport probability is used as a certification metric for realizing the Heisenberg XX model, then a single experimental measurement of $p_\\pi$ would simultaneously certify the effective Hamiltonian and the entanglement channel, a convenient diagnostic for Rydberg quantum simulators.","The degradation at large $n$ from next-nearest-neighbor couplings suggests that including NNN terms in the inversion problem—an inverse eigenvalue problem the paper mentions but does not solve—should push the optimal $n$ higher and may extend $L_c$ beyond the reported extrapolations.","The reported $L_c$ values are extrapolations of the exponential decay seen up to $L=20$; direct measurement at intermediate lengths would test whether the scaling persists or new decoherence mechanisms cut in."],"forward_implications":["For a fixed chain length, the optimal transport probability is a single-peaked function of the global detuning; approaching resonance speeds up transport but injects population into high-excitation manifolds, while moving off resonance slows transport and increases radiative loss.","The optimal principal quantum number is not the largest one; long-range interactions, specifically next-nearest-neighbor couplings, degrade the perfect refocusing beyond $n^*\\approx 71$ for $L=16$.","As the chain grows, the optimal $n^*$ drifts upward in a piecewise-continuous way while the optimal detuning stays roughly constant, so the transport time grows approximately linearly with $L$ and the optimal probability decays approximately exponentially.","At 300 K the channel supports entanglement distribution up to $L_c=54$ atoms; cooling to 4 K extends this to $L_c=90$ atoms, purely by suppressing blackbody-radiation-induced decay.","A transport probability above 50% certifies that the auxiliary and end spins share distillable entanglement, so the same measurement that demonstrates transport also demonstrates entanglement distribution."],"supporting_citations":[{"why":"introduced spin chains as quantum channels for communication, the conceptual basis for distributing entanglement across the Rydberg chain.","marker":"[3]"},{"why":"defines the perfect state transfer condition and the linear-spectrum coupling solution that the Rydberg chain is engineered to match.","marker":"[5]"},{"why":"establishes the quantum speed limit for perfect state transfer that the linear-spectrum transport time saturates.","marker":"[7]"},{"why":"supplies the second-order perturbative mapping from the Rydberg Hamiltonian to the effective XX model with flip-flop couplings.","marker":"[20]"},{"why":"provides the blackbody-radiation-dependent Rydberg decay rates used in the lifetime calculations.","marker":"[22]"},{"why":"shows how to include next-nearest-neighbor couplings in perfect transport, the benchmark for the observed degradation at large $n$.","marker":"[29]"},{"why":"supplies the atomic interaction coefficients and dipole matrix elements used to set the numerical parameters.","marker":"[31]"},{"why":"is used to compute the time evolution of the open quantum system and the transport probabilities.","marker":"[32]"},{"why":"defines the concurrence used to convert transport probability into an entanglement criterion.","marker":"[36]"},{"why":"establishes the reduction of Rydberg decay at cryogenic temperatures that underlies the 4 K estimate.","marker":"[37]"}],"fun_headline_variants":["Entanglement rides Rydberg chains at quantum speed limit","Quantum speed limit reached in Rydberg entanglement transport","Rydberg chain sends entanglement at speed limit over 50 atoms","Tuned Rydberg chains transport entanglement at speed limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative channel lengths rest on modeling radiative decay as a scalar factor $\\exp(-\\Gamma(n,\\theta)t_\\pi)$ applied after the no-decay evolution, which assumes decay removes population uniformly from the single-excitation subspace and that no population leaks back into it.","fun_headline_variants_meta":{"raw":{"variants":["Entanglement rides Rydberg chains at quantum speed limit","Quantum speed limit reached in Rydberg entanglement transport","Rydberg chain sends entanglement at speed limit over 50 atoms","Tuned Rydberg chains transport entanglement at speed limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00038,"raw_usage":{"total_tokens":2038,"prompt_tokens":988,"completion_tokens":1050,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":604,"completion_tokens_details":{"reasoning_tokens":982}},"tokens_in":604,"tokens_out":1050,"duration_ms":7686,"temperature":1.0,"reasoning_tokens":982,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:35:02.790829+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the optimal transport probability $p_\\pi^*$ in a Rydberg chain of $L\\approx 30$ atoms using the reported control parameters; if the result lies significantly below the extrapolated curve $p_\\pi^*(L)\\sim\\exp(-t_\\pi^*(L)/\\tau(n^*))$ derived from $L\\le 20$, then the decay-factorization assumption or the extrapolation to $L_c=54$ and $90$ is wrong. A more direct test would be to build a chain of $L_c$ atoms at 300 K and check whether the two-spin concurrence is actually positive.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the perfect state transfer condition and the linear-spectrum coupling solution that the Rydberg chain is engineered to match."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the second-order perturbative mapping from the Rydberg Hamiltonian to the effective XX model with flip-flop couplings."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the blackbody-radiation-dependent Rydberg decay rates used in the lifetime calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the atomic interaction coefficients and dipole matrix elements used to set the numerical parameters."},{"cited_title":"ˇSibali´ c, J","cited_arxiv_id":null,"evidence_quote":"is used to compute the time evolution of the open quantum system and the transport probabilities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the concurrence used to convert transport probability into an entanglement criterion."}],"review_version":2}