{"id":"e5baafc6-590e-449e-8bb2-45cf1399d209","arxiv_id":"2501.01133","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"First-principles calculations identify the neutral boron antisite and the positively charged nitrogen antisite as plausible origins of near-infrared and blue quantum emitters in hBN.","lead":"Using density-functional theory, the authors calculate the optical properties of antisite defects in hexagonal boron nitride and propose that two of them match observed single-photon emitters at 1.58 eV and in the blue range. If confirmed, this would give concrete microscopic models for quantum emitters in a leading 2D host material.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1.58 eV assignment relies on a static C3v/JT minimum, but the reported 2 meV HSE barrier makes the ground state a delocalized tunneling doublet, undermining the Franck-Condon ZPL.","rationale":"The reader's weakest_assumption already pointed to the 2 meV HSE barrier and the absence of convergence tests. I agree that the lack of sensitivity analysis is a genuine weakness, but the more load-bearing issue is physical: a 2 meV barrier between two equivalent out-of-plane configurations means the zero-point motion of the soft mode likely delocalizes the ground state, so the static C3v/JT minimum used for the Franck-Condon spectrum is not the correct starting point. This does not require distrusting the DFT methodology; it is an internal consistency problem between the reported barrier and the static-distortion picture. If this concern lands, the strongest claim about BN(0) as a 1.58 eV emitter is substantially weakened, although the NB(+) assignment is less affected because it does not rely on such a small barrier. The paper is honest and carefully hedged, and the computational framework is otherwise standard, so I would not reject it; the conditional verdict already given is appropriate. My concrete test would settle the question by checking the zero-point energy against the barrier and testing alpha sensitivity.","tokens_in":8031,"tokens_out":4611,"duration_ms":46316,"concrete_test":"Compute the HSE potential energy surface for BN(0) along the out-of-plane coordinate in the 256-atom supercell with a 3x3x1 k-point mesh and calculate the zero-point energy of the out-of-plane mode. If the zero-point energy exceeds the 2 meV barrier, the static Franck-Condon ZPL calculation is invalid; then recompute the spectrum from the delocalized vibronic ground state (for example, using a D3h-symmetric reference or a two-state tunneling model) and compare the resulting ZPL with 1.58 eV. As a sensitivity check, repeat with HSE alpha = 0.25 and 0.40; if the ZPL shifts by more than about 0.1 eV, the numerical assignment is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"BN(0) is assigned to the 1.58 eV emitter through a 0 K Franck-Condon calculation in which the symmetry-broken C3v/JT minimum gives ZPL = 1.58 eV and S = 1.68 (Table I). The calculation is only valid if the initial vibrational state is localized at that minimum. The paper reports a 2 meV HSE energy barrier between the two out-of-plane configurations ('The calculated energy difference is 59 meV with PBE functional while it is 2 meV with HSE functional'), and an imaginary 16 meV phonon at the planar geometry. A 2 meV barrier is below the expected zero-point energy of this soft out-of-plane mode, so the true 0 K ground state is a delocalized tunneling doublet (or the system is at most dynamically averaged), not a statically distorted C3v/C2v minimum. The Franck-Condon PL spectrum computed from the single JT-distorted minimum therefore does not represent the physical initial state, and the 1.58 eV ZPL and S = 1.68 are not well-defined. Since no convergence test or sensitivity analysis of the barrier or ZPL with respect to supercell size, k-point sampling, or HSE mixing parameter is reported, the central BN(0) assignment rests on a structural picture that may be internally inconsistent with the reported 2 meV barrier.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents DFT calculations (HSE with alpha=0.32, a 256-atom two-layer supercell, Gamma-only sampling) of native antisite defects in hexagonal boron nitride: B_N, N_B, and BNNB. Using Franck-Condon analysis, the authors compute zero-phonon-line energies, Huang-Rhys factors, and PL line shapes. They propose that the neutral boron antisite B_N(0) is a nonmagnetic single-photon emitter with a ZPL at 1.58 eV and a phonon sideband consistent with a 1.55 eV emitter observed in hBN thin films, and that the positively charged nitrogen antisite N_B(+) may explain a dim blue emitter with a ZPL at 2.63 eV. They also analyze the BNNB antisite pair and conclude it is likely not responsible for the observed blue spectra.","tokens_in":8261,"tokens_out":6328,"duration_ms":60933,"significance":"If correct, the identification of native antisite defects as quantum emitters would be a valuable step for hBN defect engineering. The paper uses standard, reproducible methods (HSE hybrid functional, Franck-Condon approximation) and provides quantitative predictions for ZPL energies, Huang-Rhys factors, formation energies, and radiative lifetimes. The N_B(+) assignment appears internally consistent. However, the central B_N(0) assignment rests on a static Jahn-Teller distorted minimum whose stability is questionable because the reported 2 meV HSE barrier is smaller than the estimated zero-point energy of the relevant out-of-plane mode. In addition, the lack of convergence tests for supercell size, k-point sampling, and the HSE mixing parameter weakens the quantitative claims. The paper's significance is therefore conditional on resolving these issues.","major_comments":[{"comment":"The reported 2 meV HSE energy barrier between the two out-of-plane configurations of B_N(0) is far smaller than the zero-point energy of the 16 meV imaginary mode (approximately 8 meV). The 0 K ground state should therefore be a delocalized tunneling doublet, not a static C3v-JT minimum. The Franck-Condon calculation that starts from the symmetry-broken minimum yields ZPL = 1.58 eV and S = 1.68 in Table I, but these quantities are not well-defined if the initial vibrational state is delocalized. The authors should either compute the vibronic spectrum explicitly including tunneling, or provide a concrete justification (for example, coupling to a bath that localizes the wavefunction) for why the static-minimum approximation is valid. As written, the central B_N(0) assignment is not supported by the calculations presented.","section":"Section on B_N(0) (paragraph beginning 'The BN is nonmagnetic...') and Table I"},{"comment":"The calculations use a single Gamma-point k-point sampling, a 256-atom supercell, and a fixed HSE mixing parameter alpha = 0.32, with no convergence tests reported. The ZPL assignments (1.58 eV and 2.63 eV) are quantitative, and the computed barrier for B_N(0) is highly functional-sensitive (59 meV with PBE versus 2 meV with HSE). The authors must demonstrate that the key results are robust with respect to supercell size, k-point sampling, and alpha. Without this, the agreement with experiment could be fortuitous, and the 2 meV barrier in particular cannot be trusted as a basis for the proposed dynamic Jahn-Teller picture.","section":"Methods paragraph (supercell, k-point, HSE parameters)"}],"minor_comments":[{"comment":"The abstract says 'solid-state single phonon emitters'; this should be 'single photon emitters'.","section":"Abstract"},{"comment":"The phrase 'Experimental observed emission' should read 'Experimentally observed emission'.","section":"Abstract and introduction"},{"comment":"The rows labeled D3h in Table I are not explained in the text; clarify whether these are unrelaxed high-symmetry geometries or constrained solutions.","section":"Table I"},{"comment":"The sentence 'the high symmetry configuration of excited state is JT unstable due to the half occupation of the e state' leaves unspecified which high-symmetry configuration is meant; please specify the point group of the excited state before the JT distortion.","section":"Discussion of B_N(0) excited state"},{"comment":"The phrase 'the data of low energy range is missing' is vague; specify which energy range of the experimental spectrum is absent and how this affects the comparison.","section":"Comparison with experiment (Fig. 2b)"},{"comment":"Reference [51] is a preprint from Researchsquare; if a peer-reviewed version is now available, it should be cited instead.","section":"Reference [51]"},{"comment":"Equation (2) uses the bulk refractive index nD = 2.1 for a two-layer hBN system; the authors should note that this is an approximation and discuss its possible effect on the radiative lifetime.","section":"Equation (2)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is generally well-written and the computational protocol is standard for the field. The main concern is the internal consistency of the B_N(0) assignment: a 2 meV barrier is far below the zero-point energy of the relevant mode, so the Franck-Condon calculation from a static distorted minimum is not justified. The lack of any convergence or sensitivity analysis compounds this problem. I would request a revision that addresses the tunneling issue, adds convergence tests, and clarifies the comparison with experiment. The paper is not ready for acceptance in its present form, but the identified problems are potentially fixable within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is the first systematic computational look at native antisite defects in hBN. It calculates ZPL energies, Huang–Rhys factors, and simulated PL spectra for the boron and nitrogen antisites in several charge states, using standard ΔSCF and Franck–Condon methods. The two headline candidates are the neutral boron antisite for the 1.6 eV NIR emitters and the positive nitrogen antisite for the dim blue emitter. The authors hedge appropriately with “might” and “tentatively,” and the comparison to experiment is honest about the limits.\n\nThe work is genuinely new: previous defect assignments in hBN focused on vacancies, carbon, and oxygen, while antisites were mostly ignored. The inclusion of the out-of-plane Jahn–Teller distortion and its effect on the optical transition is a useful contribution. The methods are standard enough to be reproducible, and the paper is refreshingly clear about what is and isn’t established.\n\nThe soft spots are real, and one is load-bearing. The paper reports a 2 meV HSE barrier between the two out-of-plane configurations of BN(0), along with an imaginary 16 meV phonon at the planar geometry. The zero-point energy of that soft mode is several meV, so the true 0 K ground state should be a delocalized tunneling doublet, not a static C3v minimum. The Franck–Condon PL spectrum computed from the single JT minimum therefore does not represent the physical initial state, and the 1.58 eV ZPL and S = 1.68 are not well-defined. The authors do not address this, and they report no convergence tests or sensitivity analysis for the barrier, the ZPL, the supercell size, or the HSE mixing parameter. That is the central weakness.\n\nThe blue emitter assignment to NB(+) is more robust because the D3h symmetry is stable and there is no JT complication, but it is based on limited experimental correlation, and the 2.63 eV ZPL sits about 0.2 eV below the commonly cited 2.8–2.9 eV blue line.\n\nOverall, this is a useful survey and a reasonable set of candidates, but the main assignment is not yet established. The paper deserves peer review because the question is important and the approach is standard; the reviewer should ask the authors to treat the dynamic JT effect properly and to provide convergence checks. I would cite it cautiously and bring it to a reading group for discussion.\n\nRecommendation: send to review, conditional on addressing the tunneling issue and the numerical convergence tests.","headline":"First systematic optical study of native antisite defects in hBN, but the main 1.58 eV assignment rests on a 2 meV tunneling barrier that likely invalidates the static Franck–Condon picture.","tokens_in":8806,"tokens_out":2333,"would_cite":true,"duration_ms":22658,"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":"A density-functional study proposes that native antisite defects, not impurities, underlie two observed single-photon emitters in hexagonal boron nitride: the neutral boron antisite with a 1.58 eV zero-phonon line and the positively…","keywords":["hexagonal boron nitride","antisite defect","single-photon emitter","zero-phonon line","Jahn-Teller distortion","density functional theory","phonon sideband","color center"],"falsifier":"Look for the predicted roughly 16 meV out-of-plane mode as a low-energy replica in the photoluminescence sideband of the 1.58 eV emitter, and check its spin dependence: the line should show no Zeeman splitting or optically detected magnetic resonance because BN(0) is nonmagnetic; either signature would rule out the assignment.","tokens_in":7796,"feed_emoji":"💎","tokens_out":7809,"duration_ms":67321,"temperature":0.7,"pith_summary":"The paper tries to establish that two native antisite defects are the microscopic sources of specific quantum emitters observed in hexagonal boron nitride (hBN). It assigns the near-infrared single-photon line around 1.58 eV to the neutral boron antisite, a boron atom on a nitrogen site, and the dim blue line at 2.63 eV to the positively charged nitrogen antisite. The decisive mechanism is an out-of-plane Jahn-Teller distortion: the boron antisite is dynamically unstable in the plane, and this soft out-of-plane motion reshapes the optical transition through a Huang-Rhys factor of 1.68. If the assignments are right, two simple intrinsic defects can explain emitters that were previously attributed to unclear vacancy or impurity centers, giving experimenters concrete atomic models to confirm or refute.","feed_headline":"Two native defects may explain hBN quantum emitters","feed_subtitle":"Assigns a 1.58 eV near-infrared line to the boron antisite and a 2.63 eV blue line to the nitrogen antisite.","key_machinery":"The central objects are the native antisite defects: BN, a boron atom replacing a nitrogen site, and NB, a nitrogen atom replacing a boron site, together with the antisite pair BNNB. The argument is carried by the dynamic Jahn-Teller mechanism: in BN(0) a roughly 16 meV imaginary phonon shows that the planar configuration is unstable, and the defect relaxes into one of two equivalent out-of-plane C3v configurations with a computed 2 meV barrier between them. Exciting the defect breaks the symmetry further to C2v, and the Franck-Condon overlap between the vibrational modes of the ground and excited states produces the zero-phonon line and the Huang-Rhys factor that are compared with experiment. The nitrogen antisite works through a pseudo-Jahn-Teller effect in which excitation restores the high-symmetry D3h geometry, which explains the absence of a permanent dipole moment.","core_discovery":"The central claim is that the neutral boron antisite BN(0) is a nonmagnetic, near-infrared single-photon source with a zero-phonon line at 1.58 eV and a phonon sideband that matches the frequently observed hBN spectra, while the positively charged nitrogen antisite NB(+) is a plausible dim blue emitter with a 2.63 eV zero-phonon line. Both defects are single substitutions, yet their optical behavior is controlled by geometry: BN(0) has two equivalent out-of-plane configurations separated by only 2 meV in the screened-hybrid calculation, and the Jahn-Teller symmetry breaking lowers the zero-phonon line from 1.74 eV to 1.58 eV and raises the Huang-Rhys factor from 1.09 to 1.68. NB(+) keeps high D3h symmetry, has no permanent dipole, and yields a 2.63 eV zero-phonon line with Huang-Rhys factor 2.13. The paper also simulates the antisite pair BNNB but concludes that its excessively strong phonon sideband rules it out as the source of the observed blue emitters.","pith_inferences":["Because the 2 meV barrier is so small, a practical check would be to recompute the BN(0) zero-phonon line with larger supercells and denser k-point sampling; if the barrier and the C3v relaxation change, the 1.58 eV value is not stable.","A direct experimental test would be resonant photoluminescence excitation on a single 1.58 eV emitter: the sideband spacing should match the predicted out-of-plane mode frequency, and the lack of a permanent dipole could be checked by Stark spectroscopy.","If out-of-plane distortion is indeed generic in hBN defects, strain or isotope engineering should tune the emission through the soft mode, offering a way to control single-photon energies beyond the electronic level structure.","The nonmagnetic prediction extends to the ground and excited states: magneto-optical measurements at the single-emitter level could separate BN(0) from the many spin-active candidates among hBN defects."],"forward_implications":["The 1.58 eV near-infrared line in hBN should be treated as a candidate for the neutral boron antisite, with its characteristic out-of-plane phonon replica at roughly 16 meV.","The dim blue emitter at 2.63 eV should be tested against the nitrogen antisite's signatures: D3h symmetry, no permanent dipole, and a Huang-Rhys factor near 2.13.","Out-of-plane relaxation must be included when modeling other planar defects in hBN, because the soft mode, not just the electronic gap, determines the emission lineshape.","The antisite pair BNNB, despite a bright transition with a 3 Debye dipole and an 11.7 ns radiative lifetime, is unlikely to be the observed blue emitter because its simulated sideband is far stronger than the measured one.","Under nitrogen-rich growth conditions, the Fermi level pinning and formation energies put both candidate antisites at low concentration, which is consistent with their appearing as rare single-photon sources."],"supporting_citations":[{"why":"Supplies the experimental 1.55 eV photoluminescence spectrum in hBN thin film whose low-energy phonon replica is matched to the simulated BN(0) spectrum.","marker":"[31]"},{"why":"Reports nonmagnetic type-II emitters near 1.6 eV, providing the observed family of near-infrared emitters that BN(0) is tentatively assigned to.","marker":"[43]"},{"why":"Reports a blue emitter with nonlinear Stark shift and very small transition dipole, the key experimental constraint used to argue for the high-symmetry NB(+) defect.","marker":"[47]"},{"why":"Provides the recently observed 2.6 eV photoluminescence peak that is robust against annealing, supporting the 2.63 eV blue emitter candidate.","marker":"[48]"},{"why":"Supplies one of the experimental blue-emitter spectra near 435 nm used for comparison with the BNNB simulation.","marker":"[13]"},{"why":"Establishes the Delta SCF approach used to compute excited-state total energies and zero-phonon-line energies.","marker":"[38]"},{"why":"Provides the Franck-Condon method used to simulate photoluminescence lineshapes and to extract Huang-Rhys factors.","marker":"[39]"},{"why":"Confirms experimentally that certain hBN quantum emitters exhibit out-of-plane distortion, supporting the geometric mechanism central to this paper.","marker":"[42]"},{"why":"Supplies the pseudo-Jahn-Teller framework used to describe the NB(0) excitation that restores the high-symmetry geometry.","marker":"[44]"}],"fun_headline_variants":["Boron antisite: a nonmagnetic single-photon source","Two native defects explain hBN quantum emitters","hBN lines at 1.58 and 2.63 eV from antisites","Single substitutions drive hBN emission colors","Out-of-plane phonons control hBN emitter spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the screened-hybrid density-functional calculations, especially the tiny 2 meV energy barrier between the two out-of-plane configurations of the neutral boron antisite, predict the zero-phonon-line energies and Huang-Rhys factors accurately enough to assign emitters.","fun_headline_variants_meta":{"raw":{"variants":["Boron antisite: a nonmagnetic single-photon source","Two native defects explain hBN quantum emitters","hBN lines at 1.58 and 2.63 eV from antisites","Single substitutions drive hBN emission colors","Out-of-plane phonons control hBN emitter spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000168,"raw_usage":{"total_tokens":1263,"prompt_tokens":952,"completion_tokens":311,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":568,"completion_tokens_details":{"reasoning_tokens":229}},"tokens_in":568,"tokens_out":311,"duration_ms":3607,"temperature":1.0,"reasoning_tokens":229,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:34:08.443143+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Look for the predicted roughly 16 meV out-of-plane mode as a low-energy replica in the photoluminescence sideband of the 1.58 eV emitter, and check its spin dependence: the line should show no Zeeman splitting or optically detected magnetic resonance because BN(0) is nonmagnetic; either signature would rule out the assignment.","supporting_citations":[{"cited_title":"Huang, M","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental 1.55 eV photoluminescence spectrum in hBN thin film whose low-energy phonon replica is matched to the simulated BN(0) spectrum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports nonmagnetic type-II emitters near 1.6 eV, providing the observed family of near-infrared emitters that BN(0) is tentatively assigned to."},{"cited_title":"Zhigulin, J","cited_arxiv_id":null,"evidence_quote":"Reports a blue emitter with nonlinear Stark shift and very small transition dipole, the key experimental constraint used to argue for the high-symmetry NB(+) defect."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the recently observed 2.6 eV photoluminescence peak that is robust against annealing, supporting the 2.63 eV blue emitter candidate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies one of the experimental blue-emitter spectra near 435 nm used for comparison with the BNNB simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the Delta SCF approach used to compute excited-state total energies and zero-phonon-line energies."},{"cited_title":"Hoese, P","cited_arxiv_id":null,"evidence_quote":"Confirms experimentally that certain hBN quantum emitters exhibit out-of-plane distortion, supporting the geometric mechanism central to this paper."},{"cited_title":"Li, J.-P","cited_arxiv_id":null,"evidence_quote":"Supplies the pseudo-Jahn-Teller framework used to describe the NB(0) excitation that restores the high-symmetry geometry."}],"review_version":1}