{"id":"42e471d7-cf5d-482a-93aa-bfc29e284ac1","arxiv_id":"2502.04165","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In an exactly solvable long-range cluster XY chain, the exponents ν and z vary with the decay exponent α as z = α-1 for 1<α<2 and z = 1 for α≥2, with ν = 1/z.","lead":"This paper exactly solves a chain of spins with long-range cluster interactions and shows that the critical exponents of its quantum phase transition vary continuously with the interaction decay rate. The result adds a new solvable model to long-range quantum criticality and sharpens expectations for Rydberg-atom quantum simulators.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The z=α-1 result hinges on the M=N/2 cutoff sequence; the paper does not prove this reproduces the thermodynamic dynamical exponent for α<2.","rationale":"The reader's weakest assumption is exactly the point I would stress: the z=α-1 result is obtained from a finite-size gap in a model whose interaction range is tied to system size, and the paper does not prove that this sequence has the same low-energy limit as the infinite-range model. I do not see an alternative internal inconsistency: the free-fermion solution is standard and the FSS numerics, while not definitive, are consistent with the analytic trend. The main gap is the missing derivation or justification of the order of limits. I would therefore keep the paper conditional rather than rejecting it: the central claim is plausible and testable. The proposed check settles the issue analytically. If the check confirms ε(π/N)∼N^{1-α}, the paper should add this argument and the conditional can be lifted; if not, the claim of continuously varying z at hc1 would need to be withdrawn or qualified. The numerical data-collapse issues noted by the reader are secondary because the analytic derivation, once the order of limits is justified, does not depend on them.","tokens_in":13459,"tokens_out":16259,"duration_ms":168381,"concrete_test":"Evaluate the exact M=∞ dispersion at hc1=-ζ(α) for 1<α<2: compute ε(k)=2[δ_k^2+ϵ_k^2]^{1/2} with δ_k=γΣ_{m=1}∞ m^{1-α} sin(mk) and ϵ_k=Σ_{m=1}∞ m^{-α}(cos mk-1), and set k=π/N. If ε(π/N) ∼ C N^{1-α} for large N with C finite and nonzero, then z=α-1 is the true thermodynamic exponent and the cutoff sequence is innocuous. If instead ε(π/N) ∼ C N^{-1}, the M=N/2 truncation misrepresents the thermodynamic limit and z must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of z=α-1 for 1<α<2 is built on Eqs. (18)-(20): the gap is written as Δ_N ∼ 2γ(Σ_{m=1}^{N/2} m^{1-α})|k| and then |k|∼1/N. For α<2 the coefficient Σ m^{1-α} ∼ N^{2-α}, so Δ_N ∼ N^{1-α}, i.e., z=α-1. The load-bearing step is the identification of this finite-size gap with the thermodynamic dynamical exponent. The linear-in-k expansion of δ_k in Eq. (18) is not valid for the infinite-range model when α<2: the velocity Σ_{m=1}∞ m^{1-α} diverges, and the true dispersion is δ_k ∼ |k|^{α-1}. If one instead takes the thermodynamic limit at fixed interaction range M, the coefficient is finite and the same argument gives z=1. The paper never proves that the M=N/2 cutoff sequence and the M=∞ critical theory yield the same low-energy exponents; the Supplemental asymptotic analysis (B.1)-(B.7) only analyzes the cutoff model. Thus the central claim that z=α-1 is the thermodynamic exponent rests on an unproven interchange of limits.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a spin-1/2 long-range cluster XY chain with algebraically decaying couplings, maps it to free fermions via the Jordan-Wigner transformation, and derives the critical fields hc1 and hc2 and the relation νz=1. The main claim is that the dynamic exponent z (and therefore ν via νz=1) varies continuously with the decay exponent α, with z=α−1 for 1<α<2 and z=1 for α≥2. The authors then present finite-size scaling analyses of correlation-function derivatives and fidelity susceptibility, using data collapse with adjustable exponents, and compute the entanglement entropy to extract a central charge c that also varies with α. The free-fermion mapping and the formulas for the critical fields are clean, but the derivation of the α-dependent z contains a technical error that undermines the central claim.","tokens_in":13748,"tokens_out":30532,"duration_ms":315755,"significance":"If correct, the result would be notable: an exactly solvable one-dimensional spin model with continuously varying critical exponents and a varying effective central charge is of interest for long-range interacting quantum simulators. The manuscript's strengths are its transparent free-fermion solution, explicit formulas for the critical fields, and the detailed numerical scaling study. However, the load-bearing step that produces z=α−1 is based on an invalid linearization of sin(mk), and the numerical verification treats the exponents being verified as fitting parameters. As it stands, the central claim is not supported and the paper would require a substantially revised derivation and a re-evaluation of the numerical evidence.","major_comments":[{"comment":"The derivation of the central result uses an invalid uniform expansion sin(mk)≈mk: for k∼1/N and m up to M=N/2, mk is of order unity. A correct evaluation gives δ_k∼N^{2−α} for α<2 (or δ_∞(k)∼k^{α−2} for 2<α<3), so the gap does not scale as N^{−(α−1)}. This invalidates z=α−1.","section":"Sec. II, Eqs. (18)–(20) and Supplemental Eqs. (A.1)–(B.7)"},{"comment":"The verification section treats the predicted exponents as fitting parameters, and the fitted values deviate from the analytic predictions by up to 14% (α=2.0, ν=1.1438 vs ν=1). The collapse therefore does not independently confirm the claimed exponents.","section":"Sec. III, Eq. (24), Fig. 4, and Table I"},{"comment":"The CFT scaling form used to extract c is not justified for α<2, where the single-particle spectrum has no conventional low-energy continuum and the finite-size gap does not close.","section":"Sec. III, Eq. (30) and Fig. 7"}],"minor_comments":[{"comment":"The determinant in Eq. (22) uses the index j in the last row without defining it; the row index should presumably be i+r−1 rather than j−1.","section":"Eq. (22)"},{"comment":"The caption mentions dashed lines indicating the location of the critical point as α→∞, but no dashed lines are visible in the figure.","section":"Fig. 1 caption"},{"comment":"The text notes that the ground-state energy diverges for α≤1, but the subsequent analysis and figures use α>1; this restriction should be stated explicitly in the model definition.","section":"Sec. II, text before Eq. (2)"},{"comment":"The phrase 'To optimize scaling behavior, we verify these critical exponents ... by adjusting fitting parameters' is self-contradictory, since verification by adjusting the exponents being verified is circular; this phrasing should be revised.","section":"Abstract"}],"recommendation":"reject","confidential_remarks":"The main result is not a local presentation issue: the derivation of z=α−1 is based on an invalid uniform linearization of sin(mk) in the finite-size gap. The error can be checked directly by evaluating the finite sum, and it affects the central claim of the paper. The numerical sections do not remedy this because the exponents are treated as adjustable fitting parameters. I recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plain: this paper has one genuinely new result and one serious presentation problem. The new result is the closed-form dynamic exponent in the long-range cluster XY chain: z=α−1 for 1<α<2 and z=1 for α≥2, with ν=1/z. The derivation is transparent — gap closes linearly in k, the coefficient Σ m^{1−α} gives the N^{2−α} factor, and the Euler–Maclaurin analysis in the Supplemental is correct. I checked the order-of-limits worry. For α<2 the infinite-range dispersion δ_k ~ k^{α−1} gives the same N^{1−α} gap exponent as the M=N/2 cutoff, so the stress-test concern is not fatal; the authors should prove or state this, but the result survives.\n\nWhat is not solid is the numerical section. Equations (24) and (28) treat ν as a free fitting parameter and then use the resulting collapse as evidence for the predicted ν. That is circular: the collapse is optimized, as the abstract admits. Worse, Table I reports ν=1.1438 at α=2 while the analytic prediction is ν=1, and the corresponding z=0.8743 conflicts with z=1. The authors call this a logarithmic correction, but log corrections do not change the leading exponent; a 14% shift needs a quantitative explanation. The Supplemental's own convergence numbers (ν→1.08 at α=2) still do not reach 1. No error bars are provided anywhere.\n\nThe central charge part is the weakest. For α<2 the low-energy dispersion is non-linear, so the model is not conformal and Eq. (30) is not justified. The reported c depends on γ, which is a red flag that the fitted slope is a finite-size artifact, not a universal quantity.\n\nWho is this for: people working on long-range spin chains and quantum criticality. The exact z(α) formula is a useful data point and likely correct. The paper deserves a serious referee, but the referee should require (i) a statement about the cutoff regularization and the α=2 log corrections, (ii) a collapse using the analytic ν without free fitting, or at least an honest discussion of the discrepancy, and (iii) reframing the c results as effective finite-size slopes. If those are fixed, this is a solid paper.","headline":"The analytic z(α)=α−1 result is solid and new, but the numerical verification is circular and Table I contradicts the analytic prediction at α=2; still worth refereeing.","tokens_in":14262,"tokens_out":4375,"would_cite":true,"duration_ms":45887,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In an exactly solvable antiferromagnetic cluster XY chain with interactions decaying as $1/r^\\alpha$, the critical exponents $\\nu$ and $z$ vary continuously with $\\alpha$, and the identity $\\nu z=1$ still holds.","keywords":["long-range interactions","cluster XY model","critical exponents","dynamic critical exponent","correlation length exponent","Jordan-Wigner transformation","fidelity susceptibility","entanglement entropy"],"falsifier":"For $1<\\alpha<2$, compute the exact gap at $h_{c1}$ from the free-fermion spectrum for a sequence of increasing $N$ with the interaction range held at $M=N/2$, and fit $z$ from $\\Delta_N \\propto N^{-z}$; if the fitted $z$ does not approach $\\alpha-1$ as $N\\to\\infty$, the claim fails. A complementary check is to construct the low-energy continuum theory by taking $N\\to\\infty$ before expanding around $k=0$: if the resulting velocity diverges and the gap instead closes as $N^{-1}$, the order of limits changes the exponent.","tokens_in":13270,"feed_emoji":"🧲","tokens_out":13593,"duration_ms":109546,"temperature":0.7,"pith_summary":"The paper studies a spin-1/2 cluster XY chain in a transverse field whose interactions decay as $J_m = m^{-\\alpha}$. Because the model maps exactly to free fermions, the finite-size gap can be computed and the dynamic critical exponent $z$ read off; the result is $z = \\alpha-1$ for $1<\\alpha<2$ and $z=1$ for $\\alpha\\ge 2$. Together with $\\nu z = 1$, this gives a correlation-length exponent $\\nu$ that also varies continuously with $\\alpha$, unlike the long-range antiferromagnetic Ising chain where the exponents are fixed. The exponents are verified by scaling collapse of the derivative of the correlation function and of fidelity susceptibility, and the central charge $c$ extracted from entanglement entropy also decreases with $\\alpha$ toward $0.5$. If correct, this provides a rare exactly solvable example in which long-range cluster interactions tune the universality class.","feed_headline":"Critical exponents vary continuously in a solvable spin chain","feed_subtitle":"For $1<\\alpha<2$ the dynamic exponent is $\\alpha-1$; for $\\alpha\\ge2$ it is $1$, making the universality class tunable.","key_machinery":"The central object is the free-fermion representation of the cluster XY Hamiltonian obtained by Jordan-Wigner transformation, whose Bogoliubov quasiparticle energies are $\\varepsilon_k = 2\\sqrt{\\delta_k^2 + \\epsilon_k^2}$ with $\\epsilon_k = \\sum_{m=1}^{M} m^{-\\alpha}\\cos(m k) + h$ and $\\delta_k = \\gamma \\sum_{m=1}^{M} m^{-\\alpha}\\sin(m k)$. The gap closes at $k=0$ or $k=\\pi$ when $h$ reaches $h_{c1} = -\\sum_{m} m^{-\\alpha}$ or $h_{c2} = -\\sum_{m} (-1)^m m^{-\\alpha}$. The dynamic exponent is extracted from the finite-size gap $\\Delta_N \\sim |\\sum_{m=1}^{N/2} m^{1-\\alpha}|/N$, where the truncation $M=N/2$ is the longest possible interaction range under periodic boundary conditions; the Euler-Maclaurin asymptotic expansion of this sum produces the piecewise formula for $z$.","core_discovery":"The central claim is that in the long-range antiferromagnetic cluster XY model $H = \\sum_{j,m} J_m[(1+\\gamma)/2 \\sigma_j^x \\sigma_{j+m}^x + (1-\\gamma)/2 \\sigma_j^y \\sigma_{j+m}^y] \\prod_{p=j+1}^{j+m-1}\\sigma_p^z - h\\sum_j \\sigma_j^z$ with $J_m = m^{-\\alpha}$, the critical exponents are continuously varying functions of $\\alpha$. The derivation starts from the size dependence of the gap at the critical field $h_{c1}$: $\\Delta_N \\sim N^{-z}$ with $N^{-z} \\sim (1/N)\\left|\\sum_{m=1}^{N/2} m^{1-\\alpha}\\right|$. Euler-Maclaurin expansion of the generalized harmonic number yields $z = \\alpha-1$ for $1<\\alpha<2$ and $z=1$ for $\\alpha\\ge2$, with logarithmic corrections at $\\alpha=2$, and the relation $\\nu z = 1$ holds throughout. The same exponents are obtained from scaling collapse of the field derivative of the farthest correlation function and of fidelity susceptibility, and the central charge $c$ from half-chain entanglement entropy varies with $\\alpha$, converging to $0.5$ as $\\alpha\\to\\infty$.","pith_inferences":["A natural next test is whether the $\\alpha$-dependent exponents persist in the true thermodynamic limit if the momentum cutoff is taken before the system size; if they do, the critical theory for $1<\\alpha<2$ is nonlocal and has a scale-dependent effective velocity.","The variation of the central charge $c$ within a free-fermion model hints that the critical theory may be a rescaled free-fermion CFT with a cutoff-dependent Fermi velocity, rather than a genuinely different conformal field theory.","The same gap-scaling derivation could be applied to other cluster-type models with Jordan-Wigner strings, predicting similarly tunable exponents and providing a family of solvable long-range models for quantum simulation.","In a Rydberg-atom array with programmable power-law interactions and cluster terms, measuring the gap scaling or correlation length as a function of $\\alpha$ would directly test whether $z=\\alpha-1$ appears in a real system."],"forward_implications":["For all $\\alpha>1$, the identity $\\nu z = 1$ holds exactly, so measuring either exponent determines the other.","In the regime $\\alpha\\ge 2$, the model falls into the standard short-range XY/Ising universality class with $\\nu=1$ and $z=1$.","In the regime $1<\\alpha<2$, tuning the interaction decay exponent continuously changes both $\\nu$ and $z$, giving a one-parameter family of critical behaviors.","The central charge $c$ obtained from entanglement entropy at criticality is also $\\alpha$-dependent, approaching $0.5$ in the large-$\\alpha$ limit.","The exponents are confirmed numerically by scaling collapse of the correlation-function derivative and fidelity susceptibility, so the prediction is directly checkable in finite-size simulations."],"supporting_citations":[{"why":"Supplies the Jordan-Wigner transformation and free-fermion diagonalization that make the model exactly solvable.","marker":"[61]"},{"why":"Provides the finite-size scaling form $N^{\\beta/\\nu} f(|h-h_c|N^{1/\\nu})$ used to collapse the correlation-function data and extract $\\nu$.","marker":"[69]"},{"why":"Gives the long-range antiferromagnetic Ising chain result against which the continuously varying exponents of the cluster XY model are contrasted.","marker":"[27]"},{"why":"Sets the general framework for critical behavior of long-range interacting systems and motivates the question of how exponents depend on $\\alpha$.","marker":"[17]"},{"why":"Introduced the cluster Ising model and its free-fermion solution, the basis for the cluster-interaction terms studied here.","marker":"[35]"},{"why":"Discusses logarithmic corrections at $\\alpha=2$ in long-range models, invoked to explain numerical deviations from $\\nu=1$ at this special point.","marker":"[62]"},{"why":"Establishes fidelity susceptibility as a scaling probe of quantum critical points, used here to confirm the exponents independently.","marker":"[74]"},{"why":"Provides the conformal-field-theory formula for block entanglement entropy, used to extract the central charge $c$.","marker":"[75]"}],"fun_headline_variants":["Tunable critical exponents in exactly solvable spin chain","Long-range cluster XY shows α-dependent critical exponents","Universality class shifts with α in solvable cluster XY","Exact critical exponents vary with interaction range in spin model","Solvable long-range cluster XY has α-tunable exponents"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation relies on identifying the finite-size gap computed with the interaction range truncated at half the system size with the true thermodynamic dynamical exponent, without proving that the thermodynamic limit and the low-momentum limit commute for $1<\\alpha<2$.","fun_headline_variants_meta":{"raw":{"variants":["Tunable critical exponents in exactly solvable spin chain","Long-range cluster XY shows α-dependent critical exponents","Universality class shifts with α in solvable cluster XY","Exact critical exponents vary with interaction range in spin model","Solvable long-range cluster XY has α-tunable exponents"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00131,"raw_usage":{"total_tokens":5368,"prompt_tokens":1004,"completion_tokens":4364,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":620,"completion_tokens_details":{"reasoning_tokens":4284}},"tokens_in":620,"tokens_out":4364,"duration_ms":34463,"temperature":1.0,"reasoning_tokens":4284,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T23:17:42.383579+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For $1<\\alpha<2$, compute the exact gap at $h_{c1}$ from the free-fermion spectrum for a sequence of increasing $N$ with the interaction range held at $M=N/2$, and fit $z$ from $\\Delta_N \\propto N^{-z}$; if the fitted $z$ does not approach $\\alpha-1$ as $N\\to\\infty$, the claim fails. A complementary check is to construct the low-energy continuum theory by taking $N\\to\\infty$ before expanding around $k=0$: if the resulting velocity diverges and the gap instead closes as $N^{-1}$, the order of limits changes the exponent.","supporting_citations":[{"cited_title":"Defenu, T","cited_arxiv_id":null,"evidence_quote":"Sets the general framework for critical behavior of long-range interacting systems and motivates the question of how exponents depend on $\\alpha$."},{"cited_title":"Smacchia, L","cited_arxiv_id":null,"evidence_quote":"Introduced the cluster Ising model and its free-fermion solution, the basis for the cluster-interaction terms studied here."},{"cited_title":"Defenu, G","cited_arxiv_id":null,"evidence_quote":"Discusses logarithmic corrections at $\\alpha=2$ in long-range models, invoked to explain numerical deviations from $\\nu=1$ at this special point."}],"review_version":1}