Scaling laws of the cavity enhancement for nitrogen-vacancy centers in diamond
We employ a fiber-based optical microcavity with high finesse to study the enhancement of phonon sideband fluorescence of nitrogen-vacancy centers in nanodiamonds. Harnessing the full tunability and open access of the resonator, we explicitly demonstrate the scaling laws of the Purcell enhancement by varying both the mode volume and the quality factor over a large range. While changes in the emission lifetime remain small in the regime of a broadband emitter, we observe an increase of the emission spectral density by up to a factor of 300. This gives a direct measure of the Purcell factor that could be achieved with this resonator and an emitter whose linewidth is narrower than the cavity linewidth. Our results show a method for the realization of wavelength-tunable narrow-band single-photon sources and demonstrate a system that has the potential to reach the strong-coupling regime.
pacs:42.50.Pq, 42.50.Ct, 42.81.Qb, 61.72.jn, 81.05.ug
The negatively charged nitrogen-vacancy (NV) center in diamond stands out due to stable single photon emission Gruber et al. (1997); Kurtsiefer et al. (2000) and exceptional spin coherence of electronic ground-state levels Balasubramanian et al. (2009). Together with spin-selective optical transitions, the NV center represents a particularly promising solid-state quantum bit with direct optical access Gruber et al. (1997); Balasubramanian et al. (2009); Maurer et al. (2012). An important milestone yet to achieve is the realization of an efficient optical interface for this emitter. Ideally it would open the way for deterministic single photon sources Su et al. (2008), non-destructive spin state detection Young et al. (2009); Robledo et al. (2011), efficient spin-photon entanglement Togan et al. (2010), and quantum coherent optical manipulation on the single photon level. A variety of optical microcavities has been studied in this respect, including photonic crystal cavities Wolters et al. (2010); Englund et al. (2010); van der Sar et al. (2011); Riedrich-Möller et al. (2012); Faraon et al. (2012), microring resonators Faraon et al. (2011); Hausmann et al. (2012), microdiscs Barclay et al. (2009a, b), Fabry-Perot cavities Dumeige et al. (2011), and plasmonic nanoresonators de Leon et al. (2012); Choy et al. (2011). Considerable success has already been achieved, e.g. by demonstrating Purcell enhancement of the zero phonon line (ZPL) at cryogenic temperature, where of the emission was channeled into a single cavity mode Faraon et al. (2011). The broad emission spectrum and the small weight of the ZPL call for particularly large Purcell factors to achieve efficient cavity coupling. Furthermore, several issues remain problematic for the techniques used so far, which all rely on solid state cavities with limited flexibility to access the cavity mode and to tune the resonance frequency.
In this work we study the cavity enhancement of the phonon sideband emission of small ensembles of NV centers in individual nanodiamonds by means of a fiber based Fabry-Perot microcavity Hunger et al. (2010). The cavity design relies on concave surface profiles machined on the end faces of optical fibers using CO laser processing Hunger et al. (2012). The technique provides excellent surface quality to support ultra-low loss mirror coatings, and offers microscopic cavity mode volumes to achieve large Purcell factors. The open resonator geometry permits us to easily access the cavity field maximum, to study various emitters with one and the same cavity, and to fully tune the cavity resonance. These features have been beneficial for recent experiments with ultracold atoms Colombe et al. (2007), trapped ions Steiner et al. (2013), molecules Toninelli et al. (2010), quantum dots Muller et al. (2009); Barbour et al. (2011); Di et al. (2012); Miguel-Sánchez et al. (2013), and NV centers in diamond Albrecht et al. (2013). Here we use these properties to explicitly demonstrate the scaling of the fluorescence enhancement for NV centers with the cavity parameters. Using direct comparison of free space and cavity enhanced emission, we determine the ideal and effective Purcell factors as a function of the mode volume and the quality factor. While narrow cavity resonances couple only to a small fraction of the broad emission spectrum, we observe an up to 300-fold enhancement of the emission spectral density, giving a direct measure of the ideal Purcell factor of the coupled system. Furthermore, we find that we can maintain a cavity finesse in the presence of rather large nanocrystals.
Setup. In our experiments we use a cavity which consists of a mirror on a concave end face of a single mode fiber and a macroscopic plane mirror, schematically shown in Fig. 1(b). The effective radius of curvature in the center of the profile machined on the fiber is m. Together, the two mirrors define a fully accessible microscopic cavity with a mode waist m and an effective cavity length as small as m, which includes the penetration of the cavity mode into the dielectric mirrors. The smallest mode volume we currently achieve in this way is m. The coatings are designed for highest reflectivity at 780 nm and are realized by multiple layer pairs of TaO and SiO. The finesse reaches a value of at 780 nm for small cavity length. It decreases continuously towards at m and drops faster for larger distances, which we attribute to diffraction loss due to the finite mirror size. The largest quality factor we thus obtain is for m.
We introduce nanodiamonds containing NV centers into the cavity by directly spin coating an aqueous colloidal solution on the large mirror with low enough concentration (m) to isolate single nanocrystals with the cavity mode. For optimal coupling between the emitter and the cavity, the mirror is terminated with a layer of SiO to shift the electric field maximum of the cavity mode slightly above the mirror surface. We employ two different samples of fluorescent nanodiamonds (FNDs) with size distributions peaking at 30 and 100 nm diameter. They were irradiated with He ions to increase the concentration of NV centers Mohan et al. (2010); Yu et al. (2005). Surface treatment with acids and thermal oxygen etching removed graphitic surface shells and other contaminants. The resulting surfaces are predominantly oxygen terminated, which favors the negative charge state of the NV center Rondin et al. (2010); Hauf et al. (2011). The emission spectrum of the NV center shows a zero phonon line at 637 nm with a room temperature linewidth of THz Fu et al. (2009), and a broad phonon sideband resulting from coupling to a phonon continuum with a distinct coupling strength maximum around a phonon frequency of THz Davies (1974). Transitions involving several phonons constitute the sideband fluorescence as schematically shown in Figure 1(a). The strength of the integrated -th sideband can be estimated by the Frank-Condon factors for displaced harmonic oscillator states, , with a Huang-Rhys factor found for NVs in bulk. In our experiments, we observe spectra peaking at a central wavelength of 690 nm with a full width at half maximum of nm [see Fig. 2(a)].
With our present mirror coatings, the NV emission spectrum is divided into two parts: For wavelengths shorter than nm, the mirrors are transmissive and the color centers emit approximately under free space conditions. For longer wavelengths, the mirrors become increasingly reflective and cavity enhanced emission can be observed. This enables us to observe free space and cavity enhanced emission within one measurement, such that the cavity enhancement can be quantified from a direct comparison. Furthermore, it strongly facilitates the search for NV centers in the cavity, since a large fraction of their emission is transmitted. High transparency at the excitation wavelength provides nearly constant excitation conditions when varying the cavity length. The transmission of the plane mirror is chosen to be a factor 2 larger than that of the fiber mirror to optimize outcoupling efficiency into the detection channel used in the experiments. Finally, the strong change of mirror reflectivity across the fluorescence spectrum permits us to study the cavity enhancement as a function of the quality factor.
The cavity setup is combined with a homebuilt confocal microscope as shown in Fig. 1(c). We use an excitation laser at 532 nm (Cobolt Samba) which is coupled to a singlemode fiber and short pass filtered (Thorlabs FES0550) to remove fiber fluorescence. The light is focused with a long working distance objective (Mitutoyo, NA) through the macroscopic mirror (diameter 0.5 in.) onto the surface supporting the FNDs. The mirror is mounted on a nanopositioning stage (PI 541.2 SL) to perform raster scans and to maximize the overlap of a nanodiamond with the cavity mode. The fluorescence is collected via the same objective, separated from the excitation light with a dichroic mirror (Thorlabs DMLP567), and spectrally filtered (Thorlabs FEL0600, FES0850, NF533-17). Thereafter, the fluorescence is either spatially filtered with a pinhole and led to a Hanbury Brown-Twiss setup with two avalanche photo diodes (Laser Components Count) behind a 50/50 beamsplitter, or coupled via a multimode fiber (Thorlabs HPSC25) to a grating spectrometer (Princeton Instruments Acton SP2500) with a CCD camera (Andor iKon-M). The cavity fiber is mounted on a three axis micropositioning mount. It can be withdrawn from the macroscopic mirror by some millimeters in order to use the setup as a conventional confocal microscope.
Characterization. In a first step we characterize FNDs on Suprasil substrates replacing the macroscopic mirror with the fiber withdrawn. We confirm the presence of negatively charged NV centers by assessing fluorescence spectra, which show no significant contribution of NV. We find that an appreciable fraction () of the fluorescent 30 nm diamonds contain single NVs by observing antibunching in the second order autocorrelation function of the photoluminescence. Comparing the average single emitter count rate with the fluorescence yield of 100-nm FNDs, we find that the latter host 30-300 NV centers per nanocrystal.
Next, we study FNDs prepared on a cavity mirror and with the fiber aligned to form a cavity. We perform confocal scans with the mirror surface facing the fiber and detect the partial transmission through the mirror as shown in Fig. 2(d). We observe strong background fluorescence within the detection spectral band, originating from fluorescence generated in the mirror and the fiber, as well as from Raman scattering in the fiber. To obtain signals much stronger than the mirror background, we use 100 nm nanodiamonds containing ensembles of NV centers. With a coating that prevents excitation light from entering the fiber, excessive background is avoidable and also permits single emitter experiments Albrecht et al. (2013). The observed point spread function for individual FNDs is fit by a Gaussian with a radius of m. This is larger than the value found when the sample faces the objective (m) due to the reduced effective numerical aperture, when focusing through a dielectric medium with refractive index . For time resolved measurements with FNDs on the mirror we find fluorescence lifetimes ranging between 16 and 23 ns. This is on the short side of the typical range for nanodiamonds and can be attributed to a lifetime reduction due to the modified density of states on the mirror.
The observed fluorescence in the measurements shown below corresponds to roughly NV centers, again estimated from the comparison with single emitter fluorescence yield. With the respective FND positioned in the focus of the microscope objective, we align the fiber to form a cavity such that the fundamental mode has optimal overlap with the FND. Figure 2(c) shows a typical spectrum obtained under such conditions. To understand the different parts of the spectrum, we show a typical free space spectrum taken on a Suprasil substrate with negligible background [Fig. 2(a)] and a calculation of the mirror transmission [Fig. 2(b)]. Multiplying the two contributions and scaling with a constant matches the modulated emission observed under cavity conditions for the short-wavelength part [Fig. 2(c)]. Within the stop band of the coating we observe enhanced emission into sharp cavity resonances. Subsequent fundamental modes at nm and nm and a series of weak higher-order transverse modes are visible. The separation of the fundamental modes gives direct access to the effective cavity length m. This length contains the part of the cavity mode penetrating into the mirror stack. From a simulation of the coating at the respective wavelengths we find that the penetration adds up to m. Evaluating the separation of the higher order transverse modes , we infer the effective radius of curvature of the fiber mirror, m, in agreement with white light interferometric measurements of the fiber surface. This determines the mode waist and volume of the plane-concave cavity, given by and , respectively. Since and , the paraxial approximation and thus the given formulas are valid to a high precision.
To characterize the background contribution, we retract the fiber and perform saturation measurements. We observe an intensity dependence that contains a saturating part and a linear contribution from background fluorescence, shown in Fig. 2(e). Fitting the integrated count rate of the transmissive part to the function , with the intensity and the fitting parameters , we find a saturation intensity GW/m for the studied FND, which agrees with typical reported values Kurtsiefer et al. (2000), and a background contribution of in the linear regime of the NV fluorescence. To determine , we measured the excitation mode waist to be m by scanning a FND through the focus and detecting the fluorescence without spatial filtering. Within the mirror stopband and for the cavity present we find a signal to background ratio of about 100 for strong cavity resonances. We correct for the background contributions in the subsequent evaluations. Taken together, the spectrum, saturation, and lifetime measurements confirm that we are indeed coupling NV emission to the cavity.
Effective Purcell Factor. The ideal Purcell factor describes the enhancement of the emission rate of an emitter whose entire fluorescence spectrum is coupled to an optical cavity resonance with quality factor and mode volume ,
Here, is the dipole matrix element, is the local electric field at the position of the emitter, and is the field maximum in the standing wave antinode. For the bad emitter regime, i.e., for broadband emitters where electron-phonon interactions distribute the transition over a spectrum with spectral width much larger than the full cavity linewidth , only a fraction of the emission couples to the cavity. This reduces the Purcell factor to
at the cavity resonance wavelength , where accounts for a detuning from the maximum of the spectrum. In this regime, the effective Purcell factor no longer depends on the quality factor of the cavity but rather on the quality factor of the emitter, .
For an experimental evaluation of the effective Purcell factor we compare background corrected spectra for three different configurations: a typical free space reference spectrum from an FND on a Suprasil substrate (sample facing the objective), and spectra and from the FND under study on the mirror (sample facing the fiber) with and without the fiber mirror being present to form a cavity, respectively. We calculate the expected free space emission spectrum of the FND under study from its emission spectrum obtained through the mirror by matching to via the scaling factor determined by a least square fit. Here, is the transmission of the macroscopic mirror which we infer from a direct transmission measurement and simulate with a transfer matrix method Hood et al. (2001).
With this at hand we can compare the integrated emission into the strongest cavity resonance with the integrated emission under free space conditions , where denotes an integration range of 3 full linewidths (FWHM) around the resonance, and the free space spectrum is integrated from 590 to 770 nm to cover the entire emission. The Purcell factor compares the emission rate into a particular cavity mode with the free space emission into the entire solid angle, such that one has to account for the probability for emitted photons to reach the detector for both cases. For the free space situation in our case, the presence of the Bragg mirror and the substrate modifies the emission pattern. We perform analytical and finite difference time domain simulations to study this aspect. Since we collect without immersion liquid, increased emission into the dielectric beyond the critical angle remains trapped in the medium. Additionally, the angular distribution within the NA of our objective remains nearly unaffected by the refraction through the mirror layers and the substrate in the transmissive spectral range. Consequently, the collection efficiency can be determined from the free space dipole emission pattern, the effective solid angle of the objective as determined from the observed point spread function, and the angle dependent transmission coefficient. Weighting with a dipole orientation distribution (see below) yields a collection efficiency . For a photon emitted into the cavity, we include the probability to leave the cavity through the respective mirror, , where and are the loss of the fiber mirror and the macroscopic mirror respectively. At the wavelength of the cavity resonance, nm, we obtain , where we use the measured values ppm, and ppm from the coating simulation. Losses on the way to the detector and detector efficiency are assumed to be equal for free space and cavity conditions and thus cancel out.
Overall, this yields an experimental value for the effective Purcell factor
For the spectrum shown in Fig. 3(a) we obtain a value . The error is estimated from systematic uncertainties of our evaluation method, mainly originating from differences between the spectra of different FNDs and from an alignment and wavelength dependent collection efficiency. Pulsed excitation measurements with and without cavity show no significant difference of the emission lifetime, as expected for the low value of .
To obtain quantitative agreement between the prediction of Eq. (2) and our measurements evaluated according to Eq. (3), we estimate the influence of the term . Over the studied spectral range, the position of the standing wave maximum shifts from the position of the emitters at the center wavelengh towards the coating for shorter wavelengths. This reduces the coupling in particular close to the edge of the stop band of the mirror, where a phase shift results from the dispersive character of the coating. For the resonance position (nm) and our coating parameters we calculate . Next, the detuning from the maximum of the NV spectrum reduces the coupling by a factor , which we read off from the reference spectrum. Finally, we estimate the effect due to the fixed but random dipole orientations. The measurements were performed with circular polarization and excitation conditions below saturation, such that the angle dependence of the excitation probability leads to a normalized distribution with being the angle between the mirror plane and the dipole axis. This reduces the observed average enhancement by a factor . Evaluating Eq. (2) for the parameters of the measurement (m) we obtain , in good agreement with the experimental value.
To demonstrate the scaling of Eq. (2) explicitly, we perform measurements as shown in Fig. 3(a) and evaluate the effective Purcell factor for the emission into the strongest cavity resonance as a function of the cavity mode volume by changing the mirror separation. A few example spectra are shown in Figure 3(b), and the evaluation of the full data set is shown in (c). The observed behavior matches the prediction of Eq. (2) very well. Notably, approaches unity only for very small mode volumes on the order of . We expect that this regime can be approached with fiber-based Fabry-Perot cavities by further miniaturization.
Ideal Purcell Factor. In contrast to the small fraction of the entire emission that is coupled to the resonator, we observe a strong enhancement of the emission spectral density on cavity resonance. In fact, comparing the spectral density for the cavity and free space case quantifies the ideal Purcell factor . This can be seen by considering a spectral element of the emission with Fourier limited width, which couples to the cavity resonantly and thereby experiences the full enhancement. With the strong variation of our coating properties across the emission spectrum, we can study the enhancement as a function of the quality factor of the cavity. By changing the mirror separation we can locate the cavity resonance at any wavelength and thereby sample quality factors spanning a range of almost three orders of magnitude. Figure 4(a) shows the stepwise tuning of a cavity resonance across one free spectral range for a short cavity corresponding to the mode with eight antinodes between the mirrors (smallest m). For comparison, we show a single spectrum for large cavity length m in Fig. 4(b), where the cavity geometry and the excitation conditions are fixed. There, several consecutive fundamental modes between 700 and 800 nm are observed within the mirror stop band, sampling different within a single measurement.
The ideal Purcell factor can be evaluated by comparing the maximal spectral density on resonance with the free space value at the same wavelength. Since the cavity linewidth is beyond the resolution of our spectrometer, we infer the peak spectral density of a Lorentzian resonance and use the cavity linewidth , where . For nm we measure with a narrowband diode laser, while for low we use broadband light and the spectrometer. A calculation of together with the measurements are shown in Fig. 4(d). The free space spectrum as well as corrections due to the collection efficiency and mirror loss are treated in analogy to the previous evaluation. Together, the ideal Purcell factor is determined experimentally by
In Fig. 4(c) the values obtained from the data sets partially shown in Fig. 4(a),(b) are compared with theory, where we again include and . Since and are evaluated at the same wavelength, the same detuning factor applies for both and thus drops out. For m and , the measurements yield a maximum value of , one of the largest values reported to date. Notably, we observe only minor changes of the finesse (or the quality factor) measured at the largest value when introducing the nanocrystal ().
Rate model. In the previous sections we have approximated the NV emission by a single transition with large effective linewidth to simplify the discussion. A more accurate treatment accounts for the individual transitions making up the emission spectrum and the respective rates involved. The coupled system can therefore be described in the framework of the dissipative Jaynes Cummings theory including excessive dephasing Auffèves et al. (2010); del Valle and Laussy (2011). We consider the dynamics governed by the coherent coupling rate , the cavity decay rate , the decay rate of the emitter , and the excessive dephasing rate . Coupling to phonons distributes the dipole matrix element between the ZPL and individual sideband transitions. To simplify the description of the sideband, a single phonon mode picture is used, where the -th phonon sideband is approximated by a Lorentzian with effective width and transition strength . In this way we can treat the ZPL () and the sideband () in the same manner. For an ensemble, inhomogeneous broadening can contribute an additional dephasing mechanism, which, however, is expected to be small compared to the phonon induced broadening in the sideband. The coupling rate for a single emitter and an individual transition is then given by
For the ensemble, we replace with to account for the orientation distribution. Finally, the energy decay rate of the cavity at the respective transition frequency is given by and MHz for typical nanodiamonds.
Here we account for all contributing transitions with respective frequencies , coupling strengths , total incoherent rates , and cavity detunings . The effective Purcell factor is then given by . One can see that becomes independent of for , supporting the previous definition in Eq. (2).
We note that the coupling of a number of emitters to the cavity does not lead to a change of the coupling rate in the present situation. This is in contrast to the case of a coherent ensemble, where a collective dipole forms that leads to an increased coupling rate . Collective enhancement in the presence of strong dephasing can be quantified by the figure of merit , which compares the enhanced emission rate under ideal conditions with the dephasing rate Temnov and Woggon (2005). While modification of the dynamics becomes apparent for , we find for our parameters, such that no collective effects are expected.
To model the measurements of our experiment, we fit and with THz to reproduce the free space spectrum , and determine all other quantities from measurements. We find THz, contributions of transitions with larger are negligible. The large difference in for the ZPL and the sidebands reflects their slightly different origin. While for the ZPL, only thermally excited low frequency phonons contribute, each sideband is made up by the coupling to a structured phonon continuum extending up to the cutoff at THz Davies (1974). Furthermore, the width of the bands increases with due to the increasing number of possible phonon decay paths. We find good agreement between the predictions of the model and the data, e.g. for the parameters of the measurement shown in Fig. 3(a) we obtain . The obtained transition strengths deviate from the expected Frank-Condon factors, which in part could be attributed to a wavelength dependent collection efficiency in our setup. This would imply an error in the evaluated Purcell factors of up to in both models and the experimental data.
The generalized description of cavity enhancement confirms that the simple picture given previously is a sufficient approximation as long as the effective dephasing dominates all other rates and masks the structure of the phonon sideband. Beyond this, the model provides insight into the different roles of transition branching and dephasing.
Conclusion and Outlook. In our experiment we have studied the scaling laws for the cavity enhancement of the fluorescence of broadband emitters with a flexible architecture, offering high finesse and large ideal Purcell factors compatible with the presence of nanodiamonds. The results indicate a method for the realization of efficient single-photon sources with narrow bandwidth and wide tunability, and suggest that fiber-based microcavities could provide a route towards the strong coupling regime in a cryogenic environment.
While NV centers in nanocrystals show larger homogeneous linewidth and stronger spectral diffusion than in clean bulk samples, dephasing rates MHz for the ZPL appear feasible. Alternatively, diamond membranes or nanomachined high quality bulk diamond could be introduced to achieve better optical and spin properties. For an estimation of the achievable coupling strength we account for the ZPL branching ratio and the splitting of the ZPL into six individual fine structure resonances. For the transitions connecting the ground state with the or levels in the excited state, an overall branching ratio remains for perfect state initialization Manson et al. (2006). Optimization of the laser machining technique is expected to achieve mirror radii of curvature m, which would result in GHz for m. Applying a feasible coating with ppm at 637 nm yields and MHz at this mirror separation. With a penetration depth m at the coating center, the mirror separation still amounts to , providing sufficient space for alignment and tunability. Together, the assumed parameters would reach the strong coupling regime, , and lead to an effective Purcell Factor . This opens the perspective for a fully quantum-coherent spin-photon interface.
We acknowledge support from L. Costa, M. Mader, M. Weber, L. Liebermeister, and D. Burchardt, and helpful discussions with R. Albrecht, C. Becher, and A. Kubanek. This work was supported by the Nanosystems Initiative Munich. T.W.H. gratefully acknowledges support by the Max Planck Foundation.
- Gruber et al. (1997) A. Gruber, A. Drabenstedt, C. Tietz, L. Fleury, J. Wrachtrup, and C. Borczyskowski, Science 276, 2012 (1997).
- Kurtsiefer et al. (2000) C. Kurtsiefer, S. Mayer, P. Zarda, and H. Weinfurter, Phys. Rev. Lett. 85, 290 (2000).
- Balasubramanian et al. (2009) G. Balasubramanian, P. Neumann, D. Twitchen, M. Markham, R. Kolesov, N. Mizuochi, J. Achard, J. Beck, J. Tissler, V. Jacques, et al., Nature Materials 8, 383 (2009).
- Maurer et al. (2012) P. C. Maurer, G. Kucsko, C. Latta, L. Jiang, N. Y. Yao, S. D. Bennett, F. Pastawski, D. Hunger, N. Chisholm, M. Markham, et al., Science 336, 1283 (2012).
- Su et al. (2008) C.-H. Su, A. D. Greentree, and L. C. L. Hollenberg, Optics Express 16, 6240 (2008).
- Young et al. (2009) A. Young, C. Y. Hu, L. Marseglia, J. P. Harrison, J. L. O’Brien, and J. G. Rarity, New J. of Phys. 11, 013007 (2009).
- Robledo et al. (2011) L. Robledo, L. Childress, H. Bernien, B. Hensen, P. Alkemade, and R. Hanson, Nature 477, 574 (2011).
- Togan et al. (2010) E. Togan, Y. Chu, A. S. Trifonov, L. Jiang, J. Maze, L. Childress, M. V. G. Dutt, A. S. Sorensen, P. R. Hemmer, A. S. Zibrov, et al., Nature 466, 730 (2010).
- Wolters et al. (2010) J. Wolters, A. W. Schell, G. Kewes, N. Nüsse, M. Schoengen, H. Döscher, T. Hannappel, B. Löchel, M. Barth, and O. Benson, Appl. Phys. Lett. 97, 141108 (2010).
- Englund et al. (2010) D. Englund, B. Shields, K. Rivoire, F. Hatami, J. Vučković, H. Park, and M. D. Lukin, Nano Letters 10, 3922 (2010).
- van der Sar et al. (2011) T. van der Sar, J. Hagemeier, W. Pfaff, E. C. Heeres, S. M. Thon, H. Kim, P. M. Pietroff, T. H. Oosterkamp, D. Bouwmeester, and R. Hanson, Appl. Phys. Lett 98, 193103 (2011).
- Riedrich-Möller et al. (2012) J. Riedrich-Möller, L. Kipfstuhl, C. Hepp, E. Neu, C. Pauly, F. Mücklich, A. Baur, M. Wandt, S. Wolff, M. Fischer, et al., Nature Nanotechnology 7, 69 (2012).
- Faraon et al. (2012) A. Faraon, C. Santori, Z. Huang, V. M. Acosta, and R. G. Beausoleil, Phys. Rev. Lett. 109, 033604 (2012).
- Faraon et al. (2011) A. Faraon, P. E. Barclay, C. Santori, K.-M. C. Fu, and R. G. Beausoleil, Nature Photonics 5, 301 (2011).
- Hausmann et al. (2012) B. Hausmann, B. Shields, Q. Quan, P. Maletinsky, M. McCutcheon, J. T. Choy, T. M. Babinec, A. Kubanek, A. Yacoby, M. D. Lukin, et al., Nano Letters 12, 1578 (2012).
- Barclay et al. (2009a) P. E. Barclay, K.-M. Fu, C. Santori, and R. G. Beausoleil, Appl. Phys. Lett 95, 191115 (2009a).
- Barclay et al. (2009b) P. E. Barclay, C. Santori, K.-M. C. Fu, R. G. Beausoleil, and O. Painter, Optics Express 10, 8081 (2009b).
- Dumeige et al. (2011) Y. Dumeige, R. AllÃ©aume, P. Grangier, F. Treussart, and J.-F. Roch, New J. of Phys. 13, 025015 (2011).
- de Leon et al. (2012) N. P. de Leon, B. J. Shields, C. L. Yu, D. E. Englund, A. V. Akimov, M. D. Lukin, and H. Park, Phys. Rev. Lett. 108, 226803 (2012).
- Choy et al. (2011) J. T. Choy, B. J. M. Hausmann, T. M. Babinec, I. Bulu, M. Khan, P. Maletinsky, A. Yacoby, and M. Loncar, Nature Photonics 5, 738 (2011).
- Hunger et al. (2010) D. Hunger, T. Steinmetz, Y. Colombe, C. Deutsch, T. W. Hänsch, and J. Reichel, New J. Phys. 12, 065038 (2010).
- Hunger et al. (2012) D. Hunger, C. Deutsch, R. J. Barbour, R. J. Warburton, and J. Reichel, AIP Advances 2, 012119 (2012).
- Colombe et al. (2007) Y. Colombe, T. Steinmetz, G. Dubois, F. Linke, D. Hunger, and J. Reichel, Nature 450, 272 (2007).
- Steiner et al. (2013) M. Steiner, H. M. Meyer, C. Deutsch, J. Reichel, and M. Köhl, Phys. Rev. Lett. 110, 043003 (2013).
- Toninelli et al. (2010) C. Toninelli, Y. Delley, T. Stöferle, A. Renn, S. Götzinger, and V. Sandoghdar, Appl. Phys. Lett. 97, 021107 (2010).
- Muller et al. (2009) A. Muller, E. B. Flagg, M. Metcalfe, J. Lawall, and G. S. Solomon, Appl. Phys. Lett. 95, 173101 (2009).
- Barbour et al. (2011) R. J. Barbour, P. A. Dalgarno, A. Curran, K. M. Nowak, H. J. Baker, D. R. Hall, N. G. Stoltz, P. M. Petroff, and R. J. Warburton, J. Appl. Phys. 110, 053107 (2011).
- Di et al. (2012) Z. Di, H. V. Jones, P. R. Dolan, S. M. Fairclough, M. B. Wincott, J. Fill, G. M. Hughes, and J. M. Smith, New J. of Phys. 14, 103048 (2012).
- Miguel-Sánchez et al. (2013) J. Miguel-Sánchez, A. Reinhard, E. Togan, T. Volz, A. Imamoğlu, B. Bresga, J. Reichel, and J. Estève, New J. Phys. 15, 045002 (2013).
- Albrecht et al. (2013) R. Albrecht, A. Bommer, C. Deutsch, J. Reichel, and C. Becher, Phys. Rev. Lett. 110, 243602 (2013).
- Mohan et al. (2010) N. Mohan, Y.-K. Tzeng, L. Yang, Y.-Y. Chen, Y. Y. Hui, C. Y. Fang, and H. C. Chang, Adv. Mater. 22, 843 (2010).
- Yu et al. (2005) S.-J. Yu, M. W. Kang, H.-C. Chang, K.-M. Chen, and Y.-C. Yu, J. Am. Chem. Soc. 127, 17604 (2005).
- Rondin et al. (2010) L. Rondin, G. Dantelle, A. Slablab, F. Grosshans, F. Treussart, P. Bergonzo, S. Perruchas, T. Gacoin, M. Chaigneau, H.-C. Chang, V. Jacques, J.-F. Roch, Phys. Rev. B 82, 115449 (2010).
- Hauf et al. (2011) M. V. Hauf, B. Grotz, B. Naydenov, M. Dankerl, S. Pezzagna, J. Meijer, F. Jelezko, J. Wrachtrup, M. Stutzmann, F. Reinhard, et al., Phys. Rev. B 83, 081304(R) (2011).
- Fu et al. (2009) K.-C. Fu, C. Santori, P.-E. Barclay, L.-J. Rogers, N.-B. Manson, and R.-G. Beausoleil, Phys. Rev. Lett. 103, 256404 (2009).
- Davies (1974) G. Davies, J. Phys. C: Solid State Phys. 7, 3797 (1974).
- Hood et al. (2001) C. J. Hood, H. J. Kimble, and J. Ye, Phys. Rev. A 64, 033804 (2001).
- Auffèves et al. (2010) A. Auffèves, D. Gerace, J.-M. Gérard, M. Franca Santos, L. C. Andreani, and J.-P. Poizat, Phys. Rev. B 81, 245419 (2010).
- del Valle and Laussy (2011) E. del Valle and F. P. Laussy, Phys. Rev. A 84, 043816 (2011).
- Temnov and Woggon (2005) V. V. Temnov and U. Woggon, Phys. Rev. Lett. 95, 243602 (2005).
- Manson et al. (2006) N. B. Manson, J. P. Harrison, and M. J. Sellars, Phys. Rev. B 74, 104303 (2006).