Laser Direct Writing and Raman Stokes-Contrast Screening of Quantum Emitter Sites in hBN

Laser Direct Writing and Raman Stokes-Contrast Screening of Quantum Emitter Sites in hBN | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Laser Direct Writing and Raman Stokes-Contrast Screening of Quantum Emitter Sites in hBN Tadas Paulauskas, Julius Janušonis, Edgaras Markauskas, Viktorija Nargelienė, and 6 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8436272/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 26 Mar, 2026 Read the published version in Discover Nano → Version 1 posted 9 You are reading this latest preprint version Abstract Laser direct writing (LDW) enables the spatially defined creation of room-temperature single-photon emitters (SPEs) in hexagonal boron nitride (hBN). However, the rapid characterization of written sites remains a bottleneck, and the available toolset for efficient screening is limited. Here, we demonstrate a streamlined LDW workflow utilizing single-shot pulses combined with a confocal screening technique that exploits the hBN E 2 g ​ Stokes line to rapidly localize and map laser-modified regions without relying a priori on defect photoluminescence (PL). This approach enables the direct correlation of site morphology with PL hotspots, revealing that the emergence of single-photon emitters coincides with a threshold regime of minimal lattice modification. Micro-Raman spectral mapping further uncovers localized compressive strain surrounding these emission sites. We classify the generated defects into two families: narrowband "red" emitters (650–750 nm) with weak phonon sidebands (PSB), and 600–650 nm emitters with stronger vibronic coupling, both exhibiting linear polarization and high single-photon purity. These results establish a practical protocol for rapid prototyping, offering a valuable addition to the characterization toolkit for scalable quantum nanophotonics. laser direct writing single-photon emitters micro-Raman hBN quantum technologies Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 1 Introduction Femtosecond laser direct writing (LDW) has emerged as a versatile technique for the maskless, scalable fabrication of nanophotonic devices and functional quantum structures [ 1 , 2 ]. By confining energy deposition to laser focal or even sub-diffraction volumes, LDW allows for the precise modification of wide-bandgap semiconductors without the need for foreign-species implantation or complex lithography. Among solid-state hosts, hexagonal boron nitride (hBN) is of particular interest for LDW-based quantum prototyping, as it hosts room-temperature single-photon sources (SPS) with brightness and stability rivaling the nitrogen-vacancy (NV) center in diamond [ 3 – 6 ]. When integrated via LDW, these emitters offer a pathway toward compact components for quantum key distribution (QKD) and quantum random number generation (QRNG) [ 7 , 8 ]. However, a common bottleneck in the LDW workflow is the rapid characterization of successfully written sites. While recent studies have demonstrated high-yield defect creation using femtosecond pulses, the process remains inherently stochastic near the optical breakdown threshold [ 9 – 12 ]. Conventional characterization relies on time-consuming photoluminescence (PL) mapping to locate active emitters within large arrays, often followed by cumbersome correlation of PL hotspots with sample morphology, thus hindering the rapid prototyping cycle required for device integration. Furthermore, while the structural origins of known hBN defects, such as the boron vacancy (V B − ​) or carbon dimer (C B ​C N ​), are established, the microscopic nature of the technologically attractive narrowband "red" emitters (650–750 nm) generated by LDW remains debated [ 17 – 21 ]. Understanding the interplay between laser-induced lattice strain, damage thresholds, and spectral purity is essential for maturing LDW into a reliable manufacturing tool. Here, we demonstrate a practical LDW and characterization workflow that addresses these challenges. We introduce an optical screening technique that exploits the intrinsic hBN E 2 g ​ Raman Stokes line to rapidly flag laser-modified regions without relying on defect fluorescence a priori [ 22 – 25 ]. By operating in the near-threshold regime, we isolate high-purity single-photon emitters where the lattice modifications are minimal yet discernible via Stokes-based screening. Spectroscopic analysis classifies these LDW-generated defects into two groups: a class of narrowband red emitters (650–750 nm) with weak phonon sidebands (PSB), and a secondary family (600–650 nm) with stronger PSB. Polarization-resolved measurements confirm linear dipole emission for both families and show co-existing unpolarized broad background in some LDW-modified sites. Furthermore, micro-Raman spectral mapping reveals localized compressive strain, which we find correlates with reduced emitter formation probability [ 26 , 27 ]. These results define a robust approach for the rapid prototyping and qualification of LDW-based quantum emitters for next-generation hBN nanophotonics. 2 Materials and Methods hBN flakes (HQ Graphene) were mechanically exfoliated onto 0.5 mm fused-silica substrates. Substrates were cleaned in acetone/isopropanol/DI water and oxygen-plasma treated prior to exfoliation. After exfoliation, samples were baked in air at 500°C for 10 min to reduce tape residues. Laser writing was performed with a 1030 nm femtosecond source (Pharos, Light Conversion, 250 fs, 602 kHz) operated in a single-shot mode using an internal pulse picker. The linearly polarized beam was attenuated to the desired pulse energy and focused through a 50×, NA 0.65 objective to a ~ 1.6 µm spot. Sites were written under ambient conditions with 5–7 µm spacing using motorized translation stages. Following writing, samples were annealed in flowing Ar (~ 200 mL/min): 20°C min⁻¹ ramp to 1000°C, 90 min dwell, then cooled to room temperature. Surface morphology was characterized by tapping-mode AFM using silicon cantilevers. Confocal fluorescence imaging and spectroscopy were performed on a custom inverted microscope using 532 nm CW excitation (typically 1 mW at the sample). Samples were scanned with a piezo stage (0.18 µm pixels, 1 ms dwell) and a 75 µm pinhole. Emission was collected through a long-pass filter (HQ545LP, Chroma) and detected using a single-photon avalanche photodiode (Tau-SPAD-50) for imaging/count-rate measurements. Spectra were acquired using a grating spectrometer with an EMCCD (Shamrock SR-303i (50 gr/mm) + Andor DU-897E-CS0-UVB), averaged over 10 acquisitions with 1 s integration each. Polarization-resolved spectra were recorded by rotating a linear polarizer in the detection path in 10° steps and averaging 10 spectra per angle (0.1 s each). \(\:{E}_{2g}\) -based screening imaging, second-order correlation, and fluorescence lifetime measurements were performed on a separate custom confocal microscope (50 µm pinhole, NA 0.90 air objective). \(\:{g}^{\left(2\right)}\left(\tau\:\right)\) measurements were carried out under CW excitation (520 nm, Integrated Optics) using a Hanbury Brown–Twiss configuration with two single-photon detectors (Excelitas SPCM-AQRH-14); timing was recorded with a time-correlated single-photon counting module Picoharp 300. Lifetimes were measured using pulsed 520 nm excitation (Picoquant LDH-D-C-520) and TCSPC detection. Details of \(\:{g}^{\left(2\right)}\left(\tau\:\right)\) processing and fitting are provided in the Supplementary Information. Raman spectroscopy and mapping were performed using 532 nm excitation with a 2400 lines/mm grating and a high-NA objective using Renishaw inVia. Raman maps were acquired with 0.5 × 0.5 µm step size and 6 s integration per point at ~ 3 mW incident power. Wavenumbers were calibrated using the Si line at 520.7 cm⁻¹. Baseline correction and peak/intensity mapping processing are described in the SI. 3 Fabrication and Raman-Based Screening of LDW Sites 3.1 Morphological Evolution and Surface Topography All samples were processed using single-shot 1030 nm femtosecond LDW pulses. Pulse energies spanned 90–10 nJ, with finer increments employed to refine the near-threshold regime (< 50 nJ) where single emitters first appear. Optical micrographs acquired pre- and post-annealing (Fig. 1 ) illustrate the morphological evolution. In the representative hBN flake, seven columns (89.5 to 35.6 nJ) remain discernible prior to annealing (Fig. 1 a). Post-anneal (Fig. 1 b), higher-energy LDW sites enlarge markedly, whereas the lowest-energy columns show negligible change. AFM topography of the central region (Fig. 1 c–d) characterizes these features as shallow craters (~ 20–40 nm deep) whose lateral dimensions scale with pulse energy (~ 2–5 µm). At the lowest energies, only few-nanometer surface steps are observed, which often vanish upon annealing. Since the flake thickness (~ 160 nm) exceeds the crater depth, the modified regions remain well-confined above the substrate. Notably, no rim pile-up or debris is evident, and some larger craters exhibit hexagonal outlines. Similar trends were observed across hBN flakes 50–500 nm thick. We hypothesize that laser-induced disorder and high local curvature accelerate hBN sublimation during the high-temperature anneal, potentially via volatile oxide formation, resulting in crater enlargement without material redeposition. 3.2 Stokes-Contrast Mapping of LDW sites Most LDW-modified regions in hBN, and indeed other morphological features such as grain boundaries or thickness variations, typically exhibit weak or negligible PL signal. Consequently, sub-bandgap excitation offers little to no contrast. Therefore, correlating defect PL with material morphology usually requires aligning pre-acquired images (e.g., optical or SEM) with subsequent scanning PL maps. Here, we introduce a practical confocal method to locate and map LDW-modified regions without relying on a priori defect PL. A key strength of this method is its compatibility with standard confocal scanning PL microscopes, requiring minimal optical layout modifications. The contrast exploits the hBN E 2 g​ Raman Stokes line (~ 560 nm for 520 nm excitation, ~ 573 nm for 532 nm), which is strongly co-polarized with the laser pump due to the in-plane isotropic hBN Raman tensor, alongside a slight leakage of p -polarized light through the long-pass dichroic mirror [ 22 – 25 ]. By setting the excitation to p -polarization with a half-wave plate, the E 2 g​ signal transmits through the dichroic while back-reflected laser light is blocked by a 550 nm long-pass filter. A simplified minimal optical layout for this E 2 g ​-based screening is shown in Fig. 2 c. This E 2 g ​-based mapping yields a bright bulk background with darker laser-written sites (Fig. 2 a), primarily because the Raman signal scales with flake thickness. Contrast is further enhanced by crater-induced depolarization and reduced collection efficiency in deformed regions. Additionally, thickness-dependent oscillations arising from pump standing waves and Stokes out-coupling effects may modulate the signal [ 24 ]. Defect PL above the dichroic cutoff (605 nm) may still appear in some spots if no short-pass filters are employed. Switching the excitation to s -polarization suppresses E 2 g ​ transmission, eliminating the bright background and revealing the defect PL directly (Fig. 2 b). If control over excitation polarization is required for defect PL analysis without the E 2 g signal, an additional long-pass filter (e.g., at 575 nm) can simply be inserted into the PL collection path. As shown in Fig. 2 a, sites created with pulse energies down to 35.6 nJ are readily identifiable via this method. Notably, below this energy threshold, neither Stokes-contrast nor defect PL is observed. This correlation remained consistent across various investigated hBN flakes, suggesting that lattice modifications required to host quantum emitters are linked to the structural changes detectable via Raman screening. A site hosting single emitter in this field of view is marked in Fig. 2 b and imaged at higher magnification in Fig. 2 d. Its spectrum (Fig. 2 e) has ZPL at 645 nm with a weak PSB near ~ 700 nm (“red” class of emitters). The E₂g peak, marked in Fig. 2 e, underpins the contrast formation in Fig. 2 a. Power dependence of the SPE follows the standard saturation model, \(\:I\left(P\right)={I}_{sat}P/(P+{P}_{sat})\) , giving \(\:{P}_{sat}\approx\:0.85\) mW and \(\:{I}_{sat}\) ≈ 2.8×10 5 cts/s, without polarization optimization and collection efficiency correction. Second-order correlation function, g ( ² ) (τ), was measured using 100 µW laser power at the objective (Fig. 2 f) and shows g ( ² ) (0) well below 0.5, indicating high single-photon purity. The noticeable photon bunching behavior is consistent with a minimal three-level system and was fitted with \(\:{g}^{\left(2\right)}\left(\tau\:\right)=1-\left(1+a\right){e}^{-{\lambda\:}_{1}\tau\:}+a{e}^{-{\lambda\:}_{2}\tau\:}\) (see SI). 3.3 Raman Spectral Imaging and Local Lattice Strain Micro-Raman spectral imaging was performed on the LDW-modified regions to further quantify lattice modifications (Fig. 3 ). The E₂ g phonon mode showed spatially varying shifts and intensity changes, albeit no signatures of a secondary hBN phase (e.g., cubic-BN) were detected. A representative spectrum from a written site is shown in Fig. 3 a. Figure 3 b illustrates a typical E₂ g peak shift at written spots relative to an unaffected reference (~ 1367 cm⁻¹). Figure 3 c shows the mapped area, with E₂ g intensity and peak-position maps in Fig. 3 d–e derived from background-removed, peak-fitted spectra (see SI). The E₂ g intensity generally dips at crater centers but can rise on rims and at the lowest-energy sites. Variations for the lowest-energy sites may point to the nonlinear thickness-dependent Raman response [ 24 ]. Notably, we observe localized positive E₂ g shifts (up to + 2 cm⁻¹) near larger crater bottoms. Interpreting the E₂ g shift, ( \(\:\varDelta\:\omega\:\) ), as a strain gauge, we estimate a compressive strain of ~ 0.1%, assuming a biaxial strain model ( \(\:{\epsilon\:}=-\varDelta\:\omega\:/2\gamma\:{\omega\:}_{0}\) ) and literature Grüneisen coefficient \(\:\gamma\:=1.04\) for multilayer hBN [ 26 ]. Under uniaxial strain the conversion factor is smaller, but the sign still indicates compression. We emphasize that the measured values are depth-averaged over the confocal volume of both pump and Stokes light, potentially underestimating the peak strain in localized LDW-affected layers. These strain patterns indicate coexisting stress fields in laser-modified hBN. Interestingly, we observe that high-intensity defect PL hotspots typically occur within the central regions of LDW sites, rather than in the strained peripheral zones (Fig. 3 b). This spatial distribution suggests that the likelihood of emitter formation during annealing is reduced at compressively strained sites. Since strain strongly modulates defect formation energies, the local compressive field may favor simple vacancies while increasing the formation energy barriers for the interstitial or complex substitutional defects required for single-photon emission in the visible spectral region. 4 Spectroscopic Characterization 4.1 Emitter Classification and Spectral Inhomogeneity We next analyze PL from LDW–fabricated single emitters. Figure 4 gathers spectra of two spectral families distinguished from typical “green” emitters by their markedly more pronounced ZPL and suppressed PSB. We further classify these families by ZPL position and PSB strength: (i) red emitters at ~ 650–750 nm with a clearly weak PSB (Fig. 4 a), and (ii) emitters at ~ 600–650 nm with a stronger PSB (Fig. 4 b). The excitation polarization was fixed for all measurements, thus, relative intensities may reflect excitation-coupling biases for differently oriented dipoles. To compare vibronic structures, we normalize each spectrum and plot them on an energy axis aligned to the ZPL (Fig. 4 c–d). Both families exhibit PSB features ~ 170–180 meV from the ZPL, consistent with coupling to high-DOS phonons of bulk hBN [ 20 , 28 ]. Regarding microscopic origin, DFT studies have proposed carbon-related centers (e.g., carbon trimers or donor–acceptor pairs) as candidates for weak-PSB red emitters, which is broadly consistent with our spectral data [ 18 , 19 ]. Despite these shared characteristics, we observe a significant spread in ZPL positions and linewidths (FWHM) across the dataset. This spectral diversity is not unique to the LDW technique but is a pervasive feature of emitters in hBN, routinely reported across various fabrication methodologies including electron irradiation, ion implantation, and thermal annealing [ 5 , 6 , 9 – 12 ]. This universality suggests that the variations stem from local environmental heterogeneity, such as strain fields and charge fluctuations from nearby defects, as well as the potential existence of multiple stable atomic configurations (e.g., distinct defect-impurity complexes or defect pairs with varying separations) that possess slightly distinct transition energies [ 28 , 29 ]. 4.2 Excited state lifetimes Excited-state lifetimes for the red emitters were quantified using pulsed 520 nm laser excitation (Fig. 5 ). The decay traces are well-described by a single-exponential model, yielding lifetimes in the ~ 2–5 ns range (Table 1). The single-exponential behavior indicates a dominant radiative decay channel characteristic of high-purity quantum emitters. For completeness, Table 1 also lists ZPL positions and Lorentzian ZPL FWHM. We observe no systematic correlation between lifetime, ZPL position, and FWHM. The variation is likely attributable to differences in the local density of optical states (LDOS) arising from the emitters' varying depths within hBN or proximity to the dielectric interface, and not necessarily intrinsic to the defect structure itself [ 30 ]. 4.3 Spectrally Resolved PL Polarization Spectrally resolved PL polarization measurements were performed on both emitter families. Representative datasets for strong-PSB (Fig. 6 a–c) and weak-PSB (Fig. 6 d–f) cases were acquired by stepping a linear polarizer in 10° increments. Both emitters exhibit clear linear-dipole behavior: the angle dependence of the ZPL-integrated intensity (background not subtracted; integration bands marked in Fig. 6 c,f) shows pronounced maxima/minima, and the PSB-integrated intensity follows the same angle without a visible offset. This alignment implies that the dominant phonon-assisted transitions inherit the ZPL dipole orientation, suggesting that electron–phonon coupling does not strongly rotate the optical selection rules for these modes. Note, however, that 10° steps can miss the true extrema, so maxima/minima need not be separated by exactly 90°. The degree of linear polarization (DoLP) for hBN emitters seldom reaches unity, as noted in the literature [ 5 , 6 , 31 – 33 ]. Several mechanisms can contribute to incomplete DoLP: (i) intrinsic PL depolarization processes, (ii) an unpolarized background that sets a floor at the polarization minimum, or (iii) experimental geometry effects, such as the projection of an out-of-plane dipole component. Our spectrally resolved measurements identify broad, ZPL-free emission bands peaking in the 580–600 nm window as the primary unpolarized background for the emitters in Fig. 6 c,f. This background is detrimental to QKD protocols because it degrades single-photon purity, likely explaining why \(\:{g}^{\left(2\right)}\left(0\right)\) values saturate near ~ 0.2 for some yellow–red hBN emitters [ 5 , 6 ]. Consequently, identifying and suppressing these broad emission defects is critical for future fabrication optimization, or suggests a strategic shift to NIR emitters and using longer-wavelength excitation to avoid this spectral window. Additionally, the incomplete suppression of ZPL and PSB features at polarization minima (Fig. 6 f) points to a combination of intrinsic depolarization and residual background (mechanisms i and ii) limiting the contrast. 4.4 Multiple-Emitter Sites and Spectral Statistics We also examined fs-written sites that host multiple emitters per site, which are common at higher pulse energies, as well as sites that show broad, ZPL-free emission. The latter appears frequently peaking at ~ 580–600 nm (Fig. 7 a), shows no polarization dependence (cf. Figure 6 ), and can persist as a broad background for polarized emitters, as discussed above. Exciting such sites with a 637 nm laser does not reproduce them next to the E₂g Stokes line (at ~ 700 nm), ruling them out as, e.g., defect-based Raman modes. A plausible origin is PL from a defect center, or a family of defects, with a large Huang–Rhys factor and strong electron–phonon coupling. Vacancy-based emitters tend to have a large Huang–Rhys factor, one such well known example is \(\:{V}_{\text{B}}^{-}\) in hBN, which exhibits a broad ZPL-free band near 800 nm [ 13 – 15 ]. At higher laser pulse energies, spectra frequently display multiple ZPLs within a single site (Fig. 7 b), typically unresolved spatially in confocal scans (examples in Fig. 7 c–d). Interestingly, some multiple-emitter sites show repeated exact ZPL energies (up to our spectral resolution), e.g., at ~ 675 and ~ 695 nm (Fig. 7 b). Given the general ZPL spread discussed earlier, such recurrence may suggest identical microscopic centers under very similar local strain. We note that in our experimental study we did not employ energies that create exceedingly large craters in hBN, as some previous studies have done. In those higher-energy cases, the craters show pronounced material pileups and folded layers at crater edges, which tend to host a high density of defects and produce an almost continuous bright PL rim. Here, in contrast, for the multiple-emitter sites at larger energies the emission originates mostly from the central spot. This can be seen in Fig. 7 c–d as larger hot-spot emission. Future studies will aim to carry out ODMR on LDW-created SPEs and such multiple-emitter hot spots, which may be interesting for magnetic sensing applications and to help deduce their atomic origin. Figure 7 e compiles ZPL statistics across flakes, incorporating both single and multiple-emitter sites. Emitters occur most frequently in the 580–650 nm range, where the average spectrum correspondingly exhibits the highest PL intensity. A secondary clustering is observed in the red band (650–725 nm), accompanied by a distinct spectral feature in the average intensity around 675–700 nm. Considering the typical spread of ZPLs in hBN, this secondary clustering likely indicates a specific "red" defect structure, or a family of defects with similar transition energies, that is preferentially generated by the LDW and annealing process. 5 Conclusions In conclusion, we have demonstrated a workflow for femtosecond laser direct writing (LDW) of quantum emitters in hBN that facilitates the rapid post-process identification of written sites. Using an optical screening modality based on the intrinsic hBN E 2 g ​ Raman Stokes line, we map laser-modified regions without relying a priori on defect photoluminescence, adding a simple and broadly implementable tool to the LDW characterization toolkit. This approach reveals that the emergence of single-photon emitters coincides with a threshold regime of minimal lattice modification, whereas higher pulse energies more frequently yield multi-emitter spectra and broadband backgrounds. The generated emitters fall into two families: narrowband red emitters (650–750 nm) with weak phonon sidebands and a 600–650 nm family with stronger vibronic coupling, both exhibiting linear polarization and high single-photon purity. However, a broadband unpolarized background centered at 580–600 nm can persist in some sites, degrading polarization contrast. Micro-Raman mapping reveals localized E 2 g ​​ shifts consistent with compressive strain, yet these depth-averaged fields are insufficient to explain the ZPL spectral diversity, pointing to atomic-scale environmental heterogeneity as the dominant factor. Together, these results establish LDW combined with Raman-based screening as a practical characterization approach for qualifying hBN quantum emitter sites, paving the way further advancements on LDW-based integrated nanophotonic devices. Declarations Supplementary information Supporting Material 1. Acknowledgments The Research Council of Lithuania for their financial support. Author contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Funding Funded by the Research Council of Lithuania, Project No. P-ITP-24-22. Data availability Data will be available on request. Ethics approval and consent to participate Not applicable. Consent to publication All authors read and approved the manuscript for publication. Competing interests The authors declare no competing interests. References Sugioka K, Cheng Y (2014) Ultrafast lasers—reliable tools for advanced materials processing. Appl Phys Rev 1:041303. https://doi.org/10.1063/1.4904320 Gattass RR, Mazur E Femtosecond laser micromachining in transparent materials. Nat Photonics 2008, 2, 219–225. https://doi.org/10.1038/nphoton.2008.47 Doherty MW, Manson NB, Delaney P, Jelezko F, Wrachtrup J, Hollenberg LC (2013) L. The nitrogen-vacancy colour centre in diamond. Phys Rep 528:1–45. https://doi.org/10.1016/j.physrep.2013.02.001 Aharonovich I, Tetienne J-P, Toth M (2016) Quantum emitters in hexagonal boron nitride. Phys Today 69(10):38–44 Tran TT, Elbadawi C, Totonjian D, Lobo CJ, Grosso G, Moon H, Englund DR, Ford MJ, Aharonovich I, Toth M (2016) Robust multicolor single-photon emission from point defects in hexagonal boron nitride. ACS Nano 10(8):7331–7338 Mendelson N et al (2017) Tunable and high-purity room temperature single-photon emission from atomic defects in hexagonal boron nitride. Nat Commun 8:15720. https://doi.org/10.1038/s41467-017-00810-2 Scarani V, Bechmann-Pasquinucci H, Cerf NJ, Dušek M, Lütkenhaus N, Peev M (2009) The security of practical quantum key distribution. Rev Mod Phys 81:1301–1350. https://doi.org/10.1103/RevModPhys.81.1301 Herrero-Collantes M, Garcia-Escartin JC (2017) Quantum random number generators. Rev Mod Phys 89:015004. https://doi.org/10.1103/RevModPhys.89.015004 Gan L, Zhang D, Zhang R, Zhang Q, Sun H, Li Y, Ning C-Z (2022) Large-scale, high-yield laser fabrication of bright and pure single-photon emitters at room temperature in hexagonal boron nitride. ACS Nano 16:14388–14396. https://doi.org/10.1021/acsnano.2c04386 Gao X, Li G, Li S et al (2021) Femtosecond laser writing of spin defects in hexagonal boron nitride. ACS Photonics 8:994–1000. https://doi.org/10.1021/acsphotonics.0c01847 Yang Y-Z, Zhu T-X, Li Z-P et al (2023) Laser direct writing of visible spin defects in hexagonal boron nitride for applications in spin-based technologies. ACS Appl Nano Mater 6:14190–14197. https://doi.org/10.1021/acsanm.3c01047 Ahmad SI, Dave A, Sarpong E et al (2023) Dielectric breakdown and sub-wavelength patterning of monolayer hexagonal boron nitride using femtosecond pulses. 2D Mater 10:045029. https://doi.org/10.1088/2053-1583/acfa0f Gottscholl A et al (2020) Initialization and readout of intrinsic spin defects in a van der Waals crystal at room temperature. Nat Mater 19:540–545 Gottscholl A et al (2021) Room-temperature coherent control of spin defects in hexagonal boron nitride. Sci Adv 7:eabf3630 Haykal A et al (2022) Decoherence of V B - spin defects in monoisotopic hexagonal boron nitride. Nat Commun 13:4347 Mackoit-Sinkevičienė M et al (2019) Carbon dimer defect as a source of the 4.1 eV luminescence in hBN. Appl Phys Lett 115:212101 Weston L, Wickramaratne D, Mackoit M, Van de Alkauskas A (2018) Walle, C. G. Native point defects and impurities in hexagonal boron nitride. Phys Rev B 97:214104 Maciaszek M, Razinkovas L, Alkauskas A (2023) Thermodynamics of carbon point defects in hexagonal boron nitride. Phys Rev B 108:014105 Auburger P, Gali A (2023) Towards ab initio identification of paramagnetic substitutional carbon defects in hexagonal boron nitride acting as quantum bits. npj Comput Mater 9:38 Fischer M, Sajid A, Iles-Smith J et al (2023) Combining experiments with the polaron model and ab initio methods towards the microscopic origin of luminescent centres in hBN. Nanoscale 15:18638–18652 Mohajerani SS, Chen S, Alaei A et al (2022) Narrowband quantum light emission from oxygen-related color centers in hexagonal boron nitride. ACS Photonics 9:2645–2654 Cai Q, Scullion D, Han Y et al (2017) Raman signature and phonon dispersion of atomically thin hexagonal boron nitride. Phys Rev B 95:235407 Stenger I, Nagler P, Neumann A et al (2017) Low-frequency Raman spectroscopy of few-atomic-layer hBN. 2D Mater 4:031003 Schaumburg F, Sleziona S, Zöllner M et al (2023) Enhanced intensity of Raman signals from hexagonal boron nitride films. Beilstein J Nanotechnol 14:154–165 Reich S, Ferrari AC, Arenal R, Loiseau A, Bello I, Robertson J (2005) Resonant Raman scattering in cubic and hexagonal boron nitride. Phys Rev B 71:205201 Androulidakis C, Koukaras EN, Poss M, Papagelis K, Galiotis C, Tawfick S (2018) Strained hexagonal boron nitride: phonon shift and Grüneisen parameter. Phys. Rev. B 97, 241414(R) Blundo E, Surrente A, Spirito D et al (2022) Vibrational properties in highly strained hexagonal boron nitride bubbles. Nano Lett 22:1150–1157 White S, Stewart C, Solntsev AS, Li C, Toth M, Kianinia M, Aharonovich I (2021) Phonon dephasing and spectral diffusion of quantum emitters in hexagonal boron nitride. Optica 8(9):1153–1158 Fournier C, Watanabe K, Taniguchi T et al (2023) Investigating the fast spectral diffusion of a quantum emitter in hBN using resonant excitation and photon correlations. Phys Rev B 107:195304 Clua-Provost T, Durand A, Fraunié J et al (2022) Impact of thickness-dependent nanophotonic effects on the optical response of color centers in hexagonal boron nitride. ACS Photonics 9:3761–3770 Jungwirth NR, Fuchs GD (2017) Optical absorption and emission mechanisms of single defects in hexagonal boron nitride. Phys Rev Lett 119:057401. https://doi.org/10.1103/PhysRevLett.119.057401 Fournier C, Watanabe K, Taniguchi T et al (2021) Position-controlled quantum emitters with reproducible emission wavelength in hexagonal boron nitride. Nat Commun 12:3779. https://doi.org/10.1038/s41467-021-24019-6 Kumar A, Peng P, Wrachtrup J et al (2024) Polarization dynamics of solid-state quantum emitters. ACS Nano 18:10745–10768 Additional Declarations No competing interests reported. Supplementary Files SupplementaryMaterial.docx Cite Share Download PDF Status: Published Journal Publication published 26 Mar, 2026 Read the published version in Discover Nano → Version 1 posted Editorial decision: Revision requested 30 Jan, 2026 Reviews received at journal 28 Jan, 2026 Reviews received at journal 25 Jan, 2026 Reviewers agreed at journal 20 Jan, 2026 Reviewers agreed at journal 19 Jan, 2026 Reviewers invited by journal 17 Jan, 2026 Editor assigned by journal 24 Dec, 2025 Submission checks completed at journal 24 Dec, 2025 First submitted to journal 23 Dec, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8436272","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":583027040,"identity":"ce9c2fce-a156-47e4-805e-5732ef4acea5","order_by":0,"name":"Tadas Paulauskas","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABAklEQVRIiWNgGAWjYHACAwjFDiIOMDDwMx8AsSxwqueBa2GGapFsSwaxJCSI12JwjIAWe/bDGz/+YLgjb3CY+ZnkjzN29sbH+I9u+MEgUYfTFp60YmkehmeGGw6zmUnz3EhO3HaMme1mD16H5RhIMzAcZtxwmMFMmuEDc4LZ/Wa223j9wv/G+OcPhsP2Gw6zf5P88aHe3riNmYAWiRwzCR6Gw4kbDvMAGTeA1rER0nLjWZk1j8Gz5JmHeYqtec4cT5xxjNnsZo+BhGQDDi3s/cmbb/6ouGPbd7x9480fx6rt+dsYn934UWHDj8sWCDA4gCGCXwMDOApHwSgYBaNgFOACAIw6UuBSDZocAAAAAElFTkSuQmCC","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":true,"prefix":"","firstName":"Tadas","middleName":"","lastName":"Paulauskas","suffix":""},{"id":583027041,"identity":"2ef003b3-ac5f-4afc-ab9a-2c56af037c8b","order_by":1,"name":"Julius Janušonis","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Julius","middleName":"","lastName":"Janušonis","suffix":""},{"id":583027042,"identity":"8684e667-cd9a-41c2-8f6e-2b4cac8b7ded","order_by":2,"name":"Edgaras Markauskas","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Edgaras","middleName":"","lastName":"Markauskas","suffix":""},{"id":583027043,"identity":"e2ab867b-a9de-4903-8b22-495588b13d2e","order_by":3,"name":"Viktorija Nargelienė","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Viktorija","middleName":"","lastName":"Nargelienė","suffix":""},{"id":583027044,"identity":"5ef588d3-768d-4cd6-a159-5908e1f704ef","order_by":4,"name":"Vakaris Šilys","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Vakaris","middleName":"","lastName":"Šilys","suffix":""},{"id":583027045,"identity":"6701fdda-7b54-4bc8-a620-7478192f7cfb","order_by":5,"name":"Ifra Bibi","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Ifra","middleName":"","lastName":"Bibi","suffix":""},{"id":583027046,"identity":"577bc821-4768-468e-8326-3a7c4365a8b9","order_by":6,"name":"Danielis Rutkauskas","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Danielis","middleName":"","lastName":"Rutkauskas","suffix":""},{"id":583027048,"identity":"33045171-e157-4326-9c59-cddf7f60a67d","order_by":7,"name":"Skirmantas Keršulis","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Skirmantas","middleName":"","lastName":"Keršulis","suffix":""},{"id":583027049,"identity":"de5d91eb-227a-4f90-914b-342c6d3a1898","order_by":8,"name":"Virginijus Bukauskas","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Virginijus","middleName":"","lastName":"Bukauskas","suffix":""},{"id":583027050,"identity":"bbf7a67d-f265-43ee-b33e-98e1157a69e8","order_by":9,"name":"Martynas Talaikis","email":"","orcid":"","institution":"Center for Physical Sciences and Technology","correspondingAuthor":false,"prefix":"","firstName":"Martynas","middleName":"","lastName":"Talaikis","suffix":""}],"badges":[],"createdAt":"2025-12-23 17:53:12","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-8436272/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-8436272/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s11671-026-04530-9","type":"published","date":"2026-03-26T16:08:51+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":101831276,"identity":"945f0301-4181-4da6-a9d3-2305b59c78a4","added_by":"auto","created_at":"2026-02-04 06:32:31","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":557149,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Optical micrograph of LDW hBN before annealing. Seven spot columns are visible with pulse energy decreasing right to left. (b) Same area after annealing; dashed boxes mark regions scanned in later figures. Scale bars 7 µm. (c) AFM topography of the yellow central box in (b). (d) Line profiles across spots in (c) (matching colors). Note, the vertical axis shows relative surface height, not absolute flake thickness.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/d01e871d80c4046ed3db35b5.png"},{"id":101881797,"identity":"bba28d72-4295-4f82-bf20-543769ccf064","added_by":"auto","created_at":"2026-02-04 15:16:34","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":614713,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Scanning confocal image of the region marked in Fig. 1b. (b) Same field of view (x = 0–30 µm) with the hBN Stokes contribution suppressed by rotating the half-wave plate. (c) Simplified optical layout for imaging and g\u003csup\u003e(\u003c/sup\u003e²\u003csup\u003e)\u003c/sup\u003e(τ) measurements. (d) Confocal image of a single emitter indicated in (b). (e) Fluorescence spectrum of the defect in (d) with ZPL at 645 nm and PSB at ~700 nm; inset shows laser power dependence, vertical axis in 10\u003csup\u003e5 \u003c/sup\u003ects/s. The hBN E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e Stokes peak is marked blue, used for contrast in (a). (f) Second-order correlation g\u003csup\u003e\u003cem\u003e(\u003c/em\u003e\u003c/sup\u003e²\u003csup\u003e)\u003c/sup\u003e(τ) of the defect with 3-level fit.\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/d99c05ac58ecbeb0d64ce0b7.png"},{"id":101831280,"identity":"a96c325a-0e21-4970-b01e-6e62cdf0877d","added_by":"auto","created_at":"2026-02-04 06:32:31","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":444873,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Representative Raman spectrum of a laser-written site. (b) Overlay of a reference (unaffected bulk) and a laser-written spot, showing the typical E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e peak redshift. (c) Optical micrograph of the area mapped in (d–e). (d) Raman map of E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e peak intensity, highlighting variations near laser-written sites. (e) Raman map of E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e peak position shift, revealing spatially varying redshifts. Reference E\u003csub\u003e2g\u003c/sub\u003e peak position for this map is the unaffected region on the left from LDW sites.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/91334ac61fcceca28af79c0f.png"},{"id":101881437,"identity":"6db61821-5d98-4ddd-9cfa-ad0b325923b9","added_by":"auto","created_at":"2026-02-04 15:12:09","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":364509,"visible":true,"origin":"","legend":"\u003cp\u003ePL spectra of individual emitters at laser-written sites, grouped by phonon-sideband (PSB) strength. (a) Weak PSB. (b) Strong PSB. (c–d) Same spectra as (a–b), normalized and plotted on an energy axis with each trace aligned to its ZPL. Panels (a–b) and panels (c–d) share the same x-axis scale. Measurements were done at room temperature.\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/7766f688f93e7adcde1cb40b.png"},{"id":101831282,"identity":"25e229ba-24f5-4d73-9d6b-6a995e934aad","added_by":"auto","created_at":"2026-02-04 06:32:31","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":165805,"visible":true,"origin":"","legend":"\u003cp\u003eTime-resolved PL intensity of emitters with weak PSB, obtained at room temperature. Decay traces are fitted with monoexponential (listed in Table 1).\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/2ba0a7fd26df765ecd39e622.png"},{"id":101881450,"identity":"9b79a97f-8e77-4389-a92c-0db610447d41","added_by":"auto","created_at":"2026-02-04 15:12:11","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":721751,"visible":true,"origin":"","legend":"\u003cp\u003eSpectrally resolved PL polarization for emitters with strong (a–c) and weak (d–f) PSB. (a, d) Polarization-dependent spectra (linear polarizer stepped by 10°), globally normalized. (b, e) ZPL and PSB intensities vs polarization angle, integrated over the bands marked in (c, f); background was not fitted/subtracted. (c, f) Spectra at angles corresponding to each emitter’s maximum and minimum ZPL intensity.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/190aa4eaa8d2e0981cab4575.png"},{"id":101880739,"identity":"593be52d-7773-44fd-ac75-784875ea3edf","added_by":"auto","created_at":"2026-02-04 15:05:44","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":733030,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Representative broad emissions peaking at ~580-600 nm from LDW sites; sharp feature at ~575 nm is E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e Stokes line. (b) Spectra illustrating multiple emitters within a single site and repeatedly observed red emitters with nearly identical ZPL positions. (c–d) Confocal PL images of written flakes showing bright emission hotspots. (e) ZPL histogram (8 nm bins, 580–800 nm range) and the average spectrum (solid line).\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/9d4c722fd7a61824c0f7bf07.png"},{"id":105754929,"identity":"c78d1c7b-d810-4567-9db7-d66a97a6ff28","added_by":"auto","created_at":"2026-03-30 16:23:10","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4408265,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/a87821c4-21f5-4088-ab4b-312dbaacd373.pdf"},{"id":101831283,"identity":"fd02335f-3283-429c-8249-b74eac259828","added_by":"auto","created_at":"2026-02-04 06:32:32","extension":"docx","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":1095987,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryMaterial.docx","url":"https://assets-eu.researchsquare.com/files/rs-8436272/v1/2735d7cd0078c6732fdb5636.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Laser Direct Writing and Raman Stokes-Contrast Screening of Quantum Emitter Sites in hBN","fulltext":[{"header":"1 Introduction","content":"\u003cp\u003eFemtosecond laser direct writing (LDW) has emerged as a versatile technique for the maskless, scalable fabrication of nanophotonic devices and functional quantum structures [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. By confining energy deposition to laser focal or even sub-diffraction volumes, LDW allows for the precise modification of wide-bandgap semiconductors without the need for foreign-species implantation or complex lithography. Among solid-state hosts, hexagonal boron nitride (hBN) is of particular interest for LDW-based quantum prototyping, as it hosts room-temperature single-photon sources (SPS) with brightness and stability rivaling the nitrogen-vacancy (NV) center in diamond [\u003cspan additionalcitationids=\"CR4 CR5\" citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. When integrated via LDW, these emitters offer a pathway toward compact components for quantum key distribution (QKD) and quantum random number generation (QRNG) [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e, \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eHowever, a common bottleneck in the LDW workflow is the rapid characterization of successfully written sites. While recent studies have demonstrated high-yield defect creation using femtosecond pulses, the process remains inherently stochastic near the optical breakdown threshold [\u003cspan additionalcitationids=\"CR10 CR11\" citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. Conventional characterization relies on time-consuming photoluminescence (PL) mapping to locate active emitters within large arrays, often followed by cumbersome correlation of PL hotspots with sample morphology, thus hindering the rapid prototyping cycle required for device integration. Furthermore, while the structural origins of known hBN defects, such as the boron vacancy (V\u003csub\u003eB\u003c/sub\u003e\u003csup\u003e\u0026minus;\u003c/sup\u003e​) or carbon dimer (C\u003csub\u003eB\u003c/sub\u003e​C\u003csub\u003eN\u003c/sub\u003e​), are established, the microscopic nature of the technologically attractive narrowband \"red\" emitters (650\u0026ndash;750 nm) generated by LDW remains debated [\u003cspan additionalcitationids=\"CR18 CR19 CR20\" citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. Understanding the interplay between laser-induced lattice strain, damage thresholds, and spectral purity is essential for maturing LDW into a reliable manufacturing tool.\u003c/p\u003e \u003cp\u003eHere, we demonstrate a practical LDW and characterization workflow that addresses these challenges. We introduce an optical screening technique that exploits the intrinsic hBN E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e​ Raman Stokes line to rapidly flag laser-modified regions without relying on defect fluorescence \u003cem\u003ea priori\u003c/em\u003e [\u003cspan additionalcitationids=\"CR23 CR24\" citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. By operating in the near-threshold regime, we isolate high-purity single-photon emitters where the lattice modifications are minimal yet discernible via Stokes-based screening. Spectroscopic analysis classifies these LDW-generated defects into two groups: a class of narrowband red emitters (650\u0026ndash;750 nm) with weak phonon sidebands (PSB), and a secondary family (600\u0026ndash;650 nm) with stronger PSB. Polarization-resolved measurements confirm linear dipole emission for both families and show co-existing unpolarized broad background in some LDW-modified sites. Furthermore, micro-Raman spectral mapping reveals localized compressive strain, which we find correlates with reduced emitter formation probability [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]. These results define a robust approach for the rapid prototyping and qualification of LDW-based quantum emitters for next-generation hBN nanophotonics.\u003c/p\u003e"},{"header":"2 Materials and Methods","content":"\u003cp\u003ehBN flakes (HQ Graphene) were mechanically exfoliated onto 0.5 mm fused-silica substrates. Substrates were cleaned in acetone/isopropanol/DI water and oxygen-plasma treated prior to exfoliation. After exfoliation, samples were baked in air at 500\u0026deg;C for 10 min to reduce tape residues. Laser writing was performed with a 1030 nm femtosecond source (Pharos, Light Conversion, 250 fs, 602 kHz) operated in a single-shot mode using an internal pulse picker. The linearly polarized beam was attenuated to the desired pulse energy and focused through a 50\u0026times;, NA 0.65 objective to a\u0026thinsp;~\u0026thinsp;1.6 \u0026micro;m spot. Sites were written under ambient conditions with 5\u0026ndash;7 \u0026micro;m spacing using motorized translation stages. Following writing, samples were annealed in flowing Ar (~\u0026thinsp;200 mL/min): 20\u0026deg;C min⁻\u0026sup1; ramp to 1000\u0026deg;C, 90 min dwell, then cooled to room temperature. Surface morphology was characterized by tapping-mode AFM using silicon cantilevers. Confocal fluorescence imaging and spectroscopy were performed on a custom inverted microscope using 532 nm CW excitation (typically 1 mW at the sample). Samples were scanned with a piezo stage (0.18 \u0026micro;m pixels, 1 ms dwell) and a 75 \u0026micro;m pinhole. Emission was collected through a long-pass filter (HQ545LP, Chroma) and detected using a single-photon avalanche photodiode (Tau-SPAD-50) for imaging/count-rate measurements. Spectra were acquired using a grating spectrometer with an EMCCD (Shamrock SR-303i (50 gr/mm)\u0026thinsp;+\u0026thinsp;Andor DU-897E-CS0-UVB), averaged over 10 acquisitions with 1 s integration each. Polarization-resolved spectra were recorded by rotating a linear polarizer in the detection path in 10\u0026deg; steps and averaging 10 spectra per angle (0.1 s each). \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{E}_{2g}\\)\u003c/span\u003e\u003c/span\u003e-based screening imaging, second-order correlation, and fluorescence lifetime measurements were performed on a separate custom confocal microscope (50 \u0026micro;m pinhole, NA 0.90 air objective). \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{g}^{\\left(2\\right)}\\left(\\tau\\:\\right)\\)\u003c/span\u003e\u003c/span\u003emeasurements were carried out under CW excitation (520 nm, Integrated Optics) using a Hanbury Brown\u0026ndash;Twiss configuration with two single-photon detectors (Excelitas SPCM-AQRH-14); timing was recorded with a time-correlated single-photon counting module Picoharp 300. Lifetimes were measured using pulsed 520 nm excitation (Picoquant LDH-D-C-520) and TCSPC detection. Details of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{g}^{\\left(2\\right)}\\left(\\tau\\:\\right)\\)\u003c/span\u003e\u003c/span\u003eprocessing and fitting are provided in the Supplementary Information. Raman spectroscopy and mapping were performed using 532 nm excitation with a 2400 lines/mm grating and a high-NA objective using Renishaw inVia. Raman maps were acquired with 0.5 \u0026times; 0.5 \u0026micro;m step size and 6 s integration per point at ~\u0026thinsp;3 mW incident power. Wavenumbers were calibrated using the Si line at 520.7 cm⁻\u0026sup1;. Baseline correction and peak/intensity mapping processing are described in the SI.\u003c/p\u003e"},{"header":"3 Fabrication and Raman-Based Screening of LDW Sites","content":"\u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Morphological Evolution and Surface Topography\u003c/h2\u003e \u003cp\u003eAll samples were processed using single-shot 1030 nm femtosecond LDW pulses. Pulse energies spanned 90\u0026ndash;10 nJ, with finer increments employed to refine the near-threshold regime (\u0026lt;\u0026thinsp;50 nJ) where single emitters first appear. Optical micrographs acquired pre- and post-annealing (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) illustrate the morphological evolution. In the representative hBN flake, seven columns (89.5 to 35.6 nJ) remain discernible prior to annealing (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea). Post-anneal (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb), higher-energy LDW sites enlarge markedly, whereas the lowest-energy columns show negligible change.\u003c/p\u003e \u003cp\u003eAFM topography of the central region (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec\u0026ndash;d) characterizes these features as shallow craters (~\u0026thinsp;20\u0026ndash;40 nm deep) whose lateral dimensions scale with pulse energy (~\u0026thinsp;2\u0026ndash;5 \u0026micro;m). At the lowest energies, only few-nanometer surface steps are observed, which often vanish upon annealing. Since the flake thickness (~\u0026thinsp;160 nm) exceeds the crater depth, the modified regions remain well-confined above the substrate. Notably, no rim pile-up or debris is evident, and some larger craters exhibit hexagonal outlines. Similar trends were observed across hBN flakes 50\u0026ndash;500 nm thick. We hypothesize that laser-induced disorder and high local curvature accelerate hBN sublimation during the high-temperature anneal, potentially via volatile oxide formation, resulting in crater enlargement without material redeposition.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003e3.2 Stokes-Contrast Mapping of LDW sites\u003c/h2\u003e \u003cp\u003eMost LDW-modified regions in hBN, and indeed other morphological features such as grain boundaries or thickness variations, typically exhibit weak or negligible PL signal. Consequently, sub-bandgap excitation offers little to no contrast. Therefore, correlating defect PL with material morphology usually requires aligning pre-acquired images (e.g., optical or SEM) with subsequent scanning PL maps. Here, we introduce a practical confocal method to locate and map LDW-modified regions without relying on a priori defect PL. A key strength of this method is its compatibility with standard confocal scanning PL microscopes, requiring minimal optical layout modifications. The contrast exploits the hBN E\u003csub\u003e2\u003cem\u003eg​\u003c/em\u003e\u003c/sub\u003e Raman Stokes line (~\u0026thinsp;560 nm for 520 nm excitation, ~\u0026thinsp;573 nm for 532 nm), which is strongly co-polarized with the laser pump due to the in-plane isotropic hBN Raman tensor, alongside a slight leakage of \u003cem\u003ep\u003c/em\u003e-polarized light through the long-pass dichroic mirror [\u003cspan additionalcitationids=\"CR23 CR24\" citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. By setting the excitation to \u003cem\u003ep\u003c/em\u003e-polarization with a half-wave plate, the E\u003csub\u003e2\u003cem\u003eg​\u003c/em\u003e\u003c/sub\u003e signal transmits through the dichroic while back-reflected laser light is blocked by a 550 nm long-pass filter. A simplified minimal optical layout for this E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e​-based screening is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec.\u003c/p\u003e \u003cp\u003eThis E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e​-based mapping yields a bright bulk background with darker laser-written sites (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea), primarily because the Raman signal scales with flake thickness. Contrast is further enhanced by crater-induced depolarization and reduced collection efficiency in deformed regions. Additionally, thickness-dependent oscillations arising from pump standing waves and Stokes out-coupling effects may modulate the signal [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. Defect PL above the dichroic cutoff (605 nm) may still appear in some spots if no short-pass filters are employed. Switching the excitation to \u003cem\u003es\u003c/em\u003e-polarization suppresses E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e ​ transmission, eliminating the bright background and revealing the defect PL directly (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb). If control over excitation polarization is required for defect PL analysis without the E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e signal, an additional long-pass filter (e.g., at 575 nm) can simply be inserted into the PL collection path.\u003c/p\u003e \u003cp\u003eAs shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea, sites created with pulse energies down to 35.6 nJ are readily identifiable via this method. Notably, below this energy threshold, neither Stokes-contrast nor defect PL is observed. This correlation remained consistent across various investigated hBN flakes, suggesting that lattice modifications required to host quantum emitters are linked to the structural changes detectable via Raman screening.\u003c/p\u003e \u003cp\u003eA site hosting single emitter in this field of view is marked in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb and imaged at higher magnification in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed. Its spectrum (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ee) has ZPL at 645 nm with a weak PSB near ~\u0026thinsp;700 nm (\u0026ldquo;red\u0026rdquo; class of emitters). The E₂g peak, marked in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ee, underpins the contrast formation in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea. Power dependence of the SPE follows the standard saturation model, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:I\\left(P\\right)={I}_{sat}P/(P+{P}_{sat})\\)\u003c/span\u003e\u003c/span\u003e, giving \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{P}_{sat}\\approx\\:0.85\\)\u003c/span\u003e\u003c/span\u003e mW and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{I}_{sat}\\)\u003c/span\u003e\u003c/span\u003e \u0026asymp; 2.8\u0026times;10\u003csup\u003e5\u003c/sup\u003e cts/s, without polarization optimization and collection efficiency correction. Second-order correlation function, g\u003csup\u003e(\u003c/sup\u003e\u0026sup2;\u003csup\u003e)\u003c/sup\u003e(τ), was measured using 100 \u0026micro;W laser power at the objective (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ef) and shows g\u003csup\u003e(\u003c/sup\u003e\u0026sup2;\u003csup\u003e)\u003c/sup\u003e(0) well below 0.5, indicating high single-photon purity. The noticeable photon bunching behavior is consistent with a minimal three-level system and was fitted with \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{g}^{\\left(2\\right)}\\left(\\tau\\:\\right)=1-\\left(1+a\\right){e}^{-{\\lambda\\:}_{1}\\tau\\:}+a{e}^{-{\\lambda\\:}_{2}\\tau\\:}\\)\u003c/span\u003e\u003c/span\u003e (see SI).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003e3.3 Raman Spectral Imaging and Local Lattice Strain\u003c/h2\u003e \u003cp\u003eMicro-Raman spectral imaging was performed on the LDW-modified regions to further quantify lattice modifications (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). The E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e phonon mode showed spatially varying shifts and intensity changes, albeit no signatures of a secondary hBN phase (e.g., cubic-BN) were detected. A representative spectrum from a written site is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ea. Figure\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb illustrates a typical E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e peak shift at written spots relative to an unaffected reference (~\u0026thinsp;1367 cm⁻\u0026sup1;). Figure\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec shows the mapped area, with E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e intensity and peak-position maps in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ed\u0026ndash;e derived from background-removed, peak-fitted spectra (see SI).\u003c/p\u003e \u003cp\u003eThe E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e intensity generally dips at crater centers but can rise on rims and at the lowest-energy sites. Variations for the lowest-energy sites may point to the nonlinear thickness-dependent Raman response [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. Notably, we observe localized positive E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e shifts (up to +\u0026thinsp;2 cm⁻\u0026sup1;) near larger crater bottoms. Interpreting the E₂\u003csub\u003e\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e shift, (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\varDelta\\:\\omega\\:\\)\u003c/span\u003e\u003c/span\u003e), as a strain gauge, we estimate a compressive strain of ~\u0026thinsp;0.1%, assuming a biaxial strain model (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\epsilon\\:}=-\\varDelta\\:\\omega\\:/2\\gamma\\:{\\omega\\:}_{0}\\)\u003c/span\u003e\u003c/span\u003e) and literature Gr\u0026uuml;neisen coefficient \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\gamma\\:=1.04\\)\u003c/span\u003e\u003c/span\u003e for multilayer hBN [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. Under uniaxial strain the conversion factor is smaller, but the sign still indicates compression. We emphasize that the measured values are depth-averaged over the confocal volume of both pump and Stokes light, potentially underestimating the peak strain in localized LDW-affected layers.\u003c/p\u003e \u003cp\u003eThese strain patterns indicate coexisting stress fields in laser-modified hBN. Interestingly, we observe that high-intensity defect PL hotspots typically occur within the central regions of LDW sites, rather than in the strained peripheral zones (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb). This spatial distribution suggests that the likelihood of emitter formation during annealing is reduced at compressively strained sites. Since strain strongly modulates defect formation energies, the local compressive field may favor simple vacancies while increasing the formation energy barriers for the interstitial or complex substitutional defects required for single-photon emission in the visible spectral region.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"4 Spectroscopic Characterization","content":"\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003e4.1 Emitter Classification and Spectral Inhomogeneity\u003c/h2\u003e \u003cp\u003eWe next analyze PL from LDW\u0026ndash;fabricated single emitters. Figure\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e gathers spectra of two spectral families distinguished from typical \u0026ldquo;green\u0026rdquo; emitters by their markedly more pronounced ZPL and suppressed PSB. We further classify these families by ZPL position and PSB strength: (i) red emitters at ~\u0026thinsp;650\u0026ndash;750 nm with a clearly weak PSB (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea), and (ii) emitters at ~\u0026thinsp;600\u0026ndash;650 nm with a stronger PSB (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb). The excitation polarization was fixed for all measurements, thus, relative intensities may reflect excitation-coupling biases for differently oriented dipoles. To compare vibronic structures, we normalize each spectrum and plot them on an energy axis aligned to the ZPL (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec\u0026ndash;d). Both families exhibit PSB features\u0026thinsp;~\u0026thinsp;170\u0026ndash;180 meV from the ZPL, consistent with coupling to high-DOS phonons of bulk hBN [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e, \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e]. Regarding microscopic origin, DFT studies have proposed carbon-related centers (e.g., carbon trimers or donor\u0026ndash;acceptor pairs) as candidates for weak-PSB red emitters, which is broadly consistent with our spectral data [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e, \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eDespite these shared characteristics, we observe a significant spread in ZPL positions and linewidths (FWHM) across the dataset. This spectral diversity is not unique to the LDW technique but is a pervasive feature of emitters in hBN, routinely reported across various fabrication methodologies including electron irradiation, ion implantation, and thermal annealing [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e, \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e, \u003cspan additionalcitationids=\"CR10 CR11\" citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. This universality suggests that the variations stem from local environmental heterogeneity, such as strain fields and charge fluctuations from nearby defects, as well as the potential existence of multiple stable atomic configurations (e.g., distinct defect-impurity complexes or defect pairs with varying separations) that possess slightly distinct transition energies [\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003e4.2 Excited state lifetimes\u003c/h2\u003e \u003cp\u003eExcited-state lifetimes for the red emitters were quantified using pulsed 520 nm laser excitation (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e). The decay traces are well-described by a single-exponential model, yielding lifetimes in the ~\u0026thinsp;2\u0026ndash;5 ns range (Table\u0026nbsp;1). The single-exponential behavior indicates a dominant radiative decay channel characteristic of high-purity quantum emitters. For completeness, Table\u0026nbsp;1 also lists ZPL positions and Lorentzian ZPL FWHM. We observe no systematic correlation between lifetime, ZPL position, and FWHM. The variation is likely attributable to differences in the local density of optical states (LDOS) arising from the emitters' varying depths within hBN or proximity to the dielectric interface, and not necessarily intrinsic to the defect structure itself [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003e4.3 Spectrally Resolved PL Polarization\u003c/h2\u003e \u003cp\u003eSpectrally resolved PL polarization measurements were performed on both emitter families. Representative datasets for strong-PSB (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ea\u0026ndash;c) and weak-PSB (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ed\u0026ndash;f) cases were acquired by stepping a linear polarizer in 10\u0026deg; increments. Both emitters exhibit clear linear-dipole behavior: the angle dependence of the ZPL-integrated intensity (background not subtracted; integration bands marked in Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ec,f) shows pronounced maxima/minima, and the PSB-integrated intensity follows the same angle without a visible offset. This alignment implies that the dominant phonon-assisted transitions inherit the ZPL dipole orientation, suggesting that electron\u0026ndash;phonon coupling does not strongly rotate the optical selection rules for these modes. Note, however, that 10\u0026deg; steps can miss the true extrema, so maxima/minima need not be separated by exactly 90\u0026deg;.\u003c/p\u003e \u003cp\u003eThe degree of linear polarization (DoLP) for hBN emitters seldom reaches unity, as noted in the literature [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e, \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e, \u003cspan additionalcitationids=\"CR32\" citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. Several mechanisms can contribute to incomplete DoLP: (i) intrinsic PL depolarization processes, (ii) an unpolarized background that sets a floor at the polarization minimum, or (iii) experimental geometry effects, such as the projection of an out-of-plane dipole component. Our spectrally resolved measurements identify broad, ZPL-free emission bands peaking in the 580\u0026ndash;600 nm window as the primary unpolarized background for the emitters in Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ec,f. This background is detrimental to QKD protocols because it degrades single-photon purity, likely explaining why \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{g}^{\\left(2\\right)}\\left(0\\right)\\)\u003c/span\u003e\u003c/span\u003e values saturate near ~\u0026thinsp;0.2 for some yellow\u0026ndash;red hBN emitters [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e, \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Consequently, identifying and suppressing these broad emission defects is critical for future fabrication optimization, or suggests a strategic shift to NIR emitters and using longer-wavelength excitation to avoid this spectral window. Additionally, the incomplete suppression of ZPL and PSB features at polarization minima (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ef) points to a combination of intrinsic depolarization and residual background (mechanisms i and ii) limiting the contrast.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003e4.4 Multiple-Emitter Sites and Spectral Statistics\u003c/h2\u003e \u003cp\u003eWe also examined fs-written sites that host multiple emitters per site, which are common at higher pulse energies, as well as sites that show broad, ZPL-free emission. The latter appears frequently peaking at ~\u0026thinsp;580\u0026ndash;600 nm (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003ea), shows no polarization dependence (cf. Figure\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e), and can persist as a broad background for polarized emitters, as discussed above. Exciting such sites with a 637 nm laser does not reproduce them next to the E₂g Stokes line (at ~\u0026thinsp;700 nm), ruling them out as, e.g., defect-based Raman modes. A plausible origin is PL from a defect center, or a family of defects, with a large Huang\u0026ndash;Rhys factor and strong electron\u0026ndash;phonon coupling. Vacancy-based emitters tend to have a large Huang\u0026ndash;Rhys factor, one such well known example is \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{V}_{\\text{B}}^{-}\\)\u003c/span\u003e\u003c/span\u003e in hBN, which exhibits a broad ZPL-free band near 800 nm [\u003cspan additionalcitationids=\"CR14\" citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eAt higher laser pulse energies, spectra frequently display multiple ZPLs within a single site (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eb), typically unresolved spatially in confocal scans (examples in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003ec\u0026ndash;d). Interestingly, some multiple-emitter sites show repeated exact ZPL energies (up to our spectral resolution), e.g., at ~\u0026thinsp;675 and ~\u0026thinsp;695 nm (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eb). Given the general ZPL spread discussed earlier, such recurrence may suggest identical microscopic centers under very similar local strain. We note that in our experimental study we did not employ energies that create exceedingly large craters in hBN, as some previous studies have done. In those higher-energy cases, the craters show pronounced material pileups and folded layers at crater edges, which tend to host a high density of defects and produce an almost continuous bright PL rim. Here, in contrast, for the multiple-emitter sites at larger energies the emission originates mostly from the central spot. This can be seen in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003ec\u0026ndash;d as larger hot-spot emission. Future studies will aim to carry out ODMR on LDW-created SPEs and such multiple-emitter hot spots, which may be interesting for magnetic sensing applications and to help deduce their atomic origin.\u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003ee compiles ZPL statistics across flakes, incorporating both single and multiple-emitter sites. Emitters occur most frequently in the 580\u0026ndash;650 nm range, where the average spectrum correspondingly exhibits the highest PL intensity. A secondary clustering is observed in the red band (650\u0026ndash;725 nm), accompanied by a distinct spectral feature in the average intensity around 675\u0026ndash;700 nm. Considering the typical spread of ZPLs in hBN, this secondary clustering likely indicates a specific \"red\" defect structure, or a family of defects with similar transition energies, that is preferentially generated by the LDW and annealing process.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"5 Conclusions","content":"\u003cp\u003eIn conclusion, we have demonstrated a workflow for femtosecond laser direct writing (LDW) of quantum emitters in hBN that facilitates the rapid post-process identification of written sites. Using an optical screening modality based on the intrinsic hBN E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e​ Raman Stokes line, we map laser-modified regions without relying \u003cem\u003ea priori\u003c/em\u003e on defect photoluminescence, adding a simple and broadly implementable tool to the LDW characterization toolkit. This approach reveals that the emergence of single-photon emitters coincides with a threshold regime of minimal lattice modification, whereas higher pulse energies more frequently yield multi-emitter spectra and broadband backgrounds. The generated emitters fall into two families: narrowband red emitters (650\u0026ndash;750 nm) with weak phonon sidebands and a 600\u0026ndash;650 nm family with stronger vibronic coupling, both exhibiting linear polarization and high single-photon purity. However, a broadband unpolarized background centered at 580\u0026ndash;600 nm can persist in some sites, degrading polarization contrast. Micro-Raman mapping reveals localized E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e​​ shifts consistent with compressive strain, yet these depth-averaged fields are insufficient to explain the ZPL spectral diversity, pointing to atomic-scale environmental heterogeneity as the dominant factor. Together, these results establish LDW combined with Raman-based screening as a practical characterization approach for qualifying hBN quantum emitter sites, paving the way further advancements on LDW-based integrated nanophotonic devices.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eSupplementary information\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSupporting Material 1.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe Research Council of Lithuania for their financial support.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFunded by the Research Council of Lithuania, Project No. P-ITP-24-22.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eData will be available on request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent to publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll authors read and approved the manuscript for publication.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eSugioka K, Cheng Y (2014) Ultrafast lasers\u0026mdash;reliable tools for advanced materials processing. Appl Phys Rev 1:041303. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1063/1.4904320\u003c/span\u003e\u003cspan address=\"10.1063/1.4904320\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGattass RR, Mazur E Femtosecond laser micromachining in transparent materials. Nat Photonics 2008, 2, 219\u0026ndash;225. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nphoton.2008.47\u003c/span\u003e\u003cspan address=\"10.1038/nphoton.2008.47\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDoherty MW, Manson NB, Delaney P, Jelezko F, Wrachtrup J, Hollenberg LC (2013) L. The nitrogen-vacancy colour centre in diamond. Phys Rep 528:1\u0026ndash;45. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.physrep.2013.02.001\u003c/span\u003e\u003cspan address=\"10.1016/j.physrep.2013.02.001\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAharonovich I, Tetienne J-P, Toth M (2016) Quantum emitters in hexagonal boron nitride. Phys Today 69(10):38\u0026ndash;44\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTran TT, Elbadawi C, Totonjian D, Lobo CJ, Grosso G, Moon H, Englund DR, Ford MJ, Aharonovich I, Toth M (2016) Robust multicolor single-photon emission from point defects in hexagonal boron nitride. ACS Nano 10(8):7331\u0026ndash;7338\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMendelson N et al (2017) Tunable and high-purity room temperature single-photon emission from atomic defects in hexagonal boron nitride. Nat Commun 8:15720. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-017-00810-2\u003c/span\u003e\u003cspan address=\"10.1038/s41467-017-00810-2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eScarani V, Bechmann-Pasquinucci H, Cerf NJ, Dušek M, L\u0026uuml;tkenhaus N, Peev M (2009) The security of practical quantum key distribution. Rev Mod Phys 81:1301\u0026ndash;1350. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/RevModPhys.81.1301\u003c/span\u003e\u003cspan address=\"10.1103/RevModPhys.81.1301\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHerrero-Collantes M, Garcia-Escartin JC (2017) Quantum random number generators. Rev Mod Phys 89:015004. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/RevModPhys.89.015004\u003c/span\u003e\u003cspan address=\"10.1103/RevModPhys.89.015004\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGan L, Zhang D, Zhang R, Zhang Q, Sun H, Li Y, Ning C-Z (2022) Large-scale, high-yield laser fabrication of bright and pure single-photon emitters at room temperature in hexagonal boron nitride. ACS Nano 16:14388\u0026ndash;14396. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acsnano.2c04386\u003c/span\u003e\u003cspan address=\"10.1021/acsnano.2c04386\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGao X, Li G, Li S et al (2021) Femtosecond laser writing of spin defects in hexagonal boron nitride. ACS Photonics 8:994\u0026ndash;1000. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acsphotonics.0c01847\u003c/span\u003e\u003cspan address=\"10.1021/acsphotonics.0c01847\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang Y-Z, Zhu T-X, Li Z-P et al (2023) Laser direct writing of visible spin defects in hexagonal boron nitride for applications in spin-based technologies. ACS Appl Nano Mater 6:14190\u0026ndash;14197. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acsanm.3c01047\u003c/span\u003e\u003cspan address=\"10.1021/acsanm.3c01047\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAhmad SI, Dave A, Sarpong E et al (2023) Dielectric breakdown and sub-wavelength patterning of monolayer hexagonal boron nitride using femtosecond pulses. 2D Mater 10:045029. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1088/2053-1583/acfa0f\u003c/span\u003e\u003cspan address=\"10.1088/2053-1583/acfa0f\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGottscholl A et al (2020) Initialization and readout of intrinsic spin defects in a van der Waals crystal at room temperature. Nat Mater 19:540\u0026ndash;545\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGottscholl A et al (2021) Room-temperature coherent control of spin defects in hexagonal boron nitride. Sci Adv 7:eabf3630\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHaykal A et al (2022) Decoherence of V\u003csub\u003eB\u003c/sub\u003e\u003csup\u003e-\u003c/sup\u003e spin defects in monoisotopic hexagonal boron nitride. Nat Commun 13:4347\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMackoit-Sinkevičienė M et al (2019) Carbon dimer defect as a source of the 4.1 eV luminescence in hBN. Appl Phys Lett 115:212101\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWeston L, Wickramaratne D, Mackoit M, Van de Alkauskas A (2018) Walle, C. G. Native point defects and impurities in hexagonal boron nitride. Phys Rev B 97:214104\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMaciaszek M, Razinkovas L, Alkauskas A (2023) Thermodynamics of carbon point defects in hexagonal boron nitride. Phys Rev B 108:014105\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAuburger P, Gali A (2023) Towards ab initio identification of paramagnetic substitutional carbon defects in hexagonal boron nitride acting as quantum bits. npj Comput Mater 9:38\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFischer M, Sajid A, Iles-Smith J et al (2023) Combining experiments with the polaron model and ab initio methods towards the microscopic origin of luminescent centres in hBN. Nanoscale 15:18638\u0026ndash;18652\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMohajerani SS, Chen S, Alaei A et al (2022) Narrowband quantum light emission from oxygen-related color centers in hexagonal boron nitride. ACS Photonics 9:2645\u0026ndash;2654\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCai Q, Scullion D, Han Y et al (2017) Raman signature and phonon dispersion of atomically thin hexagonal boron nitride. Phys Rev B 95:235407\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eStenger I, Nagler P, Neumann A et al (2017) Low-frequency Raman spectroscopy of few-atomic-layer hBN. 2D Mater 4:031003\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchaumburg F, Sleziona S, Z\u0026ouml;llner M et al (2023) Enhanced intensity of Raman signals from hexagonal boron nitride films. Beilstein J Nanotechnol 14:154\u0026ndash;165\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eReich S, Ferrari AC, Arenal R, Loiseau A, Bello I, Robertson J (2005) Resonant Raman scattering in cubic and hexagonal boron nitride. Phys Rev B 71:205201\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAndroulidakis C, Koukaras EN, Poss M, Papagelis K, Galiotis C, Tawfick S (2018) Strained hexagonal boron nitride: phonon shift and Gr\u0026uuml;neisen parameter. Phys. Rev. B 97, 241414(R)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBlundo E, Surrente A, Spirito D et al (2022) Vibrational properties in highly strained hexagonal boron nitride bubbles. Nano Lett 22:1150\u0026ndash;1157\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWhite S, Stewart C, Solntsev AS, Li C, Toth M, Kianinia M, Aharonovich I (2021) Phonon dephasing and spectral diffusion of quantum emitters in hexagonal boron nitride. Optica 8(9):1153\u0026ndash;1158\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFournier C, Watanabe K, Taniguchi T et al (2023) Investigating the fast spectral diffusion of a quantum emitter in hBN using resonant excitation and photon correlations. Phys Rev B 107:195304\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eClua-Provost T, Durand A, Frauni\u0026eacute; J et al (2022) Impact of thickness-dependent nanophotonic effects on the optical response of color centers in hexagonal boron nitride. ACS Photonics 9:3761\u0026ndash;3770\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJungwirth NR, Fuchs GD (2017) Optical absorption and emission mechanisms of single defects in hexagonal boron nitride. Phys Rev Lett 119:057401. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/PhysRevLett.119.057401\u003c/span\u003e\u003cspan address=\"10.1103/PhysRevLett.119.057401\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFournier C, Watanabe K, Taniguchi T et al (2021) Position-controlled quantum emitters with reproducible emission wavelength in hexagonal boron nitride. Nat Commun 12:3779. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-021-24019-6\u003c/span\u003e\u003cspan address=\"10.1038/s41467-021-24019-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKumar A, Peng P, Wrachtrup J et al (2024) Polarization dynamics of solid-state quantum emitters. ACS Nano 18:10745\u0026ndash;10768\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":true,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"discover-nano","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"narl","sideBox":"Learn more about [Discover Nano](https://www.springer.com/journal/11671)","snPcode":"11671","submissionUrl":"https://submission.nature.com/new-submission/11671/3","title":"Discover Nano","twitterHandle":"@SpringerOpen","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Discover Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"laser direct writing, single-photon emitters, micro-Raman, hBN, quantum technologies","lastPublishedDoi":"10.21203/rs.3.rs-8436272/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8436272/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eLaser direct writing (LDW) enables the spatially defined creation of room-temperature single-photon emitters (SPEs) in hexagonal boron nitride (hBN). However, the rapid characterization of written sites remains a bottleneck, and the available toolset for efficient screening is limited. Here, we demonstrate a streamlined LDW workflow utilizing single-shot pulses combined with a confocal screening technique that exploits the hBN E\u003csub\u003e2\u003cem\u003eg\u003c/em\u003e\u003c/sub\u003e​ Stokes line to rapidly localize and map laser-modified regions without relying \u003cem\u003ea priori\u003c/em\u003e on defect photoluminescence (PL). This approach enables the direct correlation of site morphology with PL hotspots, revealing that the emergence of single-photon emitters coincides with a threshold regime of minimal lattice modification. Micro-Raman spectral mapping further uncovers localized compressive strain surrounding these emission sites. We classify the generated defects into two families: narrowband \"red\" emitters (650\u0026ndash;750 nm) with weak phonon sidebands (PSB), and 600\u0026ndash;650 nm emitters with stronger vibronic coupling, both exhibiting linear polarization and high single-photon purity. These results establish a practical protocol for rapid prototyping, offering a valuable addition to the characterization toolkit for scalable quantum nanophotonics.\u003c/p\u003e","manuscriptTitle":"Laser Direct Writing and Raman Stokes-Contrast Screening of Quantum Emitter Sites in hBN","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2026-02-04 06:32:23","doi":"10.21203/rs.3.rs-8436272/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2026-01-30T12:07:51+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2026-01-28T17:20:26+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2026-01-25T15:38:11+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"108000281416711476001908638265576279374","date":"2026-01-21T04:08:37+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"73952796148115665202183617442425187296","date":"2026-01-19T15:25:49+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2026-01-17T07:42:25+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-12-25T03:50:04+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-12-25T03:48:58+00:00","index":"","fulltext":""},{"type":"submitted","content":"Discover Nano","date":"2025-12-23T17:44:45+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"discover-nano","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"narl","sideBox":"Learn more about [Discover Nano](https://www.springer.com/journal/11671)","snPcode":"11671","submissionUrl":"https://submission.nature.com/new-submission/11671/3","title":"Discover Nano","twitterHandle":"@SpringerOpen","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Discover Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"bc94c850-74b8-4a22-b3c0-d3acae429c6f","owner":[],"postedDate":"February 4th, 2026","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2026-03-30T16:16:53+00:00","versionOfRecord":{"articleIdentity":"rs-8436272","link":"https://doi.org/10.1186/s11671-026-04530-9","journal":{"identity":"discover-nano","isVorOnly":false,"title":"Discover Nano"},"publishedOn":"2026-03-26 16:08:51","publishedOnDateReadable":"March 26th, 2026"},"versionCreatedAt":"2026-02-04 06:32:23","video":"","vorDoi":"10.1186/s11671-026-04530-9","vorDoiUrl":"https://doi.org/10.1186/s11671-026-04530-9","workflowStages":[]},"version":"v1","identity":"rs-8436272","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8436272","identity":"rs-8436272","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}