{"id":"f34ed96d-a4d3-4d7a-9dfa-15a44b11b72c","arxiv_id":"2512.24557","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a 2D spin-1 XY model with 1/r^3 interactions, the Higgs mode disperses linearly and its Beliaev damping is strongly suppressed near the quantum critical point, making it observable.","lead":"The authors show that in a two-dimensional spin-1 magnet with 1/r^3 interactions, the Higgs amplitude mode can stay long-lived instead of being quickly damped. This suggests Rydberg-atom arrays are a promising platform for observing this elusive collective excitation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"One-loop HP truncation not justified near the 2D critical point; thermal occupation of soft NG modes may invalidate the claimed suppression of Higgs damping.","rationale":"The reader's weakest_assumption correctly identifies the one-loop truncation as the main soft spot. My analysis confirms that the low-energy phase space and the divergent thermal occupation of the NG modes are precisely the conditions under which higher-order diagrams could become important. I also note supplementary issues—the missing real-part/residue calculation—but these are secondary. The proposed two-loop/A→3Φ calculation is a direct and quantitative test of the truncation. If the correction is small, the paper's conclusion stands; if not, the conditional verdict should be downgraded. Thus my assessment does not change the reader's CONDITIONAL verdict, but it sharpens the required check.","tokens_in":13709,"tokens_out":43434,"duration_ms":373782,"concrete_test":"Derive the fourth-order Holstein-Primakoff Hamiltonian for α=3,d=2 and evaluate the A→3Φ decay rate and the two-loop self-energy for u=0.95 (Δ_h≈1.8 J z_eff) at T=0.5 J z_eff, using the same long-wavelength approximations as in Sec. IV.B. Compare the resulting total Γ/Δ_h with the one-loop value from Eq. (41). If the higher-order contribution is ≲10% of the one-loop result, the suppression scenario is robust; if it is comparable or larger, the claim of a long-lived Higgs mode is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that Γ/Δ_h remains finite and small as u→u_c (Eqs. (35)-(41))—rests on a one-loop Beliaev calculation in which the Holstein-Primakoff expansion is truncated at cubic order (Sec. III.C) and the self-energy is evaluated only from the A→ΦΦ vertex (Eq. (32)). No estimate is given for the neglected quartic terms or higher-loop diagrams. For d=2, α=3, the on-shell NG momentum scales as k0∝Δ_h^2∝(1-u^2), and the NG Bogoliubov coefficient grows as u_Φ(k0)∝k0^{-1/4}. The vertex M_{0,k0,-k0} then scales as (1-u^2)^{-1/4}, while the thermal factor in Eq. (33) is (1+2n_B)≈4T/Δ_h, which diverges as u→u_c. This large occupation of very soft NG modes enhances higher-order diagrams containing additional NG lines—e.g., the A→3Φ process from fourth-order HP terms or the two-loop sunset with two internal Φ propagators—by additional powers of T/Δ_h. The paper does not compute the quartic Hamiltonian, the real part of the self-energy, or the quasiparticle residue Z. A finite Γ/Δ_h is necessary but not sufficient for a long-lived mode: if Z→0 or if two-loop contributions dominate, the spectral function will not exhibit a sharp Higgs peak. Without a quantitative bound on these neglected terms, the conclusion that the long-range interaction makes the Higgs mode 'long-lived' is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the spin-1 XY model with a quadratic Zeeman term and power-law interactions, Eq. (1), in the XY-ferromagnetic phase near the transition to the disordered phase. Using a mean-field product state, a Holstein-Primakoff expansion truncated at cubic order, and finite-temperature Green's functions, the authors derive analytic dispersions for the Higgs and NG modes, Eqs. (24)-(25), and a one-loop Beliaev damping rate for the k=0 Higgs mode, Eqs. (33), (39), and (41). For d=2 and alpha=3, the parameter set relevant to Rydberg-atom arrays, they find a linear Higgs dispersion, a square-root NG dispersion, and a damping-to-gap ratio Gamma/Delta_h that remains finite as u->u_c, in contrast to the short-range case where it diverges. The alpha->5, d=3 limit reproduces the previously known short-range result, Eq. (42). The paper also proposes a concrete experimental protocol for exciting and detecting the Higgs mode.","tokens_in":14094,"tokens_out":18333,"duration_ms":165266,"significance":"If the one-loop calculation is quantitatively reliable, the paper gives an experimentally relevant prediction for Rydberg-atom arrays: the long-range nature of the interaction suppresses Higgs-mode damping near the critical point, making a long-lived amplitude mode observable. The analytic closed forms are a strength, and the reduction to the known short-range result is a useful consistency check. However, the central conclusion depends on an uncontrolled cubic-order, one-loop truncation near the critical point and on the numerical prefactor of Eq. (41). The significance is therefore contingent on the robustness of that approximation, which the manuscript does not currently establish.","major_comments":[{"comment":"The central claim that Gamma/Delta_h stays finite and small as u->u_c rests on a Holstein-Primakoff expansion truncated at cubic order and a one-loop self-energy containing only the A->Phi Phi process, Eq. (32). Near the critical point the on-shell NG momentum is k0 proportional to Delta_h^2 proportional to (1-u^2), the NG Bogoliubov coefficient grows as k0^{-1/4}, and the thermal factor in Eq. (33) behaves as 4T/Delta_h. The manuscript gives no estimate of the neglected quartic HP terms, of the two-loop sunset diagram with two Phi propagators, or of the quasiparticle residue Z. A finite one-loop Gamma/Delta_h is not by itself sufficient to guarantee a sharp Higgs peak if higher-order diagrams or a vanishing Z modify the spectral function. I ask for a quantitative scaling estimate, or an explicit calculation, of the leading neglected corrections in the limit u->1.","section":"Sec. III.C and III.E, Eq. (41)"},{"comment":"The step from the exact self-energy expression, Eq. (33), to the closed form, Eq. (39), is not shown. In particular, the simplification of |M_{0,k,-k}|^2 in the long-wavelength limit and the angular integration are not displayed, so the numerical prefactor of Eq. (41) cannot be checked from the text. This prefactor controls the statement that the damping is 'sufficiently smaller than unity' near the critical point. Please provide the intermediate algebra or include it in a supplemental derivation.","section":"Sec. IV.B, Eq. (39)"},{"comment":"The long-wavelength approximations for the Bogoliubov coefficients do not match the exact small-k limits of Eqs. (18)-(21). For u->1, Eq. (18) at k=0 gives u_A^2 approximately [4 sqrt(1-u^2)]^{-1}, whereas the text's u_A^2 = 1/(2 Delta_bar_h) = [8 sqrt(1-u^2)]^{-1}. Similarly, Eq. (20) with gamma_k approximately 1 - A|k| gives u_Phi^2 approximately [4 sqrt(A|k|)]^{-1}, whereas the text gives [8 sqrt(A|k|)]^{-1}. These factor-of-two discrepancies in u_A^2 and u_Phi^2 change |M_{0,k,-k}|^2 and hence the prefactor of Eq. (41). Please correct the approximations or explain the convention used.","section":"Sec. IV.B, approximations after Eq. (38)"}],"minor_comments":[{"comment":"The Fourier transform convention is not stated. The sign convention in Eq. (14) affects the momentum-conservation delta functions in Appendix A; please specify it explicitly.","section":"Eq. (14)"},{"comment":"The expansion gamma_k = 1 - A|k|^{alpha-d} is used for d < alpha < d+2. The limit alpha -> d+2, mentioned after Eq. (38), should be qualified: at alpha = d+2 logarithmic corrections can appear, and the simple power-law form requires care.","section":"Eq. (34)"},{"comment":"The statement that Eq. (42) agrees with previous works would be easier to verify if the corresponding equations in Refs. [22,29] were cited explicitly, together with the mapping z_eff -> z.","section":"Sec. IV.B, Eq. (42)"},{"comment":"The concluding paragraph mentions a possible future study of spatial inhomogeneity and Higgs bound states; this is fine, but the discussion would benefit from a sentence distinguishing mean-field u_c=1 from the numerically more accurate u_c values cited in Sec. II.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The qualitative prediction is interesting and the consistency check with the short-range limit is valuable. However, the missing derivation of Eq. (39), the apparent inconsistency in the long-wavelength Bogoliubov coefficients in Sec. IV.B, and the unquantified one-loop/truncation error near the critical point are load-bearing for the central claim. I recommend major revision rather than rejection: the authors should be asked to supply the omitted algebra, correct or clarify the asymptotic coefficients, and provide a scaling estimate of neglected higher-order corrections."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe paper does exactly what the abstract says: it takes the known HP/Bogoliubov/one-loop treatment of the spin-1 XY model and swaps the nearest-neighbor dispersion for a long-range one, γ_k ≈ 1 - A|k|^{α-d}. For d=2, α=3 the resulting dispersions are new: linear Higgs, sqrt NG. The Beliaev damping formula, Eq. (39), with the concrete d=2, α=3 limit Eq. (41) is also new. The consistency check with α→5, d=3 reproduces the previous nearest-neighbor result, which is a good sign that the algebra is right.\n\nThe physics point is genuinely interesting. In the short-range 2D case the Higgs mode is overdamped because Γ/Δ_h diverges as u→u_c. Here the one-loop calculation gives Γ/Δ_h → finite ∝ T, so the ratio stays below unity at low T. That is a concrete, falsifiable prediction for Rydberg atom arrays, and the experimental protocol in Sec. IV B is sensible.\n\nNow the soft spots, in order of size. First, the paper does not show the integration from Eq. (33) to Eq. (39). The prefactor of Eq. (41) is therefore not checkable from the text. That should have been an appendix line at minimum. Second, the Holstein-Primakoff expansion is truncated at cubic order and the self-energy is one-loop. The paper does not address the validity of this truncation near the 2D critical point. The stress-test concern about soft NG modes carrying large thermal occupation is real: on-shell k0 ~ Δ_h^2, u_Φ(k0) ~ (1-u^2)^{-1/4}, and n_B ~ T/Δ_h. Higher-order diagrams with additional NG lines could in principle contribute powers of T/Δ_h and change the ratio Γ/Δ_h. The paper doesn't estimate those. This is not a fatal flaw — the approximation is standard and the short-range limit checks out — but it means the 'long-lived' claim is established only within the one-loop approximation. Third, the mean-field treatment of the critical point ignores the known deviation of the critical exponent from 1/2, though for long-range interactions that may be less severe.\n\nWho should read this: anyone working on collective modes in long-range spin systems or Rydberg arrays. It's a solid contribution, but the missing derivation and the unquantified truncation mean it deserves a serious referee rather than immediate acceptance. I would send it to peer review.","headline":"A clean extension of the known HP/Bogoliubov machinery to long-range interactions with a concrete, checkable prediction for Rydberg arrays; the one-loop truncation and a missing derivation keep the central claim from being fully established.","tokens_in":14624,"tokens_out":4750,"would_cite":true,"duration_ms":36562,"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-dimensional spin-1 XY model with dipolar (1/r^3) interactions, the Higgs mode acquires a linear dispersion and its Beliaev damping is strongly suppressed, so the amplitude mode can be long-lived at finite temperatures near the quan","keywords":["spin-1 XY model","long-range interactions","Higgs mode","Nambu-Goldstone mode","Beliaev damping","Rydberg atom arrays","quantum phase transition","finite-temperature Green's function"],"falsifier":"A direct measurement of the zero-momentum Higgs mode's damping rate as a function of temperature and u in a Rydberg-atom array with three states, using the proposed quench protocol, would test the predicted finite Γ/Δ_h as u→1. If the damping-to-gap ratio diverges or the Higgs oscillation is not resolved, the one-loop truncation would be falsified. Alternatively, a numerically exact quantum Monte Carlo calculation of the spectral function near the critical point could check whether corrections beyond one loop alter the suppression.","tokens_in":13589,"feed_emoji":"⚛️","tokens_out":3486,"duration_ms":31313,"temperature":0.7,"pith_summary":"This paper analyzes collective excitations in a spin-1 XY model on a square lattice with long-range interactions decaying as 1/r^3, as realized in Rydberg-atom arrays. For two dimensions, it finds that the Higgs mode—a gapped amplitude fluctuation of the magnetic order—has a linear dispersion, while the Nambu-Goldstone mode has a square-root dispersion. The central result is that the damping rate of the zero-momentum Higgs mode, evaluated via the Beliaev process of decaying into two NG modes, is strongly suppressed near the quantum critical point, so the damping-to-gap ratio stays finite at low temperatures. This suggests the Higgs mode can be observed as a long-lived coherent oscillation in Rydberg experiments, in contrast to short-range interacting systems where it is overdamped.","feed_headline":"Long-range interactions make the Higgs mode long-lived","feed_subtitle":"In a 2D spin-1 XY model with 1/r^3 interactions, the Higgs mode stays underdamped near the quantum transition, a key signature for Rydberg a","key_machinery":"The Holstein-Primakoff expansion truncated at cubic order, combined with a Bogoliubov transformation and one-loop finite-temperature Green's functions. The cubic term captures the Beliaev decay of a Higgs mode into two NG modes. The crucial small-momentum expansion γ_k ≈ 1 - A|k|^{α-d} (valid for d<α<d+2) converts the lattice dispersions into anomalous power laws, and its exponent α-d controls the momentum-space integrals in the damping formula.","core_discovery":"The paper claims that in the XY-ferromagnetic ordered phase of a spin-1 XY model with a quadratic Zeeman term and long-range interactions, the dispersion relations of the Higgs and Nambu-Goldstone modes are dramatically modified when the interaction decays as 1/r^3 in two dimensions: the Higgs mode becomes linear in momentum, E_A,k ≈ Δ_h + (g_h^2/2Δ_h)|k|, and the NG mode becomes square-root-like, E_Φ,k ∝ |k|^{1/2}. Using a finite-temperature Green's-function formalism with a one-loop self-energy, the Beliaev damping of the k=0 Higgs mode is derived analytically, yielding Γ_A,0 ∝ (1-u^2) coth(βΔ_h/4) for d=2, α=3. Consequently, Γ_A,0/Δ_h converges to a finite value proportional to temperatur","pith_inferences":["The suppression of damping likely persists beyond one loop as long as the A→ΦΦ channel remains the dominant decay process, but near the critical point higher-order corrections could alter the quantitative prediction; a two-loop calculation or a numerical simulation of the spectral function would test this.","The square-root NG dispersion implies a modified density of states at low energies, which could affect spin transport or thermalization dynamics in a measurable way, perhaps through spin-diffusion experiments.","A natural extension is to include spatial inhomogeneity or finite system sizes, where the Higgs mode may hybridize with edge states or form bound states; the paper mentions this direction in its conclusion.","Because the ratio Γ/Δ_h becomes temperature-independent near u_c, a single experimental measurement of the oscillation decay over a range of temperatures could extract both the gap and the damping strength simultaneously."],"forward_implications":["The Higgs mode in Rydberg-atom arrays with dipole-dipole interactions should be observable as a long-lived oscillation after a sudden quench of the quadratic Zeeman field, with a finite damping-to-gap ratio at low temperatures.","The linear Higgs dispersion and square-root NG dispersion are distinctive signatures that can be probed via momentum-resolved spectroscopy or Bragg scattering, distinguishing long-range from short-range interactions.","The damping rate vanishes at the critical point u→1 for all temperatures within the one-loop approximation, implying undamped amplitude oscillations exactly at the transition.","The results interpolate continuously to the nearest-neighbor interaction limit as α→d+2, recovering the standard quadratic Higgs and linear NG dispersions.","The proposed excitation protocol (adiabatic preparation followed by a sudden quench) provides a concrete experimental route to test the theory."],"fun_headline_variants":["Long-range interactions suppress Higgs damping in spin-1 XY model","Linear Higgs dispersion and square-root NG mode in 2D long-range spin model","Rydberg atoms: long-range coupling tunes Higgs and NG modes","In 2D XY ferromagnet, Higgs mode goes linear, damping drops","Long-range forces reshape quantum excitations in spin model"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The calculation assumes that the Holstein-Primakoff expansion truncated at cubic order and the one-loop self-energy (retaining only the A→ΦΦ decay channel) remain accurate even close to the quantum critical point, where quantum fluctuations grow and higher-order terms or alternative decay channels could become significant.","fun_headline_variants_meta":{"raw":{"variants":["Long-range interactions suppress Higgs damping in spin-1 XY model","Linear Higgs dispersion and square-root NG mode in 2D long-range spin model","Rydberg atoms: long-range coupling tunes Higgs and NG modes","In 2D XY ferromagnet, Higgs mode goes linear, damping drops","Long-range forces reshape quantum excitations in spin model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000873,"raw_usage":{"total_tokens":3628,"prompt_tokens":768,"completion_tokens":2860,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":512,"completion_tokens_details":{"reasoning_tokens":2767}},"tokens_in":512,"tokens_out":2860,"duration_ms":18372,"temperature":1.0,"reasoning_tokens":2767,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T13:20:06.419130+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of the zero-momentum Higgs mode's damping rate as a function of temperature and u in a Rydberg-atom array with three states, using the proposed quench protocol, would test the predicted finite Γ/Δ_h as u→1. If the damping-to-gap ratio diverges or the Higgs oscillation is not resolved, the one-loop truncation would be falsified. Alternatively, a numerically exact quantum Monte Carlo calculation of the spectral function near the critical point could check whether corrections beyond one loop alter the suppression.","supporting_citations":[],"review_version":1}