{"id":"0ff0c691-0504-4154-9160-8aa21d324049","arxiv_id":"2412.01571","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In a 1D Bose-Hubbard model with hopping 1/r^α, the superfluid-Mott-insulator transition is continuous and scale-invariant for 1<α≤3, not BKT, while true long-range order in the superfluid appears only for α≤2.","lead":"This paper computes the zero-temperature phase diagram of a 1D chain of interacting bosons whose hopping amplitude decays with distance as 1/r^α, using exact large-scale quantum Monte Carlo simulations. It reports that for 1<α≤3 the superfluid-to-Mott-insulator transition is continuous and scale-invariant, belonging to a new universality class distinct from the Berezinskii-Kosterlitz-Thouless transition, and that the superfluid has long-range order only for α≤2.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Crossing of ⟨W²⟩ does not rule out BKT: BKT exhibits finite-size crossings due to the universal stiffness jump, so the paper's central inference is unsupported without a BKT baseline.","rationale":"The reader identified the lack of a BKT baseline as the weakest assumption, and I agree that a BKT comparison is essential. However, the specific flaw is more pointed: the paper's stated criterion for ruling out BKT—the existence of a crossing in ⟨W²⟩ vs t/U—is not a valid discriminator, because such crossings are a known finite-size feature of BKT transitions. This makes the central argument weaker than the reader's phrasing suggests, although the conclusion may still be true. The proposed test is decisive: applying the same analysis to the short-range BKT model would show whether the crossing criterion misclassifies a known BKT transition. If it does, the paper must switch to a drift-based or explicit BKT-fit analysis before its non-BKT claim can be accepted. The paper has substantial numerical evidence and is honest about limitations, so a conditional verdict with a demand for this BKT control is appropriate rather than outright rejection.","tokens_in":14133,"tokens_out":15813,"duration_ms":140975,"concrete_test":"Run the same QMC protocol (worm algorithm, ⟨W²⟩ measurement, L = 64–512) on the nearest-neighbor 1D Bose-Hubbard model (α=∞), whose SF-MI transition is known to be BKT, using β = L (z=1). If the ⟨W²⟩ vs t/U curves exhibit crossings similar to Fig. 3, then the presence of a crossing cannot be used to rule out BKT, directly testing the paper's central criterion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the SF-MI transition is non-BKT for 1<α≤3 rests on the statement in Section III (Fig. 3) that 'the very presence of a crossing in the ⟨W²⟩−(t/U) curves rules out the BKT universality class.' This inference is not valid. In BKT transitions, the superfluid stiffness (or the related winding-number variance) is expected to exhibit finite-size crossings near the critical point: the universal stiffness jump implies that curves for different L intersect at (or close to) the transition, and such crossings are routinely used to locate BKT transitions in the 2D XY model and in the short-range 1D Bose-Hubbard model. The correct discriminator between BKT and a power-law scale-invariant transition is not the existence of a crossing but the drift of the crossing point with L (logarithmic in BKT; power-law in a conventional transition) and/or a quantitative fit to the BKT scaling form, neither of which is provided. The paper's power-law collapse in Fig. 4 assumes the ansatz it is trying to prove; a BKT-like essential-singularity scaling could also produce an acceptable collapse over the limited L range (64–512). Additionally, the simulations use β = L^{z*} with z* = (α−1)/2 < 1, whereas a BKT transition would require β ∝ L (z=1); the smaller β may induce thermal rounding that produces crossings. Thus the evidence is insufficient to establish incompatibility with BKT.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports quantum Monte Carlo simulations (worm algorithm) of the one-dimensional Bose-Hubbard model with power-law hopping 1/r^α. For 1<α≤3 the authors claim that the superfluid-to-Mott-insulator transition at unit filling is continuous and scale invariant, with a new universality class, and that the BKT scenario is recovered only for α>3. The evidence consists of crossings of winding-number fluctuations, power-law finite-size collapses of the superfluid stiffness, an excitation spectrum E(k)∼k^{(α−1)/2}, and an analysis of single-particle correlations that is interpreted as showing long-range order for α≤2 and 'anomalous quasi-long-range order' for 2<α≤3. A grand-canonical phase diagram and correlation data in the Mott phase are also presented.","tokens_in":14386,"tokens_out":4401,"duration_ms":41552,"significance":"If substantiated, the claim of a non-BKT scale-invariant universality class for 1<α≤3 would revise the standard picture that all localization transitions in one-dimensional bosonic systems belong to the BKT class, and it would provide a benchmark for experiments with dipolar atoms, molecules, and ion chains. The paper brings large-scale QMC data (L up to 1024 for correlations, L up to 512 for critical scaling) and uses standard tools such as the worm algorithm and the Kawashima-Ito-Houdayer-Hartmann collapse quality metric. The agreement between the gap-derived ν and the collapse-derived ν for two representative α values is a genuine consistency check. However, the central discriminator against BKT is currently a single invalid logical step, and the supporting collapses use only three system sizes without a BKT baseline, so the significance of the claimed universality class is not yet established.","major_comments":[{"comment":"The statement that 'the very presence of a crossing in the ⟨W²⟩−(t/U) curves rules out the Berezinskii-Kosterlitz-Thouless universality class' is not valid. In the BKT scenario finite-size curves of the superfluid stiffness (or winding-number fluctuations) are expected to cross near the transition because of the universal stiffness jump, and such crossings are routinely used to locate BKT transitions in the 2D XY model and in the short-range 1D Bose-Hubbard model. The correct discriminator is the size dependence of the crossing point (logarithmic drift for BKT versus power-law for a conventional transition) and/or a quantitative fit of the data to a BKT scaling form; neither is provided. Since this inference is the main basis for the non-BKT central claim, it must be replaced by a concrete test.","section":"Section III, Fig. 3"},{"comment":"The power-law data collapse uses only L = 128, 256, 512 and assumes the power-law scaling form that is being tested. An essential-singularity BKT form can produce an apparently acceptable collapse over such a modest size range, especially with the freedom of three fitting parameters (ν, ζ, and t/Uc). The manuscript does not compare the collapse quality with a BKT baseline or report the drift of the crossing in the insets of Fig. 3 in a way that distinguishes logarithmic from power-law behavior. The claim of incompatibility with BKT therefore lacks the needed statistical and model-selection support.","section":"Section III, Fig. 4"},{"comment":"The simulations use inverse temperature β = L^{z*} with z* = (α−1)/2 < 1, whereas the BKT scenario one wishes to exclude has z = 1. In the competing BKT interpretation, β/L → 0 as L → ∞, so the largest simulated sizes would fall outside the genuine ground-state scaling regime and could display spurious crossings from thermal rounding. The paper needs to justify the choice β = L^{z*} by checking convergence in β for at least one α, or by performing a consistency run with β ∝ L.","section":"Section III, simulation parameters"},{"comment":"The ordering-regime claim for 2<α≤3 ('anomalous quasi-long-range order') is based on a comparison of two fits that the Supplemental Material states are 'equally effective' in describing G(ℓ). The pure power-law fit is then preferred because it 'eliminates one fitting parameter,' but this is not a valid model-selection criterion when the constant-plus-power-law form is the physically motivated expression for a long-range-ordered state. The finite-size extrapolation of the decay exponent γ in Fig. 5(e) assumes a linear 1/L dependence without independent justification. Since the abstract presents a sequence of ordering regimes as a key result, this part of the evidence needs either a stronger statistical discrimination between the two forms or a more explicit discussion of the ambiguity.","section":"Section III, Fig. 5 and Supplemental Material IB"}],"minor_comments":[{"comment":"The main text says the single-particle density matrix is analyzed 'near the critical point' but then specifies hard-core bosons at half-filling, which do not have a Mott transition; the relationship between these parameter choices should be clarified.","section":"Section III (text near Fig. 5)"},{"comment":"The statement that the pure power-law fit is 'better by eliminating one fitting parameter' contradicts the earlier sentence in the same section that both fits are equally effective; this inconsistency should be resolved.","section":"Supplemental Material, Sec. IB"},{"comment":"The notation for the transition boundary uses t/Uc, t/U_c, and t/Uc interchangeably; please standardize the spacing and subscript format in equations and figure captions.","section":"Throughout"},{"comment":"The value k=π/64 used for the Green function in the inset is stated in the caption but not in the main text; please define the momentum grid and clarify that the dispersion is evaluated at the smallest nonzero momentum.","section":"Fig. 2 inset"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's central claim is timely and would be important if true, but the current evidence does not yet support the specific claim of incompatibility with BKT. The crossing argument is logically flawed as presented, and the finite-size collapses are too limited to exclude BKT on their own. I would encourage the authors to add a direct BKT scaling fit (or an analysis of the crossing-point drift) and to address the β = L^{z*} issue. The ordering-regime classification also needs a more principled model-selection approach. I have no concerns about the authors' integrity; the issues are purely technical and appear fixable within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a serious QMC study with a big claim, and the claim is not yet supported. Send it out, but tell the authors the central discriminator needs to change.\n\nWhat's new: for the 1D Bose-Hubbard model with 1/r^α hopping, they present a phase diagram where for 1<α≤3 the constant-density SF-MI transition looks continuous with power-law scaling, and the low-energy dispersion follows E(k)~k^{(α-1)/2} with z*=(α-1)/2, matching the spin-wave prediction. They also see no long-range order for α>2, contradicting the DMRG estimate of α_c≈2.8. The QMC work is large scale: L up to 512 for stiffness, 1024 for correlations, worm algorithm, and the ν from collapse and from the gap agree. That is real substance.\n\nThe soft spot is the logic that identifies the universality class. The paper says the crossing of ⟨W²⟩ rules out BKT. That's not right. BKT also produces finite-size crossings from the universal stiffness jump; what distinguishes BKT is the logarithmic drift of the crossing and the essential-singularity form of the scaling, and neither is fitted. Their power-law collapse assumes the power-law ansatz, and with three sizes (128, 256, 512) a BKT form could plausibly collapse too. They also run β=L^{z*} with z*=(α-1)/2; if the true transition were BKT with z=1, that β would be too small and could induce thermal rounding, producing fake crossings. So the central claim is under-supported, not wrong.\n\nThe ordering-regime part is shakier. For 2<α≤3 they label 'anomalous quasi-long-range order' because a pure power-law fit is more stable than a constant+power-law fit. Both fit equally well, and preferring the simpler form is a modeling choice, not evidence. The γ→0 extrapolation for α≲2 is fine, but the anomalous QLRO category needs more support.\n\nWho's it for: people working on long-range bosons, trapped ions, dipolar gases, and anyone doing QMC scaling. The paper is a useful benchmark even if the universality-class claim is unresolved. I'd send it to peer review and ask for an explicit BKT fit and a crossing-drift analysis; that's a normal request, not a rejection.","headline":"Plausible but unproven: the non-BKT claim rests on a crossing that BKT also produces.","tokens_in":15034,"tokens_out":3435,"would_cite":false,"duration_ms":30753,"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":"For 1<α≤3, the superfluid–Mott-insulator transition of a one-dimensional Bose-Hubbard chain with power-law hopping is continuous and scale invariant, not BKT.","keywords":["Bose-Hubbard model","power-law hopping","Berezinskii-Kosterlitz-Thouless transition","superfluid-Mott insulator transition","quantum Monte Carlo","winding number","long-range order","universality class"],"falsifier":"Fit the same winding-number and superfluid-stiffness data with the BKT essential-singularity form, e.g., ξ∼exp(a/√δ) with a universal stiffness jump; if a good collapse is obtained with the available system sizes, the reported ν and z would be effective finite-size values rather than evidence for a new universality class.","tokens_in":13838,"feed_emoji":"⚛️","tokens_out":4587,"duration_ms":36928,"temperature":0.7,"pith_summary":"The paper asks what happens to the superfluid-to-Mott-insulator transition in a one-dimensional lattice gas when hopping decays as 1/r^α rather than only between neighbors. Using large-scale quantum Monte Carlo simulations, it argues that for every decay exponent 1<α≤3 the transition is continuous and scale invariant, governed by a new universality class with correlation-length exponent ν(α) and dynamical exponent z=(α−1)/2, instead of the Berezinskii-Kosterlitz-Thouless (BKT) vortex-unbinding scenario that holds for α>3. If correct, this means power-law hopping changes the nature of localization in one dimension and fixes the long-range threshold at α*=3. The same simulations show the superfluid ground state evolves from true long-range order (α≤2) to an anomalous quasi-long-range order (2<α≤3).","feed_headline":"Power-law hopping flips 1D boson transition out of BKT class","feed_subtitle":"Large-scale QMC finds a new universality class for 1<α≤3, with BKT recovered only above α=3.","key_machinery":"The load-bearing observable is the mean-square winding number ⟨$W^{2}$⟩, a scale-invariant proxy for the superfluid stiffness Y_s; a crossing of ⟨$W^{2}$⟩ curves for increasing system size identifies a continuous scale-invariant transition and directly contradicts BKT's essential-singularity scaling. Critical exponents are obtained by collapsing $L^{{−ζ/ν}}$Y_s against $L^{{1/ν}}$(t/U−(t/U)_c), using the Nelder-Mead optimization of the Kawashima-Ito-Houdayer-Hartmann quality metric, and z*=(α−1)/2 comes from fitting the low-energy Green function G(k,τ). The single-particle density matrix G(ℓ) measured on chord distance c(ℓ)=sin(πℓ/L) is used to distinguish long-range order from quasi-long-range order.","core_discovery":"For the one-dimensional Bose-Hubbard model with hopping amplitude t/|r|^α, the zero-temperature superfluid–Mott-insulator transition at unit filling is claimed to be a continuous, scale-invariant quantum phase transition for all 1<α≤3, incompatible with the BKT universality class that is recovered only for α>3. The evidence is a crossing of winding-number fluctuations with system size, power-law data collapse of the superfluid stiffness, and a sublinear dispersion E(k)∼$k^{{(α−1)/2}}$ matching spin-wave theory. The authors extract the correlation-length exponent ν from data collapse and from the gap Δ∼|t/U−(t/U)_c|^{z*ν}, and classify superfluid correlations: true long-range order for α≤2, anomalous quasi-long-range order for 2<α≤3, and conventional algebraic decay for α>3.","pith_inferences":["A direct BKT fit to the same finite-size data would sharpen the claim: if BKT scaling with logarithmic corrections also collapses the data, the reported ν would be an effective exponent rather than a new universality class.","The threshold α*=3 may be generic for power-law hoppings in one dimension, suggesting neighboring models such as long-range interacting spin chains should be re-examined for the same crossover.","The 'anomalous quasi-long-range order' label for 2<α≤3 is based on the more stable power-law fit; a larger-L study of the constant-plus-power-law alternative would decide whether a tiny condensate fraction survives.","Measuring the single-particle Green function in an ion chain with tunable α could test the dispersion prediction without needing to approach the transition."],"forward_implications":["For every 1<α≤3, the transition has a finite correlation-length exponent and power-law scaling, so standard BKT descriptions of one-dimensional bosons do not apply in this regime.","The critical hopping ratio (t/U)_c decreases as α decreases, and the Mott lobe shape changes from rounded (α<2) to pointed (α>2), matching the crossover seen in correlations.","The superfluid phase has true long-range order for α≤2 and no true long-range order for α>2; the region 2<α≤3 is identified as anomalous quasi-long-range order.","In the Mott phase, the single-particle density matrix decays as ℓ^{−α}, tying the gap to the power-law hopping exponent.","The predicted dispersion E(k)∼k^{(α−1)/2} and correlation functions are directly accessible in trapped-ion and dipolar cold-molecule experiments."],"supporting_citations":[{"why":"Provides the worm-algorithm quantum Monte Carlo method used to compute winding numbers, Green functions, and density matrices.","marker":"[28]"},{"why":"Supplies the spin-wave prediction E(k)∼k^{(α−1)/2} that the measured spectrum is compared with.","marker":"[34]"},{"why":"Gives the earlier density-matrix-renormalization-group estimate α_c≈2.8 for long-range order, which this work revises to α*=3.","marker":"[35]"},{"why":"Defines the BKT universality class for short-range one-dimensional bosons that the paper claims is excluded for α≤3.","marker":"[10]"},{"why":"Provides the nearest-neighbor critical value (t/U)_c=0.300±0.025 used as the short-range reference.","marker":"[33]"},{"why":"Supplies predictions for the long-range XY limit against which the hard-core-boson density matrices are compared.","marker":"[44]"},{"why":"Sets the weak-long-range regime d<α<α* and the threshold concept that the paper sharpens to α*=3.","marker":"[18]"}],"fun_headline_variants":["Power-law hopping rewrites 1D boson transition universality","New quantum critical class for 1D bosons with 1/r^α hopping","1D Bose-Hubbard: power-law hopping defeats BKT scaling","Long-range hopping yields continuous SF-MI transition in 1D","QMC finds distinct universality class in power-law 1D bosons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that a new universality class exists rests on finite-size data up to 512 sites being large enough to rule out BKT scaling, whose logarithmic corrections can mimic power-law crossings over this range.","fun_headline_variants_meta":{"raw":{"variants":["Power-law hopping rewrites 1D boson transition universality","New quantum critical class for 1D bosons with 1/r^α hopping","1D Bose-Hubbard: power-law hopping defeats BKT scaling","Long-range hopping yields continuous SF-MI transition in 1D","QMC finds distinct universality class in power-law 1D bosons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000325,"raw_usage":{"total_tokens":1825,"prompt_tokens":951,"completion_tokens":874,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":776}},"tokens_in":567,"tokens_out":874,"duration_ms":7344,"temperature":1.0,"reasoning_tokens":776,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:18:12.338630+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the same winding-number and superfluid-stiffness data with the BKT essential-singularity form, e.g., ξ∼exp(a/√δ) with a universal stiffness jump; if a good collapse is obtained with the available system sizes, the reported ν and z would be effective finite-size values rather than evidence for a new universality class.","supporting_citations":[{"cited_title":"Lepori and L","cited_arxiv_id":null,"evidence_quote":"Provides the worm-algorithm quantum Monte Carlo method used to compute winding numbers, Green functions, and density matrices."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the spin-wave prediction E(k)∼k^{(α−1)/2} that the measured spectrum is compared with."},{"cited_title":"Fr´ erot, P","cited_arxiv_id":null,"evidence_quote":"Gives the earlier density-matrix-renormalization-group estimate α_c≈2.8 for long-range order, which this work revises to α*=3."},{"cited_title":"Kashurnikov, A","cited_arxiv_id":null,"evidence_quote":"Provides the nearest-neighbor critical value (t/U)_c=0.300±0.025 used as the short-range reference."}],"review_version":1}