The nitrogenvacancy (NV) center in diamond is a solidstate defect qubit with favorable coherence time up to room temperature, which could be harnessed in several quantumenhanced sensor and quantum communication applications, and has a potential in quantum simulation and computing. The quantum control largely depends on the intricate details about the electronic structure and states of the NV center, the radiative and nonradiative rates between these states, and the coupling of these states to external spins, electric, magnetic, and strain fields, and temperature. This review shows how firstprinciples calculations contributed to understanding the properties of the NV center and briefly discusses the issues to be solved toward the full ab initio description of solidstate defect qubits.
We briefly introduce a prominent solidstate point defect, the nitrogenvacancy (NV) center in diamond, which acts as a quantum bit, the elementary unit of quantum information processing. A desideratum is then provided on what properties should be determined for understanding solidstate defect quantum bits. Finally, the content of this overview is shortly summarized.
Point defects may introduce levels in the fundamental bandgap of semiconductors or insulators that radically change the optical and magnetic properties of the host material. In particular, these point defects could be paramagnetic, i.e. the electron spin is greater than zero. A primary example of such a point defect is the NV center in diamond [1]. This defect consists of a nitrogen atom substituting a carbon atom near a missing carbon atom in diamond crystal, i.e. vacancy of diamond. It can accept an electron from the environment and can be negatively charged (see Figure 1A). In this charge state, multiple levels appear in the fundamental bandgap of diamond occupied by four electrons. The corresponding electron configurations constitute an S=1 ground state and an optically active S=1 excited state with below bandgap (5.4 eV) excitation energy [2], [3], [4]. As a consequence, the defect has singlephoton absorption and emission spectra in the visible region in the transparent diamond host. Thus, the NV center is also called color center in diamond. The absorption and emission spectra are broad even at cryogenic temperatures caused by the coupling of phonons to the optical transitions. The nophononline or zerophononline (ZPL) optical transition appears at 637 nm (1.945 eV) [5]. In particular, the contribution of the ZPL emission to the total emission, i.e. the DebyeWaller (DW) factor, is about 0.03, with a relatively large Stokes shift of about 0.45 eV. This makes possible to excite the NV center by green light (typically 532 nm excitation wavelength) and to detect the emitted photons in the nearinfrared region of 700 to 900 nm. The S=1 spin in the ground state can be measured by conventional electron spin resonance (ESR) techniques [6]. It is important to note that this defect always is a twospin system as the nitrogen isotopes either have I=1 (^{14}N) or I=1/2 (^{15}N) nuclear spin. As a consequence, the hyperfine interaction between the electron spin and the nitrogen nuclear spin always occurs and nuclear quadrupole coupling appears for the largely abundant ^{14}N isotope. In addition, the hyperfine signature of ^{13}C I=1/2 nuclear spins is also observable (natural abundance of 1.1%) in the corresponding ESR spectrum [3], [7], [8]. The spin Hamiltonian of the NV center in the ground state can be written as
where the gyromagnetic factors of the electron spin were first observed as isotropic g_{xx}=g_{yy}=g_{zz}=2.0028±0.0003 (Refs. [7], [9]), whereas slight anisotropy, g_{xx}=g_{yy}=2.0029(2) and g_{zz}=2.0031(2), was reported in Ref. [8], which results in the corresponding Zeeman splitting upon external magnetic field B. D is the socalled zerofield splitting (ZFS) tensor, where “zero field” refers to zero magnetic field. In this particular system, the eigenstates of this part of the Hamiltonian can be given as
A key property of the NV center is that the ESR is correlated with the spinselective fluorescence intensity [13] and twophoton ionization probability [14], where the latter occurs by subsequent absorption of two photons via the ^{3}E excited state. If the concentration of NV centers is low in diamond and the photoexcitation source is focused into a single NV center, then the ESR of the single NV center can be detected either optically [13] or electrically [15] that are called optically detected magnetic resonance (ODMR) and electrically detected magnetic resonance (EDMR), respectively. As the latter was achieved by measuring the photocurrent, it is also called photocurrentbased detected magnetic resonance (PDMR). In ODMR measurements, it was found that the fluorescence intensity between the ms=0 levels of the excited and ground states is about 30% stronger than that between the ms=±1 levels of the excited and ground states. This is also called the ODMR readout contrast. In other words, this readout process converts the ESR frequency used in electron paramagnetic resonance techniques in the microelectronvolt region into the frequency of optical photons in the electronvolt region. In addition, the number of optical photons emitted in the ODMR process is several orders of magnitude larger per center than the single microwave photon absorbed in the ESR process, which increases the detection sensitivity by a large amount. Another consequence of the ODMR process is that the optical pumping of the NV center results in almost 100% population of the ms=0 level in the ground state [6], [7], [16]. It has been recently found that a similar PDMR readout contrast is achievable for single NV centers [15]. In the PDMR readout, the readout process ends at the neutral NV defect [17]. Further optical pumping turns the neutral NV defect into the negatively charged NV defect, i.e. the NV center [18], [19], [20].
The readout and electron spinpolarization processes inherently contain flipping of the electron spin during the decay from the excited state to the ground state, which is highly selective to the ms=±1 spin state of the triplet excited state. Selection rules in the decay processes can be efficiently analyzed by group theory. The NV center in diamond exhibits C_{3}_{v} symmetry with a symmetry axis lying along the ⟨111⟩ directions of diamond [5]. By combining the defectmolecule diagram and group theory [21], [22], [23], one can predict the character of the triplet ^{3}A_{2} ground and ^{3}E excited states as well as the dark singlet states also as a function of external perturbations [24], [25]. Spinorbit interaction may connect the triplet states and singlet states resulting in nonradiative decay and spin flipping. This type of nonradiative decay is called intersystem crossing (ISC). Understanding ISC is the key for ODMR and PDMR readout.
The connection between quantum information science and NV centers is intimately bound to the ODMR and PDMR processes and readout. These readout processes make it possible to coherently manipulate single electron spins in solids by microwave fields and optical excitation [16], [26], where two states of the electron spin of the NV center realize a quantum bit, which can be read out and initialized by optical means. These readout mechanisms operate at room temperature [13], [15] and even at higher temperatures [27], and the measurement scheme can be pushed up to 1000°C with pulsed protocols [28]. The coherence time of the NV center’s electron spin in diamond with natural abundant ^{13}C isotopes can reach ≈600 μs [16], [29] and up to 2 ms in ^{12}Cenriched diamonds [30] as obtained by Hahnecho measurements, even at room temperature. By controlNOT operation [31], [32], the quantum bit information can be written from the electron spin to the nuclear spin. Singleshot readout of the quantum bit was demonstrated [33] as well as the coding of quantum information from the electron spin [34] or from the nuclear spin [35] to the polarization of the emitted photon or the photon emission itself at a given frequency [33], which realize spintophoton interfaces at cryogenic temperature. This was used to transmit quantum information over 1 km distance [36]. One can conclude that the NV center in diamond satisfies DiVincenzo’s criteria of quantum information processing [37]. In the following, we list these criteria that we amend with the NV properties: (1) scalable physical system with wellcharacterized quantum bits, i.e. the electron and nuclear spins of a single NV center; (2) ability to initialize the state of the quantum bits to a simple fiducial state, such as initialization of the ms=0 state by optical means; (3) long relevant decoherence times, much longer than the gate operation time, i.e. millisecond coherence time of the NV center; (4) a “universal” set of quantum gates, e.g. controlNOT operation, as demonstrated for the NV center; (5) a quantum bitspecific measurement capability, i.e. ODMR or PDMR readout; (6) ability to interconvert stationary and flying qubits, such as spintophoton interface of the NV center; and (7) ability to faithfully transmit flying qubits between specified locations as demonstrated between remote NV centers. Scalability is still an issue because identical NV centers require Stark shift tuning of the levels in the excited state [38], [39], [40], and despite the technological efforts on creating arrays or clusters of NV centers in diamond [41], [42], [43], [44], [45], a robust coherent coupling mechanism between multiple NV centers has not yet been demonstrated for more than three NV centers [45].
The NV center in diamond has a favorable electron spin coherence time, but this coherence time and the (spin) levels are relatively sensitive to the environment. This is undesirable for building up a quantum computer from the NV centers but can be harnessed in the measurement of magnetic [29], [46], [47], electric [48], and strain [49], [50], [51], [52] fields and temperature [53], [54], [55] at the nanoscale. The room temperature operation and favorable coherence time of the NV center paves the way toward nuclear magnetic resonance (NMR) of single molecules at ambient conditions [56], [57], [58], [59], [60], [61], [62], [63], [64]. Furthermore, the possibility of spinpolarization transfer from the electron spin toward the nuclear spins [11], [65], [66], [67] can be employed to hyperpolarize diamond particles [68], [69], [70], [71], [72], [73], [74] or external species [75], [76] attached to the diamond surface. In these quantum sensing and related applications, the NV center should reside close to the surface of diamond.
Theory proposed [77], [78] that nearsurface NV centers can be also used as a resource to carry out quantum simulations on frustrated quantum magnetism or quantum spin liquid if the surface of diamond is decorated by species with nuclear spins, e.g. I=1/2 fluorine [77] or I=1 ^{14}N [78], respectively.
NV centers can be observed in natural diamonds, but controlled preparation is required for the aforementioned quantum technology applications. The intentional production of NV defects often starts with a highquality diamond that was grown by chemical vapor deposition (CVD) with minimal contamination. In the next stage, nitrogen ions are implanted into diamond. After implantation, the diamond sample is annealed to remove the damage created by implantation and facilitates the formation of the NV center in diamond [79], [80], [81], [82], [83]. The negative charge state of the NV defect is most likely provided by substitutional nitrogen (N_{s}) donor defects that are called P1 center named after its ESR signature [84], [85], [86]. The shallow NV centers near the surface of diamond are dominantly created by this implantation technique [82], where the depth of the defects can be controlled by the energy of the bombarding ions. It was found that the coherence time and often the photostability of these shallow NV centers are compromised [82], which might be either related to the quality of the diamond surface or the quality of diamond crystal around the shallow NV centers. Indeed, the nitrogen deltadoping technique of diamond growth for creating shallow NV centers improved the coherence time [87].
In some sensing applications, large ensembles of NV centers are needed. In that case, the starting material could be a heavily nitrogendoped diamond and then vacancies can be formed by different irradiation techniques (such as electron [88], [89] or neutron irradiation [17] or implantation with inert ions [90], [91]) in that sample. Finally, annealing is applied to remove the irradiation damages and facilitate the formation of NV defects [92]. It was found that the optimization of the annealing stages is very important to achieve good coherence properties of the resulting NV centers [93].
This brief overview could lead to the impression that the NV center in diamond is probably the most studied point defect in the experiments. Still not all the signatures are well interpreted and understood purely from experimental spectra and simple models. The ab initio theory could significantly contribute to understanding the formation, photoexcitation, and photoionization processes as well as the ISC processes of NV centers. In particular, results from ab initio calculations could guide the idea of improving the coherence properties of NV centers after nitrogen implantation [94] and the notion [14], [20] and optimization [95], [96], [97] of PDMR readout of the NV center in diamond. We provide a recipe below how to achieve the full ab initio description of a solidstate defect quantum bit on the exemplary NV center in diamond.
Two major goals can be identified in the ab initio description of the NV center in diamond: (i) creation of the NV center and its interaction with the other point defects and diamond surfaces and (ii) determination of ionization energies and magnetooptical properties with the corresponding radiative and nonradiative rates, also as a function of external perturbations such as magnetic and electric fields, strain, and temperature.
Issue (i) is a very common target in point defect studies, where the formation energy or enthalpy of the point defects should be determined and the surface of the host should be studied as a function of the environment, which may influence the surface morphology and termination of the host material. However, issue (ii) requires indepth investigation of the behavior of a point defect, much deeper than usual in the community of researchers working on point defects in solids, as the thermodynamic properties of point defects are often the main target in the vast majority of such investigations. However, vacancy formation is not necessarily a quasiequilibrium process. ODMR and PDMR signals of the NV center arise from photoexcitation, which is out of thermal equilibrium of the electrons. In addition, the magnetic properties play a central role in solidstate defect quantum bits that should be described in detail in both the ground and excited states.
The list of properties for the full ab initio description of the NV center can be recognized in Figure 1B. In the ground state, the ZFS due to dipolar electron spinelectron spin interaction should be computed. The hyperfine interaction of the electron spin with the nuclear spins is also very important in the entanglement schemes [16], quantum memory [35], [98], [99], quantum error correction [100], and hyperpolarization [65] schemes as well as in understanding the decoherence of the NV center’s electron spin. For instance, by controlling the proximate nuclear spin states around the NV center, the coherence time of the electron spin has been pushed beyond 1 s [101]. The calculation of the nitrogen quadrupole constant can be important in understanding the socalled “dark state” NMR of ^{14}N [102], which can be used to identify the charge state of the NV defect [10].
In the ^{3}E excited state, the spin levels are heavily temperature dependent and show complex features as observed in the photoluminescence (PL) excitation (PLE) spectrum [103]. The ZFS occurs due to the dipolar spinspin interaction, but also spinorbit interaction takes place. As a consequence, the optically active spin triplet state should be calculated with the corresponding electron dipolar spinspin interaction and spinorbit interaction.
In the ODMR readout process, the calculation of ISC rate requires the spinorbit matrix element between the triplet states and the corresponding singlet states (see Figure 1B). This assumes that the optically inactive or dark singlet excited states and levels should be calculated too.
In the PDMR readout process, the probability of photoionization either directly from the ground state or via the real excited state is a key issue. Photoionization could be direct or via Auger process; thus, the photoionization rate of both processes should be calculated. These nonradiative rates compete with the radiative rate from the ^{3}E excited state to the ^{3}A_{2} ground state, so the latter should be also determined for the sake of a complete description.
Phonons are involved in all of these processes; thus, electronphonon coupling should be calculated to understand the DW factor of the PL spectrum and the ISC processes between the triplet and singlet states [104], [105].
One should recognize that all the interactions between electron orbitals, electron and nuclear spins, phonons, and external fields should be considered and calculated ab initio for a complete description of the operation of solidstate defect quantum bits in a realistic solid and environment. At least, two main challenges can be identified: (a) calculation of the excited states for sufficiently large models of the NV center and (b) treatment of the electronphonon coupling in the radiative and nonradiative processes. In the following, recent efforts along these directions are presented together with the other developments and results in the field.
The rest of the paper is structured as follows. The basic ab initio methods for studying solidstate defect quantum bits are presented in Section 2. We then summarize the recent methodology developments for the description of the NV center in diamond and the corresponding results in Section 3. This section starts with the basic modeling and groundstate properties in thermal equilibrium and then continuous with the treatment of excited states and optical properties, which includes the discussion of the participation of phonons in the optical transition. Next, the computational methods of the magnetic parameters are discussed, which would complete the description of the magnetooptical properties of the NV center in diamond. This is the starting point to determine the radiative and nonradiation rates as described in the next section. Finally, the simulation tools for calculating the various sources of perturbation on the aforementioned properties are presented. An outlook is provided for the next steps toward the full ab initio description of the NV center in diamond in Section 4. Finally, the paper is summarized in Section 5.
The physics of solidstate defect quantum bits is the physics of point defects in solids. As a consequence, understanding solidstate defect quantum bits means to develop and apply tools to explore the properties of point defects in solids. The most employed ab initio technique to study point defects in solids is the planewave supercell KohnSham (KS) density functional theory (DFT) calculation. Two aspects are mentioned in this statement: (i) modeling of point defects in solids and (ii) computational methodology for determining the electronic structure.
In the supercell model, the point defect is placed in a cluster of the host material that has a periodic boundary in each direction; in other words, the cluster with the defect is a unit cell (see Figure 2). If the size of the cluster is sufficiently large, then the defect can be considered as isolated. In practice, the concentration of the defect in the modeling is much higher than that in the experiments because of the limits of the computational capacity. That may lead to the dispersion of the defect levels in the fundamental bandgap, which is a clear sign of the interaction between the periodic images of the defect. We note that the kpoints in the Brillouin zone (BZ) of the supercell are folding into the BZ of the primitive cell [106]; thus, the integration of the reduced number of kpoints in the BZ of the supercell may result in converged wave functions and electron charge density. This reduced number of kpoints can be generated by the MonkhorstPack (MP) scheme [107]. For sufficiently large supercells, a single Γpoint BZ sampling suffices. We note here that only Γpoint calculation guarantees that all the symmetry operations appear for the respective wave functions and charge density, which is an important issue in the investigation of degenerate orbitals and levels. In addition, a practical advantage of Γpoint calculations is that the wave functions are real, which reduces the computational capacity and time. In recent years, the NV center is often embedded into a simple cubic supercell of diamond lattice that originally contains 512 atoms as shown in Figure 2B. In this case, Γpoint sampling is near the absolute convergent kpoint set, but a larger supercell might be needed for highly accurate calculations (e.g. spinorbit coupling in Ref. [105]). We note that the need of a relatively large diamond cluster for the accurate calculation of the NV center presently excludes to use quantum chemistry methods that are based on the extension of the HartreeFock method because of the intractable computational capacity and time or those can be applied with compromising the accuracy by the small size of the diamond cluster [108], [109]. Therefore, another approach has to be applied.
KS DFT [110] has been a very powerful method to determine the ground state of solids. The literature is very rich about KS DFT (e.g. Ref. [111]), which is not repeated here. Briefly, the total energy of the system E is a functional of the electron charge density n(r) of the interacting electron system, where n(r) can be expressed by noninteracting KS single particle wave functions, ϕ_{i} as n(r)=∑_{i}c_{i}ϕ_{i}^{2} and c_{i} is the occupation number of ϕ_{i} KS wave function. The key expression in KS DFT is the exchangecorrelation potential and functional, which is universal for a given number of electrons of the system. In theory, there exists such an exact exchangecorrelation functional. However, this functional is unknown and is approximated in practice. The simplest but unexpectedly successful approximation is the local DFT (LDA), where the exchangecorrelation functional is calculated from that of the homogeneous electron gas at each point [112]. A powerful extension of LDA was achieved by taking the generalized gradient approximation (GGA) of the electron charge density. In particular, the PerdewBurkeErnzerhof (PBE) functional is often applied in the family of GGA functionals [113]. These functionals can be extended to spinpolarized electron systems, where the functionals will also depend on the spin state. In that case, the spin density n_{s}(r)=n↑(r)−n↓(r) with n↑(r) and n↓(r) spinup and spindown densities, respectively, whereas n(r)=n↑(r)+n↓(r), in which the corresponding spinup and spindown KS wave functions and spin densities are varied independently in the sense that the wave functions are not bound to form the spin eigenstate of the system. If the final solution is not a spin eigenstate of the system, then the solution suffers from “spin contamination”. In practice, this “spin contamination” is tiny in LDA or PBE calculations. A key problem of these functionals is that they suffer from the socalled selfinteraction error, which results in too low bandgap of semiconductors or insulators. In particular, the calculated bandgaps of diamond is about 4.2 eV in LDA or PBE DFT calculations. In contrast, the electron charge and spin density of the ground state can be well calculated for many systems, including point defects in solids such as the NV center in diamond [3], [11], [66]. Nevertheless, an improvement in the applied functional or correction to the selfinteraction error was needed to calculate the ionization energies and excitation energies of the NV center in diamond and other solidstate quantum bits. This will be discussed in the next section.
It is important to note that the external potential of the system in KS DFT equations is the potential of ions, i.e. the Coulomb potential of carbon and nitrogen atoms of the NV center in diamond. This method inherently treats ions in a semiclassical fashion in the sense that the ions are classical particles, but the vibration of these classical particles in the selfconsistent adiabatic potential energy surface (APES) can be calculated at the quantum mechanical level. In other words, we apply BornOppenheimer approximation (BOA), which separates the electronic and ionic degrees of freedom. This approximation works generally well for many systems but may fail for degenerate orbitals and levels, where vibrations may effectively couple those states. We will show below that the description of the double degenerate ^{3}E excited state of the NV center requires to go beyond BOA. In the APES, the global minimum energy can be found by minimizing the quantum mechanical forces acting on the ions that can be calculated analytically in KS DFT by applying the HellmannFeynman theorem. This is called geometry optimization procedure for a given electronic configuration. The vibrations can be calculated in the quasiharmonic approximations by moving the ions out of equilibrium and fitting a parabola to the resulting energy differences around the global energy minimum. This procedure sets up a Hessian matrix that can be diagonalized to obtain the vibration (phonon) modes and eigenvectors.
In practice, ϕ_{i} KS wave functions should be expanded by known basis functions or calculated numerically on a grid. In the former case, a natural choice for threedimensional (3D) system with periodic boundary conditions is the planewave basis set, i.e. solution of Bloch states. An advantage of the planewave basis set is that the numerical convergence can be straightforwardly checked by adding moreandmore plane waves to the basis set. The disadvantage of planewave basis set is that very short wavelength (high kinetic energy) plane waves are needed to produce the strongly varying wave functions of the core orbitals (such as C 1s orbital) in a small radius near the ions, which would mean an intractably large basis set. To reduce the computational cost, the ionic Coulomb potentials are replaced so that the effect of the ionic Coulomb potential and the shielding of the core electrons are combined into soft potentials that act on the valence electrons, and the core electrons are not explicitly calculated in the KS DFT procedure. Blöchl worked out the socalled projector augmentation wave (PAW) method [114], which produces a soft potential for the valence electrons but can fully reconstruct the allelectron, meaning the core electron plus accurate valence electron, solution in the core region of ions. This is particularly important for calculating accurate hyperfine constants of point defects in solids [115]. Thus, planewave supercell calculation with the PAW method is a very powerful method for highly numerically convergent KS DFT calculations.
We shortly mention here another modeling approach that was applied to the NV center in diamond (e.g. Ref. [2]). In that approach, the diamond cluster is terminated by hydrogen atoms, i.e. molecular cluster model. As the translation symmetry disappears in this model, the inversion symmetry cannot be maintained in a spherical molecular cluster of diamond. The NV center has no inversion symmetry; thus, this property is not necessarily a disadvantage. In contrast, the surface may introduce the polarization of bonds and extra surface states that are artifacts in modeling isolated NV centers in a perfect diamond host. By applying sufficiently large clusters, the polarization of bonds may disappear at the central part of the cluster where the defect is placed. In contrast, diamond is peculiar in the sense that special surface states appear in hydrogenated diamond (see also Section 3.2). Shockley in his fundamental work [116] already predicted that (hydrogenated) diamond should have negative electron affinity. As a consequence, surface mirror image states [117] or Rydberg states [118], [119], [120] appear in hydrogenated diamond clusters that produce deep empty levels at about 1.7 eV below the conduction band of diamond. Those empty states may mix with the NV center’s empty states, which make the excitation calculation of the NV center in diamond problematic. The maximum probability of the Rydberg wave functions can be found outside of the diamond cluster [119], [120]; thus, those states can only be described by wave functions that are extended and not very much localized on the atoms. In molecular cluster models, the orbitals are often expanded by the linear combination of atomic orbitals (LCAO) that are localized around the ions. Nonorthogonal Gaussiantype orbitals [121] (GTO) are applied in the LCAO calculations as the corresponding integrals in the KS DFT equations can be computed very efficiently compared to those of Slatertype orbitals. These GTO orbitals can be chosen to be well localized that can describe the relatively localized wave functions of the NV center but are not able to describe the Rydberg states properly; thus, the empty levels from Rydberg states do not appear in the bandgap of diamond. Using this trick, a nonconverged basis set for the surface Rydberg states, the NV center in hydrogenated diamond clusters can be modeled as an isolated defect in diamond. In contrast, the valence band and conduction band edges converge very slowly as a size of the molecular cluster toward those of the perfect diamond crystal, i.e. quantum confinement effect (see Ref. [120] specialized to the NV center in diamond); thus, the molecular cluster model does not generally enable the calculation of ionization energies and thresholds. However, it can be an acceptable model for calculating the electron and spin density of the NV center and related properties. We note that allelectron basis can be applied in these GTO calculations without a significant increase of the computational time with respect to that of valenceelectron calculations; thus, the hyperfine tensors can be directly calculated in this approach [122].
In this section, we collect the recent developments on the ab initio calculation of the NV center in diamond. We start with the basic groundstate properties in bulk diamond and diamond surface and then continue with the excited state and related magnetooptical properties, which are used to calculate the corresponding decay rates and coupling parameters to external perturbations.
The formation energy of the defect
where the chemical potential of the carbon atom can be expressed from the total energy of the perfect diamond lattice [E_{tot}(C512)], whereas the chemical potential of the nitrogen atom (μ_{N}) depends on the growth conditions. It can be set to the half of the total energy of the nitrogen molecule as a reference [124]. E_{F} is the Fermi energy referenced to the valence band top E_{V}.
where
which gives the position of the Fermi level with respect to E_{V}, where the concentrations of the defects in charge states q and q+1 are equal,
This result could motivate a recent experimental effort to hinder the formation of divacancy or vacancy clusters by engineering the charge state of the single vacancies so that they are not mobile anymore [94]. In particular, in borondoped diamond, the single vacancies become positively charged or even double positively charged [124], and the formation of vacancy clusters should be significantly reduced. To this end, Fávaro de Oliveira et al. produced a thin borondoped layer near the region where nitrogen ions were implanted [94]. Indeed, this method resulted in an improved charge stability and coherence properties of NV centers in diamond [94]. We note here that atomistic molecular dynamics simulations using empirical potentials showed in this study [94] that indeed divacancies and larger vacancy aggregates are formed near the NV center after N ion implantation when the vacancies are neutral, and the proximate vacancy aggregates’ spin causes the decoherence of the NV center’s electron spin [94].
The shallow implanted NV centers below the diamond surface may suffer from the vacancy cluster formation as explained above, (e.g. Ref. [132]) but the presence of the diamond surface can also pose a problem. Typically, NV centers were implanted into the (100) diamond surface. It was found that shallow NV centers in hydrogenterminated diamond are not stable. This can be explained by the fact that hydrogenated (100) diamond has an electron affinity of about −1.3 eV [133] as predicted by Shockley [116]. When water absorbs to this diamond surface, it creates a band bending of diamond [134], which shifts down the Fermi level and converts the NV center to neutral NV defect [135], [136], [137]. As a consequence, the stability of NV centers requires diamond surface with positive electron affinity. An additional criterion is that no surfacerelated level should appear in the fundamental bandgap of diamond. Ab initio DFT supercell calculations can predict the electron affinity of surfaces of solids by comparing the calculated conduction band minimum (E_{C}) with respect to the vacuum level. In the supercell modeling of surfaces, a slab model is the only option, e.g. diamond slab with (100) bottom and top surfaces [138]. It is often desirable to choose the same type of termination at the bottom and top of the slab to avoid an artificial polarization across the slab, and the size of the vacuum region should be sufficiently large for convergent electron affinity calculations (see Figure 4A). Again, the choice of the functional is crucial. In early calculations, LDA or PBE DFT calculations were applied with too low bandgap [140]. To obtain realistic results, the calculated E_{V} was fixed and E_{C} was shifted to reproduce the experimental bandgap, i.e. scissor correction was applied to E_{C} [140]. However, HSE06 functional calculations with accurate band structure calculation showed [120] that E_{V} should shift down and E_{C} should shift up by about the same amount in the correction of PBE band edges. Thus, the trends could well produce the LDA calculations in Ref. [140], but accurate results are expected from HSE06 electronic structure calculations [120].
In the experiments, oxygenation is applied to stabilize the charge state of shallow NV centers in (100) diamond. It was found that alcohol groups (–OH) would not turn the negative electron affinity of hydrogenated diamond surface to positive electron affinity. Rather, ether bridges (C–O–C bonds) are responsible for this issue [120]. However, ether bridges will stiffen the diamond lattice with introducing defect levels in the gap of diamond, and these groups most likely roughen the diamond surface with creating optically and electrically active defects. Kaviani et al. suggested to apply smooth and wellcontrolled oxygenation techniques that would result in a mixture of C–H, C–OH, and C–O–C bonds at the diamond surface with no stress at the diamond surface but with positive electron affinity and clean bandgap (see Figure 4B). It was also found [120], in agreement with Ref. [140], that smooth fluorine termination provides a large positive electron affinity (100) diamond surface, although an empty band appears just below E_{C} according to HSE06 DFT [120]. Motivated by experimental efforts [141], nitrogenterminated diamond surface was also studied from first principles, and positive electron affinity was predicted for larger than 0.5 monolayer of N–N groups rather than C–H groups at the surface (see Figure 4C). Indeed, a recent study has found improved NV properties on nitrogenterminated diamond (100) surface over oxygenated one, although the improvement was not striking [142]. Most likely, the nitrogen termination of the (100) diamond surface was not perfect in that experiment. HSE06 DFT calculations predicted [78] that nitrogenterminated (111) diamond surface is more favorable than that on (100) diamond surface, as nitrogen can naturally replace the top C–H layer on (111) diamond without introducing any strain [78]. Thus, it is likely that (111) nitrogenterminated diamond surface may host NV centers with excellent properties for sensing applications (Figure 4D). This is appealing as preferential alignment of NV centers along the (111) axis can be realized [143], [144], [145], [146], up to 99% [147], [148], in the growth process of diamond, which is desirable for magnetometry and related applications. However, the growth of highquality (111) diamond at sufficiently high rate is still a big challenge. Alternatively, (113) diamond can be grown with 79% preferential alignment of NV centers at a considerable growth rate [149], [150]. HSE06 DFT calculations have recently obtained a surprising result that nitrogen termination is not preferred for this surface because of introducing surfacerelated bands into the bandgap of diamond, but rather oxygen termination may result in an excellent environment for hosting NV centers for quantum sensor applications [139]. On the (113) diamond surface, oxygen can form socalled epoxy bonds with the surface carbon atoms that are stable according to the ab initio simulations (Figure 4E). This bonding situation does not frustrate the top carbon layers; thus, it produces clean bandgap and positive electron affinity.
The calculated electron affinities in seeking positive electron affinity are presented in Table 1. We note that we do not provide here a comprehensive list of surface terminators creating negative electron affinities that are not relevant for NV sensor applications.
Surface  C–H  C–OH  C–O–C  Mixed  C–N  C–F 

(100)  −1.7  −0.6  +2.4  +0.5  +3.5  +3.0 
(111)  −1.6  +0.2  –  –  +3.2  +3.6 
(113)  −1.8  −0.1  +2.2  +0.5  +3.6  +3.3 
Previously, we mentioned above that the preferential alignment of NV centers have been reported during the ingrowth process of NV centers in the CVD of the (111) and (113) diamond. We note here that a single DFT study predicted [151] that the preferential alignment of NV centers might be achieved in bulk diamond by inducing 2% biaxial strain during the diamond growth at about 970°C temperature; however, this has been not verified in the experiments. The preferential alignment of NV centers in (110) and (111) diamond surfaces was studied by ab initio calculations [152], [153]. It was found that nitrogen prefers to reside on the carbon layer beneath the top carbon layer of the surface [152], and the NV centers will be preferentially aligned at the kink of the diamond terrace during the growth of the (111) diamond under hydrogenrich environments [153].
Realistic diamond surfaces may contain steps, voids, and other defects that were not considered in previous modeling. An exception is the carboxyl group (double carbonoxygen bonds) at the oxygenated (100) diamond surface, which creates levels in the bandgap diamond with localized states that can absorb light in the visible and perturbs NV measurements [120]. Thus, carboxyl groups should be eliminated from the oxygenated diamond surface. Recently, very surprising results have been reported on the (100) diamond surface that could have a direct relation to surfacerelated charge or spin noise [154]. Ab initio modeling tentatively identified stable pairs of sp^{2} C–C bonds in the void defect, i.e. single surface carbon vacancy at the top of diamond surface, which introduces an empty level about 1.4 eV above E_{V} as found by surfacesensitive analysis techniques [154]. This empty state may be filled by electrons during illumination of the diamond surface that provides a surface charge and spin. This defect is very general and relatively abundant on the (100) diamond surface that might be the source of surface noises felt by nearsurface NV centers [154]. The interaction of the NV center with a nearby acceptor defect has been recently modeled in a bulk supercell by DFT calculations [155]. The bulk model was chosen to avoid modeling problems of the slab calculations. The acceptor defect was chosen to be the neutral NV defect for simplicity [155]. A common sense among experimentalists is that the NV center is an atomiclike defect with very localized orbitals. Chou et al. showed in the study of interaction of the NV center with neutral NV acceptor [155] that the NV orbitals spread to 4 nm from the center of the NV defect in the plane of the threecarbon atoms nearest to the vacancy. If the acceptor defect is closer than about 7 nm distance from the NV center with similar wave function extension, then they can directly interact quantum mechanically without any illumination of the NV center, which can lead to a decrease in the spin coherence time of the NV center [155]. The calculated quantum mechanical tunneling rate between the defects could well explain the experimental data in diamond sample with high density of NV defects [156], which highlights the predictive power of DFT methods. Recently, it has been found in experiments [157] that shallow NV centers can ionize in the dark on experimentally relevant timescales, which can be understood as tunneling to a single local electron trap as the mechanism behind this process.
NV centers were also considered close to diamond surfaces using slab models. The modeling of the NV center in slab models faces several problems: (i) the charge correction of the negatively charged defect is not straightforward as the potential of a point charge in a surface is only conditionally convergent; thus, the handling of charged slab supercell is painstaking. (ii) Artificial polarization might appear due to the bottom of the slab (double surface). Lastly, (iii) the finite width of the slab may introduce a quantum confinement effect. Basically, the total energy of the NV center in the slab model may be converged by increasing the lateral size of the slab model [138] and by applying total energy correction techniques that were recently proposed in the literature [128], [158], [159], [160], [161], [162], [163], [164], [165], [166]. The sufficiently large lateral size is necessary to minimize the artificial interaction of the defect with its periodic images too [167] as usual in 3D bulk modeling. The other two problems may be minimized by adding the same surface termination at the bottom and top of the slab and by using a sufficiently large width of the slab [120]. According to intensive tests, about 2.2 nm width produces minute quantum confinement effect [120].
Kaviani et al. invented to do a workaround in the problem of charged slab supercell by replacing it with another but readily solvable problem [120]: a neutral slab model is used for the negatively charged NV defect (NV center) at the expense that another defect enters the slab. Namely, the N_{s} donor will donate an electron to the neutral NV acceptor defect by creating a pair of the NV center and positively charged N_{s}. If these defects are placed into the same layer of the slab and the slab has a cubiclike shape, then the dipoledipole interaction between the periodic images of the defect pairs can be minimized. Rather, the presence of ionized N_{s} near the NV center may shift the levels or split the degenerate levels of the NV center. Extensive tests showed that if the defects are placed at least 7.5 Å far from each other, then the degeneracy of the corresponding orbitals is maintained and the constant shift in the KS levels can be well monitored and corrected [120]. In a recent study, the pair of the NV center and ionized N_{s} was analyzed in detail by arriving at the same conclusion that 7.5 Å distance between the two defects suffices to avoid the splitting of the degeneracy of the corresponding orbitals, and the NV center can be well approximated as isolated [168]. These achievements make possible to directly simulate the NV center with correct total energies that is necessary for modeling direct interaction of the NV center with surface species [120].
We mention in this context that direct tunneling between the nitrogen donor and neutral NV defect has been recently discussed in detail in type 1b diamond [169], which study largely extends the original idea from Collins [170] about the electron transfer between these nearby defects. We note that one conclusion from the ab initio study of the interaction of the NV center and neutral NV acceptor defect that the rate of tunneling depends on the actual orientations of the nearby defects because of the direction of the spatial extension of the orbitals [155]. This should hold for the nitrogen donorNV defect pair too as the neutral nitrogen donor orbitals shows special spatial extension due to the giant JahnTeller distortion (e.g. Refs. [124], [171]); thus, accurate determination of the tunneling rate as a function of distance and relative orientation between the two defects requires ab initio calculations.
The quantum bit operation and readout works via optical excitation; thus, understanding the absorption and decay from the excited state is highly important. To this end, the first task is to calculate the excited states and levels. This is far from trivial. The extension of wave function of the NV center in diamond requires sufficiently large models for accurate calculations [105], [155]. In addition, if diamond band edges are also involved in the photoexcitation process, then the diamond host should be modeled at equal footing with the NV center itself. As a consequence, the highly accurate but extremely costly wave functionbased quantum chemistry methods have only limited accuracy, as these methods can be applied on small molecular cluster models [108], [109]. KS DFT, including HSE06 DFT, which can be used in large supercells in practice, has builtin limitations because of the approximations in the KS DFT functionals; in particular, the highly correlated singlet excited states, e.g. the ^{1}A_{1} state, cannot be directly calculated by KS DFT correctly. Highly correlated states may be recognized in the combination of defectmolecule picture and group theory as multideterminant states [2], [23], [24], [25]. Defect states in solids may be viewed as highly correlated electron states in contact with a bath of extended states. In a recent work, the wave function method is embedded into the DFT framework to synthesize the advantages of both methods and follow the character of the system [172]. For the chosen orbitals (typically defect orbitals), the electronelectron interaction is exactly calculated as Coulomb interaction between the electrons, i.e. configurational interaction (CI), whereas the interaction of the other electrons is treated with HSE06 DFT. The interaction between the chosen orbitals and the extended states is treated within randomphase approximation in a way that the corresponding dielectric function calculated with the constraint of extended states (cRPA) [172]. This CIcRPA method does not contain any fitting parameters, and it is transferable in the sense that the number of orbitals in the CI active space can be systematically increased. This method can be used to analyze the character of the wave functions and the features of the other methods. In particular, the position of the ^{1}A_{1} level is very sensitive to the screening of the electrons from the diamond bands [172]. As a consequence, pure quantum chemistry CI method in tiny diamond cluster models yield too high ^{1}A_{1} level [109]. In a diamond cluster Hubbard model calculation, it was shown that the ^{1}A_{1} state inherits double excitation from the lower a_{1} level to the upper e level in the gap [173], which provides an important insight about the nature of the ^{1}A_{1} state, although the limitation of the diamond cluster model resulted in a false position of the ^{1}A_{1} level in this calculation [172]. For this reason, the very popular manybody perturbation method on top of DFT calculations, called GW+BSE (see Ref. [174] and references therein), fails to properly describe ^{1}A_{1} level [175], [176] because BSE can typically describe excited states with the combination of single excitation configurations. The extended Hubbard model calculation fit to the GW calculation of the ingap defect levels provides accurate singlet and triplet levels of the NV center [172], [176]. We note here that the NV center is exceptional in the sense that the accurate description of the electronic structure of other solidstate quantum bits requires to consider resonance orbitals from the valence band in the CI active space, for which the fitting procedure is ill defined [172]. The final conclusion is that the ^{1}A_{1} level resides about 0.4 eV below the excited state ^{3}E level. We note that the electronic ^{1}E state is basically stable against distortion (it can be described as a single Slater determinant in an appropriate basis of e_{x,y} orbitals in Figure 5A [4], [177]), and it becomes slightly unstable against distortion because of the appearance of ^{1}E′ character in the ^{1}E state [172], [177]. The predominant simple character of the ^{1}E state makes the ^{1}E level almost insensitive to the choice of the computational method that lies about 0.4 eV above the groundstate ^{3}A_{2} level. The final level diagram is shown in Figure 5 in which the global energy minimum of the corresponding APES is taken into account (see below). The final picture well accounts for the measured ZPL energy between the singlet states of about 1.19 eV [178], and derivation from the measured nonradiative rates from the ^{3}E level to the ^{1}A_{1} combined with a phenomenological model on the phonon participation and density of states [104] that concluded about 0.4 eV gap between these two levels.
The CIcRPA method also justified the socalled constraint occupation DFT or ΔSCF method for calculating the ^{3}E excited state [2], [179] of which level is relatively insensitive to the size of the diamond cluster and different computational methods. In the constraint occupation DFT, the electronic configuration is set to build up the Slater determinant of the ^{3}E excited state as shown by the green arrow in Figure 5A. A big advantage of the ΔSCF method is that the forces acting on the ions can be straightforwardly calculated; thus, the APES of the ^{3}E state can be mapped by the KS DFT method. This is of high importance in understanding the nature of the interaction of electrons and phonons in the excited state and its role in the optical spectrum and nonradiative decay.
Experiments indicated [180], [181] that the dynamic JahnTeller (DJT) effect occurs in the ^{3}E excited state, which is double degenerate and has two components:
DFT HSE06 calculations yielded about 42 meV JahnTeller energy and about 9 meV barrier energy between the global minima of APES [105] that might be superior over the DFT PBE results, yielding smaller values [183]. Because of the finite barrier energy, the DJT should be solved in quadratic JahnTeller approximation. The exact solution can be expanded into series as
where the expansion is convergent with maximum of four oscillator quanta (n+m≤4), n, m⟩ is the occupation representation of e_{x} and e_{y} vibrations, respectively, and c_{nm}, d_{nm} coefficients are obtained from the solution of the E⊗e DJT electronphonon Hamiltonian [105]. We note that other systems may require a larger number of oscillator quanta (e.g. the product DJT system of neutral siliconvacancy qubit in Ref. [185]). The tilde over Ψ labels the vibronic or polaronic nature of the wave function: the electronphonon wave function cannot be factorized into electronic and phonon wave functions as assumed in BOA but will be a linear combination of such wave functions. Thus, one has to go beyond BOA in DJT systems, i.e. in strongly coupled electronphonon systems. The vibronic spectrum starts with an
Phonons may couple nondegenerate states too, which is often called the pseudoJahnTeller (PJT) effect. This occurs for the electronic ^{1}A_{1}⟩ and ^{1}E⟩, although the energy spacing between the two is more than 1 eV [177]. The dynamic coupling between these two states may be rationalized by invoking the symmetry breaking e_{x,y} phonon modes, which distort ^{1}E⟩ to ^{1}A′⟩ and ^{1}A″⟩, and ^{1}A_{1}⟩ becomes ^{1}A′⟩, so the ^{1}A′⟩ component of the distorted ^{1}E⟩ can couple to ^{1}A_{1}⟩ by Coulomb interaction. According to the DFT+CIcRPA method, the contribution of ^{1}A_{1}⟩ into ^{1}E⟩ is about 2% in the distorted geometry [172]. At high C_{3}_{v} symmetry, the contribution of ^{1}E′⟩ to ^{1}E⟩ is about 10%, which is an electronelectron correlation effect. We call the resulting multideterminant singlet state
that govern the shape of the phonon sideband in the optical spectra of singlets. We label the symmetry adapted vibrational wavefunctions, e.g. χ_{1}(A_{1})⟩=00⟩, or in general by χ_{i}(…)⟩. The g_{i} coefficients are generally tiny and can be ignored. In contrast, the nonzero d_{i} and
Similar to molecules, vibrations or phonons can contribute to optical transitions of point defects in solids. The optical transitions of molecules are often modeled by the FranckCondon theory, which applies BOA. One can sketch the APES in the electronic ground and excited states together with the phonon states and levels as shown in Figure 8. The most simple case is that the APES can be described as a parabola in both the electronic ground and excited states, so the FranckCondon theory perfectly works. Principally, the two APES parabolas may have different tangents or effective vibration modes, but they may be similar in the electronic ground and excited states. The HuangRhys theory assumes that two APES parabolas are identical, except for a shift in the minimum of the parabolas. In this case, the derivation of the phonon lineshape in the optical excitation spectrum is simplified. In the experiments, the DW factor can be measured, which has a direct relation to the HuangRhys factor S=−ln DW, which measures how many effective phonons participate in the optical transition on average. Generally, the larger the distance between the two APES parabolas, i.e. the larger movement of ions upon optical excitation, the more likely to incorporate phonons in the optical transition. The NV center has a broad absorption and emission spectrum (see Figure 8) even at low temperatures, and DW≈0.03 at cryogenic temperature, which means that about 3.5 effective phonons participate in the radiative decay.
The construction of the absorption and emission spectra in the HuangRhys approximation requires the calculation of APES in both the electronic ground and excited states as well as the vibration modes of the point defect in the electronic ground state. DFT calculations showed that quasilocal vibration modes appear at about 65 meV in which the carbon and nitrogen atoms vibrate the most around the vacancy [182], [187], and these modes appear in the absorption and PL spectra as bump features (see Figure 8). DFT ΔSCF method can be used to calculate the APES in the electronic excited state. By forcing C_{3}_{v} symmetry in the electronic excited state, Alkauskas et al. calculated the PL spectrum of the NV center in diamond [188]. They developed an embedding technique to calculate the contribution of longwavelength phonons [188], and they obtained a good agreement with the experimental data, except for a small feature in the phonon sideband of the PL spectrum (see Figure 8) as pointed out in Ref. [105]. It was shown above that DJT occurs in the ^{3}E state, which dynamically distorts the symmetry. By taking one of the parabolas in the sombrero APES, i.e. the distorted C_{1}_{h} geometry, the PL spectrum can be calculated within the HuangRhys theory, and all the features in the phonon sideband of the PL spectrum were well reproduced [105]. The ab initio theory revealed that about 10% of the phonon sideband emission is associated with the e_{x,y} phonon modes that are responsible for the dynamic distortion in the ^{3}E state [105]. This proves that the DJT feature indeed appears in the PL spectrum of the NV center. Furthermore, the DJT nature of the ^{3}E state was considered to be responsible for the anomalous temperature dependence of the ZPL width [183].
The aforementioned asymmetry in the phonon sideband in the PL and absorption spectra is attributed to the DJT nature of the excited state. We note that a characteristic double peak in the absorption spectrum was associated with a JahnTeller feature [5], but the ab initio theory (see Ref. [188] and private communication with Audrius Alkauskas) produces this double peak by the contribution of A_{1} phonons. The reason behind this asymmetry may be identified as the consequence of the dynamics of ions in DJT systems. By taking the low temperature limit for the sake of simplicity, the electron occupies the lowestenergy
The excited states and excitation spectrum of the triplets can be basically calculated by timedependent DFT (TDDFT), too [190], [191]. The usual linear response approximation was applied in the TDDFT calculations of the NV center [191], [192]. We note that the accurate calculation of the spectrum requires the proper choice of the DFT functional. In 1.4 nm diamond molecular cluster model, the PBE0 [193] hybrid density functional provides quantitatively good results within TDDFT [192]. The quantum mechanical forces acting on the ions in the electronic excited state can be calculated within the TDDFT formalism. The experimental Stokes shift of the NV center of about 0.45 eV could be well reproduced by PBE0 TDDFT calculations (see Supplementary Materials in Ref. [194]) despite the limitations of a molecular cluster model as explained above. This shows the predictive power of the TDDFT method. TDDFT calculations can be combined with molecular dynamics simulation to monitor the evolution of the ionic motions in the electronic excited state. It was found that phonons move the system from the phonon excited state to the zero point energy (lowestenergy phonon level) after photoexcitation from the electronic ground state within 50 fs, which agrees well with the observed decay time from pumpprobe PL measurements [195]. The reason behind the peculiar ultrafast motion is tentatively associated with the DJT nature of the ^{3}E excited state. The observed femtosecond electronic depolarization dynamics of the NV center [181] was associated with the nonadiabatic transitions and phononinduced electronic dephasing between the two components of the ^{3}E state based on ab initio molecular dynamics simulations.
We now discuss the PL spectrum of the singlets, in particular, the phonon sideband. The starting point is that the final state,
It was found from direct calculation of the intensities in Eq. (7) that the optical transition to the first vibronic A_{1} state of the
The absorption spectrum of the singlets should be calculated in a different manner. The PJT effect produces a barrier energy for the damped DJT APES energy surface of the
We previously discussed the orbitals of the NV center and their interaction with phonons. The group theory analysis was very powerful to identify the electronic structure of the ground and excited state orbitals. The fine structure of the NV center can be further analyzed by invoking the analysis of the C_{3}_{v} double group, which treats spinor functions. As the NV center has an even number of electrons (holes) the single group rows in the C_{3}_{v} double group appears as irreducible representation of the spin states [22], [23], [24], [25]. In particular, the ^{3}A_{2} ground state splits to A_{1}ms=0 and E ms=±1 states, whereas the ^{3}E state splits to E_{x,y} ms=0, E_{1,2}ms=±1, and A_{2} and A_{1}ms=±1. These splittings are intrinsic to the C_{3}_{v} symmetry defect, and not the external fields are responsible for their presence. Basically, spinorbit and dipolar electron spinelectron spin interactions may introduce energy spacing between spin levels. In the ^{3}A_{2} ground state, the effective angular momentum (L) of the orbital is zero; therefore, the spinorbit interaction
The spin states and levels behave in a strikingly different way as a function of temperature [103], [197]. In the ground state manifold, the ZFS D constant decreases slightly [−74.2(7) kHz/K] but well resolved [197]. In contrast, the fine structure of the ^{3}E state can be well observed at cryogenic temperatures with λ_{z}=5.33±0.03 GHz [103], [198] and D=1.42 GHz,
The spin Hamiltonian of the spinspin dipolar interaction may be written as
where r_{ij}=r_{i}−r_{j}. The 3×3 Dtensor can be diagonalized to find the spectrum and spin eigenstates. In the C_{3}_{v} symmetry, the D_{zz} component lies in the symmetry axis, which also sets the quantization axis of the electron spin along the symmetry axis of the NV center. In the ground and excited state manifold, the ZFS is D=(3/2)D_{zz} and the corresponding spin levels can be computed as
where r_{a,b}=(r_{1}−r_{2})_{a,b} and χ_{ij} is either 1 or −1 for KS i, j states of the same or different spin channels, respectively. Note that in DFT the KS states are not spin restricted as mentioned previously. Consequently, not only the unpaired KS states but also the rest of the occupied states can contribute to the spin density and the ZFS [199]. The first implementation of Eq. (9) in the DFT supercell version with LCAO basis set was reported in Ref. [199]. The first implementation into planewave supercell code was done by Ivády et al. [200], which only used pseudowave functions in the vasp code. The theory of allelectron PAW version of the Dtensor was reported by Bodrog and Gali [201]. The allelectron version was then implemented by Martijn Marsman into vasp code in 2014. Other implementations also followed later: both the pseudodensity version [202] and the allelectron versions have been recently implemented into Quantum Espresso code [203], [204] or an adaptive finiteelement 3D grid code [205].
The calculated D constant in the ground state was within 3% with respect to the lowtemperature value at ≈2.87 GHz [200] in Γpoint 512atom supercell calculation. By increasing the supercell size, this value did not change, but a larger supercell might be needed for similar defect quantum bits [206]. The ΔSCF method can be used to generate the KS states for calculating the components of the Dtensor in the ^{3}E excited state. By following the previous group theory analysis [24], [25], the calculated D=1.61 GHz,
The spinorbit interaction Hamiltonian in zeroorder approximation can be written as
where V is the nuclear potential energy, m_{e} is the electron mass, and
where λ_{⊥} and λ_{z} are the basal or nonaxial and axial parameters of the interaction, respectively. The former flips the orbital and spin state; thus, it does not play a role in the solution of the double group of C_{3}_{v} symmetry but can play a role in the ISC by coupling ms=±1 of the triplet to a singlet state [24], [25]. The latter will only shift spin levels in the ^{3}E triplet or may couple the ms=0 of the triplet to a singlet state [24], [25].
The axial spinorbit interaction strength in the ^{3}E excited state can be obtained both from total energy difference calculations and from the splitting of KS states in ΔSCF noncollinear spinorbit DFT calculations [105]. It has been shown that the two approaches gave the same result. The latter is more practical because it only calls a single calculation but should be done in the Γpoint of the BZ, as the dispersion and splitting of the defect states in low symmetry kpoints may be larger than the axial spinorbit splitting. To accurately determine small values of the axial spinorbit interaction, the calculations require high numerical convergence and accuracy as the spinorbit energy falls in the μeV regime. Finite size effect turned to be crucial for spinorbit interaction calculations [105] (see Figure 10). Thiering and Gali [105] attributed the observed finite size effect to the overlap of the defect states and used an exponential fit to eliminate the supercell size dependence of λ_{z}.
The calculated HSE06 λ_{z}=15.78 GHz is about three times larger than the data derived from PLE measurements at cryogenic temperature [103]. This difference is not a discrepancy but rather the consequence of the DJT nature of the ^{3}E state as stated in Eq. (5). Ham already suggested that the effective spinorbit energy will be reduced in such systems with a Ham reduction factor p [208]. In the particular ^{3}E state of the NV center,
The ISC between ^{1}A_{1} and ^{3}E states is associated with the perpendicular component of the spinorbit operator. Unlike the case of λ_{z}, one has to apply an approximation to calculate this, namely, that the KS wave functions building up the ^{3}E state and the ^{1}A_{1} multiplet do not change. This enables to compute λ_{⊥}ab initio using the a_{1} and e_{+,−} KS wave functions of the NV(–) in the ^{3}E excited state as
In the NV center [see Eq. (1)], the ^{14}N nuclear spin has a quadrupole splitting, C_{Q}, which can be written as
where h and e are the Planck constant and the charge of the electron and V_{zz} is the gradient of the electric field changes around the nitrogen nucleus in the zdirection (chosen symmetry axis of the defect). Q_{N} is the nuclear electric quadrupole moment of ^{14}N, which scatters between 0.0193 and 0.0208 barn in the literature [210]. V_{zz} can be calculated ab initio within the PAW framework. In this calculation, the criterion for the wave functions (planewave cutoff) and the forces should be set accordingly to achieve convergent V_{zz} (see the Supplementary Information in Ref. [10]). The calculated C_{Q} value for the NV center is −5.019±0.19 MHz, where the uncertainty comes from the uncertainty in the experimental data of Q_{N}. Nevertheless, the ab initio result is very close to the experimental data and could be used to follow the change in the quadrupole moment in NV defect as a function of the charge state [10].
The last remaining terms in Eq. (1) at zero magnetic field are related to the hyperfine interactions. Hyperfine interaction describes the interaction between the electron spin and the nucleus spins that can be generally written as
where A is the hyperfine tensor and
where n_{s}(r) is the electron spin density, r is the vector between the electron spin and nuclear spin at R, g_{n} is the nuclear gfactor, and μ_{n} is the nuclear magneton for a given nucleus n. We note that the nitrogen atom and distant carbon atoms sit in the symmetry axis by creating a hyperfine field that preserves the C_{3}_{v} symmetry. For those nuclear spins, the hyperfine tensor is diagonal with diagonal element A_{xx}=A_{yy}=A_{⊥} and A_{zz}=A_{}. These parameters can be expressed by the Fermicontact term a and a simplified dipolar coupling term b as A_{⊥}=a−b and A_{}=a+2b. The other nuclear spins reduce the symmetry; thus, all the principal values of the hyperfine tensor differ for those nuclear spins. One can recognize from the definition of the Fermicontact term that the spin density at the nucleus site should be calculated very accurately. Thus, allelectron methods should be applied: either treating the core electrons explicitly (typical for GTO codes [122]), PAW method for planewave codes [3], [11], [115], [211] or finiteelement numerical method on a grid [205].
The hyperfine constants of the NV center were first determined by planewave PAW DFT LDA method [3]. It was found that the spin density is not evenly distributed around the core of the defect but rather elongated in the plane of the threecarbon dangling bonds where the amplitude of the spin density decays like sin x/x function, where x is the distance from the vacancy (Figure 11A and B). This spin density distribution has an analog of a Friedel oscillation with the NV center as an impurity. This spin density distribution is responsible for the fact that direct quantum mechanical interaction of the NV center with acceptor defects is much stronger in the basal plane of the threecarbon dangling bonds of the NV center, i.e. in the (111) plane than that along the ⟨111⟩ direction [155]. Another important consequence of this finding is that the Fermicontact term cannot be neglected for ^{13}C spins in the (111) plane near the NV center, but some proximate ^{13}C spins have small hyperfine constants that lie in the node of the spin density oscillation. In this study [3], some ^{13}C spin quantum bits were identified that were employed in a previous entanglement measurement between the electron spin and nuclear spins of the NV center [16].
The spin density significantly redistributes in the ^{3}E excited state and is significantly localized around the nitrogen atom (Figure 11C). As a consequence, the hyperfine constants on the nitrogen spins increases up to ≈50 MHz, which broadens the resonance condition of external magnetic fields for transferring the polarization of the electron spin toward the nitrogen nuclear spin in the excited state (see Refs. [65], [66], [67]).
This type of calculations was repeated by the PAW DFT PBE method [11], [66], where the agreement between the reported hyperfine constants of ^{14}N and ^{13}C nuclear spins was very good [6], [8]. Later, it was found that this was partially fortuitous agreement [211] because of the PAW implementation with frozen core orbitals [115]. The core polarization of carbon 1s orbital significantly contributes to the Fermicontact term, in particular, where the spin density is much localized [211]. Generally, core polarization correction is required in the PAW implementation of accurate hyperfine tensor calculation [211], [212]. DFT HSE06 calculations localized the wave functions more around the vacancy than DFT PBE does, but core polarization correction will reduce the Fermicontact term of ^{13}C spins, resulting in an excellent agreement with the experimental data [211]. It was shown [211] that cancellation of errors (DFT PBE underestimation of spin density and neglect of core polarization) does not generally occur for other defect quantum bits; thus, the implementation of core polarization in the PAW formalism is crucial. In a recent study [122], nonflipping proximate ^{13}C spins in the NV center going from the ground state to the excited state were identified in an 510atom molecular cluster model by means of hybrid functional calculations of hyperfine tensors as implemented in an allelectron GTO code that can be used as longliving quantum memories.
We note that orbital hyperfine interaction may contribute to the hyperfine dipoledipole interaction in defects with an effective angular momentum of the associated orbitals (e.g. Ref. [213]), which was implemented into a Green functionbased code and applied to nickel defects in diamond [214] but not to known defect quantum bits.
By having the magnetooptical parameters in hand, one can discuss the ab initio description of ODMR and PDMR readout of the NV center in diamond. In the ODMR readout, the interplay between the radiative and ISC decay rates from the ^{3}E excited state plays a crucial role that will be first discussed below. In the PDMR readout, the direct photoionization and Auger recombination processes also enter, which we will describe next.
Radiative rates can be readily calculated from KS DFT by taking the corresponding a_{1} and e_{x,y} KS orbitals of the NV center and calculating the
where n_{r} is the refractive index of diamond, ħω is the excitation energy, c is the speed of light, and er is the optical transition dipole moment. One can realize that the ez component will be zero, and only the perpendicular components will give nonzero results as was also found in the experiments [216]. The polarization property of the absorbed and emitted photons is naturally obtained by KS DFT calculations. Finally, Eq. (15) results in about 8 ns radiative lifetime, which is relatively close to the experimental one (~12 ns that corresponds to about (2π)13 MHz rate; see Refs. [104], [217]). For the triplet optical transition, this is a very good approximation because of the relatively simple character of the electronic ground and excited states. The calculation of the optical transition between the singlet states has just become feasible by the DFT+CIcRPA method [172], where the multiplet wave functions are considered in the optical transition dipole moment calculation. The calculated radiative lifetime of the ^{1}A_{1} state is 1878 ns, which is two orders of magnitude longer than that for the triplet [172]. This ab initio result sheds light on the experimental fact that not just the competition with nonradiative decay and the low fractional population of the singlets but the inherently tiny optical transition dipole moment is responsible for the very weak infrared PL signal of the NV center [178]. This agrees well with the conclusion of absorption measurements between the singlets [196]. The very recent pumpprobe PL measurement has observed 100 ps lifetime of the ^{1}A_{1} state [218]. The combination of this recent experimental data with the ab initio radiative rate [172] and observed absorption [196] indicates that the nonradiative decay from the ^{1}A_{1} state to the ^{1}E is much faster than the radiative one. This should be verified by future ab initio calculations of the nonradiative rate between the singlets.
We note that these calculations inherently assume that the FranckCondon principles hold for the optical transitions, i.e. the optical transition dipole moment is independent from the ionic movements. This is a valid approximation for the optical transition between the triplets of the NV center in diamond, but it does not necessarily hold for all types of vibrations and solidstate defect quantum bits, in which the HerzbergTeller interaction is considerable [219]. We showed above that the fluorescence of the singlets of the NV center in diamond does not follow the FranckCondon principles either, and one should use polaronic wave functions in the calculation of the optical spectrum.
The nonradiative rates are associated with the spinorbit interaction between the triplet and singlet states mediated by phonons, i.e. ISC processes. Thus, the ISC rates from the ^{3}E spin states to the ^{1}A_{1} state and from the ^{1}E state to the ^{3}A_{2} spin states should be calculated. In the FranckCondon type of approximation, the ISC rate from the A_{1} state of the ^{3}E state may be calculated [104] as
where λ_{⊥} is given in rad/s unit, Δ is the energy spacing between A_{1}⟩ and ^{1}A_{1}⟩, and δ is the Dirac delta function. χ_{0}⟩ is the ground vibrational level of A_{1}⟩, and
where
The next step is to calculate the ISC rate from the ^{1}E state to the ^{3}A_{2} spin states in the lower branch that will also determine the lifetime of the ^{1}E state. It is known from experiments that the lifetime of ^{1}E〉 is temperature dependent and varies between about 371 and 165 ns going from cryogenic temperature to room temperature [220]. The temperature dependence could be well understood by a stimulated phonon emission process with an energy of 16.6±0.9 meV [220]. On the other hand, the ODMR contrast does not seem to vary significantly as a function of temperature. From PL decay measurements, it was derived that the decay rates toward ms=0 and ±1 favors ms=0 only by about 20%, which was a surprising result but does not contradict with the measured ODMR contrast because of the high spin selectivity in the upper branch [220]. Indeed, the nature of the ^{1}E state is complex: it is a polaronic
where the summation over all vibration wave functions of ^{3}A_{2} collapses to the number of
and
It was again found that λ_{⊥}/λ_{2}=1.2 provides consistent absolute total ISC rate
We emphasize that Eqs. (17) to (22) complete the theory of optical spinpolarization loop of the NV center in diamond where the corresponding matrix elements and coupling parameters could be mostly determined from firstprinciples calculations [105], [172], [177]. The remaining issues are discussed in Section 4.
Next, we discuss the photoionization of the NV center and related defects that are important in the PDMR readout of the NV center’s spin state. The idea of PDMR readout is based on the standard ODMR measurements of a single NV center. The phonon sideband photoexcitation of the NV center (typically 532 nm wavelength) with high laser power can lead to the unintentional ionization of the defect via the absorption of a second photon before decay from the ^{3}E state toward the ground state. The ab initio theory predicts that ≈2.7 eV requires to ionize NV center from the ground state to neutral NV defect at cryogenic temperature. The absorption of the second photon at the ^{3}E state likely occurs at the zeropoint energy as the phonon relaxation in the ^{3}E state is extremely fast [195], which is 1.945 eV above the ground state (ZPL energy of the NV center). Thus, the absorption of the second photon would promote the electron to ≈4.3 eV above the groundstate level of the NV center, which is much larger than the ionization threshold. The twophoton absorption and the corresponding ionization of the NV center will convert it to neutral NV [18]. The ab initio theory proved that this twophoton ionization is an Auger recombination process [20]. It is well known in semiconductor physics that two types of ionization of defects may occur: (i) direct ionization or (ii) Auger recombination. In Auger recombination, multiple defect levels and states may play a role. In the NV center, the direct ionization from the e_{x,y} orbital in the ^{3}E state competes with the Auger process where the highenergy electron promoted to the conduction band will fall back to the hole left in the a_{1} state, and the energy gain is used to kick out an electron from the e_{x,y} level [20]. The Auger recombination rate, i.e. the inverse process of impact ionization, can be written as a Coulomb interaction between the initial exciton state and all the possible final exciton states. The initial state i is an excitation from orbital j to a with spin ↓; here, j is the hole on the a_{1} level and a is a state close to the conduction band edge, which can be simplified as ϕ_{j}↓ϕ_{a}↓ electronhole KS wave functions. The energy of this exciton is approximated as ε_{a}↓−ε_{j}↓, where ε_{rσ} is the KS level of state rσ. The final exciton state is approximated similarly as ϕ_{k}↑ϕ_{b}↑ electronhole KS wave functions, where k denotes the hole on the double degenerate e_{x,y} orbital, b is the electron above the conduction band minimum, and the energy of the exciton is ε_{b}↑−ε_{k}↑. To obtain the rate of transition, we need to sum over all possible final states: the degenerate e_{x,y} level (index k) and the electron in the conduction band minimum (index b). The final formula may be given as (see Supplementary Materials in Ref. [20] and also Ref. [223])
We obtain a single Γ by averaging over initial states:
where ħω is the sum of two energy pulses (twophoton absorption) and V can be generally expressed as
where ϕ_{s,r,u,t} are the corresponding KS wave functions. Vörös et al. [223] implemented the computation of this rate into a planewave supercell code [203], which was applied for the NV center in diamond, and ε_{r} was simply estimated from the refraction index of diamond in this particular calculation [20]. This could be the first ab initio calculation of the Auger recombination rate for point defects in solids. The direct ionization rate can be calculated like the usual absorption rate with optical transition dipole operator just the final state is at the conduction band minimum. The calculated Auger recombination rate was 800 ps, whereas the direct ionization rate was 0.5 μs. A very important consequence of this finding is that the twophoton absorption of the NV center results in the “ground state” of neutral NV defect; thus, no PL signal of neutral NV defect is expected. The typical 532 nm excitation wavelength can excite the neutral NV defect for which the same type of twophoton absorption occurs; finally, it ends at the ground state of the NV center [20]. It was found that the twophoton absorption of the neutral NV defect is very effective at the excitation wavelength of its ZPL wavelength at 575 nm (yellow light). One can conclude that the full cycle of photoexcitation and ionization of the NV center will lead to emission of an electron and a hole, and the final state is the ground state of the NV center.
As the ejection of the electron goes through the ^{3}E excited state in which the spin states have different lifetimes, the rate of ionization (absorption probability of the second photon) should be also spin dependent. This idea was tested on a diamond device where electrodes were built on the diamond structure for the observation of the photocurrent upon the photoionization of the NV centers [14]. They observed a contrast in the photocurrent as they hit the microwave resonance of the electron spin in the ground state of ensemble NV centers [14], [96], i.e. PDMR readout. The full ab initio description of a PDMR cycle would require to calculate the Auger recombination rate for the neutral NV defect, which has not yet been reported. In contrast, the emission of the electron from the NV center and its spin selectivity is now fully understood using the ab initio theory. CIcRPA calculations showed little or no photoionization from the
In the first PDMR measurements [14], the PDMR contrast was relatively low compared to ODMR contrast because of the background current coming from mostly the P1 center in diamond. Therefore, it is highly important to understand the photoionization of the NV center and other parasitic defects to optimize the PDMR readout technique. Recent ab initio calculations have found that green light excitation much favors the ionization of the P1 center over that of the NV center; however, the relative ionization rates can be pushed toward increasing that of the NV center if blue light excitation is applied [95]. In these calculations, the valence and conduction bands were included in the computation of the absorption rates (not just the band edges); therefore, the BZ sampling should go up to 6×6×6 MP kpoint set for convergent calculations in a 512atom supercell [95]. Indeed, dualbeam photoionization of the diamond sample with an ensemble of NV centers resulted in an 3× enhanced PDMR contrast over the singlebeam photoionization scheme, where blue light was employed for direct photoionization and green light was applied to induce spinselective nonradiative decay from the ^{3}E excited state [95]. This again proves the predictive power of ab initio calculations and its role for the improvement of the control and readout of solidstate defect quantum bits. Later on, it has been found that lockin techniques with usual greenbeam excitation can result also in an enhanced PDMR readout contrast [97]. Other parasitic defects producing photoionization bands in the diamond samples could be identified by ab initio calculations [171], which can be rather used as a resource, i.e. PDMR quantum bits, which might not be visible in the ODMR readout but could be effective in the PDMR readout [171].
We note that the carrier capture crosssection rates of NV defect and/or parasitic defects in diamond can be important in understanding the atomistic processes and improving further the PDMR readout. The ab initio framework already exists to calculate these rates [224] but not yet has been applied to the NV center in diamond. The carrier capture crosssection together with the Auger recombination process may play an important role in the electroluminescence of NV defect, which only could successfully employed for the neutral NV defect [225].
Having the computed magnetooptical parameters in hand, one can calculate the coupling parameters to various external fields and temperature. The most obvious type of external perturbation is the constant magnetic field that appears in Eq. (1) known as the Zeeman effect. Generally, the interaction between the electron spin and magnetic field cannot be simply described by the gfactor of the free electron g_{e} but rather by a gtensor, where the Cartesian elements of the gtensor can be obtained from the second derivative of the relativistic manyelectron energy (e.g. Ref. [226]). The implementation of gtensor in supercell planewave codes exist (e.g. gauge including projector augmented wave approach [203], [227], [228], [229]); however, it has not been applied to the NV center in diamond but other defects (e.g. Refs. [230], [231]). For the NV center in diamond, the gtensor does not significantly deviate from that of the free electron. Under external magnetic field, the Zeeman effect appears in both the ^{3}A_{2} ground state and the ^{3}E excited state, but the gtensor in the excited state should differ at low temperatures because spinorbit interaction can contribute to the gtensor in the excited state but is almost negligible in the ground state. As a consequence, PLE measurements at cryogenic temperatures should be able to identify the difference in the gtensor in the ground and excited states in future Zeeman experiments.
When constant (positive) magnetic field is applied along the symmetry axis of the NV center with the magnitude near the Dconstant, then the ms=0 and −1 levels of the triplet state approach each other. The perpendicular component of the hyperfine constant of nitrogen spin opens a gap between these levels, socalled level anticrossing, and the electron and spin states are coupled and rotate with a rate of
In the presence of electric field and strain, the ground state spin Hamiltonian in Eq. (1) is modified. Again, there are symmetryconserving interactions that only shifts the levels, whereas a symmetrybreaking interaction may mix the spin states. Accordingly, the zcomponent of the electric field vector E only shifts the spin levels (transforms as full symmetry A_{1} representation; E_{z}) whereas the x and y components (transform as E representation) can mix the ms=±1 and ms=0, ±1 levels that can be written as
where h is the Planck constant. Eq. (26) has been known as a linear Stark shift for C_{3}_{v}S=1 EPR centers and was derived specifically for the NV center by Doherty et al. [233]. The coefficients d_{⊥}=17 Hz cm/V and d_{z}=0.35 Hz cm/V have been inferred in the experiment of Ref. [234]. However, to our knowledge, the coefficient
The firstprinciples calculations of electric field are not straightforward in the 3D supercell model as the electric field breaks the periodicity of the system in the direction of the electric field. One possible solution is to apply a slab model for modeling the NV center with a surface that is perpendicular to the direction of the electric field. In Ref. [235], the NV center in diamond was modeled by a 990atom 2nmthick diamond (111) slab, where the applied electric field was set to 0 to 0.1 V/Å along the symmetry axis of the defect in diamond slab with –OH termination. In this interval, the response was found to this perturbation to be linear on the calculated Dconstant and resulted in d_{z}=0.76 Hz cm/V (see Ref. [235]), which falls to the order of magnitude deduced from experiments. The disadvantages of this method are that (i) the calculation of the perpendicular component of the coupling parameters requires another slab model with an appropriate surface orientation and termination and (ii) the surface termination may induce “artificial” polarization around the NV center even when is buried deepest in the slab model.
Recently, an alternative approach was used to study the effect of electric fields on the properties of the NV center [236]. In this particular case, the coupling of external electric field to the optical transition was examined, which is also known as Stark shift [237]. The Stark shift is generally an undesirable effect for quantum emitters because it can cause spectral diffusion (e.g. Ref. [20]): during optical excitation of the NV center, nearby defects may be ionized that creates different electric fields around the NV center that results in a shift of ZPL of the NV center in subsequent optical excitation cycles. Another consequence is that zero (magnetic) field ODMR signals of ensemble NV centers will show a splitting due to the various charged defects around the NV centers [238]. For determining the Stark shift of the NV center, the permanent dipole moments should be calculated in the ^{3}E excited state and the ^{3}A_{2} ground state.
The alternative approach of calculating this quantity relies on the modern theory of polarization [239], [240], [241], [242] that can be applied to periodic systems. The polarization p can be written as
where M is the number of occupied bands, f is the occupation number, e is the elementary charge of the electron, and the wave functions (u_{k} ) are cell periodic and periodic in the reciprocal space. Using density functional perturbation theory, ∇_{k} u_{k}_{n}〉can be calculated from the Sternheimer equations with similar selfconsistent iterations as in the selfconsistent field DFT
By applying this theory to the NV center, the coupling parameter of the electric field to the optical transition is Δp=〈^{3}Ep^{3}E〉–〈^{3}A_{2}p^{3}A_{2}〉=0.903eÅ≈4.33 Debye, and it directs toward the NV line along the ⟨111⟩ direction. This value corresponds to 2.18 MHz cm/V, which means that a pointlike charge (e.g. positively charged N_{s}) at 50 nm distance from the NV center shifts the ZPL energy by 2.4 GHz. This result agrees well with previous experimental results [243], where Δp≈1.5 Debye was deduced for the NV center in diamond.
Basically, the advanced theory of polarization enables to calculate the d_{z}, d_{⊥}, and
The ^{3}E excited state is much more sensitive to the presence of electric field or strain as the degenerate orbitals may split in the presence of symmetry breaking fields, and they are intertwined [24], [25]. Therefore, the effect of strain is first discussed before discussing the effect of electric field in the ^{3}E excited state.
We note that the effect of hydrostatic pressure on the spin levels in the ground state was first studied in experiments [244] and the ab initio theory [200] before determining the spinstress coupling parameters. In the theoretical study, the lattice constant of the 512atom supercell was varied to induce the appropriate pressure on the NV center and the Dconstant was calculated within pseudowave function approximation [200]. The experimentally deduced linear coupling coefficient is 14.58 MHz/GPa, which was measured at room temperature [244]. The results of ab initio simulations are valid at T=0 K and yielded 10.3 MHz/GPa (see Ref. [200]). It was found in the ab initio simulations that the local relaxation of ions and adaption of the wave functions to the new ionic positions upon external pressure plays a crucial role in the final response; thus, ab initio simulations are essential in the calculation of the spinstress parameters. We further note that the pressure shifts the energy of the ZPL optical transition too by 5.5 meV/GPa [245]. DFT PBE ΔSCF calculations yielded 5.75 meV/GPa shift [246], in good agreement with the experimental data. Furthermore, they predicted that the emission in the phonon sideband will be significantly suppressed and redistributed toward the ZPL emission under giant (100 GPa) pressure [246]. This has not yet been confirmed in experiments to our knowledge. Next, we turn to the spinstrain coupling parameters.
The spinstrain coupling in the ground state of the NV center is more complicated than spinelectric field coupling as strain is a tensor and can result in more coupling coefficients than the electric field does. Indeed, the correct spinstrain coupling coefficients have been only recently derived, which contains six independent real couplingstrength parameters, h_{41}, h_{43}, h_{25}, h_{26}, h_{15}, and h_{16}, and has the following form [237]:
where ε_{ij}=(∂u_{i}/∂x_{j}+∂u_{j}/∂x_{i})/2 denotes the strain tensor and u(r) is the displacement field. Similar to Eq. (26), subscripts 0, 1, and 2 here refer to the difference in the electron spin quantum numbers ms connected by the corresponding Hamiltonian. The spinstress Hamiltonian H_{σ} is analogous to Eq. (29a), with substitutions ε↦σ and h↦g. The conversion from h to g can be calculated using the stiffness tensor (C) of diamond, which connects the external stress to the internal strain in diamond (e.g. Ref. [237]) with the parameters C_{11}=1076 GPa, C_{12}=125 GPa, and C_{44}=576 GPa in the cubic reference frame.
The six spinstrain couplingstrength parameters of Eq. (29a) were determined using numerical DFT calculations [237]. To model the structure subject to mechanical strain, described by the strain tensor ε , the cubic supercell was deformed to a parallelepiped, whose edge vectors are obtained by transforming the undeformed edge vectors with the matrix 1+ε in the cubic reference frame and allow the atomic positions to relax. For each strain configuration, the elements of the 3×3 ZFS matrix D were calculated using the VASP implementation by Martijn Marsman with the PAW formalism [201].
We illustrate this methodology to obtain the six spinstrain couplingstrength coefficients with the example of h_{16}. To determine h_{16}, the supercell using a strain tensor is deformed whose only nonvanishing element is ε_{yz}, and the D matrix is obtained from the calculation. Due to Eq. (29a), the chosen strain configuration implies that the Hamiltonian has the form
This, together with the above definition of the D matrix, yields
The numerical error of our DFT calculations was estimated using a linear fit to the data points that were taken from 11 equidistant values of ε_{yz} between −0.01 and +0.01 (see Ref. [237]). The final results for all the coupling parameters are summarized in Table 2.
Parameter (MHz/strain)  Parameter (MHz/GPa)  

h_{43}  2300±200  g_{43}  2.4±0.2 
h_{41}  −6420±90  g_{41}  −5.17±0.07 
h_{25}  −2600±80  g_{25}  −2.17±0.07 
h_{26}  −2830±70  g_{26}  −2.58±0.06 
h_{15}  −5700±200  g_{15}  3.6±0.1 
h_{16}  −19,660±90  g_{16}  18.98±0.09 
In Table 3, we compare the numerical DFT results of Table 2 to the experimental results of Ref. [247]. In Ref. [247], four of the six independent spinstress couplingstrength parameters of the spinstress interaction Hamiltonian were measured. Ref. [247] defined these four spinstress couplingstrength parameters, denoted as a_{1}, a_{2}, b, and c, in a “hybrid” representation, where the spinstress Hamiltonian is expressed in terms of the NV frame components of the spin vector (S_{x}, S_{y}, and S_{z}) and the cubic frame components of the stress tensor (σ_{XX}, σ_{XY}, etc.). To be able to make a comparison between the DFT results and the experimental ones, the notations of Ref. [247] were taken and d, e,
Parameters  Relation  DFT (MHz/GPa)  Experimental results (MHz/GPa) 

a_{1} 

−2.66±0.07  −4.4±0.2 
a_{2} 

2.51±0.06  3.7±0.2 
b 

1.94±0.02  2.3±0.3 
c 

−2.83±0.03  −3.5±0.3 
d 

−0.12±0.01  – 
e 

0.66±0.01  – 
where
The relations between the hybrid representation parameters (a_{1}, a_{2}, b, c, d, and e) and the NV frame parameters (g_{41}, etc.) are given in the first two columns of Table 3. Importantly, H_{σ}_{0} and H_{σ}_{2} are identical to the spinstress Hamiltonian in Eqs. (1) and (2) of Ref. [247]. It was found (Table 3) that the d and e parameters are in the same order of magnitude as the other four parameters; thus, this full sixparameter description is required to deduce accurate parameters from experimental data. From the calculated spinstress coupling coefficients in the ground state, the ODMR contrast and basic sample dependent data (such as coherence time of the NV center’s spin), the sensitivity of strain quantum sensors can be deduced based on Hahn echo measurements at a given direction of the applied stress [248]. Next, we turn to the spinstrain interaction in the optically allowed excited state.
The ^{3}E excited state can be written in the {A_{1}, A_{2}, E_{x}, E_{y}, E_{1}, E_{2}} manifold, and the strain will affect these as [24], [25]
where
In Ref. [24], DFT PBE simulation on molecular cluster model was applied together with electric fields to estimate the d_{z} and d_{⊥} coefficients. The effect of the electric field was separated to “ionic” effect, i.e. piezo effect, when strain appears due to the presence of electric field, and “electronic” effect, when the cloud of the electrons and the corresponding potential changes due to the applied electric fields (see Ref. [24] and Supplementary Materials in Ref. [235]). In these simulations, the change in the geometry and splitting of the degenerate orbitals as a function of the applied electric field can be used to derive the order of magnitude for d_{z} and d_{⊥}. It was found that hd_{⊥}≈6.6 GHz (MV/m)^{−1} and hd_{⊥}≈0.6 GHz (MV/m)^{−1} (see Appendices C and D in Ref. [24]).
As mentioned above, the electric field is intertwined with strain via the piezo effect; thus, electric field may be applied to suppress the effect of strain [cf. Eqs. (32) and (33)] and tune two individual NV centers’ ZPL energies into the same position [39], [40]. It was also found that if a symmetrybreaking strain is present and a symmetrybreaking electric field is applied to the NV center at the same time, then the optical transition will show quadratic effect as a function of the strength of the electric field with anticrossing feature approaching the zero electric field; however, linear behavior shows up at very low strain fields [24], [38]. This explains the early measurements on the electric field dependence of the ZPL lines where single NV centers showed often linear characteristics upon the applied electric fields, but also quadratic features were observed for other single NV centers [243]. It was also shown that the optical transition from the A_{2} state toward the ms=±1 ground state will preserve the circular polarization at low strain fields, which was successfully used in entanglement schemes [34]. Here, we note that, as we have already discussed above, phonons can couple A_{1} state into A_{2}. This coupling is very small at cryogenic temperature [105] but it becomes substantial at T=10 K and above as can be inferred from the deduced ISC rates from the ms=±1 ^{3}E states toward the ^{1}A_{1} state [104]. Thus, the aforementioned conclusion about the strain dependence of the circular polarization in the optical transition between the A_{2} state and the ms=±1 ^{3}A_{2} state is valid at very low temperatures.
The temperature dependence of the lifetime of the
First, we discuss the basic spin Hamiltonian of the ^{3}E excited state. The orthorombic component of the D constant appears in the presence of the perpendicular component of the strain
where D_{⊥}=(D_{xx}−D_{yy})/2 and R(T)=(1−e^{−}^{hε}^{⊥/}^{k}^{B}^{T})/(1+e^{−}^{hε}^{⊥/}^{k}^{B}^{T}) is the temperature reduction factor. In the limit of zero strain, D_{⊥} disappears. At finite nonzero strain, D_{⊥} has the largest value in the limit of T=0 K and approaches zero in the limit of k_{B}T≫hε_{⊥}. The temperature dependence of the D component in the ^{3}E state is similar to that in the ^{3}A_{2} ground state (e.g. Ref. [250]), which will be discussed below. Similar to D_{⊥}, the spinorbit splitting in the ^{3}E state is reduced by increasing the temperature due to the Raman scattering of acoustic E phonons that couples E_{+} and E_{−} states [250]. This model was applied to understand the temperature dependence of ODMR and ZPL linewidths [251], and they found consistent results for the observed temperature dependence of the ZPL linewidth [180] with assuming quadratic interaction with acoustic A_{1} phonons in the optical transition [251], in contrast to the pure DJT model within a single effective phonon solution [183]. Here, we note that the DW factor was found also temperature dependent in this study [250], which may be used for alloptical thermometry with the NV center in nanodiamonds.
The aforementioned mechanisms are responsible for the temperature dependence of ISC rates between the ms=±1 states and the ^{1}A_{1} state [104]. We note that all these findings are based on the electronic solution of the ^{3}E state. In contrast, the strong electronphonon coupling due to DJT results in already strongly coupled vibronic solution even at T=0 K in the ^{3}E state. This implies that the phonon scattering between the vibronic levels and states should be considered, which can result in additional scattering channels, e.g. the participation of acoustic A_{1} phonons in the scattering between the
The temperature dependence of the Dconstant in the ground state can neither be described by a simple model in which the change in the Dconstant is attributed to the lattice constants increase at elevated temperatures because this effect is minor in diamond as confirmed in ab initio simulations [200]. Rather, dynamical effects with electronphonon coupling should be responsible for this effect as explained by a phenomenological approach in Ref. [244]. All these findings imply that the ab initio calculation of the temperature dependence of the magnetooptical properties in quasidegenerate states requires the calculation of single phonon emission and absorption, Raman process, or even more complicated scattering processes between the corresponding orbitals and spin states.
A recent example for the ab initio simulation of the phonon scattering due to electronphonon coupling between twolevel system is the calculation of electron spin flipping time (T_{1}) or rate (1/T_{1}) in Ref. [252]. It was found previously in experiments [253] that the rate goes as ~T or as a constant at 4<T<50 K depending on the spin density in the diamond sample and the applied magnetic fields, as ~exp(1/T) at T<50<200 K and as ~T^{5} at T>200 K (see Refs. [253], [254], [255], [256]). The different temperature regimes can be identified as single phonon absorption and emission, Orbachtype process, and secondorder Raman scattering process, respectively [253], [255]. In highpurity diamond samples with no external magnetic fields, the T_{1} exceeded 8 h for NV electron spin at 15 mK temperature [256], which is determined by phononic vacuum fluctuations.
Ab initio simulations computed T_{1} time of the ^{3}A_{2} ground state spin at the very low temperature regime where the single phonon absorption and emission is relevant [252]. At zero perturbation with perfect C_{3}_{v} symmetry, the system can be described as a twolevel system where the upper ms=±1 level is double degenerate with an energy gap of the Dconstant, which separates it from the lower ms=0 level. The Dtensor depends on the distance of the electron spins localized around the vacancy. During the vibration of atoms, the distance between the atoms dynamically change so the Dtensor dynamically changes. This effect may result in a spin flip because of the spin raising and lowering operators that appear in the dipolar electron spinelectron spin operator (e.g. S_{i}S_{j} in Eq. (8) and Ref. [257]). In the ab initio simulation, the derivative of the Dtensor elements is computed as a function of the ion coordinates along the normal coordinates of the phonon modes for each phonon mode, which effectively changes the distance between the electron spins localized on the carbon dangling bonds in the NV center [252]. To accelerate the calculation of Dtensor components, the total spin density was maximally localized on the e_{x,y} orbitals by Wannier localization [258]. The final rate at T → 0 K limit depends on the phonon density of states that may depend on the type of defects and defects’ concentrations, which shift the phonon density of states toward lower energies [252]. The final calculated Γ_{0} rates fall between 2×10^{−5} and 3×10^{−5} s^{−1} depending on the choice of the phonon density states, which are close to the experimental one at 3.47(16) 10^{−5} s^{−1} in Ref. [256]. The temperature dependence can be written as Γ_{0}×{1+3/[exp(ħD/k_{B}T)–1]}, which goes linearly with temperature because ħD is only 138 mK.
The ab initio theory on NV quantum bit in diamond is close to complete, but some issues have to be solved. For the full description of the ISC rates, the APES should be calculated also for the highly correlated singlet states. This requires to extend the CIcRPA or similar methods for calculating the quantum mechanical forces acting on the ions. This may reveal the issue of the critical spinorbit matrix element associated with the perpendicular component of the spinorbit coupling between the triplet and singlet states. An accurate approach in the calculation of spin density matrix might provide consistent dipolar electron spinelectron spin coupling parameters, which can be important for improving the calculation of the Dtensor components in the electronic excited state.
In general, the calculation of temperature dependence of the magnetooptical properties should still be further elaborated. Although the JahnTeller theory developed by Bersurker and Ham can be successfully applied to calculate the effective electronphonon coupling for degenerate states, one has to go beyond the single effective phonon approach for studying the accurate temperature dependence together with the calculation of Raman scattering process between the resultant states, an analog to the calculations in Ref. [252]. One possible route along this line is the application of the manybody perturbation theory of the electronphonon coupling [259], [260] that was successfully applied to the simulation of the ionization potential of diamondoids [261].
We further note that the spinorbit energy is much smaller than the electronphonon energy for the NV center in diamond; thus, it can be treated as a perturbation. In other defect quantum bits, the spinorbit and electronphonon energies are in the same order of magnitude; thus, the spinorbit and electronphonon coupling should be solved simultaneously [209]. This is highly nontrivial with going beyond an effective single phonon approach in the electronphonon calculations. Spinorbit interaction may contribute to the ZFS either in the first order (e.g. Ref. [209]) or the second order, which should be considered for defects with relatively large spinorbit energies.
Despite the persisting issues raised above, recent progress in the description of solidstate defect quantum bits has significantly contributed to the optimization of quantum bit operation and prediction of key properties that could be harnessed in experiments. We note that the full description of quantum bit cannot be separated from the description of the environment, which means the interaction with other defects’ charge and spin that are present either in bulk or close to the surface, beside other external perturbations such as strain, electric and magnetic fields, and temperature. Thus, the description of charge and spin dynamics caused by these species should be connected to the full characterization of the target defect quantum bit, to compute the readout contrasts, the longitudinal spin relaxation time, and the coherence time. By computing these key quantities, the sensitivity of quantum sensing protocols and optimization of quantum control can be designed and may guide future experimental studies. Only this comprehensive approach makes the ab initio search for alternative solidstate defect quantum bits reliable and powerful that might be superior for a given quantum technology application.
In this paper, we reviewed the ab initio theory on the NV center in diamond, which is an exemplary solidstate point defect quantum bit. We summarized the methods to calculate the key magnetooptical properties of this quantum bit and derived the decay rates and coupling parameters to different types of external perturbation and temperature. We briefly mentioned the missing pieces for the complete ab initio description of the operation of this quantum bit. The complete description of solidstate quantum bits would make possible to convert the phenomenological description of the spin control and quantum sensing to fully ab initio solution or, at least, a solution based on ab initio wave functions, densities, and density matrices.
A.G. acknowledges the National Excellence Program of QuantumCoherent Materials Project (Hungarian NKFIH grant no. KKP129866, Funder Id: http://dx.doi.org/10.13039/501100003549), the EU QuantERA QMagine Project (grant no. 127889, Funder Id: http://dx.doi.org/10.13039/501100003549), the EU H2020 Quantum Technology flagship project ASTERIQS (grant no. 820394, Funder Id: http://dx.doi.org/10.13039/100010661), and the National Quantum Technology Program (grant No. 20171.2.1NKP201700001, Funder Id: http://dx.doi.org/10.13039/501100012550).
[1] du Preez L. Ph.D. dissertation. University of Witwatersrand, 1965. Search in Google Scholar
[2] Goss JP, Jones R, Breuer SJ, Briddon PR, Öberg S. The twelveline 1.682 eV luminescence center in diamond and the vacancysilicon complex. Phys Rev Lett 1996;77:3041–4.1006211610.1103/PhysRevLett.77.3041 Search in Google Scholar
[3] Gali A, Fyta M, Kaxiras E. Ab initio supercell calculations on nitrogenvacancy center in diamond: electronic structure and hyperfine tensors. Phys Rev B 2008;77:155206.10.1103/PhysRevB.77.155206 Search in Google Scholar
[4] Larsson JA, Delaney P. Electronic structure of the nitrogenvacancy center in diamond from firstprinciples theory. Phys Rev B 2008;77:165201.10.1103/PhysRevB.77.165201 Search in Google Scholar
[5] Davies G, Hamer MF. Optical studies of the 1.945 eV vibronic band in diamond. Proc R Soc Lond Ser A 1976;348:285–98.10.1098/rspa.1976.0039 Search in Google Scholar
[6] Loubser JHN, van Wyk JP. Diamond Research (London), Vol. 11–15. London: Industrial Diamond Information Bureau, 1977:4–7. Search in Google Scholar
[7] Loubser JHN, van Wyk JA. Electron spin resonance in the study of diamond. Rep Prog Phys 1978;41:1201.10.1088/00344885/41/8/002 Search in Google Scholar
[8] Felton S, Edmonds AM, Newton ME, et al. Hyperfine interaction in the ground state of the negatively charged nitrogen vacancy center in diamond. Phys Rev B 2009;79:075203.10.1103/PhysRevB.79.075203 Search in Google Scholar
[9] He XF, Manson NB, Fisk PTH. Paramagnetic resonance of photoexcited NV defects in diamond. II. Hyperfine interaction with the 14N nucleus. Phys Rev B 1993;47:8816.10.1103/PhysRevB.47.8816 Search in Google Scholar
[10] Pfender M, Aslam N, Simon P, et al. Protecting a diamond quantum memory by charge state control. Nano Lett 2017;17:5931–7.10.1021/acs.nanolett.7b0179628872881 Search in Google Scholar
[11] Smeltzer B, Childress L, Gali A. 13C hyperfine interactions in the nitrogenvacancy centre in diamond. N J Phys 2011;13:025021.10.1088/13672630/13/2/025021 Search in Google Scholar
[12] Ivády V, Gali A, Abrikosov IA. HybridDFT + Vw method for band structure calculation of semiconducting transition metal compounds: the case of cerium dioxide. J Phys Condens Matter 2017;29:454002.10.1088/1361648X/aa8b93 Search in Google Scholar
[13] Gruber A, Dräbenstedt A, Tietz C, Fleury L, Wrachtrup J, von Borczyskowski C. Scanning confocal optical microscopy and magnetic resonance on single defect centers. Science 1997;276:2012–4.10.1126/science.276.5321.2012 Search in Google Scholar
[14] Bourgeois E, Jarmola A, Siyushev P, et al. Photoelectric detection of electron spin resonance of nitrogenvacancy centres in diamond. Nat Commun 2015;6:8577.10.1038/ncomms957726486014 Search in Google Scholar
[15] Siyushev P, Nesladek M, Bourgeois E, et al. Photoelectrical imaging and coherent spinstate readout of single nitrogenvacancy centers in diamond. Science 2019;363:728–31.10.1126/science.aav278930765564 Search in Google Scholar
[16] Childress L, Gurudev Dutt MV, Taylor JM, et al. Coherent dynamics of coupled electron and nuclear spin qubits in diamond. Science 2006;314:281–5.10.1126/science.113187116973839 Search in Google Scholar
[17] Mita Y. Change of absorption spectra in typeIb diamond with heavy neutron irradiation. Phys Rev B 1996;53:11360.10.1103/PhysRevB.53.11360 Search in Google Scholar
[18] Gaebel T, Domhan M, Wittmann C, et al. Photochromism in single nitrogenvacancy defect in diamond. Appl Phys B 2006;82:243–6.10.1007/s0034000520562 Search in Google Scholar
[19] Aslam N, Waldherr G, Neumann P, Jelezko F, Wrachtrup J. Photoinduced ionization dynamics of the nitrogen vacancy defect in diamond investigated by singleshot charge state detection. N J Phys 2013;15:013064.10.1088/13672630/15/1/013064 Search in Google Scholar
[20] Siyushev P, Pinto H, Vörös M, Gali A, Jelezko F, Wrachtrup J. Optically controlled switching of the charge state of a single nitrogenvacancy center in diamond at cryogenic temperatures. Phys Rev Lett 2013;110:167402.2367963710.1103/PhysRevLett.110.167402 Search in Google Scholar
[21] Coulson CA, Kearsley MJ. Colour centres in irradiated diamonds. I. Proc R Soc Lond Ser A Math Phys Sci 1957;241:433–54.10.1098/rspa.1957.0138 Search in Google Scholar
[22] Lenef A, Rand SC. Electronic structure of the NV center in diamond: theory. Phys Rev B 1996;53:13441.10.1103/PhysRevB.53.13441 Search in Google Scholar
[23] Manson NB, Harrison JP, Sellars MJ. Nitrogenvacancy center in diamond: model of the electronic structure and associated dynamics. Phys Rev B 2006;74:104303.10.1103/PhysRevB.74.104303 Search in Google Scholar
[24] Maze JR, Gali A, Togan E, et al. Properties of nitrogenvacancy centers in diamond: the group theoretic approach. N J Phys 2011;13:025025.10.1088/13672630/13/2/025025 Search in Google Scholar
[25] Doherty MW, Manson NB, Delaney P, Hollenberg LCL. The negatively charged nitrogenvacancy centre in diamond: the electronic solution. N J Phys 2011;13:025019.10.1088/13672630/13/2/025019 Search in Google Scholar
[26] Jelezko F, Gaebel T, Popa I, Gruber A, Wrachtrup J. Observation of coherent oscillations in a single electron spin. Phys Rev Lett 2004;92:076401.10.1103/PhysRevLett.92.076401 Search in Google Scholar
[27] Toyli DM, Christle DJ, Alkauskas A, Buckley BB, Van de Walle CG, Awschalom DD. Measurement and control of single nitrogenvacancy center spins above 600 K. Phys Rev X 2012;2:031001. Search in Google Scholar
[28] Liu GQ, Feng X, Wang N, Li Q, Liu RB. Coherent quantum control of nitrogenvacancy center spins near 1000 kelvin. Nat Commun 2019;10:1344.3090297210.1038/s41467019093272 Search in Google Scholar
[29] Maze JR, Stanwix PL, Hodges JS, et al. Nanoscale magnetic sensing with an individual electronic spin in diamond. Nature 2008;455:644–7.1883327510.1038/nature07279 Search in Google Scholar
[30] Balasubramanian G, Neumann P, Twitchen D, et al. Ultralong spin coherence time in isotopically engineered diamond. Nat Mater 2009;8:383–7.1934997010.1038/nmat2420 Search in Google Scholar
[31] Jelezko F, Gaebel T, Popa I, Dunham M, Gruber A, Wrachtrup J. Observation of coherent oscillation of a single nuclear spin and realization of a twoqubit conditional quantum gate. Phys Rev Lett 2004;93:130501.10.1103/PhysRevLett.93.13050115524692 Search in Google Scholar
[32] Neumann P, Beck J, Steiner M, et al. Singleshot readout of a single nuclear spin. Science 2010;329:542–4.10.1126/science.118907520595582 Search in Google Scholar
[33] Robledo L, Childress L, Bernien H, Hensen B, Alkemade PFA, Hanson R. Highfidelity projective readout of a solidstate spin quantum register. Nature 2011;477:574–8.10.1038/nature1040121937989 Search in Google Scholar
[34] Togan E, Chu Y, Trifonov AS, et al. Quantum entanglement between an optical photon and a solidstate spin qubit. Nature 2010;466:730–4.10.1038/nature09256 Search in Google Scholar
[35] Yang S, Wang Y, Rao DDB, et al. Highfidelity transfer and storage of photon states in a single nuclear spin. Nat Photonics 2016;10:507–11.10.1038/nphoton.2016.103 Search in Google Scholar
[36] Hensen B, Bernien H, Dréau AE, et al. Loopholefree Bell inequality violation using electron spins separated by 1.3 kilometres. Nature 2015;526:682–6.10.1038/nature1575926503041 Search in Google Scholar
[37] DiVincenzo DP. The physical implementation of quantum computation. Fortsch Phys 2000;48:771–83.10.1002/15213978(200009)48:9/11<771::AIDPROP771>3.0.CO;2E Search in Google Scholar
[38] Tamarat P, Manson NB, Harrison JP, et al. Spinflip and spinconserving optical transitions of the nitrogenvacancy centre in diamond. N J Phys 2008;10:045004.10.1088/13672630/10/4/045004 Search in Google Scholar
[39] Bassett LC, Heremans FJ, Yale CG, Buckley BB, Awschalom DD. Electrical tuning of single nitrogenvacancy center optical transitions enhanced by photoinduced fields. Phys Rev Lett 2011;107:266403.2224317110.1103/PhysRevLett.107.266403 Search in Google Scholar
[40] Pfaff W, Hensen BJ, Bernien H, et al. Quantum information. Unconditional quantum teleportation between distant solidstate quantum bits. Science 2014;345:532–5.2508269610.1126/science.1253512 Search in Google Scholar
[41] Balasubramanian G, Chan IY, Kolesov R, et al. Nanoscale imaging magnetometry with diamond spins under ambient conditions. Nature 2008;455:648–51.1883327610.1038/nature07278 Search in Google Scholar
[42] Degen CL. Scanning magnetic field microscope with a diamond singlespin sensor. Appl Phys Lett 2008;92:243111.10.1063/1.2943282 Search in Google Scholar
[43] Dolde F, Fedder H, Doherty MW, et al. Electricfield sensing using single diamond spins. Nat Phys 2011;7:459.10.1038/nphys1969 Search in Google Scholar
[44] Teissier J, Barfuss A, Appel P, Neu E, Maletinsky P. Strain coupling of a nitrogenvacancy center spin to a diamond mechanical oscillator. Phys Rev Lett 2014;113:020503.10.1103/PhysRevLett.113.02050325062153 Search in Google Scholar
[45] Barfuss A, Teissier J, Neu E, Nunnenkamp A, Maletinsky P. Strong mechanical driving of a single electron spin. Nat Phys 2015;11:820–4.10.1038/nphys3411 Search in Google Scholar
[46] Balasubramanian G, Chan IY, Kolesov R, et al. Nanoscale imaging magnetometry with diamond spins under ambient conditions. Nature 2008;455:648–51.1883327610.1038/nature07278 Search in Google Scholar
[47] Degen CL. Scanning magnetic field microscope with a diamond singlespin sensor. Appl Phys Lett 2008;92:243111.10.1063/1.2943282 Search in Google Scholar
[48] Dolde F, Fedder H, Doherty MW, et al. Electricfield sensing using single diamond spins. Nat Phys 2011;7:459.10.1038/nphys1969 Search in Google Scholar
[49] Teissier J, Barfuss A, Appel P, Neu E, Maletinsky P. Strain coupling of a nitrogenvacancy center spin to a diamond mechanical oscillator. Phys Rev Lett 2014;113:020503.10.1103/PhysRevLett.113.02050325062153 Search in Google Scholar
[50] Barfuss A, Teissier J, Neu E, Nunnenkamp A, Maletinsky P. Strong mechanical driving of a single electron spin. Nat Phys 2015;11:820–4.10.1038/nphys3411 Search in Google Scholar
[51] MacQuarrie ER, Gosavi TA, Bhave SA, Fuchs GD. Continuous dynamical decoupling of a single diamond nitrogenvacancy center spin with a mechanical resonator. Phys Rev B 2015;92:224419.10.1103/PhysRevB.92.224419 Search in Google Scholar
[52] Golter DA, Oo T, Amezcua M, Stewart KA, Wang H. Optomechanical quantum control of a nitrogenvacancy center in diamond. Phys Rev Lett 2016;116:143602.2710470910.1103/PhysRevLett.116.143602 Search in Google Scholar
[53] Kucsko G, Maurer PC, Yao NY, et al. Nanometrescale thermometry in a living cell. Nature 2013;500:54–8.10.1038/nature12373 Search in Google Scholar
[54] Toyli DM, de las Casas CF, Christle DJ, Dobrovitski VV, Awschalom DD. Fluorescence thermometry enhanced by the quantum coherence of single spins in diamond. Proc Natl Acad Sci USA 2013;110:8417–21.10.1073/pnas.1306825110 Search in Google Scholar
[55] Neumann P, Jakobi I, Dolde F, et al. Highprecision nanoscale temperature sensing using single defects in diamond. Nano Lett 2013;13:2738–42.10.1021/nl401216y23721106 Search in Google Scholar
[56] Staudacher T, Shi F, Pezzagna S, et al. Nuclear magnetic resonance spectroscopy on a (5nanometer)3 Sample volume. Science 2013;339:561–3.2337200910.1126/science.1231675 Search in Google Scholar
[57] Mamin HJ, Kim M, Sherwood MH, et al. Nanoscale nuclear magnetic resonance with a nitrogenvacancy spin sensor. Science 2013;339:557–60.10.1126/science.123154023372008 Search in Google Scholar
[58] DeVience SJ, Pham LM, Lovchinsky I, et al. Nanoscale NMR spectroscopy and imaging of multiple nuclear species. Nat Nanotechnol 2015;10:129–34.10.1038/nnano.2014.31325559712 Search in Google Scholar
[59] Häberle T, SchmidLorch D, Reinhard F, Wrachtrup J. Nanoscale nuclear magnetic imaging with chemical contrast. Nat Nanotechnol 2015;10:125–8.10.1038/nnano.2014.29925559711 Search in Google Scholar
[60] Rugar D, Mamin HJ, Sherwood MH, et al. Proton magnetic resonance imaging using a nitrogenvacancy spin sensor. Nat Nanotechnol 2015;10:120–4.10.1038/nnano.2014.28825531089 Search in Google Scholar
[61] Boss JM, Cujia KS, Zopes J, Degen CL. Quantum sensing with arbitrary frequency resolution. Science 2017;356: 837–40.10.1126/science.aam700928546209 Search in Google Scholar
[62] Schmitt S, Gefen T, Störner FM, et al. Submillihertz magnetic spectroscopy performed with a nanoscale quantum sensor. Science 2017;356:832–7.10.1126/science.aam553228546208 Search in Google Scholar
[63] Aslam N, Pfender M, Neumann P, et al. Nanoscale nuclear magnetic resonance with chemical resolution. Science 2017;357:67–71.2857245310.1126/science.aam8697 Search in Google Scholar
[64] Glenn DR, Bucher DB, Lee J, et al. Highresolution magnetic resonance spectroscopy using a solidstate spin sensor. Nature 2018;555:351–4.10.1038/nature2578129542693 Search in Google Scholar
[65] Jacques V, Neumann P, Beck J, et al. Dynamic polarization of single nuclear spins by optical pumping of nitrogenvacancy color centers in diamond at room temperature. Phys Rev Lett 2009;102:057403.10.1103/PhysRevLett.102.05740319257552 Search in Google Scholar
[66] Gali A. Identification of individual 13C isotopes of nitrogenvacancy center in diamond by combining the polarization studies of nuclear spins and firstprinciples calculations. Phys Rev B 2009;80:241204.10.1103/PhysRevB.80.241204 Search in Google Scholar
[67] Ivády V, Szász K, Falk AL, et al. Theoretical model of dynamic spin polarization of nuclei coupled to paramagnetic point defects in diamond and silicon carbide. Phys Rev B 2015;92:115206.10.1103/PhysRevB.92.115206 Search in Google Scholar
[68] Chen Q, Schwarz I, Jelezko F, Retzker A, Plenio MB. Optical hyperpolarization of 13C nuclear spins in nanodiamond ensembles. Phys Rev B 2015;92:184420.10.1103/PhysRevB.92.184420 Search in Google Scholar
[69] Álvarez GA, Bretschneider CO, Fischer R, et al. Local and bulk 13C hyperpolarization in nitrogenvacancycentred diamonds at variable fields and orientations. Nat Commun 2015;6:8456.10.1038/ncomms9456 Search in Google Scholar
[70] King JP, Jeong K, Vassiliou CC, et al. Roomtemperature in situ nuclear spin hyperpolarization from optically pumped nitrogen vacancy centres in diamond. Nat Commun 2015;6:8965.10.1038/ncomms996526639147 Search in Google Scholar
[71] Scheuer J, Schwartz I, Chen Q, et al. Optically induced dynamic nuclear spin polarisation in diamond. N J Phys 2016;18:013040.10.1088/13672630/18/1/013040 Search in Google Scholar
[72] Wunderlich R, Kohlrautz J, Abel B, Haase J, Meijer J. Optically induced cross relaxation via nitrogenrelated defects for bulk diamond 13C hyperpolarization. Phys Rev B 2017;96:220407.10.1103/PhysRevB.96.220407 Search in Google Scholar
[73] Ajoy A, Liu K, Nazaryan R, et al. Orientationindependent room temperature optical 13C hyperpolarization in powdered diamond. Sci Adv 2018;4:eaar5492.2979578310.1126/sciadv.aar5492 Search in Google Scholar
[74] Schwartz I, Scheuer J, Tratzmiller B, et al. Robust optical polarization of nuclear spin baths using Hamiltonian engineering of nitrogenvacancy center quantum dynamics. Sci Adv 2018;4:eaat8978.3018206010.1126/sciadv.aat8978 Search in Google Scholar
[75] Abrams D, Trusheim ME, Englund DR, Shattuck MD, Meriles CA. Dynamic nuclear spin polarization of liquids and gases in contact with nanostructured diamond. Nano Lett 2014;14:2471–8.2475475510.1021/nl500147b Search in Google Scholar
[76] FernándezAcebal P, Rosolio O, Scheuer J, et al. Toward hyperpolarization of oil molecules via single nitrogen vacancy centers in diamond. Nano Lett 2018;18:1882–7.2947008910.1021/acs.nanolett.7b05175 Search in Google Scholar
[77] Cai J, Retzker A, Jelezko F, Plenio MB. Towards a largescale quantum simulator on diamond surface at room temperature. Nat Phys 2013;9:168.10.1038/nphys2519 Search in Google Scholar
[78] Chou JP, Retzker A, Gali A. Nitrogenterminated diamond (111) surface for roomtemperature quantum sensing and simulation. Nano Lett 2017;17:2294–8.10.1021/acs.nanolett.6b0502328339209 Search in Google Scholar
[79] Meijer J, Burchard B, Domhan M, et al. Generation of single color centers by focused nitrogen implantation. Appl Phys Lett 2005;87:261909.10.1063/1.2103389 Search in Google Scholar
[80] Rabeau JR, Reichart P, Tamanyan G, et al. Implantation of labelled single nitrogen vacancy centers in diamond using 15N. Appl Phys Lett 2006;88:023113.10.1063/1.2158700 Search in Google Scholar
[81] Pezzagna S, Naydenov B, Jelezko F, Wrachtrup J, Meijer J. Creation efficiency of nitrogenvacancy centres in diamond. N J Phys 2010;12:065017.10.1088/13672630/12/6/065017 Search in Google Scholar
[82] OforiOkai BK, Pezzagna S, Chang K, et al. Spin properties of very shallow nitrogen vacancy defects in diamond. Phys Rev B 2012;86:081406.10.1103/PhysRevB.86.081406 Search in Google Scholar
[83] Van Dam SB, Walsh M, Degen MJ, et al. Optical coherence of diamond nitrogenvacancy centers formed by ion implantation and annealing. Phys Rev B 2019;99:161203.10.1103/PhysRevB.99.161203 Search in Google Scholar
[84] Smith WV, Sorokin PP, Gelles IL, Lasher GJ. Electronspin resonance of nitrogen donors in diamond. Phys Rev 1959;115:1546.10.1103/PhysRev.115.1546 Search in Google Scholar
[85] Loubser JHN, Preez LD. New lines in the electron spin resonance spectrum of substitutional nitrogen donors in diamond. Br J Appl Phys 1965;16:457.10.1088/05083443/16/4/307 Search in Google Scholar
[86] Farrer RG. On the substitutional nitrogen donor in diamond. Solid State Commun 1969;7:685–8.10.1016/00381098(69)905936 Search in Google Scholar
[87] Ohno K, Heremans FJ, Bassett LC, et al. Engineering shallow spins in diamond with nitrogen deltadoping. Appl Phys Lett 2012;101:082413.10.1063/1.4748280 Search in Google Scholar
[88] Collins AT, Kiflawi I. The annealing of radiation damage in type Ia diamond. J Phys Condensed Matter 2009;21:364209.10.1088/09538984/21/36/364209 Search in Google Scholar
[89] Acosta VM, Bauch E, Ledbetter MP, et al. Diamonds with a high density of nitrogenvacancy centers for magnetometry applications. Phys Rev B 2009;80:115202.10.1103/PhysRevB.80.115202 Search in Google Scholar
[90] Waldermann F, Olivero P, Nunn J, et al. Creating diamond color centers for quantum optical applications. Diamond Relat Mater 2007;16:1887–95.10.1016/j.diamond.2007.09.009 Search in Google Scholar
[91] Huang Z, Li WD, Santori C, et al. Diamond nitrogenvacancy centers created by scanning focused helium ion beam and annealing. Appl Phys Lett 2013;103:081906.10.1063/1.4819339 Search in Google Scholar
[92] Mainwood A. Nitrogen and nitrogenvacancy complexes and their formation in diamond. Phys Rev B 1994;49:7934–40.10.1103/PhysRevB.49.7934 Search in Google Scholar
[93] Capelli M, Heffernan A, Ohshima T, et al. Increased nitrogenvacancy centre creation yield in diamond through electron beam irradiation at high temperature. Carbon 2019;143:714–9.10.1016/j.carbon.2018.11.051 Search in Google Scholar
[94] Fávaro de Oliveira F, Antonov D, Wang Y, et al. Tailoring spin defects in diamond by lattice charging. Nat Commun 2017;8:15409.2851358110.1038/ncomms15409 Search in Google Scholar
[95] Bourgeois E, Londero E, Buczak K, et al. Enhanced photoelectric detection of NV magnetic resonances in diamond under dualbeam excitation. Phys Rev B 2017;95:041402.10.1103/PhysRevB.95.041402 Search in Google Scholar
[96] Hrubesch FM, Braunbeck G, Stutzmann M, Reinhard F, Brandt MS. Efficient electrical spin readout of NV− centers in diamond. Phys Rev Lett 2017;118:037601.10.1103/PhysRevLett.118.03760128157351 Search in Google Scholar
[97] Gulka M, Bourgeois E, Hruby J, et al. Pulsed photoelectric coherent manipulation and detection of N−V center spins in diamond. Phys Rev Appl 2017;7:044032.10.1103/PhysRevApplied.7.044032 Search in Google Scholar
[98] Fuchs GD, Burkard G, Klimov PV, Awschalom DD. A quantum memory intrinsic to single nitrogen–vacancy centres in diamond. Nat Phys 2011;7:789–93.10.1038/nphys2026 Search in Google Scholar
[99] Maurer PC, Kucsko G, Latta C, et al. Roomtemperature quantum bit memory exceeding one second. Science 2012;336:1283–6.10.1126/science.122051322679092 Search in Google Scholar
[100] Waldherr G, Wang Y, Zaiser S, et al. Quantum error correction in a solidstate hybrid spin register. Nature 2014;506:204–7.10.1038/nature12919 Search in Google Scholar
[101] Abobeih MH, Cramer J, Bakker MA, et al. Onesecond coherence for a single electron spin coupled to a multiqubit nuclearspin environment. Nat Commun 2018;9:2552.2995932610.1038/s4146701804916z Search in Google Scholar
[102] Waldherr G, Beck J, Steiner M, et al. Dark states of single nitrogenvacancy centers in diamond unraveled by single shot NMR. Phys Rev Lett 2011;106:157601.10.1103/PhysRevLett.106.15760121568615 Search in Google Scholar
[103] Batalov A, Jacques V, Kaiser F, et al. Low temperature studies of the excitedstate structure of negatively charged nitrogenvacancy color centers in diamond. Phys Rev Lett 2009;102:195506.1951897310.1103/PhysRevLett.102.195506 Search in Google Scholar
[104] Goldman ML, Sipahigil A, Doherty MW, et al. Phononinduced population dynamics and intersystem crossing in nitrogenvacancy centers. Phys Rev Lett 2015;114:145502.10.1103/PhysRevLett.114.14550225910136 Search in Google Scholar
[105] Thiering G, Gali A. Ab initio calculation of spinorbit coupling for an NV center in diamond exhibiting dynamic JahnTeller effect. Phys Rev B 2017;96:081115.10.1103/PhysRevB.96.081115 Search in Google Scholar
[106] Evarestov RA. Use of representative points of the Brillouin zone for the selfconsistent calculations of solids in the large unit cell approach. Phys Stat Solidi B 1975;72:569–78.10.1002/pssb.2220720213 Search in Google Scholar
[107] Monkhorst HJ, Pack JK. Special points for Brillouinzone integrations. Phys Rev B 1976;13:5188.10.1103/PhysRevB.13.5188 Search in Google Scholar
[108] Zyubin AS, Mebel AM, Hayashi M, Chang HC, Lin SH. Quantum chemical modeling of photoadsorption properties of the nitrogenvacancy point defect in diamond. J Comput Chem 2009;30:119–31.10.1002/jcc.2104218548526 Search in Google Scholar
[109] Delaney P, Greer JC, Larsson JA. Spinpolarization mechanisms of the nitrogenvacancy center in diamond. Nano Lett 2010;10:610–4.2008527110.1021/nl903646p Search in Google Scholar
[110] Kohn W, Sham LJ. Selfconsistent equations including exchange and correlation effects. Phys Rev Lett 1965;140:A1133. Search in Google Scholar
[111] Parr RG, Yang W. Densityfunctional theory of atoms and molecules. Cambridge: Oxford University Press, 1989. Search in Google Scholar
[112] Ceperley DM, Alder BJ. Ground state of the electron gas by a stochastic method. Phys Rev Lett 1980;45:566.10.1103/PhysRevLett.45.566 Search in Google Scholar
[113] Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys Rev Lett 1996;77: 3865–8.10.1103/PhysRevLett.77.386510062328 Search in Google Scholar
[114] Blöchl PE. Projector augmentedwave method. Phys Rev B 1994;50:17953.10.1103/PhysRevB.50.17953 Search in Google Scholar
[115] Blöchl PE. Firstprinciples calculations of defects in oxygendeficient silica exposed to hydrogen. Phys Rev B 2000;62:6158.10.1103/PhysRevB.62.6158 Search in Google Scholar
[116] Shockley W. On the surface states associated with a periodic potential. Phys Rev 1939;56:317.10.1103/PhysRev.56.317 Search in Google Scholar
[117] Cole MW, Cohen MH. Imagepotentialinduced surface bands in insulators. Phys Rev Lett 1969;23:1238.10.1103/PhysRevLett.23.1238 Search in Google Scholar
[118] Mulliken RS. The Rydberg states of molecules.1a Parts IV1b. J Am Chem Soc 1964;86:3183–97.10.1021/ja01070a001 Search in Google Scholar
[119] Vörös M, Gali A. Optical absorption of diamond nanocrystals from ab initio densityfunctional calculations. Phys Rev B 2009;80:161411.10.1103/PhysRevB.80.161411 Search in Google Scholar
[120] Kaviani M, Deák P, Aradi B, Frauenheim T, Chou JP, Gali A. Proper surface termination for luminescent nearsurface NV centers in diamond. Nano Lett 2014;14:4772–7.2505462110.1021/nl501927y Search in Google Scholar
[121] Boys SF, Egerton AC. Opening remarks. Proc R Soc Lond Ser A Math Phys Sci 1950;200:542. Search in Google Scholar
[122] Nizovtsev AP, Kilin SY, Pushkarchuk AL, et al. Nonflipping 13C spins near an NV center in diamond: hyperfine and spatial characteristics by density functional theory simulation of the C510[NV]H252 cluster. N J Phys 2018;20:023022.10.1088/13672630/aaa910 Search in Google Scholar
[123] Northrup JE, Felice RD, Neugebauer J. Energetics of H and NH2 on GaN(1010) and implications for the origin of nanopipe defects. Phys Rev B 1996;56:R4325. Search in Google Scholar
[124] Deák P, Aradi B, Kaviani M, Frauenheim T, Gali A. Formation of NV centers in diamond: a theoretical study based on calculated transitions and migration of nitrogen and vacancy related defects. Phys Rev B 2014;89:075203.10.1103/PhysRevB.89.075203 Search in Google Scholar
[125] Makov G, Payne MC. Periodic boundary conditions in ab initio calculations. Phys Rev B 1995;51:4014–22.10.1103/PhysRevB.51.4014 Search in Google Scholar
[126] Lany S, Zunger A. Assessment of correction methods for the bandgap problem and for finitesize effects in supercell defect calculations: case studies for ZnO and GaAs. Phys Rev B 2008;78:235104.10.1103/PhysRevB.78.235104 Search in Google Scholar
[127] Freysoldt C, Neugebauer J, Van de Walle CG. Fully ab initio finitesize corrections for chargeddefect supercell calculations. Phys Rev Lett 2009;102:016402.10.1103/PhysRevLett.102.01640219257218 Search in Google Scholar
[128] Komsa HP, Rantala TT, Pasquarello A. Finitesize supercell correction schemes for charged defect calculations. Phys Rev B 2012;86:045112.10.1103/PhysRevB.86.045112 Search in Google Scholar
[129] Deák P, Aradi B, Frauenheim T, Janzén E, Gali A. Accurate defect levels obtained from the HSE06 rangeseparated hybrid functional. Phys Rev B 2010;81:153203.10.1103/PhysRevB.81.153203 Search in Google Scholar
[130] Krukau AV, Vydrov OA, Izmaylov AF, Scuseria GE. Influence of the exchange screening parameter on the performance of screened hybrid functionals. J Chem Phys 2006;125:224106.10.1063/1.240466317176133 Search in Google Scholar
[131] Weber JR, Koehl WF, Varley JB, et al. Quantum computing with defects. Proc Natl Acad Sci USA 2010;107:8513–8.10.1073/pnas.1003052107 Search in Google Scholar
[132] Dhomkar S, Jayakumar H, Zangara PR, Meriles CA. Charge dynamics in nearsurface, variabledensity ensembles of nitrogenvacancy centers in diamond. Nano Lett 2018;18:4046.2973361610.1021/acs.nanolett.8b01739 Search in Google Scholar
[133] Maier F, Ristein J, Ley L. Electron affinity of plasmahydrogenated and chemically oxidized diamond (100) surfaces. Phys Rev B 2001;64:165411.10.1103/PhysRevB.64.165411 Search in Google Scholar
[134] Grotz B, Hauf MV, Dankerl M, et al. Charge state manipulation of qubits in diamond. Nat Commun 2012;3:729.2239562010.1038/ncomms1729 Search in Google Scholar
[135] Fu KMC, Santori C, Barclay PE, Beausoleil RG. Conversion of neutral nitrogenvacancy centers to negatively charged nitrogenvacancy centers through selective oxidation. Appl Phys Lett 2010;96:121907.10.1063/1.3364135 Search in Google Scholar
[136] Rondin L, Dantelle G, Slablab A, et al. Surfaceinduced charge state conversion of nitrogenvacancy defects in nanodiamonds. Phys Rev B 2010;82:115449.10.1103/PhysRevB.82.115449 Search in Google Scholar
[137] Hauf MV, Grotz B, Naydenov B, et al. Chemical control of the charge state of nitrogenvacancy centers in diamond. Phys Rev B 2011;83:081304.10.1103/PhysRevB.83.081304 Search in Google Scholar
[138] Chou JP, Gali A. Nitrogenvacancy diamond sensor: novel diamond surfaces from ab initio simulations. MRS Commun 2017;7:551–62.10.1557/mrc.2017.75 Search in Google Scholar
[139] Li S, Chou JP, Wei J, Sun M, Hu A, Gali A. Oxygenated (113) diamond surface for nitrogenvacancy quantum sensors with preferential alignment and long coherence time from first principles. Carbon 2019;145:273–80.10.1016/j.carbon.2019.01.016 Search in Google Scholar
[140] Tiwari AK, Goss JP, Briddon PR, et al. Calculated electron affinity and stability of halogenterminated diamond. Phys Rev B 2011;84:245305.10.1103/PhysRevB.84.245305 Search in Google Scholar
[141] Stacey A, O’Donnell KM, Chou JP, et al. Nitrogen terminated diamond. Adv Mater Interfaces 2015;2:1500079.10.1002/admi.201500079 Search in Google Scholar
[142] Kawai S, Yamano H, Sonoda T, et al. Nitrogenterminated diamond surface for nanoscale NMR by shallow nitrogenvacancy centers. J Phys Chem C 2019;123:3594–604.10.1021/acs.jpcc.8b11274 Search in Google Scholar
[143] Edmonds AM, D’HaenensJohansson UFS, Cruddace RJ, et al. Production of oriented nitrogenvacancy color centers in synthetic diamond. Phys Rev B 2012;86:035201.10.1103/PhysRevB.86.035201 Search in Google Scholar
[144] Michl J, Teraji T, Zaiser S, et al. Perfect alignment and preferential orientation of nitrogenvacancy centers during chemical vapor deposition diamond growth on (111) surfaces. Appl Phys Lett 2014;104:102407.10.1063/1.4868128 Search in Google Scholar
[145] Lesik M, Tetienne JP, Tallaire A, et al. Perfect preferential orientation of nitrogenvacancy defects in a synthetic diamond sample. Appl Phys Lett 2014;104:113107.10.1063/1.4869103 Search in Google Scholar
[146] Tahara K, Ozawa H, Iwasaki T, Mizuochi N, Hatano M. Quantifying selective alignment of ensemble nitrogenvacancy centers in (111) diamond. Appl Phys Lett 2015;107:193110.10.1063/1.4935709 Search in Google Scholar
[147] Ozawa H, Tahara K, Ishiwata H, Hatano M, Iwasaki T. Formation of perfectly aligned nitrogenvacancycenter ensembles in chemicalvapordepositiongrown diamond (111). Appl Phys Express 2017;10:045501.10.7567/APEX.10.045501 Search in Google Scholar
[148] Osterkamp C, Mangold M, Lang J, et al. Engineering preferentiallyaligned nitrogenvacancy centre ensembles in CVD grown diamond. Sci Rep 2019;9:5786.3096248010.1038/s41598019423147 Search in Google Scholar
[149] Lesik M, Plays T, Tallaire A, et al. Preferential orientation of NV defects in CVD diamond films grown on (113)oriented substrates. Diamond Relat Mater 2015;56:47–53.10.1016/j.diamond.2015.05.003 Search in Google Scholar
[150] Chouaieb S, Martínez LJ, Akhtar W, et al. Optimizing synthetic diamond samples for quantum sensing technologies by tuning the growth temperature. Diamond Relat Mater 2019;96:85–9.10.1016/j.diamond.2019.04.022 Search in Google Scholar
[151] Karin T, Dunham S, Fu KM. Alignment of the diamond nitrogen vacancy center by strain engineering. Appl Phys Lett 2014;105:053106.10.1063/1.4892544 Search in Google Scholar
[152] Atumi MK, Goss JP, Briddon PR, Rayson MJ. Atomistic modeling of the polarization of nitrogen centers in diamond due to growth surface orientation. Phys Rev B 2013;88:245301.10.1103/PhysRevB.88.245301 Search in Google Scholar
[153] Miyazaki T, Miyamoto Y, Makino T, et al. Atomistic mechanism of perfect alignment of nitrogenvacancy centers in diamond. Appl Phys Lett 2014;105:261601.10.1063/1.4904988 Search in Google Scholar
[154] Stacey A, Dontschuk N, Chou JP, et al. Evidence for primal sp2 defects at the diamond surface. Adv Mater Interfaces 2019;6:1801449.10.1002/admi.201801449 Search in Google Scholar
[155] Chou JP, Bodrog Z, Gali A. Firstprinciples study of charge diffusion between proximate solidstate qubits and its implications on sensor applications. Phys Rev Lett 2018;120:136401.10.1103/PhysRevLett.120.13640129694166 Search in Google Scholar
[156] Choi J, Choi S, Kucsko G, et al. Depolarization dynamics in a strongly interacting solidstate spin ensemble. Phys Rev Lett 2017;118:093601.10.1103/PhysRevLett.118.093601 Search in Google Scholar
[157] Bluvstein D, Zhang Z, Jayich ACB. Identifying and mitigating charge instabilities in shallow diamond nitrogenvacancy centers. Phys Rev Lett 2019;122:076101.10.1103/PhysRevLett.122.07610130848640 Search in Google Scholar
[158] Lozovoi AY, Alavi A, Kohanoff J, LyndenBell RM. Ab initio simulation of charged slabs at constant chemical potential. J Chem Phys 2001;115:1661.10.1063/1.1379327 Search in Google Scholar
[159] Lozovoi AY, Alavi A. Reconstruction of charged surfaces: general trends and a case study of Pt(110) and Au(110). Phys Rev B 2003;68:245416.10.1103/PhysRevB.68.245416 Search in Google Scholar
[160] Scivetti I, Persson M. The electrostatic interaction of an external charged system with a metal surface: a simplified density functional theory approach. J Phys Condensed Matter 2013;25:355006.10.1088/09538984/25/35/355006 Search in Google Scholar
[161] Wang D, Han D, Li XB, et al. Determination of formation and ionization energies of charged defects in twodimensional materials. Phys Rev Lett 2015;114:196801.2602418910.1103/PhysRevLett.114.196801 Search in Google Scholar
[162] Vinichenko D, Sensoy MG, Friend CM, Kaxiras E. Accurate formation energies of charged defects in solids: a systematic approach. Phys Rev B 2017;95:235310.10.1103/PhysRevB.95.235310 Search in Google Scholar
[163] Bal KM, Neyts EC. Modelling molecular adsorption on charged or polarized surfaces: a critical flaw in common approaches. Phys Chem Chem Phys 2018;20:8456.10.1039/C7CP08209F29557427 Search in Google Scholar
[164] Freysoldt C, Neugebauer J. Firstprinciples calculations for charged defects at surfaces, interfaces, and twodimensional materials in the presence of electric fields. Phys Rev B 2018;97:205425.10.1103/PhysRevB.97.205425 Search in Google Scholar
[165] Tahini HA, Tan X, Smith SC. Fermi level determination for charged systems via recursive density of states integration. J Phys Chem Lett 2018;9:4014–9.10.1021/acs.jpclett.8b0163129968476 Search in Google Scholar
[166] Smart TJ, Wu F, Govoni M, Ping Y. Fundamental principles for calculating charged defect ionization energies in ultrathin twodimensional materials. Phys Rev Mater 2018;2:124002.10.1103/PhysRevMaterials.2.124002 Search in Google Scholar
[167] Pinto H, Jones R, Palmer DW, et al. Firstprinciples studies of the effect of (001) surface terminations on the electronic properties of the negatively charged nitrogenvacancy defect in diamond. Phys Rev B 2012;86:045313.10.1103/PhysRevB.86.045313 Search in Google Scholar
[168] Löfgren R, Pawar R, Öberg S, Larsson JA. Charged dopants in neutral supercells through substitutional donor (acceptor): nitrogen donor charging of the nitrogenvacancy center in diamond. N J Phys 2018;20:023002.10.1088/13672630/aaa382 Search in Google Scholar
[169] Manson NB, Hedges M, Barson MSJ, et al. NV−–N+ pair centre in 1b diamond. N J Phys 2018;20:113037.10.1088/13672630/aaec58 Search in Google Scholar
[170] Collins AT. The Fermi level in diamond. J Phys Condensed Matter 2002;14:3743.10.1088/09538984/14/14/307 Search in Google Scholar
[171] Londero E, Bourgeois E, Nesladek M, Gali A. Identification of nickelvacancy defects by combining experimental and ab initio simulated photocurrent spectra. Phys Rev B 2018;97:241202.10.1103/PhysRevB.97.241202 Search in Google Scholar
[172] Bockstedte M, Schütz F, Garratt T, Ivády V, Gali A. Ab initio description of highly correlated states in defects for realizing quantum bits. npj Quant Mater 2018;3:31.10.1038/s4153501801036 Search in Google Scholar
[173] Ranjbar A, Babamoradi M, Heidari Saani M, Vesaghi MA, Esfarjani K, Kawazoe Y. Manyelectron states of nitrogenvacancy centers in diamond and spin density calculations. Phys Rev B 2011;84:165212.10.1103/PhysRevB.84.165212 Search in Google Scholar
[174] Onida G, Reining L, Rubio A. Electronic excitations: densityfunctional versus manybody Green’sfunction approaches. Rev Mod Phys 2002;74:601.10.1103/RevModPhys.74.601 Search in Google Scholar
[175] Ma Y, Rohlfing M, Gali A. Excited states of the negatively charged nitrogenvacancy color center in diamond. Phys Rev B 2010;81:041204.10.1103/PhysRevB.81.041204 Search in Google Scholar
[176] Choi S, Jain M, Louie SG. Mechanism for optical initialization of spin in NV− center in diamond. Phys Rev B 2012;86:041202.10.1103/PhysRevB.86.041202 Search in Google Scholar
[177] Thiering G, Gali A. Theory of the optical spinpolarization loop of the nitrogenvacancy center in diamond. Phys Rev B 2018;98:085207.10.1103/PhysRevB.98.085207 Search in Google Scholar
[178] Rogers LJ, Armstrong S, Sellars MJ, Manson NB. Infrared emission of the NV centre in diamond: zeeman and uniaxial stress studies. N J Phys 2008;10:103024.10.1088/13672630/10/10/103024 Search in Google Scholar
[179] Gali A, Janzén E, Deák P, Kresse G, Kaxiras E. Theory of spinconserving excitation of the N−V− center in diamond. Phys Rev Lett 2009;103:186404.10.1103/PhysRevLett.103.18640419905820 Search in Google Scholar
[180] Fu KMC, Santori C, Barclay PE, Rogers LJ, Manson NB, Beausoleil RG. Observation of the dynamic JahnTeller effect in the excited states of nitrogenvacancy centers in diamond. Phys Rev Lett 2009;103:256404.2036627010.1103/PhysRevLett.103.256404 Search in Google Scholar
[181] Ulbricht R, Dong S, Chang IY, et al. JahnTellerinduced femtosecond electronic depolarization dynamics of the nitrogenvacancy defect in diamond. Nat Commun 2016;7:13510.2784893810.1038/ncomms13510 Search in Google Scholar
[182] Zhang J, Wang CZ, Zhu ZZ, Dobrovitski VV. Vibrational modes and lattice distortion of a nitrogenvacancy center in diamond from firstprinciples calculations. Phys Rev B 2011;84:035211.10.1103/PhysRevB.84.035211 Search in Google Scholar
[183] Abtew TA, Sun YY, Shih BC, Dev P, Zhang SB, Zhang P. Dynamic JahnTeller effect in the NV− center in diamond. Phys Rev Lett 2011;107:146403.10.1103/PhysRevLett.107.14640322107219 Search in Google Scholar
[184] Bersurker I. The JahnTeller effect. Cambridge University Press, 2006. Search in Google Scholar
[185] Thiering G, Gali A. The (e_{g}⊗e_{u})⊗E_{g} product Jahn–Teller effect in the neutral groupIV vacancy quantum bits in diamond. npj Comput Mater 2019;5:1. Search in Google Scholar
[186] Rogers LJ, Doherty MW, Barson MSJ, Onoda S, Ohshima T, Manson NB. Singlet levels of the NV−centre in diamond. N J Phys 2015;17:013048.10.1088/13672630/17/1/013048 Search in Google Scholar
[187] Gali A, Simon T, Lowther JE. Singlet levels of the NV−centre in diamond. N J Phys 2011;13:025016. Search in Google Scholar
[188] Alkauskas A, Buckley BB, Awschalom DD, de Walle CGV. Firstprinciples theory of the luminescence lineshape for the triplet transition in diamond NV centres. N J Phys 2014;16:073026.10.1088/13672630/16/7/073026 Search in Google Scholar
[189] Thiering G, Gali A. Characterization of oxygen defects in diamond by means of density functional theory calculations. Phys Rev B 2016;94:125202.10.1103/PhysRevB.94.125202 Search in Google Scholar
[190] Runge E, Gross EKU. Densityfunctional theory for timedependent systems. Phys Rev Lett 1984;52:997.10.1103/PhysRevLett.52.997 Search in Google Scholar
[191] Casida ME. Timedependent density functional response theory for molecules. Recent Advances in Density Functional Theory. Singapore: World Scientific, 1995:155–92. Search in Google Scholar
[192] Gali A. Timedependent density functional study on the excitation spectrum of point defects in semiconductors. Phys Stat Solidi B 2011;248:1337.10.1002/pssb.201046254 Search in Google Scholar
[193] Perdew JP, Ernzerhof M, Burke K. Rationale for mixing exact exchange with density functional approximations. J Chem Phys 1996;105:9982–85.10.1063/1.472933 Search in Google Scholar
[194] Vlasov II, Shiryaev AA, Rendler T, et al. Molecularsized fluorescent nanodiamonds. Nat Nanotechnol 2014;9:54–8.10.1038/nnano.2013.25524317283 Search in Google Scholar
[195] Ulbricht R, Dong S, Gali A, Meng S, Loh ZH. Vibrational relaxation dynamics of the nitrogenvacancy center in diamond. Phys Rev B 2018;97:220302.10.1103/PhysRevB.97.220302 Search in Google Scholar
[196] Kehayias P, Doherty MW, English D, et al. Infrared absorption band and vibronic structure of the nitrogenvacancy center in diamond. Phys Rev B 2013;88:165202.10.1103/PhysRevB.88.165202 Search in Google Scholar
[197] Acosta VM, Bauch E, Ledbetter MP, Waxman A, Bouchard LS, Budker D. Temperature dependence of the nitrogenvacancy magnetic resonance in diamond. Phys Rev Lett 2010;104:070801.2036686810.1103/PhysRevLett.104.070801 Search in Google Scholar
[198] Bassett LC, Heremans FJ, Christle DJ, Yale CG, Burkard G, Buckley BB, Awschalom DD. Ultrafast optical control of orbital and spin dynamics in a solidstate defect. Science 2014;345:1333–7.10.1126/science.1255541 Search in Google Scholar
[199] Rayson MJ, Briddon PR. First principles method for the calculation of zerofield splitting tensors in periodic systems. Phys Rev B 2008;77:035119.10.1103/PhysRevB.77.035119 Search in Google Scholar
[200] Ivády V, Simon T, Maze JR, Abrikosov IA, Gali A. Pressure and temperature dependence of the zerofield splitting in the ground state of NV centers in diamond: a firstprinciples study. Phys Rev B 2014;90:235205.10.1103/PhysRevB.90.235205 Search in Google Scholar
[201] Bodrog Z, Gali A. The spin–spin zerofield splitting tensor in the projectoraugmentedwave method. J Phys Condensed Matter 2014;26:015305.10.1088/09538984/26/1/015305 Search in Google Scholar
[202] Seo H, Ma H, Govoni M, Galli G. Designing defectbased qubit candidates in widegap binary semiconductors for solidstate quantum technologies. Phys Rev Mater 2017;1:075002.10.1103/PhysRevMaterials.1.075002 Search in Google Scholar
[203] Giannozzi P, Baroni S, Bonini N, et al. QUANTUM ESPRESSO: a modular and opensource software project for quantum simulations of materials. J Phys Condensed Matter 2009;21:395502.10.1088/09538984/21/39/395502 Search in Google Scholar
[204] Biktagirov T, Schmidt WG, Gerstmann U. Calculation of spinspin zerofield splitting within periodic boundary conditions: towards allelectron accuracy. Phys Rev B 2018;97:115135.10.1103/PhysRevB.97.115135 Search in Google Scholar
[205] Ghosh K, Ma H, Gavini V, Galli G. Allelectron density functional calculations for electron and nuclear spin interactions in molecules and solids. Phys Rev Mater 2019;3:043801.10.1103/PhysRevMaterials.3.043801 Search in Google Scholar
[206] Davidsson J, Ivády V, Armiento R, Son NT, Gali A, Abrikosov IA. First principles predictions of magnetooptical data for semiconductor point defect identification: the case of divacancy defects in 4H–SiC. N J Phys 2018;20:023035.10.1088/13672630/aaa752 Search in Google Scholar
[207] Fuchs GD, Dobrovitski VV, Hanson R, et al. Excitedstate spectroscopy using single spin manipulation in diamond. Phys Rev Lett 2008;101:117601.1885133210.1103/PhysRevLett.101.117601 Search in Google Scholar
[208] Ham FS. Dynamical JahnTeller effect in paramagnetic resonance spectra: orbital reduction factors and partial quenching of spinorbit interaction. Phys Rev 1965;138:A1727.10.1103/PhysRev.138.A1727 Search in Google Scholar
[209] Thiering G, Gali A. Ab initio magnetooptical spectrum of groupiv vacancy color centers in diamond. Phys Rev X 2018;8:021063. Search in Google Scholar
[210] Stone NJ. Table of nuclear magnetic dipole and electric quadrupole moments. Atomic Data Nucl Data Tables 2005;90:75.10.1016/j.adt.2005.04.001 Search in Google Scholar
[211] Szász K, Hornos T, Marsman M, Gali A. Hyperfine coupling of point defects in semiconductors by hybrid density functional calculations: the role of core spin polarization. Phys Rev B 2013;88:075202.10.1103/PhysRevB.88.075202 Search in Google Scholar
[212] Yazyev OV, Tavernelli I, Helm L, Röthlisberger U. Core spinpolarization correction in pseudopotentialbased electronic structure calculations. Phys Rev B 2005;71:115110.10.1103/PhysRevB.71.115110 Search in Google Scholar
[213] Blügel S, Akai H, Zeller R, Dederichs PH. Hyperfine fields of 3D and 4D impurities in nickel. Phys Rev B 1987;35:3271–83.10.1103/PhysRevB.35.3271 Search in Google Scholar
[214] Larico R, Justo JF, Machado WVM, Assali LVC. Electronic properties and hyperfine fields of nickelrelated complexes in diamond. Phys Rev B 2009;79:115202.10.1103/PhysRevB.79.115202 Search in Google Scholar
[215] Weisskopf V, Wigner E. Berechnung der natürlichen Linienbreite auf Grund der Diracschen Lichttheorie. Z Phys 1930;63:54–73.10.1007/BF01336768 Search in Google Scholar
[216] Alegre TPM, Santori C, MedeirosRibeiro G, Beausoleil RG. Polarizationselective excitation of nitrogen vacancy centers in diamond. Phys Rev B 2007;76:165205.10.1103/PhysRevB.76.165205 Search in Google Scholar
[217] Batalov A, Zierl C, Gaebel T, et al. Temporal coherence of photons emitted by single nitrogenvacancy defect centers in diamond using optical Rabioscillations. Phys Rev Lett 2008;100:077401.1835259410.1103/PhysRevLett.100.077401 Search in Google Scholar
[218] Ulbricht R, Loh ZH. Excitedstate lifetime of the NV− infrared transition in diamond. Phys Rev B 2018;98:094309.10.1103/PhysRevB.98.094309 Search in Google Scholar
[219] Norambuena A, Reyes SA, MejíaLopéz J, Gali A, Maze JR. Microscopic modeling of the effect of phonons on the optical properties of solidstate emitters. Phys Rev B 2016;94:134305.10.1103/PhysRevB.94.134305 Search in Google Scholar
[220] Robledo L, Bernien H, van der Sar T, Hanson R. Spin dynamics in the optical cycle of single nitrogenvacancy centres in diamond. N J Phys 2011;13:025013.10.1088/13672630/13/2/025013 Search in Google Scholar
[221] Kalb N, Humphreys PC, Slim JJ, Hanson R. Dephasing mechanisms of diamondbased nuclearspin memories for quantum networks. Phys Rev A 2018;97:062330.10.1103/PhysRevA.97.062330 Search in Google Scholar
[222] Doherty MW, Manson NB, Delaney P, Jelezko F, Wrachtrup J, Hollenberg LCL. The nitrogenvacancy colour centre in diamond. Phys Rep 2013;528:1–45.10.1016/j.physrep.2013.02.001 Search in Google Scholar
[223] Vörös M, Rocca D, Galli G, Zimanyi GT, Gali A. Increasing impact ionization rates in Si nanoparticles through surface engineering: a density functional study. Phys Rev B 2013;87:155402.10.1103/PhysRevB.87.155402 Search in Google Scholar
[224] Alkauskas A, Yan Q, Van de Walle CG. Firstprinciples theory of nonradiative carrier capture via multiphonon emission. Phys Rev B 2014;90:075202.10.1103/PhysRevB.90.075202 Search in Google Scholar
[225] Mizuochi N, Makino T, Kato H, et al. Electrically driven singlephoton source at room temperature in diamond. Nat Photonics 2012;6:299–303.10.1038/nphoton.2012.75 Search in Google Scholar
[226] Schreckenbach G, Ziegler T. Calculation of the GTensor of electron paramagnetic resonance spectroscopy using gaugeincluding atomic orbitals and density functional theory. J Phys Chem A 1997;101:3388–99.10.1021/jp963060t Search in Google Scholar
[227] Pickard CJ, Mauri F. Allelectron magnetic response with pseudopotentials: NMR chemical shifts. Phys Rev B 2001;63:245101.10.1103/PhysRevB.63.245101 Search in Google Scholar
[228] Pickard CJ, Mauri F. FirstPrinciples theory of the EPR g tensor in solids: defects in quartz. Phys Rev Lett 2002;88:086403.10.1103/PhysRevLett.88.086403 Search in Google Scholar
[229] Kadantsev ES, Ziegler T. Implementation of a DFTbased method for the calculation of the zeeman gTensor in periodic systems with the use of numerical and SlaterType atomic orbitals. J Phys Chem A 2009;113:1327–34.1917364010.1021/jp805466c Search in Google Scholar
[230] Von Bardeleben HJ, Cantin JL, Rauls E, Gerstmann U. Identification and magnetooptical properties of the NV center in 4H−SiC. Phys Rev B 2015;92:064104.10.1103/PhysRevB.92.064104 Search in Google Scholar
[231] Von Bardeleben HJ, Cantin JL, Csóré A, Gali A, Rauls E, Gerstmann U. NV centers in 3C, 4H, and 6H silicon carbide: a variable platform for solidstate qubits and nanosensors. Phys Rev B 2016;94:121202.10.1103/PhysRevB.94.121202 Search in Google Scholar
[232] Gali A. Theory of the neutral nitrogenvacancy center in diamond and its application to the realization of a qubit. Phys Rev B 2009;79:235210.10.1103/PhysRevB.79.235210 Search in Google Scholar
[233] Doherty MW, Dolde F, Fedder H, et al. Theory of the groundstate spin of the NV− center in diamond. Phys Rev B 2012;85:205203.10.1103/PhysRevB.85.205203 Search in Google Scholar
[234] van Oort E, Glasbeek M. Electricfieldinduced modulation of spin echoes of NV centers in diamond. Chem Phys Lett 1990;168:529–32.10.1016/00092614(90)85665Y Search in Google Scholar
[235] Falk AL, Klimov PV, Buckley BB, et al. Electrically and mechanically tunable electron spins in silicon carbide color centers. Phys Rev Lett 2014;112:187601.10.1103/PhysRevLett.112.18760124856721 Search in Google Scholar
[236] Udvarhelyi P, Nagy R, Kaiser F, Lee SY, Wrachtrup J, Gali A. Spectrally stable defect qubits with no inversion symmetry for robust spintophoton interface. Phys Rev Appl 2019;11:044022.10.1103/PhysRevApplied.11.044022 Search in Google Scholar
[237] Udvarhelyi P, Shkolnikov VO, Gali A, Burkard G, Pályi A. Spinstrain interaction in nitrogenvacancy centers in diamond. Phys Rev B 2018;98:075201.10.1103/PhysRevB.98.075201 Search in Google Scholar
[238] Mittiga T, Hsieh S, Zu C, et al. Imaging the local charge environment of nitrogenvacancy centers in diamond. Phys Rev Lett 2018;121:246402.3060873210.1103/PhysRevLett.121.246402 Search in Google Scholar
[239] KingSmith RD, Vanderbilt D. Theory of polarization of crystalline solids. Phys Rev B 1993;47:1651–4.10.1103/PhysRevB.47.1651 Search in Google Scholar
[240] Vanderbilt D, KingSmith RD. Electric polarization as a bulk quantity and its relation to surface charge. Phys Rev B 1993;48:4442–55.10.1103/PhysRevB.48.4442 Search in Google Scholar
[241] Resta R. Macroscopic polarization in crystalline dielectrics: the geometric phase approach. Rev Mod Phys 1994;66:899.10.1103/RevModPhys.66.899 Search in Google Scholar
[242] Gajdoš M, Hummer K, Kresse G, Furthmüller J, Bechstedt F. Linear optical properties in the projectoraugmented wave methodology. Phys Rev B 2006;73:045112.10.1103/PhysRevB.73.045112 Search in Google Scholar
[243] Tamarat P, Gaebel T, Rabeau JR, et al. Stark shift control of single optical centers in diamond. Phys Rev Lett 2006;97:083002.10.1103/PhysRevLett.97.08300217026299 Search in Google Scholar
[244] Doherty MW, Acosta VM, Jarmola A, et al. Temperature shifts of the resonances of the NV− center in diamond. Phys Rev B 2014;90:041201.10.1103/PhysRevB.90.041201 Search in Google Scholar
[245] Kobayashi M, Nisida Y. High pressure effects on photoluminescence spectra of color centers in diamond. Jpn J Appl Phys 1993;32:279–81.10.7567/JJAPS.32S1.279 Search in Google Scholar
[246] Deng B, Zhang RQ, Shi XQ. New insight into the spinconserving excitation of the negatively charged nitrogenvacancy center in diamond. Sci Rep 2014;4:5144.24888367 Search in Google Scholar
[247] Barson MSJ, Peddibhotla P, Ovartchaiyapong P, et al. Nanomechanical sensing using spins in diamond. Nano Lett 2017;17:1496–503.10.1021/acs.nanolett.6b0454428146361 Search in Google Scholar
[248] Udvarhelyi P, Gali A. Ab Initio spinstrain coupling parameters of divacancy qubits in silicon carbide. Phys Rev Appl 2018;10:054010.10.1103/PhysRevApplied.10.054010 Search in Google Scholar
[249] Rogers LJ, McMurtrie RL, Sellars MJ, Manson NB. Alloptical thermometry and thermal properties of the optically detected spin resonances of the NV(−) center in nanodiamond. N J Phys 2009;11:063007. Search in Google Scholar
[250] Plakhotnik T, Doherty MW, Cole JH, Chapman R, Manson NB. Alloptical thermometry and thermal properties of the optically detected spin resonances of the NV− center in nanodiamond. Nano Lett 2014;14:4989–96.2510203510.1021/nl501841d Search in Google Scholar
[251] Plakhotnik T, Doherty MW, Manson NB. Electronphonon processes of the nitrogenvacancy center in diamond. Phys Rev B 2015;92. Search in Google Scholar
[252] Gugler J, Astner T, Angerer A, Schmiedmayer J, Majer J, Mohn P. Ab initio calculation of the spin lattice relaxation time T1 for nitrogenvacancy centers in diamond. Phys Rev B 2018;98:214442.10.1103/PhysRevB.98.214442 Search in Google Scholar
[253] Jarmola A, Acosta VM, Jensen K, Chemerisov S, Budker D. Temperature and magneticfielddependent longitudinal spin relaxation in nitrogenvacancy ensembles in diamond. Phys Rev Lett 2012;108:197601.10.1103/PhysRevLett.108.19760123003089 Search in Google Scholar
[254] Takahashi S, Hanson R, van Tol J, Sherwin MS, Awschalom DD. Quenching spin decoherence in diamond through spin bath polarization. Phys Rev Lett 2008;101:047601.1876436510.1103/PhysRevLett.101.047601 Search in Google Scholar
[255] Norambuena A, Muñoz E, Dinani HT, et al. Spinlattice relaxation of individual solidstate spins. Phys Rev B 2018;97:094304.10.1103/PhysRevB.97.094304 Search in Google Scholar
[256] Astner T, Gugler J, Angerer A, et al. Solidstate electron spin lifetime limited by phononic vacuum modes. Nat Mater 2018;17:313–17.10.1038/s415630170008y29434307 Search in Google Scholar
[257] Waller I. Über die Magnetisierung von paramagnetischen Kristallen in Wechselfeldern. Z Phys 1932;79:370–88.10.1007/BF01349398 Search in Google Scholar
[258] Marzari N, Mostofi AA, Yates JR, Souza I, Vanderbilt D. Maximally localized Wannier functions: theory and applications. Rev Mod Phys 2012;84:1419.10.1103/RevModPhys.84.1419 Search in Google Scholar
[259] Marini A. Ab Initiofinitetemperature excitons. Phys Rev Lett 2008;101:106405.10.1103/PhysRevLett.101.10640518851235 Search in Google Scholar
[260] Cannuccia E, Marini A. Effect of the quantum zeropoint atomic motion on the optical and electronic properties of diamond and transpolyacetylene. Phys Rev Lett 2011;107:255501.2224308910.1103/PhysRevLett.107.255501 Search in Google Scholar
[261] Gali A, Demján T, Vörös M, Thiering G, Cannuccia E, Marini A. Electronvibration coupling induced renormalization in the photoemission spectrum of diamondoids. Nat Commun 2016;7:11327.2710334010.1038/ncomms11327 Search in Google Scholar