Phonon-assisted emission and absorption of individual color centers in hexagonal boron nitride
Daniel Wigger Robert Schmidt Osvaldo Del Pozo-Zamudio Johann A. Preuß Philipp Tonndorf Robert Schneider Paul Steeger Johannes Kern Yashar Khodaei Jaroslaw Sperling Steffen Michaelis de Vasconcellos Rudolf Bratschitsch Tilmann Kuhn


Phonon coupling in hBN color centers]Phonon-assisted emission and absorption of individual
 color centers in hexagonal boron nitride





^1Institute of Solid State Theory, University of Münster, 48149 Münster, Germany
^2Institute of Physics and Center for Nanotechnology, University of Münster, 48149 Münster, Germany
^3Instituto de Investigación en Comunicación Optica, Universidad Autónoma de San Luis Potosí, 78210 San Luis Potosi, México
^4Terahertz Technology and Photonics, Hübner GmbH & Co K.G., 34123 Kassel, Germany

d.wigger@wwu.de




Defect centers in hexagonal boron nitride represent room-temperature single-photon sources in a layered van der Waals material. These light emitters appear with a wide range of transition energies ranging over the entire visible spectrum, which renders the identification of the underlying atomic structure challenging. In addition to their eminent properties as quantum light emitters, the coupling to phonons is remarkable. Their photoluminescence exhibits significant side band emission well separated from the zero phonon line (ZPL) and an asymmetric broadening of the ZPL itself. In this combined theoretical and experimental study we show that the phonon side bands can be well described in terms of the coupling to bulk longitudinal optical (LO) phonons. To describe the ZPL asymmetry we show that in addition to the coupling to longitudinal acoustic (LA) phonons also the coupling to local mode oscillations of the defect center with respect to the entire host crystal has to be considered. By studying the influence of the emitter's wave function dimensions on the phonon side bands we find reasonable values for size of the wave function and the deformation potentials. We perform photoluminescence excitation measurements to demonstrate that the excitation of the emitters is most efficient by LO-phonon assisted absorption.



[
    Daniel Wigger^1, Robert Schmidt^2, Osvaldo Del Pozo-Zamudio^2,3, Johann A. Preuß^2, Philipp Tonndorf^2, Robert Schneider^2, Paul Steeger^2, Johannes Kern^2, Yashar Khodaei^4, Jaroslaw Sperling^4, Steffen Michaelis de Vasconcellos^2, Rudolf Bratschitsch^2 and Tilmann Kuhn^1
    December 30, 2023
=====================================================================================================================================================================================================================================================================================





§ INTRODUCTION


Single-photon emitters are at the heart of many promising quantum technologies such as quantum computing and quantum cryptography. Although there are various solid-state emitters of single photons, there is still no system available, which simultaneously meets all requirements <cit.>. For example, self-assembled semiconductor quantum dots <cit.> represent a mature technology platform, but presently yield optimum performance only at cryogenic temperatures. Another class of single-photon emitters are defect centers in insulators, such as color centers in diamond <cit.>. These defect centers exhibit prominent single-photon emission also at room temperature but suffer from a high refractive index and challenging processability of the host crystal. Recently, a new class of single-photon emitters in atomically thin semiconductors has been discovered <cit.> and has been deterministically positioned on the nanoscale by strain engineering <cit.>. Defect centers in hexagonal boron-nitride (hBN) combine features from these classes <cit.>. They have the characteristics of an atomically sized defect center and at the same time share the advantages of a layered structure, i.e., the ultimate limit of miniaturization due to their atomic thickness and high mechanical robustness.

In this work, we investigate the photoluminescence (PL) and photon absorption of defect centers in hBN nanocrystals. By developing a theoretical model to calculate the PL spectrum taking into account the coupling to bulk longitudinal optical (LO) and longitudinal acoustic (LA) phonons and to the oscillation of a local mode, we explain measured emission spectra achieving an excellent agreement between experiment and theory. In addition, we verify by means of photoluminescence excitation (PLE) spectroscopy that coupling to LO phonons provides an efficient way of exciting individual single-photon emitters. This then allows for the isolation of a desired single-photon emitter from the wide range of emission energies of various defect centers in hBN.





§ RESULTS



 §.§ Emission spectra


Localized single-photon emitters are investigated in hBN nanopowder deposited on a Si/SiO_2 substrate (see Supplementary Material for details about the structure). Figure <ref> presents typical measured room-temperature PL spectra of six different localized light emitters in hBN. The transition energy of the emitters E_ ZPL, which is called the zero phonon line (ZPL) (labeled with numbers (1) to (6) in Fig. <ref>(a)) varies over a large range between 1.6 eV and 2.5 eV <cit.>. The energy of the excited states of the emitters has already been extensively studied, both theoretically via DFT calculations <cit.> and experimentally <cit.>. Although different types of atomic defects such as N- or B-vacancies or substitutions with carbon or oxygen might exist, the distribution of transition energies is remarkably homogeneous. Therefore, it is unlikely that each emission energy in hBN stems from a different type of defect. One promising proposal for the measured wide range of transition energies is that at least part of the energy spread is due to a Stark shift <cit.> resulting from an electric field which arises from a trapped charge near the color center <cit.>. Another suggested origin is a local strain distribution in the crystal <cit.>. Our results suggest that both effects might contribute to the broad energy spread of the ZPL energies.



The narrow ZPLs of the emitters are accompanied by characteristic phonon side bands at lower energies. To compare the phonon side bands of the different emitters more easily, the spectra are plotted against the detuning from the ZPL energy in Fig. <ref>(b). Two side bands well separated from the ZPL between 0.2 eV and 0.15 eV below E_ ZPL are visible, which can be attributed to two longitudinal optical (LO) phonons <cit.>. Also the ZPL itself is always more or less asymmetric with a steeper high energy side, which hints towards the contribution of low-energy (LE) phonon modes with energies below 50 meV. As will be shown later, we also assume two types of LE modes to be important. Comparing the phonon side bands of the different emitters in Fig. <ref>(b) we find significant variations. Not only the relative height of the LO phonon side bands compared to the ZPL intensity, but also the detailed structure, i.e., the relative heights of the LO_1 and LO_2 peak with respect to each other, vary from emitter to emitter. Having a closer look at the LE phonon modes near the ZPL, also this spectral phonon feature has different shapes and strengths for different emitters as highlighted in Fig. <ref>(b). In the following sections we will study these phonon side bands and develop a detailed understanding of their origin. In this context we will discuss which arguments support the identification of the LE phonons as longitudinal acoustic (LA) phonons and local mode oscillations of the defect with respect to the entire crystal.



 §.§ Phenomenological analysis

The properties of the LE and LO phonon side bands change significantly among the emitters. Therefore, we first try gaining a broad, quantitative overview of the different features. For this purpose we use a phenomenological model to fit the room temperature PL spectra of 165 different emitters. Figure <ref>(c) shows an exemplary spectrum (blue) and its phenomenological fit (orange). It is obvious that two LO phonon modes with energies at E_ LO,1≈ 165 meV and E_ LO,2≈ 195 meV have to be taken into account (red line). These energies are also supported by Raman measurements and ab initio calculations <cit.>. While trying different models for the LE phonons we found that at least two LE modes are required to achieve reasonable fits for all measured spectra. The best agreement was found for E_ LE_1 = 14 meV and E_ LE_2= 30 meV. We will later show that these energies can be well identified as representative for LA phonons in the case of LE_1 and for a local mode in the case of LE_2. But for now we will stick to the label LE phonons for the low energy side bands. The applied fit function for the PL intensity reads

    I(E) = I_ ZPL(E) + I_ LO(E) + I_ LE(E) ,

where we include a ZPL, single-phonon processes for the LO phonons and single- and two-phonon processes for the LE phonons. We also take single-phonon absorption processes for the LE modes into account. For all peaks we assume a Lorentzian shape. From each fit we retrieve the integrated peak weights A_ j (j=  ZPL, LO_1, LO_2, LE_1, LE_2) and respective peak widths γ_ k (k=  ZPL, LE_1, LE_2), where we assume that the LO side bands have the same broadening as the ZPL, which is a good approximation for dispersionless LO modes. Details about the fit function can be found in the Supplementary Material. In Fig. <ref>(c) we additionally show the two LE side bands (yellow for LE_1 and violet for LE_2) in phonon emission (solid) and absorption (dotted). The dashed turquoise line depicts the contribution from two phonon processes. Note that the absorption and two-phonon side bands are strictly related to the corresponding single-phonon emission side bands and therefore do not introduce additional fit parameters. The fit slightly deviates from the measurement for energies larger than the ZPL energy. The reason is, that a longpass filter is used close to the transition energy to block the exciting laser. This often reduces the measured high energy side of the ZPL. Therefore, we only used energies smaller than E_ ZPL for the fits.

For the large data set of 165 emitters we identify correlations between the different spectral features. We start by comparing the LO phonon side bands to the ZPL. For this purpose, Fig. <ref>(a) shows the LO contribution, i.e., the ratio A_ LO/A_ ZPL (A_ LO=A_ LO_1+A_ LO_2) between the LO side band and the ZPL weight, as a function of the ZPL energy. The colors of the dots represent the energy of the driving laser. The photon energies of the three lasers are marked by the dashed lines. We find transition energies almost equally distributed between 1.6 eV and 2.2 eV <cit.>. Even a few emitters appear above 2.4 eV. Larger transition energies could not be measured, because of the longpass filter blocking the exciting laser. Having a look at the values of the LO contribution, i.e., the vertical distribution of the dots, we find that they have a clear tendency following a positive slope. While for transition energies below 1.9 eV most of the points lie between 0 and 0.2, they appear between 0.1 and 0.45 for energies around 2.1 eV. This shows that the emitters with smaller transition energies couple on average less efficiently to the optical phonons. This is in good agreement with the assumption that the shift of the ZPL energy is associated with the Stark effect. A static electric field acting on the emitter's dipole should increase the distance between the different charge centers. In consequence this should also influence the coupling to the polar LO phonons. The additional spread of LO side band strengths for a fixed energy, especially around 2.1 eV, could be due to the effect of local strain. Distortions of the lattice constants in the vicinity of the defect center should not only influence the distance between the charge centers, but also change the size of the emitter's wave function and therefore the effective distance between the charge centers. Later on we will discuss the influence of these properties on the phonon coupling in more detail.



To shed more light on the properties of the LO phonons, Fig. <ref>(b) shows the energies of the LO phonon side bands as a function of the ZPL energy. We find two clusters of data points, which we associate to the two LO phonon modes near the Γ point. The energies of LO_1 (yellow) appear around 165 meV, those of LO_2 (violet) around 195 meV, which agrees well with calculated band structures and Raman measurements <cit.>. As indicated by the dotted lines, we find that the phonon energies have a slight trend to higher values for larger ZPL energies. As we identify these energies with the phonon energies of the hBN bulk crystal, this correlation suggests that the origin of the energy shift of the ZPL should also influence the LO phonon energies. This could be the case for a local strain field, which naturally changes the phonon dispersion. But also the dielectric environment might alter the LO phonon energies, as it affects in particular the LO-TO splitting <cit.>.

Figure <ref>(c) focusses on the LE phonon side bands by plotting the ratio between the first side band weight A_ LE_1 and the ZPL weight A_ ZPL as a function of the ZPL energy. The dots are scattered over a wide range from 0 to almost 1 for any transition energy. However, the majority of the points is equally distributed between 0.1 and 0.7. This finding does not provide any strong correlation between the coupling to the LE_1 (LA) phonons and the transition energy. The color and size of the dots represents the LE_1 contribution to the entire LE side band, i.e., A_ LE_1/A_ LE, where A_ LE=A_ LE_1+A_ LE_2 (small dots are blue, large dots red). Also this quantity shows no correlation in the dot pattern. Small blue and large red dots are found everywhere. This suggests that the coupling to LE phonons is not strongly affected by the parameters that are expected to govern the transition energy, i.e., the distance of the charge centers. However, we find a wide spread of dots on the vertical axis, i.e., a non-negligible variation of coupling strengths for a given transition energy. This finding will be traced back to variations of the emitter's size as discussed in more detail later.

Finally, Fig. <ref>(d) presents the ZPL width as a function of ZPL energy. We find a slight trend to smaller line widths at small transition energies, which is a similar to what has been found in Fig. <ref>(a) for the strength of the LO side bands. This is in line with established models which associate the ZPL width with second-order phonon coupling mechanisms <cit.>. The mechanism suggested in Ref. <cit.> relies on the scattering of acoustic phonons with the polaron formed by optical phonons. Consequently, a weaker coupling to LO phonons should result in a smaller lattice displacement of the polaron and therefore in a smaller scattering rate for the acoustic phonons. This could lead to the narrowing of the ZPL we found here. Another reason for this trend could be a stronger spectral wandering induced by trapped charges which are also considered as the origin of the Stark shift of the ZPL energies <cit.>.



 §.§ Theoretical model

The phenomenological model used to fit the 165 spectra in Sec. <ref> already includes some assumptions about the expected phonon modes, but it does not allow to retrieve further information about the emitters themselves. To support our assumptions about the mode energies we use a microscopically motivated model to reproduce the measured PL spectra. We also use this model to find possible origins for the correlations and non-correlations of the different phonon features and the ZPL energy found in Fig. <ref>. As already discussed in detail, our model should include the coupling to LO phonons to reproduce the two side bands around 165 meV and 195 meV below the ZPL. Also the coupling to acoustic phonons, which typically leads to a low-energy broadening of the ZPL, should be present in these systems. From defect centers in diamond, such as the extensively investigated NV^- center, it is well known that the prominent phonon side band is dominated by the coupling to local mode oscillations <cit.>. These modes describe displacements of the defect atoms with respect to the entire host crystal. Due to the similarity of the atomic defect structures, we expect that also in hBN local modes might have a significant influence on the PL spectrum.

We model the PL intensity spectrum of an emitter by calculating the optical susceptibility for a two-level system coupled to phonons, which  in the time domain is given by <cit.>


    χ(t) = iθ(t) exp[iω_0 t + Φ(t)-t/T_2] ,


where θ(t) denotes the Heaviside step function, ħω_0=E_ ZPL is the polaron-shifted transition energy of the emitter, i.e., the energy of the zero phonon line,  T_2 the dephasing time and Φ describes the phonon-induced dephasing with


    Φ(t) = ∑_ j∫_0^∞𝒥_ j(ω_ ph)/πω_ ph^2{(ħω_ ph/2k_ BT)[cos(ω_ ph t)-1]-isin(ω_ ph t)} dω_ ph .


This phonon influence is determined by the temperature T and the phonon spectral density 𝒥_ j at the phonon energy ħω_ ph for the phonon mode j. Note that such a model was commonly used to simulate the phonon coupling of F-centers <cit.>. In this section T=300 K is considered to allow comparison with the experimental results in Sec. <ref>. In the case of coupling to bulk phonons the spectral density is given by <cit.>


    𝒥_ j(ω_ ph) = ∑_ q|g_ j( q)|^2 δ[ω_ ph-ω_ j( q)] .


Here, g_ j( q) denotes the coupling matrix element for the coupling to a phonon with wave vector q in the phonon branch j and ω_ j( q) is the corresponding phonon dispersion relation. With Eq. eq:chi the absorption spectrum α(ω) of the emitter is simply given by the imaginary part of the Fourier transform


    α(ω) =   Im [ ∫_-∞^∞χ(t)  e^-iω t  dt]


From this expression the emission spectrum I(ω) is then retrieved by inverting the absorption α with respect to the ZPL at ω_0, i.e., 


    I(ω) = α(-ω + 2ω_0) .


We account for the coupling of the defect to LO_1 and LO_2 phonons via the Fröhlich-interaction


    g_ LO_j( q) = √(e^2ω_ LO_j( q)/2ε_0 V ħ[1/ε_∞( q)-1/ε_ s( q)]) 1/q [ F_ e( q) -F_ h( q)]


and for the coupling to LA phonons via the deformation potential coupling


    g_ LA( q) = q/√(2ϱħω_ LA( q)V) [ D_ eF_ e( q) - D_ hF_ h( q)] ,


where V is a normalization volume, ϱ is the mass density, ε_ s and ε_∞ are the static and high frequency dielectric functions, respectively, and D_ e and D_ h are the deformation potentials for electrons and holes, respectively. The form factors F_ e,h( q) are determined by the wave function of the excited state of the emitter as we explain in the following.

Another common coupling mechanism between charges and LA phonons is via the piezoelectric effect. In Ref. <cit.> it was shown that the piezoelectric constants vanish with an increasing number of hBN layers. As we are investigating nanopowder samples, we are dealing with small bulk crystals with large layer numbers. The thickness of the investigated hBN crystals is continuously distributed between 40 nm and 100 nm as exemplarily shown in the Supplementary Material. In Ref. <cit.> it was demonstrated that due to inversion symmetry the piezoelectric constants vanish for even layer numbers. If this coupling mechanism to LA phonons played an important role, it should only appear for half of the studied defects, i.e., those located in samples with odd layer numbers. This effect should lead to a separation of data points for crystals of even and odd layer number in Fig. <ref>(a) but one does not find such a feature. Therefore, we conclude that the piezoelectric coupling to LA phonons does not play an important role in our samples and assume that this coupling mechanism can be neglected.

One central aspect of hBN is its pronounced anisotropy arising from its layered structure. In our model we take this into account by distinguishing in-plane and out-of-plane directions. We assume that the dipole of the emitter lies in the plane of one hBN layer. Therefore, we consider a wave function of the defect's excited state consisting of two differently charged centers with an in-plane distance d as schematically shown in Fig. <ref>(a). The negative charge will be referred to as electron (e) and the positive one as hole (h). Each of the charge centers has a Gaussian wave function with an in-plane localization length a and an out-of-plane localization length b. This form of the wave function can be seen as an approximation for a wave function typically calculated from ab initio theory. The advantage of our approach is the flexibility to easily study the influence of changes of geometrical quantities. With this form of the wave function the form factors in Eqs. eq:LO and eq:LA read


    F_ e( q)    =   exp[-(q_r a/2)^2-(q_z b/2)^2 -  i/2 d· q]     and
    
    F_ h( q)    =   exp[-(q_r a/2)^2-(q_z b/2)^2 + i/2 d· q] ,


where the phonon wave vector q is split into an in-plane q_r and an out-of-plane q_z component and d is the vector connecting the two charge centers. After integrating |g_ j( q)|^2 in Eq. eq:J over the in-plane angle of the phonon wave vector, the influence of the wave function geometry on the coupling strengths is given by


    |g_ LO_j(q_r,q_z)|^2 ∼exp{-1/2[(q_r a)^2 + (q_z b)^2]}[1-J_0(dq_r)]

and

    |g_ LA(q_r,q_z)|^2 ∼exp{-1/2[(q_r a)^2 + (q_z b)^2]}[D_ e^2+D_ h^2-2D_ eD_ hJ_0(dq_r)] ,


where J_0 is the Bessel function of first kind and of zeroth order. Due to these form factors, coupling to phonons is essentially restricted to the wave vector region q_r ≲ a^-1, q_z ≲ b^-1.



For the Fröhlich coupling, different dielectric constants for the in-plane and the out-of-plane direction are considered. The energies of the phonons, which enter the model, are extracted from the simulated dispersion relation shown in Fig. <ref>(b) adapted from Ref. <cit.>; we assume two LO modes (LO_1 and LO_2) with constant dispersions. For the LA phonons we consider linear dispersions, where we again distinguish between the in-plane and the out-of-plane direction. More details about the model parameters can be found in the Supplementary Material. As will be shown below, the 14 meV LE_1 mode introduced in Sec. <ref> can be seen as a reasonable approximation for the side band stemming from the LA phonons. In Ref. <cit.> it was calculated that the C_ BN_ V center has the strongest coupling to an in-plane breathing mode of the defect center, which has an energy of 30 meV. Also other studies suggest this atomic structure as promising candidate for the hBN color center <cit.>. This motivates us to interpret the LE_2 mode in the phenomenological model in Sec. <ref> as representative of such a breathing mode. The coupling to this local mode is modeled by a Lorentzian-like distribution of the phonon spectral density <cit.>
 
    𝒥_ LOC(ω_ ph) = 1/π16 g^2 ω_ LOCΔω_ ph/(ω_ ph^2-ω_ LOC^2)^2+4 Δ^2ω_ ph^2 .

The parameters entering this coupling are the local mode energy ħω_ LOC=E_ LOC, the coupling strength g and the width Δ. For the considered parameters the spectral shape of the side band due to the local mode agrees well with the ordinary Lorentzian shape considered for the phenomenological fits in Sec. <ref>. We show in the Supplementary Material that our model for the local mode coupling reproduces PL spectra of different defect centers in diamond very well, showing the accuracy of the model. To reproduce the finite width of the ZPL, we assume a constant dephasing rate 1/T_2, providing us with a dephasing time T_2. We note that this dephasing time includes pure dephasing, homogeneous broadening mechanisms as well as inhomogeneous processes, e.g., spectral wandering <cit.>.

While the phonon dispersion relations and dielectric constants are well known for hBN, the deformation potentials are only hardly known. Simulations indicate that the relevant values are in the range of some hundred meV to a few eV <cit.>. We therefore will try finding reasonable values for the deformation potentials D_ e and D_ h for the LA phonon coupling of the considered defect centers in hBN. The other unknown parameters in our model are the in-plane and out-of-plane localization lengths a and b, respectively, and the charge separation d in the defect state. To be able to extract information on these parameters from the measured spectra, we will analyze their role for the coupling to LO and LA phonon modes, respectively.



For the coupling to LO phonons the only unknown quantities are the geometrical parameters in-plane localization length a, out-of-plane localization length b, and charge distance d. To investigate their influence, we first concentrate on the coupling to LO phonons. In Fig. <ref>(a) we plot the contribution of the LO single-phonon side band, i.e., the fraction between the weight of the first LO phonon side bands of LO_1 and LO_2 and the ZPL weight, as a function of the charge distance d for different localization lengths a between 0.1 nm (blue) and 1 nm (red). The distance d is varied on the same scale. This should be a reasonable range for these quantities as the in-plane lattice constant of hBN is 0.4 nm and the extension of an atomistic defect state should be of similar size. Under the assumption that the Stark effect is responsible for the different transition energies of the emitters, we additionally plot the measured data points from <ref>(a) as gray dots into Figs. <ref>(a). The relation between the energies of the experimental values of Fig. <ref>(a) (gray dots in Fig. <ref>(a)) and the charge distances d in Fig. <ref>(a) is obtained by the fits of the measured spectra in Sec. <ref>, which yield values between d = 0.2 nm and d=0.5 nm. By the given choice of the ZPL energy axis the majority of the measurement points appear in this range. The simulated LO contributions vanish for vanishing in-plane distance d because in this limit there is no dipole moment associated with the defect. For a given in-plane localization length a the LO contribution grows with increasing d. When increasing the localization length a the curves rise slower, because the Fröhlich coupling strength depends on the effective dipole strength of the wave function. To achieve the same strength for larger localization lengths a the charge separation has to be accordingly larger. The increase of the LO contribution with increasing distance nicely resembles the distribution of the measured dot pattern, which supports the interpretation that the shift of the ZPL is at least partially related to a Stark shift. To further analyze the influence of the LO phonons, we show exemplary normalized PL spectra for different values of a and d in Fig. <ref>(b). The colors agree with the corresponding curves in (a). At the bottom for a=0.5 nm the two single-phonon LO side bands are well resolved around -0.2 eV, also the three two-phonon side bands between -0.3 eV and -0.4 eV can be seen. As expected the amplitude of the side bands increases for larger in-plane distance d. At the top of Fig. <ref>(b) the same spectra are shown for a=0.1 nm. Here, the LO side bands grow very rapidly when increasing d, they even become larger than the ZPL. For the largest considered d the two-phonon side bands are even larger than the single-phonon side bands. Indeed, from Eqs. eq:Phi and eq:J it can be shown that, as long as successive side bands do not spectrally overlap, the contributions of single-phonon (A_ LO) and two-phonon (A_ 2LO) side bands with respect to the ZPL are given by

    A_ LO/A_ ZPL = ∑_ j, q| g_ LO_j( q)/ω_ LO_j|^2  ,     A_ 2LO/A_ ZPL =1/2( A_ LO/A_ ZPL)^2  .

Thus, as soon as A_ LO>2A_ ZPL the two-phonon sideband will exceed the single-phonon one. However, in the measured PL spectra of hBN color centers this regime where the side bands exceed the ZPL is not reached. Therefore, we can conclude that the localization length should lie between a=0.3 nm and 1 nm and the charge distance between d=0 nm and 1 nm. 

After finding reasonable values for the in-plane size of the wave function (a and d), in Figs. <ref>(c) we focus on the values for the deformation potentials D_ e and D_ h. When considering the special case of a vanishing charge distance, d=0, the deformation potential coupling strength in Eq. eq:gLA is proportional to (D_ e-D_ h)^2. So we will use this limit to find reasonable values for the difference of the deformation potentials. To provide a quantitative measure, we calculate the LA contribution to the entire spectrum via

    1- Im[ χ_ LA(t→∞; ω_0=0; T_2→∞) ] .

This value is plotted as a function of |D_ e-D_ h| in the inset in Fig. <ref>(c) for wave function dimensions of a=0.4 nm, b=0.2 nm and d=0 nm. For D_ e=D_ h the LA contribution vanishes and no LA side band appears in the spectrum, while the side band dominates the entire spectrum for large values of |D_ e-D_ h|. In Fig. <ref>(c) we found that the LA contribution in the measured spectra also for small transition energies, i.e., small d, yields values between zero and one. Therefore, reasonable values for the difference of the deformation potentials should lead to an LA contribution of around 0.5. We find that this is fulfilled for |D_ e-D_ h|≈ 0.6 eV. However, the deformation potentials not only determine the LA contribution but are also responsible for the shape of the spectra. To demonstrate this influence, we show exemplary spectra for different values of |D_ e-D_ h| still for a vanishing dipole moment d=0 in Fig. <ref>(c). We find that with increasing |D_ e-D_ h| (blue to red) the asymmetry of the ZPL peak increases till the entire spectrum is dominated by the LA side band for the 1 eV case (red). Here the ZPL is only a rather small peak on the broad phonon background. From here on we fix |D_ e-D_ h|=0.6 eV.

As already mentioned, we assume that the coupling to the LE_1 mode phenomenologically describes the coupling to LA phonons. One might wonder why the coupling to a discrete mode with a Lorentzian lineshape approximates reasonably well the coupling to the continuum of LA phonons. To confirm that this is indeed a valid approximation, the black solid line in Fig. <ref>(d) shows a simulated room-temperature PL spectrum, where the coupling to the continuum of LA phonons is fully taken into account. This spectrum was then fitted by the phenomenological model of Sec. <ref> yielding the orange line. The contributions of the two LE modes (yellow and violet) and the ZPL (green) in the phenomenological model are plotted separately as before. We find an excellent agreement between the simulated spectrum including only LA phonons and the phenomenological fit. Furthermore we find that only the LE_1 mode has a relevant contribution in the fit. This shows that indeed the LE_1 mode is a reasonable approximation for the contribution of the LA phonons to the spectrum.

In the next step we investigate the influence of the sum of the deformation potentials. For a non-vanishing charge distance d the Bessel function in Eq. eq:gLA does not vanish leading to the fact that the phonon coupling depends on D_ e and D_ h separately. In Fig. <ref>(a) we show simulated PL spectra for a=0.4 nm, b=0.2 nm and d=0.5 nm and different values of D_ e+D_ h between 0 eV (blue) and 1 eV (red). We find that the weight of the LA side bands close to the ZPL strongly increases with the sum of the deformation potentials. It is also nicely visible that the LO phonon side bands are additionally broadened by the LA phonons, demonstrating that multi-phonon contributions are well reproduced by the model. In Fig. <ref>(c) we found that the LA contribution (LE_1 mode) to the spectrum exhibits no significant correlation with the transition energy, i.e., with d. This finding should be reproduced by the choice of D_ e+D_ h. Therefore, in Fig. <ref>(b) we plot the LA contribution as a function of the charge distance d for different values of D_ e+D_ h between 0 eV (blue) and 4 eV(red) (a=0.4 nm and b=0.2 nm). All curves start at the same value of slightly below 0.6 for d=0 nm. We find that for deformation potential sums of 1 eV or larger the LA contribution grows significantly with increasing d. One would therefore expect a clear correlation between these two quantities. This makes us conclude that D_ e+D_ h should be smaller than 1 eV. In the special case of D_ e+D_ h=0.6  eV=|D_ e-D_ h| the LA contribution is a straight line, because in this case either D_ h=0 eV or D_ e=0 eV, which makes the charge distance irrelevant. However, we do not expect one of the  deformation potentials, either for electrons or for holes, to vanish in hBN. This would mean that either electrons or holes do not scatter with LA phonons via the deformation potential coupling, which makes us disregard this value. We found good results for D_ e+D_ h=0.2 eV (D_ e=0.4 eV, D_ h=-0.2 eV) as will be demonstrated below in Sec. <ref> when directly comparing the simulations with measured PL spectra.



After having selected reasonable values for the deformation potentials, we now study the influence of the other two geometrical parameters, the in-plane localization length a and the out-of-plane localization length b of the wave function. Figures <ref>(c,d) focus on the in-plane localization length a. In Fig. <ref>(c) exemplary spectra for d=0.5 nm, b=0.2 nm and different values of a between 0.2 nm (blue) and 1 nm (red) are shown. While the LO phonon side bands increase significantly when reducing a, as discussed earlier, the LA phonon side band next to the ZPL does not change significantly. To provide a quantitative measure for the influence of a on the LA coupling in Fig. <ref>(d) we plot the LA contribution as a function of a for different values of d. All curves slightly decrease from values between 0.6 and 0.7 to 0.4 for a=1 nm. This spread can give rise to part of the spread of LA contributions found in Fig. <ref>(c).

Finally, the influence of the out-of-plane localization length b is considered. In Fig. <ref>(e) we again show PL spectra for a=0.4 nm, d=0.5 nm and different values of b between 0.2 nm (blue) and 1 nm (red). Two pronounced features are found in the spectra: (i) the LA phonon side band increases with decreasing b; (ii) the two LO single-phonon side bands LO_1 and LO_2 slightly increase with decreasing b, but not by the same rate. While LO_1 and LO_2 are almost equally strong for b=1 nm (red), LO_1 is significantly stronger for 0.2 nm (blue). These two features are quantified in Fig. <ref>(f). The solid lines show the LA contribution and the dashed lines the ratio between the peak intensities of the two LO single-phonon side bands I_ LO_2/I_ LO_1 as functions of b for different values of a. The LA contribution rapidly drops from values up to 0.9 for b=0 nm to 0.2 for b=1 nm. This spread of LA phonon side band contributions agrees well with the measured values found in Fig. <ref>(c). The fraction of the LO side band peaks increases for growing b. This variation of the shape of the LO side band was also found in the measured spectra in Fig. <ref>. The reason for this variation of the LO side bands is that the larger LO phonon energy is not present for the out-of-plane phonon wave vectors, as can be seen in Fig. <ref>(b). So b mainly influences the LO_1 phonon side band. A larger out-of-plane localization length b in Eq. eq:gLO leads to a smaller range of out-of-plane phonon wave vectors that is excited by the emitter, which leads to a weaker effective coupling strength. Therefore the fraction I_ LO_2/I_ LO_1 reduces for growing b. The out-of-plane lattice constant of hBN is 0.66 nm, which makes the considered values for b a reasonable range for this parameter.

This parameter study shows that slight variations of the geometrical parameters of the defect's wave function in the range of reasonable dimensions lead to strong changes in the PL spectra. In the next step, we apply this model to directly reproduce measured PL spectra by varying the geometrical parameters. We also take the local mode for the low-energy coupling in the range of 30 meV into account and vary its coupling parameters.



 §.§ Comparison between experiment and theory




Figure <ref> shows PL measurements of different defect centers. While the spectra in Figs. <ref>(a,b) are taken at room temperature, Figs. <ref>(c,d) are measured at cryogenic temperatures of about 7 K. Figure <ref>(e) depicts a temperature series taken on the same emitter. The temperature decreases from top to bottom from 300 K to 20 K. In addition to the experiments (colored lines) we show corresponding simulations as black lines. All parameters for the simulations are listed in Tab. <ref>.



Comparing the different spectra, we find a convincing agreement between the measured and the calculated spectra. The theoretical approach used in Ref. <cit.>, which considers the phonon spectral function as a free parameter, also led to good agreements with measured spectra. However, the authors did not draw conclusions on the emitters' wave functions and details of the electron-phonon coupling matrix elements. From the parameters for the dimensions of the wave functions in our model listed in Tab. <ref> we find that all sizes (a, b, and d) are in the range between 0.1 nm and 0.5 nm. This is a reasonable parameter range, because it shows that the assumed emitter wave function spreads over a few unit cells at most and is mainly restricted to the layer in which the defect is located. While the LO phonon energies scatter slightly around the already discussed values of E_ LO_1≈ 165 meV and E_ LO_2≈ 195 meV, the energy of the local mode is rather fixed. This might have two reasons. On the one hand, the side band in the PL spectrum partly overlaps with the LA contribution at room temperature, which makes a precise determination of the peak position less easy. On the other hand, strain distributions could have less influence on the energetic position of the mode, but rather change the coupling strength g and the width Δ. We want to remark that, the low temperature spectrum in Fig. <ref>(d) shows a pronounced peak at E≈ 2.05 eV, i.e., 35 meV below the ZPL. The dashed line includes this peak as a local mode in the simulation. However, this spectral feature could also stem from another nearby emitter. Usually, PL spectra at cryogenic temperatures exhibit a large number of narrow lines in these samples, which makes a clear identification of the lines difficult.

To reproduce the entire temperature series in Fig. <ref>(e) we fix the dimensions of the wave function and the local mode coupling strength for T=300 K as given in Tab. <ref>. For each temperature we have to additionally determine the local mode width Δ and the dephasing time T_2. We find that both the ZPL and the local mode side band, which remains as narrow peak slightly above E=2.1 eV for T=20 K, become narrower with lower temperatures. This can already be seen when comparing the room-temperature spectra in Figs. <ref>(a,b) to the low-T ones in Figs. <ref>(c,d). In the Supplementary Material we discuss another temperature series of PL spectra taken on a different emitter, where we additionally measured the PL lifetime. In agreement with Ref. <cit.> we find that the life time does not depend on the temperature and that it is in the range of a few nanoseconds <cit.>. This indicates that the ZPL line width at room temperature is dominated by the pure dephasing mechanism and spectral wandering, which we combine in the fit parameter of the dephasing time T_2. For decreasing temperatures, the ZPL width shrinks significantly as also found in Refs. <cit.>. However, the values for low temperatures between 5 K and 50 K are underestimated because the measured line width is on the order of the limited resolution of the spectrometer (100 µeV). Hence, we are not able to draw reliable conclusions about the low-temperature T_2 times from our data. The evolution of the two low-energy phonon side bands with decreasing temperature nicely supports our assumptions of their nature. The local mode results in an isolated peak at the same energy. This is not expected for phonons with continuous spectrum, i.e., for LA modes. The LA phonons lead to a low-energy side band directly attached to the ZPL, as it is well known from semiconductor QDs <cit.>. In the spectra in Fig. <ref>(e) we find that the spectral shape of the LO_1 phonon side band is additionally broadened to the high energy side compared to the simulation. When moving from high to low T, all spectral features become narrower, which makes it obvious that this broadening evolves into another maximum just below E=2 eV. The fits for the low-temperature spectra in Figs. <ref>(c,d) show that our model reproduces the LO phonon side bands, despite our rough approximations for the phonon dispersion when assuming constant LO phonon energies. However, it was shown that strain strongly influences the emitters' spectra <cit.> and that it is highly likely that the emitters appear near the surface of the sample <cit.>. Local strain or the proximity to the surface will alter the phonon band structure and could thereby lead to variations of the spectral shape of the LO side bands. This could be a reason for the atypical broadening of the LO_1 side band for this emitter.




 §.§ Photoluminescence excitation spectra


We now focus again on the distribution of ZPL energies in Fig. <ref>(a). We find that many of the green dots cluster around E_ ZPL = 2.15 eV, while the red dots cluster around 1.8 eV. These energies are approximately 150–200 meV below the respective photon energy of the laser used for excitation (dashed lines), which is the range of LO phonon energies. This suggests that LO-phonon assisted absorption could efficiently drive the emitters. Therefore, we plot a histogram of detunings between the exciting laser energy E_ L and the ZPL energy E_ ZPL in Fig. <ref>(a), where the colors correspond to the colors in Fig. <ref>(a). In total, nearly 50% of all 165 investigated emitters appear between 150 meV and 200 meV below the laser energy. The remaining emitters are almost equally distributed over the remaining 800 meV. This agrees with the findings in Ref. <cit.>.



To confirm that the emitters are efficiently excited by LO-phonon assisted absorption, we have performed PLE measurements at room temperature to obtain information about the absorption properties of the emitters at different photon energies. The results are shown in Fig. <ref>(b). We chose an emitter which has a strong LO phonon sideband in the PL spectrum, as shown with the orange line. At the high-energy side of the ZPL at 2.03 eV the spectrum is again cut off by a longpass filter. We measure PL spectra for a wide range of different excitation energies from 2.09 eV to 2.3 eV as shown in the Supplementary Material. After fitting and subtracting the background and Raman lines as described in the Supplementary Material, we extract the PLE data by integrating the PL spectra either over the ZPL, i.e., from 1.95 eV to 2.05 eV (blue shaded area), or over the LO phonon sideband, i.e., from 1.75 eV to 1.9 eV (yellow shaded area). The PLE data for the ZPL and for the LO phonon sideband are shown in blue and yellow, respectively. We find that both PLE data sets exhibit a pronounced maximum between 2.2 eV and 2.3 eV. This energy range nicely agrees with the phonon sidebands of the mirrored PL spectrum (green line). This observation demonstrates that the excitation of the emitter is very efficient via the LO-phonon assisted process.

Our study shows that the LO-phonon assisted excitation of the excited defect state provides an optimized way to drive the system. Therefore we can conclude that it is conceivable that we select a specific ZPL energy ħω_0, when pumping the system with energies between ħω_0+160 meV and ħω_0+200 meV, i.e, via a LO-phonon assisted absorption. This provides a strategy of isolating single-photon emitters with a desired emission energy from the wide range of possible transitions.



§ CONCLUSIONS


We have shown that the coupling to different phonon modes plays a crucial role for localized light emitters in hBN. On the one hand, coupling to bulk phonons leads to an asymmetric broadening of the ZPL in the case of LA modes and to the appearance of prominent sidebands well separated from the ZPL in the case of LO modes. On the other hand the coupling to a local mode oscillation contributes to the broadening of the ZPL at room temperature, while the respective side band is discernible at cryogenic temperatures. By fitting measured PL spectra with our theoretical model we were able to extract parameters for the wave function geometries of the emitters. We have shown that the distance between the positive and the negative charge center may be connected to the energy of the ZPL, supporting the assignment of the ZPL energy spread with the Stark effect caused by nearby charges. It was also possible to find reasonable values for the deformation potential of electrons and holes in hBN. By measuring PLE spectra, we have demonstrated that these sidebands are also present in absorption and lead to an efficient absorption via LO-phonon assisted transitions. Finally we have shown that it is possible to preferentially select emitters with a desired transition from the wide range of emission energies possible in hBN by aiming at the LO-phonon assisted absorption. Our results lead to a deeper understanding of the fundamental properties of color centers in hBN and pave the way to a tailored control of the excitation process.





We gratefully acknowledge support by the research group FOR 1493 of the Deutsche Forschungsgemeinschaft. We thank Michael Rohlfing for fruitful discussions.



§ REFERENCES


10
url<#>1#1urlprefixURL 



aharonovich2016sol
Aharonovich I, Englund D and Toth M 2016 Nature Photon. 10
  631–641

santori2001tri
Santori C, Pelton M, Solomon G, Dale Y and Yamamoto Y 2001 Phys. Rev.  Lett. 86 1502

michler2003sin
edited by Peter Michler 2003 Single quantum dots: Fundamentals,
  applications and new concepts (Topics in Applied Physics) vol 90 (Springer
  Verlag, Berlin)

zaitsev2000vib
Zaitsev A M 2000 Phys. Rev. B 61 12909

jelezko2006sin
Jelezko F and Wrachtrup J 2006 Phys. Status Solidi A 203
  3207–3225

beha2012opt
Beha K, Batalov A, Manson N B, Bratschitsch R and Leitenstorfer A 2012 Phys. Rev. Lett. 109 097404

tonndorf2015sin
Tonndorf P, Schmidt R, Schneider R, Kern J, Buscema M, Steele G A,
  Castellanos-Gomez A, van der Zant H S J, Michaelis de Vasconcellos S and
  Bratschitsch R 2015 Optica 2 347–352

koperski2015sin
Koperski M, Nogajewski K, Arora A, Cherkez V, Mallet P, Veuillen J Y, Marcus J,
  Kossacki P and Potemski M 2015 Nature Nanotech. 10 503–506

srivastava2015opt
Srivastava A, Sidler M, Allain A V, Lembke D S, Kis A and Imamoğlu A 2015
  Nature Nanotech. 10 491–496

he2015sin
He Y M, Clark G, Schaibley J R, He Y, Chen M C, Wei Y J, Ding X, Zhang Q, Yao
  W, Xu X, Lu C Y and Pan J W 2015 Nature Nanotech. 10 497–502

chakraborty2015vol
Chakraborty C, Kinnischtzke L, Goodfellow K M, Beams R and Vamivakas A N 2015
  Nature Nanotech. 10 507–511

tonndorf2017sin
Tonndorf P, Schwarz S, Kern J, Niehues I, Del Pozo-Zamudio O, Dmitriev A I,
  Bakhtinov A P, Borisenko D N, Kolesnikov N N, Tartakovskii A I, Michaelis de
  Vasconcellos S and Bratschitsch R 2017 2D Materials 4 021010

palacios2016ato
Palacios-Berraquero C, Barbone M, Kara D M, Chen X, Goykhman I, Yoon D, Ott
  A K, Beitner J, Watanabe K, Taniguchi T, Ferrari A C and Atatüre M 2016
  Nature Comm. 7

kern2016nan
Kern J, Niehues I, Tonndorf P, Schmidt R, Wigger D, Schneider R, Stiehm T,
  Michaelis de Vasconcellos S, Reiter D E, Kuhn T and Bratschitsch R 2016 Adv. Mater. 28 7101–7105

branny2017det
Branny A, Kumar S, Proux R and Gerardot B D 2017 Nature Comm. 8
  15053

palacios2017lar
Palacios-Berraquero C, Kara D M, Montblanch A R P, Barbone M, Latawiec P, Yoon
  D, Ott A K, Loncar M, Ferrari A C and Atatüre M 2017 Nature Comm.
  8 15093

tran2016quaI
Tran T T, Bray K, Ford M J, Toth M and Aharonovich I 2016 Nature
  Nanotech. 11 37–41

tran2016quaII
Tran T T, Zachreson C, Berhane A M, Bray K, Sandstrom R G, Li L H, Taniguchi T,
  Watanabe K, Aharonovich I and Toth M 2016 Phys. Rev. Appl. 5
  034005

tran2016rob
Tran T T, Elbadawi C, Totonjian D, Lobo C J, Grosso G, Moon H, Englund D R,
  Ford M J, Aharonovich I and Toth M 2016 ACS Nano 10 7331–7338

vuong2016pho
Vuong T Q P, Cassabois G, Valvin P, Ouerghi A, Chassagneux Y, Voisin C and Gil
  B 2016 Phys. Rev. Lett. 117 097402

martinez2016eff
Martínez L J, Pelini T, Waselowski V, Maze J R, Gil B, Cassabois G and
  Jacques V 2016 Phys. Rev. B 94 121405

jungwirth2016temp
Jungwirth N R, Calderon B, Ji Y, Spencer M G, Flatté M E and Fuchs G D 2016
  Nano Lett. 16 6052–6057

shotan2016pho
Shotan Z, Jayakumar H, Considine C R, Mackoit M, Fedder H, Wrachtrup J,
  Alkauskas A, Doherty M W, Menon V M and Meriles C A 2016 ACS
  Photonics 3 2490–2496

exarhos2017opt
Exarhos A L, Hopper D A, Grote R R, Alkauskas A and Bassett L C 2017 ACS
  nano 11 3328–3336

li2017non
Li X, Shepard G D, Cupo A, Camporeale N, Shayan K, Luo Y, Meunier V and Strauf
  S 2017 ACS Nano 11 6652–6660

jungwirth2017opt
Jungwirth N R and Fuchs G D 2017 Phys. Rev. Lett. 119 057401

sontheimer2017pho
Sontheimer B, Braun M, Nikolay N, Sadzak N, Aharonovich I and Benson O 2017
  Phys. Rev. B 96 121202

grosso2017tun
Grosso G, Moon H, Lienhard B, Ali S, Efetov D K, Furchi M M, Jarillo-Herrero P,
  Ford M J, Aharonovich I and Englund D 2017 Nat. commun. 8 705

tran2017res
Tran T T, Kianinia M, Nguyen M, Kim S, Xu Z Q, Kubanek A, Toth M and
  Aharonovich I 2017 ACS Photonics 5 295–300

chejanovsky2017qua
Chejanovsky N, Kim Y, Zappe A, Stuhlhofer B, Taniguchi T, Watanabe K, Dasari D,
  Finkler A, Smet J H and Wrachtrup J 2017 Sci. Rep. 7 14758

noh2018sta
Noh G, Choi D, Kim J H, Im D G, Kim Y H, Seo H and Lee J 2018 Nano
  Lett. 18 4710–4715

xu2018sin
Xu Z Q, Elbadawi C, Tran T T, Kianinia M, Li X, Liu D, Hoffman T B, Nguyen M,
  Kim S, Edgar J H, Wu X, Song L, Ali S, Ford M, Toth M and Aharonovich I 2018
  Nanoscale 10 7957–7965

grinyaev2002dee
Grinyaev S N, Konusov F V and Lopatin V V 2002 Phys. Solid State 44 286–293

tawfik2017fir
Tawfik S A, Ali S, Fronzi M, Kianinia M, Tran T T, Stampfl C, Aharonovich I,
  Toth M and Ford M J 2017 Nanoscale 9 13575–13582

abdi2018col
Abdi M, Chou J P, Gali A and Plenio M B 2018 ACS Nano 5
  1967–1976

proscia2018nea
Proscia N V, Shotan Z, Jayakumar H, Reddy P, Cohen C, Dollar M, Alkauskas A,
  Doherty M, Meriles C A and Menon V M 2018 Optica 5 1128–1134

bourrellier2014nan
Bourrellier R, Amato M, Galvão Tizei L H, Giorgetti C, Gloter
  A, Heggie M I, March K, Stéphan O, Reining L, Kociak M and Zobelli A 2014
  ACS Photonics 1 857–862

serrano2007vib
Serrano J, Bosak A, Arenal R, Krisch M, Watanabe K, Taniguchi T, Kanda H, Rubio
  A and Wirtz L 2007 Phys. Rev. Lett. 98 095503

sohier2017bre
Sohier T, Gibertini M, Calandra M, Mauri F and Marzari N 2017 Nano
  Lett. 17 3758–3763

machnikowski2006cha
Machnikowski P 2006 Phys. Rev. Lett. 96 140405

huxter2013vib
Huxter V M, Oliver T A A, Budker D and Fleming G R 2013 Nature Phys.
  9 744–749

kehayias2013inf
Kehayias P, Doherty M W, English D, Fischer R, Jarmola A, Jensen K, Leefer N,
  Hemmer P, Manson N B and Budker D 2013 Phys. Rev. B 88
  165202

duke1965pho
Duke C B and Mahan G D 1965 Phys. Rev. 139 A1965

mahan2013man
Mahan G D 2013 Many-particle physics (Springer Science & Business
  Media)

schmitt1987the
Schmitt-Rink S, Miller D A B and Chemla D S 1987 Phys. Rev. B
  35 8113

krummheuer2002the
Krummheuer B, Axt V M and Kuhn T 2002 Phys. Rev. B 65
  195313

reiter2014the
Reiter D E, Kuhn T, Gl'́assl M and Axt V M 2014 J. Phys. Condens.  Matter 26 423203

michel2011pho
Michel K H and Verberck B 2011 Phys. Rev. B 83 115328

sajid2018def
Sajid A, Reimers J R and Ford M J 2018 Phys. Rev. B 97 064101

lopez2018tow
López-Morales G I, Proscia N V, López G E, Meriles C A and Menon V M
  2018 arXiv preprint arXiv:1811.05924

garg1985eff
Garg A, Onuchic J N and Ambegaokar V 1985 J. Chem. Phys. 83
  4491–4503

thorwart2000con
Thorwart M, Hartmann L, Goychuk I and Hänggi P 2000 J. Mod. Opt.
  47 2905–2919

norambuena2016mic
Norambuena A, Reyes S A, Mejía-Lopéz J, Gali A and Maze J R 2016 Phys. Rev. B 94 134305

wiktor2016abs
Wiktor J and Pasquarello A 2016 Phys. Rev. B 94 245411

jakubczyk2016imp
Jakubczyk T, Delmonte V, Fischbach S, Wigger D, Reiter D E, Mermillod Q,
  Schnauber P, Kaganskiy A, Schulze J H, Strittmatter A, Langbein W, Kuhn T,
  Reitzenstein S and Kasprzak J 2016 ACS Photonics 3 2461–2466















Supplementary Material





§ EXPERIMENTAL DETAILS




 §.§ Sample preparation


Single-photon emitters (SPE) in hexagonal boron nitride are observed in hBN nano-powder (Sigma-Aldrich, grain size <150 nm). The powder is micromechanically exfoliated using scotch tape and stamped onto a Silicon substrate with a 80 nm or 270 nm thick thermal oxide layer on top. Figure <ref> shows exemplary images of the nano powder sample. Panels (a) and (b) show a reflectance and PL image of the same sample area, respectively. The dark spots in (a) are clusters of the hBN powder crystals. The PL map in (b) shows bright spots at the positions of the clusters, where each spot stems from one or more defect centers. The SEM image in (c) shows the detailed stacked disc-like structure of the individual hBN crystals.





 §.§ Photoluminescence and photoluminescence excitation spectroscopy

The photoluminescence experiments are carried out in a homebuilt confocal microscope. For excitation at 2.76 eV (450 nm) a diode laser and at 2.33 eV (532 nm) a frequency-doubled Nd:YAG laser is used. For recording the PLE spectra, the emitters are excited with a tunable laser light source based on a continuous-wave (cw) optical parametric oscillator (OPO). The turn-key system (from Hübner Photonics) covers the wavelength range of 450 nm – 650 nm (2.76 eV – 1.91 eV).

For the room-temperature measurements, the laser is focussed to the diffraction limit by an objective lens (numerical aperture NA=0.9), resulting in a focus size of ≈ 400 nm. The excitation power is 100 µW and kept constant for all excitation energies.

For the temperature-dependent measurements, the sample resides in a flow-cryostat cooled by liquid Helium. In this experiment, the laser is focussed on the sample using an objective lens with a NA of 0.75. In both cases, the same objective lens is used to collect the photoluminescence of the single emitters and the PL is analyzed in a spectrometer with an attached liquid nitrogen-cooled CCD camera. A longpass filter at 2.06 eV is used to remove stray light of the laser from the PL signal.




§ THEORY





 §.§ Phenomenological model


As explained in the main text we use a phenomenological fit function to extract relative spectral contributions from the phonon side bands of 165 different emitters.

The full fit function consists of three parts with

    I(E) = I_ ZPL(E) + I_ LO(E) + I_ LE(E) ,

where the zero phonon line (ZPL) is given by

    I_ ZPL(E) =  A_ ZPLγ_ ZPL/(E-E_ ZPL)^2 + γ_ ZPL^2 .

The longitudinal optical (LO) phonon side bands are assumed to have the same width as the ZPL, which leads to the function

    I_ LO(E)    =    I_ LO_1(E)+ I_ LO_2(E)
       =   A_ LO_1γ_ ZPL/[E-(E_ ZPL - E_ LO_1 )]^2 + γ_ ZPL^2
       +   A_ LO_2γ_ ZPL/[E-(E_ ZPL - E_ LO_2 )]^2 + γ_ ZPL^2 ,

where E_ LO_1,2 are the LO phonon energies.

Because the considered low energy (LE) phonon contributions have energies of E_ LE_1=14 meV and E_ LE_2=30 meV they are in the range of the thermal energy at room temperature k_ BT≈ 25 meV. Therefore we take phonon emission (I_ LE_j(E)) and absorption (I_ LE_j^ abs(E)) processes into account. Additionally we consider two phonon processes (I_ LE_j^2 ph(E)) for these modes and combination of LO and LE processes. Therefore, the LE phonon side bands read

    I_ LE(E) = ∑_j=1,2( I_ LE_j(E) + I_ LE_j^ abs(E) + I_ LE_j^2 ph(E) ) ,

The single LE phonon lines read

    I_ LE_j(E)    =   A_ LE_jγ_ LE_j/[E-(E_ ZPL - E_ LE_j )]^2 + γ_ LE_j^2 ,
    
    I_ LE_j^ abs(E)    =   n_j/n_j+1A_ LE_jγ_ LE_j/[E-(E_ ZPL + E_ LE_j )]^2 + γ_ LE_j^2 ,

where n_j is the thermal occupation of LE mode j from the Bose-Einstein distribution. to reduce the number of fit parameters we assume that the weights of the two phonon processes scale like the single LE phonon processes. This leads to the fit function

    I_ LE_j^2 ph(E)    =   ∑_k=1,2{A_ LO_k/A_ ZPL[  I_ LE_j(E-E_ LO_k) + I_ LE_j^ abs(E-E_ LO_k) ] 
                +  A_ LE_k/A_ ZPL[  I_ LE_j(E-E_ LE_k) + I_ LE_j^ abs(E-E_ LE_k) ] } .





 §.§ Model parameters

Motivated by Ref. <cit.>, for the LO phonon energies we consider two constant energies and for the LA phonons linear dispersions. We distinguish between the in-plane and the out-of-plane direction with


    ħω_ LO_1(q_r) = E_ LO_1 ,        ħω_ LO_2(q_r) = E_ LO_2 ,        ω_ LA(q_r) = c_r q_r
    ħω_ LO_1(q_z) = E_ LO_1 ,         ħω_ LO_2(q_z) = E_ LO_1 ,        ω_ LA(q_z) = c_z q_z


The sound velocities are chosen to c_z=3.44 nm/ps <cit.> for the out-of-plane direction and for the in-plane direction we consider the mean value of the two high symmetry directions with c_r=16 nm/ps, which we extract from the dispersion relations in Ref. <cit.>.
Also the dielectric constants are different in the two distinct lattice directions with <cit.>


    ε_∞(q_r) = 4.95 ,           ε_ s(q_r) = 7.04 ,
    ε_∞(q_z) = 4.10 ,           ε_ s(q_z) = 5.09


We interpolate between the two given directions via

    ω_ j( q)    =   cos^2(φ_q) ω_ j(q_r) + sin^2(φ_q)ω_ j(q_z) 
    ε( q)    =   cos^2(φ_q) ε(q_r) + sin^2(φ_q)ε(q_z)

where φ_q =  atan(q_z/q_r) is the angle of q with respect to the hBN-plane.




§ TEMPERATURE-DEPENDENT PHOTOLUMINESCENCE MEASUREMENT



To confirm the conclusions from the temperature-dependent measurements drawn in the main text, in Fig. <ref> we show a second data set for a different localized emitter.



In Fig. <ref>(a) the measured PL spectra (color) are shown for decreasing temperatures from top to bottom. The corresponding simulations are shown as black lines. The parameters used are given in the plot next to each curve. Here only the dephasing time T_2 is adjusted for each temperature, while the other parameters are kept constant. We did not find strong indications for a contribution of the local mode, leading to a coupling strength of g=0 in the simulations. The deformation potential couplings are the same as in the main text (D_ e=0.4 eV and D_ h=-0.2 eV). We again find that the ZPL strongly narrows when reducing the temperature.

In addition we measured the lifetime of the PL signal for different temperatures. These measurements are carried out in the same setup as the low-temperature PL measurements (see Sec. 1.2). However, the emitters are excited by a tunable femtosecond fiber laser system (≈ 250 fs pulse length, 40 MHz repetition rate) at an energy of 2.16 eV <cit.>. The photoluminescence is detected using a Picoquant PDM Series single-photon counting module (timing accuracy 50 ps) and the PL decay measurement is performed with a Becker & Hickl SPC-130-EM time-correlated single-photon counting card. The excitation power is kept constant for all measurements at 150 µW. The results are shown in Fig. <ref>(b). We extract the lifetime T_1 by fitting the data with a convolution of the instrument response function and a single exponential decay (black curves). We find that the lifetime is almost constant T_1≈ 2.3 ns for all considered temperatures. These values are on the same time scale as those of defect centers in diamond <cit.>. 

In Refs. <cit.> it was shown, that the homogenous linewidth of the emitters at low temperatures is in the sub-µeV range and the corresponding T_2  in the ns range. In this case, the dephasing is dominated by the lifetime of the excited state, i.e., additional dephasing effects such as spectral wandering or pure dephasing vanish. This supports our observed trend of a strong narrowing of the ZPL with decreasing temperature. However, the values for low temperatures (5 K and 50 K) are underestimated because the measured line width is on the order of the limited resolution of the spectrometer (100 µeV). Hence, we are not able to draw reliable conclusions about the low-temperature T_2 times from our data. At room temperature, we find that the radiative lifetime is four orders of magnitude longer than the dephasing time. Therefore, we can conclude that the dephasing due to the radiative lifetime is negligible compared to other dephasing mechanisms.




§ COLOR CENTERS IN DIAMOND


Figure <ref> shows PL spectra of the nitrogen vacancy (NV^-) and the H3 color center in diamond in orange, respectively. The single-crystal diamond plate used for this experiment has a lateral size of 3× 3 mm^2 and a thickness of 1 mm and is produced by a high-pressure, high-temperature process (element6) and contains approximately 200 ppm nitrogen. For the excitation of the NV^- center (H3-center) we used a 532 nm laser diode (405 nm laser diode) with a power of 100 µW at the sample position. The focus diameter is 1.2 µm for 532 nm excitation and 1.4 µm for 405 nm excitation. Both spectra were measured at cryogenic temperatures. The blue curves represent the calculated spectra, where the coupling to a single local mode was considered. The model parameters with the best agreement with the experiment are listed in Tab. <ref>.









§ PLE MEASUREMENTS


In Fig. <ref>(a) PL spectra for different excitation energies are shown. The ZPL maxima and the phonon sidebands reside on a background signal. For excitation energies between 2.09 eV and 2.22 eV, Raman lines appear, which shift with the laser energy. The Raman shifts of the three lines are 63 meV, 118 meV, and 168 meV. We attributed the first two to the Si substrate <cit.> and the third one to the LO_1 energy <cit.>. To extract the intensity of the ZPL and the phonon sidebands, we fit the background by an exponential multiplied by an error function to reproduce the edge of the filter. The Raman lines are modeled by Gaussians. The background- and Raman-line-corrected spectra are shown in Fig. <ref>(b).

For small excitation energies around 2.1 eV, which are very close to the ZPL energy, The PL spectra are very weak. Therefore both PLE signals in Fig. 7 in the main text drop significantly. This surprising observation indicates that in the energetic range where excitation processes - including the local mode - should occur, the emitter is not excited. Looking at the PL spectra in Fig. <ref>(a) we see that for small excitation energies the ZPL overlaps with a strong Raman line, which might spoil the evaluation of the PLE data. Apart from that, the deviation between the PLE spectrum and the mirrored emission spectrum remains unclear and deserves additional investigation which is however beyond the scope of this paper. 






§ REFERENCES


10
url<#>1#1urlprefixURL 



serrano2007vib
Serrano J, Bosak A, Arenal R, Krisch M, Watanabe K, Taniguchi T, Kanda H, Rubio
  A and Wirtz L 2007 Phys. Rev. Lett. 98 095503

bosak2006ela
Bosak A, Serrano J, Krisch M, Watanabe K, Taniguchi T and Kanda H 2006 Phys. Rev. B 73 041402

michel2009the
Michel K H and Verberck B 2009 Phys. Rev. B 80 224301

geick1966nor
Geick R, Perry C H and Rupprecht G 1966 Phys. Rev. 146 543

sotier2009fem
Sotier F, Thomay T, Hanke T, Korger J, Mahapatra S, Frey A, Brunner K,
  Bratschitsch R and Leitenstorfer A 2009 Nature Phys. 5 352

batalov2008tem
Batalov A, Zierl C, Gaebel T, Neumann P, Chan I Y, Balasubramanian G, Hemmer
  P R, Jelezko F and Wrachtrup J 2008 Phys. Rev. Lett. 100
  077401

dietrich2017nar
Dietrich A, Bürk M, Steiger E S, Antoniuk L, Tran T T, Nguyen M,
  Aharonovich I, Jelezko F and Kubanek A 2017 arXiv:1712.06938

tran2017res
Tran T T, Kianinia M, Nguyen M, Kim S, Xu Z Q, Kubanek A, Toth M and
  Aharonovich I 2017 ACS Photonics 5 295–300

uchinokura1972ram
Uchinokura K, Sekine T and Matsuura E 1972 Solid State Commun. 11
  47–49

reich2005res
Reich S, Ferrari A C, Arenal R, Loiseau A, Bello I and Robertson J 2005 Phys. Rev. B 71 205201





