{"id":"3264715f-0239-46c2-937d-0cc0b573231f","arxiv_id":"2506.16973","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Gap-protected countertwisting with locally interacting spins reaches Heisenberg-limited metrological gain in time t_H ∝ N^{γ/d} log N, sublinear in N for 2D dipolar and 3D van der Waals systems.","lead":"This paper shows that spin squeezing can reach the ultimate Heisenberg sensitivity limit using only local, power-law interactions between spins, rather than the global all-to-all couplings usually required. The concrete result is a predicted time scaling, t_H ∝ N^{1/2} log N for 2D dipolar arrays, which makes near-Heisenberg-limited sensing plausibly reachable in molecules, Rydberg atoms, and solid-state spins.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted t_H scaling assumes the quadratic spin-wave theory remains valid up to χt ≈ log N, where the k=0 occupation is O(N); the omitted four-magnon coupling of the k=0 pump to k=kc modes is O(J), exceeding the linear gap, and the numerics do not yet confirm convergence (γ≈1.3 vs 1 at N=4096).","rationale":"After reading the full manuscript and supplement, the strongest claim is the sublinear time to near-Heisenberg-limited sensitivity with local interactions. For this to hold, the k≠0 modes must remain unexcited for the entire squeezing time. The paper's derivation (Eqs. 2-3) establishes this only within the quadratic, low-occupation approximation. The reader's weakest assumption—validity of the linear spin-wave expansion up to optimal time—is exactly where the argument is thinnest. My stress-test sharpens it: at χt ≈ log N, the k=0 occupation is O(N), so the leading nonlinear vertices are not suppressed by any small parameter; the pump-to-mode-kc coupling is O(J) while the protective gap is O(J N^{−γ/d}). Without an explicit estimate of the resulting scattering rate (or a demonstration of kinematic suppression), the analytic prediction of t_H is incomplete. The DTWA simulations are a useful check but semiclassical, and the published fit γ≈1.3 at N=4096 (Fig. S3b) shows the prediction γ=1 is not yet confirmed at the largest accessible sizes. These considerations support the reader's CONDITIONAL verdict; they do not overturn the likely correctness of the proposal, which is corroborated by exact diagonalization at N=16, spectral checks, and the all-to-all baseline. Hence no change to the verdict.","tokens_in":23379,"tokens_out":29235,"duration_ms":317087,"concrete_test":"Perform a Floquet stability analysis of the time-dependent uniform spin-wave condensate: linearize the full classical (or exact HP, including quartic terms) equations of motion around the uniform trajectory with amplitude α_0(t)=sinh(χt), and compute the Floquet/Lyapunov exponent Γ_kc of the smallest nonzero momentum mode. If Γ_kc ≥ χ before χt ≈ log N, the linear stability condition (Eq. 3) is insufficient and the predicted t_H scaling fails. A complementary numerical check: initialize DTWA at N=4096 with a small coherent amplitude at momentum kc, evolve to χt ≈ log N at χ=χ_c, and measure the amplification of that mode; if it grows by more than a factor of order one, the k=0 pump nonlinearly excites k≠0 modes, invalidating the two-mode model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central scaling law t_H ∼ N^{γ/d} log N rests on the quadratic Holstein-Primakoff Hamiltonian (Eq. 2) and the stability bound χ < χ_c = 1 − 1/f_kc (Eq. 3). These are early-time conditions: |ω_k/χ_k| > 1 is evaluated for infinitesimal spin waves. Reaching a Heisenberg sensitivity N(Δφ)² = c/N requires evolving to χt ≈ log N, at which the k=0 magnon occupation is n_0 ∼ N (the TACT optimum). The local boson occupation n_0/N is then O(1), so the quartic terms omitted from Eq. 2 are not small. The four-magnon process 2(k=0) → k+(−k) has effective coupling (J/N)×a_0² ∼ J, while the mode frequency that should suppress it is ω_kc ∼ J(1−1/f_kc) ∼ J N^{−γ/d}; for large N the coupling exceeds the protective gap by N^{γ/d}. Neither the paper nor the supplement estimates this nonlinear scattering rate or shows it is kinematically suppressed. DTWA, which partly includes these nonlinearities, is only semiclassical; its largest-size fit gives γ≈1.3 for (2,3) instead of the predicted γ=1 (Fig. S3b), so the asymptotic scaling is not empirically converged. Thus the load-bearing assumption is that nonlinearities do not drive a modulational instability of the k=0 condensate before t_H; this is asserted, not derived.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a gap-protected countertwisting Hamiltonian, Eq. (1), with power-law Heisenberg and countertwisting interactions, and analyzes its spin-squeezing dynamics. Using a Holstein-Primakoff reduction to the quadratic spin-wave Hamiltonian Eq. (2), the authors derive the critical squeezing rate χ_c = 1 − 1/f_{k_c} (Eq. 3) and predict a time t_H ∼ N^{γ/d} log N to approach Heisenberg-limited sensitivity N(Δφ)² = c/N, with constants c_sq ≈ 3.9 and c_QFI ≈ 1.56 inherited from the all-to-all countertwisting model. The claims are supported by DTWA simulations up to N = 4096, benchmarked against exact diagonalization at small N, together with an analysis of Floquet implementations and robustness to disorder and density fluctuations.","tokens_in":23574,"tokens_out":4625,"duration_ms":50709,"significance":"If the central scaling claim holds, the paper would constitute an important step for spin squeezing with local and power-law interactions: it identifies a concrete Hamiltonian whose metrological gain reaches the Heisenberg limit in sublinear time for two-dimensional dipolar and three-dimensional van der Waals systems, and it provides a transparent, non-circular derivation of the allowable squeezing rate from the momentum-space couplings f_k. The paper is also commendable for computing the all-to-all constants c_sq and c_QFI from the collective model and then comparing them with local-interaction simulations, rather than extracting them as fit parameters. The main weakness is that the asymptotic time-to-Heisenberg scaling relies on the validity of quadratic spin-wave theory up to occupation numbers of order N, and the numerical evidence at the largest simulated sizes has not yet converged to the predicted exponent. These issues are load-bearing but, in my view, addressable within the manuscript's scope.","major_comments":[{"comment":"The derivation of χ_c and of the time scaling t_H ∼ N^{γ/d} log N is based on the stability condition |ω_k/χ_k| > 1, which is evaluated for infinitesimally populated spin-wave modes. However, reaching the claimed sensitivity requires evolving to χt_opt ∼ log N, where the k = 0 mode occupation is n_0 ∼ N. At that point the quartic terms omitted from Eq. (2) are not small: the four-magnon process 2(k=0) → k+(−k) has an effective coupling of order J, while the protective gap at the nearest competing mode is ω_{k_c} ∼ J(1−1/f_{k_c}) ∼ J N^{−γ/d}. For large N this nonlinear coupling exceeds the linear gap, and the manuscript neither estimates the associated scattering rate nor shows that the process is kinematically suppressed. Because the central claim is precisely that the dynamics remains confined to the k = 0 mode up to t_H, this missing estimate is load-bearing.","section":"Eqs. (2)-(3); Supplemental Sec. II.1; Fig. S3(b)"},{"comment":"The numerical support for the asymptotic scaling is not yet convergent. At (d, α) = (2, 3), the fit of the QFI scaling N(Δφ_QFI)² ∝ (χ_c/χ)^{−2/γ} yields γ ≈ 1.3 at the largest simulated system size N = 4096, whereas the predicted value is γ = 1. The supplement states that the fitted values 'approach' the prediction, but the trend is still approximately 30% away at the largest size, and DTWA is a semiclassical method that cannot by itself exclude a quantum contribution to the exponent. This is precisely the kind of deviation that the omitted nonlinear terms could produce, so the empirical case for the claimed asymptotic time scaling is incomplete.","section":"Supplemental Sec. III.1, Fig. S3(b); main text Fig. 2(b-c)"},{"comment":"The paper uses F = 4(ΔS_max)² for the quantum Fisher information and extracts the constant c_QFI ≈ 1.56 from DTWA simulations, but DTWA operates with second moments of a semiclassical distribution. For the strongly non-Gaussian states reached near the Heisenberg limit, this expression is not guaranteed to equal the true QFI, and the supplement's echo protocol is the concrete proposal for accessing this sensitivity. However, the echo analysis is performed only within the same DTWA framework and is not benchmarked against exact diagonalization for the local-interaction case. A short discussion of the validity of F = 4(ΔS_max)² for the states generated by Eq. (1) would strengthen the claim that the reported QFI values are rigorous metrological bounds.","section":"Main text 'In approaching the Heisenberg limit'; Supplemental Sec. IV"}],"minor_comments":[{"comment":"The axis annotation '40¯6' appears to be a typo for 4096; please correct it.","section":"Fig. 2(a)"},{"comment":"The open-boundary generalization of the critical rate uses the two largest eigenvalues f̃ and f̃′ of the real-space coupling matrix, but the main text does not define these quantities; adding a sentence tying f̃ and f̃′ to f_{k_c} would help the reader.","section":"Eq. (3) and Supplemental Sec. II.3"},{"comment":"The abstract uses 'local interactions' while the model includes power-law interactions ranging down to all-to-all α = 0; please clarify in the abstract what is meant by 'local' (for example, 'short-range or power-law interacting').","section":"Abstract and introduction"},{"comment":"The horizontal axis label in the main-text figure appears as 'πtott' but should read J_tot t; also, the analytical model of Eq. (S13) is first used in the main text, so citing the supplement section at that point would be helpful.","section":"Fig. 4(a) and Supplemental Sec. V"},{"comment":"Several quantitative statements, including the scalings in Fig. 5(a) and the three-regime model in Fig. 4, are deferred to the supplement via a single reference [49]; please add section-specific pointers to the supplement so that the reader can locate the derivations.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper's central idea is attractive and the derivation of χ_c is not circular: it follows from the Hamiltonian couplings through Eq. (3), and the all-to-all constants are computed independently and then compared with local simulations. The main risk is the unestimated nonlinear instability of the k=0 condensate at the optimal time; if the authors can supply a quantitative estimate of four-magnon scattering, or additional numerical evidence at larger N and a direct check that the k=0 occupation remains dominant, the manuscript would be publishable. I recommend major revision rather than rejection because the issue is a missing derivation/support for a key assumption, not an internal contradiction in the presented results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Best to give you my read. The paper proposes a gap-protected countertwisting Hamiltonian with power-law interactions and claims that, for a sufficiently slow squeezing rate, one recovers the all-to-all Heisenberg scaling in time t_H ∝ N^{γ/d} log N, sublinear for 2D dipolar and 3D van der Waals systems. The result is genuinely new: prior gap-protected squeezing work was one-axis twisting (XXZ), and local countertwisting had only been studied at mean-field level. So the question is good and the answer is non-obvious.\n\nWhat I like: the spin-wave derivation of the critical rate χ_c = 1 − 1/f_kc is clean and parameter-free; the constants c_sq and c_QFI are inherited from the all-to-all model, not fitted. The DTWA simulations up to N=4096 show a reasonable collapse when rescaled by χ_c. The comparison with XXZ in Fig. 4 is informative, and the Floquet-engineering section adds practical value. The robustness-to-fluctuations argument (fixed squeezed quadrature versus rotating) is simple and convincing.\n\nThe soft spot is the one the stress-test flags. The quadratic Holstein-Primakoff expansion is validated at early times and against spectral gaps, but the protocol reaches the Heisenberg limit at χt ≈ log N, when the k=0 magnon occupation is O(N). At that point the local boson density is O(1), so the quartic terms are not small. The paper does not estimate the rate of four-magnon processes that couple the k=0 pump to k=kc modes, nor does it show they are kinematically suppressed. This is a missing argument rather than a demonstrated failure, but it is load-bearing for the claimed time scaling. The DTWA fits do not settle it: the fitted exponent in Fig. S3b is γ≈1.3 at N=4096, trending toward the predicted 1 but not converged. Also, no code or data is shipped, and the experimental section assumes coherent evolution without a decoherence error budget.\n\nNet: this deserves a serious referee. The derivation is honest and the numerics are consistent with the prediction, but the asymptotic scaling rests on an unproven stability assumption. A referee should ask for either a quantitative estimate of the quartic couplings or a more direct large-N calculation. I would send it to peer review and would cite it once the nonlinearity question is addressed.","headline":"A clean, parameter-free derivation of a sublinear time to Heisenberg scaling for local countertwisting, with a legitimate but unresolved question about nonlinear corrections at the optimal time.","tokens_in":24372,"tokens_out":6442,"would_cite":true,"duration_ms":66785,"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":"Local spins reach near-Heisenberg sensitivity sublinearly in time","keywords":["spin squeezing","Heisenberg limit","two-axis countertwisting","gap protection","power-law interactions","quantum Fisher information","spin-wave theory","quantum metrology"],"falsifier":"At the largest simulated size $N=4096$, the fitted exponent extracted from the quantum Fisher information is still about $\\gamma\\approx 1.3$ rather than the predicted 1; if simulations at $N=10^4$ to $10^5$ show the fitted exponent plateauing above 1, the claimed sublinear time scaling fails.","tokens_in":2336,"feed_emoji":"⚛️","tokens_out":3286,"duration_ms":108917,"temperature":0.7,"pith_summary":"This paper argues that a local version of the two-axis countertwisting Hamiltonian, protected by an energy gap from a Heisenberg coupling, can squeeze a collection of spins close to the Heisenberg limit of phase sensitivity even though interactions decay with distance. The central result is a rate condition: if the countertwisting rate $\\chi$ is kept slower than a critical value $\\chi_c$ that shrinks as $N^{-\\gamma/d}$, the dynamics stays in the collective zero-momentum mode and the metrological gain recovers the all-to-all scaling $G\\propto N$. For 2D dipolar ($\\alpha=3$) and 3D van der Waals ($\\alpha=6$) interactions, the time to reach that scaling grows only sublinearly with particle number, $t_H\\sim N^{1/2}\\log N$ and $t_H\\sim N^{2/3}\\log N$. If true, this removes the need for all-to-all interactions for near-Heisenberg quantum-enhanced sensing and gives specific near-term platforms, including molecules, Rydberg atoms, and solid-state spins, a concrete protocol.","feed_headline":"Local spins reach near-Heisenberg sensitivity sublinearly in time","feed_subtitle":"A gap-protected countertwisting protocol recovers all-to-all sensitivity for 2D dipolar and 3D van der Waals systems.","key_machinery":"The central object is a linear spin-wave (Holstein-Primakoff) expansion of the XYZ Hamiltonian around the fully polarized state, which reduces the early-time dynamics to coupled bosonic modes: $H\\approx\\sum_k[\\chi_k(a^\\dagger_{-k}a^\\dagger_k+a_k a_{-k})/2+\\omega_k a^\\dagger_k a_k]$. For $H_{\\mathrm{gct}}$, the zero-momentum mode is resonant ($\\omega_0=0$) and squeezes exponentially, while each nonzero mode $k$ is detuned by $\\omega_k=f_k-1$ and parametrically driven at rate $\\chi_k=\\chi f_k$; requiring $|\\omega_k/\\chi_k|>1$ yields the critical rate $\\chi_c=1-1/f_{k_c}$. This stability condition, together with the dispersion $\\omega_k\\propto k^\\gamma$ ($\\gamma=\\min(2,\\alpha-d)$), turns the correlation length $L_c\\propto\\chi^{-1/\\gamma}$ into a prediction for metrological gain $G\\propto L_c^d\\propto\\chi^{-d/\\gamma}$, and converts the rate bound into the time scaling $t_H\\sim N^{\\gamma/d}\\log N$. The mechanism is gap protection: the Heisenberg coupling $s_i\\cdot s_j$ energetically suppresses leakage out of the collective $S=N/2$ manifold, letting local interactions mimic permutation-symmetric entanglement.","core_discovery":"The paper's central claim is that the gap-protected countertwisting Hamiltonian $H_{\\mathrm{gct}}=\\sum_{i\\neq j}J(r_{ij})[s_i\\cdot s_j+(\\chi/2)(s_i^+s_j^++s_i^-s_j^-)]$ with power-law couplings $J\\propto r^{-\\alpha}$ reproduces, for sufficiently slow countertwisting, the Heisenberg scaling $G\\propto N$ of the all-to-all countertwisting model. The Heisenberg term is the load-bearing protector: it opens an energy gap between manifolds of different total spin, constraining the dynamics to the permutation-symmetric $k=0$ mode while countertwisting parametrically amplifies that mode. Spin-wave theory sets the critical countertwisting rate at $\\chi_c=1-1/f_{k_c}\\propto N^{-\\gamma/d}$ with $\\gamma=\\min(2,\\alpha-d)$; below this rate the optimal squeezing parameter approaches $N(\\Delta\\varphi_{\\mathrm{sq}})^2=c_{\\mathrm{sq}}/N$ with $c_{\\mathrm{sq}}\\approx 3.9$, and the quantum Fisher information gives $N(\\Delta\\varphi_{\\mathrm{QFI}})^2=c_{\\mathrm{QFI}}/N$ with $c_{\\mathrm{QFI}}\\approx 1.56$. The time to reach Heisenberg scaling is therefore $t_H\\sim N^{\\gamma/d}\\log N$, which is sublinear for the experimentally relevant cases $(d,\\alpha)=(2,3)$ and $(3,6)$. Numerical simulations based on the discrete truncated Wigner approximation up to $N=4096$, benchmarked against exact diagonalization at $N=16$, support the scaling predictions and the predicted robustness to disorder and density fluctuations.","pith_inferences":["Beyond the paper: because the distance-dependent plateau in pairwise correlations directly encodes the dispersion exponent $\\gamma$, a real-space measurement of $C_{\\max}(r)$ in a dipolar platform could independently confirm the critical-rate prediction $\\chi_c\\propto N^{-\\gamma/d}$ without requiring very large atom numbers.","Beyond the paper: the fixed squeezed quadrature of countertwisting suggests the protocol could be adapted to parameter estimation in the presence of slowly varying fields, where one-axis twisting protocols degrade because their squeezed axis rotates with the unknown parameter.","A testable extension the authors do not pursue is a time-dependent rate $\\chi(t)$, starting fast for early exponential squeezing and slowing near the critical rate to grow the correlation plateau; whether this outcompetes the optimal fixed-rate schedule is open.","The paper leaves open whether $t_H$ can be improved to the Lieb-Robinson speed limits; comparing the achieved exponents $\\gamma/d$ (1/2 for 2D dipolar, 2/3 for 3D van der Waals, 2 for 1D short-range) against those bounds is a concrete quantitative check."],"forward_implications":["For 2D dipolar ($\\alpha=3$) and 3D van der Waals ($\\alpha=6$) interactions, squeezing at the critical rate reaches near-Heisenberg sensitivity within $t_H\\sim N^{1/2}\\log N$ and $t_H\\sim N^{2/3}\\log N$, respectively.","If practical constraints cap the total interaction time, a smaller gain $G$ can still be obtained by squeezing faster: $\\chi\\sim G^{-\\gamma/d}$ for time $t\\sim G^{\\gamma/d}\\log G$.","The protocol inherits countertwisting's robustness to fluctuations in interaction strength: under fractional density fluctuations $\\epsilon$, the optimal squeezing degrades as $N^{-(1-\\epsilon)}$ rather than saturating like one-axis twisting, which hits a floor near $\\epsilon^2$.","Because countertwisting and Heisenberg spreading are distinct terms, the Hamiltonian can be generated by Floquet engineering from native Ising or XY interactions with global $\\pi/2$ pulses, with accessible rates $0\\leq\\chi\\leq\\chi_{\\max}$ depending on the native anisotropy.","An echo protocol that reverses the sign of interactions after a small rotation accesses the quantum-Fisher-information-bound sensitivity, extending the metrological gain beyond what the squeezing parameter alone shows."],"supporting_citations":[{"why":"It supplies the two-axis countertwisting Hamiltonian and the all-to-all Heisenberg-scaling baseline that this paper generalizes.","marker":"[4]"},{"why":"It introduces the use of Heisenberg interactions to gap the different total-spin manifolds, which is the gap-protection mechanism used throughout.","marker":"[20]"},{"why":"It provides the dispersion law $\\gamma(\\alpha,d)=\\min(2,\\alpha-d)$ for power-law interactions used to derive $\\chi_c\\propto N^{-\\gamma/d}$.","marker":"[21]"},{"why":"It shows gap-protected spin squeezing under XXZ models with power-law interactions, the theoretical basis this paper extends to countertwisting.","marker":"[24]"},{"why":"It supplies the phase-transition comparison for gap-protected one-axis twisting used to benchmark $H_{\\mathrm{gct}}$ against $H_{\\mathrm{XXZ}}$.","marker":"[28]"},{"why":"It demonstrates dipolar spin squeezing in Rydberg atom arrays, one experimental context for implementing the proposed protocol.","marker":"[31]"},{"why":"It demonstrates spin squeezing in nitrogen-vacancy-center dipolar systems and provides a platform for gap-protected countertwisting.","marker":"[35]"},{"why":"It demonstrates Floquet engineering of countertwisting with dipolar molecules, making the proposed Hamiltonian experimentally available.","marker":"[38]"},{"why":"It supplies the discrete truncated Wigner approximation used for the large-$N$ numerical simulations up to $N=4096$.","marker":"[48]"},{"why":"It contains the benchmark, spin-wave derivations, and scaling fits on which the predictions rely.","marker":"[49]"}],"fun_headline_variants":["Local spins reach Heisenberg limit sublinearly","Countertwisting local spins: Heisenberg-limited sensitivity","Sublinear-time Heisenberg-limited spin squeezing","Gap-protected protocol for Heisenberg-limited metrology","Robust local interactions approach Heisenberg bound"],"cache_read_input_tokens":26112,"weakest_assumption_plain":"The main calculation assumes that the linear spin-wave description of the early dynamics stays accurate up to the optimal squeezing time, so that nonlinear corrections and leakage out of the protected total-spin manifold can be neglected at large particle number.","fun_headline_variants_meta":{"raw":{"variants":["Local spins reach Heisenberg limit sublinearly","Countertwisting local spins: Heisenberg-limited sensitivity","Sublinear-time Heisenberg-limited spin squeezing","Gap-protected protocol for Heisenberg-limited metrology","Robust local interactions approach Heisenberg bound"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000432,"raw_usage":{"total_tokens":2242,"prompt_tokens":1024,"completion_tokens":1218,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":640,"completion_tokens_details":{"reasoning_tokens":1146}},"tokens_in":640,"tokens_out":1218,"duration_ms":11139,"temperature":1.0,"reasoning_tokens":1146,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:17:18.284535+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At the largest simulated size $N=4096$, the fitted exponent extracted from the quantum Fisher information is still about $\\gamma\\approx 1.3$ rather than the predicted 1; if simulations at $N=10^4$ to $10^5$ show the fitted exponent plateauing above 1, the claimed sublinear time scaling fails.","supporting_citations":[{"cited_title":"Block, B","cited_arxiv_id":null,"evidence_quote":"It supplies the phase-transition comparison for gap-protected one-axis twisting used to benchmark $H_{\\mathrm{gct}}$ against $H_{\\mathrm{XXZ}}$."},{"cited_title":"Bornet, G","cited_arxiv_id":null,"evidence_quote":"It demonstrates dipolar spin squeezing in Rydberg atom arrays, one experimental context for implementing the proposed protocol."},{"cited_title":"[65, 66], for additional details and supporting calculations","cited_arxiv_id":null,"evidence_quote":"It contains the benchmark, spin-wave derivations, and scaling fits on which the predictions rely."}],"review_version":1}