Fast non-Abelian geometric gates via transitionless quantum driving.

A practical quantum computer must be capable of performing high fidelity quantum gates on a set of quantum bits (qubits). In the presence of noise, the realization of such gates poses daunting challenges. Geometric phases, which possess intrinsic noise-tolerant features, hold the promise for performing robust quantum computation. In particular, quantum holonomies, i.e., non-Abelian geometric phases, naturally lead to universal quantum computation due to their non-commutativity. Although quantum gates based on adiabatic holonomies have already been proposed, the slow evolution eventually compromises qubit coherence and computational power. Here, we propose a general approach to speed up an implementation of adiabatic holonomic gates by using transitionless driving techniques and show how such a universal set of fast geometric quantum gates in a superconducting circuit architecture can be obtained in an all-geometric approach. Compared with standard non-adiabatic holonomic quantum computation, the holonomies obtained in our approach tends asymptotically to those of the adiabatic approach in the long run-time limit and thus might open up a new horizon for realizing a practical quantum computer.

Fast and robust quantum gates play a central role in realizing a practical quantum computer. While the robustness offers resilience to certain errors such as parameter fluctuations, the fast implementation of designated quantum gates increases computational speed, which in turn decreases environment-induced errors. A possible approach towards robust quantum computation is to implement quantum gates by means of different types of geometric phases1234; an approach known as holonomic quantum computation (HQC)5678910. Such geometric gates depend solely on the path of a system evolution, rather than its dynamical details.

Universal quantum computation based purely on geometric means has been proposed in the adiabatic regime, resulting in a precise control of a quantum-mechanical system7. Despite the appealing features, the adiabatic evolution is associated with long run time, which increases the exposure to detrimental decoherence and noise. However, this drawback can be eliminated by using non-adiabatic HQC schemes based on Abelian89 or non-Abelian geometric phases10. The latter has been developed further in refs 11, 12, 13, 14 experimentally demonstrated in refs 15, 16, 17, 18 and its robustness to a variety of errors has been studied in refs 19,20.

Adiabatic processes can also be carried out swiftly by employing transitionless quantum driving algorithm (TQDA) if the quantum system consists of non-degenerate subspaces21. This is also known as adiabatic shortcut in the literature222324252627. A key notion of TQDA is to seek a transitionless Hamiltonian so that the system evolves exactly along the same adiabatic passage of a given target Hamiltonian, but at any desired rate. This is achieved with the aid of an additional Hamiltonian that suppresses the energy level fluctuations caused by the changes in the system parameters.

In this report, we generalize TQDA to degenerate subspaces, where non-Abelian geometric phases are acquired after a cyclic evolution. With the help of the generalized TQDA, we propose a universal set of non-adiabatic holonomic single- and two-qubit gates. Specifically, non-Abelian geometric phases or quantum holonomies are acquired by a degenerate subspace after a cyclic evolution. TQDA-based geometric phases are realized via non-adiabatic evolution, dictated by an additional transition-suppressing Hamiltonian. We further simplify the transitionless Hamiltonian by selectively choosing geodesic path segments forming a loop in the system parameter space. This effectively reduces the complexity of the control parameters, which in turn simplifies the physical implementation of the TQDA-gates. Finally, we show how to realize the effective Hamiltonians in superconducting tunable coupling transmons.

Results

Degenerate TQDA

Let us consider a quantum-mechanical system characterized by a Hamiltonian that depends on a set of control parameters . For an eigenvalue of the Hamiltonian, we assume there exists a set of eigenstates that span the proper subspace of Hilbert space . When the parameter is varied from an initial value to a final value , the adiabatic theorem entails that does not mix with other subspaces in when the corresponding run time tends to infinity. The states satisfying the Schrödinger equation along the path in parameter space are related to as with being unitary matrices. By substituting this relation into the Schrödinger equation, we arrive at

up to a global dynamical phase. Here, is a time-ordering operator, the connection is an anti-Hermitian matrix, and we have assumed the initial condition , being the unit matrix. We focus on a cyclic adiabatic evolution, i.e., . A non-trivial transformation relating the initial and final states is a quantum gate. This transformation operator, be it Abelian or non-Abelian , is a quantum holonomy.

The transitionless Hamiltonian that exactly generates the adiabatic time evolution governed by the unitary operator , can be obtained by reverse engineering21 yielding

Since we have made no assumptions regarding the dimensionality of the subspaces, this result generalizes Berry’s TQDA21 to the degenerate case. From this point onwards, we denote as the target Hamiltonian. One notes that the expectation values and coincide, i.e., the dynamical phases acquired are the same, although the eigenstates of are not eigenstates of the transitionless Hamiltonian H. In cyclic evolution, a non-Abelian geometric phase, induced by the matrix-valued connection appearing in the additional transition-suppressing Hamiltonian, is picked up by each subspace. In the following, we demonstrate how to arrive at desired quantum gates from non-Abelian geometric phases under the degenerate TQDA.

Fast holonomic single-qubit gates

To implement our proposed holonomic gates, we consider a generic system consisting of four bare energy eigenstates , and , coupled by external oscillating fields in a tripod configuration7. Here, are qubit states; and are excited and auxiliary states, respectively. By assuming that all fields are resonant with the transitions, the target Hamiltonian of such a system is written as

where each serves as a time-dependent control parameter for the transition in a frame that rotates with the fields. has two dark and two bright eigenstates, and the former eigenstates provide a proper subspace to encode a qubit. To obtain our desired holonomic gates, we evolve the subspace along some loops in parameter space. In the standard adiabatic scheme7, this task is accomplished by changing the control parameters slowly so that transitions between the dark and bright subspaces become negligible. In the following, we show explicitly how to speed up this process with the help of degenerate TQDA. Since the qubit information is stored in the dark subspace, there is no dynamical contribution. Therefore, the desired quantum holonomy is acquired at the end of the evolution in the same manner as in ref. 7.

We first construct a phase-shift gate , where the phase depends solely on the path in the parameter space defined by the f’s. This holonomic gate is implemented adiabatically with the Λ-like Hamiltonian , where the state is decoupled . The Hamiltonian has a parameter-dependent dark eigenstate , and two bright eigenstates. The dark state initially coincides with by choosing at . After completing a cyclic adiabatic evolution, this qubit state acquires the Berry phase factor , with S being the surface enclosed by the path shown in Fig. 1(a). Thus, is minus half the solid angle enclosed on the parameter sphere with polar angles θ and ϕ.

Figure 1

Schematic diagram for the orange slice and geodesic triangle schemes.

The spheres in (a,c) represent the parameter space, where the control parameters θ and are changed following a specific path. When we set , the starting point in the parameter space is along the z axis. The parameters θ and are polar angle and azimuthal angle in the parameter space, respectively. (a) Schematic diagram for the orange slice scheme. The state is evolved by changing the parameters along the path A-B-C-D-A. The opening angle is . (b) The level structure for the first and third steps in achieving . Transitions between all three levels are needed, which makes it a Δ model Hamiltonian. (c) Schematic diagram for the geodesic triangular scheme. The manifold, consisting of and , is evolved by changing the parameters along the path A-B-C-A, while the opening angle is . (d) The level structure for the last step in achieving . Transitions between all four levels are needed.

The shortcut to this adiabatic process is realized by adding an extra Hamiltonian term to the target Hamiltonian . It is evident that relies both on the structure and rate of change of . In order to simplify , we propose evolving the system in three steps forming an “orange slice” path on the parameter sphere along which the connection A vanishes (see Fig. 1(a,b) and Supplementary information for details). First, while changing θ from 0 to π. During this step, terms containing and vanish, and becomes a Δ-like Hamiltonian, in which the Rabi frequency between and is . Therefore, a transition is added to the system caused by the temporal change in θ. In the second step, we keep , but vary from 0 to . Here, the total Hamiltonian is just a transition between and with detuning . Thus, the detuning is caused by the temporal change in the parameter . Finally, we keep , while θ is decreased to 0. During this final step, is again a Δ-like Hamiltonian. When , has returned to its initial form. As shown in Fig. 1(a), the acquired geometric phase is equal to minus the azimuthal angle , i.e., minus half the solid angle enclosed by the orange slice shaped loop. We note that both the detuning in the second step and the transition in the first and third steps are caused by the fast parameter change. Thus, is recovered in the adiabatic limit when the rate of the parameter change tends to zero. The steps of the control Hamiltonian needed to implement are shown in Table 1.

Table 1

Transitions needed in the “orange slice” scheme and the “geodesic triangle” scheme.

	Step1	Step2	Step3
UP
UB

The “orange slice” scheme and the “geodesic triangle” scheme are proposed to achieve the single qubit gates and , respectively. Each scheme contains three steps and needed transitions for each step are listed. The symbol denotes the transition between the states and .

Similarly, an adiabatic gate , can be realized by utilizing the full tripod structure with , , and . This Hamiltonian has two degenerate parameter-dependent dark eigenstates, which we take as and , together with two non-degenerate bright eigenstates. Similar to the adiabatic phase-shift gate, our qubit is encoded within the dark subspace. The desired holonomic gate is obtained after traversing a loop on the parameter sphere with polar angles θ and , resulting in as the solid angle swept7.

Shortcut to is realized by choosing a geodesic triangle path on the corresponding parameter sphere, again in order to force the connection A to vanish (see Fig 1(c,d) and Supplementary information for details). The purpose of this choice is to simplify the implementation of the gate, just as the choice of an “orange slice” path simplified the implementation of . Initially, we set , which implies and , followed by increasing θ to , while keeping constant. During this step, exhibits a Δ-like structure involving the states , and . Thereafter, we keep θ unchanged but rotate to . has again a Δ-like structure, but now involving , and . Finally, θ is decreased to zero along the geodesic curve and then is tuned to zero. The level structure for the last step is depicted in Fig. 1(d). Unlike the “orange slice” scheme, no detuning emerges in the “geodesic triangle” scheme. Therefore, all the photon-assisted transitions are resonant. The steps of the control Hamiltonian needed to implement are shown in Table. 1. We may note that all transitions between the four levels , and are needed to realize an arbitrary single-qubit gate by combining and .

Fast holonomic two-qubit gate

Since the combination of the gates and allows the realization of an arbitrary single-qubit gate, it remains to construct an entangling two-qubit holonomic gate to realize universal quantum computation2829. A holonomic two-qubit gate can be realized in the TQDA scenario, by controlling suitable coupling parameters between two four-level tripod systems. Specifically, we note that there are four auxiliary states, from the two four-level systems (see Methods and Supplementary information for details). The key idea is to construct a Δ-like structure to combine two auxiliary states from and one of the two-qubit states. The target Hamiltonian of the two-qubit gate can be taken as   By using the same “orange slice” scheme as in the realization of , a geometric phase is acquired by the state while the other states remain unchanged. This amounts to applying an entangling gate on the two-qubit space, where is the geometric phase obtained after the cyclic evolution.

The single- and two-qubit gates in superconducting circuits

To realize our non-adiabatic schemes experimentally, we consider a superconducting tunable coupling transmon (TCT)30 (see Fig. 2(a)), in which we assume that the transition frequency of the system and system-resonator coupling strength can be tuned independently. In the symmetric configuration, where , , the system Hamiltonian is described by3031

Figure 2

Diagram for superconducting circuits.

(a) Circuit model of a tunable coupling transmon (TCT) system, which is used to realize our single-qubit gates in absence of a transmission line resonator. (b) A proposed setup to realize the two-qubit gate with two TCTs mediated via a resonator and (c) its equivalent circuit model, where the two TCTs are capacitively coupled to a coplanar resonator.

where are the charging energies of the upper and lower islands, while , , are the Josephson energies (see Methods for definition of the other parameters in ). In the limit of large , the energy level splitting is observed with increase in as shown in Fig. 3(a), giving rise to the four energy levels and needed to implement our proposal.

Figure 3

Energy spectrum and transition matrix for the superconducting circuits.

(a) Energy spectrum of the TCT as a function of for . We confine ourselves to the eigenstates and . (b) Moduli of the unnormalized transition matrix elements between states and as a function of for the labelled eigenstates in (b), while and .

When a time-dependent microwave field impinges onto our four-level TCT (see Figs 2(a) and 3(a) and Methods), microwave-assisted transitions occur. To investigate allowed energy level transitions, we consider a situation where there is a minimum flux entering the upper SQUID loop , where a small change of flux is present in the lower one. In this case, the perturbation Hamiltonian with small becomes   with the microwave amplitude, a carrier frequency of the incoming microwave pulse that connects two desired eigenstates and , and . In Fig. 3(b), we plot the transition matrix elements for the eigenstates of interest versus while and . From this result, we conclude that a particular energy value of should be tuned via the external magnetic flux, threading the lower SQUID loop, in order to give rise to the desired energy level transitions. If we confine ourselves to the selected four eigenstates, the Hamiltonian of the superconducting circuit in a rotating frame can be expressed as , where is a complex Rabi frequency. When we invoke the “orange slice” and “geodesic triangle” schemes, gives rise to and , respectively. We note that since our schemes relax the adiabatic constraint, the holonomic gates can be obtained at any speed. However, from a practical perspective, the desired gate time depends on details of the specific experimental setup such as the energy gaps and how fast we can vary the input microwave pulses.

To implement our two-qubit gate, we consider a superconducting architecture with two TCTs placed inside a coplanar resonator (see Fig. 2(b,c)). The direct coupling between the two transmons (labeled as l and r) is mediated by a cavity in the resonant regime32. With the built-in tunability of the transmons, we realize the Hamiltonian as follow. We assume that the transmon levels are anharmonic and each level transition can be tuned at will by matching with the resonator frequency via ac-Stark-shift fields33, alternatively by invoking the built-in tunability of the transmons energy levels30 or by tuning the resonator frequnecy34. For instance, to achieve the transition between and , we switch on vacuum Rabi couplings to transmon l, and to transmon r, where the g’s are the coupling strengths between the respective transition and the single mode resonator, , , and is the creation (annihilation) operator of the bosonic resonator mode, respectively. All the other required transitions can also be realized in the same manner (see Supplementary information for details).

Discussion and Conclusions

We would like to put forward some remarks on the holonomic gates based on TQDA proposed here and the nonadiabatic holonomic gates proposed in ref. 10. On one hand, we notice that both schemes are genuinely non-adiabatic in that they can be performed at any desired rate. On the other hand, the holonomies obtained in the proposed TQDA approach tends asymptotically to those of the adiabatic approach in the long run-time limit27; a feature which is not shared by the scheme in ref. 10. As a consequence, the TQDA-based gates have an immediate geometric meaning in terms of solid angles enclosed in parameter space of the underlying target Hamiltonian, while no such explicit geometric meaning is present in the other scheme. This also makes our scheme explicitly robust to parameter fluctuations that preserve the solid angle. In addition, our scheme has the advantage over standard adiabatic scheme in that non-adiabatic transitions, which inevitably degrades the control in any realistic finite-time implementation of the adiabatic gates19, are suppressed using TQDA.

In conclusion, we have proposed a non-Abelian generalization of the TQDA and show how it can be used for holonomic quantum computation. The key feature of this algorithm is to realize quantum holonomies that can be performed at arbitrary rate by applying an additional transition-suppressing Hamiltonian so that the path of the original eigensubspaces as well as the purely geometric nature of the resulting matrix-valued phases are preserved. Non-adiabatic TQDA-based single- and two-qubit gates are proposed in a four-level transmon and two four-level transmons coupled with a cavity, respectively. With high controllability3536 and scalability3738 of the superconducting circuit, our scheme can also be incorporated within a more general landscape where the transmons are embedded in a resonator lattice3940414243. This scheme therefore opens up new experimental horizons towards a robust all-geometric high-speed large-scale quantum computation architecture.

Methods

Superconducting tunable coupling transmon

The four level system can be provided by a superconducting tunable coupling transmon (TCT)30 (see Fig. 2(a)). To obtain the effective Hamiltonian , we define , being the superconducting magnetic flux quantum) as the phase differences on the upper and lower superconducting quantum interference device (SQUID) loops, and as the dimensionless gate voltages. The tunable Josephson energies , where and the frustration parameters . are external magnetic fluxes threading the upper and lower superconducting quantum interference devices (SQUIDs). In , is the respective associated charge. are the charging energies of the upper and lower islands, while represents the interaction energy between them. By letting , the effective capacitors and are defined by and , respectively.

Additional Information

How to cite this article: Zhang, J. et al. Fast non-Abelian geometric gates via transitionless quantum driving. Sci. Rep. 5, 18414; doi: 10.1038/srep18414 (2015).