The nitrogen-vacancy center (NV) in diamond, with its exceptional spin coherence and convenience in optical spin initialization and readout, is increasingly used both as a quantum sensor and as a building block for quantum networks. Employing photonic structures for maximizing the photon collection efficiency in these applications typically leads to broadened optical linewidths for the emitters, which are commonly created via nitrogen ion implantation. With studies showing that only native nitrogen atoms contribute to optically coherent NVs, a natural conclusion is to either avoid implantation completely or substitute nitrogen implantation by an alternative approach to vacancy creation. Here, we demonstrate that implantation of carbon ions yields a comparable density of NVs as implantation of nitrogen ions and that it results in NV populations with narrow optical linewidths and low charge-noise levels even in thin diamond microstructures. We measure a median NV linewidth of 150 MHz for structures thinner than 5 μm, with no trend of increasing linewidths down to the thinnest measured structure of 1.9 μm. We propose a modified NV creation procedure in which the implantation is carried out after instead of before the diamond fabrication processes and confirm our results in multiple samples implanted with different ion energies and fluences.
The negatively charged nitrogen vacancy center (NV) in diamond is a widely applicable system for both quantum communication1 and quantum sensing.2 It features an optically addressable electron spin with long spin coherence both at low- and at room temperature,3–5 hyperfine-coupled nuclear spins that can be used as long-term storage qubits,6 and an electronic structure that enables sensing of magnetic, electric, or mechanical fields as well as temperature.7 Nanoscale probing of magnetism in 2D materials,8 imaging of domain walls in antiferromagnets,9 and measurement of current flow in graphene10 have benefited from this sensitivity, while demonstrations of entanglement between distant NVs11 and implementations of few-node quantum networks12 build on the NV's electronic, nuclear, as well as optical coherence. Because these applications rely on optical initialization and interrogation, an efficient optical interface is critical. For NVs deep in the bulk, the solution has so far been to employ solid immersion lenses,13,14 which improve the collection efficiency by up to a factor of 10 over unstructured diamond but do not increase the weak NV zero-phonon line (ZPL) fraction. Approaches based on optical resonators, such as open microcavities,15–21 shorten the radiative lifetime, increase the fraction of photons emitted into the ZPL, and improve the collection efficiency, but require more extensive diamond microfabrication. For quantum sensing, the common approach to increase the collection efficiency is to use waveguiding photonic structures as scanning probes,22,23 which also requires a considerable amount of diamond microstructuring.
On account of its permanent electric dipole moment, the NV's excited state is sensitive to charge noise and, thus, prone to optical linewidth broadening, especially in microfabricated diamond.24 Lifetime-limited optical linewidths, ∼13 MHz, have, therefore, been measured only on native NVs in bulk diamond.25 It is important to differentiate between two cases of reported linewidths: the ones measured under repeated off-resonant repumping of the charge- and spin state, which we refer to as extrinsically broadened linewidths, and the ones measured between repumping pulses, which we refer to as dephasing linewidths. In the former case, charge noise triggered by the repump pulses induces spectral jumps and results in a broadened linewidth, while in the latter case this contribution is avoided. Recent work demonstrated close to lifetime-limited dephasing linewidths (<60 MHz) for implanted NVs when using conditional charge state repumping.26 Under repeated repumping, most implanted centers in bulk instead exhibit extrinsically broadened linewidths on the order of 100 MHz, about ten times the lifetime limit.21,27 For micro- and nano-fabricated diamond, the sensitivity to fluctuating electric fields results in a further deterioration of the optical quality of the NVs: remaining lattice damage from implantation as well as surface and subsurface charge traps introduced during etching conspire such that extrinsically broadened linewidths reaching up to several GHz are observed.21,26,28–30 Furthermore, studies have shown that NVs created from implanted nitrogen have worse optical properties compared to NVs created from diamond-native nitrogen and irradiation-induced vacancies, in that they form a separate, broad-linewidth distribution.31,32 This has led to a number of studies exploring different NV creation methods where nitrogen ion implantation is avoided, including laser writing33,34 and electron irradiation.24 Such approaches have resulted in improved NV optical coherence in bulk and, in the case of electron irradiation, also in 3.8 μm-thick diamond membranes, but at the cost of loss of depth control (electron irradiation) and the need for highly specialized setups (laser writing). In general, NV optical linewidths have been shown to increase drastically for structure thicknesses less than 3.8 μm, with only a few cases of narrow lines reported in thinner structures.32,35
In this study, we explore a refined method of NV creation, in which we avoid implantation of nitrogen ions and instead use 12C+ implantation to create vacancies.36 The NVs are formed by recombination of the created vacancies with native nitrogen in the diamond lattice. Using carbon as a substitute for nitrogen retains the ability to create NVs with nanoscale depth resolution (unlike electron irradiation or optical writing), without creating a broad-linewidth population (unlike nitrogen implantation). Using 12C+, which has a mass comparable to 15N+, for the implantation further ensures that the depth distribution of vacancies is very similar to the one created by implantation of nitrogen. Since the lattice consists mainly of 12C, this further ensures that no heteroatoms are introduced into the lattice. In contrast to prior work, we opt for implantation post-fabrication (IPF), implanting the ions after fabrication of the microstructures, following the approach used by Kasperczyk et al.32 The rationale behind this reversal is to avoid exposing already formed NVs to the potentially deleterious effects of fabrication, in particular, the aggressive dry etching steps and electron beam lithography.37,38
We fabricate our samples from electronic grade (100) single-crystal diamond ([N] < 5 ppb, [B] < 1 ppb, Element Six) and follow a previously developed fabrication procedure.21,39,40 The main steps of the process flow are illustrated in Fig. 1(a). We create <5 μm-thin microstructures out of ∼50 μm bulk diamond via electron beam lithography and inductively coupled plasma reactive ion etching (ICP-RIE). The etching is carried out through the aperture (typically 1 × 1 mm2) of a quartz mask until the structures within the aperture become free-standing. A combination of the etching process, mask geometry and mask position results in a thickness gradient: the structures in the center of the aperture are up to several micrometers thinner compared to the structures at its edges. The thickness gradient varies in its extent from sample to sample. Once the microstructures have been created, a thin chromium/gold layer is evaporated at a slight angle onto the back surface to reduce charging during implantation. Finally, 12C+ is implanted into the microstructures from the front side at an angle of 7°, and NVs are created after stripping the metal layer and performing a multi-step annealing process similar to the one described in Ref. 39.
Key aspects of the carbon implantation post-fabrication procedure. (a) Main fabrication steps: (i) A bulk diamond is used as a starting material. (ii) Thin microstructures are patterned on the front side of the diamond via lithography and etching. The structures are then released by a back-side etch through a quartz mask. (iii) The microstructures are implanted with 12C+ from the front side. (b) SRIM-simulated density of implanted 12C+ ions (top) and resulting vacancies (bottom) for an implantation energy of 12 (squares) and 50 keV (circles), as a function of the depth. A fluence of 1 × 1010 and 5 × 108 ions/cm2 was taken for the implantation energy of 12 and 50 keV, respectively. The depths corresponding to distribution maxima are stated together with the difference in depth compared to implanting 15N+ with the same energy and fluence. The 50 keV data were multiplied by a factor of 10 for better visibility.
Key aspects of the carbon implantation post-fabrication procedure. (a) Main fabrication steps: (i) A bulk diamond is used as a starting material. (ii) Thin microstructures are patterned on the front side of the diamond via lithography and etching. The structures are then released by a back-side etch through a quartz mask. (iii) The microstructures are implanted with 12C+ from the front side. (b) SRIM-simulated density of implanted 12C+ ions (top) and resulting vacancies (bottom) for an implantation energy of 12 (squares) and 50 keV (circles), as a function of the depth. A fluence of 1 × 1010 and 5 × 108 ions/cm2 was taken for the implantation energy of 12 and 50 keV, respectively. The depths corresponding to distribution maxima are stated together with the difference in depth compared to implanting 15N+ with the same energy and fluence. The 50 keV data were multiplied by a factor of 10 for better visibility.
As can be seen in Fig. 1(b), Stopping and Range of Ions in Matter (SRIM)-simulations for 12C+ show that the ion- and vacancy distributions after implantation differ only very marginally from implantation with 15N+. For both 12 and 50 keV implantation energies, the difference in implantation depth of ions (vacancies) inside the diamond, defined as with being the depth corresponding to the maximum of the distributions, results in a relative depth difference below 20% (15%). We choose the implantation energies 12 and 50 keV for two different purposes: the former to create shallow NVs suitable for optomechanics experiments requiring coherent, near-surface NVs in cantilevers,41 and the latter to create deeper NVs located at a vacuum electric field antinode in an open microcavity.21 As can be seen in Fig. 1(b), we expect the majority of the resulting NVs to be distributed between 15.3 and 20.4 nm and between 64.6 and 79.9 nm from the top surface of the diamond, respectively. According to our simulations, the expected vacancy yield per implanted ion is virtually the same for carbon and nitrogen when implanting with 12 (50) keV: 74 (175) vacancies/ion for 15N+ and 68 (151) vacancies/ion for 12C+.
Following this procedure, we fabricate three samples, denoted sample A, B, and C. Samples A and B have square, 20 × 20 μm2 microstructures and sample C has rectangular, 5–10 μm wide and 50–100 μm long microstructures. We create NVs at different depths via carbon IPF and study the linewidths of the resulting NVs as a function of the microstructure thickness. Sample A and B are both 12C+-implanted with an energy of 55 keV and a fluence of ions/cm2, while sample C is implanted with an energy of 12 keV and a fluence of 1 × 1010 ions/cm2.
We image the emitters using a home-built confocal microscope, where we use a 532 nm diode-pumped solid-state laser for samples A and B and a 515 nm diode laser module for sample C for off-resonant excitation while detecting the resulting photoluminescence (PL) after filtering by a long-pass filter (Semrock, 594 nm RazorEdge) on a single-photon counting photodetector (Excelitas). Figure 2 shows brightfield images as well as confocal scans of representative regions on the three samples.
Sample overview. (a)–(c) Brightfield microscope images of samples A–C, respectively, with labels indicating the regions which were sampled for linewidth measurements. The areas shown in the images correspond to the apertures in the quartz masks used in the etching process. (d–f) Representative confocal fluorescence image of NVs in regions A2, B1, and C1, respectively.
Sample overview. (a)–(c) Brightfield microscope images of samples A–C, respectively, with labels indicating the regions which were sampled for linewidth measurements. The areas shown in the images correspond to the apertures in the quartz masks used in the etching process. (d–f) Representative confocal fluorescence image of NVs in regions A2, B1, and C1, respectively.
We record NV photoluminescence (PL) spectra on a cryogenically cooled CCD camera coupled to a grating spectrometer (Princeton Instruments). An example spectrum is shown in Fig. 3(a). The thickness of the measured microstructures is determined either from the spectra by fitting the amplitude modulation of the NV's phonon sideband (PSB) arising from interference in the diamond or by direct measurement with a commercial laser confocal microscope (Keyence VK-X1100). We choose measurement sites on the three samples covering a range of different thicknesses. The chosen regions are indicated in Figs. 2(a)–2(c): in sample A, thicknesses span 1.9–4.6 μm and in sample B 1.9–4.9 μm. The free-standing microstructures in sample C are 2.5 μm thick and we also measure in a neighboring bulk region with a thickness of approximately 50 μm (not indicated).
Measurement examples. (a) PL of an NV in sample A. The PSB has been fitted to extract the wavelengths that give constructive interference within the microstructure (dotted lines), in turn providing the structure thickness (5.4 μm). The inset illustrates the etaloning effect allowing for thickness determination. (b) PLE of the same NV as in (a). The inset shows the PLE measurement sequence.
Measurement examples. (a) PL of an NV in sample A. The PSB has been fitted to extract the wavelengths that give constructive interference within the microstructure (dotted lines), in turn providing the structure thickness (5.4 μm). The inset illustrates the etaloning effect allowing for thickness determination. (b) PLE of the same NV as in (a). The inset shows the PLE measurement sequence.
We characterize the NV ZPL optical linewidths through photoluminescence excitation (PLE) measurements in either a liquid helium bath (CryoVac) or a closed-loop cold-finger (attocube attoDRY800) cryostat. We resonantly excite the NVs by sweeping a tunable narrow-linewidth external cavity diode laser (New Focus Velocity TLB-6704) locked to a wavemeter (HighFinesse WS-U) across the ZPL and collect the photons emitted into the PSB with a high-numerical aperture objective (Partec 50×, 0.82 NA, or Olympus LMPLFLN 100×, 0.8 NA). We minimize power broadening in the resulting PLE spectra by measuring at the lowest possible resonant laser power that still gives an adequate signal-to-noise ratio. Typical resonant laser powers are on the order of 100 nW, which leads to a power broadening of up to three times the homogeneous linewidth yet remains marginal compared to our measured linewidths.34 The measurement sequence consists of a 1–2 μs green initialization pulse to reset the charge- and spin state, followed by a resonant pulse (637 nm ± detuning, 5–10 μs) during which PSB counts are accumulated. The sequence is repeated at a 100 kHz rate for a duration of 10 ms for each laser detuning. An example PLE measurement from the average of 100 scans is shown in Fig. 3(b). The residual PLE background arises via leakage of photons generated during the green repumping pulse into the counting window. We note that by systematically using an off-resonant green initialization pulse, our experiment yields a measurement of the extrinsically broadened linewidths, i.e., linewidths including spectral jumps caused by the repump pulses.
Figure 4 summarizes the measured PLE linewidths in samples A, B, and C. We measure a median (mean) Gaussian-fitted full width at half maximum (FWHM) linewidth of 143 (227 MHz), 138 (181 MHz), and 304 MHz (691 MHz) for the three samples, respectively. The data are well-described by lognormal distributions.32 We argue that the median is a better figure-of-merit than the mean as the linewidth distributions are asymmetric. A few broad-linewidth NVs in the tail of the distribution will increase the mean but barely affect the amount of preselection needed in any application requiring narrow-linewidth NVs. In this case, the median is more informative, providing the linewidth value for which half of the population is equal or narrower. Figure 4(d) shows the corresponding empirical cumulative distribution functions, indicating the fraction of measured linewidths below or equal to a certain value.
Linewidth statistics. (a)–(c) Histograms over the measured linewidths for sample A (implantation energy of 50 keV), B (implantation energy of 50 keV), and C (implantation energy of 12 keV), respectively. The solid red lines represent lognormal fits and the median M is stated for each distribution. (d) Empirical cumulative distribution functions for all samples, plotted on a logarithmic scale for better data visibility. The dashed lines represent a 95% confidence interval.
Linewidth statistics. (a)–(c) Histograms over the measured linewidths for sample A (implantation energy of 50 keV), B (implantation energy of 50 keV), and C (implantation energy of 12 keV), respectively. The solid red lines represent lognormal fits and the median M is stated for each distribution. (d) Empirical cumulative distribution functions for all samples, plotted on a logarithmic scale for better data visibility. The dashed lines represent a 95% confidence interval.
In order to provide a quantitative threshold under which NV linewidths can be described as narrow, we base ourselves on the projected indistinguishability between two ZPL photons—a cornerstone for remote entanglement generation.42 Two-photon interference with high (0.9) visibility can still be obtained despite broadening at the cost of temporal filtering, down to the typical ∼300 ps timing resolution of silicon single photon counting diodes.43–45 In our case, we cannot discriminate between broadening due to spectral jumps and broadening due to pure dephasing. Considering the latter as the sole contribution leads to the worst-case scenario, in which temporal selection is technically possible for measured linewidths 150 MHz. In samples A and B, which we group as the two samples with the same implantation parameters, 54% of the measured linewidths are below 150 MHz. In sample C, only 26% of the linewidths are below 150 MHz, and the median linewidth is about two times higher. This difference can be explained by the shallower implantation depth of sample C and a correspondingly larger contribution of surface-related charge noise to the linewidth broadening.26 For the three samples combined, including bulk, we measure 48% of linewidths below 150 MHz.
Figure 5 demonstrates the statistical changes in linewidths as a function of the thickness of the measured structure. We see no general trend of increasing mean (not shown) or median linewidth with decreasing structure thickness, in contrast to previous studies.21,24,32 Sample C shows broader lines than the other two samples regardless of the sampled thickness, which can again be attributed to the shallower NV depth in this sample. Considering all the sampled microstructures, spanning 1.9–4.9 μm in thickness, 52% of the NVs have a linewidth below 150 MHz. No thinner regions were created on samples A and C. On sample B, the fraction of NVs exhibiting optically addressable ZPLs dropped significantly for structures less than 1.9 μm thick, which precluded further studies. We note that narrow, <250 MHz, extrinsically broadened linewidths have previously not been measured consistently in diamond microstructures thinner than 3.8 μm.24,32,35 IPF with nitrogen as employed by Kasperczyk et al. described a few cases of narrow linewidths in thin structures; as we show here, going one step further and substituting nitrogen for carbon avoids the formation of broader-linewidth NVs from implanted nitrogen and, therefore, results in only the native-nitrogen, narrow-linewidth NV distribution.
Median linewidth as a function of diamond thickness. The error bars correspond to the geometrical standard deviation of a lognormal fit to each distribution. The data points are labeled according to the sample and corresponding microstructure, in order of increasing thickness within each sample. The inset shows a close-up of the median linewidths of the four thinnest microstructures.
Median linewidth as a function of diamond thickness. The error bars correspond to the geometrical standard deviation of a lognormal fit to each distribution. The data points are labeled according to the sample and corresponding microstructure, in order of increasing thickness within each sample. The inset shows a close-up of the median linewidths of the four thinnest microstructures.
We speculate that a possible explanation for the positive effect of IPF on NV optical linewidths is that the dynamics of NV formation are affected by preexisting lattice damage. Performing etching and electron beam lithography on structures with existing emitters imposes lattice damage on the complete structure, while implantation after the microfabrication might preferentially create emitters on sites with less lattice disorder. Dedicated studies are, however, required to single out the effect of IPF from carbon substitution alone.
In conclusion, we demonstrate that carbon IPF creates NVs with on average lower charge-noise and reduced spectral diffusion compared to other NV creation approaches. We show that substituting nitrogen with carbon during implantation results in similar NV densities for moderate implantation fluences that allow optical isolation of single emitters. We measure narrow linewidths with high probability: 48% of the measured NVs have an extrinsically broadened linewidth below 150 MHz, and this applies across the full measured range of diamond microstructure thicknesses down to 1.9 μm. We show that by implanting an ion species other than nitrogen, the creation of NVs displaying extreme linewidth broadening is avoided, confirming earlier conjectures.31,32 Moreover, we demonstrate that our method is easy to implement, is reproducible, and provides an advantage over conventional fabrication methods for different samples and implantation parameters.
Our method should prove useful for any application relying on optically coherent NVs, both in bulk and in microstructured diamond. As an example, coupling the emission of such NVs (∼150 MHz linewidth, ∼2 μm thick microstructure) to an open microcavity with a finesse of a few thousand would increase the fraction of photons emitted into the ZPL by more than an order of magnitude.21 The success probability of spin–spin entanglement between remote NVs would, therefore, increase by over two orders of magnitude compared to the state-of-the-art when using the Barrett–Kok protocol.11,24,46 This would dramatically improve the generation rate of entangled qubit pairs and remove an important roadblock toward quantum networks beyond few-node prototypes.
The authors thank Tomasz Jakubczyk for fruitful discussions. They acknowledge the financial support from the National Centre of Competence in Research (NCCR) Quantum Science and Technology (QSIT), a competence center funded by the Swiss National Science Foundation (SNF), and by SNF project No. 188521. J.A.Z. acknowledges the financial support from the Swiss Nanoscience Institute (SNI). A.C. acknowledges the financial support from the Quantum Science and Technologies at the European Campus (QUSTEC) project of the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie Grant Agreement No. 847471.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Viktoria Yurgens: Data curation (lead); Formal analysis (lead); Investigation (lead); Software (equal); Visualization (lead); Writing – original draft (lead); Writing – review & editing (lead). Andrea Corazza: Data curation (equal); Investigation (equal); Writing – review & editing (equal). Josh A. Zuber: Data curation (equal); Formal analysis (equal); Investigation (equal); Writing – review & editing (equal). Marie Gruet: Formal analysis (equal); Software (equal); Writing – review & editing (equal). Mark Kasperczyk: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Writing – review & editing (equal). Brendan J. Shields: Conceptualization (equal); Investigation (equal); Writing – review & editing (equal). Richard J. Warburton: Conceptualization (lead); Funding acquisition (lead); Project administration (lead); Resources (lead); Supervision (lead); Validation (lead); Writing – review & editing (equal). Yannik Fontana: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Project administration (equal); Software (equal); Supervision (equal); Validation (equal); Writing – review & editing (equal). Patrick Maletinsky: Conceptualization (lead); Funding acquisition (lead); Methodology (equal); Project administration (lead); Resources (lead); Supervision (lead); Validation (lead); Writing – review & editing (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.