{"id":"b1c93c65-63cf-44e8-bba0-6b8bae1f422b","arxiv_id":"2506.14358","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Measurements show double-quantum spin relaxation in boron vacancies in hBN grows faster with temperature than single-quantum relaxation and is attributed to a high-energy phonon mode.","lead":"Researchers measured how fast the electron spin of a boron vacancy in hexagonal boron nitride (hBN) relaxes as temperature rises from 293 to 393 K, looking at two different spin transitions. They find the double-quantum transition relaxes much faster at high temperature and may set the limit for quantum sensing in this 2D material.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The high-energy-mode attribution in Eqs. (2)-(3) is an untested four-parameter fit to only five temperatures; without the reported Supporting Information, the \"five times at 400 K\" extrapolation is not established.","rationale":"The experimental work is careful: the pulse sequences in Fig. 2 implement the standard three-level rate-equation readouts, and the extracted room-temperature rates Ω = 33.26 ± 1.83 kHz and γ = 81.60 ± 9.12 kHz are consistent with the stated exponentials e^{-3Ωτ} and e^{-(2γ+Ω)τ}. I therefore do not put the burden on the raw measurement. The soft spot is the second half of the central claim, the quantitative mode attribution. Equations (2)-(3) are an ansatz borrowed from NV centers; no microscopic derivation for V_B^- in hBN is given, and the only support is that the resulting curves look reasonable. With four free parameters per rate and five temperatures, the fit is a curve through points; it cannot by itself identify the 165.75 meV mode as the cause. The missing Supporting Information removes the reader's ability to check this; the extrapolation to 400 K is a prediction of that unverified model. This is a correctable but important limitation, so the conditional verdict is appropriate; the concern does not overturn the measured temperature trends.","tokens_in":9744,"tokens_out":10923,"duration_ms":119125,"concrete_test":"Recover the raw F1(τ), F2(τ) decays (or author-provided decay curves) for all five temperatures and refit Ω(T), γ(T) with three nested models: (i) the full three-mode Eqs. (2)-(3); (ii) a two-mode version omitting the 165.75 meV term; (iii) a single effective mode plus a linear-in-T term. Report corrected AIC and 95% confidence intervals for each fitted coefficient. If the 165.75 meV coupling B_3 (or A_3) is not significant (ΔAIC < 2 or confidence interval includes zero), the paper's central attribution to the high-energy mode fails; if model (i) is decisively favored and B_3 is significantly positive, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II proposes Ω(T) = Σ_i A_i n_i(n_i+1) + A_S and γ(T) = Σ_i B_i n_i(n_i+1) + B_S (Eqs. 2-3), imported by analogy from the NV-in-diamond model [37]. The conclusion that the 165.75 meV mode dominates high-temperature relaxation depends on these fits being both valid and unique. The data set is only five temperatures (293, 318, 343, 368, 393 K) per spot; over this range the three Bose factors n_i(n_i+1) are smooth and strongly correlated, so the faster rise of γ can be absorbed by the high-energy term without being forced by the data. The fit coefficients, their uncertainties, and reduced χ² are not reported, and the Supporting Information (referenced for the derivation and for the 'higher energy, greater coupling coefficient' statement) is absent. The claim \"At 400 K, γ can reach the values five times of Ω\" is therefore a model extrapolation beyond the measured 393 K, not a measurement. Because the model collapses the full phonon spectral function onto three δ-peaks and ignores first-order processes and temperature-dependent couplings, the high-energy attribution is not falsified by the presented fit.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports measurements of the single-quantum (Ω) and double-quantum (γ) spin relaxation rates of negatively charged boron vacancies (V_B−) in hexagonal boron nitride over 293–393 K, using standard optical/microwave pulse sequences on two defect ensembles (Spot 1 and Spot 2). Both rates increase with temperature, with γ rising faster than Ω. The authors extrapolate that γ reaches five times Ω above 400 K. They fit the temperature dependence with a second-order spin-phonon model containing three effective phonon modes at 23.48, 77.39, and 165.75 meV, whose energies are taken from a DFT phonon density of states, plus sample-dependent constants. From the fit they conclude that the higher-energy phonon mode dominates high-temperature relaxation.","tokens_in":10047,"tokens_out":4370,"duration_ms":45284,"significance":"If the attribution holds, this is the first temperature-dependent double-quantum relaxation dataset for V_B− in hBN and identifies a regime where double-quantum relaxation may dominate spin-phonon decoherence. The experimental method is conventional and the two-spot reproducibility is a strength. However, the central mechanism claim currently rests on an imported NV-in-diamond model, an under-reported fit, and an extrapolation beyond the measured temperature range. The dataset itself is valuable and publishable, but the quantitative attribution needs additional transparency and falsifiability before it can be accepted as established.","major_comments":[{"comment":"The fits of Ω(T) and γ(T) use four parameters each (A1,A2,A3,A_S and B1,B2,B3,B_S) against only five temperatures per spot, but the manuscript gives no fitted parameter values, no parameter uncertainties, no residuals, and no reduced chi-square. Over 293–393 K the three Bose factors n_i(n_i+1) are smooth and strongly correlated, so the decomposition into three effective modes is not shown to be unique. Without fit-quality metrics, the statement in Section II that 'the higher energy phonon mode, the greater the coupling coefficient' and the resulting high-energy-mode attribution are not supported. The Supporting Information Section 4, cited for this result, is not available to the reader.","section":"Section II, Eqs. (2)-(3), Fig. 3(d)"},{"comment":"The claim 'At 400 K, γ can reach the values five times of Ω' is an extrapolation, not a measurement: the highest measured temperature is 393 K, and the model is only fitted over 293–393 K. The factor of five depends on the assumed model and on fitted coefficients whose uncertainties are not reported. This sentence should be explicitly labeled as a model extrapolation, or the claim should be supported by data at or above 400 K.","section":"Section II and Abstract"},{"comment":"The second-order spin-phonon model from NV centers in diamond [37] is transferred to V_B− in hBN by analogy, without independent validation. The model assumes only second-order processes with temperature-independent coupling coefficients and a three-peak phonon spectral function; it neglects first-order processes, local modes, temperature-dependent couplings, and magnetic noise. Since the high-energy-mode attribution is conditional on this specific model being valid, the manuscript should include at least a sensitivity test or a comparison with a simpler alternative (e.g., a one-mode model or a model allowing first-order terms) to show that the neglected contributions do not change the conclusion.","section":"Section II, 'Analogize to V_B spins'"},{"comment":"The phonon energies used in Eqs. (2)-(3) are taken from the DFT phonon density of states of a defective 3×3 monolayer hBN, while the experiment is performed on a mechanically exfoliated flake from bulk hBN. The manuscript does not discuss whether the monolayer phonon spectrum is representative of the multilayer/bulk sample. Because the phonon energies enter directly into the Bose factors and therefore control the fitted decomposition, this approximation needs an explicit justification or a sensitivity analysis.","section":"Fig. 3(c) and Section II"}],"minor_comments":[{"comment":"The manuscript contains numerous grammatical errors (e.g., 'is limited spin-phonon interactions' in the abstract, 'drives spin-lattice relaxation ... play a crucial role' in the Introduction); a careful proofread is needed.","section":"Throughout"},{"comment":"The figure does not show error bars for the relaxation rates, although the room-temperature values are quoted with uncertainties; without error bars it is difficult to judge whether the faster rise of γ relative to Ω is statistically significant.","section":"Fig. 3(d)"},{"comment":"The phrase 'first-principle calcutations' is a typo for 'first-principle calculations'.","section":"Fig. 3(c) caption"},{"comment":"The manuscript repeatedly cites Supporting Information Sections 2 and 4 for essential derivations and fitting details, but the Supporting Information is not included with the arXiv submission; the published version must include it for the claims to be verifiable.","section":"References to Supporting Information"},{"comment":"The spin sublevel notation is inconsistent ('ms' vs 'm_s'); it should be uniformly typeset as 'm_s' throughout the text and figures.","section":"Notation"}],"recommendation":"major_revision","confidential_remarks":"I recommend that the editor request the missing Supporting Information and the full fit results (parameter values, uncertainties, and goodness-of-fit) as part of the revision. If the fit metrics do not support a unique high-energy-mode decomposition, the conclusion should be weakened to a qualitative statement about γ growing faster than Ω. The novelty relative to Ref. [36], which already measured temperature-dependent single-quantum relaxation of V_B−, is mainly in the double-quantum channel; the authors should make this distinction explicit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe genuinely new thing here is the first temperature-dependent measurement of the double-quantum relaxation rate gamma for V_B- in hBN, from 293 to 393 K. The pulse sequences are standard, the exponential fits look clean, and two independent spots give consistent trends. The qualitative result — gamma rises faster with temperature than the single-quantum rate Omega — is credible from the raw data alone. If it holds, it changes how T1 limits are estimated for hBN-based sensors.\n\nThe soft spot is the theoretical attribution. Equations (2)-(3) import the NV-in-diamond second-order model, with three effective phonon modes, and the fit has eight free parameters: four coupling constants per rate, including sample offsets. There are only five temperatures per spot, and the three Bose factors n_i(n_i+1) are smooth and strongly correlated over this range. The fit cannot uniquely pin the 165.75 meV mode as the dominant high-temperature channel; it just absorbs the faster rise of gamma. No fit coefficient uncertainties or reduced chi-squared are reported, and the Supporting Information that supposedly contains the derivation, the fit details, and the \"higher energy, greater coupling coefficient\" statement is not included in the preprint. So the \"five times at 400 K\" line is a model extrapolation 7 K beyond the last data point, not a measurement, and the high-energy-mode attribution is not independently tested.\n\nI don't think this sinks the paper. The dataset stands on its own as a new measurement. The issue is framing: the abstract and conclusion lean harder on the model than the evidence supports. The load-bearing interpretive claim needs either more data (wider temperature range, more temperatures) or a genuinely predictive test, such as computing the coupling coefficients from first principles rather than fitting them.\n\nWho is this for? People working on spin-phonon limits in hBN and 2D van der Waals spin defects. It deserves a serious referee, not a desk reject: the experiment appears competently executed and the new observable is worth reporting. A referee should require the SI, the fit uncertainties, and a rewording of the high-temperature claim so it doesn't outrun the data.\n\nCandidly, [your name]","headline":"Useful new gamma(T) dataset for V_B- in hBN, but the model attribution outruns the five-temperature fit.","tokens_in":10630,"tokens_out":3115,"would_cite":true,"duration_ms":30596,"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":"For boron vacancies in hexagonal boron nitride, the double-quantum spin relaxation rate grows faster with temperature than the single-quantum rate, reaching five times it at 400 K and dominating high-temperature spin-phonon decoherence.","keywords":["hexagonal boron nitride","boron vacancy","double-quantum relaxation","single-quantum relaxation","spin-phonon interaction","quantum sensing","temperature dependence","spin coherence"],"falsifier":"Extend the relaxation measurements from 393 K to 450 K or higher on the same defect arrays. The model predicts $\\gamma$ will keep growing with the thermal occupation factor of the 165.75 meV mode and remain several times $\\Omega$; if $\\gamma$ saturates, crosses $\\Omega$, or the ratio drops, the claimed high-temperature dominance and its attribution to the high-energy phonon mode would be falsified. An independent check is to compute the coupling coefficients $A_i$ and $B_i$ from first principles and see whether the 165.75 meV mode actually has the largest values.","tokens_in":9580,"feed_emoji":"⚛️","tokens_out":11542,"duration_ms":107772,"temperature":0.7,"pith_summary":"This paper sets out to measure how the spin relaxation of negatively charged boron vacancies in hexagonal boron nitride (hBN) changes with temperature, separately for single-quantum transitions between the $m_s=0$ and $m_s=\\pm1$ levels and double-quantum transitions between $m_s=-1$ and $m_s=+1$. Over 293 to 393 K, both relaxation rates rise, but the double-quantum rate $\\gamma$ rises much faster than the single-quantum rate $\\Omega$, and the authors argue that by 400 K $\\gamma$ is about five times $\\Omega$. They attribute this behavior to second-order spin-phonon interactions with three effective phonon modes, with the highest-energy mode responsible for the rapid growth. If correct, the result identifies the dominant decoherence channel at high temperature and gives a concrete target for extending the coherence time of hBN spin sensors.","feed_headline":"In hBN, double-quantum spin relaxation dominates at high temperature","feed_subtitle":"Double-quantum spin decay hits five times the single-quantum rate at 400 K, setting the coherence limit for hBN sensing.","key_machinery":"The machinery is a second-order spin-phonon relaxation model, adapted from nitrogen-vacancy centers in diamond, plus a three-mode phonon spectral input. The relaxation rates are written as $\\Omega(T)=\\sum_{i=1,2,3} A_i n_i(n_i+1)+A_S$ and $\\gamma(T)=\\sum_{i=1,2,3} B_i n_i(n_i+1)+B_S$, where $n_i=(e^{\\hbar\\omega_i/k_B T}-1)^{-1}$ is the thermal occupation number of the $i$-th effective phonon mode, $A_i,B_i$ are temperature-independent coupling coefficients, and $A_S,B_S$ absorb sample-specific constants. The phonon energies come from a density-functional-theory phonon density of states, whose three peaks at 23.48, 77.39 and 165.75 meV are taken as the modes that dominate relaxation. The experimental identification is carried by two fluorescence pulse sequences whose decays are governed by $3\\Omega$ and $2\\gamma+\\Omega$, allowing the two rates to be separated.","core_discovery":"The central claim is that the double-quantum relaxation rate $\\gamma$ of the $V_\\mathrm{B}^-$ center in hBN grows with temperature faster than the single-quantum rate $\\Omega$, and that this growth is governed by second-order spin-phonon interactions with three effective phonon modes. Using two pulse sequences, the authors isolate decays at rates $3\\Omega$ and $2\\gamma+\\Omega$, and from fits they extract $\\Omega$ and $\\gamma$ from 293 to 393 K in two defect ensembles. They find both rates increase with temperature, with $\\gamma$ rising rapidly enough that at 400 K it reaches five times $\\Omega$. Applying a model of the form $\\Omega(T)=\\sum_i A_i n_i(n_i+1)+A_S$ and $\\gamma(T)=\\sum_i B_i n_i(n_i+1)+B_S$, with thermal occupation numbers $n_i$ for phonon energies 23.48, 77.39, and 165.75 meV, the authors reproduce the measurements and find that the higher-energy phonon mode has the larger coupling coefficients. They conclude that double-quantum relaxation may dominate the spin-phonon decoherence channel at high temperature, so the total spin-lattice relaxation time obeys $1/T_1 = 3\\Omega + \\gamma$ with $\\gamma$ setting the limit.","pith_inferences":["Beyond the paper: the fitted coupling coefficients are not uniquely pinned by the data unless the model's functional form is right; an independent first-principles calculation of the coefficients, or a measurement at temperatures far above 393 K, is needed to confirm that the 165.75 meV mode is the true driver.","Beyond the paper: a direct extension above 400 K would turn the extrapolated five-fold ratio into a measured fact; if $\\gamma$ continues to follow the thermal occupation factor of the high-energy mode, the claimed dominance is strengthened, and if not, the model would need revision.","Beyond the paper: if double-quantum relaxation dominates, coherence stored in the $m_s=\\pm1$ subspace is especially fragile; strategies that reshape the phonon environment, such as strain engineering or encapsulation, become a natural path to longer coherence times.","Beyond the paper: the sample-related constants $A_S$ and $B_S$ may be absorbing non-phonon relaxation such as magnetic noise; comparing samples with different defect densities or isotopic compositions would separate intrinsic phonon effects from ensemble artifacts."],"forward_implications":["At high temperature, the $|m_s=-1\\rangle$ to $|m_s=+1\\rangle$ double-quantum channel, not the single-quantum channel, is the main spin-phonon decoherence path for the boron vacancy in hBN.","Because the total spin-lattice relaxation time obeys $1/T_1 = 3\\Omega + \\gamma$, a double-quantum rate several times $\\Omega$ makes $\\gamma$ the quantity that sets the high-temperature $T_1$ limit.","The highest-energy effective phonon mode (165.75 meV) carries the largest coupling coefficients, so high-temperature relaxation is driven by high-frequency lattice vibrations rather than by the lowest phonon mode.","The same second-order relaxation model, applied to other spin defects in hBN, would predict which quantum channel limits their coherence at elevated temperatures."],"supporting_citations":[{"why":"Supplies the second-order spin-phonon relaxation model and its functional form, which the paper adapts to the boron vacancy.","marker":"[37]"},{"why":"Provides the triple-peak phonon vibration spectral function and earlier temperature-dependent spin-phonon data for hBN boron vacancies.","marker":"[36]"},{"why":"Supplies the double-quantum measurement pulse sequences and the relation between total relaxation and the two rates used to extract $\\Omega$ and $\\gamma$.","marker":"[40]"},{"why":"Supports the phonon mode energies used as effective modes in the relaxation model through spin-phonon calculations for boron-vacancy centers in boron nitride.","marker":"[44]"}],"fun_headline_variants":["Double-quantum spin relaxation outpaces single in hBN at high T","hBN spin decay: double-quantum rate dominates above 400 K","At high T, double-quantum relaxation governs hBN spin coherence","Double-quantum spin relaxation sets the limit in hBN at high T"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the boron vacancy relaxes by the same two-phonon process used for nitrogen-vacancy centers in diamond, with three fixed vibration modes and temperature-independent coupling strengths, and that no other relaxation mechanism such as magnetic noise or local defects contributes significantly.","fun_headline_variants_meta":{"raw":{"variants":["Double-quantum spin relaxation outpaces single in hBN at high T","hBN spin decay: double-quantum rate dominates above 400 K","At high T, double-quantum relaxation governs hBN spin coherence","Double-quantum spin relaxation sets the limit in hBN at high T"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000564,"raw_usage":{"total_tokens":2705,"prompt_tokens":1005,"completion_tokens":1700,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":621,"completion_tokens_details":{"reasoning_tokens":1618}},"tokens_in":621,"tokens_out":1700,"duration_ms":11823,"temperature":1.0,"reasoning_tokens":1618,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:17:03.774860+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Extend the relaxation measurements from 393 K to 450 K or higher on the same defect arrays. The model predicts $\\gamma$ will keep growing with the thermal occupation factor of the 165.75 meV mode and remain several times $\\Omega$; if $\\gamma$ saturates, crosses $\\Omega$, or the ratio drops, the claimed high-temperature dominance and its attribution to the high-energy phonon mode would be falsified. An independent check is to compute the coupling coefficients $A_i$ and $B_i$ from first principles and see whether the 165.75 meV mode actually has the largest values.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the triple-peak phonon vibration spectral function and earlier temperature-dependent spin-phonon data for hBN boron vacancies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the double-quantum measurement pulse sequences and the relation between total relaxation and the two rates used to extract $\\Omega$ and $\\gamma$."},{"cited_title":"Estaji, I","cited_arxiv_id":null,"evidence_quote":"Supports the phonon mode energies used as effective modes in the relaxation model through spin-phonon calculations for boron-vacancy centers in boron nitride."}],"review_version":1}