{"id":"13f5732b-dd4b-4326-a267-b976c917d5f5","arxiv_id":"2509.07528","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A dissipative Rydberg lattice with several atoms per site is reported to switch between uniform and antiferromagnetic density-wave excitation order as the on-site interaction is tuned, with a claimed weak-universality scaling signature.","lead":"This paper studies an optical lattice where each site holds several Rydberg atoms, and reports that tuning the on-site interaction switches the laser-driven excitation pattern between a uniform state and an alternating, antiferromagnetic density wave. The interest: on-site interactions are tunable in experiments, so this may offer a practical knob for ordered states in dissipative Rydberg lattices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Two-sublattice ansatz is imposed before solving the lattice equations; the full mean-field model may have other steady states, so the central phase diagram is unverified.","rationale":"The reader's weakest assumption and my concern coincide: the two-sublattice reduction is the main load-bearing assumption. I considered the detuning-shift equivalence (Eq. 11) as an alternative; it is real and limits novelty, but it does not make the phase structure false, and the unequal-population section gives U a non-uniform role. The two-sublattice ansatz is more load-bearing because if it fails, the phase diagram, the oscillatory phase, and the scaling analysis built on theta all lose their foundation. The Pulser run is too small and uses a different dissipator to serve as a decisive check. The proposed full-lattice mean-field simulation directly tests the ansatz at the same level of approximation as the paper's own equations and would settle whether the phase classification is an artifact. No change to the CONDITIONAL verdict is needed; the condition is exactly this verification.","tokens_in":18263,"tokens_out":6803,"duration_ms":87423,"concrete_test":"Integrate the full site-resolved mean-field equations (4)-(5) on a 16x16 square lattice (or a 1D ring with L=12) without the two-sublattice ansatz, using random initial conditions and the parameters of Fig. 2 (Delta=1, Omega=3, V=8, N=4, for U in {-0.3, 0.6, 2.0} and around the left critical point). After long times, compute the spatial Fourier spectrum and local two-point correlations of {omega_j}. If the only stable long-time patterns are uniform or period-2, the ansatz is validated; if period-3/4, striped, or domain-wall states dominate, the central phase diagram is incomplete and theta is not a sufficient order parameter.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's phase classification is built on the two-sublattice reduction in Sec. II (after Eq. 5): all sites in each sublattice are assigned a common value omega_1 or omega_2, and the order parameter theta=|omega_1 - omega_2|/2 (Eq. 10) presupposes that any non-uniform order is exactly bipartite. This reduction is not derived from the full site-resolved mean-field equations (4)-(5); it is an ansatz. The full equations can in principle support period-3 or period-4 density waves, striped states, spirals, or local domain walls. None of these are captured by the two-mode ODEs (6)-(9), and theta would either miss them or report a misleading average. The only non-mean-field check is the 8-site 1D Pulser simulation in Sec. V, which uses a depolarizing noise model rather than the spontaneous-emission dissipator of Eq. 3, and whose order parameter is not defined via a full spatial pattern analysis. Thus the central claim that the steady state is 'either uniform or antiferromagnetic-like' is not established for the 2D/3D dissipative lattice to which mean-field theory is claimed to apply; the phase diagram in Fig. 2 may be an artifact of the truncation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a zero-hopping Bose-Hubbard-type chain of Rydberg 'superatoms' with two internal states, spontaneous emission, and nearest-neighbor Rydberg interactions, described by the Lindblad master equation. A mean-field decoupling reduces the lattice dynamics to two coupled ODEs for the two sublattice magnetizations (Eqs. 6-9), assuming a bipartite sublattice ansatz. The paper reports uniform/paramagnetic, antiferromagnetic/density-wave, and (mean-field-only) oscillatory phases as a function of the on-site interaction U, using the order parameter theta = |omega_1 - omega_2|/2 (Eq. 10). It analyzes fixed-point stability, semiclassical Monte Carlo correlation functions, and an 8-site Pulser simulation with depolarizing noise. It also proposes empirical scaling forms near the transition and at maximal antiferromagnetic order, with effective exponents beta = 0.025V + 0.45 and n = 0.017V + 0.585, interpreted as a signature of weak universality.","tokens_in":18453,"tokens_out":6334,"duration_ms":67649,"significance":"If the bipartite mean-field reduction is valid in 2D/3D, the paper provides a simple demonstration that density-wave order appears in a window of effective detuning in a multi-boson dissipative Rydberg lattice. The comparison with the 8-site Pulser simulation gives a genuine, though limited, cross-check for the density-wave phase, and the fixed-point analysis is a useful diagnostic. Credit is due for honestly stating that the oscillatory phase is only seen in mean field, and for making the numerical simulation reproducible via the open-source Pulser package. However, the paper's central framing overstates the novelty: the U-dependence is equivalent to a detuning shift through Eq. 11, the two-sublattice ansatz is imposed rather than derived, the quantum simulation uses depolarizing noise rather than the spontaneous-emission dissipator of Eq. 3, and the weak-universality claim rests on empirical fits with no theoretical underpinning. The core observation is plausible, but the manuscript requires substantial revision to support its advertised conclusions.","major_comments":[{"comment":"The bipartite ansatz, in which all sites in each sublattice share a single omega_1 or omega_2, is imposed before solving the lattice equations. The site-resolved mean-field equations (4)-(5) can in principle admit longer-period density waves, stripes, or domain-wall states. The order parameter theta = |omega_1 - omega_2|/2 presupposes exactly two uniform sublattices, so the phase diagram in Fig. 2 may be an artifact of this truncation. This is load-bearing for the central phase classification and should be tested, e.g., by a linear stability analysis of the full site-resolved equations or by allowing a period-3/period-4 ansatz.","section":"Sec. II, after Eq. (5); Eq. (10)"},{"comment":"The paper itself shows that the Hamiltonian of Eq. 1 is dynamically equivalent to Eq. 11 with Delta_m = Delta + U(N_j - 1). Consequently, tuning U is exactly equivalent to tuning the detuning in the single-species dissipative Rydberg lattice studied in Refs. [29,30]. The claim in Sec. VI that 'it is the on-site interaction U that serves as the driving parameter' and that this is 'fundamentally different' from earlier detuning-driven transitions is therefore not supported. The manuscript should either reframe the results as a multi-boson realization of the known detuning-driven density-wave transition, or identify a genuinely occupancy-dependent effect beyond the linear shift in Delta_m.","section":"Sec. IV, Eq. (11); Sec. VI"},{"comment":"The 'weak universality' claim is not established. The scaling forms are empirical fits whose normalization points (U_T, U_m, theta_m, u_1m) are taken from the same mean-field data, so the observed collapse is partly built in. The authors admit that beta and n are not conventional critical exponents and that Eq. 18 is inaccurate near (u1,theta)=(0,0). A signature of weak universality requires more than a linear dependence of a fitted exponent on V; it requires a theoretical justification or at least a demonstration of universal scaling ratios or correlation-length behavior. This claim should be substantially weakened or removed unless additional evidence is provided.","section":"Sec. VI, Eqs. (18)-(19) and Figs. 17-21"},{"comment":"The Pulser simulation uses a depolarizing noise model rather than the spontaneous-emission dissipator of Eq. 3, and it treats only eight sites in one dimension. This cannot quantitatively validate the mean-field phase diagram claimed for 2D/3D, nor can it address the oscillatory phase. The statement that the oscillatory region is absent because 'one cannot capture its signature with such a small number of sites' is not a substitute for a controlled check (e.g., a tensor-network simulation of the actual Lindblad master equation). The comparison in Fig. 7 and Fig. 8 is indicative but not a full validation of the model.","section":"Sec. V"},{"comment":"The semiclassical Monte Carlo procedure starts from mean-field omega values and then rearranges configurations using an ad hoc transition probability P taken from the single-atom result of Ref. [47]. This rearrangement is not derived from the master equation of the multi-boson system. The correlation functions obtained from this procedure therefore do not provide an independent confirmation of the mean-field phases; they are, at best, a consistency check. This limitation should be stated clearly, and the correlations should be compared with an exact small-system calculation where possible.","section":"Sec. IV A, Eq. (12)"}],"minor_comments":[{"comment":"The sign convention for U = U_GG - U_EG and its relation to Delta_m is confusing. Since U can be negative, the text should explicitly state how a negative U shifts the effective detuning, especially when discussing 'large absolute values' of U.","section":"Sec. II, Eq. (11)"},{"comment":"The caption lists values V=5,6,8,10,12,14 but does not label which color corresponds to which V in the text. Add an explicit legend or color list.","section":"Fig. 2 caption"},{"comment":"The quantity P is called a 'transition probability' but is not normalized over outcomes. Clarify whether P is a rate or a probability per unit time, and specify the parameter regime in which P < 1 is guaranteed.","section":"Sec. IV A, Eq. (12)"},{"comment":"Several typos and grammatical issues need correction, e.g., 'possess' for 'possess', 'principle quantum number' for 'principal quantum number', 'rols' for 'roles', and the double comma in the Fig. 13 caption.","section":"Throughout"},{"comment":"The text says 'we could not find traditional critical exponents here' but later uses the word 'criticality' in the section title. This inconsistency should be resolved; the scaling forms should be clearly distinguished from true critical exponents.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The manuscript occupies a borderline niche: the density-wave observation in both mean-field and the 8-site Pulser simulation is credible, but the novelty is undercut by Eq. 11, which the authors themselves state. I would ask the authors to compare their phase diagram directly with the known detuning-driven phase diagram of Refs. [29,30] and to either remove the weak-universality claim or support it with a proper scaling analysis. If the paper is revised in that direction, it could be a useful contribution; in its present form, the abstract and Sec. VI overclaim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read arXiv:2509.07528. The density-wave order is real: it shows up in the mean-field equations and in the 8-site Pulser simulation, and the fixed-point analysis is careful. The authors also do the honest thing of flagging that the oscillatory phase is mean-field only. That part is solid.\n\nThe trouble is the central novelty. For equal sublattice populations, Eq. 11 folds U into Delta_m = Delta + U(N-1), and U appears nowhere else. So the U-driven phase diagram is a reparametrization of the detuning-driven dissipative Rydberg lattice studied in Refs. [29,30]. The paper never compares to the single-atom limit or to the detuning-driven case, and the title's emphasis on occupancy and on-site interaction overstates what is actually new. The unequal-population case produces a staggered detuning, which is also a known knob from the staggered-detuning literature. So the 'roles of occupancy' is mostly an effective detuning story.\n\nThe weak-universality claim is the shakiest part. The exponents beta and n are fitted for each V, and the scaling form in Eq. 18 is admitted to be inaccurate in a window around u1=0, which is exactly the critical point. That means the 'scale-free' description does not actually cover the critical regime. Calling that a signature of weak universality is a stretch; at best it is an empirical collapse away from criticality. A referee should ask the authors to either make the claim quantitative about where the scaling holds or drop it.\n\nThe numerics have a mismatch: the Pulser simulation uses a depolarizing noise model, not the spontaneous-emission dissipator of Eq. 3, and the text does not discuss that choice. The AFM order appearing in both is still meaningful, but the simulation is not a check of the same open quantum model. Also, the correlation baseline correction is post hoc, and no code or data are provided. Minor, but standard for this kind of claim.\n\nOne structural worry from the stress-test: the two-sublattice ansatz is imposed before solving the lattice. For a density-wave order it is a natural guess, and the 8-site chain supports it, but the full site-resolved mean-field equations could in principle host longer-period or striped states. The paper doesn't rule those out. This makes the phase diagram potentially incomplete, not necessarily wrong.\n\nOverall, this is a competent paper with a sincere cross-support for the density-wave phase, but the claimed novelty and the universality signature are inflated. I would send it to review, because a careful referee can help the authors reframe the contribution and tighten or remove the scaling claim. I wouldn't cite it as a new mechanism in my own work, but it's a useful case study of how an exact reparametrization can masquerade as a new control knob.","headline":"The density-wave order is real and well supported, but the U-driven phase diagram is an exact reparametrization of the known detuning-driven case, and the weak-universality claim rests on scaling forms that miss the critical point.","tokens_in":19128,"tokens_out":3467,"would_cite":false,"duration_ms":34786,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["67.85.-d","42.50.Gy","3.67.-a","03.67.Bg","05.40.-a","05.70.Fh","05.90.+m"],"model":"deepseek-v4-flash","headline":"In a dissipative Rydberg lattice with several atoms per site, the on-site interaction strength U controls whether the steady state is uniform or an antiferromagnet-like density wave, with an oscillatory mean-field subphase and V-dependent s","keywords":["Rydberg lattice","dissipative open quantum system","mean-field approximation","density-wave order","antiferromagnetic order","on-site interaction","weak universality","Bose-Hubbard model"],"falsifier":"Take a larger (16-24 site) 1D or 2D dissipative Rydberg Bose-Hubbard simulation with fixed Omega, Delta, and V, sweep U through the predicted density-wave window, and compute the steady-state odd-even spatial correlation S_r and order parameter theta. If S_r never alternates in sign and theta stays flat near zero in the predicted U range, the density-wave claim fails. For the oscillatory phase, look for non-decaying irregular population oscillations with no nearby stable fixed point in a system larger than eight sites; if none appear at the predicted U, that phase is a mean-field artefact. In","tokens_in":17949,"feed_emoji":"⚛️","tokens_out":8413,"duration_ms":85839,"temperature":0.7,"pith_summary":"The paper studies a dissipative optical lattice in which several bosonic Rydberg atoms occupy each site (a zero-hopping Bose-Hubbard-like chain) and asks what steady-state pattern of Rydberg excitations forms under laser driving and spontaneous decay. Its central claim is that the on-site interaction strength U, equivalently the modified detuning Delta_m = Delta + U(N-1), controls whether the lattice is uniformly excited or develops an antiferromagnet-like staggered density wave of excitations. Mean-field dynamics also reveals an oscillatory subphase of the density-wave state, although the small one-dimensional quantum simulation does not reproduce it. The paper further claims that unequal atom numbers on alternating sublattices stabilize and widen the density-wave region, and that the order parameter near the transition obeys scaling forms whose effective exponents grow linearly with the Rydberg interaction V, a signature of weak universality. If correct, this gives an experimentally accessible knob for producing ordered, potentially entanglement-relevant states in a dissipative many-body Rydberg system.","feed_headline":"Tuning on-site repulsion flips a Rydberg lattice to a density wave","feed_subtitle":"On-site interaction picks the pattern of Rydberg excitations, opening a tunable route to antiferromagnetic-like order.","key_machinery":"The central machinery is the superatom plus two-sublattice mean-field ansatz: each multi-atom site is treated as a two-level system with at most one Rydberg excitation, and the lattice is assumed to split into two sublattices, each described by a single excitation imbalance omega_i. The order parameter theta = |omega_1 - omega_2|/2 then carries the entire phase classification. A companion fixed-point analysis of the mean-field equations identifies which branches are stable, explaining the uniform, density-wave, and oscillatory phases; the semi-classical Monte Carlo reconstruction uses the two-site excitation probability derived from the effective detuning Delta_m = Delta + U(N-1).","core_discovery":"By writing each multi-atom site as a superatom with at most one Rydberg excitation, and applying a two-sublattice mean-field ansatz, the authors derive four coupled equations for the excitation imbalances and coherences on the two sublattices. The order parameter theta = |omega_1 - omega_2|/2 separates the steady states: theta near zero for the uniform phase, positive for a staggered density-wave phase, within which a fluctuating, unstable oscillatory subclass appears when no stable fixed point exists. Fixed-point analysis and correlation functions support this picture: the uniform phase has nearly zero spatial correlations, while the density-wave phase shows alternating-sign correlations th","pith_inferences":["If the oscillatory phase is real, it should reappear as non-decaying staggered population oscillations in larger 1D or 2D simulations; the eight-site run may simply be too small to host it, so its absence there does not settle the question.","The fitted forms beta = 0.025V + 0.45 and n = 0.017V + 0.585 imply that the Rydberg interaction itself, not just U, could be used to dial the effective scaling behaviour; this is a testable prediction for higher-dimensional tensor-network or quantum-trajectory simulations.","Unequal sublattice populations create an alternating effective detuning Delta + U(N_j - 1), so selective loading effectively engineers a staggered-field Rydberg system, an avenue the paper notes but does not exploit for time-reversal or localization effects."],"forward_implications":["Sweeping the on-site interaction U (equivalently the modified detuning Delta_m = Delta + U(N-1)) switches the steady state between uniform and antiferromagnet-like density-wave order in the Rydberg excitation distribution.","The density-wave region is wider when the two sublattices are loaded with different atom numbers, and the peak order parameter is marginally larger, so selective loading is a practical way to stabilize the ordered state.","Within the density-wave region at the mean-field level there is an oscillatory subclass with no stable fixed points and weak alternating-sign correlations; the small one-dimensional quantum simulation does not capture it.","Near the uniform-to-density-wave transition and near the point of maximum order, the order parameter follows scale-free forms whose effective exponents beta and n increase linearly with the Rydberg interaction V, which the authors read as a signature of weak universality.","In a small one-dimensional lattice the order parameter peaks near Delta_m = 0 and the transition is a smooth crossover rather than a sharp one, with the ordered window widening as the dissipation rate gamma is lowered."],"supporting_citations":[{"why":"Introduces the antiferromagnetic density-wave phase in a driven nonequilibrium Rydberg sequence that this paper extends to multi-boson sites.","marker":"[29]"},{"why":"Reports antiferromagnetic long-range order in dissipative Rydberg lattices and supplies the density-wave classification used here.","marker":"[30]"},{"why":"Shows dissipative Rydberg pumping can prepare entangled and antiferromagnetic states, motivating why the ordered phase matters.","marker":"[31]"},{"why":"Provides the superatom description and the two-site excitation probability used in the semi-classical Monte Carlo correlation analysis.","marker":"[47]"},{"why":"Supplies the open-source simulation package used for the direct eight-site quantum numerical check.","marker":"[54]"},{"why":"Supplies the Ising-universality comparison against which the paper contrasts its own scaling behaviour.","marker":"[55]"},{"why":"Documents continuously varying effective exponents, the comparison class for the paper's V-linear scaling exponents.","marker":"[60]"},{"why":"Provides the hidden-superuniversality idea invoked to interpret the V-linear effective exponents.","marker":"[63]"},{"why":"Establishes Feshbach-resonance tuning as the experimental route for varying the on-site interaction U.","marker":"[71]"}],"fun_headline_variants":["On-site repulsion flips Rydberg lattice to density wave","Rydberg lattice order set by on-site interaction","Occupancy and interaction dictate Rydberg phases","Tuning on-site repulsion creates density-wave order","Dissipative Rydberg lattice: density wave from on-site repulsion"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The phase classification assumes the steady state always splits into two uniform sublattices, each sharing one excitation value; if the true state has longer-period waves, domain walls, or coexistence, or if mean-field is wrong in the relevant dimensions, the phase diagram, oscillatory region, and scaling fits lose their basis.","fun_headline_variants_meta":{"raw":{"variants":["On-site repulsion flips Rydberg lattice to density wave","Rydberg lattice order set by on-site interaction","Occupancy and interaction dictate Rydberg phases","Tuning on-site repulsion creates density-wave order","Dissipative Rydberg lattice: density wave from on-site repulsion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0006,"raw_usage":{"total_tokens":2619,"prompt_tokens":701,"completion_tokens":1918,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":1835}},"tokens_in":445,"tokens_out":1918,"duration_ms":14993,"temperature":1.0,"reasoning_tokens":1835,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T22:03:23.070840+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a larger (16-24 site) 1D or 2D dissipative Rydberg Bose-Hubbard simulation with fixed Omega, Delta, and V, sweep U through the predicted density-wave window, and compute the steady-state odd-even spatial correlation S_r and order parameter theta. If S_r never alternates in sign and theta stays flat near zero in the predicted U range, the density-wave claim fails. For the oscillatory phase, look for non-decaying irregular population oscillations with no nearby stable fixed point in a system larger than eight sites; if none appear at the predicted U, that phase is a mean-field artefact. In","supporting_citations":[{"cited_title":"Comparat and P","cited_arxiv_id":null,"evidence_quote":"Introduces the antiferromagnetic density-wave phase in a driven nonequilibrium Rydberg sequence that this paper extends to multi-boson sites."},{"cited_title":"Kübler, J","cited_arxiv_id":null,"evidence_quote":"Reports antiferromagnetic long-range order in dissipative Rydberg lattices and supplies the density-wave classification used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows dissipative Rydberg pumping can prepare entangled and antiferromagnetic states, motivating why the ordered phase matters."},{"cited_title":"Greiner, O","cited_arxiv_id":null,"evidence_quote":"Provides the superatom description and the two-site excitation probability used in the semi-classical Monte Carlo correlation analysis."},{"cited_title":"Samajdar, W","cited_arxiv_id":null,"evidence_quote":"Supplies the open-source simulation package used for the direct eight-site quantum numerical check."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Ising-universality comparison against which the paper contrasts its own scaling behaviour."},{"cited_title":"Numerical investigation of quantum phases and phase transitions in a two-leg ladder of Rydberg atoms","cited_arxiv_id":"2411.05494","evidence_quote":"Documents continuously varying effective exponents, the comparison class for the paper's V-linear scaling exponents."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the hidden-superuniversality idea invoked to interpret the V-linear effective exponents."},{"cited_title":"Sarkar, M","cited_arxiv_id":null,"evidence_quote":"Establishes Feshbach-resonance tuning as the experimental route for varying the on-site interaction U."}],"review_version":1}