# Universal non-adiabatic holonomic gates in quantum dots and single-molecule magnets

###### Abstract

Geometric manipulation of a quantum system offers a method for fast, universal, and robust quantum information processing. Here, we propose a scheme for universal all-geometric quantum computation using non-adiabatic quantum holonomies. We propose three different realizations of the scheme based on an unconventional use of quantum dot and single-molecule magnet devices, which offer promising scalability and robust efficiency.

###### pacs:

03.65.Vf, 03.67.Lx, 73.63.Kv, 75.50.Xx^{†}

^{†}: New J. Phys.

## 1 Introduction

Holonomic quantum computation (HQC), originally conceived by Zanardi and Rasetti [1], has become one of the key approaches to perform robust quantum computation. The idea of HQC is based on the Wilczek-Zee holonomy [2] that generalizes the Berry phase [3] to non-Abelian (non-commuting) geometric phases accompanying adiabatic evolution. It was shown in Ref. [1] that adiabatic quantum holonomies generically allow for universal quantum computation. Different physical settings for adiabatic HQC have been proposed [4, 5, 6].

Recently, universal HQC based on Anandan’s non-adiabatic non-Abelian geometric phase [7] has been proposed [8]. This scheme allows for high-speed implementations of quantum gates; a feature that can be particularly useful in solid-state devices, where qubits typically have short coherence time. Non-adiabatic HQC can be made robust to a wide class of errors, such as decay, dephasing, and fluctuations in control parameters [9], and can be realized in decoherence free subspaces and noiseless subsystems insensitive to various collective errors [10, 11]. First experimental demonstrations of the scheme for universal non-adiabatic HQC proposed in Ref. [8] have been carried out for a superconducting qubit [12] and in an NMR system [13].

Here, we develop a setting for conditional holonomic gates in a four-level configuration. We illustrate how non-adiabatic holonomic gates in our scheme can be implemented in solid state devices consisting of few-electron quantum dots [14] and single-molecule magnets (SSMs) [15], which are promising candidates for quantum computation [16, 17, 18, 19]. By utilizing the scalability of these systems, the geometric nature of the proposed holonomic gates is a feature that may allow for robust and fast manipulations of a large number of qubits.

The outline of the paper is as follows. In the next section, we demonstrate one- and two-qubit quantum gates based on non-adiabatic quantum holonomies. The scheme is based on a generic four-level model with off-diagonal couplings. Realizations of the model in various kinds of solid-state devices are identified and examined in Sec. 3. The paper ends with the conclusions.

## 2 Non-adiabatic quantum holonomy

The general structure of our system is described by a four-dimensional effective state space whose dynamics is governed by a Hamiltonian of the form

(1) |

in the ordered orthonormal basis . We assume that is a complex-valued and time-independent matrix that satisfies , where the latter ensures that is invertible. The Hamiltonian is turned on and off as described by the time-dependent scaling function . We demonstrate in the following how can be used for universal non-adiabatic HQC and how it can be realized in systems of four coupled quantum dots or in triangular molecular antiferromagnets.

The Hamiltonian in Eq. (1) naturally splits the state space into two subspaces spanned by and . Non-adiabatic holonomy transformations are realized if (i) the Hamiltonian vanishes on each evolving subspace, and (ii) the projection operators and move in a cyclic fashion. Condition (i) is satisfied since , , where the first equality follows from the fact that commutes with the time evolution operator . For condition (ii), we first note that since is assumed to be invertible, there is a unique singular value decomposition , where and are the positive and unitary parts, respectively. By using the singular value decomposition, the time evolution operator can be expressed as

(2) |

where . The evolution of the ’s is thus cyclic provided there is a such that . Under this condition the holonomy matrices for the two subspaces read

(3) |

where with . and are the paths of the two subspaces spanned by and , respectively. These paths reside in the Grassmannian manifold , being the space of all two-dimensional subspaces of the four-dimensional state space. One can see from Eq. (3) that the holonomy is non-trivial if and only if .

Next, we encode the basis states in into the computational two-qubit states , , , and . Running the Hamiltonian along certain time intervals for appropriate choices of provides a realization of holonomic conditional gates of the form

(4) |

where the holonomies and act on the target qubit conditionalized on the states and , respectively, of the control qubit. One can see that by switching the roles of control and target qubits, which can be done by swapping the encoding of states and , holonomic conditional two-qubit gate of the form can be carried out.

We now verify that is sufficient for implementing a universal set of one- and two-qubit gates. The key observation is that any triplet and can be realized since is arbitrary up to the condition . First, by choosing such that , it follows that and reduces to the one-qubit holonomic gate with the Pauli operators acting on the target qubit and a unit vector. Any given , characterized by the spherical polar angles , can for instance be realized by choosing and , such that and . By applying sequentially two such gates with different ’s corresponding to and yields , which is an arbitrary SU(2) transformation acting on the target qubit. Thus, our setup allows for arbitrary holonomic one-qubit gates. Secondly, by choosing with and implements the holonomic conditional gate that may entangle the qubits. Any entangling two-qubit gate is universal when assisted by arbitrary single-qubit transformations [20]. In other words, our scheme is universal.

Conceptually, the proposed holonomic conditional gate can be viewed as a non-Abelian version of the Abelian geometric phase shift gate developed in Ref. [21, 22, 23]. In the Abelian case, evolution is chosen such that two orthogonal qubit states and , say, pick up zero dynamical phases (for instance, by parallel transporting and ), which implies that and , being the Aharonov-Anandan (AA) geometric phase [24] with the solid angle enclosed on the Bloch sphere. This defines the AA geometric phase gate , . In our non-Abelian scheme, conditional two-qubit gate is implemented if there exist pair of loops and in the Grassmannian manifold of the target qubit along which there is no dynamical phase. Note that the orthogonal subspaces and correspond to the state space of the target qubit when the control qubit is in the state or , respectively. Zero dynamical contribution to the evolution of these two orthonormal subspaces is a signature of the geometric nature of our holonomic conditional two-qubit gate, just as zero dynamical phase in the Abelian schemes [21, 22, 23] is a signature of the geometric nature of the above AA gates.

## 3 Physical implementations

### 3.1 Coupled quantum dots

#### 3.1.1 Tight-binding model

As possible physical implementations of the holonomic gates, we here discuss two settings involving four coupled quantum dots arranged in closed ring geometry. The first setting involves the well-known two-dimensional electron gas supported, e.g., by a GaAs/AlGaAs heterostructure, depleted by a convenient set of top gate electrodes to form four quantum dots where electrons are confined. Devices of this type can be routinely fabricated nowadays with state of the art electron beam lithography [25]. In this system, electrons can tunnel from one dot to another with a tunnel coupling controlled individually by inter-dot gates. We also assume that the local potential, and hence the available single-particle energy levels in each quantum dot, can be individually manipulated with back gates. HQC is realized when this four-quantum-dot device contains just a single electron. If only one single-particle orbital in each quantum dot is relevant, the system can be described by a four-site tight-binding model, forming a closed ring. The second-quantized Hamiltonian reads

(5) |

where () are electron creation (annihilation) operators for Wannier orbitals localized at dot-site , with on-site energy and spin quantum number ; we assume periodic boundary conditions . The are nearest-neighbor hopping parameters, which can be controlled by gate voltages.

In the presence of an external magnetic field , the Hamiltonian in Eq. (5) acquires the Zeeman term , where the Bohr magneton has been absorbed into . Furthermore, the magnetic field generates a gauge-invariant modification of the tight-binding hopping parameters via the Peierls phase factors [26, 27, 28]

(6) |

where is the integral of the vector potential along the hopping path from site to site .

We fix the on-site energies by adjusting the voltage on the back gates of the device. We now assume that hopping parameters can be turned on and off with the same time-dependent scaling function by sweeping the top gates. In the invariant Fock subspace spanned by the four (spin-polarized) one-electron states , the time dependent Hamiltonian in this ordered basis has the form of Eq. (1), with the time-independent submatrix

(7) |

where we have put . Note that for a given magnetic field we must choose the complex-valued hopping parameters so that their phases add up to the total magnetic flux , where is the area enclosed by the four-dot ring. A change of these phases corresponds to a gauge transformation under which the holonomies of the two subspaces transform into , where with indicating the value of at site . This shows the holonomies associated with this field-dependent tight-binding Hamiltonian are gauge covariant quantities.

In order to take advantage of non-adiabaticity to control decoherence effects, the typical switching time should be as short as possible. At the same time, the action of the Hamiltonian must be ”adiabatic” enough to prevent transitions to higher orbital levels not included in Eq. (5). This requirement is expressed by the condition , where is the single-particle level spacing. For a quantum dot with a size of 100 nm, typically , which would permit to be as fast as 10 ps. This time is considerably shorter than the phase coherence time for charge qubits in quantum dots ( 1-10 ns [25]).

#### 3.1.2 Spin model

As the second setting to realize our scheme, we start again with the four-dot system introduced above, but we consider the case where four electrons are confined in the device. Now, electron-electron interactions must be included and we can do so by adding to the Hamiltonian of Eq. (5) an onsite Hubbard repulsion term , where . In the limit double occupancy in each dot is suppressed and the low-energy physics of the system can be described by the spin Hamiltonian [29]

(8) | |||||

where is a spin vector operator localized at site . In Eq. (8), the first term is an isotropic Heisenberg interaction, with antiferromagnetic exchange coupling , the second term is an anisotropic exchange described by second rank symmetric tensors , and the last term is an antisymmetric Dzialoshinski-Moriya (DM) interaction. The general exchange coupling arises from the interplay between the electron-electron repulsion and spin-orbit interaction. For the geometry of four coplanar quantum dots in the plane only components of the DM exchange vectors are relevant. In an in-plane electric field, the parameters in Eq. (8) can be modified so that the effective Hamiltonian of the spin-model is given by

(9) |

Here, we have with and , where is a spin-dependent hopping term describing the spin-orbit interaction.

Let be the eigenvector of the component of the spin operator at site corresponding to the eigenvalue . The four-dimensional subspace with , spanned by the ordered basis is an invariant subspace [30]. To realize the HQC scheme, the parameters and have to be turned on and off with the same time-dependent function . Under this condition, the Hamiltonian in the ordered basis has exactly the form of Eq. (1) with the time-independent submatrix

(10) |

where we have put and . It is important to stress that since both and are proportional to the inter-dot tunneling constants, they can be efficiently manipulated with time-dependent gate voltages.

### 3.2 Single-molecule magnet

The spin Hamiltonian in Eq. (8) also describes the magnetic properties of a class of single-molecule magnets (SSMs) composed of antiferromagnetic spin rings [15]. Such systems would in principle be more scalable, compact, and reproducible than the devices based on quantum dots, since their exchange constants and spin coherence properties are controlled by the molecular structure of the molecule, which can be chemically engineered [31].

The systems that we are interested in are triangular SMMs such as [32, 33]. Their magnetic core consists of three spins positioned at the vertexes of an equilateral triangle and coupled by an antiferromagnetic Heisenberg exchange interaction. The ground state (GS) manifold of these frustrated SMMs is given by two degenerate total spin doublets () of opposite spin chirality . The energy gap to the first excited quadruplet is typically in the order of 1 meV. A spin-orbit induced DM interaction lifts the degeneracy between the two chiral doublets with a splitting [34, 35]. What makes these triangular SMMs interesting for quantum manipulation is that, due to the lack of inversion symmetry, an electric field in the -plane of the molecule couples the two GS doublets of opposite chirality [34, 29, 36]. In the four-dimensional GS manifold spanned by the states , the effective low-energy spin Hamiltonian in the presence of external electric and magnetic fields is [34]

(11) |

where and are respectively the chirality and spin operators. Due to the symmetry of the molecule . The parameter is the effective electric dipole coupling, and it gives the strength of the coupling between the two states with opposite chirality brought about by the electric field. In Cu SMMs is not small [36], and for typical electric fields generated by scanning tunneling microscope ( kV/cm) the spin-chirality manipulation (Rabi) time is ps [34, 36].

To implement our holonomic scheme in this setting, we assume and are turned on and off by a common scaling function and that is negligibly small. Under these conditions, the Hamiltonian in Eq. (11) takes the form

(12) |

in the ordered basis , , , of the ground state subspace. Here, and , where we have put and , . This expression of the Hamiltonian shows that triangular SMMs such as Cu can be used to implement efficient and robust gates. Although the restricted form of associated with this system does not allow for universal computation, we would still be able to obtain some of the most important gates such as all the Pauli-matrices and the controlled-Z gate in a fashion that is highly resistant to errors. The crucial challenge relies on the ability of generating electric and magnetic fields characterized by time-dependence function . Fig. 1 shows the fidelity of the Pauli-Y gate performed with a triangular SMM, as a function of the spin-orbit induced splitting and operation time . Here, we have adopted the definition [37] , where is the ideal gate and is the gate carried out by the perturbative Hamiltonian in Eq. (12). The fidelity determines the Hilbert-Schmidt distance between and . For typical values of meV [35] the fidelity is over for ps.

It is worth mentioning that the spin model described by the Hamiltonian in Eq. (9) can be realized in other settings, such as magnetic impurities on the surface of topological insulators [38], cold atoms in a two-dimensional optical lattice [39], and nuclear spin qubits interacting via a semiconductor-heterojunction [40]. These systems may provide alternative platforms to implementing our holonomic gates.

### 3.3 Scalability

The four-level configuration presented above can be used in a scalable pattern necessary for universal quantum computation on any number of qubits. This can be achieved by implementing the Hamiltonian

(13) |

where and label the qubits. Any gate acting on qubits and can be implemented by turning on and off while the other terms are kept off. can be realized with the first two examples using the techniques introduced in Ref. [41] to design scalable architectures in quantum dot systems. This Hamiltonian would also be realized by arranging triangular SMMs on a surface in a way that each of them can be manipulated individually.

## 4 Conclusions

In conclusion, we have introduced a scheme to implement high-speed universal holonomic quantum gates that considerably differs from the other setups proposed for HQC. It allows for implementation of conditional gates using non-adiabatic quantum holonomies [7] acting on qubit pairs encoded in a four level structure. We make use of the two-qubit system as a unit cell for scalable all-geometric quantum computation. We propose three different realizations of our scheme based on an unconventional use of quantum dot and single-molecule magnet devices. We emphasize that the physical realizations of HQC presented here differ substantially from previous proposals of implementation of charge and spin qubits in quantum dots and single-molecule magnets. Finally, we would like to point out that the four level configuration proposed here is quite a general scheme that can be realized with other physical systems actively considered for quantum computer implementation [38, 39, 40].

## Acknowledgments

This work was supported by School of Computer Science, Physics and Mathematics at Linnaeus University (Sweden). E.S. thanks National Research Foundation and the Ministry of Education (Singapore) for support.

## References

- [1] Zanardi P and Rasetti M 1999 Holonomic quantum computation Phys. Lett. A 264 94
- [2] Wilczek F and Zee A 1984 Appearance of gauge structure in simple dynamical systems Phys. Rev. Lett. 52 2111
- [3] Berry M V 1984 Quantal phase factors accompanying adiabatic changes Proc. Roy. Soc. London Ser. A 392 45
- [4] Duan L M, Cirac J I and Zoller P 2001 Geometric manipulation of trapped ions for quantum computation Science 292 1695
- [5] Faoro L, Siewert J and Fazio R 2003 Non-Abelian holonomies, charge pumping, and quantum computation with Josephson junctions Phys. Rev. Lett. 90 028301
- [6] Solinas P, Zanardi P, Zanghì N and Rossi F 2003 Semiconductor-based geometrical quantum gates Phys. Rev. B 67 121307
- [7] Anandan J 1988 Non-adiabatic non-Abelian geometric phase Phys. Lett. A 133 171
- [8] Sjöqvist E, Tong D M, Andersson L M, Hessmo B and Singh K 2012 Non-adiabatic holonomic quantum computation New J. Phys. 14 103035
- [9] Johansson M, Sjöqvist E, Andersson L M, Ericsson M, Hessmo B, Singh K and Tong D M 2012 Robustness of nonadiabatic holonomic gates Phys. Rev. A 86 062322
- [10] Xu G F, Zhang J, Tong D M, Sjöqvist E and Kwek L C 2012 Nonadiabatic holonomic quantum computation in decoherence-free subspaces Phys. Rev. Lett. 109 170501
- [11] Zhang J, Kwek L C, Sjöqvist E, Tong D M and Zanardi P 2013 Quantum computation in noiseless subsystems with fast non-Abelian holonomies arXiv:1308.1919 [quant-ph]
- [12] Abdumalikov A A, Fink J M, Juliusson K, Pechal M, Berger S, Wallraff A and Filipp S 2013 Experimental realization of non-Abelian non-adiabatic geometric gates Nature 496 482
- [13] Feng G, Xu G and Long G 2013 Experimental realization of nonadiabatic holonomic quantum computation Phys. Rev. Lett. 110 190501
- [14] Kouwenhoven L P, Austing D G and Tarucha S 2001 Few-electron quantum dots Rep. Prog. Phys. 63 701
- [15] Gatteschi D, Sessoli R and Villain J 2006 Molecular Nanomagnets (Oxford University Press, Oxford)
- [16] Leuenberger M N and Loss D 2001 Quantum computing in molecular magnets Nature 410 789
- [17] Tejada J, Chudnovsky E M, del Barco E, Hernandez J M and Spiller T P 2001 Magnetic qubits as hardware for quantum computers Nanotechnology 12 181
- [18] Hanson R, Kouwenhoven L P, Petta J R, Tarucha S and Vandersypen L M K 2007 Spins in few-electron quantum dots Rev. Mod. Phys. 79 1217
- [19] Żac R A, Röthlisberger B, Chesi S and Loss D 2010 Quantum computing with electron spins in quantum dots Riv. Nuovo Cimento 033 345
- [20] Bremner M J, Dawson C M, Dodd J L, Gilchrist A, HarrowA W, Mortimer D, Nielsen M A and Osborne T J 2002 Practical scheme for quantum computation with any two-qubit entangling gate Phys. Rev. Lett. 89 247902
- [21] Wang X B and Keiji M 2001 Phys Rev. Lett. 87 097901
- [22] Zhu S-L and Wang Z D 2002 Implementation of universal quantum gates based on nonadiabatic geometric phases Phys. Rev. Lett. 89 097902
- [23] Zhu S-L and Wang Z D 2003 Universal quantum gates based on a pair of orthogonal cyclic states: application to NMR systems Phys. Rev. A 67 022319
- [24] Aharonov Y and Anandan J 1987 Phase change during a cyclic quantum evolution Phys. Rev. Lett. 58 1593
- [25] Petersson K D, Petta J R, Lu H and Gossard A C 2010 Quantum coherence in a one-electron semiconductor charge qubit Phys. Rev. Lett. 105 246804
- [26] Graf M and Vogl P 1995 Electromagnetic fields and dielectric response in empirical tight-binding theory Phys. Rev. B 51 4940
- [27] Ismail-Beigi S, Chang E K and Louie S G 2001 Coupling of nonlocal potentials to electromagnetic fields Phys. Rev. Lett. 87 087402
- [28] Cehovin A, Canali C M and MacDonald A H 2004 Orbital and spin contributions to the g tensors in metal nanoparticles Phys. Rev. B 69 045411
- [29] Trif M, Troiani F, Stepanenko D and Loss D 2010 Spin electric effects in molecular antiferromagnets Phys. Rev. B 82 045429
- [30] Note that this subspace involves the excited states with of the spin Hamiltonian.
- [31] Wedge C J, Timco G A, Spielberg E T, George R E, Tuna F, Rigby S, McInnes E J L, Winpenny R E P, Blundell S J and Ardavan A 2012 Chemical engineering of molecular qubits Phys. Rev. Lett. 108 107204
- [32] Kortz U, Al-Kassem N K, Savelieff M G, Al Kadi N A and Sadakane M 2001 Synthesis and characterization of copper-, zinc-, manganese-, and cobalt-substituted dimeric heteropolyanions, [(-XWO) M (HO)] (n = 12, X = As, Sb, M = Cu, Zn; n = 10, X = Se, Te, M = Cu) and [(-As W O) WO (HO) M (H O)] (M = Zn, Mn, Co) Inorg. Chem. 40 4742
- [33] Choi K-Y, Matsuda Y H, Nojiri H, Kortz U, Hussain F, Stowe A C, Ramsey C and Dalal N S 2006 Observation of a Half step magnetization in the Cu-type triangular spin ring Phys. Rev. Lett. 96 107202
- [34] Trif M, Troiani F, Stepanenko D and D. Loss D 2008 Spin-electric coupling in molecular magnets Phys. Rev. Lett. 101 217201
- [35] Nossa J F, Islam M F, Canali C M and Pederson M R 2012 First-principles studies of spin-orbit and Dzyaloshinskii-Moriya interactions in the Cu single-molecule magnet Phys. Rev. B 85 085427
- [36] Islam M F, Nossa J F, Canali C M and Pederson M R 2010 First-principles study of spin-electric coupling in a Cu single molecular magnet Phys. Rev. B 82 155446
- [37] Our fidelity is the linear function of the scalar cost function introduced in [Cabrera R, Shir O M, Wu R and Rabitz H 2011 Fidelity between unitary operators and the generation of robust gates against off-resonance perturbations J. Phys. A: Math. Theor. 44 095302].
- [38] Zhu J-J, Yao D-X, Zhang S-C and Chang K 2011 Electrically controllable surface magnetism on the surface of topological insulators Phys. Rev. Lett. 106 097201
- [39] Radić J, Di Ciolo A, Sun K and Galitski V 2012 Exotic quantum spin models in spin-orbit-coupled Mott insulators Phys. Rev. Lett. 109 085303
- [40] Mozyrsky D, Privman V and Glasser M L 2001 Indirect interaction of solid-state qubits via two-dimensional electron gas Phys. Rev. Lett. 86 5112
- [41] Trifunovic L, Dial O, Trif M, Wootton J R, Abebe R, Yacoby A and Loss D 2012 Long-distance spin-spin coupling via floating gates Phys. Rev. X 2 011006