## Abstract

The nitrogen-vacancy (NV) color center in diamond has rapidly emerged as an important solid-state system for quantum information processing. Whereas individual spin registers have been used to implement small-scale diamond quantum computing, the realization of a large-scale device requires the development of an on-chip quantum bus for transporting information between distant qubits. Here, we propose a method for coherent quantum transport of an electron and its spin state between distant NV centers. Transport is achieved by the implementation of spatial stimulated adiabatic Raman passage through the optical control of the NV center charge states and the confined conduction states of a diamond nanostructure. Our models show that, for two NV centers in a diamond nanowire, high-fidelity transport can be achieved over distances of order hundreds of nanometers in timescales of order hundreds of nanoseconds. Spatial adiabatic passage is therefore a promising option for realizing an on-chip spin quantum bus.

## 1 Introduction

Diamond-based quantum information processing (QIP) is possible due to the remarkable properties of the nitrogen-vacancy (NV) defect center. The spin states of individual NV centers can be optically addressed for initialization and read-out [1] and possess the longest room-temperature spin-coherence time of any solid-state defect [2]. Cryogenic temperatures further enhance these spin properties by up to four orders of magnitude [3] and enable spin-photon entanglement through coherent optical transitions [4]. This includes the adiabatic control of the NV spin using optical pulses [5]. These properties have been used to operate individual spin registers for small-scale quantum computing, including demonstrations of error correction [6], simple algorithms [7], and simulations [8]. Spin registers have been realized through clusters of NV centers as well as single NV centers coordinated with clusters of paramagnetic defects such as ^{13}C atoms and substitutional nitrogen (N_{S}) centers [9], [10], [11], [12], [13], [14], [15].

Individual spin registers are ultimately limited by the number of paramagnetic defects addressable by a single NV center (5–10 strongly coupled qubits [16]). The development of large-scale QIP therefore requires the use of quantum buses to entangle a network of interconnected clusters. At cryogenic temperatures, optical quantum buses have been used to entangle spin registers separated by macroscopic distances [4], [17], [18]. Unfortunately, this option is challenging for scalable QIP as optical buses are difficult to compactly incorporate on chip while maintaining low losses. One scalable suggestion is intracluster spin transport along a chain of paramagnetic defects [19], [20]. Despite the promise of room-temperature operation, engineering the spin chain is difficult with existing fabrication techniques. Recently, Doherty et al. [21] have proposed entanglement using semiclassical electron transport between two NV–^{14}N_{S} pairs embedded in a diamond nanowire. This scheme is severely limited as semiclassical transport in a nanowire is prone to erroneous capture from

Consequently, here we propose an on-chip spin quantum bus for diamond QIP using spatial stimulated adiabatic Raman passage (STIRAP). Spatial STIRAP involves the transport of massive particles between spatially distinct locations using conventional STIRAP-type coupling. It was originally explored in the context of electronic transfer in double quantum dots [24], [25] and for general methods of spatial adiabatic passage, which promise robust quantum information transport [26], [27]. The conventional STIRAP technique is a method of population transfer between the two lowest levels of a Λ-type atomic system [28]. By applying optically coupling fields in the so-called counterintuitive direction, where the unpopulated transition is coupled before the populated transition, the system can be maintained in the optical dark state, which suppresses population in the central (excited) state. This suppression of population in the central state, which will ideally have zero population in the adiabatic limit, leads to the suppression of the effects of spontaneous emission. In the context of NV centers, STIRAP has been proposed [29] and demonstrated [30] for single qubit control; however, we are not aware of any proposals for the spatial transport of electrons between NV centers of the form considered here.

In our proposal for spatial STIRAP, the two disconnected states of the Λ scheme correspond to an electron occupying one of two distant NV centers embedded in a diamond nanostructure. As depicted in Figure 1A, we envision mediating adiabatic passage through states of the conduction band minimum, which are discretized by the confining potential of the nanostructure. While all nanostructures lead to discretization, we will consider nanowires as an archetypal structure. We have identified two different designs for realizing a nanowire confining potential, via the potential of a diamond surface or an electrostatic potential applied by electrodes in a bulk structure. The schematics for the transport scheme in these designs are presented in Figure 1. Although STIRAP is performed similarly in both cases, the susceptibility to different decoherence mechanisms varies. We also predict that quantum transport will maintain coherence of the electron’s spin. Spin transport is known theoretically to be excellent in diamond due to an indirect band-gap, low spin-orbit coupling, low electron-phonon (e-p) scattering, and low background impurity spins [31]. Here, under STIRAP we explicitly exploit the low background impurities and low spin-orbit coupling of the conduction band to achieve spin coherent transport.

This proposal for quantum transport is fundamentally different from existing schemes for spatial STIRAP. As coupling is mediated optically rather than through the manipulation of tunneling amplitudes [27], the spatially separated states are truly disconnected and there is no occupation of the intermediary space during transport. In this sense, it is colloquially akin to the teleportation of a massive particle. Furthermore, to our knowledge, this is the first spatial STIRAP proposal that includes coherent spin transport.

In Section 1, we present a proposal for implementing STIRAP in a diamond nanowire including a spin initialization and a measurement protocol. These spin protocols play the essential role of identifying the transported electrons (thereby distinguishing them from background sources) and will be the means for which spin entanglement is mediated. We identify a feasibility condition for transport, Δ*E _{c}*/

*ħ*≫Ω≫Γ, where Δ

*E*

_{c}defines the energy gap between the first and second conduction states, Ω is the Rabi frequency of the optical transition between defect and conduction states, and Γ is the optical decoherence rate. As a conservative estimate, these constraints can be considered satisfied when Δ

*E*/

_{c}*ħ*≳10Ω≳10Γ. In Section 2, we calculate Δ

*E*

_{c}by modeling the effect of the nanowire confining potential on the diamond conduction band. Effective mass theory (EMT) is used to derive a level scheme and determine the energies of discretized conduction states. In Section 3, we present the first

*ab initio*calculations for the photoionization mechanism of the NV center. This is used to evaluate values of Ω that can be practically achieved. We then identify optical decoherence mechanisms that contribute to Γ. Finally, in Section 4, we evaluate Ω, Δ

*E*

_{c}, and Γ at cryogenic temperatures for both nanowire designs introduced previously. The dimensions of the wires are optimized to maximize the Rabi frequency relative to decoherence and to satisfy the STIRAP feasibility constraints.

## 2 Quantum transport in diamond nanowires

STIRAP is a method for evolving a three-state system from an initial state |1⟩ to a final state |3⟩ [28], [32]. This is achieved through coupling to a mutual intermediary state, |2⟩, although the direct transition |1⟩→|3⟩ may be forbidden. To do so, a pump laser pulse with Rabi frequency Ω* _{P}*(

*t*) is tuned to the |1⟩→|2⟩ transition, whereas a Stokes laser pulse with Rabi frequency Ω

*(*

_{S}*t*) is tuned to the |2⟩→|3⟩ transition. Applying the Stokes laser pulse before the pump laser pulse with some overlapping time τ results in an adiabatic evolution of the state from |1⟩→|3⟩ without any population of |2⟩. The requirement for adiabaticity is given by Ω

*≫1, where*

_{τ}We propose using STIRAP for quantum transport of an electron between two NV centers. In Figure 1, NVs are embedded within the nanowire confining potential and labeled as *A* and *B*. This allows us to define a three-level system given by

where |NV^{−}⟩ and |NV^{0}⟩ correspond to the ^{3}A_{2} and ^{2}E electronic ground states of the negative and neutral charge states of the NV centers and |*ψ*⟩ is a nanowire conduction band minimum state. A schematic of the STIRAP level scheme is presented in Figure 2. The adiabatic passage between states |1⟩ and |3⟩ therefore constitutes the electron transport from NV *A* to NV *B* with no occupation of the intermediary space. The degeneracy of states |1⟩ and |3⟩ may be lifted by applying a potential difference Φ between each end of the length *L* wire. For a defect separation distance *s*, this induces an energy splitting of Δ*E _{d}*=

*e*Φ

*s*/

*L*.

Let us now consider the transportation of the electron spin during STIRAP. This requires a protocol for spin initialization of both |NV^{−}⟩* _{A}* and |NV

^{0}⟩

*, which we now present in turn. First, suppose we wish to transport the spin-state |↑⟩. We restrict our attention to photoionization by absorption of a single photon from an optical pulse and do not consider the two-photon ionization process of the NV*

_{B}^{−}[33]. As the ionizing pulse cannot differentiate between the two electrons occupying the

*e*-orbitals of |NV

^{−}⟩

*, both must be initialized into |↑⟩. This can be achieved through the optical polarization of |NV*

_{A}^{−}⟩

*into the*

_{A}*m*=0 triplet state followed by the application of a microwave control pulse. The resulting spin-state can be expressed in the two

_{s}*e*-orbital basis as |+1⟩=|↑⟩|↑⟩. Hence, through this protocol, we can guarantee the ionization of a spin-up electron during STIRAP.

Now, consider the initialization of the |NV^{0}⟩* _{B}* spin-state. After STIRAP, electron transportation can only be validated if the read-out of the |NV

^{−}⟩

*spin-state agrees with the spin initialized at |NV*

_{B}^{−}⟩

*. Otherwise, it is impossible to discern if the charging of |NV*

_{A}^{0}⟩

*into |NV*

_{B}^{−}⟩

*was instead due to the erroneous electron capture. It is therefore essential for the read-out that |NV*

_{B}^{0}⟩

*is initialized into either the |↑⟩ or |↓⟩ spin eigenstates. To see how this could be achieved, consider that the |NV*

_{B}^{−}⟩

*spin-state after transport (assuming |NV*

_{B}^{0}⟩

*was initialized as |↑⟩) is given by |↑⟩|↑⟩=|+1⟩. Similarly, we obtain the*

_{B}*m*=0 triplet-state when |NV

_{s}^{0}⟩

*is initialized as |↓⟩. Transport can then be confirmed by performing read-out of the spin state at NV*

_{B}*B*.

The spin initialization of NV *B* can be achieved as follows. Consider NV *B* initially prepared into the negative charge state with spin projection *m _{s}*=+1. The application of a >2.6 eV photoionizing pulse will then eject an electron with |↑⟩ spin, leaving NV

*B*in the neutral charge state with |↑⟩ spin. Equivalently, the |↓⟩ spin state can be prepared for NV

*B*prepared in the

*m*=−1 state. Although not proven, it is assumed that photoionization is a spin-conserving process. Even if this is not the case, one way to ameliorate the risk of nonconserving spin transitions is to polarize the N nuclear spin when NV

_{s}*B*is in the negative charge state. After the ionizing pulse, this nuclear spin polarization can be swapped back onto the electron spin of the neutral charge state via microwave control. Such a scheme is now feasible given that recent progress has identified the NV

^{0}spin Hamiltonian [34].

Transportation also requires that NV^{0} spin relaxation is long with respect to the STIRAP operation time. As of yet, the spin relaxation of the NV^{0} ground state has not been directly observed due to the apparent absence of an electron paramagnetic resonance (EPR) signal [35]. Although this could potentially imply a high rate of spin relaxation, the new work by Barson et al. indicated that this is likely not the case. Rather, the apparent absence of an EPR signal is due to strain broadening in ensembles [34]. Drawing analogy between the similar electronic groundstates of the NV^{0} and Si-V^{−} centers, we expect spin relaxation to be not overly fast. This is because within the Si-V^{−} center, e-p scattering is known to cause spin dephasing but not relaxation [36]. Consequently, at cryogenic temperatures, NV *B* should maintain its spin projection throughout transport.

During transport, spin decoherence mechanisms can be neglected. This remains true as long as the spin-quantization axis of the conduction states match those of the NV center and there is no rapid spin dephasing of these conduction states. This is investigated in Supplementary Section B, where we present symmetry arguments within the closure approximation. We prove that if the NV is coaxial to the wire, then the spin-quantization axis is conserved throughout STIRAP. Conveniently, NV centers can be preferentially aligned with a (111) wire axis during fabrication with chemical vapor deposition (CVD) [37]. For other nanowire directions, the spin-quantization axis of the conduction state can be aligned with the NV center via the application of a magnetic field. The magnitude of this field must be sufficiently large relative to the component of the spin-orbit interaction orthogonal to (111).

There are two further requirements for the successful implementation of STIRAP. First, coherent transport demands that there be no off-resonant coupling to states other than the conduction band minimum. Hence, the energy splitting between the first and second conduction states Δ*E*_{c} must be significantly greater than the Rabi frequency, Δ*E _{c}*/

*ħ*≫Ω. Second, the optical decoherence rate Γ sets a bound on Ω, as it induces an unknown phase on the electron wave function as the transition energy fluctuates. Provided that the Rabi frequency is large (Ω≫Γ), this dephasing will be sufficiently small to achieve high-fidelity transport. The following two sections address each of these feasibility requirements in turn. We first investigate the effects of wire confinement on the conduction band before evaluating the Rabi frequency and identifying sources of optical decoherence.

## 3 Diamond nanowire conduction band

The confining potential of the nanowire leads to the discretization of the conduction band suitable for STIRAP. These confinement effects are well described through EMT [38], as detailed in Supplementary Section A. As shown in Figure 1, the nanowire was modeled as a square prism of monocrystalline diamond with dimensions of *w*×*w*×*L* (*w*<*L*). The two NV centers are separated by a distance *s* and located on the longitudinal. Consider that bulk diamond possesses six equivalent conduction band minima located at the *k*-points *X* points in each orthogonal axis of the fcc Brillouin zone. The dispersion relationship at each

In the effective mass framework, the *n*th conduction band state for the *i*th valley of the nanowire is given by

where *F _{n,i}* is an envelope function determined by the confining potential and

*u*is the bulk-diamond Bloch function. For simplicity, we approximate the confining potential as an infinite square well. Although more realistic confining potentials such as a finite well would result in perturbations to the electronic structure, the key features are well described within our approximate model. This leads to the following envelope function

_{i}where *n*=(*n _{x}*,

*n*,

_{y}*n*) and

_{z}*V*is the volume of the diamond unit cell.

_{c}The true eigenstates of the nanowire conduction band are linear combinations of the valley wavefunctions (1), which respect to its *D*_{4h} symmetry. As presented in Supplementary Section A, symmetry analysis reveals four degenerate conduction band minima states formed from the valleys with _{1g}, B_{1g}, and E_{u}). Similarly, the valley states with _{1g} and *A*_{2u} symmetry. For *L*≫*w*, the energy of the twofold degenerate states exceeds that of the fourfold degenerate states, so only the latter are relevant for STIRAP. For the simple case of a wire aligned along a (100) axis, the eigenenergies of the fourfold conduction band minima states are given by

where *m*_{⊥} and *m*_{||} are the transverse and longitudinal effective masses of an electron in bulk diamond. For *L*≫*w*, this corresponds to

In Section 4, we use Equation (4) to assess the feasibility condition Δ*E _{c}*/

*ħ*≫Ω for the two different nanowire designs.

EMT is limited as it cannot account for coupling between valley states via the wire confining potential. Known as the valley-orbit interaction, this coupling leads to a fine splitting of the conduction band minima states [40]. Although the resulting eigenstates may be determined by symmetry, the magnitude of the splitting is notoriously difficult to quantify but assumed to be small [41], [42], [43], [44]. A proper treatment would require detailed *ab initio* calculations, which are left as an opportunity for future work. The final electronic consideration is splitting induced by the applied electric field. As explored in Supplementary Section A, this Stark effect does not influence the states of the conduction band minimum. Although it does produce a small splitting of the A_{1g} and *A*_{2u} states, for *L*≫*w*, this poses no impediment to the performance of STIRAP.

## 4 Optical processes in diamond nanowires

Assessing the feasibility of STIRAP requires comparing attainable photoionization Rabi frequencies to optical decoherence mechanisms. In this section we first present calculations of the photoionization Rabi frequency including Franck-Condon effects. We then identify and characterize three optical decoherence mechanisms; erroneous capture by

### 4.1 Photoionization Rabi frequency

The conduction dipole moment was first obtained for the photoionization of bulk NV centers using *ab initio* calculations. Density functional theory was performed using the VASP plane-wave code with a 512-atom diamond supercell [45], [46], [47], [48] and Heyd-Scuseria-Ernzerhof functionals [49]. The energy of the ^{3}A_{2} ground state of the NV^{−} center and the ^{2}E ground state of the NV^{0} center (plus an ionized electron) were calculated as a function of generalized atomic coordinates. As detailed in Supplementary Section C, the calculations produced an ionization energy of *ω*=2.6 eV (close to experimental values [33]) and a Huang-Rhys factor of *S*=1.39. The direct calculation of the optical matrix elements yields a transition dipole moment of *d*_{bulk}=0.085*e* Å, normalized to the volume of the 512-atom supercell.

The bulk photoionization dipole moment can be used to determine the corresponding moment in a nanowire. Noting that this moment is identical for both the pump and Stokes pulse, we calculate in Supplementary Section C that

where *A*. For defects centered at opposite ends of the wire, the normalization function in Equation (5) can be simplified as

where *V*_{sc}=2.837 nm^{3} is the volume of the 512-atom supercell. Assuming a laser power *P* with radial spot-size *r*, the effective photoionization Rabi frequency within a diamond nanowire is

where *n _{D}* is the index of refraction in diamond. The achievable spot size and power are related to each other depending on the technical aspects of the optics employed. A conservative estimate says that a diffraction-limited spot size

*r*=200 nm can be attained with a lasing power of

*P*=100 mW. These parameters will be used in Section 4 to evaluate the photoionization Rabi frequency.

### 4.2 Decoherence mechanisms

One source of optical decoherence is the erroneous capture of conduction band electrons by _{S} defects are formed during CVD growth due to N-based gases present in the plasma [50]. These defects act as electron donors, readily ionizing to ^{2} is the *k*_{B} is the Boltzmann constant, *T* is the temperature, *m** is the isotropic mean of the effective mass components, *ρ*_{e} is the probability density of the conduction band states [21].

We may conservatively estimate that *ρ_{e}L*

^{2}

*w*<1; thus, the capture rate is independent of volume. Assuming that

*T*=4 K and

_{cap}≲2.5 MHz. This places a constraint on the photoionization Rabi frequency for STIRAP in any nanostructure. In this analysis, we have assumed passivation of all surface traps by N

_{S}donors. If this is not the case, surface defects must be considered as another source of electron capture. Other defects, such as vacancy clusters introduced during fabrication, may also act as electron traps. However, they will not be dominant in number within a single wire in comparison to N

_{S}defects.

A second source of decoherence can be attributed to the e-p scattering of the conduction band states. We consider only first-order processes involving absorption and emission of acoustic phonons with higher-energy conduction states. Higher-order interactions can be considered negligible at 4 K. The e-p scattering rate, Γ_{ep}, may be determined through Fermi’s golden rule with the Hamiltonian

where Ξ* _{d}*=8.7 eV is the diamond deformation potential [51] and

_{ep}are specific to the nanowire design and are presented in the following section.

A fundamental source of optical decoherence is that due to spontaneous emission. The spontaneous emission rate for transitions |2⟩→|1⟩ and |2⟩→|3⟩ may be estimated using conventional expressions as [52]

note that this decoherence rate is dependent on wire geometry through the normalization function *F* as per Equation (5). We now evaluate the Rabi frequency and total decoherence rate within both nanowire designs – surface confinement and electrostatic confinement – and optimize wire dimensions to satisfy the STIRAP feasibility conditions.

## 5 Optimization of nanowire design

### 5.1 Surface confinement

As displayed in Figure 1A, electronic confinement may be realized through the fabrication of freestanding nanowires [53], [54] or nanopillars [55] using reactive ion etching. To avoid injection of spurious electrons, it is necessary that the electrodes affixed to the wire ends are insulating and separated from the diamond surface with a high resistivity contact. NV centers may be introduced into wires through either N ion-implantation [56], [57] or N δ-doping [58]. A pair of NV centers suitable for STIRAP can then be identified by confocal microscopy. Ideally, the NVs will be coaligned along a (111) axis and situated at depths greater than ~50 nm (this is necessary for obtaining optimal bulk-like properties [59]).

We now derive an expression for the decoherence rate due to e-p scattering. In Supplementary Section D, we present calculations of Γ_{ep} by approximating *m*=(*m _{x}*,

*m*,

_{y}*m*), which are given by

_{z}where *c _{l}* is the longitudinal speed of sound in diamond. The scattering rate can then be calculated as

where *M _{n,m}* is the overlap integral between the electron and phonon wavefunctions, ρ

_{C}is the density of diamond, and

*n*is the Bose-Einstein distribution. The density of states,

_{B}*ρ*, is assumed to be a Lorentzian, given by

where *γ*≈*ω _{m}*/

*Q*for

*Q*is the quality factor. For diamond microcantilevers,

*Q*>10

^{6}has been observed [61].

We now optimize the dimensions of a surface-confined nanowire with respect to the STIRAP feasibility condition Δ*E _{c}*/

*ħ*≫Ω≫Γ. In Figure 3A, we present the photoionization Rabi frequency derived in Equation (7) as a function of wire dimension. In Figure 3B, we compare the Rabi frequencies to the total decoherence rate Γ=Γ

_{cap}+ Γ

_{SE}+ Γ

_{ep}for varying nanowire dimensions. We assume that

*Q*=10

^{6},

*T*=4 K and that each NV center is positioned 100 nm from its respective wire end. For simplicity, we have only considered wires aligned along a (100) axis, although similar results likely hold for other wire directions.

Figure 3B indicates that, in general, decoherence effects are minimized for smaller wire dimensions. This can be attributed to the normalization of the Rabi frequency, which scales inversely proportional with wire volume as evident in Figure 3A [cf. Equation (7)]. The sporadic variations in Ω/Γ can be attributed to the resonance between electronic and photonic levels. Consulting Equation (11), certain wire dimensions produce an increased density of resonant states, which amplifies the rate of e-p scattering. Off resonance, the dominant decoherence mechanism is electron capture. A range of wire dimensions with *L*<0.6 μm and *w*<0.2 μm satisfy the STIRAP feasibility condition with at least Ω>10Γ, reaching a maximum of Ω≈20Γ for short wires 0.3 μm long. Note that for all wire dimensions considered in Figure 3B, Δ*E _{c}*/

*ħ*>10

^{3}Ω; therefore, off-resonant coupling does not pose an impediment to STIRAP.

### 5.2 Electrostatic confinement

An alternative means for realizing a nanowire confining potential is through electrodes in a bulk structure. Consider the diamond substrate presented in Figure 1B, where two coaligned NV centers can be identified through confocal microscopy. To generate the confining potential, a square nanoelectrode with width *w* is affixed to the surface directly above the NV centers. On the opposing surface of the diamond substrate, an electrode plate is fixed and grounded. Applying a positive potential to the square electrode will then generate a confining potential. This localizes the electron density near the surface into a wire-like geometry. The effective length of this wire may be tuned through the magnitude of the potential. Conveniently, the potential also generates an electric field gradient that can be used to energetically distinguish each NV center during STIRAP. These electrostatic properties are demonstrated through simulations using COMSOL Multiphysics software presented in Supplementary Section E.

In contrast to surface-confined wires, electron scattering within electrostatically confined wires is due to interactions with bulk phonons. The calculation of Γ_{ep} using Fermi’s golden rule and the Hamiltonian (8) is presented in Supplementary Section D. We find that

where *G _{n}* is a complicated expression for the overlap integral.

A further benefit of the electrostatic confinement design is that decoherence due to _{S} defects distant from the confining potential can be δ-doped into the substrate during CVD growth [58]. The sacrificial defects donate electrons to the surface traps until charge conservation occurs, reducing the density of _{cap} by 95%.

In Figure 3C, we present the dimensional dependence of Ω/Γ for an electrostatically confined nanowire aligned along the (100) axis. We have that *G _{n}*∝

*w*

^{−}^{4}in Equation (13); therefore, e-p scattering dominates the Rabi frequency for wire widths less than approximately 0.4 μm. For larger wire dimensions, both e-p scattering and spontaneous emission are negligible and the dominant source of decoherence is due to

*w*<0.9 μm and 0.75 μm<

*L*<1.25 μm, which yields Ω>10Γ. The ratio Ω/Γ is optimized at ≈25 for wire dimensions of

*L*=0.75 μm and

*w*=0.6 μm. Furthermore, for all wire dimensions, Δ

*E*/

_{c}*ħ*≫10

^{3}Ω; hence, crosstalk with higher-energy states in the conduction band is negligible.

## 6 Conclusion

In this paper, we propose spatial STIRAP to realize an on-chip spin quantum bus for scalable diamond QIP devices. Our scheme considers coherent quantum transport of an electron and its spin state between two NV centers embedded in a diamond nanowire. This is achieved by implementing STIRAP for the optical control of the NV center charge states and the confined nanowire conduction states. In contrast to existing spatial STIRAP protocols that manipulate tunneling amplitudes between distant sites, our proposal realizes quantum transport without any occupation of the intermediary space. Furthermore, to our knowledge, this is the first implementation of spatial STIRAP including spin-state transport.

A proposal for our scheme was first presented along with a protocol for initializing the NV spin-states. We identified a feasibility condition for implementing STIRAP, Δ*E _{c}*/

*ħ*≫Ω≫Γ, requiring the photoionization Rabi frequency to exceed optical decoherence rates and be limited by coupling with higher-energy conduction states. For the latter, EMT was employed to model the effects of the nanowire confining potential on the diamond conduction band. We then presented the first

*ab initio*calculations of the NV photoionization Rabi frequency. This allowed the STIRAP feasibility condition to be evaluated in two nanowire designs, which achieved electronic confinement through either a surface structure or electrodes. By optimizing the Rabi frequency relative to the decoherence rate for varying wire dimension, we conclude that STIRAP is feasible in both designs in timescales on the order of hundreds of nanoseconds.

The first experimental stages should perform photoionization spectroscopy to assess the spectral linewidth of discretized states in the nanowire conduction band. Ultimately, this will determine whether deterministic photoionization to the conduction band minimum is feasible. A second key experiment is to directly the measure the spin-relaxation time of the NV^{0} center to confirm that it is sufficiently long to enact the STIRAP protocol. This study has also introduced multiple avenues for future theoretical research. For example, a detailed *ab initio* model of the nanowire electronic structure would allow for the precise calculation of valley-orbit and spin-orbit interactions. This would require the development of new computational tools for simulating mesoscopic diamond structures.

**Funding source: **Australian Research Council

**Award Identifier / Grant number: **DE170100169

**Award Identifier / Grant number: **DP170103098

**Award Identifier / Grant number: **FT160100357

**Award Identifier / Grant number: **CE170100026

**Funding source: **National Science Foundation

**Award Identifier / Grant number: **NSF-1547830

**Award Identifier / Grant number: **NSF-1619896

**Award Identifier / Grant number: **NSF-1914945

**Funding source: **Research and Innovation Programme

**Award Identifier / Grant number: **820394

**Funding statement: **We acknowledge funding from the Australian Research Council (DE170100169, DP170103098, FT160100357, and CE170100026). C.A.M. acknowledges the support from the National Science Foundation through grants NSF-1547830, NSF-1619896, and NSF-1914945 and the Research Corporation for Science Advancement through a FRED Award. A.A. acknowledges the funding from the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement no. 820394 (project Asteriqs). This research was undertaken with the assistance of resources and services from the National Computational Infrastructure, which is supported by the Australian Government.

## References

[1] Doherty MW, Manson NB, Delaney P, et al. The nitrogen-vacancy colour centre in diamond. Phys Rep 2013;528:1–45.10.1016/j.physrep.2013.02.001Search in Google Scholar

[2] Balasubramanian G, Neumann P, Twitchen D, et al. Ultralong spin coherence time in isotopically engineered diamond. Nat Mater 2009;8:383–7.10.1038/nmat2420Search in Google Scholar PubMed

[3] Jarmola A, Acosta VM, Jensen K, et al. Temperature- and magnetic-field-dependent longitudinal spin relaxation in nitrogen-vacancy ensembles in diamond. Phys Rev Lett 2012;108:197601.10.1103/PhysRevLett.108.197601Search in Google Scholar PubMed

[4] Togan E, Chu Y, Trifonov AS, et al. Quantum entanglement between an optical photon and a solid-state spin qubit. Nature 2010;466:730–4.10.1038/nature09256Search in Google Scholar PubMed

[5] Yale CG, Buckley BB, Christle DJ, et al. All-optical control of a solid-state spin using coherent dark states. Proc Natl Acad Sci 2013;110:7595–600.10.1073/pnas.1305920110Search in Google Scholar PubMed PubMed Central

[6] Taminiau TH, Cramer J, van der Sar T, et al. Universal control and error correction in multi-qubit spin registers in diamond. Nat Nanotechnol 2014;9:171–6.10.1038/nnano.2014.2Search in Google Scholar PubMed

[7] Shi F, Rong X, Xu N, et al. Room-temperature implementation of the Deutsch-Jozsa algorithm with a single electronic spin in diamond. Phys Rev Lett 2010;105.10.1103/PhysRevLett.105.040504Search in Google Scholar PubMed

[8] Wang Y, Dolde F, Biamonte J, et al. Quantum simulation of helium hydride cation in a solid-state spin register. ACS Nano 2015;9:7769–74.10.1021/acsnano.5b01651Search in Google Scholar PubMed

[9] Gaebel T, Domhan M, Popa I, et al. Room-temperature coherent coupling of single spins in diamond. Nat Phys 2006;2:408–13.10.1038/nphys318Search in Google Scholar

[10] Hanson R, Mendoza FM, Epstein RJ, et al. Polarization and readout of coupled single spins in diamond. Phys Rev Lett 2006;97:087601.10.1103/PhysRevLett.97.087601Search in Google Scholar PubMed

[11] Neumann P, Mizuochi N, Rempp F, et al. Multipartite entanglement among single spins in diamond. Science 2008;320:1326–9.10.1126/science.1157233Search in Google Scholar PubMed

[12] Neumann P, Kolesov R, Naydenov B, et al. Quantum register based on coupled electron spins in a room-temperature solid. Nat Phys 2010;6:249–53.10.1038/nphys1536Search in Google Scholar

[13] Dolde F, Jakobi I, Naydenov B, et al. Room-temperature entanglement between single defect spins in diamond. Nat Phys 2013;9:139–43.10.1038/nphys2545Search in Google Scholar

[14] Dolde F, Bergholm V, Wang Y, et al. High-fidelity spin entanglement using optimal control. Nat Commun 2014;5:3371.10.1038/ncomms4371Search in Google Scholar PubMed

[15] Waldherr G, Wang Y, Zaiser S, et al. Quantum error correction in a solid-state hybrid spin register. Nature 2014;506:204–7.10.1038/nature12919Search in Google Scholar PubMed

[16] Doherty M. The state of diamond quantum computing. Aust Phys 2017;110:7595–600.Search in Google Scholar

[17] Bernien H, Hensen B, Pfaff W, et al. Heralded entanglement between solid-state qubits separated by three metres. Nature 2013;497:86–90.10.1038/nature12016Search in Google Scholar PubMed

[18] Hensen B, Bernien H, Dreaú AE, et al. Loophole-free Bell inequality violation using electron spins separated by 1.3 kilometres. Nature 2015;526:682–6.10.1038/nature15759Search in Google Scholar PubMed

[19] Yao NY, Jiang L, Gorshkov AV, et al. Scalable architecture for a room temperature solid-state quantum information processor. Nat Commun 2012;3:800.10.1038/ncomms1788Search in Google Scholar PubMed

[20] Nikolopoulos GM, Jex I. Quantum state transfer and network engineering. 2014:1–250.10.1007/978-3-642-39937-4Search in Google Scholar

[21] Doherty MW, Meriles CA, Alkauskas A, et al. Towards a room-temperature spin quantum bus in diamond via electron photoionization, transport, and capture. Phys Rev X 2016;6:041035.10.1103/PhysRevX.6.041035Search in Google Scholar

[22] Stacey A, Dontschuk N, Chou JP, et al. Evidence for primal sp2 defects at the diamond surface: candidates for electron trapping and noise sources. Adv Mater Interfaces 2019;6:201801449.10.1002/admi.201801449Search in Google Scholar

[23] Kaviani M, Deák P, Aradi B, et al. Proper surface termination for luminescent near-surface NV centers in diamond. Nano Lett 2014;14:4772–7.10.1021/nl501927ySearch in Google Scholar PubMed

[24] Brandes T, Vorrath T. Adiabatic transfer of electrons in coupled quantum dots. Phys Rev B Condensed Matter Mater Phys 2002;66:753411–12.10.1103/PhysRevB.66.075341Search in Google Scholar

[25] Vitanov NV, Rangelov AA, Shore BW, et al. Stimulated Raman adiabatic passage in physics, chemistry, and beyond. Rev Mod Phys 2017;89:015006.10.1103/RevModPhys.89.015006Search in Google Scholar

[26] Greentree AD, Cole JH, Hamilton AR, et al. Coherent electronic transfer in quantum dot systems using adiabatic passage. Phys Rev B Condensed Matter Mater Phys 2004;70:1–6.10.1103/PhysRevB.70.235317Search in Google Scholar

[27] Menchon-Enrich R, Benseny A, Ahufinger V, et al. Spatial adiabatic passage: a review of recent progress. Rep Prog Phys 2016;79:074401.10.1088/0034-4885/79/7/074401Search in Google Scholar PubMed

[28] Bergmann K, Theuer H, Shore BW. Coherent population transfer among quantum states of atoms and molecules. Rev Mod Phys 2002;70:1003–25.10.1103/RevModPhys.70.1003Search in Google Scholar

[29] Coto R, Jacques V, Hétet G, et al. Stimulated Raman adiabatic control of a nuclear spin in diamond. Phys Rev B 2017;96:1003–25.10.1103/PhysRevB.96.085420Search in Google Scholar

[30] Zhou BB, Baksic A, Ribeiro H, et al. Accelerated quantum control using superadiabatic dynamics in a solid-state lambda system. Nat Phys 2017;13:330–4.10.1038/nphys3967Search in Google Scholar

[31] Restrepo OD, Windl W. Full first-principles theory of spin relaxation in group-IV materials. Phys Rev Lett 2012;109:166604.10.1103/PhysRevLett.109.166604Search in Google Scholar PubMed

[32] Ivanov PA, Vitanov NV, Bergmann K. Effect of dephasing on stimulated Raman adiabatic passage. Phys Rev A Atomic Mol Opt Phys 2004;70:063409.10.1103/PhysRevA.70.063409Search in Google Scholar

[33] Aslam N, Waldherr G, Neumann P, et al. Photo-induced ionization dynamics of the nitrogen vacancy defect in diamond investigated by single-shot charge state detection. N J Phys 2013;15:013064.10.1088/1367-2630/15/1/013064Search in Google Scholar

[34] Barson MSJ, Krausz E, Manson NB, Doherty MW. The fine structure of the neutral nitrogen-vacancy center in diamond. Nanophotonics 2019;8:1985–91.10.1515/nanoph-2019-0142Search in Google Scholar

[35] Felton S, Edmonds AM, Newton ME, et al. Electron paramagnetic resonance studies of the neutral nitrogen vacancy in diamond. Phys Rev B Condensed Matter Mater Phys 2008;77:081201.10.1103/PhysRevB.77.081201Search in Google Scholar

[36] Jahnke KD, Sipahigil A, Binder JM, et al. Electron-phonon processes of the silicon-vacancy centre in diamond. N J Phys 2015;17:043011.10.1088/1367-2630/17/4/043011Search in Google Scholar

[37] Michl J, Teraji T, Zaiser S, et al. Perfect alignment and preferential orientation of nitrogen-vacancy centers during chemical vapor deposition diamond growth on (111) surfaces. Appl Phys Lett 2014;104:102407.10.1063/1.4868128Search in Google Scholar

[38] Stoneham AM. Theory of defects in solids. Acta Crystallogr Section A 1976;32:527.10.1107/S0567739476001162Search in Google Scholar

[39] Madelung O. Semiconductors: data handbook. Microelectron J 2004;35:685.10.1016/j.mejo.2004.04.001Search in Google Scholar

[40] Murphy-Armando F, Fagas G, Greer JC. Deformation potentials and electron-phonon coupling in silicon nanowires. Nano Lett 2010;10:869–73.10.1021/nl9034384Search in Google Scholar PubMed

[41] Resca L, Resta R. Large binding due to dispersive screening and Bloch function interference in many-valley semiconductors. Solid State Commun 1979;29:275–7.10.1016/0038-1098(79)91055-XSearch in Google Scholar

[42] Resta R. A note on the many-valley effective mass theory. J Phys C Solid State Phys 1977;10.10.1088/0022-3719/10/8/004Search in Google Scholar

[43] Luttinger JM, Kohn W. Motion of electrons and holes in perturbed periodic fields. Phys Rev 1955;97:869–83.10.1103/PhysRev.97.869Search in Google Scholar

[44] Altarelli M, Hsu WY, Sabatini RA. Donor binding energies in multivalley semiconductors. J Phys C Solid State Phys 1977;10.10.1088/0022-3719/10/21/004Search in Google Scholar

[45] Kresse G, Furthmüller J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput Mater Sci 1996;6:15–50.10.1016/0927-0256(96)00008-0Search in Google Scholar

[46] Kresse G, Hafner J. *Ab initio* molecular-dynamics simulation of the liquid-metalamorphous-semiconductor transition in germanium. Phys Rev B 1994;49:14251–69.10.1103/PhysRevB.49.14251Search in Google Scholar

[47] Kresse G, Furthmüller J. Efficient iterative schemes for *ab initio* total-energy calculations using a plane-wave basis set. Phys Rev B Condensed Matter Mater Phys 1996;54:11169–86.10.1103/PhysRevB.54.11169Search in Google Scholar

[48] Kresse G, Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys Rev B Condensed Matter Mater Phys 1999;59:1758–75.10.1103/PhysRevB.59.1758Search in Google Scholar

[49] Heyd J, Scuseria GE, Ernzerhof M. Hybrid functionals based on a screened Coulomb potential. J Chem Phys 2003;118:8207–15.10.1063/1.1564060Search in Google Scholar

[50] Samlenski R, Haug C, Brenn R, et al. Incorporation of nitrogen in chemical vapor deposition diamond. Appl Phys Lett 1995;67:2798.10.1063/1.114788Search in Google Scholar

[51] Nava F, Canali C, Jacoboni C, et al. Electron effective masses and lattice scattering in natural diamond. Solid State Commun 1980;33:475–7.10.1016/0038-1098(80)90447-0Search in Google Scholar

[52] Fox M. Quantum optics: an introduction. 1st ed. Oxford: Oxford University Press, 2006.Search in Google Scholar

[53] Babinec TM, Hausmann BJ, Khan M, et al. A diamond nanowire single-photon source. Nat Nanotechnol 2010;5:195–9.10.1038/nnano.2010.6Search in Google Scholar PubMed

[54] Burek MJ, de Leon NP, Shields BJ, et al. Free-standing mechanical and photonic nanostructures in single-crystal diamond. Nano Lett 2012;12:6084–9.10.1021/nl302541eSearch in Google Scholar PubMed

[55] Momenzadeh SA, Stöhr RJ, de Oliveira FF, et al. Nanoengineered diamond waveguide as a robust bright platform for nanomagnetometry using shallow nitrogen vacancy centers. Nano Lett 2015;15:165–9.10.1021/nl503326tSearch in Google Scholar PubMed

[56] Bayn I, Chen EH, Trusheim ME, et al. Generation of ensembles of individually resolvable nitrogen vacancies using nanometer-scale apertures in ultrahigh-aspect ratio planar implantation masks. Nano Lett 2015;15:1751–8.10.1021/nl504441mSearch in Google Scholar PubMed

[57] Pezzagna S, Wildanger D, Mazarov P, et al. Nanoscale engineering and optical addressing of single spins in diamond. Small 2010;6:2117–21.10.1002/smll.201000902Search in Google Scholar PubMed

[58] Ohashi K, Rosskopf T, Watanabe H, et al. Negatively charged NVC in a 5nm thin ^{12}C diamond film. Nano Lett 2013;13:4733–8.10.1021/nl402286vSearch in Google Scholar
PubMed

[59] Ohno K, Joseph Heremans F, Bassett LC, et al. Engineering shallow spins in diamond with nitrogen delta-doping. Appl Phys Lett 2012;101.10.1063/1.4748280Search in Google Scholar

[60] Nishiguchi N, Ando Y, Wybourne MN. Acoustic phonon modes of rectangular quantum wires. J Phys Condensed Matter 1997;9:5751–64.10.1088/0953-8984/9/27/007Search in Google Scholar

[61] Tao Y, Boss JM, Moores BA, et al. Single-crystal diamond nanomechanical resonators with quality factors exceeding one million. Nat Commun 2014;5:3638.10.1038/ncomms4638Search in Google Scholar PubMed

## Supplementary Material

The online version of this article offers supplementary material (https://doi.org/10.1515/nanoph-2019-0144).

**Received:**2019-05-16

**Revised:**2019-08-04

**Accepted:**2019-08-05

**Published Online:**2019-08-30

©2019 Marcus W. Doherty et al., published by De Gruyter, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 Public License.