{"id":"0b382d88-ef32-44b7-8741-4196cdb26c30","arxiv_id":"2505.05645","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"The dynamical critical exponent z of a fractional-order quantum Ising chain is claimed to equal q/2, with simulations showing continuously tunable entanglement light cones.","lead":"A chain of quantum spins with unusual fractional long-range couplings is claimed to spread entanglement faster or slower depending on a single number q, the fractional order. The paper predicts the spread exponent z=q/2 and reports matrix-product-state simulations for q from 0 to 2.5.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's own RG caveat undermines the claimed continuous tunability; the z=q/2 observation is not tested against the predicted q>2/3 crossover.","rationale":"The reader's verdict identifies the central issue: the numerical claim of tunable z(q) conflicts with the two-loop RG result cited in the paper itself, and the paper's only response is an untested finite-size hypothesis. My stress-test confirms that this is the load-bearing weakness. The derivation of z=q/2 in Section 4 has an apparent algebraic gap: Eq. (34) is asserted from Eqs. (27) and (33) without providing the step, and the q=2 check fails unless the hypergeometric modulus identity in Eq. (38) holds; the paper does not prove it. However, the primary rejection basis is the numerical evidence vs. the cited RG prediction: the model's own Section 5 states that for q>2/3 the flow drives z to 1, which is exactly what the data show at q>~2.5, and the authors admit the apparent mean-field window q up to 2 is likely a finite-size effect. Since the claim 'continuously tunable exponent z(q) for 0<q<2.5' requires the observed z=q/2 to be the asymptotic value, and the paper provides no test at larger L, the claim is unsubstantiated. The recommendation of REJECT is appropriate; the paper could become CONDITIONAL if the authors provide larger-scale data and a complete derivation, but as written the central claim does not hold. I agree with the reader's weakest_assumption, though I would also flag the analytic gap as a secondary independent issue. No ad hominem is involved; the critique is on the argument and evidence. The concrete test I propose is one specific check that would settle the finite-size vs. genuine-tunability question: rerun the two main extractions at L=400-1000 or use an infinite-MPS method. This is feasible with the exponential-sum MPO already implemented. If the z(q) curve is stable with L, the central claim survives; if it bends toward 1 for q>2/3, the claim collapses. This directly addresses the most load-bearing uncertainty.","tokens_in":12949,"tokens_out":2245,"duration_ms":21212,"concrete_test":"Repeat the finite-size gap scaling and entanglement-front extraction at q=1.5, 2.0, 2.2, and 2.4 for L=400, 600, and 1000 (or use iTPS/infinite-MPS with momentum-space methods) and compare the fitted z with the two-loop RG prediction. If z drifts toward 1 as L grows for any q>2/3, the observed z=q/2 is a crossover artifact and the tunability claim collapses; if z remains q/2 within error bars, the discrepancy with Ref. [46] is real and needs explanation. Additionally, independently evaluate the hypergeometric identity in Eqs. (38) and (34): check numerically whether |Gamma(q+1)/Gamma(q/2+1)^2 * 2F1(1,-q/2;(q+2)/2;e^{ik})|^2 - ((1/2)*binomial(q,q/2))^2 is indeed O(|k|^{2q}) for q in (0,2.5).","verdict_should_be":"REJECT","load_bearing_attack":"The central claim is that the fractional multiscale Ising model realizes a continuously tunable dynamical exponent z(q)=q/2 for 0<q<2.5. The numerical evidence, however, contradicts the two-loop RG result the authors themselves cite (Ref. [46], Maghrebi et al.), which states that for couplings J(r)~r^{-(1+q)} in 1D the mean-field relation z=q/2 holds only for q<2/3, while for q>2/3 fluctuations bend z toward 1. Section 5 explicitly acknowledges that the extracted exponents deviate from z=q/2 in the fluctuation-dominated regime and \"hypothesise[s] that these discrepancies arise from finite-size effects\" without testing larger L. The observed agreement with z=q/2 up to q~2 is therefore a potential crossover artifact: at q=2.2 the fit is 'subtle' and at q>~2.5 the system reverts to z=1, which is precisely the behavior expected if the asymptotic exponent is always 1 for q>2/3 and the apparent z=q/2 is a pre-asymptotic finite-size effect. The analytic derivation of z=q/2 is also not self-contained: Eq. (34) is presented as following from (27) and (33), but the intermediate step is not shown, and the q=2 sanity check Ek~|k| at gc fails unless the hypergeometric identity in Eq. (38) is verified; the modulus-squared subtraction in Eq. (38) is written without proof of the required asymptotic cancellation. Because the paper's own cited literature predicts the observed deviations, and because the finite-size escape hatch is explicitly untested, the central claim of a tunable z(q) confirmed across 0<q<2.5 is not supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a one-dimensional transverse-field Ising model with couplings generated by a discretized Riesz fractional derivative, J(r) ~ r^{-(1+q)} plus subleading corrections. It claims, via a truncated Jordan-Wigner mean-field calculation, that the dynamical critical exponent is z = q/2 (Eq. 41), and that MPS/TDVP simulations for 0 < q < 2.5 confirm a continuously tunable exponent: sublinear light cones for q < 2, a faint superlinear cone for 2 < q < 2.5, and a return to ballistic z = 1 for q ≳ 2.5. The numerical pipeline uses exponential-sum MPO fitting, finite-size gap scaling up to L = 200, and bond-entropy wavefront analysis.","tokens_in":13307,"tokens_out":11266,"duration_ms":118908,"significance":"If correct, the model would provide a one-parameter family of quantum light cones ranging from faster-than-linear to local, with a direct connection to Lévy statistics, and would be an interesting platform for near-term quantum simulators. The exponential-sum MPO representation and the dual extraction of z from gap scaling and wavefronts are reasonable technical contributions, and the comparison between these two estimators is a sensible cross-check. However, the central analytic derivation is internally inconsistent, and the numerical claim is not robust against the crossover predicted by the two-loop RG result the authors themselves cite. As it stands, the significance is not established.","major_comments":[{"comment":"Equation (27) does not follow from the preceding definition of ξ_k in Eq. (23). With ξ_k = 2g - 2J0 Σ J(r) cos(kr), the Bogoliubov dispersion is E_k^2/4 = (C_k - g)^2 + S_k^2 = C_k^2 + S_k^2 - 2g C_k + g^2, not C_k^2 + S_k^2 + 2g C_k + g^2 as written. Because of this sign error, at q = 2 and g_c = 1 the dispersion from Eq. (27) is 4|cos(k/2)|, which does not close at k = 0; the corrected sign gives 4|sin(k/2)| and closes at k = 0. Moreover, substituting Eq. (33) into Eq. (27) does not yield Eq. (34): the substitution produces cross terms linear in Re F that are absent from Eq. (34), and at q = 2 Eq. (34) gives E ∝ 3|sin(k/2)|, a factor mismatch with the direct evaluation. Thus the derivation of z = q/2 in Eq. (41) is not supported.","section":"§4, Eqs. (27), (34)"},{"comment":"The key asymptotic step is the claim in Eq. (38) that |F(q,k)|^2 − (B/2)^2 ∝ |k|^{2q}, where B = binomial(q, q/2). This cancellation is asserted after invoking the Euler transformation in Eq. (37), but the required cancellation of the constant and linear terms is not demonstrated. Since this subtraction is the entire reason that the |sin(k/2)|^q term is claimed to dominate, the conclusion E_k ∝ |k|^{q/2} is not established. A direct check at q = 2 contradicts the formula, as noted in the previous comment.","section":"§4, Eqs. (36)–(38)"},{"comment":"The numerical evidence does not support the abstract's claim of a confirmed continuously tunable exponent z(q) for 0 < q < 2.5. The paper cites Ref. [46], whose two-loop RG gives z = q/2 only for q < 2/3 and drives z toward 1 for q > 2/3. The data in Fig. 6 show z ≈ q/2 up to about q ≈ 2 and then a drift back toward z = 1, and Section 5 states that the authors 'hypothesise that these discrepancies arise from finite-size effects' without simulating larger systems. With L ≤ 200, the apparent tunable plateau may be a pre-asymptotic finite-size artifact, and the observed behavior at q ≈ 2.2–2.5 is exactly the crossover the cited RG predicts. The claimed superlinear cone is described as 'subtle', and the return to z = 1 at q ≳ 2.5 is not tested against the expected fluctuation-driven crossover. This is load-bearing because the central claim is continuous tunability.","section":"§5, Fig. 6 and the RG comparison"},{"comment":"The truncation of Jordan-Wigner strings to unity in Eq. (18) is an uncontrolled approximation. The argument that higher-order vertices carry higher canonical dimension is a mean-field/RG statement, and in one dimension the cited two-loop RG shows that fluctuations do change the exponent. Therefore the analytic result z = q/2 should be regarded only as a mean-field prediction for the model, not as a derivation for the spin chain itself; the numerical comparison cannot be used to confirm the mean-field exponent unless the thermodynamic-limit convergence of the numerics is established.","section":"§4, Jordan-Wigner truncation"}],"minor_comments":[{"comment":"The binomial coefficient with non-integer q and integer r, (q choose q/2 + r), is nonstandard; please state explicitly that it is the generalized binomial coefficient Γ(q+1)/(Γ(q/2+r+1)Γ(q/2−r+1)) and define J(0) when it first appears in Eq. (29).","section":"§2, Eq. (3)"},{"comment":"The green markers in Fig. 6 are described as the exponent for the 'asymptotic power law of the fractional interactions', but the fitting procedure for this quantity is not described in the text; please add a sentence or a reference to the relevant equation.","section":"§3.2 and §5, Fig. 6"},{"comment":"The relationship between the hypergeometric sum F(q,k) and the physical sums C_k and S_k should be stated more carefully, including the sign of S_k: with J(r) = (−1)^{r+1} binom(q, q/2+r), the sine sum J0 Σ J(r) sin(kr) equals −J0 Im F, not J0 Im F as written in Eq. (33). The spectrum is insensitive to this sign because S_k appears squared, but the formula as written is inconsistent with Eq. (3).","section":"§4, Eq. (32)"},{"comment":"The finite-size fits use several free parameters (a, b, ω, z, and similarly for g_c) on systems L = 10 to L = 200; please report the number of system sizes, the fit ranges, and the resulting parameter uncertainties so that the robustness of the extracted z(q) can be assessed.","section":"§3.2, Eq. (10)"}],"recommendation":"reject","confidential_remarks":"The decisive issue is that the manuscript's own Section 5 caveat is in tension with the abstract's central claim. The analytic derivation of z = q/2 contains an algebraic sign error and an unproven hypergeometric cancellation, and the numerical plateau z ≈ q/2 for q > 2/3 is explicitly attributed to possible finite-size effects without a test at larger L. These are load-bearing problems that the current manuscript does not resolve; a resubmission would need to correct the derivation, verify Eq. (38), and provide either larger-L data or a quantitative account of the crossover to z = 1."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me give you the short version: this paper has a genuinely interesting model and a real technical contribution, but the central claim—that the fractional Ising chain realizes a continuously tunable dynamical exponent z=q/2—is not supported as stated.\n\nThe model is a 1D transverse-field Ising chain whose couplings come from the Riesz fractional derivative, J(r)~r^{-(1+q)} with sign alternation for q>2. That kernel is a nice way to interpolate between short-range, power-law, and frustrated regimes, and the paper is the first to study the q>2 sign-alternating window. The q<2 segment is a known power-law result, but the fractional formulation plus the q>2 physics is new. The MPS work is serious: they fit the nonlocal kernel with a sum of exponentials to build an MPO, and they extract z two independent ways (finite-size gap scaling and bond-entropy wavefronts) that agree to within a few percent. That is real evidence of careful numerics, and the exponential-fit method is worth borrowing.\n\nThe problems are where the stress-test note points. The analytic step from (27) and (33) to (34) is not shown, and when I tried to close it, the q=2 limit does not come out right unless the hypergeometric subtraction in (38) is verified—and that subtraction is asserted, not proven. Since q=2 is the one limit where you know the answer, this is a serious gap in the derivation of the headline result.\n\nThe bigger issue is numerical. The paper itself cites Maghrebi et al. (Ref. [46]), who found in 1D that the mean-field relation z=q/2 holds only for q<2/3; for larger q fluctuations drive z toward 1. The numerics here show z tracking q/2 all the way to q~2, then reverting to z≈1 around q≳2.5. The authors call the deviation finite-size effects but do not test larger chains. At L=200, the observed curve is exactly the shape you would get from the RG crossover, so the “confirmation” over 0<q<2.5 is not established. And there is a terminology problem: for q<2, z=q/2<1 gives t^{1/z}=t^{2/q}, a superlinear cone, yet the paper calls it sublinear. That inversion will confuse readers and undermines confidence in the presentation.\n\nWho is this for? It is for people working on long-range quantum simulators and fractional generalizations of lattice models. It deserves a serious referee—the model and the q>2 physics are worth engaging with—but the paper needs major revision: verify the derivation, run longer chains, and either reconcile with the RG prediction or explicitly frame the observed z(q) as a finite-size crossover. I would not cite the central claim in its current form.","headline":"Genuinely novel model and solid numerics, but the central z=q/2 claim is undermined by an unverified derivation and an untested finite-size escape hatch.","tokens_in":13891,"tokens_out":4183,"would_cite":false,"duration_ms":44745,"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":"The paper claims that the dynamical critical exponent of the fractional multiscale transverse-field Ising model is z=q/2, giving a continuously tunable family of entanglement light cones.","keywords":["fractional derivative","Lévy flights","quantum Ising model","long-range interactions","dynamical critical exponent","light cone","matrix product states","entanglement entropy"],"falsifier":"Perform the same gap and wavefront analysis on chains several times longer than L=200 at fixed q in, say, 0.7<q<2.5: if the measured z bends toward 1 as L grows, the continuously tunable z=q/2 is a finite-size crossover, not the asymptotic behavior. A trapped-ion or Rydberg implementation that measures an entanglement-front exponent consistent with z=1 for q>2 at accessible times would similarly contradict the superlinear-cone claim.","tokens_in":12683,"feed_emoji":"⚛️","tokens_out":7986,"duration_ms":79742,"temperature":0.7,"pith_summary":"The paper claims that a one-dimensional transverse-field Ising chain with couplings generated by a fractional (Riesz) derivative has a dynamical critical exponent set directly by the fractional order q: z=q/2. If true, one control knob sweeps the model through three dynamical regimes—sublinear entanglement light cones for q<2, a faint superlinear cone with z<1 for 2<q<2.5, and the standard ballistic z=1 regime for q≳2.5. The claim matters because it ties quantum information propagation to Lévy-flight statistics and because the required fractional interactions are within reach of current trapped-ion, Rydberg, and cavity-mediated quantum simulators. The authors support it with a mean-field Jordan-Wigner calculation plus matrix-product-state time evolution on chains up to 200 sites.","feed_headline":"Fractional coupling sets spin-chain light-cone speed to z=q/2","feed_subtitle":"One parameter bridges sublinear, superlinear, and linear entanglement spreading in a fractional Ising model.","key_machinery":"The central object is the Riesz fractional derivative, a Fourier multiplier -|k|^q that discretizes into the sign-alternating binomial coupling J(r)=(-1)^{r+1} binom(q, q/2+r). The analytical argument runs through a truncated Jordan-Wigner transformation, a Fourier transform, and a Bogoliubov diagonalization; the crux is Eq. (30), which turns the cosine sum over the coupling into |2 sin(k/2)|^q, so the excitation gap inherits a |k|^{q/2} power law. On the numerical side, the workhorse is an exponential-sum decomposition of J(r) into roughly 10-14 terms, which keeps matrix-product-operator bond dimension small and makes time-dependent variational principle evolution with nonlocal couplings practical. Dynamical exponents are extracted two independent ways: finite-size gap scaling with a subleading correction term, and bond-entropy light-cone contours fit to |j-j*| ~ $t^{{1/z}}$.","core_discovery":"The authors establish that at the quantum critical point g_c = (1/2) binom(q, q/2), the Bogoliubov quasiparticle dispersion of the fractional multiscale Ising model scales as E_k ∝ |k|^{q/2} in the long-wavelength limit, yielding z=q/2. Their truncated Jordan-Wigner mapping replaces the nonlocal Jordan-Wigner strings by unity, Fourier transforms the bilinear Hamiltonian, and uses a closed-form hypergeometric identity to evaluate the coupling sums; an Euler transform then isolates the leading |k|^q term in the gap. The numerical simulations find a continuously tunable exponent over 0<q<2.5, with z close to q/2 up to q≈2, a superlinear cone for 2<q<2.5, and a return to z=1 for larger q. The paper explicitly notes that two-loop renormalization-group results predict z=q/2 only in the strictly long-range window q<2/3, and that the apparent agreement at larger q is hypothesized to come from finite-size effects.","pith_inferences":["A direct test of the finite-size hypothesis is to compute the gap and wavefront exponents on chains several times longer than 200 sites; if z drifts toward 1 for q>2/3 at larger L, the continuously tunable z becomes a pre-asymptotic crossover rather than the asymptotic universality class.","The correspondence with Lévy flights suggests a quantitative link between the entanglement-front exponent 1/z and the classical transport exponent 2/q; this could be tested by engineering a single-site perturbation in a trapped-ion chain and measuring the front shape, not just the scaling exponent.","Because the mean-field derivation relies on dropping the Jordan-Wigner strings, a non-perturbative check—such as a fermionic simulation that keeps string terms or an independent entanglement-scaling calculation—would isolate whether z=q/2 is an artifact of the truncation.","In d≥2 the paper expects the frustrated superlinear regime to be more pronounced; extending the exponential-sum matrix-product-operator method to two-dimensional cylinder geometries would turn that expectation into a concrete prediction."],"forward_implications":["For q<2, entanglement spreads sublinearly, so information propagates faster than in the nearest-neighbour model while still displaying a well-defined scaling front.","For 2<q<2.5, the sign-alternating frustrated tail produces a superlinear cone with z<1, outside the universality class of simple power-law models.","At q≳2.5 the model effectively becomes local again, recovering z=1 and linear light-cone propagation.","The same model reproduces the standard transverse-field Ising critical point at q=2 and a ballistic regime, making the fractional model a unifying interpolation rather than a separate theory.","If the mean-field exponent holds beyond numerical finite sizes, the model gives a one-parameter family of light cones that can be probed in trapped-ion or Rydberg simulators by engineering near-fractional kernels."],"supporting_citations":[{"why":"Supplies the closed-form identity that turns the binomial coupling sum into |2 sin(k/2)|^q, the step that gives the |k|^{q/2} dispersion.","marker":"[37]"},{"why":"Two-loop renormalization-group study of power-law long-range Ising chains; provides the q<2/3 baseline prediction that the paper's larger-q agreement is compared against.","marker":"[46]"},{"why":"Introduces the exponential-sum decomposition of long-range couplings that makes the matrix-product-operator representation and time-dependent variational principle evolution feasible.","marker":"[34]"},{"why":"Provides the ground-state optimization and energy-variance convergence used to locate the critical point and compute finite-size gaps.","marker":"[40]"},{"why":"Exact solution of the long-range Kitaev chain that supports the truncated Jordan-Wigner quadratic Hamiltonian.","marker":"[44]"},{"why":"Renormalization-group argument that higher-order Jordan-Wigner string terms are irrelevant, justifying the mean-field dynamical exponent.","marker":"[42]"},{"why":"Establishes the linear light-cone baseline for short-range models that the fractional model is shown to depart from.","marker":"[1]"}],"fun_headline_variants":["Fractional derivative order dials spin-chain light cone speed","Quantum light cone exponent set by fractional order: z=q/2","One parameter tunes entanglement speed: fractional Ising model","Sublinear to superlinear: fractional order controls spin light cone","Levy-like spreading in spin chains tied to fractional order q"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that chains of 200 sites already expose the thermodynamic-limit dynamical exponent; the paper's own Section 5 states that accepted two-loop theory predicts z=q/2 only for q<2/3 and that the observed agreement for larger q is hypothesized, not tested, to be a finite-size effect.","fun_headline_variants_meta":{"raw":{"variants":["Fractional derivative order dials spin-chain light cone speed","Quantum light cone exponent set by fractional order: z=q/2","One parameter tunes entanglement speed: fractional Ising model","Sublinear to superlinear: fractional order controls spin light cone","Levy-like spreading in spin chains tied to fractional order q"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000998,"raw_usage":{"total_tokens":4314,"prompt_tokens":1125,"completion_tokens":3189,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":741,"completion_tokens_details":{"reasoning_tokens":3105}},"tokens_in":741,"tokens_out":3189,"duration_ms":26882,"temperature":1.0,"reasoning_tokens":3105,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:02:02.394400+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform the same gap and wavefront analysis on chains several times longer than L=200 at fixed q in, say, 0.7<q<2.5: if the measured z bends toward 1 as L grows, the continuously tunable z=q/2 is a finite-size crossover, not the asymptotic behavior. A trapped-ion or Rydberg implementation that measures an entanglement-front exponent consistent with z=1 for q>2 at accessible times would similarly contradict the superlinear-cone claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the closed-form identity that turns the binomial coupling sum into |2 sin(k/2)|^q, the step that gives the |k|^{q/2} dispersion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Two-loop renormalization-group study of power-law long-range Ising chains; provides the q<2/3 baseline prediction that the paper's larger-q agreement is compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the exponential-sum decomposition of long-range couplings that makes the matrix-product-operator representation and time-dependent variational principle evolution feasible."},{"cited_title":"Classical and Quantum Phase Transitions in Multiscale Media: Universality and Critical Exponents in the Fractional Ising Model","cited_arxiv_id":"2501.14134","evidence_quote":"Provides the ground-state optimization and energy-variance convergence used to locate the critical point and compute finite-size gaps."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Exact solution of the long-range Kitaev chain that supports the truncated Jordan-Wigner quadratic Hamiltonian."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Renormalization-group argument that higher-order Jordan-Wigner string terms are irrelevant, justifying the mean-field dynamical exponent."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the linear light-cone baseline for short-range models that the fractional model is shown to depart from."}],"review_version":1}