{"id":"8ddf80ae-2d3e-48cd-8142-83e4dd3a6728","arxiv_id":"2502.05490","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"DFT calculations predict that the symmetry of hBN defects controls whether their emission energy shifts linearly or quadratically under an electric field, and that the dielectric environment strongly changes the shift magnitude.","lead":"This paper uses density functional theory to compute how electric fields shift the emission energy of seven candidate atomic defects in hexagonal boron nitride. It finds that the defect's local symmetry decides whether the shift is linear or quadratic in field strength, which could help identify which defects produce the observed single-photon emitters.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed CNCB3 match to experiment relies on ε=3.52, but the paper's own local-field slope fixes ε≈2.69; using the measured slope removes the agreement.","rationale":"I considered two candidate load-bearing concerns: the dielectric-constant convention (reader's choice) and the slab-thickness dependence of the polarizability. The latter is explicitly acknowledged only for CB, not for CNCB3, and is therefore not the decisive issue for the paper's main experimental match. The former is decisive: the paper's quantitative assignment of CNCB3 to the experimentally measured ~150 ų emitter hinges on multiplying the fitted polarizability by ε² with ε chosen as 3.52. However, the paper's own Fig. 4c/d computes the local-field screening factor as 0.37, which fixes ε=2.69. Adopting ε=3.52 is therefore inconsistent with the computed local field and is introduced post hoc to match experiment. This is not a minor error-bar issue; it removes the quantitative evidence for the CNCB3 assignment. The qualitative symmetry rule—D3h → quadratic, Cs with out-of-plane distortion → linear, C2v → quasi-quadratic—does not depend on the rescaling and is well grounded in group theory and the computed geometries. Hence the correct verdict remains CONDITIONAL: the symmetry classification is likely correct, but the specific defect identifications and numerical values require a consistent local-field calculation or direct experimental corroboration. The reader's CONDITIONAL verdict is therefore unchanged.","tokens_in":9959,"tokens_out":11998,"duration_ms":118794,"concrete_test":"Recompute Table I using the screening factor derived directly from the moving-averaged electrostatic potential of the actual 3-layer defect slab (not the pristine 5/9-layer slope), or equivalently re-fit the Stark shifts of CNCB3 with the local field set by the measured slope 0.37 (ε=2.69). If the resulting polarizability is not within the experimental uncertainty of 150 ų, the claimed match is an artifact of the ε=3.52 rescaling. A complementary check: repeat the CNCB3 Stark-shift calculation in a 5-layer or 9-layer slab and verify convergence of the extracted polarizability; if it does not converge near the experimental value, the 3-layer-based assignment is unreliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the identification of CNCB3 as the quadratic emitter with polarizability close to the experimental ~150 ų. This match is only obtained if one adopts ε=3.52 in Table I. But the paper's own local-field calculation (Fig. 4c,d) measures the slope of the moving-averaged electrostatic potential as 0.0185 eV/Å for an applied field of 0.05 eV/Å, i.e., a screening factor of 0.37, which fixes ε=2.69 via slope=ε⁻¹. The ε=3.52 value corresponds to a screening factor of 0.284, contradicting the computed local field; it is merely a thickness convention chosen to bring CNCB3's polarizability from 96.5 to 161.3 ų, bracketing the experimental 150 ų. Since the paper makes no independent determination of the correct local-field factor for the 3-layer defect slab, the claimed agreement is a post-hoc rescaling, not a prediction. If the measured slope is used, CNCB3's polarizability is 96.5 ų, a 36% underestimate, and the dipole moments for the linear defects are ~1.15 D instead of ~1.5 D. The symmetry-based classification itself (linear vs quadratic) is independent of this rescaling and remains plausible, but the specific assignment of CNCB3 to the ~2 eV emitter and the quoted numerical values are not supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles (HSE-DFT with ΔSCF) calculations of the Stark shifts of eight candidate defects in hexagonal boron nitride, using 3-layer slab models with an applied out-of-plane electric field. The central claim is that the local symmetry of the defect controls the shape of the Stark shift: centrosymmetric D3h defects (CB, CNCB3) exhibit quadratic shifts, Cs defects with out-of-plane distortion (CBON, BN, B_DB) exhibit linear shifts, and C2v defects exhibit quasi-quadratic shifts. The authors fit dipole-moment and polarizability changes from the computed energy shifts, discuss the dielectric screening that converts the applied slab field to the local field, and compare the resulting values with experimental Stark data on ~2 eV emitters, proposing CNCB3 as a candidate for the quadratic emitter.","tokens_in":10191,"tokens_out":6120,"duration_ms":54088,"significance":"The symmetry-based classification is the paper's strongest contribution: it is simple, consistent with the calculated relaxed geometries, and makes a falsifiable experimental prediction that the linear-versus-quadratic character of the Stark shift can be used to identify the symmetry class of individual hBN emitters. The paper also deserves credit for explicitly analyzing the sensitivity of the results to the thickness convention used to define the slab dielectric constant (Eq. 1, Fig. 4, Table I), rather than hiding this dependence. However, because the quantitative match to experiment for the leading candidate CNCB3 is obtained only by adopting a dielectric constant that contradicts the paper's own local-field slope, the numerical predictions and the specific defect assignment are not yet robust. The classification itself, being largely symmetry-based, survives this concern.","major_comments":[{"comment":"The claimed agreement between the CNCB3 polarizability and the experimental ~150 Å3 value is obtained only by adopting εh,⊥=3.52, but the local-field slope computed in Fig. 4(c,d) fixes εh,⊥≈2.69: the moving-averaged electrostatic potential gives 0.0185 eV/Å for an applied field of 0.05 eV/Å, a screening factor of 0.37. With the paper's own computed screening factor, the fitted CNCB3 polarizability is 96.5 Å3, a 36% underestimate, and the linear-defect dipole moments are ~1.15 D rather than ~1.5 D. The manuscript should therefore either (i) determine the local-field factor for the 3-layer defect slab from first principles, consistently with the 5- and 9-layer calculation, or (ii) present the CNCB3 match not as a prediction but as an illustration of the dielectric-convention sensitivity, and state clearly that the quantitative assignment is convention-dependent.","section":"Section III, Eq. (1), Fig. 4(c,d), Table I"},{"comment":"The statement that the calculated linear dipole moments are \"consistent with experiment value from -0.9 to 0.9 D\" is not supported by the numbers: with εh,⊥=2.69, CBON (1.15 D) and B_DB (-1.14 D) lie outside that range, and with εh,⊥=3.52 all three values (1.49, 1.87, -1.49 D) exceed it substantially. This should be reworded or supplemented with a justification for comparing to a broader experimental distribution, since this comparison is part of the linear-defect identification.","section":"Section III, Fig. 5b, Table I"},{"comment":"The classification of CNCB3 as a quadratic emitter and the fitted polarizability rely on fixing D3h symmetry by excluding the dynamic Jahn-Teller effect, as the text states (\"the dynamic JT effect is not included to fix the symmetry\"). This is a load-bearing assumption for the paper's main candidate assignment: if the JT effect is active, the degeneracy of the e'' state is lifted, inversion can be broken, and a linear Stark component would appear. The authors should quantify the JT stabilization energy (e.g., from a symmetry-broken calculation including electron-phonon coupling) or otherwise justify that the D3h approximation is valid at the relevant energies and timescales, and show that the quadratic Stark behavior is robust to this approximation.","section":"Section III, CNCB3 discussion and Fig. 2"}],"minor_comments":[{"comment":"The sentence \"The projector augmented wave (PAW) potentials ... is used\" contains a subject-verb agreement error; it should be \"are used.\"","section":"Section II"},{"comment":"The polarizability values are written as bare numbers with a trailing \"3\" (e.g., \"53.0 3\", \"95.6 3\"); the unit \"ų\" should be inserted consistently in the text and the table.","section":"Section III and Table I"},{"comment":"The origin of the adopted εh,⊥=3.52 is not explained: Fig. 4b states that the thickness-rescaled value can reach 3.25, and the text does not show how 3.52 follows from any specific thickness convention; please provide the thickness leading to 3.52 and relate it to the physical interlayer distance.","section":"Section III"},{"comment":"The statement that the CNCB3 ZPL of 2.04 eV is \"not far\" from 1.88 eV should be quantified, since the difference is 0.16 eV; please report the experimental line width or the range of reported ZPL values used for the comparison.","section":"Section III"},{"comment":"Several language issues appear, including \"especically\" in the Introduction, \"experimental date 0.24 D\" for \"experimental data\", \"external magnetic field is a effective to flip\" for \"is an effective way to flip\", and \"nevertheless its might be not the case\" for \"it might not be the case\"; these should be corrected.","section":"Section III"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about the dielectric-convention dependence, but the abstract and conclusion present the CNCB3 match as a quantitative identification, which currently rests on an ε=3.52 choice that the paper's own local-field calculation does not justify. A revision that either computes the local field in the 3-layer defect slab or clearly separates the symmetry-based prediction from the convention-dependent numerical values would make the paper suitable for publication. The symmetry classification seems worth preserving even if the quantitative assignment is softened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The durable part of this paper is the symmetry classification: D3h defects give quadratic Stark shifts, Cs defects with out-of-plane distortion give linear shifts, and C2v defects give quasi-quadratic shifts. That is standard defect physics, but the paper applies it systematically to seven hBN defect models and tabulates fitted dipole and polarizability changes. If that survives review, it gives experimentalists a quick fingerprint for unknown emitters around 2 eV.\n\nThe soft spot is the dielectric rescaling. The paper's own local-field calculation in Fig. 4(c,d) fixes epsilon around 2.69, but it adopts epsilon=3.52 to bring CNCB3's polarizability close to the experimental 150 cubic angstroms. That is a post-hoc rescaling, not a prediction; with epsilon=2.69 the polarizability is 96.5 cubic angstroms, about a 36% underestimate, and the dipole moments shift by roughly 30% as well. The paper is honest about the thickness dependence, but that honesty cuts the other way: it shows the quantitative predictions are convention-dependent. There are no error bars, no code or data, and the text acknowledges that alpha for CB changes from 53 to 256 cubic angstroms depending on slab thickness. That is a convergence red flag.\n\nNone of this kills the qualitative rule, which is robust to the rescaling. But the specific assignment of CNCB3 to the ~2 eV emitter and the quoted numerical values are not supported by the evidence presented. The conclusion's 'exact defect structures' language overstates what the body shows.\n\nFor what it is worth, the Berry-phase comparison for BN is a nice addition, even if the result does not match. The authors clearly know the literature and are trying to engage with experiment.\n\nWho is this for? People working on hBN color centers who want a starting point for defect assignment, and DFT groups who want a cautionary tale about local-field conventions. It deserves a serious referee, but the referee should push for a defensible local-field factor, convergence tests, and data deposition.\n\nI would send it to review rather than desk reject, with the expectation of a substantial revision.","headline":"The symmetry classification of Stark shifts is the durable contribution; the CNCB3 match to experiment is a post-hoc dielectric rescaling, not a prediction.","tokens_in":10737,"tokens_out":2487,"would_cite":false,"duration_ms":23841,"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":"Defect symmetry dictates linear vs quadratic Stark shifts in hBN","keywords":["hexagonal boron nitride","Stark shift","single-photon emitters","point-group symmetry","zero-phonon line","dielectric screening","density functional theory","quantum defects"],"falsifier":"Measure the Stark shift of a single hBN emitter whose defect structure has been independently identified, for example by electron microscopy or spin resonance, and check whether the shift is linear for a Cs-symmetry defect and purely quadratic for a D3h-symmetry defect, while also reproducing the predicted polarizability magnitude within the factor-of-two dielectric uncertainty.","tokens_in":1734,"feed_emoji":"⚡","tokens_out":1715,"duration_ms":58714,"temperature":0.7,"pith_summary":"This paper seeks to explain why single-photon emitters in hexagonal boron nitride show different Stark shifts, some linear in the applied electric field and some quadratic, by attributing the behavior to the local symmetry of the emitting defect. It argues that centrosymmetric D3h defects (CB and CNCB3) produce strong quadratic shifts, that defects with Cs symmetry and out-of-plane distortion (CBON, BN, and B_DB) produce linear shifts with dipole changes around one debye, and that C2v defects give weaker quasi-quadratic shifts. If correct, measuring the shape and magnitude of the Stark shift around 2 eV would let experimenters infer the symmetry class, and sometimes the specific structure, of an emitter whose microscopic identity is otherwise unknown. The paper also maintains that the local dielectric environment substantially rescales the extracted dipole and polarizability values, so the same defect can look quantitatively different in different samples.","feed_headline":"Defect symmetry decides linear vs quadratic Stark shift in hBN","feed_subtitle":"Calculations tie ~2 eV single-photon emitters to specific defect structures through their electric-field response.","key_machinery":"The load-bearing objects are the defect point-group symmetries (D3h, Cs, and C2v) together with the out-of-plane distortion that breaks mirror symmetry and creates a permanent dipole along the stacking direction. These determine whether the zero-phonon-line shift, written as $\\Delta ZPL = -\\Delta\\mu_z E_z - \\tfrac12 E_z \\Delta\\alpha_z E_z$, is dominated by the linear term (permanent dipole change) or the quadratic term (polarizability change). The calculations use slab models with the defect embedded in the central layer, a moving average of the electrostatic potential to extract the effective local field, and a dielectric rescaling to convert the applied slab field into the field the defect actually experiences.","core_discovery":"The central claim is a symmetry-to-response mapping: the point-group symmetry of a defect in hexagonal boron nitride determines whether its zero-phonon-line Stark shift is linear or quadratic, and the fitted coefficients identify the defect. For the defect set studied, the calculations yield quadratic transition polarizabilities of about 53 and 96 cubic ångströms for the D3h defects CB and CNCB3, linear dipole changes of about 1.1 to 1.4 debye, with sign encoding the direction of out-of-plane distortion, for the Cs defects CBON, BN, and B_DB, and smaller quasi-quadratic polarizabilities for C2v defects. The VNCB defect is singled out as a candidate for the experimentally observed V-shaped Stark response because it distorts out of plane in the ground state and relaxes to a planar configuration in the excited state. The authors present the mapping as a step toward identifying unknown emitters near 2 eV and toward using defects as local dielectric sensors.","pith_inferences":["If the symmetry-to-Stark-shift mapping holds, Stark spectroscopy could serve as a rapid symmetry assay for unknown emitters, and intermediate or mixed behavior would signal either several emitting defects in one spot or a field-induced symmetry breaking such as a dynamic Jahn-Teller effect.","The strong dependence of extracted parameters on the dielectric rescaling suggests that deliberately measuring the same emitter in flakes of different thickness could turn the present uncertainty into a probe of the local screening length.","A testable extension would be to apply an in-plane electric field to a C2v defect: the field should break the remaining mirror symmetry and convert the quasi-quadratic shift into a linear shift, a prediction that could be checked with currently available gated hBN devices."],"forward_implications":["Linear Stark shifts around 2 eV point to noncentrosymmetric defects with out-of-plane distortion, and the sign of the slope indicates the direction of the distortion.","Quadratic Stark shifts point to centrosymmetric D3h defects, with CB and CNCB3 distinguishable by their fitted polarizabilities of about 53 versus 96 cubic ångströms (or about 98 versus 161 cubic ångströms under the larger dielectric constant).","C2v defects should show weaker quasi-quadratic shifts, so their small quadratic response is a marker of lower symmetry without a permanent out-of-plane dipole.","The choice of local dielectric constant changes predicted polarizabilities by up to roughly a factor of two, so quantitative comparison with experiment requires knowing or measuring the local screening environment.","The experimentally observed V-shaped Stark response can be produced by a defect whose ground state is distorted out of plane but whose excited state is planar, linking a distinctive line shape to a specific structure."],"supporting_citations":[{"why":"Provides the experimental Stark-tuning data for hBN single-photon emitters whose linear and quadratic coefficients the paper compares against.","marker":"[13]"},{"why":"Reports the room-temperature giant Stark effect in a van der Waals material, serving as the baseline for the predicted 30 meV/(V/nm) response of BN.","marker":"[14]"},{"why":"Reports very large and reversible Stark-shift tuning in layered hBN, supplying an experimental polarizability scale for comparison.","marker":"[15]"},{"why":"Supplies the first-principles method for fitting the Stark-shift dipole and polarizability of a defect and for evaluating the effective local field.","marker":"[42]"},{"why":"Provides the moving-average electrostatic-potential approach used to extract the effective local field inside the slab.","marker":"[43]"},{"why":"Supplies the dielectric-constant formalism and the fixed-layer-thickness definition that underlie the dielectric rescaling in Eq. (1).","marker":"[41]"},{"why":"Identifies symmetric carbon tetramers (CNCB3) as spin qubits in hBN, giving the paper its D3h defect candidate with an ~2 eV transition.","marker":"[34]"},{"why":"Proposes the CBON defect with a matching hyperfine constant, giving the paper its Cs-symmetry candidate with an out-of-plane distortion.","marker":"[35]"},{"why":"Proposes dangling bonds in hBN as single-photon emitters, supplying the B_DB defect model used in the linear-Stark-shift group.","marker":"[36]"},{"why":"Provides the defect-state assignments for VNCB in hBN, including the out-of-plane ground-state configuration used for the V-shape Stark analysis.","marker":"[19]"}],"fun_headline_variants":["Defect symmetry sets linear or quadratic Stark shift in hBN","hBN defect symmetry controls Stark shift response type","Symmetry rules Stark shift in hBN single-photon emitters","How defect symmetry shapes hBN's Stark shift","Stark shift in hBN defects: symmetry is key"],"cache_read_input_tokens":12800,"weakest_assumption_plain":"The quantitative predictions rest on the rescaling that converts the electric field applied to the slab into the effective field felt by the defect, and the paper's own dielectric constant varies from 2.69 to 3.52 depending on how the layer thickness is defined, which changes the inferred dipole changes and polarizabilities by up to a factor of about two.","fun_headline_variants_meta":{"raw":{"variants":["Defect symmetry sets linear or quadratic Stark shift in hBN","hBN defect symmetry controls Stark shift response type","Symmetry rules Stark shift in hBN single-photon emitters","How defect symmetry shapes hBN's Stark shift","Stark shift in hBN defects: symmetry is key"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000582,"raw_usage":{"total_tokens":2712,"prompt_tokens":891,"completion_tokens":1821,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":507,"completion_tokens_details":{"reasoning_tokens":1741}},"tokens_in":507,"tokens_out":1821,"duration_ms":12116,"temperature":1.0,"reasoning_tokens":1741,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T19:06:57.985213+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Stark shift of a single hBN emitter whose defect structure has been independently identified, for example by electron microscopy or spin resonance, and check whether the shift is linear for a Cs-symmetry defect and purely quadratic for a D3h-symmetry defect, while also reproducing the predicted polarizability magnitude within the factor-of-two dielectric uncertainty.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental Stark-tuning data for hBN single-photon emitters whose linear and quadratic coefficients the paper compares against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the room-temperature giant Stark effect in a van der Waals material, serving as the baseline for the predicted 30 meV/(V/nm) response of BN."},{"cited_title":"Nikolay, N","cited_arxiv_id":null,"evidence_quote":"Reports very large and reversible Stark-shift tuning in layered hBN, supplying an experimental polarizability scale for comparison."},{"cited_title":"Alaerts, Y","cited_arxiv_id":null,"evidence_quote":"Supplies the first-principles method for fitting the Stark-shift dipole and polarizability of a defect and for evaluating the effective local field."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the moving-average electrostatic-potential approach used to extract the effective local field inside the slab."},{"cited_title":"Laturia, M","cited_arxiv_id":null,"evidence_quote":"Supplies the dielectric-constant formalism and the fixed-layer-thickness definition that underlie the dielectric rescaling in Eq. (1)."},{"cited_title":"Barcza, and V","cited_arxiv_id":null,"evidence_quote":"Identifies symmetric carbon tetramers (CNCB3) as spin qubits in hBN, giving the paper its D3h defect candidate with an ~2 eV transition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proposes dangling bonds in hBN as single-photon emitters, supplying the B_DB defect model used in the linear-Stark-shift group."},{"cited_title":"Sajid, J","cited_arxiv_id":null,"evidence_quote":"Provides the defect-state assignments for VNCB in hBN, including the out-of-plane ground-state configuration used for the V-shape Stark analysis."}],"review_version":1}