{"id":"6e5a323a-4f92-4b75-974d-a4fb7aa83377","arxiv_id":"2412.02005","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Individual silicon-vacancy centers can be characterized by extracting ground and excited state strain/Jahn-Teller coupling rates from a photoluminescence excitation splitting and a coherent population trapping resonance.","lead":"This paper reports an experimental way to measure the energy-level structure of a single silicon-vacancy color center in diamond under a magnetic field. Two strain/Jahn-Teller coupling rates are extracted from two optical measurements, a step toward making these centers more predictable for quantum information uses.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The extracted coupling rates depend on assumed literature lambda_so values; if an individual center's lambda_so differs, all four reported rates shift, and the post-hoc unequal-f_e fit does not remove this sensitivity.","rationale":"The paper's central method is original and the transverse-field extraction is internally consistent: the analytic two-level structure in Eq. (6) follows from the model Hamiltonian, and the use of CPT plus PLE splitting gives an overdetermined, physically motivated route to the two coupling rates. I do not find a fatal flaw in the derivation or a circular argument in the extraction itself. The most load-bearing concern is the input lambda_so values: all four reported rates scale with these assumed parameters, and the model checks that use the same parameters cannot independently validate them. The reader's weakest_assumption identifies the same premise. I agree partially because the B-field dependence in Figs. 2 and 3 provides some support for the model, but that support is conditional on the very lambda_so values in question. The post-hoc unequal-f_e fit is a weaker verification and should not be presented as an independent confirmation. The proposed zero-field A-D spectrum test would settle whether the literature lambda_so values are adequate for individual centers. Given the current evidence, the conditional verdict remains appropriate: the method is promising and likely correct, but the quantitative numbers need an independent lambda_so check or explicit uncertainty propagation before they can be taken as definitive.","tokens_in":10498,"tokens_out":30504,"duration_ms":298711,"concrete_test":"For each of SiV1 and SiV2, record the zero-field PLE spectrum of the four A-D transitions. Under the model, the ground and excited doublet separations are sqrt(lambda_g^2 + Gamma_g^2) and sqrt(lambda_e^2 + Gamma_e^2). Using the paper's extracted Gamma_g and Gamma_e, solve for lambda_g and lambda_e from these zero-field line separations and compare with 45 and 257 GHz. If the inferred lambda values differ by more than the fitting uncertainty, the reported Gamma_s values are not robust; if they match, the assumed literature lambda_so values are independently supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 4, the measured CPT splitting and the C2-C3 splitting are converted into Gamma_s^g and Gamma_s^e using lambda_so^g = 45 GHz and lambda_so^e = 257 GHz from Ref. [29]. These literature values are not measured for the two implanted centers. The conversion uses Eq. (6), and for the ground state Gamma_s^g/lambda_so^g is about 0.7, so the exact relation between the measured splitting and Gamma_s is nonlinear and as sensitive to lambda_so as it is to Gamma_s. An error in lambda_so therefore shifts the extracted rates by a comparable fraction. The verification curves in Figs. 2 and 3 use the same literature lambda_so values as fixed inputs, so agreement with those curves tests the Hamiltonian structure only conditionally on the assumed lambda_so; the subsequent introduction of unequal quenching factors f_e = 0.13 and 0.17 is a post-hoc adjustment that provides additional flexibility but does not independently pin down lambda_so. If the large strain the paper infers modifies the effective spin-orbit splitting for an individual defect, the central quantitative results would change.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a combined model of static strain, dynamic Jahn-Teller coupling, and magnetic field for individual silicon-vacancy (SiV) centers, introducing an effective strain/JT coupling rate for the ground and excited doublets. Using an analytical diagonalization valid for a transverse magnetic field, the authors extract these rates for two SiV centers from two measurements: the frequency separation between the spin-conserved C2 and C3 transitions in a photoluminescence excitation spectrum, and the coherent-population-trapping resonance position of the ground-state doublet. They then compare calculated and measured magnetic-field amplitude and direction dependences, and introduce unequal ground/excited orbital quenching factors (0.1 vs 0.13 or 0.17) to improve agreement. The central claim is that this provides a general, two-measurement characterization route for individual SiV centers.","tokens_in":10770,"tokens_out":2941,"duration_ms":33844,"significance":"If the extraction is robust, the paper offers a practically useful method for characterizing the energy-level structure of individual SiV centers, a key requirement for optical spin control, quantum spin-mechanics, and cavity QED applications. The analytical diagonalization in Eq. (6) is a useful contribution, and the experimental data cover field-amplitude and field-direction dependences for two centers. The paper is explicit that the two extracted coupling rates are obtained from two measured frequencies with no additional fitting parameters, and it makes a falsifiable prediction of the field dependences shown in Figs. 2 and 3. However, the significance is tempered by the small number of centers, the absence of reported uncertainties, and the sensitivity of the quantitative extraction to literature spin-orbit parameters that are not measured for the individual centers.","major_comments":[{"comment":"The extracted values Γ_s^g = 33, 46 GHz and Γ_s^e = 77, 129 GHz are computed using λ_so^g = 45 GHz and λ_so^e = 257 GHz taken from Ref. [29] for a different sample. Because Γ_s^g is comparable to λ_so^g (ratios of 0.73 and 1.02), the relation between the measured splittings and Γ_s is not in the perturbative limit, and an error in λ_so shifts the extracted Γ_s by a comparable fraction. This sensitivity is load-bearing because the same assumed λ_so values are then used in all the verification curves in Figs. 2, 3, and 5. The manuscript needs a sensitivity analysis showing how the extracted rates and the agreement in Figs. 2 and 3 change when λ_so^g and λ_so^e vary within their plausible ranges, or an experimental constraint on λ_so for the measured centers.","section":"Sec. 4, Eq. (6)"},{"comment":"The unequal orbital quenching factors f_e = 0.13 and 0.17 are introduced after observing systematic deviations, and the manuscript reports no uncertainty or independent corroboration for these values. Since the abstract and conclusion claim that the work 'reveals contributions from unequal orbital magnetic coupling,' the post-hoc adjustment needs a quantitative justification: for example, a fit with confidence intervals, a chi-square comparison with the equal-f model, or a statement of the range of f_e values that are consistent with the data. Otherwise the claim of unequal quenching factors is weaker than the central coupling-rate extraction.","section":"Sec. 4, Figs. 2 and 3"},{"comment":"No experimental uncertainties are reported for the CPT resonance positions, the spin-conserved splittings, or the extracted Γ_s values. The scatter of data points around the calculated lines in Figs. 2 and 3 is not quantified, making it difficult to assess the claimed 'overall good agreements' and the systematic deviations that motivate f_e. Please add error bars (fit uncertainties, linewidths, magnetic-field calibration) and, if possible, report the fit parameters and their covariance for the Lorentzian and CPT fits.","section":"Secs. 3 and 4, Figs. 2-4"}],"minor_comments":[{"comment":"The equation numbering is inconsistent: there are two equations labeled (2), one for the magnetic-field Hamiltonian and one for the strain Hamiltonian. Please renumber the equations throughout the manuscript.","section":"Sec. 2, Eqs. (2) and (3)"},{"comment":"The notation 'X-Y plane' and 'z-Z plane' is confusing because x,y,z are used for the SiV internal axes and X,Y,Z for the lab frame. Please define the lab frame and SiV frame explicitly and use distinct symbols (e.g., lowercase for SiV axes and uppercase for lab axes) in all figures and text.","section":"Sec. 3, Fig. 2 caption and text"},{"comment":"The caption says 'dotted lines' and 'solid lines' but does not identify which curve corresponds to which C1-C4 transition. Please add a legend or explicit curve labeling so the figure can be read without reference to the main text.","section":"Sec. 4, Fig. 5 caption"},{"comment":"The statement that 'there is no other fitting parameter used in these calculations' is slightly overstated, since f_e is adjusted later in the same section. Please clarify the sequence: the spin-conserved splitting calculations use the extracted Γ_s values as fixed inputs, and f_e is a subsequent adjustment.","section":"Sec. 4, paragraph on quenching factors"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely problem and the proposed two-measurement extraction is elegant. My main concern is quantitative robustness: the extracted coupling rates are directly proportional to assumed λ_so values that are not measured for the two implanted centers, and the ratios Γ_s/λ_so are not small, so the central numbers could shift substantially. The post-hoc f_e adjustment also needs a more rigorous statistical treatment. These issues are fixable within the manuscript's scope and do not, in my view, invalidate the approach. I would support publication after the authors provide a sensitivity analysis and uncertainty estimates."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper has a genuinely useful idea—compress static strain and Jahn-Teller coupling into two effective rates and extract both from two standard measurements under a transverse field. That is new relative to prior work, which handled strain and JT separately or ignored one of them. The analytic eigenenergies in Eq. (6) are simple enough to make the extraction practical, and the field-amplitude and field-direction dependences in Figs. 2 and 3 do track the model. The two centers studied give concrete numbers (ground rates 33 and 46 GHz; excited rates 77 and 129 GHz) that the SiV community will want to compare against.\n\nThe soft spots are real but not disqualifying. The conversion from measured splittings to rates assumes lambda_so^g=45 GHz and lambda_so^e=257 GHz from Ref. [29], and those values are not measured for these two centers. The relation is nonlinear and as sensitive to lambda_so as to the measured splittings, so the absolute rates carry an unquantified systematic error. The paper does say 'estimated,' but it never bounds the sensitivity. Also, the unequal excited-state quenching factors (0.13 and 0.17) are introduced after seeing systematic deviations in the equal-quenching fits. That is honest—the prose says so—but it's a post-hoc adjustment, not a prediction, and with only two centers the generalization is unclear. On top of that, no uncertainties or raw data are reported, which makes it hard to assess fit quality.\n\nStill, the central method holds up. The extraction protocol is the contribution, and it works as a protocol. A careful referee should ask for an uncertainty analysis and a discussion of how much lambda_so can vary from center to center, plus ideally a couple more centers. That is a reasonable bar, not a rejection.\n\nI would bring this to a reading group and I would send it to review. It is the kind of paper that gives the SiV field a practical tool, even if the numbers will shift with better-known parameters later.","headline":"A practical two-measurement protocol for SiV strain/JT coupling, with unquantified systematic error from assumed spin-orbit splittings.","tokens_in":11284,"tokens_out":2570,"would_cite":true,"duration_ms":189304,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two optical measurements, a CPT resonance and a two-line PLE splitting, directly yield the ground- and excited-state strain/Jahn-Teller coupling rates of an individual diamond silicon-vacancy center.","keywords":["silicon vacancy center","diamond defect","Jahn-Teller coupling","static strain","energy level structure","coherent population trapping","photoluminescence excitation","transverse magnetic field"],"falsifier":"Measure the same two centers with an independent probe of the ground doublet splitting, for example by resolving the weak spin-flip C4 transition directly in a high-signal PLE spectrum and comparing its position with the CPT resonance, and check that the implied $\\tilde{\\epsilon}_g$ matches the value extracted by the paper's two-measurement route. A direct zero-field measurement of the ground- and excited-state doublet splittings on those same centers would also test the assumed $\\lambda_{\\rm so}$ values.","tokens_in":10311,"feed_emoji":"💎","tokens_out":8828,"duration_ms":80470,"temperature":0.7,"pith_summary":"The paper develops a way to characterize the energy-level structure of an individual silicon-vacancy (SiV) center in diamond without large fitting sets. It shows that in a magnetic field normal to the SiV axis, two measured quantities, the coherent population trapping (CPT) resonance between the two lower ground states and the frequency separation between the two spin-conserved C transitions in a photoluminescence excitation (PLE) spectrum, fix two effective coupling rates, one for the ground doublet and one for the excited doublet. These rates merge static strain and dynamic Jahn-Teller coupling, which share the same symmetry and cannot be separated in these experiments. For two measured centers the authors obtain ground-state rates of 33 and 46 GHz and excited-state rates of 77 and 129 GHz, and they use the field-amplitude and field-direction dependence to uncover unequal orbital quenching factors for ground and excited states. If the approach is right, a quick two-laser characterization can turn the large, seemingly random center-to-center variation in SiV spin splittings into a concrete diagnostic of local strain.","feed_headline":"Two measurements yield a diamond qubit's key level parameters","feed_subtitle":"A CPT resonance and one PLE line splitting give the ground- and excited-state strain/JT rates of a single SiV center.","key_machinery":"The load-bearing object is the effective strain/Jahn-Teller coupling rate $\\tilde{\\epsilon}=\\sqrt{(\\tilde{\\epsilon}_x)^2+(\\tilde{\\epsilon}_y)^2}$, defined separately for the ground and excited doublets, with $\\tilde{\\epsilon}_x$ and $\\tilde{\\epsilon}_y$ formed by adding the static strain components $\\beta,\\gamma$ to the Jahn-Teller coupling rates $\\delta_x,\\delta_y$ of the same symmetry. Because strain and Jahn-Teller terms share symmetry, only the combined rate enters the energy levels, and in a transverse magnetic field the Hamiltonian is analytically diagonalizable. The eigenenergy formula (Eq. 6) then shows that the splitting of each doublet is set by the product of the magnetic coupling $2\\gamma_s B_\\perp$ and $\\tilde{\\epsilon}$, divided by the spin-orbit splitting $\\lambda_{\\rm so}$; this direct proportionality is what lets one measurement fix $\\tilde{\\epsilon}_g$ and one additional measurement fix $\\tilde{\\epsilon}_e$. The same machinery predicts the full field-angle and field-amplitude dependence, which the authors use to verify the model and to expose the unequal orbital quenching factors.","core_discovery":"The paper's central claim is that the combined effect of static strain and dynamic Jahn-Teller coupling on a SiV center can be condensed into two scalar coupling rates, $\\tilde{\\epsilon}_g$ for the ground doublet and $\\tilde{\\epsilon}_e$ for the excited doublet, and that these rates are directly readable from the spectrum. For a transverse magnetic field the authors diagonalize the full Hamiltonian analytically; the resulting eigenenergies (Eq. 6) make the field-induced splitting of each doublet approximately $\\Delta E\\approx 2\\gamma_s B_\\perp\\tilde{\\epsilon}/\\lambda_{\\rm so}$. The CPT resonance gives $\\Delta E_g$, and the separation between the two spin-conserved C transitions gives $\\Delta E_e$, so both rates follow without numerical fitting. Applying this to two implanted centers yields $\\tilde{\\epsilon}_g = 33, 46$ GHz and $\\tilde{\\epsilon}_e = 77, 129$ GHz. The same model, when compared with the measured dependence on field amplitude and direction, also shows that the ground and excited doublets have unequal orbital quenching factors (0.1 for ground; 0.13 and 0.17 for excited), resolving systematic discrepancies that an equal-factor model leaves behind.","pith_inferences":["One natural extension is to use the same two-measurement protocol as a strain survey: scanning a sample center by center would map the local strain distribution that causes the large variations in SiV spin splittings.","The finding of unequal quenching factors suggests that the quenching factor may not be a universal constant for all SiV centers; a multi-center study could test whether $f_e/f_g$ varies systematically with implantation parameters or local strain.","Because Eq. 6 is analytical, real-time tracking of the spin-conserved splitting could in principle be used to monitor time-dependent strain fluctuations at a single center, turning the characterization method into a strain sensor."],"forward_implications":["For any individual SiV center, only two measurements are needed to determine the two strain/JT coupling rates: the CPT resonance position and the spin-conserved splitting of the C transition.","Because the model gives analytical eigenenergies for arbitrary transverse-field orientation, the measured angular dependence is predicted from the extracted rates with no additional fitting parameters.","The extracted rates for SiV1 and SiV2 (ground 33/46 GHz, excited 77/129 GHz) indicate strain contributions well above the earlier low-strain JT values of 15 and 76 GHz, showing that local strain, not just Jahn-Teller coupling, dominates these implanted centers.","Unequal ground and excited orbital quenching factors (0.1 vs 0.13/0.17) imply that treatments of SiV levels assuming a single quenching factor of 0.1 will deviate when the magnetic field is along or near the SiV axis."],"supporting_citations":[{"why":"It supplies the assumed spin-orbit splittings of 45 and 257 GHz and the earlier maximum Jahn-Teller coupling rates (15 and 76 GHz) that the paper's extracted rates are compared against.","marker":"[29]"},{"why":"It provides the symmetry-adapted strain Hamiltonian with components alpha, beta, gamma that the paper merges with the Jahn-Teller coupling to define the effective rates.","marker":"[6]"},{"why":"It establishes the coherent-population-trapping dark resonance used to measure the ground-state doublet splitting in the two-measurement protocol.","marker":"[26]"},{"why":"It describes the low-temperature confocal magneto-optical setup that enables the PLE and CPT measurements on individual centers.","marker":"[30]"}],"fun_headline_variants":["Two measurements decode SiV center's energy levels","CPT and PLE line split give SiV's strain-JT rates","Off-axis field exposes SiV orbital coupling asymmetry","SiV energy levels from one CPT and one PLE split","Unequal orbital quenching in SiV exposed by off-axis field"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The extraction's weakest point is that it takes the spin-orbit splittings $\\lambda_{\\rm so}^g=45$ GHz and $\\lambda_{\\rm so}^e=257$ GHz and the ground-state orbital quenching factor $f_g=0.1$ from earlier literature and applies them to each individual center; if a measured center's actual spin-orbit or quenching values differ, the quoted coupling rates shift accordingly.","fun_headline_variants_meta":{"raw":{"variants":["Two measurements decode SiV center's energy levels","CPT and PLE line split give SiV's strain-JT rates","Off-axis field exposes SiV orbital coupling asymmetry","SiV energy levels from one CPT and one PLE split","Unequal orbital quenching in SiV exposed by off-axis field"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00077,"raw_usage":{"total_tokens":3429,"prompt_tokens":979,"completion_tokens":2450,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":2367}},"tokens_in":595,"tokens_out":2450,"duration_ms":17845,"temperature":1.0,"reasoning_tokens":2367,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:56:08.367419+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the same two centers with an independent probe of the ground doublet splitting, for example by resolving the weak spin-flip C4 transition directly in a high-signal PLE spectrum and comparing its position with the CPT resonance, and check that the implied $\\tilde{\\epsilon}_g$ matches the value extracted by the paper's two-measurement route. A direct zero-field measurement of the ground- and excited-state doublet splittings on those same centers would also test the assumed $\\lambda_{\\rm so}$ values.","supporting_citations":[{"cited_title":"Maity, L","cited_arxiv_id":null,"evidence_quote":"It supplies the assumed spin-orbit splittings of 45 and 257 GHz and the earlier maximum Jahn-Teller coupling rates (15 and 76 GHz) that the paper's extracted rates are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the symmetry-adapted strain Hamiltonian with components alpha, beta, gamma that the paper merges with the Jahn-Teller coupling to define the effective rates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It establishes the coherent-population-trapping dark resonance used to measure the ground-state doublet splitting in the two-measurement protocol."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It describes the low-temperature confocal magneto-optical setup that enables the PLE and CPT measurements on individual centers."}],"review_version":1}