{"id":"3332720d-901f-40f0-ad28-a28a7bde1405","arxiv_id":"2501.19165","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"NV-center T1 noise magnetometry should show a sharp relaxation enhancement when its frequency matches the energy gap between a vacancy-bound Majorana zero mode and a low-energy hybridized Majorana mode in a non-Abelian Kitaev spin liquid.","lead":"The paper predicts that a nitrogen-vacancy center in diamond can detect the exotic Majorana fermions bound to vacancies in a Kitaev quantum spin liquid, by measuring how its spin relaxation time changes with an applied magnetic field. If correct, this gives experimentalists a new, non-invasive probe for the long-sought non-Abelian phase in Kitaev materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted 1/T1 resonance positions depend on flux-gap denominators Δ, Δ1, Δ2, Δ3 in the effective Hamiltonian (Supp. Eqs. S23–S38), which are never computed or stated; without them the central quantitative claim is not reproducible.","rationale":"The reader's weakest assumption is the correct one to stress. The paper is a theoretical proposal whose novelty is the specific resonant 1/T1 signature; the quantitative position and height of that resonance are the central claim. The analytic Majorana mapping and the projection-operator algebra in Supplementary Section III are standard and, as far as the text shows, internally consistent. The weak point is not the formalism but the input parameters to the effective Hamiltonian. The authors explicitly distinguish Δ1, Δ2, Δ3 in the perturbative coefficients, yet give no way to obtain them, and no code or parameter file is provided. Since the resonances in Fig. 2(b) are sharp, even 20–30% variations in these denominators can move peaks by more than the announced linewidth and could change which field window is identified as the non-Abelian signature. This is an addressable, not disqualifying, problem: computing the denominators is a finite exact-diagonalization calculation in the Kitaev sector, so the conditional verdict stands and could be upgraded to accept if the check is done and the results are robust.","tokens_in":23792,"tokens_out":5030,"duration_ms":54645,"concrete_test":"Compute Δ, Δ1, Δ2, Δ3 for the exact vacancy configuration of Fig. 2 (bound-flux pair on the 40×40 lattice, vacancies at (32,27) A-site and (9,14) B-site) by diagonalizing the non-Zeeman Kitaev Majorana Hamiltonian in the relevant Z2 flux sectors, using the same method that gives the 0.065J bulk flux gap. Insert the computed values into Eqs. (S23)–(S26) and (S35)–(S38), regenerate the low-energy spectrum and the 1/T1(EZ) curve of Fig. 2, and compare the resonance field and peak amplitude with the published result. If the peak shifts by more than ~20% in EZ or disappears, the quantitative central claim needs revision; if it moves negligibly, the unspecified denominators are not load-bearing. A complementary check is exact diagonalization of the full spin model on a small cluster with the same two vacancies to test the effective Hamiltonian's spectrum directly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative prediction — a specific field-tunable 1/T1 peak whose position is set by the energy difference between a vacancy-bound MZM and a low-energy hybridized mode — is computed from the effective Hamiltonian of Supplementary Section II. Every term that determines the hybridized mode energies, Eqs. (S23)–(S26) and (S35)–(S38), is proportional to inverse powers of local flux-gap denominators Δ, Δ1, Δ2, Δ3. These denominators are never computed or stated anywhere in the paper or supplement. The only related number given is the bulk flux gap, ~0.065J (Ref. 67), which applies to flipping a gauge field far from vacancies. The processes in Fig. S1(b)–(d) flip different sets of bonds around a vacancy, so Δ1, Δ2, Δ3 should differ from the bulk gap and from each other; the paper even keeps them distinct, but then does not provide values. Because the resonance condition is an eigenvalue of the effective hopping matrix T, an O(1) uncertainty in a denominator translates directly into an O(1) uncertainty in the predicted resonance field and in the peak shape of 1/T1(EZ). This does not invalidate the mechanism, but it means the headline quantitative signature, as presented, is neither reproducible nor independently checkable without the missing input parameters or a code release.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes using T1-based noise magnetometry with an NV center qubit to detect Majorana zero modes bound to vacancies in a non-Abelian Kitaev quantum spin liquid. In the Majorana representation, a Zeeman field releases the dangling b-Majorana fermions adjacent to a vacancy, which hybridize with the itinerant c-Majorana sector and form low-energy modes; the authors show that the NV relaxation rate 1/T1 is resonantly enhanced when the NV working frequency Omega = Omega_0 + (g_q/g_s) EZ matches the energy difference between such a hybridized mode and a vacancy-bound MZM. The paper derives the 1/T1 formula and the required spin correlation functions, constructs a third-order perturbative effective Majorana Hamiltonian for the vacancy sector (Supplementary Section II), and computes the field-dependent 1/T1 spectrum and its spatial distribution for one vacancy-pair configuration on a 40x40 lattice, with additional data for NV polarization, temperature, and field direction in the supplement. The authors estimate signal increments of order 30 Hz at d=30a and 3 kHz at d=10a at roughly 1 T and 1 K, and discuss the experimental challenges.","tokens_in":24104,"tokens_out":22155,"duration_ms":201833,"significance":"If the mechanism proposed here works in practice, NV noise magnetometry would be a genuinely new, non-invasive, local probe of vacancy-bound Majorana zero modes in Kitaev materials, complementary to the thermal Hall and tunneling techniques emphasized in the introduction. The paper makes falsifiable qualitative predictions that are already useful: the resonance is field-tunable (Fig. 2), and the maximal signal is displaced from the vacancy cores (Fig. 3), a feature the authors verify by comparing the approximation in Eq. (6) with the full expression in Eq. (1). The derivation of the projected spin correlation function in Supplementary Section III (Eqs. (S51)-(S77)) is non-trivial, and the 1/T1 peaks in Fig. 2(b) emerge from the spectrum of the effective Hamiltonian rather than being fitted to the target result, so circularity is not a concern. The main caveat is quantitative: the flux-gap denominators that set the scale of the effective Hamiltonian are never stated, and the perturbation expansion is used near the edge of its validity, so the predicted resonance fields must be regarded as provisional until those inputs are provided or benchmarked numerically.","major_comments":[{"comment":"The effective hopping matrix T that determines the hybridized-mode energies in Fig. 2(a), and therefore the resonance field in Fig. 2(b), is built from the couplings J~', J~'', kappa', kappa'', t^(1), and t^(2) in Eqs. (S23)-(S26) and (S35)-(S38), each of which contains inverse powers of the flux-gap denominators Delta_1, Delta_2, Delta_3 (and of the unspecified Delta_in in Eq. (S36)). These local flux-gap energies are never computed or stated in the paper or the supplement; the only related number given is the bulk flux gap of about 0.065J (main text, Ref. [67]), which applies to flux flips far from vacancies and may differ from the vacancy-adjacent values that the paper itself keeps distinct. Because the predicted 1/T1 peak position is an eigenvalue of the effective matrix T, the central quantitative claim of the paper is not reproducible or independently checkable without these inputs. I request that the authors report the computed values of Delta_1, Delta_2, Delta_3 (and Delta_in) for the bound-flux background, state the conversion between the perturbative field h and the plotted EZ, and show how the resonance fields in Fig. 2(b) shift under O(1) variations of these denominators.","section":"Supplementary Section II, Eqs. (S23)-(S38)"},{"comment":"The effective Hamiltonian in Supplement Section II is quoted as accurate to order h^3, but at the resonance fields shown in Fig. 2(a) the expansion parameter is not small: the ratio of the t^(1) to t^(0) couplings is of order h/Delta, which lies between roughly 0.3 and 0.6 at EZ ~ 0.02-0.03J with Delta ~ 0.065J. The hybridized-mode energies are therefore subject to possibly large corrections from the omitted orders, so the quantitative accuracy of the predicted peak positions is not established. I recommend benchmarking the effective Hamiltonian against exact diagonalization or a higher-order calculation for the same vacancy configuration on a small cluster; if that is not feasible, the quantitative claims in the Results section and in Supplementary Section V should be reframed as order-of-magnitude estimates. This concern does not affect the qualitative mechanism, namely a low-energy mode whose energy grows with the Zeeman field and crosses the NV working frequency.","section":"Supplementary Section II.D; Fig. 2(a)"},{"comment":"All numerical results (Figs. 2, 3, S3-S6) are obtained for one vacancy configuration on one supercell: vacancies at the A-sublattice of unit cell (32,27) and the B-sublattice of unit cell (9,14) on a 40x40 lattice with the Zeeman field along the a-axis. The abstract and the summary claim a general pathway for identifying non-Abelian KQSLs, but no robustness scan over vacancy separation, sublattice arrangement, or supercell size is provided. This matters for the experimental estimates, because the spatial footprint (over one-third of the supercell) is a property of the wavefunction geometry of this particular configuration. I would like to see at least one additional configuration (for example, two vacancies on the same sublattice, or a significantly different separation) to confirm that the field-tunable resonance and the size of the footprint are robust features rather than artifacts of the chosen geometry.","section":"Results and discussion; Figs. 2-3"}],"minor_comments":[{"comment":"The title contains the typo 'magnetormetry' for 'magnetometry', and the same misspelling pattern appears as 'Kiteav' and 'Kiatev' throughout the main text and supplement (for example, 'Kiteav honeycomb model' in the introduction, 'Kiteav lattice' in the Fig. 1 caption, and 'KITAEV'/'Kiatev' in Supplementary Section II).","section":"Title and abstract"},{"comment":"The resonance condition used to interpret Fig. 2 (the NV frequency being matched by the energy difference between a hybridized mode and an MZM, with Omega = Omega_0 + (g_q/g_s) EZ) is stated only in words; an explicit equation, together with the conversion between EZ and the perturbative field h of the supplement, would greatly improve reproducibility.","section":"Main text, Eq. (1) and Fig. 2"},{"comment":"In the Fig. 3(d) caption, 'are obtained rom Eq. (1)' should read 'are obtained from Eq. (1)'.","section":"Fig. 3(d) caption"},{"comment":"In Supplementary Section V, the quoted increments (1/T1 larger than 30 Hz at d=30a and 3 kHz at d=10a over one-third of the supercell) should specify the baseline against which the increment is measured, and a short sensitivity statement for the assumed g-factors (g_q = 2, g_s = 2.5) and the J range (2.8-6.8 meV) would help the reader assess the robustness of the estimates.","section":"Supplementary Section V, Eq. (S80)"},{"comment":"In Fig. S4, the three temperature curves are not labeled in the caption; please state which curve corresponds to k_B T = 0.01J, 0.02J, and 0.05J, since the claim that bulk-mode contributions remain small at roughly 3 K relies on distinguishing them.","section":"Fig. S4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and its novelty is adequate: noise magnetometry of vacancy-bound Majorana modes is distinct from the authors' earlier PRL (Takahashi et al., PRL 131, 236701 (2023)), which uses non-local spin correlations, and from the tunneling probes cited in the introduction. My main concern is the one flagged in the report: the missing flux-gap denominators make the headline quantitative prediction unreproducible. The authors are well positioned to fix this, since the same numerical machinery that produced Fig. 2 can produce the flux-flip energy costs; a benchmark of the effective Hamiltonian against exact diagonalization would additionally resolve my concern about the expansion parameter. Should the benchmark show large deviations, the paper's claims could still be reframed as qualitative or order-of-magnitude predictions without losing its core message. I have no concerns about the overlap with prior work, which is properly cited."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: this paper has a genuinely new proposal — using NV center T1 noise magnetometry to see vacancy-bound Majorana zero modes in a non-Abelian Kitaev quantum spin liquid. The predicted field-tunable resonance between a MZM and a low-energy hybridized dangling-Majorana mode is a concrete, falsifiable signature. If it works, it is a real alternative to tunneling probes. The core derivation of 1/T1 from the noise spectrum and the Majorana correlation functions is careful, and the supplement does a lot of honest legwork with the projection operator and Wick's theorem. They also talk straight about experimental difficulty: NV placement, field alignment, low temperature.\n\nThe soft spot is exactly where the stress-test note points. The effective Hamiltonian in Supplement II uses three vacancy-specific flux-gap denominators Δ1, Δ2, Δ3 (plus the bulk Δ) in every perturbation term, Eqs. (S23)–(S38). These numbers are never computed or stated anywhere. The resonance condition is an eigenvalue of the effective hopping matrix, so an O(1) uncertainty in a denominator gives O(1) uncertainty in the predicted 1/T1 peak field and peak shape. That means the headline quantitative signature, as presented, is not reproducible or independently checkable. This doesn't kill the mechanism — the effect is robust to the exact values — but it does mean a reader cannot verify the numbers without either the missing inputs or a code release.\n\nMinor issues: only one vacancy configuration is shown, and there is no clean-system baseline in the same calculation (bulk modes are stated to be negligible, which is convincing but not quantified). The estimate of 1/T1 enhancement relies on a simplified far-field formula, but that part is clearly stated and reasonable.\n\nWho should read this: people working on Kitaev materials and quantum spin liquid detection; also anyone interested in theory-paper reproducibility, because the missing denominators are a textbook example of a checkable claim that isn't checkable. I'd send it to a serious referee; the idea deserves airing and the gaps are fixable. I'd want the referee to demand values for Δ1..Δ3 or a code release before acceptance.","headline":"New NV-magnetometry route to vacancy-bound Majoranas in Kitaev spin liquids, but the quantitative resonance positions hinge on unspecified flux-gap denominators that must be pinned down.","tokens_in":24594,"tokens_out":2385,"would_cite":true,"duration_ms":22652,"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":"This paper claims that the relaxation time of a nitrogen-vacancy (NV) center qubit above a vacancy-doped Kitaev spin liquid carries a sharp, field-tunable signature of Majorana zero modes bound to the vacancies.","keywords":["Kitaev quantum spin liquid","Majorana zero modes","NV center magnetometry","T1 relaxation","vacancy defects","non-Abelian anyons","noise spectroscopy","α-RuCl3"],"falsifier":"Compute the flux-gap energy denominators $\\Delta_1$, $\\Delta_2$, $\\Delta_3$ of the supplementary perturbation theory from the microscopic Kitaev parameters and check whether the resulting hybridized-mode spectrum reproduces the resonance fields of Fig. 2; if the true denominators shift the predicted $1/T_1$ peak to a Zeeman field where no enhanced relaxation is observed in a vacancy-doped sample, the central claim fails. Alternatively, a fixed-frequency noise-magnetometry scan over field that shows no $1/T_1$ peak anywhere in the perturbative window would falsify the prediction directly.","tokens_in":23638,"feed_emoji":"💎","tokens_out":12674,"duration_ms":102533,"temperature":0.7,"pith_summary":"The paper proposes that a nitrogen-vacancy (NV) center qubit placed above a vacancy-doped Kitaev spin liquid can act as a detector for the non-Abelian phase. The argument is that a Zeeman field releases dangling Majorana fermions at lattice vacancies, which hybridize into low-energy modes together with Majorana zero modes bound to the vacancies, and that the NV relaxation rate $1/T_1$ jumps whenever the NV working frequency matches the energy difference between a hybridized mode and a zero mode. Because an external field tunes that energy difference, the signature is a field-swept resonance rather than a fixed background feature. This matters because detecting isolated Majorana zero modes is regarded as a key step toward fault-tolerant topological quantum computation, and current experimental probes of Kitaev materials rely mostly on tunneling and thermal-transport measurements.","feed_headline":"Diamond qubit T1 spike flags Majorana modes in Kitaev liquids","feed_subtitle":"Sweeping the magnetic field should make the NV relaxation rate peak at each vacancy-bound Majorana resonance.","key_machinery":"The load-bearing object is a quadratic effective Majorana Hamiltonian (Supplementary Eq. S34), obtained by treating the Zeeman term as a perturbation to order $h^3$ in the Kitaev honeycomb model. Away from vacancies the field renormalizes the $c$-Majorana nearest-neighbor hopping and generates next-nearest-neighbor hopping; on the sites bordering a vacancy it couples the released dangling $b$-Majoranas to the itinerant $c$-Majoranas, which is what produces the MZMs and the low-energy hybridized modes. This Hamiltonian enters the $1/T_1$ formula of Eq. (1), a Fermi golden-rule expression in which the relaxation rate is a lattice sum of the spin correlation matrix evaluated at the NV frequency, weighted by a dipole-orientation matrix; the spatial approximation of Eq. (6) then reduces the signal to a $1/r^6$ falloff set by the principal components of the correlation and orientation matrices at each vacancy, which is why the maximum of $1/T_1$ sits away from the vacancy positions.","core_discovery":"On its own terms, the paper shows that the $T_1$-based noise spectrum of an NV center is a sensitive, field-tunable probe of Majorana zero modes (MZMs) trapped at vacancies in a non-Abelian Kitaev quantum spin liquid (KQSL). In a site-diluted Kitaev honeycomb model at the isotropic point, each vacancy leaves dangling $b$-Majorana fermions on its neighboring sites; a Zeeman field releases these from the local $Z_2$ gauge structure and hybridizes them with the itinerant $c$-Majorana fermions, producing MZMs bound to the vacancies together with low-energy hybridized modes whose wave functions overlap strongly with the vacancy-adjacent sites. The spin correlations of these modes enter the NV relaxation rate through a Fermi golden-rule expression, so $1/T_1$ increases sharply whenever the NV working frequency matches the energy difference between a hybridized mode and an MZM, and the external Zeeman field sweeps that difference through resonance. Because spin correlations away from the vacancies are gated by the large flux gap, the predicted $1/T_1$ enhancement stands out from other fluctuations in the spin liquid and should be observable in candidate materials such as $\\alpha$-RuCl$_3$.","pith_inferences":["The mechanism should generalize beyond pairs of vacancies: any defect that binds Ising anyons in a Kitaev spin liquid (a possibility the paper itself notes) should produce an analogous field-tunable $1/T_1$ resonance, so the same measurement could map which defects in a real material actually host zero modes.","The offset of the maximum $1/T_1$ away from the vacancy positions is a distinctive fingerprint: a signal from an ordinary local magnetic impurity would peak directly above the defect, so imaging the relaxation-rate pattern could distinguish fractionalized Majorana physics from impurity magnetism.","A practical test is an exfoliated Kitaev-material flake bonded to diamond with controlled vacancies: holding the NV frequency fixed and sweeping the in-plane field should give a $1/T_1$ peak at the predicted resonance field, and the absence of any such peak across the accessible field range would rule the scenario out."],"forward_implications":["NV noise magnetometry becomes a field-tunable diagnostic for the non-Abelian phase in Kitaev materials, supplementing tunneling and thermal-Hall measurements.","Holding the NV frequency fixed and sweeping the Zeeman field should produce a $1/T_1$ peak whenever the field brings the MZM-hybridized-mode energy difference through resonance.","Because spin correlations away from vacancies are gated by the flux gap, the predicted $T_1$ signature should survive at temperatures up to roughly 3 K and against other spin fluctuations in the spin liquid.","Rotating the in-plane Zeeman field away from the $a$-axis can lower the field needed for the first $1/T_1$ resonance to about $0.01J$, easing experimental requirements."],"supporting_citations":[{"why":"Supplies the exactly solvable honeycomb model, the Majorana representation of the spins, and the non-Abelian phase whose Majorana zero modes the paper aims to detect.","marker":"[1]"},{"why":"Shows that vacancies in the non-Abelian Kitaev spin liquid trap Ising anyons in the bound-flux sector, the starting state for the calculation.","marker":"[47]"},{"why":"Establishes the vacancy-induced low-energy density of states in the Kitaev spin liquid that underlies the site-dilution picture.","marker":"[54]"},{"why":"The NV-center noise-magnetometry framework, $T_1$ relaxometry, that the paper adapts to Kitaev spin liquids.","marker":"[56]"},{"why":"Supplies the $T_1$-based magnetic noise spectroscopy method for two-dimensional systems that the paper extends to the Kitaev model.","marker":"[57]"},{"why":"Contains the perturbative effective Hamiltonian to order $h^3$ and the spin-correlation derivation (Eq. 5) that produce the predicted $1/T_1$ peaks.","marker":"[65]"},{"why":"Gives the flux gap of about $0.065J$ below which spin correlations away from vacancies vanish, the basis for the claimed protection of the signal.","marker":"[67]"}],"fun_headline_variants":["NV T1 peak flags vacancy-bound Majorana modes in Kitaev liquid","Field-swept T1 resonance exposes non-Abelian Kitaev phase","Noise magnetometry reveals Majorana zero modes in Kitaev spin liquid","Zeeman-tuned T1 spikes probe Majoranas in Kitaev quantum spin liquid","Vacancy Majorana modes light up NV T1 spectrum in Kitaev liquids"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the perturbative effective Hamiltonian in Supplementary Section II is quantitatively accurate, because the flux-gap energy denominators $\\Delta$, $\\Delta_1$, $\\Delta_2$, $\\Delta_3$ appearing in Eqs. (S23)-(S38) are never computed or stated in the paper, yet they set the hybridized-mode energies and therefore the Zeeman-field values at which the predicted $1/T_1$ peaks occur.","fun_headline_variants_meta":{"raw":{"variants":["NV T1 peak flags vacancy-bound Majorana modes in Kitaev liquid","Field-swept T1 resonance exposes non-Abelian Kitaev phase","Noise magnetometry reveals Majorana zero modes in Kitaev spin liquid","Zeeman-tuned T1 spikes probe Majoranas in Kitaev quantum spin liquid","Vacancy Majorana modes light up NV T1 spectrum in Kitaev liquids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000302,"raw_usage":{"total_tokens":1758,"prompt_tokens":984,"completion_tokens":774,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":600,"completion_tokens_details":{"reasoning_tokens":674}},"tokens_in":600,"tokens_out":774,"duration_ms":7615,"temperature":1.0,"reasoning_tokens":674,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T21:04:31.897455+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the flux-gap energy denominators $\\Delta_1$, $\\Delta_2$, $\\Delta_3$ of the supplementary perturbation theory from the microscopic Kitaev parameters and check whether the resulting hybridized-mode spectrum reproduces the resonance fields of Fig. 2; if the true denominators shift the predicted $1/T_1$ peak to a Zeeman field where no enhanced relaxation is observed in a vacancy-doped sample, the central claim fails. Alternatively, a fixed-frequency noise-magnetometry scan over field that shows no $1/T_1$ peak anywhere in the perturbative window would falsify the prediction directly.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that vacancies in the non-Abelian Kitaev spin liquid trap Ising anyons in the bound-flux sector, the starting state for the calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the vacancy-induced low-energy density of states in the Kitaev spin liquid that underlies the site-dilution picture."},{"cited_title":"Agarwal, R","cited_arxiv_id":null,"evidence_quote":"Supplies the $T_1$-based magnetic noise spectroscopy method for two-dimensional systems that the paper extends to the Kitaev model."},{"cited_title":"Baskaran, S","cited_arxiv_id":null,"evidence_quote":"Gives the flux gap of about $0.065J$ below which spin correlations away from vacancies vanish, the basis for the claimed protection of the signal."}],"review_version":1}