{"id":"9cc1c128-98bd-4634-87cd-96782ab5c0e0","arxiv_id":"2505.23409","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A self-consistent parton mean-field derivation, with no fitting, reproduces the chiral spin liquid ansatz for a Rydberg model and yields the twofold ground-state degeneracy that exact diagonalization missed.","lead":"By solving a parton mean-field theory self-consistently, the authors derive the parameters of a chiral spin liquid in a Rydberg lattice from the microscopic Hamiltonian rather than by fitting. Their projected state matches exact diagonalization and shows the twofold ground-state degeneracy expected for a topological spin liquid.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The self-consistent solution at g=0.7 (real χ1 with eχ↑=eχ↓) contradicts the spin-symmetry condition Eq. (19) used to initialize it, so the PSG class assignment and CSL interpretation rest on an unexplained, load-bearing assumption.","rationale":"The paper's strongest contribution is a parameter-free self-consistent derivation of the parton amplitudes and the overlap with exact-diagonalization states, which is genuine evidence. It also honestly states limitations around phase transitions and finite-size overlaps. The load-bearing weakness is the unresolved incompatibility between Eq. (19) and the reported solution; because the PSG classification step is exactly what upgrades the ED overlap into a statement about topological order, an unexamined symmetry violation there propagates to the central claim. I do not see a second issue that is more central: the 10^-2 degeneracy test is weak but secondary, and the lack of code/data affects reproducibility rather than the logical argument. A direct spin-symmetry test of the projected wavefunction would settle the matter. I therefore keep the CONDITIONAL verdict pending that check.","tokens_in":9154,"tokens_out":10786,"duration_ms":120446,"concrete_test":"At g=0.7 on the cluster used for Fig. 5 (e.g., 24c), construct the Gutzwiller-projected wavefunction from the reported converged parameters and compute ⟨S_tot^2⟩ (equivalently, verify ⟨S^x_i S^x_j⟩ = ⟨S^y_i S^y_j⟩ = ⟨S^z_i S^z_j⟩ on all bonds). If the projected state is not an SU(2) singlet, the solution violates the spin-rotation condition Eq. (19) at the physical level, and the identification with PSG class 1 of Ref. [19] would need to be re-examined under the actual reduced symmetry group. This test distinguishes the two possible resolutions: projection restores SU(2) (concern resolved) versus projection does not (central claim unsupported).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The weakest load-bearing step is in Sec. III: after stating Wen's spin-rotation condition, Eq. (19), eχ↑ = −eχ†↓, and initializing the ↓ correlations with it, the self-consistency loop is run independently for ↑ and ↓ and is reported to converge to eχ1↑ = eχ1↓ = −0.238 (real). These two statements are incompatible for a real nonzero χ1: Eq. (19) would require χ1,↑ = −χ1,↓* = χ1,↓ only if χ1 = 0 or purely imaginary. The paper neither resolves this nor proves that the converged point is the physical saddle point. This matters because the subsequent identification of V = vσ0, W/v ≈ −0.26(1+i)σ0 with ansatz class 1 of Ref. [19] uses the PSG classification, which is only valid for the symmetry structure assumed there. If the actual mean-field ansatz breaks the symmetry encoded in Eq. (19), the PSG assignment and the topological-degeneracy argument do not follow without a new classification for the reduced (non-SU(2)) symmetry of the Rydberg model. The ED overlap is real evidence but cannot distinguish a valid CSL from a symmetry-broken trial state on small clusters.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a fermionic parton mean-field theory for Rydberg excitations on a honeycomb lattice with nearest-neighbor hopping and density-dependent, complex next-nearest-neighbor hopping. Starting from the microscopic hard-core boson Hamiltonian, the authors decouple the interaction in terms of spinon hopping amplitudes, solve the resulting mean-field equations self-consistently, and obtain real nearest-neighbor and complex next-nearest-neighbor spinon hoppings without fitting to numerics. They identify these amplitudes with PSG ansatz class 1 of Ref. [19], report large Gutzwiller-projected wavefunction overlaps with exact-diagonalization ground states in the interval 0.4 ≲ g ≲ 0.9, and claim a twofold topological degeneracy under twisted boundary conditions as evidence for a chiral spin liquid. The central difficulty is that the reported converged solution appears incompatible with the spin-symmetry condition Eq. (19) used to initialize the iteration, and the degeneracy evidence is supported only by a loose numerical tolerance.","tokens_in":9339,"tokens_out":4190,"duration_ms":50672,"significance":"If the self-consistent construction is valid, the paper would provide a parameter-free, microscopic parton mean-field description of a chiral spin liquid in a Rydberg platform, with explicit analytical expressions for the spinon hoppings and testable predictions for chirality, spin correlations, and wavefunction overlaps. This would be a valuable step beyond the variational PSG fitting of Ref. [19]. The absence of any fitted parameters and the direct connection to the microscopic Hamiltonian are genuine strengths, as is the reproduction of ED spin-chirality and overlap data without tuning. However, the symmetry inconsistency identified in the manuscript undermines the PSG classification and the associated spin-rotation invariance, and the topological-degeneracy test is not yet quantitative enough to support the central claim. The paper's significance is therefore conditional on resolving these load-bearing issues.","major_comments":[{"comment":"The text states the spin-symmetry condition eχ↑ = −eχ†↓ and initializes the ↓ correlations with it, but the self-consistent iteration is then run independently for ↑ and ↓ and is reported to converge to eχ1,↑ = eχ1,↓ = −0.238 and eχ2,↑ = eχ2,↓ = 0.019 + i0.280. These values are incompatible with Eq. (19): for real nonzero eχ1, Eq. (19) requires eχ1,↑ = −eχ1,↓, not equality, and for eχ2 it requires eχ2,↑ = −(eχ2,↓)*, which for eχ2,↑ = eχ2,↓ = 0.019 + i0.280 is 0.019 + i0.280 = −0.019 + i0.280, a contradiction. Since Eq. (19) is the link to Wen's spin-liquid construction and to the PSG classification of Ref. [19], the converged fixed point does not, as written, satisfy the symmetry structure used to classify it. The authors need to derive a corrected symmetry condition for their independent-↑/↓ decoupling or prove that the equal-spin fixed point is nevertheless the physical saddle point, e.g., by showing it is invariant under the projective symmetry group of the original Hamiltonian.","section":"Sec. III, Eq. (19) and Fig. 2"},{"comment":"The twofold topological degeneracy is inferred from two zero eigenvalues of the 4×4 overlap matrix 'up to a tolerance of 10^-2', but no scale is given for this tolerance. The relevant test is whether the two small eigenvalues are separated from the nonzero eigenvalues by a clear gap, and how this separation behaves with system size. A tolerance of 10^-2 could be comparable to typical overlap magnitudes on the 24-site lattices used in the ED comparison. The claim of topological degeneracy requires either a quantitative eigenvalue hierarchy or finite-size scaling of the overlap matrix; as it stands, the evidence is not robust enough to support the central CSL conclusion.","section":"Sec. IV, ground-state degeneracy discussion"},{"comment":"The self-consistency loop is initialized with random values in (0,1), but the paper provides no basin-of-attraction study, no demonstration that the iteration converges to a unique fixed point independent of the starting values, and no comparison of the free energy of the converged solution with other possible fixed points. The conclusion states that 'the self-consistent solution yields a unique mean-field Hamiltonian', but that uniqueness is not established. If multiple fixed points exist with comparable energies, the physical selection rule must be given. This is load-bearing because the PSG assignment and the subsequent CSL interpretation rely on the converged solution being the correct physical saddle point.","section":"Sec. III, iteration scheme and closure"}],"minor_comments":[{"comment":"The antisymmetric tensor εαβ is used without definition; please specify its orientation and explain more concretely why setting η=0 is exact within the half-filled particle-number-conserved subspace, beyond the brief statement given.","section":"Sec. III, Eq. (5)"},{"comment":"The caption mentions 'different lattice shapes', but the specific shapes, sizes, and boundary conditions are not listed. Please identify the finite lattices used for the ED and mean-field overlaps, since the text notes that phase-transition points depend heavily on the ED cluster shape.","section":"Sec. IV, Fig. 5"},{"comment":"The overlap formula uses a single |ψ_ED⟩, but the text discusses ED ground states in the context of degeneracy; please clarify whether the ED ground state is numerically degenerate and how the overlap is computed in that case.","section":"Sec. IV, Eq. (21)"},{"comment":"The sentence 'The mean-field approximation to Hamiltonian (1), obtained by replacing the number operators by their expectation value at half filling' could be misinterpreted, since the self-consistent calculation later keeps the density-dependent phase fluctuations beyond the simple replacement; please rephrase to distinguish this initial illustration from the full parton treatment.","section":"Sec. II, model Hamiltonian"},{"comment":"The paper would benefit from a short statement on reproducibility, e.g., whether the self-consistent iteration and the twisted-boundary-condition diagonalization are implemented in publicly available code; this is not required but would strengthen confidence in the numerical claims.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The reader's report and my own reading agree on the main weakness: the converged solution contradicts the symmetry condition Eq. (19), and the degeneracy check is too loosely quantified. I view this as fixable in a revision if the authors can either reconcile the symmetry condition with their decoupling, or appropriately qualify the PSG classification claim while preserving the honest ED-overlap evidence. The paper is within the scope of the journal and the self-consistent approach is a genuine step forward if the symmetry issue is resolved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe headline is twofold: this paper does something genuinely useful and something genuinely shaky. The useful part is deriving the parton mean-field parameters for the Rydberg honeycomb model of Ref. [18] self-consistently from the microscopic Hamiltonian, rather than fitting them as in Ref. [19]. That is a real step forward, and the resulting overlap with ED wavefunctions in g ∈ [0.4, 0.9] is a credible check. The twofold degeneracy on a torus, while not conclusive, is a new observation that small-cluster ED could not provide.\n\nThe shaky part is an internal contradiction at the center of the argument. Sec. III states the spin-symmetry condition eχ↑ = −eχ†↓ (Eq. 19), correctly translated from Wen's χ↑ = χ↓, and uses it to initialize the ↓ correlations. The self-consistency loop is then run independently for the two species and converges to eχ1,↑ = eχ1,↓ = −0.238, both real. That satisfies Eq. (19) only if χ is zero or purely imaginary. The paper calls the convergence \"remarkable\" but never notes the conflict. This matters because the next step—identifying V = vσ0 and W/v ≈ −0.26(1+i)σ0 with ansatz class 1 of Ref. [19]—relies on the PSG classification, which assumes the symmetry encoded in Eq. (19). If the true solution breaks that symmetry, the classification may not apply, and the degeneracy claim does not follow. The authors need to explain why the symmetry condition is not required, or show that a corrected calculation still lands in the same universality class.\n\nOther issues are secondary. The degeneracy test via overlap-matrix eigenvalues \"zero up to 10^-2\" is loose; a system-size or robustness check would help. No code or data, so reproducibility is limited. The uniqueness of the saddle point is asserted rather than proved, but that is common in this literature and not a major flaw. To give credit where it is earned: the ED overlap and spin-chirality comparisons are appropriate evidence, honestly presented, and the authors explicitly note that their mean-field cannot predict phase transitions.\n\nThis paper is for people working on parton mean-field descriptions of spin liquids and on Rydberg-based quantum simulators. It deserves a serious referee: the method is a genuine contribution and the open problem is important. But the referee should be instructed to focus on the Eq. (19) discrepancy; if it is not resolved, the central classification claim should not stand.\n\nMy recommendation: send to peer review with that red flag clearly raised.\n\nBest.","headline":"A useful self-consistent parton derivation for the Rydberg CSL model, but the converged solution violates its own spin-symmetry condition, leaving the PSG classification unsupported.","tokens_in":9944,"tokens_out":5556,"would_cite":false,"duration_ms":59894,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A parameter-free parton mean-field solution places the Rydberg model's ground state in a chiral spin liquid class and reproduces exact diagonalization results.","keywords":["Rydberg atoms","quantum spin liquid","chiral spin liquid","parton mean-field theory","projective symmetry group","density-dependent Peierls phase","honeycomb lattice","Gutzwiller projection"],"falsifier":"Compute the energy landscape of the mean-field Hamiltonian over the full space of complex $\\chi_{1,\\uparrow}$, $\\chi_{1,\\downarrow}$, $\\chi_{2,\\uparrow}$, $\\chi_{2,\\downarrow}$ (relaxing the spin-symmetry condition) and check whether the reported fixed point is the global minimum; alternatively, test spin-rotation invariance of the Gutzwiller-projected ground state by evaluating the static spin-structure factor $S(k)$ for residual magnetic order, or run ED on larger clusters (e.g., 36 sites) to see if the twofold degeneracy persists and matches the parton prediction.","tokens_in":8840,"feed_emoji":"🌀","tokens_out":8441,"duration_ms":70259,"temperature":0.7,"pith_summary":"The paper aims to show that the Rydberg Hamiltonian on a honeycomb lattice with density-dependent Peierls phases has a chiral spin liquid ground state, by deriving the spin-liquid wavefunction from first principles rather than by fitting ansatz parameters to numerics. The authors fermionize the hard-core boson model, perform a parton mean-field decoupling, and solve the resulting self-consistency equations analytically in momentum space. The self-consistent solution yields hopping amplitudes $V = v \\sigma_0$ and $W/v \\approx -0.26(1+i) \\sigma_0$ at $g=0.7$, which places it in the same projective symmetry group class as the previously fitted ansatz. Its Gutzwiller-projected ground state shows large overlap with exact diagonalization states for $0.4 \\lesssim g \\lesssim 0.9$ and exhibits a twofold topological degeneracy on the torus, the hallmark of a chiral spin liquid.","feed_headline":"Parton theory derives Rydberg spin-liquid state with no fitting","feed_subtitle":"Self-consistent mean-field solution matches exact numerics and shows the twofold degeneracy of a chiral spin liquid.","key_machinery":"The central object is the fermionic parton representation of the hard-core boson operators, with spinon operators $f_{i,\\alpha}$ replacing the bosons; the mean-field decoupling of the sixth-order term leaves hopping amplitudes $\\chi_{1,\\alpha}$ (nearest neighbor) and $\\chi_{2,\\alpha}$ (next-nearest neighbor) to be determined. The self-consistency equations (16) and (17) express these amplitudes as Brillouin-zone integrals over the lower-band eigenvector components $(a(k), c(k))$ and the lattice phase factors $\\varphi_1(k)$, $\\varphi_2(k)$, evaluated for both spin species independently. The Gutzwiller projection then maps the parton ground state back to the physical bosonic Hilbert space. This machinery is what converts the microscopic Hamiltonian into the concrete prediction $V = v \\sigma_0$, $W/v \\approx -0.26(1+i) \\sigma_0$.","core_discovery":"The paper's central claim is that the self-consistent parton mean-field solution of the Rydberg Hamiltonian (1), obtained without any fitting to numerics, produces hopping amplitudes $V = v \\sigma_0$ and $W/v \\approx -0.26(1+i) \\sigma_0$ at $g=0.7$ (eq. 20), placing it in PSG ansatz class 1 of Ref. [19]. The Gutzwiller-projected ground state of this mean-field Hamiltonian has large overlap with exact-diagonalization ground states in the window $0.4 \\lesssim g \\lesssim 0.9$, and shows a twofold topological degeneracy under twisted boundary conditions that ED could not resolve. The paper presents these results as evidence that the model hosts a chiral spin liquid, and as a microscopic derivation of the previously ad-hoc PSG ansatz.","pith_inferences":["The reported solution violates the spin-symmetry condition (19) that it is supposed to satisfy ($\\chi_1$ real rather than purely imaginary), so the spin-rotation invariance of the resulting chiral spin liquid remains an open question that the paper does not resolve.","A natural next test would be to compute the topological entanglement entropy or modular matrices from the projected parton wavefunction, which would confirm the topological order independently of the degeneracy count.","The self-consistent method could be extended to other lattice geometries or to finite doping to see whether the chiral spin liquid survives away from half filling."],"forward_implications":["The self-consistent mean-field solution places the Hamiltonian in PSG ansatz class 1 of Ref. [19] without any fitting, giving a microscopic origin to the previously ad-hoc hopping amplitudes.","The twofold topological degeneracy on the torus, which exact diagonalization could not resolve, emerges naturally in the parton description and is a hallmark of a chiral spin liquid.","The overlap with exact-diagonalization ground states is large for $0.4 \\lesssim g \\lesssim 0.9$, indicating that the projected parton wavefunction captures the physics of the frustrated regime.","Because the self-consistent solution is computed analytically in momentum space, the approach applies to arbitrary system sizes and can access the thermodynamic limit."],"supporting_citations":[{"why":"Supplies the microscopic Rydberg Hamiltonian and the exact-diagonalization evidence for a chiral spin liquid, including the absence of torus degeneracy that the parton approach restores.","marker":"[18]"},{"why":"Provides the projective symmetry group classification and the fitted ansatz (no. 1) whose form the self-consistent solution reproduces without fitting.","marker":"[19]"},{"why":"Establishes the PSG classification scheme for chiral spin liquids used to place the mean-field Hamiltonian in the correct symmetry class.","marker":"[20]"},{"why":"Introduces the fermionic parton construction and the spin-symmetry condition (eq. 19) that the self-consistency iteration is intended to respect.","marker":"[21]"},{"why":"Details the twisted-boundary-condition method used to extract the twofold topological degeneracy of the projected wavefunction.","marker":"[30]"}],"fun_headline_variants":["Parton mean-field derives Rydberg CSL without fitting","No-fit parton ansatz matches Rydberg spin liquid","Self-consistent partons yield chiral Rydberg liquid","Parton theory captures Rydberg spin liquid without fits","Self-consistent partons derive Rydberg CSL without fitting"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the self-consistently converged mean-field fixed point is the physically correct ground state; the paper notes the spin-symmetry condition (19) from the projective construction, but its reported solution (real $\\chi_1 = -0.238$ with $\\chi_{1,\\uparrow} = \\chi_{1,\\downarrow}$) is incompatible with that condition, which would require purely imaginary $\\chi_1$, and the paper does not prove this fixed point is the global saddle point nor explain the discrepancy.","fun_headline_variants_meta":{"raw":{"variants":["Parton mean-field derives Rydberg CSL without fitting","No-fit parton ansatz matches Rydberg spin liquid","Self-consistent partons yield chiral Rydberg liquid","Parton theory captures Rydberg spin liquid without fits","Self-consistent partons derive Rydberg CSL without fitting"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000547,"raw_usage":{"total_tokens":2572,"prompt_tokens":862,"completion_tokens":1710,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":478,"completion_tokens_details":{"reasoning_tokens":1625}},"tokens_in":478,"tokens_out":1710,"duration_ms":13433,"temperature":1.0,"reasoning_tokens":1625,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:46:41.813494+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the energy landscape of the mean-field Hamiltonian over the full space of complex $\\chi_{1,\\uparrow}$, $\\chi_{1,\\downarrow}$, $\\chi_{2,\\uparrow}$, $\\chi_{2,\\downarrow}$ (relaxing the spin-symmetry condition) and check whether the reported fixed point is the global minimum; alternatively, test spin-rotation invariance of the Gutzwiller-projected ground state by evaluating the static spin-structure factor $S(k)$ for residual magnetic order, or run ED on larger clusters (e.g., 36 sites) to see if the twofold degeneracy persists and matches the parton prediction.","supporting_citations":[{"cited_title":"Ohler, M","cited_arxiv_id":null,"evidence_quote":"Supplies the microscopic Rydberg Hamiltonian and the exact-diagonalization evidence for a chiral spin liquid, including the absence of torus degeneracy that the parton approach restores."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the projective symmetry group classification and the fitted ansatz (no. 1) whose form the self-consistent solution reproduces without fitting."},{"cited_title":"Bieri, C","cited_arxiv_id":null,"evidence_quote":"Establishes the PSG classification scheme for chiral spin liquids used to place the mean-field Hamiltonian in the correct symmetry class."},{"cited_title":"Mei and X.-G","cited_arxiv_id":null,"evidence_quote":"Details the twisted-boundary-condition method used to extract the twofold topological degeneracy of the projected wavefunction."}],"review_version":1}