{"id":"eeec99fc-4344-492d-86f2-5cac1ff21115","arxiv_id":"2608.07697","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A frequency-domain continuous-wave protocol extracts NV spin-lattice relaxation times from lock-in ODMR response, validated on bulk and nanodiamond samples and applied to Mn2+ sensing.","lead":"This paper introduces a continuous-wave, frequency-domain method that extracts NV diamond spin-lattice relaxation times from lock-in ODMR response, replacing slow pulsed sequences with a fast measurement. If it holds up, it gives nanodiamond bio-sensing and long-T1 qubit characterization a much faster route to quantitative relaxometry.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (5)-(6) assume additive linear rates; if cross-terms or saturation leak into the extrapolation window, the zero-power intercept is not the dark T1.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the additive linear extrapolation in Eq. (6) is asserted without a derivation, and if cross-terms or saturation leak into the fitted window, the extracted intercept is not the dark T1. This is the single most load-bearing issue because the entire method reduces to that extrapolation: the frequency response only yields Teff, and T1 is obtained by the linear zero-power limit. The empirical validation against pulsed TDR is genuine supporting evidence, but it is limited to two bulk samples and the IN3x3 comparison is not clean (the paper itself attributes the discrepancy to a suboptimal π-pulse). The abstract's unsupported '60 μs' endpoint and the 'parameter-free' overstatement are secondary. A conditional verdict is appropriate: the method is plausible and internally consistent, but the central formula needs either a full derivation or a two-dimensional cross-term test. Our read does not change the reader's CONDITIONAL verdict, so the recommended verdict remains unchanged.","tokens_in":13978,"tokens_out":13867,"duration_ms":137117,"concrete_test":"Build a rate-equation model of the NV- ground/excited triplet and metastable singlet, and compute the slowest eigenvalue as a function of P_MW and I. If the expansion contains a P_MW×I term, Eq. (6)'s additive form is invalid. Experimentally, measure 1/Teff vs I at multiple MW powers and fit 1/Teff = 1/T1 + αP_MW + βI + γP_MW I; a nonzero γ at 2σ in the fitted window would require a 2D extrapolation. Compare the corrected T1 with pulsed TDR.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central extraction of T1 rests on Eq. (6), 1/Teff ≈ 1/T1 + αP_MW + βI, which is asserted without a derivation: the preceding Bloch equation (1) yields only the MW term 2P, and the optical term Γp is added as 'a rate-equation treatment yields, to leading order' (text before Eq. (5)). No microscopic model is given, and the identification of Γp with the stimulated-emission rate σλI/hν is physically questionable—NV optical polarization proceeds via spin-selective intersystem crossing, not stimulated emission. If the full rate dynamics contain a P_MW×I cross-term (e.g., MW-induced changes in the effective optical pumping cycle, or light-induced broadening modifying the MW drive), then a linear extrapolation in I at fixed P_MW gives an intercept 1/T1 + f(P_MW). The FDR vs TDR agreement in two bulk samples is supportive, but the IN3x3 values (4.75±0.14 ms vs 4.15±0.30 ms) differ by about two combined errors, and the paper's explanation (suboptimal π-pulse) is plausible but not verified. Thus the method's central claim would fail if the additive linear model is not exact in the fitted window.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces a continuous-wave frequency-domain relaxometry (FDR) protocol for measuring the spin-lattice relaxation time T1 of NV centers. Microwave amplitude is modulated at frequency ωm and the lock-in X and Y components of the ODMR signal are recorded; the response is modeled as a first-order low-pass filter with time constant T_eff. The authors claim that extrapolating 1/T_eff to zero microwave power and zero laser intensity yields the intrinsic T1, and they validate this against pulsed time-domain relaxometry (TDR) in two bulk diamond samples, measure a cryogenic T1 of about 200 ms in a DNV sample, and apply the method to nanodiamond ensembles of three sizes, including detection of Mn2+ ions in water. The paper also reports substantial measurement speed-ups and lower optical intensities relative to confocal pulsed relaxometry.","tokens_in":14241,"tokens_out":4975,"duration_ms":49054,"significance":"If the central claim holds, FDR would be a genuinely useful addition to quantum relaxometry: it avoids fast pulse switching, works in the frequency domain, and promises access to T1 values that are difficult to measure with pulsed protocols. The bulk validation, the consistency between the X and Y components, the phonon-like temperature dependence, and the demonstrated speed advantage in nanodiamonds are concrete strengths. However, the method's quantitative foundation is currently incomplete: the key additive-rate relation is asserted rather than derived, the saturated-regime analysis used for the sensing application is not validated, and the abstract claims a dynamic range that does not appear in the data. These issues are load-bearing for the paper's central claims, so the manuscript needs substantial revision before the results can be fully trusted.","major_comments":[{"comment":"The central extrapolation rests on Eq. (6), 1/T_eff ≈ 1/T1 + αP_MW + βI, but the text only states that \"a rate-equation treatment yields, to leading order\" and does not present that treatment. Equation (1) is a two-level Bloch equation containing only the microwave-induced rate 2P and T1; the optical contribution Γp is then introduced as σλI/(hν) (stimulated emission), which is not the standard NV optical polarization mechanism (spin-selective intersystem crossing). Because a P_MW×I cross-term or a nonlinearity within the fitted window would change the zero-power intercept, this missing derivation is load-bearing. The authors should supply a full rate-equation derivation, justify the form of Γp for NV centers, and ideally test the linearity assumption explicitly over the extrapolation window.","section":"Theoretical background of FDR, Eqs. (5)-(6)"},{"comment":"The abstract claims T1 measurements \"directly demonstrated from 60 μs to 200 ms,\" implying a range of more than three orders of magnitude. However, the main text and the Supplementary Information report no measurement at 60 μs; the shortest quoted values are about 0.7 ms for ND50 and the cryogenic value is about 200 ms. The dynamic-range claim should either be supported with data (for example, a fast-relaxing sample or a fitted low-frequency response) or removed and reworded to reflect what was actually measured.","section":"Abstract and Results"},{"comment":"The Mn2+ sensing experiments rely on one-shot measurements at +15 dBm, where the text explicitly states the system is \"not in the linear approximation regime anymore\" and the Y-signal shape deviates from Eq. (S1). The paper then asserts that \"the position of the signal maximum still correctly reflects the relaxation contribution and can be used in relative measurements.\" This assertion is neither derived nor validated against an independent reference. Since the quantitative sensitivity factors (about 5 for ND70 and 6.3 for ND50) depend on this saturated-regime assumption, the authors should validate it, for instance by comparing Y-maximum shifts with a full T_eff extraction or with pulsed TDR on the same samples.","section":"Supplementary Note 2 and Application of FDR to nanodiamond samples"},{"comment":"The claimed \"excellent agreement\" between FDR and TDR is only partial. For IN3x3, FDR gives T1 = (4.75 ± 0.14) ms and TDR gives T1 = (4.15 ± 0.30) ms, which differ by about two combined errors. The explanation that this reflects suboptimal π-pulse settings is plausible but is not verified. Given that the validation rests on only two bulk samples, the discrepancy should either be resolved by improved TDR calibration in IN3x3 or be presented more cautiously as a systematic uncertainty in the comparison.","section":"Validation of FDR in bulk samples, Fig. 2(f)"},{"comment":"The nanodiamond T1 values (1.6 ms, 0.8 ms, and 0.7 ms for ND100, ND70, and ND50) are reported without uncertainties or error bars, even though the text notes that \"errors in the T1 values significantly increase at low MW and laser powers.\" Quantitative relaxometry-based sensing and the size-dependent comparisons require error estimates for these extrapolated values, especially since the method is proposed as a quantitative metrological tool.","section":"Application of FDR to nanodiamond samples, Fig. 4"}],"minor_comments":[{"comment":"The text says the effective relaxation time can be determined \"from equation (7)\", but the manuscript contains no numbered equation (7); the reference should be corrected to the appropriate equation or a new equation should be numbered.","section":"Validation of FDR in bulk samples"},{"comment":"The quantity Γp = σλI/(hν) is called the polarization rate of stimulated emission. Even after a full derivation is supplied, the wording should distinguish the net optical pumping rate from the stimulated-emission cross-section contribution, as the two are not the same in NV centers.","section":"Theoretical background of FDR, Eq. (5)"},{"comment":"The conclusion describes the protocol as \"parameter-free,\" which is inconsistent with the need for the slopes α and β in Eq. (6) and the associated extrapolation procedure; the phrasing should be revised.","section":"Conclusion"},{"comment":"There are several typographical errors that should be corrected: \"timeT1\" in the introduction, \"estimsted\" in Supplementary Note 2, \"MV power\" in Supplementary Note 3, and \"spectum\" in the Supplementary Figure 1 caption.","section":"Throughout manuscript"},{"comment":"The text states that the linear approximation to zero laser power is shown by solid red lines, but some displayed points are in the sub-linear regime; the manuscript should clarify which points were included in each linear fit and how the fit range was chosen.","section":"Validation of FDR in bulk samples, Figs. 2(c,d)"},{"comment":"In Eq. (1), n0 is defined as the population immediately after optical polarization, whereas T1 relaxation normally refers to approach to thermal equilibrium; the definition should be made consistent throughout the derivation.","section":"Theoretical background of FDR, Eq. (1)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and addresses a topic of clear interest to the quantum-sensing and NV-center community. The main concern is that the core rate-equation model is asserted rather than derived, and the saturated-regime measurements used for the nanodiamond sensing claims lack validation. I also recommend that the authors be asked to reconcile the abstract's 60 μs claim with the actual data. There is no indication of any ethical or attribution problem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid empirical methods paper that does something genuinely new — extracting ground-state T1 from the frequency response of cw-ODMR under microwave amplitude modulation. The validation against pulsed TDR in two bulk diamonds is real, and the X/Y component consistency is a nice check. The method deserves referee time, but the theoretical underpinning is thinner than it looks.\n\nWhat's new: frequency-resolved lock-in spectroscopy has been used for photoluminescence lifetimes, and cw-ODMR modulation is standard for magnetometry, but using the low-frequency roll-off to pull out the spin-lattice relaxation time, then extrapolating to zero MW and laser power, is a new application. The two-bulk-sample agreement with pulsed TDR (4.75 vs 4.15 ms for IN3x3, 5.55 vs 6.27 for DNV) supports the central claim, as does the cryogenic behavior and the nanodiamond relative sensing with water and Mn2+. The speedup claims — minutes vs hours for long T1 — look credible. No data or code shipped, though the paper is transparent enough about the fitting.\n\nSoft spots, in proportion. The biggest one is Eq. (5)–(6): the paper says 'a rate-equation treatment yields, to leading order' 1/Teff ≈ 1/T1 + 2P + Γp, but never shows the treatment. The identification of Γp with the stimulated-emission rate is physically odd — NV optical polarization runs through intersystem crossing, not stimulated emission. And the assumed additivity of MW and optical rates is exactly what could break if there are cross-terms. The data are reassuring: the DNV intercepts are stable across three MW powers, and the IN3x3 discrepancy is plausibly explained by the short T2* spoiling the pi-pulse. So I'm not claiming the method fails — but the derivation needs to actually be worked out before the paper is fully convincing. Minor: the abstract's '60 µs' lower endpoint doesn't appear in the body; and 'parameter-free' overstates things, since slopes and intercepts are fitted. The nanodiamond +15 dBm sensing operates outside the linear regime, though that's disclosed and used only for relative contrast.\n\nBottom line: the idea is new, the empirical case is solid, and the limitations are mostly acknowledged. The missing rate-equation derivation is fixable without invalidating the results. I'd send this to peer review — the referees should ask for the theory, not reject on the current evidence. I'd bring it to a reading group.","headline":"Genuinely new cw-ODMR method for T1 relaxometry with real validation, but the missing rate-equation derivation for the load-bearing extrapolation needs to be supplied before the theory is fully convincing.","tokens_in":14786,"tokens_out":2295,"would_cite":true,"duration_ms":22541,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["76.30.Mi","76.70.Hb"],"model":"deepseek-v4-flash","headline":"The paper claims that $T_1$ relaxation times of nitrogen-vacancy centers can be extracted from the frequency response of cw-ODMR under low-frequency microwave amplitude modulation, turning relaxometry from a pulsed time-domain measurement…","keywords":["quantum relaxometry","nitrogen-vacancy centers","T1 relaxation","cw-ODMR","frequency-domain spectroscopy","lock-in detection","nanodiamond sensing","spin-lattice relaxation"],"falsifier":"On a bulk NV sample with independently calibrated pulsed $T_1$, measure $T_\\mathrm{eff}$ on a grid of microwave and laser powers, then halve the lowest powers and remeasure; if the linear extrapolation intercept shifts by more than the combined uncertainty, the additivity model fails. A numerical solution of the Bloch equations with the actual square-wave amplitude modulation should reproduce the measured $X$ and $Y$ curves with a single low-pass mode; needing extra modes or a power-dependent $T_1$ to fit would disprove the single-$T_\\mathrm{eff}$ description.","tokens_in":13721,"feed_emoji":"⏱️","tokens_out":9723,"duration_ms":85003,"temperature":0.7,"pith_summary":"The paper seeks to establish that the spin-lattice relaxation time $T_1$ of nitrogen-vacancy centers in diamond can be measured from the frequency response of a continuous-wave optically detected magnetic resonance (cw-ODMR) signal under low-frequency microwave amplitude modulation, with no pulsed sequence and no time-resolved detection. The central move is to treat the spin system as a first-order low-pass filter: the lock-in in-phase and quadrature components roll off with an effective time constant $T_\\mathrm{eff}$, and $1/T_\\mathrm{eff}$ is argued to equal $1/T_1$ plus terms linear in microwave and laser power, so extrapolating to zero power recovers the dark $T_1$. If correct, this turns relaxometry into a fast frequency-domain measurement that spans 60 $\\mu$s to 200 ms in one protocol, works on inhomogeneous nanodiamond ensembles where $\\pi$ pulses are impractical, and cuts acquisition time by orders of magnitude compared with pulsed relaxometry.","feed_headline":"Frequency sweep replaces pulsed T1 measurement","feed_subtitle":"Spin-lattice relaxation from 60 µs to 200 ms comes from a lock-in response, no pulse sequences required.","key_machinery":"The central object is the effective longitudinal time constant $T_\\mathrm{eff}$ of the driven spin system, extracted from lock-in detection of cw-ODMR. The identity $1/T_\\mathrm{eff}\\approx 1/T_1+2P+\\Gamma_p$, with $P\\propto P_\\mathrm{MW}$ and $\\Gamma_p\\propto I$, is the load-bearing relation: it converts a set of continuous-wave measurements at different microwave and laser powers into a linear extrapolation whose intercept is the intrinsic relaxation rate $1/T_1$. The frequency response is read either from the half-decay point of the in-phase component $X(\\omega_m)$ or from the maximum of the quadrature component $Y(\\omega_m)$, which occurs at $f_m=1/(2\\pi T_\\mathrm{eff})$.","core_discovery":"Under resonant continuous-wave drive, the NV ground-state transition behaves as an effective two-level system whose population response to microwave amplitude modulation is a first-order low-pass filter. The lock-in components are $X(\\omega_m)\\propto 1/(1+(\\omega_m T_\\mathrm{eff})^2)$ and $Y(\\omega_m)\\propto -\\omega_m T_\\mathrm{eff}/(1+(\\omega_m T_\\mathrm{eff})^2)$, so $T_\\mathrm{eff}$ is fixed by the roll-off or by the quadrature maximum at $f_m=1/(2\\pi T_\\mathrm{eff})$. To leading order, $1/T_\\mathrm{eff}\\approx 1/T_1+2P+\\Gamma_p$, with $P\\propto P_\\mathrm{MW}$ and $\\Gamma_p\\propto I$, which justifies a two-dimensional linear extrapolation to zero microwave and laser power for the intrinsic dark $T_1$. The paper validates the protocol on bulk ensembles, obtaining $T_1=(4.75\\pm0.14)$ ms and $(5.55\\pm0.20)$ ms, consistent with pulsed time-domain relaxometry, and a cryogenic value near 200 ms; it then reports intrinsic nanodiamond $T_1$ values of 1.6 ms, 0.8 ms, and 0.7 ms and uses relative $T_1$ changes to detect micromolar Mn$^{2+}$ ions.","pith_inferences":["If the linear-additivity relation survives wider tests, the same frequency-domain readout should transfer to other long-lived ground-state spin defects, such as divacancy or transition-metal defects in silicon carbide, because the low-pass response is generic to driven two-level systems.","The nanodiamond results suggest that some previously reported short $T_1$ values for small nanodiamonds may reflect pulse imperfections rather than intrinsic relaxation; a direct test is to measure the same particles with FDR and an optimized pulsed sequence.","The water-induced $T_1$ increase seen in small nanodiamonds implies that hydration state, not only surface chemistry, determines effective relaxation in biological settings, so FDR's speed could enable real-time hydration or osmolarity studies."],"forward_implications":["$T_1$ can be measured without fast optical and microwave switching, $\\pi$-pulse calibration, or dark evolution intervals, and one frequency sweep covers relaxation times from 60 $\\mu$s to 200 ms.","In bulk samples, FDR and pulsed time-domain relaxometry agree within error, with two to four times smaller uncertainty in about six times less acquisition time.","Nanodiamond ensembles with randomly oriented NV axes, where reliable $\\pi$-pulses are impractical, become measurable; reported intrinsic $T_1$ values of 1.6, 0.8, and 0.7 ms are among the longest reported for such particles.","Relaxometry-based sensing of paramagnetic species can be done by relative $T_1$ changes in one-shot measurements: ND70 shows about a fivefold $T_1$ reduction at 500 $\\mu$M Mn$^{2+}$ with a high-SNR acquisition in roughly 2 minutes.","At cryogenic temperatures, where pulsed measurements are impractical, FDR gives $T_1$ near 200 ms."],"supporting_citations":[{"why":"Supplies the frequency-resolved spectroscopy framework that the FDR method transfers to ground-state spin relaxation.","marker":"[17]"},{"why":"Provides the Bloch-equation population dynamics from which the first-order low-pass response and $T_\\mathrm{eff}$ are derived.","marker":"[20]"},{"why":"Gives the stimulated emission cross-section making the optical polarization rate $\\Gamma_p\\propto I$ in the extrapolation model.","marker":"[21]"},{"why":"Defines the pulsed $T_1$ protocol whose values are used as the benchmark for bulk FDR validation.","marker":"[7]"},{"why":"Identifies the phonon contribution that explains the low-temperature plateau assigned to spin-lattice relaxation.","marker":"[22]"},{"why":"Supplies the in-solution nanodiamond relaxometry comparison for measurement speed and method limitations.","marker":"[25]"},{"why":"Provides the state-of-the-art Gd$^{3+}$ relaxometry sensitivity benchmark and the confocal laser-intensity values compared in Supplementary Table 1.","marker":"[31]"}],"fun_headline_variants":["CW relaxometry reads T1 from lock-in roll-off","No pulses: T1 from CW frequency response","Frequency-domain T1 from 60 us to 200 ms","Lock-in T1 without pulsed sequences","CW relaxometry spans 60 us to 200 ms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that Eq. (5) holds to leading order in the fitted window, so the measured rate satisfies $1/T_\\mathrm{eff}\\approx 1/T_1+2P+\\Gamma_p$ with $P$ linear in microwave power and $\\Gamma_p$ linear in laser intensity; the paper asserts this rate-equation result without deriving it, and if these contributions do not add linearly, the zero-power intercept is not the true dark $T_1$.","fun_headline_variants_meta":{"raw":{"variants":["CW relaxometry reads T1 from lock-in roll-off","No pulses: T1 from CW frequency response","Frequency-domain T1 from 60 us to 200 ms","Lock-in T1 without pulsed sequences","CW relaxometry spans 60 us to 200 ms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002149,"raw_usage":{"total_tokens":8382,"prompt_tokens":1041,"completion_tokens":7341,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":7265}},"tokens_in":657,"tokens_out":7341,"duration_ms":50927,"temperature":1.0,"reasoning_tokens":7265,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:23:36.103748+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On a bulk NV sample with independently calibrated pulsed $T_1$, measure $T_\\mathrm{eff}$ on a grid of microwave and laser powers, then halve the lowest powers and remeasure; if the linear extrapolation intercept shifts by more than the combined uncertainty, the additivity model fails. A numerical solution of the Bloch equations with the actual square-wave amplitude modulation should reproduce the measured $X$ and $Y$ curves with a single low-pass mode; needing extra modes or a power-dependent $T_1$ to fit would disprove the single-$T_\\mathrm{eff}$ description.","supporting_citations":[{"cited_title":"Philosophical Magazine B 50(5), 579–597 (1984)","cited_arxiv_id":null,"evidence_quote":"Supplies the frequency-resolved spectroscopy framework that the FDR method transfers to ground-state spin relaxation."},{"cited_title":"Principles of Magnetic Resonance, pp","cited_arxiv_id":null,"evidence_quote":"Provides the Bloch-equation population dynamics from which the first-order low-pass response and $T_\\mathrm{eff}$ are derived."},{"cited_title":"Nature communications8(1), 14000 (2017)","cited_arxiv_id":null,"evidence_quote":"Gives the stimulated emission cross-section making the optical polarization rate $\\Gamma_p\\propto I$ in the extrapolation model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the pulsed $T_1$ protocol whose values are used as the benchmark for bulk FDR validation."},{"cited_title":"Physical Review Letters130(25), 256903 (2023) 15","cited_arxiv_id":null,"evidence_quote":"Identifies the phonon contribution that explains the low-temperature plateau assigned to spin-lattice relaxation."},{"cited_title":"Physical Review Applied20(3), 034018 (2023)","cited_arxiv_id":null,"evidence_quote":"Supplies the in-solution nanodiamond relaxometry comparison for measurement speed and method limitations."},{"cited_title":"ACS sensors5(12), 3862–3869 (2020)","cited_arxiv_id":null,"evidence_quote":"Provides the state-of-the-art Gd$^{3+}$ relaxometry sensitivity benchmark and the confocal laser-intensity values compared in Supplementary Table 1."}],"review_version":1}