{"id":"5637e863-253c-4f67-b11e-d657750cb2e6","arxiv_id":"2412.11201","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A DMRG study of a two-leg Rydberg ladder with ax=2ay reveals Z_p^+ and Z_2p density-wave phases, two floating phases, and Ising, Potts, Ashkin-Teller, and chiral transitions between them.","lead":"This paper uses DMRG simulations to map the ground-state phase diagram of a two-leg ladder of Rydberg atoms in which the rung spacing is twice the leg spacing. It identifies several density-wave ordered phases, two families of floating phases, and the universality classes of the transitions between them, providing a map for future quantum simulation experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Floating phases FL and FL+ rest on boundary-dominated Friedel oscillations and structure-factor peaks, not on demonstrated algebraic decay; the claimed QLRO regions and their BKT/PT boundaries are therefore not yet established.","rationale":"The reader's conditional verdict already identifies the floating-phase characterization as a load-bearing weakness, and my stress-test converges on that point rather than on the interaction truncation. The 21-rung interaction cutoff is long enough that the omitted van der Waals tail is small on the plotted parameter range, so it is not the strongest threat to the central claim. The more serious issue is that Sec. IIIF asserts quasi-long-range order and BKT/PT transitions without the scaling analysis that would certify them. This does not overturn the reader's conditional verdict; it reinforces the need for additional evidence before the floating phases and their transition lines are treated as established.","tokens_in":18643,"tokens_out":2954,"duration_ms":34014,"concrete_test":"Compute C_m(i,r)=<m_i m_{i+r}> and C_B(i,r)=<B_i B_{i+r}> for L=289,397,541,613 at representative points (Delta/Omega,Rb/a)=(9.9,2.536) for FL and (9.9,2.552) for FL+, with truncation-error and bond-dimension convergence checks. Fit each to A cos(qr+phi)/r^eta and to A cos(qr+phi) exp(-r/xi)/sqrt(r). If eta is approximately constant across L and xi diverges with L, the QLRO claim is supported; if xi saturates or the exponent fit degrades with L, FL and FL+ must be relabeled and the BKT/PT boundary assignments should be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central phase diagram in Sec. IIIA includes two new floating phases, FL and FL+, and Sec. IIIF characterizes them at fixed L=289 via edge-enhanced profiles of <m_i> and <B_i>, fits to A sin(k i + phi), and incommensurate peaks of Sm(k) and SB(k). This evidence does not distinguish algebraic quasi-long-range order from an exponentially correlated incommensurate state with a moderately large correlation length: Appendix A explicitly shows that a structure-factor peak can coincide with the oscillation wave vector even when correlations decay exponentially, so the prominent peak heights in Fig. 13(c) alone do not imply QLRO. The Friedel oscillations shown in Figs. 13(d,e) are strongest near the edges and decay toward the center, which is the expected pattern for open-boundary critical systems, but similar edge-dominated profiles can appear in gapped phases near boundaries. No exponential-vs-algebraic fit of C_m(i,r) or C_B(i,r), no scaling of peak heights with L, and no Luttinger-parameter or central-charge estimate are provided. Since the transitions out of the floating phases into the disordered and crystalline phases are assigned to BKT and PT universality classes, and since the existence of an FL phase near Z6 is left as an open question, the classification of these regions is load-bearing for the claimed phase diagram. If FL/FL+ are actually short-range ordered incommensurate states, the diagram and its transition assignments change qualitatively.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript uses DMRG to map the ground-state phase diagram of a two-leg Rydberg ladder with lattice spacing ax=2ay, described by a blockade-model Hamiltonian with van der Waals interactions truncated to 21 consecutive rungs. The authors identify Z_p^+ density-wave orders that preserve top-bottom reflection symmetry, Z_{2p} orders that break it, Ising transitions between them, Potts and Ashkin-Teller CFT points on the commensurate lines, chiral transitions away from these points, and two floating phases FL and FL+ characterized by quasi-long-range incommensurate bond-order or density-difference correlations. Critical exponents ν and z are extracted from Binder-cumulant and energy-gap collapses, and the Kibble-Zurek exponent is predicted.","tokens_in":18985,"tokens_out":10378,"duration_ms":91164,"significance":"If the floating-phase classification is confirmed, the paper provides a detailed phase diagram for a two-leg Rydberg ladder with strong interactions and proposes concrete experimental signatures, such as a Kibble-Zurek exponent maximum at the CFT points. The numerical work is carefully described: DMRG truncation errors below 10^-10, explicit sweep convergence criteria, and system sizes up to L=613 are reported. The CFT identifications are supported by consistency of ν with Potts (0.838 vs 5/6) and Ashkin-Teller (0.8076) values and by z close to 1, as well as by the observed maximum of ν at the CFT points. The main weakness is that the quasi-long-range nature of FL and FL+ is not directly demonstrated.","major_comments":[{"comment":"The identification of FL and FL+ as quasi-long-range ordered phases rests on finite-size Friedel-oscillation profiles and structure-factor peak heights at fixed L=289. This evidence does not distinguish algebraic decay from exponential decay with a moderate correlation length; Appendix A explicitly shows that a structure-factor peak can coincide with the oscillation wavevector even for exponentially decaying correlations. The authors do not provide fits of C_m(i,r) or C_B(i,r) to algebraic versus exponential forms, do not show scaling of the peak heights with L, and do not extract a Luttinger parameter or central charge. Because the BKT and PT transition lines surrounding FL/FL+ in Fig. 1(a) are part of the central phase-diagram claim, the QLRO classification must be backed by direct correlation-function scaling or an equivalent finite-size analysis.","section":"Sec. IIIF, Fig. 13, Appendix A"},{"comment":"The central claim that the intersections of the Z_3^+ and Z_4^+ boundaries with the commensurate lines are Potts and Ashkin-Teller CFT points rests on the reported exponents, but no error bars or systematic-uncertainty estimates are given for ν, z, or the CFT-point coordinates (2.7901, 2.3749) and (3.3584, 3.4307). The CFT points are obtained from intersections of spline interpolations of seven or ten boundary points and extrapolated constant-k lines, all of which carry finite-size uncertainties. Without a sensitivity analysis (for example, varying the L ranges used in the Binder collapses, bootstrapping over system sizes, or testing the sensitivity to the 1/L correction parameter b in U'_4), the 0.6% discrepancy between ν=0.838 and 5/6 cannot be assessed.","section":"Secs. IIID and IIIE, Figs. 8 and 12"},{"comment":"The Hamiltonian retains van der Waals interactions only within any set of 21 consecutive rungs, and no convergence check with respect to this truncation is provided. Although the 1/r^6 tail beyond 21 rungs is small for the blockade radii studied, the phase boundaries and CFT points are determined by intersections of curves whose locations could shift slightly under a different truncation. A comparison with a longer truncation (for example, 31 rungs) or with the full interaction for at least one representative path in the phase diagram would establish that the reported diagram is not an artifact of this cutoff.","section":"Secs. IIA and IIF"}],"minor_comments":[{"comment":"The phrase 'within any set of consecutive twenty-one rungs' is ambiguous; please specify whether all pairs with rung-index difference up to 20 are included.","section":"Sec. IIA"},{"comment":"The notation Z_p for even p and Z_{2p} for the reflection-broken counterpart of Z_p^+ is potentially confusing because p is used with two meanings; consider a table defining each label and its order parameter.","section":"Fig. 1(a) and Sec. IIIA"},{"comment":"The statement that z 'continuously increases ... until it reaches the Lifshitz point with z=2' is not supported by data shown in Fig. 8; either add data near the Lifshitz point or soften the claim.","section":"Sec. IIID"},{"comment":"The sentence 'The prominent heights of the incommensurate peaks indicate the presence of floating phases' should be replaced by a quantitative comparison with disordered-phase peak heights at the same L and normalization, since Appendix A shows that peak position alone is not sufficient.","section":"Sec. IIIF"},{"comment":"The expression for the Fourier transform at η=2 appears to have a typo: the term 4(k-q)^2-4π|k-q|+4π^2/3 should be checked for missing parentheses or factors.","section":"Appendix A"},{"comment":"The energy gap ΔE in Eq. (7) is not defined explicitly; please state that it is the difference between the ground state and the first excited state in the DMRG spectrum.","section":"Sec. IIE"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid numerical DMRG study and likely within scope for a condensed-matter physics journal. The main required revision is a direct demonstration of quasi-long-range order in the floating phases; without it, the most novel portion of the phase diagram is not established. The comparison with Refs. [25] and [74] is appropriate, and I see no concerns about novelty or attribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid numerical mapping of a quasi-1D Rydberg ladder geometry that hasn't been done before. The authors retain interactions over 21 consecutive rungs, mirroring experimental conditions, and identify Z_p^+ and Z_2p phases with a two-step symmetry-breaking scenario. The strongest part is the transition analysis: Binder cumulant collapses give nu=0.838 at the Z_3+ boundary, matching the Potts value 5/6 within 0.6%, and nu=0.8076 at the Z_4+ boundary, inside the Ashkin-Teller window. The extracted z values sit near 1 at the CFT points and increase away from them. That internal consistency with known universality classes is real evidence that the phase diagram is essentially correct where it matters.\n\nWhat's genuinely new: the ax=2ay aspect ratio with a large blockade radius and full 21-rung interactions had not been mapped previously; earlier work used the complementary ay=2ax geometry, a small blockade radius, or a square ladder. The identification of Z_p^+ bond-order phases and two distinct floating-phase channels (bond-order vs density-difference) is a useful addition.\n\nThe soft spots are mostly concentrated in the floating-phase section. The FL and FL+ regions are characterized at fixed L=289 via edge-enhanced Friedel oscillations and structure-factor peaks. That evidence does not distinguish algebraic quasi-long-range order from an exponentially correlated incommensurate state with a moderately large correlation length. The paper does not provide correlation-function decay fits, finite-size scaling of peak heights, or a Luttinger parameter estimate, and the authors' own Appendix A shows that a structure-factor peak can sit at the oscillation wave vector even when correlations decay exponentially. The BKT/PT assignments for the floating-phase boundaries are borrowed from Ref. [25], a different ladder aspect ratio, rather than demonstrated for this model. This matters because the FL/FL+ regions are a central new claim. To the authors' credit, they are explicit about open questions, including the unresolved existence of the FL phase near the Z6 lobe.\n\nMinor issues: no error bars on exponents, a free 1/L correction in the modified Binder cumulant, and no code or data deposit. The 21-rung interaction truncation could shift boundaries, though it is a deliberate modeling choice tied to experiment.\n\nBottom line: the crystalline orders and the Potts/AT/chiral transitions look trustworthy, and this will be a useful reference for the Rydberg-ladder community. The floating-phase classification needs stronger evidence before I would treat it as established. I would send this to a serious referee, asking specifically for correlation-function decay analysis and L-scaling of the floating-phase signatures, and for more cautious wording on QLRO.","headline":"A credible DMRG phase diagram for the ax=2ay Rydberg ladder with solid Potts/AT/chiral transition analysis, but the floating-phase classification is suggestive rather than proven.","tokens_in":19490,"tokens_out":2256,"would_cite":true,"duration_ms":21710,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In a two-leg Rydberg ladder, increasing laser detuning first breaks translational symmetry into Z_p^+ bond-order waves, then breaks top-bottom reflection through an Ising transition into Z_2p phases, with two floating phases completing…","keywords":["Rydberg ladder","phase diagram","density wave order","bond-order wave","floating phase","chiral phase transition","Ising transition","DMRG"],"falsifier":"Recompute the phase diagram with interactions kept over 30, 40, or 60 consecutive rungs and check whether the Z_3^+ lobe boundary and the claimed CFT point at (Δ/Ω, Rb/a) = (2.7901, 2.3749) move outside the reported error bars, and independently fit the correlation functions C_B(i,r) and C_m(i,r) in the two floating phases at L = 613 to power-law decay; if the fitted exponents are not algebraic or the oscillation amplitudes vanish with system size, the FL and FL+ classification would not hold.","tokens_in":18482,"feed_emoji":"⚛️","tokens_out":4553,"duration_ms":41786,"temperature":0.7,"pith_summary":"The paper maps the ground-state phase diagram of a two-leg Rydberg ladder with lattice spacings ax = 2ay using the density matrix renormalization group. Its central claim is that detuning drives a two-step symmetry-breaking sequence: the disordered phase first develops a rung bond-order wave that breaks translational Z_p symmetry while preserving top-bottom reflection, and then a further increase in detuning breaks reflection symmetry through an Ising transition, producing staggered Z_2p density waves. The paper also identifies two distinct floating phases, one carrying quasi-long-range bond order and one carrying quasi-long-range density-difference order, and locates conformal points of the three-state Potts and Ashkin-Teller universality classes on the commensurate lines. A sympathetic reader cares because the result shows how lattice geometry controls which symmetry breaks first in programmable Rydberg arrays, and because it predicts specific critical exponents and Kibble-Zurek signatures that experiments could probe.","feed_headline":"Rydberg ladders order in two steps: bond wave first, reflection later","feed_subtitle":"Simulations map the ax = 2ay ladder where bond-order waves come first and Ising transitions break reflection later.","key_machinery":"The central objects are the rung bond-order operator B_i = a†_{i,1}a_{i,2} + a†_{i,2}a_{i,1} and the rung density-difference operator m_i = n_{i,2} − n_{i,1}, whose correlations distinguish the Z_p^+ phases from the Z_2p phases and define the two floating phases. The argument is carried by structure factors peaked at the incommensurate wave vector, Binder cumulant data collapse for critical points and exponents, entanglement entropy to locate phase boundaries, and finite-size extrapolation of commensurate lines to locate the CFT points where the chiral transition lines cross them.","core_discovery":"The paper establishes that, in the ax = 2ay Rydberg ladder, the crystalline orders accessible from the disordered phase appear in a defined order: intermediate detuning produces Z_p^+ phases in which every pth rung carries a positive bond order, an entangled state of the form (|r,g⟩ + |g,r⟩)/√2 that preserves the top-bottom reflection symmetry; larger detuning then drives an Ising transition that breaks reflection symmetry, doubling the period to Z_2p phases with staggered occupations of the top or bottom leg. Along the commensurate lines with k = 2π/3 and k = π/2, the direct melting of the Z_3^+ and Z_4^+ orders into the disordered phase occurs through a three-state Potts CFT point and an Ashkin-Teller CFT point respectively, flanked by non-CFT chiral transition lines; the correlation-length exponent ν peaks and the dynamical exponent z dips at these CFT points, and the Kibble-Zurek exponent reaches about 0.46 and 0.45 there. Between crystalline orders the paper finds two floating phases: FL+, with quasi-long-range incommensurate bond-order correlations, and FL, with quasi-long-range incommensurate rung density-difference correlations, separated from the disordered phase by BKT transitions and from crystalline order by Pokrovsky-Talapov transitions.","pith_inferences":["If the two-step symmetry-breaking sequence is generic, similar bond-order-first ordering should appear in other quasi-1D arrays with strong intra-rung coupling, such as zigzag or triangular ladders, whenever the intra-rung interaction dominates the inter-rung interaction.","The twenty-one-rung interaction cutoff is a natural source of error: extending the range to 30 or more rungs could shift the lobe boundaries and the reported CFT coordinates, so an interaction-range extrapolation would be a direct robustness test.","The existence of two floating phases suggests a general rule that each crystalline order melts into a floating phase preserving its short-distance order type, so the number of floating phases should track the number of order-parameter species rather than merely the number of commensurate periods.","A concrete experimental probe of the Z_3^+ bond-order wave could use site-resolved readout of the rung correlators after a short evolution, looking for the period-3 oscillation in ⟨B_i⟩ with zero staggered density difference, which would distinguish it from the Z_6 phase without needing Bragg scattering."],"forward_implications":["The two-step ordering sequence gives a concrete prediction: in an ax = 2ay ladder, increasing detuning should first produce a symmetric bond-order wave with period p, and only at higher detuning should the top-bottom reflection symmetry break into a staggered period-2p density wave.","The Z_3^+ and Z_4^+ lobes should each host exactly one conformal point, of three-state Potts and Ashkin-Teller type respectively, and the Kibble-Zurek exponent should be maximal at those points, giving an experimentally accessible signature.","Two distinct floating phases should intervene between crystalline orders: FL+ adjacent to symmetric Z_p^+ order and FL adjacent to reflection-broken Z_p order, with BKT and Pokrovsky-Talapov transitions on their respective boundaries.","The aspect ratio of the ladder is decisive: unlike the ay = 2ax geometry studied in earlier work, the ax = 2ay geometry supports symmetric Z_p^+ orders, so geometry, not just blockade radius, determines which symmetry breaks first.","The phases and critical points found here should be realizable in Rydberg tweezer arrays with interactions truncated to twenty-one consecutive rungs, making the predicted exponents and order patterns testable in current quantum simulators."],"supporting_citations":[{"why":"Provides the experimental and numerical baseline for a two-leg Rydberg ladder with ay = 2ax, including the observed floating phase and the absence of Z_p^+ order, which the present aspect-ratio comparison builds on.","marker":"[25]"},{"why":"Supplies the earlier effective-Hamiltonian mapping of the lower part of the ax = 2ay phase diagram and numerical evidence for the Ising transitions among Z_2, Z_2 × Z_2, and Z_2^+ orders.","marker":"[38]"},{"why":"Defines the chiral universality class used to classify the non-CFT transitions along the Z_3^+ and Z_4^+ boundaries away from the CFT points.","marker":"[42]"},{"why":"Supplies the one-dimensional Rydberg chain phase diagram with Ising, Potts, chiral, and floating phases that the ladder results are compared against and extend.","marker":"[17]"},{"why":"Provides the framework for commensurate-to-incommensurate melting transitions and the role of floating phases for periodicities 3 and higher in Rydberg chains.","marker":"[15]"},{"why":"Establishes the floating-phase description of one-dimensional Rydberg Ising chains, which the FL and FL+ phases generalize to the two-leg ladder.","marker":"[24]"},{"why":"Shows Z_p^± bond-order density waves in a square ladder with short interaction cutoffs, providing the contrast that motivates retaining longer-range interactions in the present work.","marker":"[39]"},{"why":"Reports a numerical study of the ax = ay ladder with the same crystalline orders labeled differently, used by the paper as a comparison for how lattice geometry changes the phase diagram.","marker":"[74]"},{"why":"Demonstrates Kibble-Zurek measurements in Rydberg arrays and supports the paper's claim that the calculated Kibble-Zurek exponent peaks are experimentally testable.","marker":"[3]"}],"fun_headline_variants":["Rydberg ladder orders in two steps: bond wave, then reflection","Two-leg Rydberg ladder: bond orders before reflection break","Rydberg ladder phases: Potts, Ashkin-Teller, and floating waves","Potts and Ashkin-Teller points mark Rydberg ladder transitions","Rydberg ladder hosts floating phases between ordered states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The phase diagram rests on truncating van der Waals interactions to twenty-one consecutive rungs and on classifying the floating phases from Friedel oscillations and structure-factor peak heights in ladders up to 613 rungs, without an explicit demonstration of algebraic decay; if the cutoff shifts the lobe boundaries or CFT points, or if those oscillations are finite-size artifacts, the reported diagram would change.","fun_headline_variants_meta":{"raw":{"variants":["Rydberg ladder orders in two steps: bond wave, then reflection","Two-leg Rydberg ladder: bond orders before reflection break","Rydberg ladder phases: Potts, Ashkin-Teller, and floating waves","Potts and Ashkin-Teller points mark Rydberg ladder transitions","Rydberg ladder hosts floating phases between ordered states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001243,"raw_usage":{"total_tokens":5215,"prompt_tokens":1176,"completion_tokens":4039,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":792,"completion_tokens_details":{"reasoning_tokens":3946}},"tokens_in":792,"tokens_out":4039,"duration_ms":24699,"temperature":1.0,"reasoning_tokens":3946,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:11:20.326499+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the phase diagram with interactions kept over 30, 40, or 60 consecutive rungs and check whether the Z_3^+ lobe boundary and the claimed CFT point at (Δ/Ω, Rb/a) = (2.7901, 2.3749) move outside the reported error bars, and independently fit the correlation functions C_B(i,r) and C_m(i,r) in the two floating phases at L = 613 to power-law decay; if the fitted exponents are not algebraic or the oscillation amplitudes vanish with system size, the FL and FL+ classification would not hold.","supporting_citations":[{"cited_title":"Zhang, S.-W","cited_arxiv_id":null,"evidence_quote":"Supplies the earlier effective-Hamiltonian mapping of the lower part of the ax = 2ay phase diagram and numerical evidence for the Ising transitions among Z_2, Z_2 × Z_2, and Z_2^+ orders."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the one-dimensional Rydberg chain phase diagram with Ising, Potts, chiral, and floating phases that the ladder results are compared against and extend."},{"cited_title":"Chepiga and F","cited_arxiv_id":null,"evidence_quote":"Provides the framework for commensurate-to-incommensurate melting transitions and the role of floating phases for periodicities 3 and higher in Rydberg chains."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports a numerical study of the ax = ay ladder with the same crystalline orders labeled differently, used by the paper as a comparison for how lattice geometry changes the phase diagram."}],"review_version":1}