# Experimental Realization of Universal Geometric Quantum Gates with Solid-State Spins

Experimental realization of a universal set of quantum logic gates is the central requirement for implementation of a quantum computer. An all-geometric approach to quantum computation 1 (); 2 () offered a paradigm for implementation where all the quantum gates are achieved based on the Berry phases 3 () and their non-abelian extensions, the holonomies 5 (), from geometric transformation of quantum states in the Hilbert space 5a (). Apart from its fundamental interest and rich mathematical structure, the geometric approach has some built-in noise-resilient features 1 (); 2 (); 6 (); 7 (). On the experimental side, geometric phases and holonomies have been observed using nuclear magnetic resonance with thermal ensembles of liquid molecules 8 (); 10 (), however, such systems are known to be non-scalable for quantum computing 11 (). There are proposals to implement geometric quantum computation in scalable experimental platforms such as trapped ions 12 (), superconducting qubits 13 (), or quantum dots 14 (), and a recent experiment has realized geometric single-bit gates with the superconducting system 15 (). Here, we report the first experimental realization of a universal set of geometric quantum gates with solid-state spins of the diamond defects. The diamond defects provide a scalable experimental platform 26 (); 16b (); 16 () with the potential for room-temperature quantum computing 16b (); 16 (); 17 (); 18b (), which has attracted strong interest in recent years 18 (). Based on advance of coherent control in this system 26 (); 16b (); 16 (); 17 (); 18b (); 18 (), our experiment shows that all-geometric and potentially robust quantum computation can be realized with solid-state spin qubits.

Under adiabatic cyclic evolution, a non-degenerate eigenstate of a quantum system acquires a phase factor, which has a dynamical component proportional to the time integral of the eigenenergy and a geometric component determined by the global property of the evolution path. This geometric phase, first discovered by Berry 3 (), has found connection with many important physics phenomena 19 (). If the system has degenerate eigenstates, the Berry’s phase is replaced by an geometric unitary operator acting on the degenerate subspace, termed as holonomy from the differential geometry. The holonomies are in general non-commutable with each other. In the proposal of geometric quantum computation 1 (); 2 (), such holonomies are exploited to realize a universal set of quantum gates, compositions of which then can fulfill arbitrary quantum computation tasks. As holonomies are determined by global geometric properties, geometric computation is more robust to certain control errors 1 (); 2 (); 6 (); 7 (). Implementation of geometric quantum computation has been proposed in several qubit systems 12 (); 13 (); 14 (), however, it remains experimentally challenging to realize a universal set of gates all by holonomies, because of the requirements of slow adiabatic evolution and a complicated level structure.

In the recent proposal of non-adiabatic geometric quantum computation 6 (); 20 (), universal quantum gates are constructed fully by geometric means without requirement of the adiabatic condition, thereby combining speed with universality. Under a cyclic evolution of the system Hamiltonian (with ), let () denote instantaneous orthonormal bases (moving frames) which coincide with the basisvectors of the computational space at with . The evolution operator on the basis states has two contributions: a dynamic part and a fully geometric part 6 (). If the parallel-transport condition is satisfied for any at any time , the dynamic contribution becomes identically zero, and is given by

(1) |

where denotes the time-ordered integration and represents the connection matrix 6 (). The form of is identical to the Wilczek-Zee holonomy in the adiabatic case 5 (); 6 ().

Our experiment realizes a universal set of quantum gates all by use of the nonadiabatic holonomies 6 (). Single-bit gates, together with entangling controlled-NOT (CNOT) operation, are universal for quantum computation. Our realization is based on control of electron and nuclear spins in a diamond nitrogen-vacancy (NV) center that form effectively a quantum register 18 (). To realize the single-bit geometric gates, we manipulate the electron spin states of a NV center (Fig. 1a) in a synthetic diamond at room temperature (see Methods for description of the experimental setup). The NV center has a spin-triplet ground state. We take the Zeeman components and as the qubit basis states and use as an ancillary level for geometric manipulation of the qubit. The spin state is initialized through optical pumping to the level and read out by distinguishing different fluorescence levels of the states under illumination of a short green laser pulse 18 () (see Methods for calibration of fluorescence levels of different states). We apply a magnetic field of G along the NV axis using a permanent magnet. Under this field, the nearby nuclear spins are polarized by optical pumping 22 (), enhancing the coherence time of the electron spin.

The transitions from the qubit states to the ancillary level are coupled by microwave pulses controlled through an arbitrary waveform generator (AWG), with Rabi frequencies , , respectively (Fig. 1b). We vary the amplitude but fix the ratio to be constant. The Hamiltonian for the coupling between these three levels takes the form

(2) |

where is the Planck constant divided by and denotes the Hermitian conjugate. Define the bright state as and the dark state as . When makes a cyclic evolution with , the bright state evolves as , where , while the dark state remains unchanged. After a cyclic evolution with , the bright (dark) state picks up a geometric phase of (), respectively, as illustrated in Fig. 1c. We take the moving frame as , , which makes a cyclic evolution with . For this evolution, one can easily check that the condition is always satisfied, so there is no dynamic contribution to the evolution operator 6 (). Using the expression (1), we find the holonomy

(3) |

under the computational basis .

We evolve the Rabi frequencies along three different loops, with the parameters chosen respectively as , , . The three geometric gates resulting from these cyclic evolutions are denoted by the NOT gate , the rotation gate , and the Hadamard gate , respectively. The combination of the gates and gives the well-known -gate , which, together with the Hadamard gate , make a universal set of single-bit gates. To characterize these geometric gates, we use quantum process tomography by preparing and measuring the qubit in different bases 23b (), with the time sequence shown in Fig. 1d. The matrix elements for each process are shown in Fig. 2a-2c, which are compared with the corresponding elements of the ideal gates. From the process tomography (see Methods), we find the process fidelity , , respectively for the , , and gates. The major contribution to the infidelity actually comes from the state preparation and detection error in quantum process tomography. To measure the intrinsic gate error, we concatenate a series of gates and examine the fidelity decay as the number of gates increases 18b (). As an example, we show in Fig. 2d the fidelity decay by concatenating the NOT gates. From the data, we find the intrinsic error per gate is about . This can be compared with the error rate for the dynamic NOT gate using optimized pulses by the same method of measurement 18b (). The achieved high fidelity indicates that the geometric manipulation is indeed resilient to control errors.

To realize the geometric quantum CNOT gate, we use one nearby nuclear spin as the control qubit (with the basisvectors ) and the NV center electron spin as the target qubit 22a (). Both the electron spin and the nuclear spin are polarized through optical pumping under the G magnetic field, which is confirmed by the optically detected magnetic resonance (ODMR) spectroscopy shown in Fig. 3b. The spins are interacting with each other through hyperfine and dipole couplings, and the resultant level configuration is shown in Fig. 3a. By applying state-selective microwave (MW) and radio-frequency (RF) pulses, we can couple different levels. In particular, with the MW0 and MW1 pulses with Rabi frequencies , , we have the following coupling Hamiltonian

(4) |

where we have fixed the ratio . Under a cyclic evolution of with , we find the holonomy using the formula (1), where denotes the unit matrix. This achieves exactly the quantum CNOT gate.

To characterize the geometric CNOT gate, we apply the gate to the qubit basis states as well as their superpositions, and measure the fidelity of the final states compared with the ideal outputs through the quantum state tomography 23b (). The superposition of the nuclear spin states required for state preparation and measurement is generated through RF pulses, which take longer times compared with the microwave pulses due to the much smaller magnetic moment of the nuclear spin. The electron spin decoherence is significant during the slow RF pulses. To correct that, we apply a Hahn spin echo in the middle of the whole operation with the time sequence shown in Fig. 3c. The measured state fidelities are listed in Fig. 4a under typical input states. As a hallmark of the entangling operation, the geometric CNOT gate generates entanglement from the initial product state. As an example, for the input state (unnormalized), the matrix elements of the output density operator are shown in Fig. 4b, with a measured entanglement fidelity of and concurrence of , which unambiguously confirms entanglement 11 ().

Our experimental realization of a universal set of holonomic gates with individual spins paves the way for all-geometric quantum computation in a solid-state system. The electron and nuclear spins of different NV centers can be wired up quantum mechanically to form a scalable network of qubits through, e.g., the direct dipole interaction 16b (); 17 (), the spin-chain assisted coupling by the nitrogen dopants 16 (); 23 (), or the photon-mediated coupling 26 (); 24 (); 25 (). The technique employed here for geometric realization of universal gates may also find applications in other scalable experimental systems, such as trapped ions or superconducting qubits. The geometric phase is closely related to the topological phase 27b (); 27 (), and the demonstration of gates all by holonomies is an important step towards realization of topological computation 27 (), the most robust way of quantum computing.

## References

- (1) Zanardi, P. & Rasetti, M. Holonomic quantum computation. Phys. Lett. A 264, 94–99 (1999).
- (2) Pachos, J., Zanardi, P. & Rasetti, M. Non-Abelian Berry connections for quantum computation. Phys. Rev. A 61,010305(R) (2000).
- (3) Berry, M. V. Quantal phase-factors accompanying adiabatic changes. Proc. R. Soc. Lond. A 392, 45–57 (1984).
- (4) Wilczek, F. & Zee, A. Appearance of gauge structure in simple dynamical systems. Phys. Rev. Lett. 52, 2111–2114 (1984).
- (5) Lloyd, S. Computation from geometry. Science 292, 1669 (2001).
- (6) Sjoqvist, E. et al. Non-adiabatic holonomic quantum computation. N. J. Phys. 14, 103035 (2012).
- (7) Johansson, M. et al. Robustness of non-adiabatic holonomic gates. Phys. Rev. A 86, 062322 (2012).
- (8) Jones, J. A., Vedral, V., Ekert, A. & Castagnoli, G. Geometric quantum computation using nuclear magnetic resonance. Nature 403, 869-871 (1999).
- (9) Feng, G., Xu, G. & Long G. Experimental Realization of Nonadiabatic Holonomic Quantum Computation. Phys. Rev. Lett. 110, 190501 (2013).
- (10) Nielsen, M. A. and Chuang, I. L. Quantum computation and quantum information. (Cambridge Univ. Press, 2010).
- (11) Duan, L. M., Cirac, J. I. & Zoller, P. Geometric manipulation of trapped ions for quantum computation. Science 292, 1695–1697 (2001).
- (12) Falci G. et al., Detection of geometric phases in superconducting nanocircuits. Nature 407, 355-358 (2000).
- (13) Solinas, P., Zanardi, P., Zanghi, N. & Rossi, F. Holonomic quantum gates: a semiconductor-based implementation. Phys. Rev. A 67, 062315 (2003).
- (14) Abdumalikov, A. A. et al. Experimental realization of non-Abelian non-adiabatic geometric gates. Nature 496, 482 (2013).
- (15) Pfaff, W. et al. Unconditional quantum teleportation between distant solid-state quantum bits. Science 10.1126, 1253512 (2014).
- (16) Neumann, P. et al. Scalable quantum register based on coupled electron spins in a room temperature solid. Nature Physics 6, 249-253 (2010).
- (17) Yao, N. Y. et al. Scalable architecture for a room temperature solid-state quantum information processor. Nature Commun. 3, 800 (2012).
- (18) Dolde, F. et al. Room-temperature entanglement between single defect spins in diamond. Nature Physics 9, 139-143 (2013).
- (19) Dolde, F. et al. High-fidelity spin entanglement using optimal control. Nature Comm. 5, 3371 (2014).
- (20) Doherty, M. W. et al. The nitrogen-vacancy colour centre in diamond. Physics Reports 528, 1-45 (2013).
- (21) Shapere, A. & Wilczek F. Geometric phases in physics (World Scientific, Singapore, 1989).
- (22) Zhu, S. L. & Wang, Z. D. Implementation of universal quantum gates based on nonadiabatic geometric phases. Phys. Rev. Lett. 89, 097902 (2002).
- (23) Jacques, V. et al. Dynamic polarization of single nuclear spins by optical pumping of nitrogen-vacancy color centers in diamond at room temperature. Phys. Rev. Lett. 102, 057403 (2009).
- (24) Jelezko, F. et al. Observation of coherent oscillation of a single nuclear spin and realization of a two-qubit conditional quantum gate. Phys. Rev. Lett. 93, 130501 (2004).
- (25) White, A. G. et al. Measuring two-qubit gates. J. Opt. Soc. Am. B, 24, 172-183 (2007).
- (26) Yao, N. Y. et al. Robust quantum state transfer in random unpolarized spin chains. Phys. Rev. Lett. 106, 040505 (2011).
- (27) Jiang, L., Taylor, J. M., Sorensen, A. S. & Lukin, M. D. Distributed Quantum Computation Based-on Small Quantum Registers. Phys. Rev. A 76, 062323 (2007).
- (28) Duan, L. M. & Monroe, C. Quantum networks with trapped ions. Rev. Mod. Phys. 82, 1209 (2010).
- (29) Loredo, J. C., Broome, M. A., Smith, D. H., & White, A. G. Observation of Entanglement-Dependent Two-Particle Holonomic Phase. Phys. Rev. Lett. 112, 143603 (2014).
- (30) Pachos, J. K. Introduction to Topological Quantum Computation (Cambridge Univ. Press, 2012).

Acknowledgements We thank M. Lukin’s group for helpful discussions. This work was supported by the National Basic Research Program of China 2011CBA00300 (2011CBA00302) and the quantum information project from the Ministry of Education of China. LMD acknowledges in addition support from the IARPA MUSIQC program, the AFOSR and the ARO MURI program.

Author Contributions L.M.D. conceived the experiment and supervised the project. C.Z., W.B.W., L.H., W.G.Z., C.Y.D., F.W. carried out the experiment. L.M.D. and C.Z. wrote the manuscript.

Author Information Reprints and permissions information is available at www.nature.com/reprints. The authors declare no competing financial interests. Correspondence and requests for materials should be addressed to L.M.D. (lmduan@umich.edu).

## Appendix A Methods

### a.1 Experimental setup

We use a home-built confocal microscopy, with an oil-immersed objective lens (), to address and detect single NV center in a type IIa single-crystal synthetic diamond sample (Element Six). A nm diode laser, controlled by an acoustic optical modulator (AOM), is used for spin state initialization and detection. We collect fluorescence photons (wavelength ranging from nm) into a single-mode fiber and detect them by the single-photon counting modular (SPCM), with a counting rate kHz and a signal-to-noise ratio . The diamond sample is mounted on a -axis closed-loop Piezo for sub-micrometer resolution scanning. An impedance-matched gold coplanar waveguide (CPW) with m gap, deposited on a cover-glass, is used for delivery of radio-frequency (RF) and microwave (MW) signals to the NV center.

In our experiment, we find a single NV center with a proximal of MHz hyperfine strength (Fig. 1). To polarize the nearby nuclear spins ( and the host ), we apply a magnetic field of G along the NV axis using a permanent magnet. Under this field, the electron spin levels and become almost degenerate in the optically excited state (called the esLAC, the electron spin level anti-crossing 22 ()), which facilitates electron-spin nuclear-spin flip-flop process during optical pumping. The spin flip-flop process leads to polarization of the nitrogen nuclear spin on the NV site and the nearby nuclear spins after green laser illumination 22 (). The Zeeman energy from the G magnetic field shifts the energy difference between electron spin states and () from the zero-field splitting MHz to MHz ( MHz) and the nuclear spin hyperfine splitting from MHz to MHz ( MHz) for () levels. Due to the large splitting of levels, we apply two independent MW sources (Rohde-Schwarz), locked by a MHz reference Rubidium clock, to address each transition. To adjust the frequency and phase of the MW pulses, we mix each MW output with an AWG (Tektronix, MHz sample rate). RF signals for nuclear spin manipulation are generated directly by another analog channel of the AWG. All the MW and RF signals are amplified by independent amplifiers, combined through a home-made circuit, and delivered to the CPW. The digital markers of the AWG are used to control the pulse sequence (including laser and SPCM) with a timing resolution of ns.

For each experimental cycle, we start the sequence with laser illumination to polarize NV electron spin and nearby nuclear spins and end it with a laser pulse for spin state detection. We collect signal photons for ns right after the detection laser rises, and another ns for reference later. With a photon collection rate of kHz, we have an average of photon counts per cycle. For measurement of each data, we repeat the experimental cycle at least times, resulting in a total photon counts of . The error bars of our data account for the statistical error associated with the photon counting. To calculate the error bar of each data, we use Monte Carlo simulation by assuming a Poissonian distribution for the photon counts. For each simulation trial, we calculate the value of each data. Then, by sampling over all the trails according to the Poissonian distribution, we get statistics of the data (including its mean value and standard deviation, the error bar).

### a.2 Calibration of fluorescence levels for different states

Due to the esLAC that induces spin flip-flop during the detection and the imperfect initial polarization of the electron and nuclear spins, each spin component ( ) may fluorescent at different levels. Note that the spins are dominantly in the state after the optical pumping. To calibrate the fluorescence level of each state, we therefore associate the detected fluorescence level right after the optical pumping with the state . With MW or RF -pulses (the -pulses are calibrated through Rabi oscillations), we can make a complete transfer between and any other spin component. For instance, with a -pulse between and right after the optical pumping, we associate the detected fluorescence level with the state. In this way, the characteristic fluorescence level of each component can be calibrated. With the calibrated fluorescence level for each spin component, we then read out the system state after the geometric gates through quantum state tomography 23b ().

### a.3 Quantum Process tomography

A quantum process can be described by a completely positive map acting on an arbitrary initial state , transferring it to . In quantum process tomography (QPT), we choose a fixed set of basis operators so that the map is identified with a process matrix . We experimentally measure this process matrix through the maximum likelihood technique 23b (). For single-bit QPT, we set the basis operators as , , , and choose four different initial states , , , and . We reconstruct the corresponding final density operators through the standard quantum state tomography and use them to calculate the process matrix . This process matrix is compared with the ideal one by calculating the process fidelity . The process fidelity also determines the average gate fidelity by the formula 23b (), where is defined as the fidelity averaged over all possible input states with equal weight and is the dimension of the state space (with for a single qubit).