{"id":"799e7df1-ff6c-4d5b-8477-2a78cfa3347a","arxiv_id":"2412.10542","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In the 3D Kitaev model on the hyperhoneycomb lattice, local spin correlations develop a persistent oscillation below the thermal transition, and the fitted amplitude of this oscillation tracks the Wilson loop order.","lead":"This paper calculates how local spin correlations evolve in time in a three-dimensional Kitaev spin liquid model, using quantum Monte Carlo. It finds a clear oscillation that appears only below the thermal transition, which it calls a dynamic order and proposes as a detection signature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Dynamic-order claim rests on a four-parameter fit to N=4 QMC data with no error bars or finite-size scaling; the T=0 parton comparison does not establish an order parameter at finite T.","rationale":"The reader's weakest assumption identifies exactly the same soft spot: the QMC evidence is a single finite-size system with a fixed Wilson-loop criterion for Tw and no finite-size scaling or statistical error propagation. My analysis agrees that this is the most load-bearing limitation, because the central claim is not merely that an oscillation exists in the spin-liquid phase (which the T=0 parton calculation supports) but that Adyn(T) acts as a dynamic order parameter for the thermal transition. The four-parameter fit is the bridge between raw QMC data and the order-parameter claim, and without error bars or larger-N data that bridge is not secure.\n\nI credit the paper for its exact real-time formula and for the sign-ambiguity resolution, which are nontrivial and likely correct. The observation of an oscillation below Tw is also a defensible numerical finding. However, the step from that observation to 'order parameter of the thermal transition' requires the thermodynamic-limit test described above. Since the reader already reached CONDITIONAL for these same reasons, I do not think a new verdict is needed; the recommendation is UNCHANGED, with the concrete finite-size and bootstrap analysis as the condition for strengthening to ACCEPT.","tokens_in":15344,"tokens_out":4730,"duration_ms":48962,"concrete_test":"Run the identical QMC for N=5 and N=6 with periodic boundary conditions, using at least 10^5 measurement sweeps after 10^4 thermalization sweeps. For each N, fix Tw(N) with the same |<W>|=0.03 criterion, then fit Re<S^z_j(t)S^z_j(0)> to Eq. (6) over a fixed time window (e.g., 0 < 3Jt < 40) using bootstrap resampling over QMC sweeps. Compute Adyn(T/Tw) at T/Tw = 0.9, 1.0, 1.05, and 1.1 with confidence intervals. If Adyn at T/Tw=1.05 or 1.1 is significantly nonzero at N=5 or N=6, or if the ratio Adyn(T<Tw)/Adyn(T>Tw) does not grow with N, the sharp-onset/order-parameter claim is not supported. In addition, check whether Tw(N) extrapolates so that the Adyn onset coincides with the thermodynamic transition rather than with the finite-size Wilson-loop threshold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the identification of the fitted amplitude Adyn(T) in Eq. (6) as an order parameter for the thermal transition. That claim requires Adyn(T)>0 below Tw and Adyn(T)=0 above Tw in the thermodynamic limit, with a sharp onset at the true transition. The evidence is QMC at a single main size, N=4, with about 6000 measurement steps after 3000 thermalization steps, no statistical error bars, and no finite-size scaling. Tw itself is set by the arbitrary criterion |<W>|=0.03; the authors note that this Wilson-loop onset can differ from the specific-heat peak for finite N and that the two coincide only in the thermodynamic limit.\n\nA four-parameter damped-cosine fit to a finite-time, noisy, monotonically decaying correlation function can return a nonzero Adyn even without a true long-lived oscillation, and with no residuals, fit-range sensitivity, or bootstrap errors shown, the statement that Adyn 'precisely follows' <W> is not quantitatively supported. The parton comparison in SM Sec. I is performed at T=0, using N=5-8 (plus a single N=200 example), and it only demonstrates that an oscillation exists in the zero-temperature spin-liquid ground state; it does not test the finite-temperature onset or the vanishing of Adyn above Tw. The paper also states that the oscillation period varies with system size, which further suggests that the fitted frequency and amplitude may track finite-size Majorana levels rather than a thermodynamic order parameter.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Kitaev model on the hyper-honeycomb lattice at nonzero temperature and computes the real-time local spin correlation function using a Majorana representation combined with quantum Monte Carlo. The authors report that the correlation function develops a damped oscillatory component below a temperature Tw, which they identify with the thermal transition of the Z2 quantum spin liquid using the Wilson loop expectation value. They extract a fitted amplitude Adyn(T) from a four-parameter damped-cosine ansatz and propose that Adyn(T) acts as a \"dynamic order parameter\" for the featureless thermal transition. A spin structure factor calculation and a zero-temperature parton calculation are presented in support of the interpretation.","tokens_in":15716,"tokens_out":4252,"duration_ms":41969,"significance":"If established, the proposed dynamic order would be a practically valuable observable for detecting featureless thermal transitions in quantum spin liquids, since it is local, dynamical, and accessible in principle to neutron scattering and related probes. The paper also contains a genuine technical contribution: the exact finite-temperature Majorana determinant expression for the real-time correlation function and the Pfaffian-based resolution of the square-root sign ambiguity are nontrivial and enable QMC evaluation without analytic continuation. It is also fair that the Adyn-versus-Wilson-loop agreement is an empirical numerical comparison, not a circular definition. However, the central claim that Adyn(T) is an order parameter for the thermodynamic transition is currently supported only by a four-parameter fit to N=4 QMC data with no statistical uncertainties, no finite-size scaling, and no finite-temperature parton comparison, so the significance is not yet established at the level claimed.","major_comments":[{"comment":"The central claim that Adyn(T) behaves as an order parameter rests on fits of the four-parameter damped-cosine ansatz in Eq. (6) to QMC time traces at N=4, but the paper reports no statistical uncertainty in Adyn(T), no fit-range sensitivity, no residuals, and no comparison with a null model such as a purely exponential decay without an oscillatory component. A four-parameter cosine can absorb noise in a short, noisy, monotonically decaying time series and return a nonzero Adyn even when no long-lived oscillation exists. Please provide bootstrap or jackknife error bars on Adyn(T), a test of the sensitivity of Adyn to the fitting window, and a goodness-of-fit comparison against a non-oscillatory model before asserting that Adyn 'precisely follows' the Wilson loop tendency.","section":"Dynamic order, Eq. (6)"},{"comment":"The identification of Tw via the criterion |<W>| = 0.03 for N=4, together with the authors' statement that this onset can differ from the specific-heat peak for finite size and coincides with it only in the thermodynamic limit, means that the apparent sharp onset of Adyn around Tw has not been shown to locate the thermodynamic transition. The authors state that calculations were performed for N=3, 4, and 5, but the main results and the Adyn extraction are shown only for N=4. A finite-size study of Adyn(T) for these sizes, with the same fitting procedure, and an analysis of how Tw converges to the specific-heat peak as N increases are needed before Adyn can be called an order parameter. The authors also note that the oscillation period varies with system size, which further suggests that the fitted frequency and amplitude may track finite-size Majorana levels rather than a thermodynamic order parameter.","section":"Dynamic order, Tw identification"},{"comment":"The comparison with the parton analysis is performed at T=0, whereas the dynamic-order claim concerns a finite-temperature transition. The T=0 parton calculation, with N=5-8 and N=200, demonstrates that an oscillation exists in the zero-temperature spin-liquid ground state, but it does not test whether the oscillation appears only below Tw or whether Adyn vanishes above Tw. To support the finite-temperature dynamic-order scenario, the manuscript needs either a finite-temperature parton calculation near Tw or a clear finite-size scaling analysis of the QMC Adyn(T); the current 'adiabatic connection' argument is not sufficient because the transition itself is a finite-temperature phenomenon.","section":"Supplementary Material Sec. I, Fig. 4"}],"minor_comments":[{"comment":"The figure caption refers to the dynamic order as 'Ap', while Eq. (6) and the text define 'Adyn'; please unify the notation.","section":"Fig. 1 caption"},{"comment":"The Wilson loop operator W is cited only as 'defined in [6]'; for self-containedness, please define it explicitly in the main text and state how its expectation value distinguishes the confined and deconfined phases.","section":"Model and Dynamic order"},{"comment":"The time-averaging cutoff T0 in Eq. (4) is easily confused with the temperature T; consider using a different symbol such as tau or t_cutoff.","section":"Eq. (4)"},{"comment":"The reported peak positions depend on the chosen broadening eta (0.09 in Fig. 2 and 0.01 in Fig. 3); please state how the peak frequencies shift with eta and why the chosen values do not affect the qualitative conclusion.","section":"Fig. 2 and Fig. 3"},{"comment":"In the text after Eq. (7), 'all the positive eigenvalues of the matrix Aj,alpha' appears to contain a typo; the matrix should be A or A' depending on the context.","section":"Supplementary Material Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The technical core, especially the exact Majorana determinant formula and the Pfaffian sign resolution, is solid and likely publishable. My concern is the gap between the evidence and the headline claim of a 'dynamic order parameter': the load-bearing evidence is a single-size fit with no error bars and no finite-size scaling. This is fixable within the manuscript's scope by adding statistical analysis, finite-size data, and a finite-temperature comparison, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nThe genuinely new thing here is a finite-temperature real-time spin correlation calculation for the 3D hyperhoneycomb Kitaev model, using an exact Majorana determinant formula that avoids numerical analytic continuation. The authors also resolve the Pfaffian sign ambiguity, a technical step that matters. The reported onset of oscillation in the local spin correlation below Tw is directly visible in Fig. 2, not something manufactured by the fit. That is a solid observation.\n\nThe paper's weakness is the leap from that observation to the claim that Adyn(T) acts as an order parameter for the thermal transition. The evidence for the Adyn(T) vs ⟨W⟩ correspondence is a four-parameter damped-cosine fit applied to one main system size, N=4, with roughly 6000 Monte Carlo steps and no error bars, no residuals, and no fit-range sensitivity analysis. A four-parameter fit to a finite-time, noisy, monotonically decaying correlation can return a nonzero damped-cosine amplitude even without true long-lived oscillation. So the \"precisely follows\" statement is not yet quantitatively supported. Tw itself is set by |⟨W⟩|=0.03, which is an arbitrary threshold; the authors themselves note the finite-N mismatch with the specific-heat peak. The parton analysis in the SM is done at T=0, so it shows the oscillation exists in the ground state, but it does not test whether Adyn vanishes above Tw. The statement that the oscillation period varies with system size further suggests the fitted frequency may track finite-size Majorana levels rather than a thermodynamic order parameter.\n\nThat said, these are addressable. The qualitative physics—a sharp growth of coherent spin dynamics entering the spin-liquid phase—is almost certainly real. The dynamic order concept is interesting and experimentally relevant. What is missing is statistical and finite-size support, not a fundamental error in approach. A serious referee could ask for error bars, bootstrap or residual analysis, N=3,4,5 comparison with finite-size scaling, and a clearer definition of Tw. I hope the authors can supply that.\n\nThis paper deserves a serious referee. I would not desk-reject it. If it is my call, I'd send it to review with a request for strengthening the statistical analysis, and would expect heavy revision before acceptance.","headline":"Dynamic-order claim needs stronger finite-size evidence, but the exact real-time correlation method and the visible oscillation onset make this a worthwhile paper for a serious referee.","tokens_in":16216,"tokens_out":2999,"would_cite":true,"duration_ms":25081,"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 a coherent oscillation in the local spin correlation function of the three-dimensional Kitaev model marks the quantum spin liquid phase at non-zero temperatures, and that its fitted amplitude acts as an order…","keywords":["quantum spin liquid","dynamic order","Kitaev model","hyper-honeycomb lattice","quantum Monte Carlo","spin correlation function","Wilson loop","thermal phase transition"],"falsifier":"A finite-size scaling study of the fitted amplitude Adyn(T) across N=4,5,6,7,8 with the same Wilson-loop criterion would settle whether the sharp onset of oscillation survives in the thermodynamic limit; if the fitted Adyn(T) does not sharpen toward a step function as N grows, or remains nonzero above Tw in the thermodynamic limit, the dynamic-order claim would be falsified.","tokens_in":15126,"feed_emoji":"🌀","tokens_out":8635,"duration_ms":74369,"temperature":0.7,"pith_summary":"This paper aims to establish that a quantum spin liquid, a state lacking any symmetry-breaking order parameter, can be identified at non-zero temperatures by a sharp onset of coherent oscillation in its local spin correlation function. Using quantum Monte Carlo on the Kitaev model on the hyper-honeycomb lattice, the authors find that below a thermal transition temperature, the fitted amplitude Adyn(T) of the dominant oscillatory mode closely tracks the Wilson loop expectation value that theoretically defines the spin liquid phase. They propose that Adyn(T) acts as an order parameter for the otherwise featureless transition between the spin liquid and a trivial paramagnet, making the quantum spin liquid detectable through experimentally accessible spin dynamics such as neutron scattering or NMR.","feed_headline":"A local spin oscillation tags the quantum spin liquid phase","feed_subtitle":"Its fitted amplitude tracks the topological Wilson loop, turning a featureless transition into a measurable dynamic order.","key_machinery":"The key object is the exact Majorana-fermion expression for the dynamical local spin correlation function, obtained by extending the exact two-dimensional Kitaev-model result of reference [29] to the hyper-honeycomb lattice. This expression rewrites ⟨Szj(t)Szj(0)⟩ as a sum over Z2 gauge-field configurations of ratios of determinants of the free-Majorana matrix, allowing quantum Monte Carlo evaluation of real-time dynamics without analytic continuation. A four-parameter fit of the correlation function to Astat + Adyn e−t/tφ cos(φt) then extracts the dynamic-order amplitude Adyn(T); the paper also resolves the sign ambiguity of the square-root determinants using Pfaffian properties, which were previously a technical obstacle.","core_discovery":"Below the thermal transition temperature Tw, the local spin correlation function ⟨Szj(t)Szj(0)⟩ develops a persistent oscillation with a temperature-independent frequency, while above Tw it decays monotonically. Fitting these correlations to a damped cosine form, the extracted amplitude Adyn(T) vanishes at Tw and precisely follows the tendency of the non-local Wilson loop ⟨Ŵ⟩, which detects the condensation of string-type excitations. The authors therefore argue that Adyn(T) plays the role of an order parameter for the thermal transition and name it the dynamic order, providing a local and experimentally accessible probe of the quantum spin liquid phase.","pith_inferences":["If the correspondence between Adyn(T) and ⟨Ŵ⟩ holds beyond the fitted system sizes, local spin dynamics may encode the thermal behavior of non-local flux excitations, suggesting that similar local dynamic markers could exist in other topologically ordered states.","In two-dimensional Kitaev spin liquids, where the thermal transition becomes a crossover, the same analysis may reveal a smeared but still detectable oscillation onset, giving a finite-temperature diagnostic in a regime where no true phase transition exists.","The single-moded cosine fit may under-report the dynamic order at lower temperatures; a multi-mode fit consistent with the spin structure factor's flat dispersions could reveal additional structure tied to fractionalized excitations."],"forward_implications":["A local, experimentally accessible spin response can identify the quantum spin liquid phase at finite temperature, replacing non-local probes like the Wilson loop.","The spin structure factor below Tw exhibits temperature-independent flat dispersions, giving a spectroscopic fingerprint that distinguishes the spin liquid from a trivial paramagnet.","The dynamic order can detect featureless thermal transitions without symmetry breaking, as reported in candidate materials such as β-Li2IrO3 and Sr2VO3FeAs.","The exact QMC method for dynamical spin correlations without analytic continuation can be extended to other exactly solvable spin models."],"supporting_citations":[{"why":"Introduces the exactly solvable Kitaev model and its Majorana-fermion representation, which the present calculation builds on.","marker":"[4]"},{"why":"Establishes the thermal transition in the hyper-honeycomb Kitaev model and defines the Wilson loop operator used to locate Tw.","marker":"[6]"},{"why":"Provides the zero-temperature parton analysis of the hyper-honeycomb spin liquid that the finite-temperature QMC results are compared against for the thermodynamic limit.","marker":"[11]"},{"why":"Supplies the method for calculating dynamical spin correlations in the Kitaev model via Majorana fermions, extended here to finite temperature.","marker":"[12]"},{"why":"Gives the exact expression for the dynamical local spin correlation function in the two-dimensional Kitaev model that the authors extend to the hyper-honeycomb lattice.","marker":"[29]"}],"fun_headline_variants":["Local spin oscillation reveals quantum spin liquid order","Dynamic order from local spin dynamics in spin liquid","Spin oscillation amplitude orders quantum spin liquid","Local spin rhythm marks the quantum spin liquid phase","Quantum spin liquid detected via dynamic spin order"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the quantum Monte Carlo data for a small finite system (N=4, with the transition temperature fixed by the arbitrary criterion |⟨Ŵ⟩|=0.03) represent the thermodynamic limit, since no finite-size scaling or statistical error analysis is provided to show the sharp oscillation onset persists as N approaches infinity.","fun_headline_variants_meta":{"raw":{"variants":["Local spin oscillation reveals quantum spin liquid order","Dynamic order from local spin dynamics in spin liquid","Spin oscillation amplitude orders quantum spin liquid","Local spin rhythm marks the quantum spin liquid phase","Quantum spin liquid detected via dynamic spin order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00033,"raw_usage":{"total_tokens":1774,"prompt_tokens":817,"completion_tokens":957,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":433,"completion_tokens_details":{"reasoning_tokens":890}},"tokens_in":433,"tokens_out":957,"duration_ms":7169,"temperature":1.0,"reasoning_tokens":890,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:51:45.032407+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A finite-size scaling study of the fitted amplitude Adyn(T) across N=4,5,6,7,8 with the same Wilson-loop criterion would settle whether the sharp onset of oscillation survives in the thermodynamic limit; if the fitted Adyn(T) does not sharpen toward a step function as N grows, or remains nonzero above Tw in the thermodynamic limit, the dynamic-order claim would be falsified.","supporting_citations":[{"cited_title":"Udagawa, Physical Review B 98, 220404 (2018)","cited_arxiv_id":null,"evidence_quote":"Gives the exact expression for the dynamical local spin correlation function in the two-dimensional Kitaev model that the authors extend to the hyper-honeycomb lattice."}],"review_version":1}