{"id":"b9145e48-dd72-43e3-a4f4-e34efa6db54c","arxiv_id":"1908.06829","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In 4H-SiC, the silicon vacancy spin relaxation time T1 stays near 300 microseconds across irradiation types and fluences, while the spin coherence time T2 decreases with emitter density according to an empirical two-parameter scaling formula.","lead":"This paper measures how long the spin of a silicon vacancy defect in silicon carbide stays coherent after the crystal is hit by neutrons, electrons, or protons. It finds that the spin relaxation time stays nearly constant across a wide range of defect densities, while the coherence time shrinks as defect density grows, and it proposes a simple two-parameter formula to describe the trade-off.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Universal scaling rests on the untested assumption that 1/T2 decoherence is linear in PL-derived emitter density; since κ varies by irradiation type and neutron T2 data are lower bounds, Eq. (1) is currently a parametrization rather than an established predictive law.","rationale":"The reader's weakest-assumption diagnosis—that 1/T2^(i) proportional to NV is assumed without independent evidence—is the same concern I find most load-bearing. The paper's own numbers strengthen the concern: κ varies by more than an order of magnitude between electron and neutron irradiation and even between the two neutron reactors, so the 'same equation with two independent parameters' is only a flexible fitting form. The self-admitted lower-bound nature of neutron T2 data adds further uncertainty to the extracted κ values, though this particular artifact would, if corrected, likely move neutron points toward longer T2 and possibly toward the electron curve. My examination therefore does not change the conditional verdict. The qualitative conclusions—T1 robust to irradiation type, T2 degrades with radiation damage, electron irradiation is gentler than neutron irradiation—are supported by the data, and the model in Eq. (1) is a reasonable phenomenological starting point. What is not supported is the stronger 'universal scaling' claim, which should be softened or backed by an independent measure of total defect density. Since the reader already recommended CONDITIONAL, I recommend no verdict change.","tokens_in":15376,"tokens_out":5771,"duration_ms":65439,"concrete_test":"Perform electron paramagnetic resonance (EPR) or another bulk spin-counting measurement on the same electron- and neutron-irradiated 4H-SiC samples used for Fig. 5b to obtain the total paramagnetic defect density N_tot. Then fit 1/T2 - 1/T2_pristine versus N_tot across all irradiation types with a single slope. If a single-slope linear fit fails while per-type fits to NV succeed, the apparent universal scaling is an artifact of the PL proxy; if the single-slope fit succeeds, the linearity assumption is validated and the universal form is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that T2 follows T2 = T2_pristine / (1 + T2_pristine κ NV) universally, with the same equation and two independent parameters describing all irradiation data. The load-bearing step is the proportionality 1/T2^(i) = κ NV, stated in the Discussion: 'We assume that 1/T2^(i) is proportional to the number of all irradiation induced defects, which in turn scales with NV.' No independent measurement of the total irradiation-induced defect density is provided; NV is calibrated from photoluminescence of optically active VSi centers. This proxy is known to break down at high fluence, where PL drops because non-radiative defects dominate while total damage continues to increase—indeed, the authors report a deviation from Eq. (1) for NV > 7×10^15 cm^-3. Furthermore, κ is not universal: the fits give κ = 0.8×10^-11 s^-1 cm^3 for electrons, κ = 3.8×10^-11 and 25.2×10^-11 s^-1 cm^3 for the two neutron reactors, and the authors concede κ should depend on ion type, energy, and depth. Thus Eq. (1) with a different κ per condition is a definition of κ, not a predictive universal scaling. The neutron T2 values are additionally self-admitted lower bounds due to microwave field inhomogeneity, so even the reported κ values are not firmly established. The design rule 'know NV and irradiation type, predict T2' therefore relies on an untested monotonic relationship between PL-derived NV and the total defect spin bath.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a systematic pulsed ODMR study of T1 and T2 for silicon vacancies in 4H-SiC after neutron, electron, and proton irradiation over wide fluence ranges, including annealing experiments. The central claim is that T2 follows a universal scaling with the optically active emitter density NV, T2 = T2_pristine / (1 + T2_pristine kappa NV), with T2_pristine and kappa as the only parameters, and that T1 remains near 300 microseconds up to NV ≈ 7x10^15 cm^-3. The data include new measurements and literature points, and the paper proposes this scaling as a design rule for choosing irradiation conditions.","tokens_in":15712,"tokens_out":5091,"duration_ms":44192,"significance":"If the proposed scaling were established, it would provide a practically useful rule for estimating T2 from NV and irradiation type. The paper's dataset is valuable: it compares multiple irradiation particles and energies, demonstrates T1 robustness, and shows annealing can partially recover T2. However, the central universal-scaling claim currently rests on fitting Eq. (1) with irradiation-type-dependent kappa values and an untested proportionality between decoherence rate and the PL-derived NV, so the paper's significance is diminished until the scaling is substantiated with independent calibrations and uncertainty-aware fits.","major_comments":[{"comment":"The claimed universal scaling is not supported because kappa is not a single parameter across the data. The fits use kappa = 0.8x10^-11, 3.8x10^-11, and 25.2x10^-11 s^-1 cm^3 for electron, BER II neutron, and TRIGA neutron irradiations, respectively. With a different kappa per irradiation condition, Eq. (1) becomes a parametrization, and the concluding statement that 'all data can be well described by the same equation with two independent parameters' is misleading because the two parameters are re-fit per data set. To claim universality, the authors would need to show kappa is determined by known irradiation parameters (e.g., particle type, energy, or displacement cross-section) and to test the scaling form with at least one parameter-free prediction.","section":"Discussion, Eq. (1); Conclusion"},{"comment":"This assumption is the load-bearing premise of Eq. (1), but it is not tested. NV is calibrated from PL intensity of optically active VSi centers (Results, Fig. 2), not from the total irradiation-induced defect concentration. The authors themselves note that PL drops at high fluence because non-radiative defects dominate while total damage increases; indeed, they report a deviation from Eq. (1) for NV > 7x10^15 cm^-3. In the absence of an independent measure of total defect density, the linear relation 1/T2^(i) = kappa NV is an unverified proxy, and the model's claimed predictive power is not established. The manuscript should either provide independent evidence for this proportionality (e.g., from EPR or absorption measurements) or explicitly limit the claim to the regime where PL-derived NV is monotonic in fluence.","section":"Discussion: 'We assume that 1/T2^(i) is proportional to the number of all irradiation induced defects...'"},{"comment":"The authors state that coherence times of neutron-irradiated samples 'should be seen as a lower bound' because of microwave-field inhomogeneity. If T2 is a lower bound, then the fitted kappa values for neutron data (3.8x10^-11 and 25.2x10^-11 s^-1 cm^3) are not reliable: the true T2 may be longer and kappa correspondingly smaller. Fitting lower-bound data to extract kappa and then comparing these kappa values with electron data biases the conclusion that neutron irradiation produces a larger decoherence per emitter. These points should either be corrected with homogeneous-pulse measurements or excluded from quantitative kappa extraction, with the qualitative trend stated separately.","section":"Discussion, neutron-irradiated samples"},{"comment":"The manuscript does not report uncertainties for the fitted parameters T2_pristine and kappa, and many T2 data points in Fig. 5b lack error bars. Without these, the statement that the data are 'well described' by Eq. (1) cannot be assessed, and the comparison of kappa values across irradiation types is not meaningful. Please report standard errors or confidence intervals for all fitted parameters, and show fit residuals or include the number of data points and goodness-of-fit metrics. The same applies to the T1 = 300 microseconds plateau, for which only representative error bars are shown.","section":"Results, Fig. 5; Methods"}],"minor_comments":[{"comment":"The abbreviation 'ZFL' in 'zero-phonon line (ZFL)' should be 'ZPL'.","section":"Introduction and Results"},{"comment":"'varyable delay time' should be 'variable delay time'.","section":"Fig. 4c caption"},{"comment":"'inrease' should be 'increase'.","section":"Discussion"},{"comment":"'red close symbols' should be 'red closed symbols'.","section":"Fig. 5 caption"},{"comment":"The axis labels appear garbled as 'N#(cm&')' and similar; please ensure correct typesetting of the emitter-density and fluence units.","section":"Fig. 2 axes"},{"comment":"The NV calibration is described only by reference [25]; please add a short description of the calibration procedure and its uncertainty in the Methods section.","section":"Results, NV calibration"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's central claim is overstated relative to the evidence. The data are valuable, and the paper could be publishable after a major revision that reframes Eq. (1) as a phenomenological fit with irradiation-dependent kappa, provides uncertainty estimates, and clearly delineates the regime of validity. I recommend major revision rather than rejection because the experimental dataset is broad and the annealing results are a useful contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know about this paper is that the measurement campaign is the real contribution: a broad, careful pulsed-ODMR comparison of T1 and T2 for silicon vacancies in 4H-SiC across neutron, electron, and proton irradiation over a wide fluence range, plus annealing data. The qualitative findings are solid: T1 stays near 300 µs up to NV ≈ 7×10^15 cm^-3 regardless of irradiation type, T2 degrades with damage, electron energy barely matters, and annealing partially recovers neutron-irradiated coherence. If you work on SiC defect spins, this is a useful reference dataset and a fair design-rule summary for choosing irradiation conditions.\n\nWhere the paper gets soft is the central 'universal scaling' claim. Equation (1), T2 = T2_pristine / (1 + T2_pristine κ NV), is presented as describing all data with two independent parameters, but κ is fitted separately for electrons, the two neutron reactors, and the literature proton points (0.8, 3.8, 25.2, and presumably larger values, all in units of 10^-11 s^-1 cm^3). The authors concede κ should depend on ion type, energy, and depth. With a different κ per condition, Eq. (1) is a parametrization, not a predictive law. The stress-test note is right about the load-bearing assumption: the decoherence rate is assumed proportional to PL-derived NV, but NV is only the optically active VSi density, not the total irradiation-induced spin bath. At high fluence the PL drops because non-radiative defects dominate, and the authors themselves see the fit deviate for NV > 7×10^15 cm^-3.\n\nOther soft spots, all minor-to-moderate: the neutron T2 values are self-admitted lower bounds due to microwave-field inhomogeneity; fit parameters come without uncertainties; many Fig. 5 points lack error bars; the NV calibration is described only by reference; and the data are not deposited. None of this kills the paper, but it does mean the 'universal scaling' headline overclaims. The authors are mostly honest about the limitations in the Discussion, so this is not a case of hidden circularity.\n\nWho is this for? SiC quantum-device experimentalists and anyone choosing between electron, neutron, or proton irradiation for making VSi ensembles. The paper deserves a serious referee: the dataset is new and the qualitative conclusions are reproducible in spirit. I would send it to review, with the request that the authors either soften 'universal' to 'phenomenological' with per-irradiation κ values, or provide an independent validation set, fit uncertainties, and tabulated data.\n\nRecommendation: engage with it; accept for review; expect heavy revision on the scaling narrative.","headline":"A genuinely useful systematic dataset for SiC emitter coherence, dressed up with a 'universal scaling' that is really a per-irradiation fit formula.","tokens_in":16334,"tokens_out":1585,"would_cite":true,"duration_ms":19247,"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 establishes that irradiation-induced spin decoherence in 4H-SiC follows one universal scaling with emitter density, $T_2 = T_2^{\\mathrm{pristine}}/(1+T_2^{\\mathrm{pristine}}\\kappa N_V)$, across electron, neutron, and proton…","keywords":["silicon carbide","silicon vacancy","spin coherence","irradiation-induced defects","optically detected magnetic resonance","4H-SiC","quantum sensing","spin-lattice relaxation"],"falsifier":"A decisive test is to measure $T_2$ in electron-irradiated samples with the same emitter density $N_V\\approx 10^{15}\\,\\mathrm{cm}^{-3}$ produced at 1 MeV and at 25 MeV: the formula with a single $\\kappa$ predicts identical $T_2$, so any resolvable difference would falsify the energy-independence claim. A second target: for one irradiation type, vary fluence across several $N_V$ values and check that $T_2^{\\mathrm{pristine}}/T_2-1$ is exactly linear in $N_V$ with zero intercept; curvature would show that the irradiation decoherence rate is not simply proportional to $N_V$.","tokens_in":15113,"feed_emoji":"🧲","tokens_out":12554,"duration_ms":109688,"temperature":0.7,"pith_summary":"The paper sets out to turn irradiation damage in silicon carbide from an uncontrolled nuisance into a predictable design parameter for quantum devices. Using pulsed optically detected magnetic resonance on silicon vacancies in 4H-SiC, the authors argue that the spin coherence time $T_2$ is governed by one universal law: $T_2 = T_2^{\\mathrm{pristine}}/(1+T_2^{\\mathrm{pristine}}\\kappa N_V)$, where $N_V$ is the density of optically active emitters and $\\kappa$ is a fixed decoherence rate per emitter for a given irradiation type. They find that electron energy from 1 to 25 MeV does not matter, that the spin-lattice relaxation time $T_1$ remains near 300 $\\mu$s up to $N_V\\approx 7\\times 10^{15}\\,\\mathrm{cm}^{-3}$ for all irradiation types, and that annealing can partially recover $T_2$ in heavily neutron-damaged samples. This matters because it gives a predictive rule for choosing irradiation particle, fluence, and annealing to hit a target coherence time in SiC-based qubits and sensors.","feed_headline":"One formula predicts spin coherence in irradiated silicon carbide","feed_subtitle":"Neutron, electron, and proton damage all shrink T2 along the same curve, making qubit tuning predictable.","key_machinery":"The load-bearing object is the negatively charged silicon vacancy in 4H-SiC (the V2 center, an $S=3/2$ defect spin read out by pulsed optically detected magnetic resonance, ODMR), and the load-bearing identity is the two-parameter formula $T_2 = T_2^{\\mathrm{pristine}}/(1+T_2^{\\mathrm{pristine}}\\kappa N_V)$. The mechanism behind it is additive spin-bath decoherence: separate channels from nuclear spins, paramagnetic impurities, residual defects, and irradiation-induced defects each contribute a rate, and the irradiation channel is taken to be proportional to the optically active emitter density $N_V$. The formula does the work of collapsing data taken with neutrons, electrons, and protons of several energies onto one family of curves, with $\\kappa$ parameterizing how much collateral damage each irradiation type produces per useful vacancy.","core_discovery":"On the paper's own terms, the discovery is a common scaling curve for irradiation-induced decoherence of the V2 silicon-vacancy spin in 4H-SiC. The decoherence rate is written as $1/T_2 = 1/T_2^{(n)}+1/T_2^{(p)}+1/T_2^{(r)}+1/T_2^{(i)}$, with the nuclear contribution $1/T_2^{(n)}$ suppressed by a 15 mT field, and the irradiation contribution assumed to grow linearly with the density of optically active vacancies: $1/T_2^{(i)}=\\kappa N_V$. The resulting expression $T_2 = T_2^{\\mathrm{pristine}}/(1+T_2^{\\mathrm{pristine}}\\kappa N_V)$ describes all measured electron, neutron, and proton data with two independent parameters; electron irradiation gives $\\kappa=0.8\\times 10^{-11}\\,\\mathrm{s^{-1}}\\,\\mathrm{cm^3}$, neutron irradiation gives larger $\\kappa$ values ($3.8$ and $25.2\\times 10^{-11}\\,\\mathrm{s^{-1}}\\,\\mathrm{cm^3}$ in the two reactors used), and proton data approach the same law. $T_1$ stays essentially constant at about 300 $\\mu$s up to $N_V=7\\times 10^{15}\\,\\mathrm{cm}^{-3}$, showing that spin-lattice relaxation is an intrinsic spin-phonon (Raman) property rather than an irradiation effect. Annealing a heavily neutron-irradiated sample raises $T_2$ by a factor of about 2.5 without removing the VSi centers, which the authors attribute to healing of other paramagnetic defects.","pith_inferences":["If the proportionality is exact, $\\kappa$ becomes a direct measure of collateral damage per useful emitter, so comparing $\\kappa$ across irradiation types should correlate with computed displacement damage and could be predicted from stopping-power simulations.","The same collapse plot could be applied to divacancy and other spin defects in SiC: if their decoherence is also dominated by irradiation-induced paramagnetic defects, their $T_2$ should fall on the same functional form with a defect-specific $\\kappa$.","The annealing results suggest a separate route to high coherence: irradiate densely and then selectively heal non-vacancy defects below the vacancy-annealing temperature; testing intermediate temperatures between $500$ and $700\\,^{\\circ}\\mathrm{C}$ would map the trade-off between recovered $T_2$ and lost $N_V$."],"forward_implications":["Once $\\kappa$ for a given beam line or reactor is known, Eq. (1) predicts $T_2$ at any fluence, so irradiation recipes can be designed to hit a target coherence time without trial-and-error.","Electron irradiation is the gentlest route to dense coherent ensembles, while neutron irradiation buys higher volume density at the price of shorter $T_2$; annealing recovers part of that loss.","$T_1$'s insensitivity up to $7\\times 10^{15}\\,\\mathrm{cm}^{-3}$ means the spin-lattice limit will not bottleneck quantum operations in that density range, leaving $T_2$ as the controlling figure of merit.","The same two-parameter curve should organize published and future data for other irradiation particles and energies once $N_V$ and $T_2^{\\mathrm{pristine}}$ are fixed, making the result a benchmark for device-oriented studies."],"supporting_citations":[{"why":"Supplies the nuclear-spin-bath decoherence model and the field-suppression effect that put $1/T_2^{(n)}$ inside $T_2^{\\mathrm{pristine}}$.","marker":"[6]"},{"why":"Gives the room-temperature spin-echo $T_2=81\\ \\mu$s and ESEEM analysis used as a comparison point for electron-irradiated HPSI wafers.","marker":"[7]"},{"why":"Provides the single-defect $T_1\\approx 500\\ \\mu$s and $T_2\\approx 100\\ \\mu$s and the coherent-control protocol that the pulsed ODMR measurements build on.","marker":"[27]"},{"why":"Provides the sample with known emitter density used to calibrate photoluminescence intensity to $N_V$.","marker":"[25]"},{"why":"Supplies the low-fluence proton-irradiation data point ($N_V\\approx 3\\times 10^{13}\\,\\mathrm{cm}^{-3}$, $T_2=42\\ \\mu$s) on the universal curve.","marker":"[47]"},{"why":"Supplies the higher-density proton-irradiation data point ($N_V\\approx 3\\times 10^{14}\\,\\mathrm{cm}^{-3}$, $T_2=22\\ \\mu$s) used for comparison with Eq. (1).","marker":"[50]"},{"why":"Defines how neutron fluence is counted (only epithermal and fast neutrons) and describes displacement damage in neutron-irradiated SiC.","marker":"[54]"},{"why":"Establishes the Raman spin-phonon origin of $T_1$ and the pulsed-ODMR methodology that underpin both measured times.","marker":"[9]"},{"why":"Provides the ensemble-coherence framework for decoherence by paramagnetic impurities, supporting the additive-rate form of Eq. (1).","marker":"[56]"}],"fun_headline_variants":["Universal curve links irradiation dose to spin coherence in SiC","Irradiation damage scales spin decoherence in SiC qubits","One equation predicts spin coherence for any irradiation type","T2 scales with vacancy density across all irradiation sources","SiC spin coherence: irradiation effect reduced to one curve"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The assumption the whole scaling rests on is that each additional optically active silicon vacancy adds exactly the same decoherence rate $1/T_2^{(i)}=\\kappa N_V$, regardless of irradiation type, particle energy, or local damage structure; if that proportionality breaks, Eq. (1) is only a curve fit and the claimed universality collapses.","fun_headline_variants_meta":{"raw":{"variants":["Universal curve links irradiation dose to spin coherence in SiC","Irradiation damage scales spin decoherence in SiC qubits","One equation predicts spin coherence for any irradiation type","T2 scales with vacancy density across all irradiation sources","SiC spin coherence: irradiation effect reduced to one curve"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000216,"raw_usage":{"total_tokens":1503,"prompt_tokens":1085,"completion_tokens":418,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":701,"completion_tokens_details":{"reasoning_tokens":338}},"tokens_in":701,"tokens_out":418,"duration_ms":4427,"temperature":1.0,"reasoning_tokens":338,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:32:53.799781+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test is to measure $T_2$ in electron-irradiated samples with the same emitter density $N_V\\approx 10^{15}\\,\\mathrm{cm}^{-3}$ produced at 1 MeV and at 25 MeV: the formula with a single $\\kappa$ predicts identical $T_2$, so any resolvable difference would falsify the energy-independence claim. A second target: for one irradiation type, vary fluence across several $N_V$ values and check that $T_2^{\\mathrm{pristine}}/T_2-1$ is exactly linear in $N_V$ with zero intercept; curvature would show that the irradiation decoherence rate is not simply proportional to $N_V$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the nuclear-spin-bath decoherence model and the field-suppression effect that put $1/T_2^{(n)}$ inside $T_2^{\\mathrm{pristine}}$."},{"cited_title":"Widmann, S.-Y","cited_arxiv_id":null,"evidence_quote":"Provides the single-defect $T_1\\approx 500\\ \\mu$s and $T_2\\approx 100\\ \\mu$s and the coherent-control protocol that the pulsed ODMR measurements build on."},{"cited_title":"Fuchs, B","cited_arxiv_id":null,"evidence_quote":"Provides the sample with known emitter density used to calibrate photoluminescence intensity to $N_V$."},{"cited_title":"Kraus, D","cited_arxiv_id":null,"evidence_quote":"Supplies the low-fluence proton-irradiation data point ($N_V\\approx 3\\times 10^{13}\\,\\mathrm{cm}^{-3}$, $T_2=42\\ \\mu$s) on the universal curve."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the higher-density proton-irradiation data point ($N_V\\approx 3\\times 10^{14}\\,\\mathrm{cm}^{-3}$, $T_2=22\\ \\mu$s) used for comparison with Eq. (1)."},{"cited_title":"Wendler, T","cited_arxiv_id":null,"evidence_quote":"Defines how neutron fluence is counted (only epithermal and fast neutrons) and describes displacement damage in neutron-irradiated SiC."},{"cited_title":"Simin, H","cited_arxiv_id":null,"evidence_quote":"Establishes the Raman spin-phonon origin of $T_1$ and the pulsed-ODMR methodology that underpin both measured times."}],"review_version":1}