Magnetic ordering of nitrogen-vacancy centers in diamond via resonator-mediated coupling

Nitrogen-vacancy centers in diamond, being a promising candidate for quantum information processing, may also be an ideal platform for simulating many-body physics. However, it is difficult to realize interactions between nitrogen-vacancy centers strong enough to form a macroscopically ordered phase under realistic temperatures. Here we propose a scheme to realize long-range ferromagnetic Ising interactions between distant nitrogen-vacancy centers by using a mechanical resonator as a medium. Since the critical temperature in the long-range Ising model is proportional to the number of spins, a ferromagnetic order can be formed at a temperature of tens of millikelvin for a sample with $\sim10^4$ nitrogen-vacancy centers. This method may provide a new platform for studying many-body physics using qubit systems.

Introduction. The negatively charged nitrogen-vacancy (NV) centers in high-purity diamond have been considered as a promising candidate for solid state quantum information processing due to their long coherence time [1][2][3] and high feasibility in initialization, control, and readout of their spin states [4]. Simulation of many-body physics using NV centers, in analogue to cold atom physics [5,6], has been proposed [7][8][9][10]. To realize phase transitions in NV center qubit systems under realistic temperatures, however, sufficiently strong interactions between NV centers located sufficiently close (<30 nm) is till highly challenging [7,8,11]. A new opportunity is to use resonators as mediators [7,12,13], which has potential of coupling NV centers at distance.
In this letter, we propose to realize long-range coupling between many separated NV centers via a mechanical resonator. A remarkable feature of the long-range interacting system is: the critical temperature for ferromagnetic phase transitions is proportional to the number of spins, so a ferromagnetic order could be formed at a temperature of tens of millikelvin (mK) for a sample with ∼ 10 4 NV centers, even though the mediated coupling between two NV centers is less than 200 kHz(∼ 10 µK).
Microscopic model. The essential idea can be understood by considering a prototype system shown in Fig. 1(a). A mechanical resonator with frequency ω r is attached with a magnet. The NV sample is placed right under the magnet with a distance d. The oscillation of the mechanical resonator generates a time dependent magnetic field that causes Zeeman shift of the spins in the sample. In contrast to the schemes of direct dipolar coupling between NV centers, our scheme does not require short distance between NV centers. We use a diamond sample where the orientation axes of the NV color centers can lie along one of the four possible crystallographic axes in diamonds, [111] ( Fig. 1(b)) which can be grown by chemical vapor deposition method [14,15]. The NV centers have a spin triplet ground state (S = 1) with a large zero field splitting ∆ = 2.87 GHz. We apply a static magnetic field B NV along the [111] direction of the sample so that two Zeeman states of the electron spin qubits with magnetic quantum number |1 and |0 are nearly degenerate, δ ≡ (∆ − γ e B NV ) ∼ 0 ( Fig. 1(c)).
We make use of the two Zeeman states of the electron spin |1 and |0 and define a pseudo-spin by σ z = |1 1| − |0 0| and the corresponding spin flip operators σ † = |1 0| and σ − = |0 1|. In each NV center there is a 14 N nuclear spin which interacts with the on-site NV center spin. The Hamiltonian of the NV center near the degenerate point 1 MHz is the 14 N nuclear spin quadrupole splitting [16], γ N is the gyromagnetic ratio of the nitrogen nuclear spin and A ≈ 2 MHz is the hyperfine coupling between the electronic spin and the 14 N nuclear spin [16].
The motion of the mechanical resonator is described by the Hamiltonian H r = ω r b † b, with ω r as the fundamental vibration mode of the resonator and b as the corresponding annihilation operator. For example a silicon nitride string resonator has dimensions (325 × 0.35 × 0.1)µm with ω r = 2π × 1.0 MHz and Q = 1.3 × 10 6 at room temperature [17].
The magnetic field felt by the electronic spin of the NV centers can be approximated by a magnetic dipole [18], , where d j is the distance between the equilibrium position of the resonator and j−th NV center, G m is the magnetic field gradient at the position of the NV centers and x is the amplitude of mechanical resonator oscillation (x ∼ 10 −12 m at temperature of mK). This magnetic field will induce a Zeeman shift to the NV spins with Hamiltonian H z = j η j σ x j (b + b † ) and η j = γ e |G m |a 0 , where γ e is the gyromagnetic ratio for electron and a 0 = √ /2mω r is the amplitude of zero point fluctuations for a resonator of mass m.
Due to the distance distribution of the NV centers to the magnet, the coupling of the NV centers to the mechanical resonator have a distribution. A magnetic tip with size of 100 nm produces a magnetic gradient G m ∼ 7.8 × 10 6 T/m at a distance 25 nm away from the tip [19]. In such case η/2π could reach 200 kHz. For a sample with typical distance of NV centers ∼ 20 nm, the neighbor NV centers have a direct dipolar interaction with strength ∼ 5 kHz (∼ 0.25 µK), much less than the critical temperature to be discussed later. For the sake of simplicity, in the following we assume the coupling between the mechanical resonator and NV centers to be uniform and neglect the direct dipolar coupling between NV centers. Without these assumptions, however, the results in this paper would only be quantitatively affected. The coupling constant η/2π ∼ 200 kHz considerably exceeds both the electronic spin coherence time (milliseconds) of the NV centers and the intrinsic damping rate of the mechanical resonator, γ = ω r /Q of high-Q mechanical resonator [17]. We neglect the interaction between the mechanical resonator and the 14 N nuclear spins since they are three orders of magnitude smaller than that between the NV center and the mechanical resonator. Summarizing the three terms described above, the total microscopic Hamiltonian including the spin qubits in the NV centers, nitrogen nuclear spins, the mechanical resonator and the interaction among them is, Effective long-range Ising model. To obtain the effective interactions between NV centers, we calculate the partition function of the hybrid system Z = Tr[e −βH ]. We shall neglect the Zeeman energy of the 14 N nuclear spin since its strength (∼ 0.3 MHz) is much smaller than the zero field splitting ∼ 5 MHz and the hyperfine interaction ∼ 3 MHz. The hyperfine interaction between the electronic spin and the 14 N nuclear spin term AI · σ can be separated into the diagonal term AI x σ x and the flip-flop term A(I y σ y + I z σ z ). Let us first neglect the flip-flop term and the qubit splitting (δ = 0). By partial trace of the phonon bath and the nuclear spins, the partition function can be factorised as where Z phonon is the partition function of an independent phonon bath, Z14 N is the partition function of the nuclear spins and the effective Hamiltonian for the NV electron spins is with h(β) = β −1 ln[(1 + 2 cosh(βA x ))/3] being an effective magnetic field resulting from the on-site interaction between the NV center and the 14 N nuclear spin. The physics for the effective long-range interactions between the NV centers is that an NV center spin flips by virtually emitting a phonon, and then the phonon is virtually absorbed by another NV center spin far away with a spin flip. Therefore two NV spins at distance can flip-flop via virtual exchange of phonons, which leads to long-range σ x i σ x j coupling. In the temperature range we are interested in(∼ 10 mK), the effective field, h(β) ∼ 0.3µK, is negligible. Therefore the distant NV centers interact through an effective long-range ferromagnetic Ising interaction with the Hamiltonian where the ferromagnetic coupling strength J = 2η 2 /ω r . If the flip-flop term of the hyperfine interaction is taken into account, it would contribute an effective transverse field to the electronic spins, which, however, is ∼ 0.3µK and hence negligible. Ferromagnetic ordering. The long-range Ising model can be solved by mapping the ensemble onto to a single large spin. With the notation S x = N j=1 σ x j /2, the Hamiltonian of the long-range Ising model can be written as The spin quantum number S takes values 0, 1, 2, · · · , N/2 for even N and 1/2, 3/2, · · · , N/2 for odd N. Moreover, the spin degeneracy in each large spin subspace is D In the large N limit [20] Z ≈ e −β/2 N where ϕ = (1/2 + X) ln(1/2 + X) + (1/2 − X) ln(1/2 − X) − 2βX 2 andβ = βN J. The integration in the partition function can be evaluated by the saddle point approximation method. Forβ < 1, the saddle point appears at X = 0, and forβ > 1, there are two symmetric saddle point located respectively within (−1/2, 0) and (0, 1/2). Thereforeβ = 1 is the critical point for the ferromagnetic phase transition in the long-range Ising model. So the critical temperature for the long-range ferromagnetic Ising model, k B T c = N J , is proportional to the number of spins and the coupling strength [21].
We present the scaled free energy of the long-range Ising model as a function of the spin magnetization M x ≡ S x /N for different temperatures in Fig. 2. One can see that at high temperature ( Fig. 2(a)) the free energy minimum is situated at M x = 0 and a paramagnetic state is stable, while at low temperature ( Fig. 2(b)) the free energy minimum is doubly degenerated with M x = 1 or M x = −1, and a ferromagnetic state is preferred. For the hybrid system, the phase transition temperature can be controlled through the detuning from the degeneracy point between the NV spin states |1 and |0 . If the NV qubits have a small splitting δ = ∆ − γB NV , the resonator mediated interaction between the NV centers becomes The critical temperature becomes Here the magnetic field is tuned so that the NV spin states |1 and |0 are near degenerate. Hence the splitting δ is much smaller than the critical temperature. In Fig. 3 we plot the phase diagrams of the NV centers as a function of splitting for N = 10000. In a diamond sample with NV concentration ∼ 50 ppb, coupling between a resonator and such many NV centers is possible. A ferromagnetic order can be formed at the temperature of 50 mK. In Fig. 4(a), we present the numerical and the mean field results of the order parameter as a function of reduced temperature, k B T/N J. One can see that the order parameter vanishes at and above the critical temperature k B T/N J = 1 and also the finite size numerical results approach to that of the mean field predictions reasonably well as the number of spins increases. The specific heat shown in Fig. 4(b) has a finite jump which signals a second order phase transitions occurs.
Conclusion. In summary, we propose a scheme to realize ferromagnetic ordering of distant nitrogen-vacancy centers by using a mechanical resonator to mediate long-range Ising-type interaction. The critical temperature for the ferromagnetic phase transition in the long-range Ising model is proportional to the number of spins, so the ferromagnetic order could be formed at the temperature of tens of millikelvin for a sample with ten thousand nitrogen-vacancy centers. In addition, it may also be possible to use a superconducting resonator [22,23] as a medium to realize the long-range ferromagnetic coupling between the NV centers. Since the interactions between the NV centers mediated by a superconducting resonator is usually small compared to that by a mechanical resonator, high density NV samples are required to observe the magnetic ordering in the NV centers.
We are grateful to P. Bertet for useful discussions.