{"id":"7a1535a1-369d-4b72-9733-83c06390df76","arxiv_id":"2607.03322","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"ADAQ high-throughput DFT of >8000 c-BN defects, refined by HSE, reassigns the GC-2 line to ONV_B^- and nominates C_N V_B^-, C_B V_N^+ and Na_B^- as bright spin-1/2 emitters.","lead":"High-throughput DFT screening of >8000 point defects in cubic boron nitride identifies several bright color-center candidates, including a re-attribution of the experimental GC-2 line to the negative oxygen-vacancy complex. The work supplies a public defect database and concrete experimental tests for quantum emitters and sensors in an underexplored wide-bandgap host.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified to the central GC-2 reassignment; residual uncertainties are standard for the field and already flagged by the reader.","rationale":"The reader correctly isolates the GC-2 reassignment as the strongest claim and the one-phonon Debye-Waller filter as the weakest methodological assumption. That assumption affects only the ranking of secondary candidates (C_NV_B^-, Na_B^-, etc.) and is already acknowledged as a typical field limitation. Because the central spectroscopic argument for ONV_B^- does not rely on the one-phonon approximation and is supported by direct HSE total-energy differences that match experiment once α is adjusted to the experimental gap, no load-bearing concern arises that would alter the CONDITIONAL verdict. The concrete test above simply tightens the already-solid GC-2 identification without changing the overall assessment.","tokens_in":17612,"tokens_out":515,"duration_ms":5019,"concrete_test":"Recompute the ONV_B^- ZPL with HSE (α=33 %) in a 1000-atom supercell (or with a 2×2×2 k-mesh) and extract the full luminescence lineshape including multi-phonon contributions; if the ZPL remains within 50 meV of 1.63 eV and the Debye-Waller factor stays >5 %, the GC-2 assignment is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s strongest claim—that ONV_B^- (HSE α=33 % ZPL = 1.600 eV) is a more plausible origin of the experimental GC-2 line at 1.63 eV than the previously proposed neutral ONV_B—rests on transparent hybrid-functional total-energy differences and a clear Jahn-Teller analysis (Sec. III A, Table IV). The low-symmetry JT state of the neutral defect lies ~0.3 eV below the observed energy and requires a two-phonon process to populate the high-symmetry configuration, while the negative charge state matches both energy and Debye-Waller factor without such caveats. No internal inconsistency appears in the electronic-structure diagrams (Figs. 4–5), the charge-transition levels, or the comparison with prior work. The one-phonon ΔQ ranking used for the broader candidate list (reader’s weakest assumption) is secondary to this spectroscopic reassignment and does not undermine it.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript presents a high-throughput DFT screening (ADAQ/PBE) of >8000 point-defect complexes (s/p-element substitutions, vacancies, interstitials) in c-BN, filtering for thermodynamic stability on the defect hull, ZPL > 0.5 eV, high TDM, low mass-weighted geometry change ΔQ (Eq. 1), and non-zero spin. A short-list is refined with HSE06 (α = 25 % and 33 %), focusing on optical properties, Jahn-Teller distortions and symmetry. The central spectroscopic claim is that the negatively charged oxygen-vacancy complex ONV_B^- (HSE α = 33 % ZPL = 1.600 eV) is a more plausible microscopic origin of the experimental GC-2 line at 1.63 eV than the previously proposed neutral ONV_B, whose low-symmetry JT state lies ~0.3 eV lower and requires multi-phonon population of the high-symmetry configuration. Additional bright candidates (C_NV_B^-, C_BV_N^+, Na_B^-) are identified and characterized.","tokens_in":17821,"tokens_out":1321,"duration_ms":19651,"significance":"c-BN remains underexplored relative to diamond and h-BN for quantum defects; a systematic combinatorial screen that produces an open database and several concrete, experimentally testable candidates is therefore valuable. The reassignment of GC-2 rests on transparent hybrid-functional total-energy differences, explicit JT analysis (Figs. 4-5, Table IV) and convergence checks (supercell size, k-points, α), which are strengths. The work also supplies falsifiable experimental proposals (oxygen-concentration dependence, ODMR of D/E, doping-dependent PL) and releases the full defect data set, both of which raise the impact beyond a pure computational catalogue.","major_comments":[{"comment":"Sec. II C and Eq. (1): the one-phonon Huang-Rhys/Debye-Waller ranking constructed solely from mass-weighted ΔQ is used as a hard filter (ΔQ < 1 amu^{1/2} Å) that reduces 384 ZPL-qualified defects to the 104 candidates of Table II and ultimately to the HSE short-list of Table III. Multi-phonon side-band structure and non-radiative rates are not computed; for several short-listed defects (especially Ca_NV_B and C_BV_N^+) the one-phonon DW factors are already 10^{-9}-10^{-13}, so the ranking that justified their inclusion (or exclusion of others) is not robust. A short sensitivity discussion or a note that the GC-2 reassignment itself does not rely on this filter would strengthen the claim.","section":"Sec. II C, Eq. (1), Fig. 1, Tables II-III"},{"comment":"Sec. III A / Fig. 2 / Table IV: the assertion that ONV_B^- is the more plausible GC-2 centre is based on HSE ZPL matching, yet all formation energies and charge-transition levels (including the ε(+/0) = 1.42 eV quoted for the hull) are reported only at the PBE level. Because the PBE gap is 4.45 eV versus the experimental 6.4 eV, the Fermi-level window of stability for the negative charge state can shift by several hundred meV under HSE. Reporting at least the HSE charge-transition levels (or a brief statement that the negative state remains on the hull for mid-gap Fermi energies) is needed to close this gap in the central claim.","section":"Sec. III A, Fig. 2, Table IV"}],"minor_comments":[{"comment":"Throughout: inconsistent spacing and notation for defect labels (ONVB / O NVB / O_NVB / ONV_B). Adopt a single, standard subscript style (e.g., O_NV_B) and apply it uniformly, including in the abstract and Table I.","section":"Abstract, Tables I-III, Sec. III"},{"comment":"Fig. 1 caption and text: 'eciency' contains a non-ASCII character; replace with 'efficiency'.","section":"Fig. 1"},{"comment":"Table III footnote a and Sec. III C: the statement that the HSE result for C_NV_B^- is 'the third lowest excitation' is clear only after reading the full paragraph; a one-sentence clarification in the table caption would help.","section":"Table III, Sec. III C"},{"comment":"Appendix A: the VASP version dependence (regular 5.4.4 vs bug-fixed) that causes the excited state of C_BVB to relax back to the ground state should be flagged more prominently, as it affects reproducibility of the 1.009 eV ZPL.","section":"Appendix A"},{"comment":"Sec. II B: the plane-wave cut-off (600 eV PBE, 520 eV HSE) and stopping criteria are given, but a brief note on residual forces or total-energy convergence for the 512-atom cells would reassure readers that the ΔQ values are numerically stable.","section":"Sec. II B"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid, well-executed high-throughput study whose central GC-2 reassignment is internally consistent and already more carefully documented than many comparable defect papers. The two major points above are fixable with modest additional calculations or clarifying text; I do not see them as grounds for major revision or rejection. Fit for a materials/condensed-matter journal is good."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real news is the combinatorial scale and the GC-2 reassignment. They ran ADAQ on >8000 s/p-element defects in a 512-atom cell, filtered on hull stability, ZPL >0.5 eV, TDM and ΔQ, then did careful HSE (α=25 % and 33 %) on a short list. The public database is useful; the new HSE numbers for ONV_B^-, C_N V_B^-, C_B V_N^+ and Na_B^- are not in the literature.\n\nWhat they do well: the JT analysis on neutral ONV_B is transparent. Low-symmetry excited state sits ~0.3 eV below the high-symmetry one (E_JT ≈ 0.21 eV), so the old 1.60 eV match to GC-2 is unlikely. Negative charge with α=33 % lands at 1.600 eV and needs no two-phonon rescue. They also give concrete experimental tests (oxygen concentration series, ODMR D/E, doping-dependent PL). Carbon and sodium candidates look bright enough to be worth a look, and the C_B V_B appendix reproduces the earlier special-k-point result once they use the same settings.\n\nSoft spots are ordinary for this genre. Debye-Waller factors rest on the one-phonon ΔQ approximation; multi-phonon and non-radiative rates could reorder the lower-ranked candidates, but that does not touch the GC-2 claim. Charge-transition levels are still PBE-level, and they keep the PBE lattice constant for the HSE runs. None of this is hidden.\n\nThis is for people who actually grow or measure c-BN defects, or who maintain defect databases. It is not a methods paper and not a theory-of-everything. The math and citation pattern look solid; self-cites are to the workflow they are using. I would send it to referees without hesitation and would cite the GC-2 reassignment and the Na_B^- / carbon numbers myself. Worth bringing to the next reading group if anyone is working on BN or high-throughput defect screening.","headline":"Solid high-throughput map of c-BN color centers that cleanly reassigns GC-2 to ONV_B^- and flags a few new bright spin-1/2 candidates; residual one-phonon ranking is a standard field limitation, not a load-bearing flaw.","tokens_in":18508,"tokens_out":577,"would_cite":true,"duration_ms":5505,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"High-throughput screening of over 8000 defects in cubic boron nitride identifies the negatively charged oxygen-vacancy center as the likely source of the GC-2 emission line and flags carbon and sodium defects as additional bright emitters.","keywords":["cubic boron nitride","color centers","point defects","zero-phonon line","hybrid DFT","quantum emitters","oxygen-vacancy complex","high-throughput screening"],"falsifier":"Measure photoluminescence spectra and zero-field-splitting parameters on c-BN samples whose oxygen concentration and Fermi level are deliberately varied; if the intensity of the 1.63 eV GC-2 line does not scale with oxygen content or the measured D and E parameters disagree with the calculated values for ONV_B^-, the assignment is ruled out.","tokens_in":18449,"feed_emoji":"💎","tokens_out":1016,"duration_ms":8255,"temperature":0.7,"pith_summary":"Cubic boron nitride is a wide-bandgap crystal whose optical defects could serve as single-photon sources and quantum sensors, yet most measured emission lines remain unassigned. This paper systematically generates and screens more than eight thousand point-defect complexes containing s- and p-block impurities with density-functional theory, then re-examines the most promising candidates with hybrid functionals. The central result is that the negatively charged oxygen-vacancy complex better matches the experimental GC-2 zero-phonon line at 1.63 eV than the neutral charge state previously proposed, while two carbon-related defects and a sodium substitutional also emerge as bright, spin-bearing emitters in technologically useful wavelength windows. A sympathetic reader cares because an identified microscopic origin immediately suggests doping and spectroscopy experiments that can confirm or refute the assignment and because the same database supplies a short-list of new color-center candidates for quantum technologies.","feed_headline":"Oxygen-vacancy defect matches the GC-2 line in cubic BN","feed_subtitle":"High-throughput screen of 8000 defects also flags carbon and sodium centers as bright emitters","key_machinery":"The Automatic Defect Analysis and Qualification (ADAQ) high-throughput workflow that generates, relaxes and ranks thousands of substitutional, vacancy and interstitial complexes by formation energy on the defect hull, zero-phonon line, transition dipole moment and mass-weighted geometry difference ΔQ, followed by hybrid-functional refinement of the short-listed candidates.","core_discovery":"After filtering a database of more than 8000 calculated defects for thermodynamic stability, zero-phonon lines above 0.5 eV, high transition dipole moments and low geometry change, hybrid-functional calculations show that the negatively charged oxygen-vacancy complex ONV_B^- yields a zero-phonon line of 1.600 eV (exact-exchange mixing 33 percent) and is therefore a more plausible microscopic origin of the experimental GC-2 line than the previously suggested neutral ONV_B. The same calculations identify the carbon defects C_NV_B^- and C_BV_N^+ and the sodium defect Na_B^- as additional bright emitters.","pith_inferences":["If the Jahn-Teller stabilization of the neutral oxygen-vacancy excited state is as large as calculated, earlier literature that matched only the high-symmetry geometry to GC-2 will need systematic re-examination.","Controllable p- and n-doping of c-BN, still experimentally difficult, becomes a high-priority materials goal because it would allow charge-state control of the newly proposed emitters.","The same screening pipeline can be reapplied to other wide-gap hosts once their experimental zero-phonon lines are catalogued."],"forward_implications":["The GC-2 line can be assigned to ONV_B^- once oxygen-concentration and doping-dependent photoluminescence and ODMR measurements are performed.","Carbon-rich c-BN is predicted to host two spin-1/2 emitters near 1.03 eV and 1.74 eV that are accessible with existing growth methods.","Sodium doping is predicted to produce a bright near-infrared spin-1/2 emitter at 1.30 eV.","The public defect database supplies a ready short-list of further candidates for hybrid-functional or experimental follow-up."],"fun_headline_variants":["Oxygen-vacancy complex better matches GC-2 line in c-BN","Hybrid DFT links ONVB- to experimental GC-2 emission","Screen of 8000 defects flags ONVB- for GC-2 origin","Carbon and sodium defects join as bright c-BN emitters","Negatively charged oxygen-vacancy explains GC-2 line"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The ranking of which defects emit brightly enough to be useful rests on a simple one-phonon estimate of the Debye-Waller factor from geometry change; multi-phonon and non-radiative processes are not calculated and could reorder the candidates.","fun_headline_variants_meta":{"raw":{"variants":["Oxygen-vacancy complex better matches GC-2 line in c-BN","Hybrid DFT links ONVB- to experimental GC-2 emission","Screen of 8000 defects flags ONVB- for GC-2 origin","Carbon and sodium defects join as bright c-BN emitters","Negatively charged oxygen-vacancy explains GC-2 line"]},"model":"grok-4.5","effort":"low","cost_usd":0.007626,"raw_usage":{"total_tokens":1881,"prompt_tokens":819,"num_sources_used":0,"completion_tokens":93,"cost_in_usd_ticks":76260000,"prompt_tokens_details":{"text_tokens":819,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":969,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":819,"tokens_out":93,"duration_ms":6898,"temperature":1.0,"reasoning_tokens":969,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-12T03:16:40.591925+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Measure photoluminescence spectra and zero-field-splitting parameters on c-BN samples whose oxygen concentration and Fermi level are deliberately varied; if the intensity of the 1.63 eV GC-2 line does not scale with oxygen content or the measured D and E parameters disagree with the calculated values for ONV_B^-, the assignment is ruled out.","supporting_citations":[],"review_version":1}