Frontiers journals are at the top of citation and impact metrics

Original Research ARTICLE

Front. ICT, 15 September 2016 | https://doi.org/10.3389/fict.2016.00019

Energetic Cost of Superadiabatic Quantum Computation

  • 1Instituto de Física, Universidade Federal Fluminense, Niterói, Brazil
  • 2Information Sciences Institute, University of Southern California, Marina del Rey, CA, USA
  • 3Ming Hsieh Department of Electrical Engineering, Center for Quantum Information Science and Technology, University of Southern California, Los Angeles, CA, USA

We discuss the energetic cost of superadiabatic models of quantum computation. Specifically, we investigate the energy–time complementarity in general transitionless controlled evolutions and in shortcuts to the adiabatic quantum search over an unstructured list. We show that the additional energy resources required by superadiabaticity for arbitrary controlled evolutions can be minimized by using probabilistic dynamics, so that the optimal success probability is fixed by the choice of the evolution time. In the case of analog quantum search, we show that the superadiabatic approach induces a non-oracular counter-diabatic Hamiltonian, with the same energy–time complexity as equivalent adiabatic implementations.

1. Introduction

Shortcuts to adiabatic passage (Demirplak and Rice, 2003, 2005; Berry, 2009; Torrontegui et al., 2013) provide a remarkable mechanism for speeding up quantum tasks, which can be achieved through the use of a counter-diabatic assistant driving. These techniques have been introduced to mimic the transitionless adiabatic dynamics, but with the usual constraint on the adiabatic runtime lifted. Transitionless quantum driving has been applied to a number of quantum information protocols, such as population transfer (Chen et al., 2014a; Lu et al., 2014a) and entanglement generation (Chen et al., 2014b, 2015a,b; Lu et al., 2014b). In the context of many-body systems, realizable settings have been investigated for assisted evolutions in quantum critical phenomena (del Campo et al., 2012; del Campo, 2013; Saberi et al., 2014). More recently, counter-diabatic approaches have been proposed for fast implementation of individual unitaries in quantum circuits, leading to universal superadiabatic schemes of quantum computing (QC) via local Hamiltonians (Santos and Sarandy, 2015; Santos et al., 2016). Such methods may be potentially relevant to accelerating the implementation of n-qubit controlled gates in digitized proposals of adiabatic quantum computing [see Kieferová and Wiebe (2014), Martinis and Geller (2014), Barends et al. (2015), and Hen (2015)].

The superadiabatic speedup is intrinsically connected with an increase of the energy resources demanded by the quantum computer (Santos and Sarandy, 2015; Santos et al., 2016), which in turn implies a rather versatile computational cost that is controlled by the energetic capacity available to the physical apparatus. Here, we show that this energy–time complementarity can be exploited in quantum information processing. First, we consider controlled evolutions (CE) as a mechanism to implement superadiabatic universal QC (Santos and Sarandy, 2015), which generalizes the original adiabatic approach introduced in Hen (2015). We then show that, within the superadiabatic scenario, the energetic cost can be minimized by replacing the deterministic realization of quantum gates for probabilistic implementations based on a probability distribution of a binary random variable described by an angle parameter. By doing so, the energy expense can be minimized by adjusting the probability distribution, provided the choice of the evolution time of the computational process. Second, we analyze the effects of the energy–time complementarity in analog quantum search (Grover, 1997), where the oracular approach designed by the local adiabatic Grover algorithm is known to be optimal (van Dam et al., 2001; Roland and Cerf, 2002). In this case, we show that the superadiabatic approach naturally requires an unphysical non-oracular counter-diabatic Hamiltonian, with the energy–time complexity equivalent to non-oracular adiabatic implementations.

The paper is organized as follows. In Section 2, we describe the adiabatic implementation of quantum gates via CE and several adiabatic quantum search approaches. We then provide their superadiabatic versions and introduce the metric for energetic cost used in our work. In Section 3, we investigate the energy complexity of the superadiabatic realizations of both quantum gates via CE and analog quantum search. In particular, we consider the properties of the probabilistic model of QC through CE and the consequences of the energy–time complementarity for the search problem. In Section 4, we present our conclusions and future perspectives.

2. Materials and Methods

Our aim in this Section is to discuss adiabatic implementations of QC, their superadiabatic generalizations, and the energetic cost measure adopted in this work.

2.1. Quantum Gates by Adiabatic Controlled Evolutions

Let us begin by using adiabatic CE (Hen, 2015) to implement n-controlled gates (Santos and Sarandy, 2015). To this end, we will consider the adiabatic evolution of a composite system 𝒯 𝒜 associated with a Hilbert space ℋ𝒯 ⊗ ℋ𝒜, where 𝒯 denotes a target subsystem containing n + 1 qubits and 𝒜 denotes an auxiliary subsystem containing a single qubit. We will use the first n qubits of 𝒯 as the control register of the n-controlled gate, while the last qubit will play the role of its target register. Then, a rotation of the target qubit of an angle ϕ around a direction n^ in the Bloch sphere will be performed when the state of the control register is |11⋯ 1⟩. We will adopt here the decimal representation |11⋯ 1⟩ ≡ |N − 1⟩, with N = 2n. An n-controlled rotation over a single qubit can be adibatically implemented by preparing the auxiliary qubit in the initial state |0⟩, with the adiabatic Hamiltonian given by (Santos and Sarandy, 2015)

H(s)=𝟙PN1,nH0(s)+PN1,nHϕ(s),

where Pk,n±=|kk||n^±n^±| is the set of all orthogonal projectors on the subspace 𝒯 and |n^±n^±|=12(𝟙±n^σ), with σ=σx,σy,σz. The Hamiltonians H0(s) and Hϕ(s) are given by

Hξ(s)=ωσzcosθ(s)+sinθ(s)[σxcosξ+σysinξ],

where θ (s) = θ0s, θ0 is a constant angle, ξ = {0, ϕ}, and s denotes the normalized time s = t/τ, with τ the total evolution time. The system is prepared in the initial state |Ψ(0)⟩ = |ψn⟩ ⊗ |0⟩, where

|ψn=m=0N1ϵ=±γm,ϵ|m,n^ϵ.

Then, by adiabatic evolution, the system will evolve to the final state |Ψ(1)⟩ given by

|Ψ(1)=𝟙PN1,n^ψn|E00(1)+PN1,n^ψn|Eϕ0(1),

where |Eξ0(1)=cosθ020+eiξsinθ021 is the ground state of Hξ(1). Then, equivalently, we can write

|Ψ(1)=cosθ02ψn|0+sinθ02|ψnrot1,

with

|ψnrot=k=0N2ϵ=±γk,ϵk,n^ϵ+N1[γN1,+n^++eiϕγN1,n^].

The rotated state |ψnrot is the target of the n-controlled gate. However, note that |Ψ(1)⟩ in equation (5) is an entangled state. Thus, a measurement must be performed on the auxiliary system, where the action of the gate will be considered successful if 𝒜 is measured in the state |1⟩, which occurs with probability sin2 (θ0/2). On the other hand, if the outcome of a measurement on 𝒜 yields |0⟩, the adiabatic evolution should be restarted through the Hamiltonian in equation (1), as the state of the system is projected onto the initial state |Ψ(0)⟩. Naturally, by choosing θ0 = π, we deterministically ensure the success of the computation. However, as we will show, deterministic evolutions may demand more energy resources than probabilistic processes when transitionless drivings are considered. In particular, observe also that the scheme presented here allows for the implementation of arbitrary n-controlled gates, which lead to versatile sets of universal gates, e.g., single qubit rotations and controlled-NOT operations (Nielsen and Chuang, 2000).

2.2. Adiabatic Quantum Search

Instead of adiabatic implementations of quantum circuits, we can also consider the original approach of adiabatic QC (Farhi et al., 2001), where a single annealing process is performed using energy penalties attributed to quantum states that violate the solutions of an optimization problem. Here, we employ this method to analyze three possible adiabatic implementations of quantum search over an unstructured list. An adiabatic QC approach for the quantum search through Grover’s algorithm (Grover, 1997) was first proposed in Farhi et al. (2000) and improved by using local adiabaticity (van Dam et al., 2001; Roland and Cerf, 2002), where the adiabatic evolution is required for each local time interval, instead of being globally applied as in the original proposal. In both cases, the search for a marked element in an unstructured list of N = 2n elements (labeled by n qubits) can be achieved by employing a Hamiltonian of the form

H0(s)=f(s)(𝟙|++|)+g(s)(𝟙|mm|),

where |m⟩ is the marked state, s is the normalized time (0 ≤ s ≤ 1), |+=1Ni=0N1|i, and f (0) = g(1) = 1 and f (1) = g(0) = 0. The eigenspectrum of this Hamiltonian can be exactly derived [see Das et al. (2003) and Orus and Latorre (2004)]. In particular, the two lowest eigenstates can be written as

|E±(s)=N±(s)|m+b±(s)|ϕ,

where the normalization constant is N±(s)=11+(N1)b±(s)2, |ϕ=im|i, and

b±(s)=1E±(s)f(s)N¯,

with N¯=11N, and the corresponding energies E±(s) given by

E±(s)=f(s)+g(s)±[f(s)+g(s)]24f(s)g(s)N¯2.

The other higher-energy eigenstates form an (N − 2)-fold degenerate eigenspace, whose energy is given by

Edeg=f(s)+g(s).

In order to explicitly provide the eigenstates |Edegk (k = 1, ⋯ , N − 2) associated with the eigenenergy Edeg, we write

|Edegk=n=0N1cnk|n.

Then, from the eigenvalue equation for H0(s), it directly follows that the set {cnk} is just required to satisfy the constraints n=0N1cnk=0 and cmk=0. As a consequence, the states |Edegk can be suitably chosen as time-independent vectors.

By imposing a local adiabatic evolution (van Dam et al., 2001; Roland and Cerf, 2002), i.e., by requiring adiabaticity at each infinitesimal time interval, the runtime is minimized for the path [see also Kieferová and Wiebe (2014)]

f(s)=1g(s),g(s)=N1tanarctanN112s2N1.

This results in a quadratic speedup over the classical search, i.e., we obtain the time complexity O(N) expected by the Grover quantum search (van Dam et al., 2001; Roland and Cerf, 2002).

It is possible to reduce the time complexity of the Grover quantum search by transferring the algorithmic cost to other physical resources. The second implementation of the adiabatic Grover search considered here has been introduced in Das et al. (2003) and Wen and Qiu (2008). It is also based on the Hamiltonian in equation (7) to perform the evolution, but requiring that the functions f (s) and g(s) satisfy

f(s)=1s+N(1s)s,
g(s)=s+N(1s)s.

This implementation achieves the solution at constant time complexity O(1). As is apparent from equations (14) to (15), the original time resource has been transferred to the coupling strengths f (s) and g(s) and, as discussed in detail in the next Section, will be reflected in the energy scaling required by the system.

The two previous versions of the adiabatic Grover’s algorithm are based on oracular Hamiltonians, which we take here to be operators able to recognize the correct answer of a problem (Nielsen and Chuang, 2000). This is indeed the case if one chooses a Hamiltonian composed of an operator Om in the form Om = 𝟙 − |m⟩ ⟨m|. The action of Om in the computational basis {|i⟩} is

Om|i=(𝟙|mm|)|i=0(i=m),|i(im),

so that this operator recognizes the marked state, providing no hint about its identity if acting upon any other state. Adiabatic versions of the quantum search have also been proposed via non-oracular Hamiltonians. Our third implementation of Grover’s algorithm is based on the non-linear non-oracular (NLNO) Hamiltonian proposed in Wen et al. (2009). In this work, the time-dependent Hamiltonian in equation (7) is replaced for

H0(s)=f(s)(𝟙|++|)+g(s)(𝟙|mm|)+h(s)(|+m|+|m+|),

where h(0) = h(1) = 0. The Hamiltonian in equation (17) contains an operator O¯m=|+m|+|m+|. The action of O¯m in the computational basis {|i⟩} is

O¯m|i=(|+m|+|m+|)|i=1N|m+|+(i=m),1N|m(im).

Observe that equation (18) implies that O¯m cannot exactly recognize a marked element, even though it could effectively recover the marked state for N ≫ 1 with a single operation over the uniform superposition provided by the state |+⟩. Naturally, the non-oracular form of the Hamiltonian involves all the individual computational states, requiring, therefore, much more than the capacity of the Hamiltonian to recognize the marked element. This is an obviously artificial approach, whose discussion here is kept just for comparison with the superadiabatic scenario. Assuming a restricted feasibility of such a Hamiltonian, we proceed by looking at its eigenspectrum. The ground and first excited states have the same structure as in equation (8), with

b±(s)=N¯f(s)+2h(s)NE±(s)N¯f(s)h(s)N.

The two lowest energy levels are given by

E±(s)=12f(s)+g(s)+2h(s)N±f(s)+g(s)24f(s)g(s)N¯+4h2(s)4h(s)Nf(s)+g(s).

As before, the higher-energy states form an (N − 2)-fold degenerate subspace, with energy given by f(s) + g(s). As shown in Wen et al. (2009), this formulation also shows constant time complexity O(1), which can be obtained by choosing a suitable interpolation, such as

f(s)=1s,g(s)=s,h(s)=s(1s).

2.3. Speeding up Adiabaticity through Superadiabatic Evolutions

The performance of adiabatic QC is dictated by a long total evolution time compared to the inverse of a power of the energy gap (Messiah, 1962; Teufel, 2003; Sarandy et al., 2004; Jansen et al., 2007). However, the adiabatic evolution can be sped up through shortcuts to adiabaticity via counter-diabatic Hamiltonians (Demirplak and Rice, 2003, 2005; Berry, 2009). The fundamental idea underlying these shortcuts to adiabaticity is to add a new contribution HCD(t), called counter-diabatic Hamiltonian, to the original adiabatic Hamiltonian H(t). This term is constructed such that it allows the mimicking of the adiabatic evolution, however, without any constraint on the total time of evolution. The total composite Hamiltonian is

HSA(t)=H(t)+HCD(t),

which is called superadiabatic Hamiltonian. In particular, it is possible to show that the counter-diabatic term reads (Berry, 2009).

HCD(t)=iℏnn˙(t)n(t)+n˙(t)|n(t)n(t)n(t),

where |n(t)⟩ is the eigenstate of H(t) associated with the energy En(t). The goal of the counter-diabatic term HCD(t) in the Hamiltonian HSA(t) is exactly to eliminate the diabatic contributions of H(t). Thus, if the system is initially prepared in the ground state of H(0), then the system will deterministically evolve to the instantaneous ground state of the Hamiltonian H(t) with no constraints over the evolution time. Note that, in general, one would need to be able to explicitly calculate all the eigenstates of H(t) to derive a shortcut to adiabaticity using the counter-diabatic driving. However, this may not be a hard requirement in the case of superadiabatic versions of circuit implementations, where one-qubit rotations and two-qubit entangling gates are enough to achieve QC universality (Nielsen and Chuang, 2000). In particular, as we shall see for this case, HCD(t) can be realized through a simple time-independent operator.

2.4. Energetic Cost of Quantum Evolutions

To quantify the expense of energy in a quantum evolution driven by a Hamiltonian H(t), we adopt as the cost measure the average norm of H(t) computed for a total time of evolution τ. This yields (Kieferová and Wiebe, 2014; Santos and Sarandy, 2015; Zheng et al., 2015; Santos et al., 2016)

Σ(τ)=1τ0τH(t)dt=01H(s)ds,

where s = t/τ is the parameterized time and the norm here is defined by the Frobenius norm (Hilbert–Schmidt norm) A=TrAA. Naturally, other norms can be adopted as, for instance, the spectral norm A2=λmaxAA, where λmax denotes the maximum eigenvalue of [A A]. For the Hamiltonians investigated in this work, these norms will imply into a cost simply related by a constant D1/2, with D denoting the dimension of corresponding the Hilbert space. The Frobenius norm as well as arbitrary superadiabatic evolutions with total evolution time τ, the energetic cost can be written as

ΣSA(τ)=1τ0τTrHSA2(t)dt=1τ0τTrH2(t)+HCD2(t)dt,

where we have used that Tr({H(t), HCD(t)}) = 0 (Santos and Sarandy, 2015). This explicitly shows that a superadiabatic evolution has an energetic cost larger than its corresponding adiabatic evolution. By evaluating the trace in equation (25), we obtain

ΣSA(τ)=01mEm2(s)+2μm(s)τ2ds,

where μm(s) = ⟨∂sm(s)|∂sm(s)⟩ − |⟨m(s)|∂sm(s)⟩|2 and {Em(s)} is the energy spectrum of the adiabatic Hamiltonian H(t), with {|m(s)⟩} denoting its eigenbasis. Notice that the adiabatic limit is recovered when taking τ → ∞. Thus, the speedup obtained by the superadiabatic dynamics is limited by the energetic cost of the evolution. Indeed, this energy–time complementarity can be formally discussed through the quantum speed limit (Deffner and Lutz, 2013), which suggests that the superadiabatic evolution time is compatible with arbitrarily short time intervals (implying into corresponding arbitrarily large energies) (Santos and Sarandy, 2015), while the adiabatic evolution time obeys the lower bound τAd ∝ 1/ωn, with ω associated with the energy gap and n ∈ ℕ+ (Messiah, 1962; Teufel, 2003; Sarandy et al., 2004; Jansen et al., 2007).

3. Results

We now consider the performance of adiabatic and superadiabatic quantum computation, focusing on their time–energy complexity. This will be investigated both for the universal model of superadiabatic controlled gates and for the superadiabatic implementations of the Grover search.

3.1. Quantum Gates by Superadiabatic Controlled Evolutions

Let us begin by discussing the superadiabatic model of universal QC via CE implemented by shortcuts to adiabaticity (Santos and Sarandy, 2015). To this end, let us first write the complete set of eigenstates of H(t) as (Santos and Sarandy, 2015)

|E0mϵk(s)=|m,n^ϵ|E0k(s),
|E0(N1)+k(s)=|N1,n^+|E0k(s),
|Eϕ(N1)k(s)=|N1,n^|Eϕk(s),

where m ∈ {0, ⋯ , N − 2}, ϵ, k ∈ { ± }, and

|Eξ+(s)=sinθ0s20+eiξcosθ0s21,
|Eξ(s)=cosθ0s20+eiξsinθ0s21,

with ξ ∈ {0, ϕ} and {|Eξ±(s)} denoting the set of eigenstates of each adiabatic Hamiltonian Hξ(s), as provided by equation (2). Thus, by using the equation (22), we can show that the superadiabatic Hamiltonian is given by (Santos and Sarandy, 2015)

HSA(s)=1PN1,n^H0SA(s)+PN1,n^HϕSA(s),

where each term HξSA(s) corresponds to the superadiabatic Hamiltonian associated with the adiabatic Hamiltonian Hξ(s), i.e., HξSA(s)=Hξ(s)+HξCD, with

HξCD=θ02τσycosξσxsinξ

being the (time-independent) counter-diabatic contribution to achieve the evolution at total time τ (Santos and Sarandy, 2015).

3.2. Energy–Time Complementarity in the CE Model of Quantum Gates

Let us now consider equation (26) to investigate the time–energy complementarity relationship in both adiabatic and superadiabatic CE models of universal quantum gates. To this end, we need the set of eigenvalues and eigenstates of the adiabatic Hamiltonian in equation (1), which are given by Eqs 27–31. The spectrum of H(s) has (2N)-degenerate levels, with {|E0mϵ+(s),|E0(N1)++(s),|Eϕ(N1)+(s)⟩} and {|E0mϵ(s),|E0(N1)+(s),|Eϕ(N1)(s)⟩} associated with the levels E+ = ℏω and E = −ℏω, respectively. So, by using Eqs 27–31, we can show that μlm(s)=θ024τ2. In addition, the energetic cost to implement any gate controlled by n qubits is ΣSA(τ,n)=2n2ΣSAsing(τ) (Santos and Sarandy, 2015), where ΣSAsing is the energetic cost to implement any single qubit unitary transformation, with

ΣSAsingωτ,θ0=2ω1+θ024ωτ2.

A similar result can be obtained from the spectral norm, with energetic cost given by ΣSAsingωτ,θ0|2=(12)ΣSAsingωτ,θ0, since the Hilbert space has dimension D = 4 in this case. Note that the energetic cost is independent of the parameter θ0 in the adiabatic limit ωτ → ∞. Therefore, the best computational adiabatic strategy is to set θ = π, which deterministically ensures the implementation of the gate with probability one. On the other hand, probabilistic quantum computation can be energetically favored in the superadiabatic regime. Indeed, from equation (5), we can see that, by setting 0 < θ0 < π, the implementation of the quantum gate is achieved with a non-vanishing probability. Thus, we can investigate whether or not it is possible to find out a specific value of θ0 such that the energetic cost is better in average than the deterministic choice θ0 = π. To address this point, let us define the quantity

N=1sin2θ02,

which is the average number of evolutions for a successful computation. So, the average energetic cost to implement a probabilistic evolution is

Σ¯=NΣ,

where Σ is the cost of a single evolution. Without loss of generality, we will consider the cost of single gates, since similar arguments apply for the cost of n-qubit controlled gates. So, by performing superadiabatic probabilistic quantum computing, the average energetic cost is given by

Σ¯SAsingωτ,θ0=NΣSAsingωτ,θ0=2ωcsc2θ021+θ024ωτ2.

The function Σ¯SAsingωτ,θ0 as θ0 → 0 and exhibits a minimum in the interval 0 < θ0 < π as a function of ωτ. Indeed, the angle θ0min that minimizes Σ¯SAsingωτ,θ0 grows monotonically with ωτ, with θ0minπ as ωτ → ∞ (adiabatic limit). Then, optimizing Σ¯SAsingωτ,θ0 for θ0, we obtain

θ0Σ¯SAsingωτ,θ0=ηθ0,ωτθ04ωτ2+θ02cotanθ02=0,

where we have defined the function

ηθ0,ωτ=csc2θ022ωτ21+θ024ωτ2.

Note that η(θ0, ωτ) is non-vanishing in the whole interval Iθ00,π. Thus, to obtain the critical angle θ0min in Σ¯SAsingωτ,θ0, we use equation (38) to note that ωτ satisfies

ωτ=θ0min2tanθ0min2θ0min,

where we can see a dependence of θ0min on the choice of ωτ. In addition, note that θ0min is such that tanθ0min2  θ0min, since the quantity ωτ is required to be real and positive. The probabilistic advantage is plotted in Figure 1, where it is shown that the optimal value for θ0 is a continuous function of ωτ, being distinct of the deterministic implementation θ0 = π. In the inset, we show the global minimum of the average energy for ωτ = 0.01, which occurs for θ0 < π. In particular, θ0min moves away from π as ωτ is lowered, i.e., in the strong superadiabatic regime. As ωτ shifts toward the adiabatic limit, we find that θ0minπ. The optimization of the energy cost is shown in the lower inset, where we define the fraction of energy required by the optimized probabilistic model as a function of ωτ as

Σrelωτ=Σ¯SAsing(ωτ,θ0min)ΣSAsingωτ,π.
FIGURE 1
www.frontiersin.org

Figure 1. Optimal value θ0min for the angle parameter θ0 as a function of ωτ, with ωτ in logarithmic scale. The points are obtained from equation (40), with the curve denoting the numerical fit. Upper inset: average energy in units of ℏω as a function of θ0 for ωτ = 0.01. The results are obtained from equation (37). Lower inset: fraction Σrel(ωτ) of energy required by the optimized probabilistic model as a function of ωτ, with data in logarithmic scale. The points are obtained from equation (41), with the curve denoting the numerical fit.

Notice that Σrel(ωτ) decreases in the superadiabatic regime, implying into a large reduction of the energetic cost for small values of ωτ. On the other hand, Σrel(ωτ) →1 in the adiabatic limit, since θ0minπ.

3.3. Superadiabatic Quantum Search

Here, we derive a superadiabatic Hamiltonian HSA(s) for the oracular quantum search governed by the adiabatic Hamiltonian H0(s) in equation (7). We will adopt linear interpolation, with f (s) = 1 − s and g(s) = s, as in Farhi et al. (2000) and write HSA(s) = H0(s) + HCD(s). In order to determine the counter-diabatic Hamiltonian HCD(s), we observe that, since H0(s) has real eigenstates, we use that n˙(s)|n(s)=0 in equation (23), which implies that

HCD(s)=iℏτξ=±|E˙ξ(s)Eξ(s)|,

where the energies |E±(s)⟩ are given by equation (8) and

|E˙±(s)=(N1)b±b˙±(1+(N1)b±2)32|m+b˙±(1+(N1)b±2)32|ϕ.

Note that only the ground and first excited states contribute to HCD(s), since the higher-energy degenerate sector {|Edegk} is composed by time-independent eigenvectors [see equation (12)]. Note also that the counter-diabatic Hamiltonian will naturally be non-oracular [see equation (18)], with contributions from operators, such as |ϕ⟩ ⟨m| and |m⟩ ⟨ϕ|. This is the reason behind the time complexity O(1) for the superadiabatic Hamiltonian. Naturally, such a result leads to an artificial approach. In a more physical scenario, superadiabaticity could be applied to the quantum search via the direct implementation of the Grover quantum circuit, through the controlled evolution approach discussed in Section 3.1.

3.4. Energy–Time Complementarity in the Quantum Search

Let us now analyze the time–energy complementarity relationship in the adiabatic and superadiabatic versions of the Grover search. In the adiabatic regime, the energetic cost can be computed from equation (26) and using τ → ∞. Therefore, the adiabatic cost can be written as

Σad=01dsmEm2(s)=01dsE+(s)2+E(s)2+(N2)Edeg(s)2.

Let us initially consider the oracular Hamiltonian H0(s) in equation (7), whose eigenvalues are given by equations (10) and (11). By considering the case of local adiabatic evolution provided by the interpolation in equation (13) and by taking N ≫ 1, we obtain E±(s) ∼ Edeg(s) ∼ O(1), which implies from equation (44) into an energetic cost ΣadLA that scales as O(N). On the other hand, in the superenergetic version of the quantum search, we adopt the interpolation in Eqs 14 and 15. Then, by taking N ≫ 1, we obtain now E±(s)Edeg(s)O(N), which implies into an energetic cost ΣadSE that scales as O(N). This higher energetic cost is a consequence of the complementarity between energy and time, which arises to compensate the constant time complexity O(1) of the superenergetic version. Naturally, the composite energy–time complexity is kept constant for both cases. This overall complexity is reduced by taking non-oracular artificial Hamiltonians. In the case of the adiabatic NLNO model, we use the Hamiltonian in equation (17), whose ground state and first excited state energies are given now by equation (20), with the higher energies kept as in equation (11). Its energetic cost ΣadNO can also be computed from equation (44) by considering the interpolation in equation (21) and by taking N ≫ 1. Then, we obtain E±(s) ∼ Edeg(s) ∼ O(1), which yields ΣadNO scaling as O(N).

For the superadiabatic algorithm, equation (26) must be used. Without loss of generality, we set energy units such that /τ = 1. We find that for N ≫ 1, the value of μ±(s) in equation (26) evaluate to

μ±(s)=E˙±(s)|E˙±(s)=(N1)b˙±(s)2(1+(N1)b±(s)2)2,

which in turn gives the superadiabatic search energetic cost ΣSA of order O(N), just reproducing the scaling of the NLNO adiabatic search. Similar results can be obtained if one chooses the spectral norm in the energetic cost, up to a common scaling factor D12=N related to the dimension of the Hilbert space. These results are summarized in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Energy–time complexity for several versions of oracular and non-oracular Hamiltonians for the Grover quantum search.

4. Discussion

We have discussed the energetic cost of shortcuts to adiabaticity and their consequences in quantum information processing. Specifically, we considered both the superadiabatic universal gate model via CE and the superadiabatic analog quantum search. For the gate model, we have shown that, differently from the adiabatic scenario, superadiabatic probabilistic gate implementations are energetically favorable with respect to deterministic gate implementations. This implies that the additional energy resources required by superadiabatic evolutions can be minimized by a suitable probabilistic model. Indeed, probabilistic evolutions have recently been considered in similar applications for QC. In particular, they have been used to cancel errors in adiabatic processes (Kieferová and Wiebe, 2014) and as a technique to decompose unitary operations (Paetznick and Svore, 2014; Bocharov et al., 2015). Here, we have shown a new aspect of probabilistic QC, which corresponds to an advantage in the energy balance for superadiabatic dynamics while keeping its performance for a fixed evolution time. For analog quantum search, we have shown that the superadiabatic approach induces a non-oracular counter-diabatic Hamiltonian, with energy–time complexity equivalent to non-oracular adiabatic implementations. This explicitly shows that the Grover optimality is robust against transitionless drivings, which is reflected by a fixed energy–time scaling of the Hamiltonian.

Implications of probabilistic superadiabatic QC under decoherence is a further challenge of immediate interest. In a quantum open-systems scenario, there is a compromise between the time required by adiabaticity and the decoherence time of the quantum device. Therefore, a superadiabatic implementation may provide a direction to obtain an optimal running time for the quantum algorithm while keeping an inherent protection against decoherence. In this context, it is our interest to understand to what extent decoherence can affect the optimal angle θ0min, investigating in particular if it can be robust against classes of decohering processes. Concerning specifically the Grover search, it would be interesting to understand whether superadiabatic implementations are equivalent to arbitrary non-oracular adiabatic Hamiltonians, as suggested in our present discussion. Moreover, the behavior of correlations, such as entanglement and the investigation of experimental proposals in the superadiabatic scenario are also topics under investigation.

Author Contributions

All the authors equally contributed to the conception of the work, development of the main results, and writing of the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

MS thanks Daniel Lidar for his hospitality at the University of Southern California. IC, AS, and MS acknowledge support from CNPq/Brazil and the Brazilian National Institute for Science and Technology of Quantum Information (INCT-IQ).

References

Barends, R., Shabani, A., Lamata, L., Kelly, J., Mezzacapo, A., Las Heras, U., et al. (2016). Digitized adiabatic quantum computing with a superconducting circuit. Nature 534, 222–226. doi: 10.1038/nature17658

CrossRef Full Text | Google Scholar

Berry, M. V. (2009). Transitionless quantum driving. J. Phys. A: Math. Theor. 42, 365303. doi:10.1088/1751-8113/42/36/365303

CrossRef Full Text | Google Scholar

Bocharov, A., Roetteler, M., and Svore, K. M. (2015). Efficient synthesis of universal repeat-until-success quantum circuits. Phys. Rev. Lett. 114, 080502. doi:10.1103/PhysRevLett.114.080502

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y. H., Xia, Y., Chen, Q. Q., and Song, J. (2014a). Efficient shortcuts to adiabatic passage for fast population transfer in multiparticle systems. Phys. Rev. A 89, 033856. doi:10.1103/PhysRevA.89.033856

CrossRef Full Text | Google Scholar

Chen, Y.-H., Xia, Y., Chen, Q. Q., and Song, J. (2014b). Shortcuts to adiabatic passage for multiparticles in distant cavities: applications to fast and noise-resistant quantum population transfer, entangled states preparation and transition. Laser Phys. Lett. 11, 115201. doi:10.1088/1612-2011/11/11/115201

CrossRef Full Text | Google Scholar

Chen, Y.-H., Xia, Y., Chen, Q. Q., and Song, J. (2015a). Fast and noise-resistant implementation of quantum phase gates and creation of quantum entangled states. Phys. Rev. A 91, 012325. doi:10.1103/PhysRevA.91.012325

CrossRef Full Text | Google Scholar

Chen, Y.-H., Xia, Y., Song, J., and Chen, Q.-Q. (2015b). Shortcuts to adiabatic passage for fast generation of Greenberger-Horne-Zeilinger states by transitionless quantum driving. Sci. Rep. 5, 15616. doi:10.1038/srep15616

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, S., Kobes, R., and Kunstatter, G. (2003). Energy and efficiency of adiabatic quantum search algorithms. J. Phys. A: Math. Gen. 36, 2839. doi:10.1088/0305-4470/36/11/313

CrossRef Full Text | Google Scholar

Deffner, S., and Lutz, E. (2013). Energy time uncertainty relation for driven quantum systems. J. Phys. A: Math. Theor. 46, 335302. doi:10.1088/1751-8113/46/33/335302

CrossRef Full Text | Google Scholar

del Campo, A. (2013). Shortcuts to adiabaticity by counter-diabatic driving. Phys. Rev. Lett. 111, 100502. doi:10.1103/PhysRevLett.111.100502

CrossRef Full Text | Google Scholar

del Campo, A., Rams, M. M., and Zurek, W. H. (2012). Assisted finite-rate adiabatic passage across a quantum critical point: exact solution for the quantum Ising model. Phys. Rev. Lett. 109, 115703. doi:10.1103/PhysRevLett.109.115703

PubMed Abstract | CrossRef Full Text | Google Scholar

Demirplak, M., and Rice, S. A. (2003). Adiabatic population transfer with control fields. J. Phys. Chem. A 107, 9937. doi:10.1021/jp030708a

CrossRef Full Text | Google Scholar

Demirplak, M., and Rice, S. A. (2005). Assisted adiabatic passage revisited. J. Phys. Chem. B 109, 6838. doi:10.1021/jp040647w

PubMed Abstract | CrossRef Full Text | Google Scholar

Farhi, E., Goldstone, J., Gutmann, S., Lapan, J., Lundgren, A., and Preda, D. (2001). A quantum adiabatic evolution algorithm applied to random instances of an NP-complete problem. Science 292, 472. doi:10.1126/science.1057726

PubMed Abstract | CrossRef Full Text | Google Scholar

Farhi, E., Goldstone, J., Gutmann, S., and Sipser, M. (2000). Quantum computation by adiabatic evolution. MIT-CTP-2936 (arXiv:quant-ph/0001106).

Google Scholar

Grover, L. K. (1997). Quantum mechanics helps in searching for a needle in a haystack. Phys. Rev. Lett. 79, 325. doi:10.1103/PhysRevLett.79.325

CrossRef Full Text | Google Scholar

Hen, I. (2015). Quantum gates with controlled adiabatic evolutions. Phys. Rev. A 91, 022309. doi:10.1103/PhysRevA.91.022309

CrossRef Full Text | Google Scholar

Jansen, S., Ruskai, M.-B., and Seiler, R. (2007). Bounds for the adiabatic approximation with applications to quantum computation. J. Math. Phys. 48, 102111. doi:10.1063/1.2798382

CrossRef Full Text | Google Scholar

Kieferová, M., and Wiebe, N. (2014). On the power of coherently controlled quantum adiabatic evolutions. New J. Phys. 16, 123034. doi:10.1088/1367-2630/16/12/123034

CrossRef Full Text | Google Scholar

Lu, M., Xia, Y., Shen, L. T., and Song, J. (2014a). An effective shortcut to adiabatic passage for fast quantum state transfer in a cavity quantum electronic dynamics system. Laser Phys. 24, 105201. doi:10.1088/1054-660X/24/10/105201

CrossRef Full Text | Google Scholar

Lu, M., Xia, Y., Shen, L. T., Song, J., and An, N. B. (2014b). Shortcuts to adiabatic passage for population transfer and maximum entanglement creation between two atoms in a cavity. Phys. Rev. A 89, 012326. doi:10.1103/PhysRevA.89.012326

CrossRef Full Text | Google Scholar

Martinis, J. M., and Geller, M. R. (2014). Fast adiabatic qubit gates using only σz control. Phys. Rev. A 90, 022307. doi:10.1103/PhysRevA.90.022307

CrossRef Full Text | Google Scholar

Messiah, A. (1962). Quantum Mechanics. North-Holland, Amsterdam: North-Holland Pub. Co.

Google Scholar

Nielsen, M. A., and Chuang, I. L. (2000). Quantum Computation and Quantum Information. Cambridge: Cambridge University Press.

Google Scholar

Orus, R., and Latorre, J. I. (2004). Universality of entanglement and quantum computation complexity. Phys. Rev. A 69, 052308. doi:10.1103/PhysRevA.69.052308

CrossRef Full Text | Google Scholar

Paetznick, A., and Svore, K. M. (2014). Repeat-until-success: non-deterministic decomposition of single-qubit unitaries. Quantum Inf. Comput. 14, 1277.

Google Scholar

Roland, J., and Cerf, N. J. (2002). Quantum search by local adiabatic evolution. Phys. Rev. A 65, 042308. doi:10.1103/PhysRevA.65.042308

CrossRef Full Text | Google Scholar

Saberi, H., Opatrny, T., Molmer, K., and del Campo, A. (2014). Adiabatic tracking of quantum many-body dynamics. Phys. Rev. A 90, 060301(R). doi:10.1103/PhysRevA.90.060301

CrossRef Full Text | Google Scholar

Santos, A. C., and Sarandy, M. S. (2015). Superadiabatic controlled evolutions and universal quantum computation. Sci. Rep. 5, 15775. doi:10.1038/srep15775

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, A. C., Silva, R. D., and Sarandy, M. S. (2016). Shortcut to adiabatic gate teleportation. Phys. Rev. A 93, 012311. doi:10.1103/PhysRevA.93.012311

CrossRef Full Text | Google Scholar

Sarandy, M. S., Wu, L.-A., and Lidar, D. (2004). Consistency of the adiabatic theorem. Quantum Inf. Process. 3, 331. doi:10.1007/s11128-004-7712-7

CrossRef Full Text | Google Scholar

Teufel, S. (2003). “Adiabatic perturbation theory in quantum dynamics,” in Lecture Notes in Mathematics, Vol. 1821 (Berlin Heidelberg: Springer-Verlag).

Google Scholar

Torrontegui, E., Ibáñez, S., Martínez-Garaot, S., Modugno, M., del Campo, A., Guéry-Odelin, D., et al. (2013). Shortcuts to adiabaticity. Adv. Atom. Mol. Opt. Phys. 62, 117. doi:10.1016/B978-0-12-408090-4.00002-5

CrossRef Full Text | Google Scholar

van Dam, W., Mosca, M., and Vazirani, U. (2001). “How powerful is adiabatic quantum computation?,” in Proceedings of the 42nd Annual Symposium on Foundations of Computer Science, Las Vegas, 279.

Google Scholar

Wen, J., Huang, Y., and Qiu, D. (2009). Entanglement proprieties of adiabatic quantum algorithms. Int. J. Quantum Inform. 7, 1531. doi:10.1142/S0219749909006000

CrossRef Full Text | Google Scholar

Wen, J., and Qiu, D. (2008). Entanglement in adiabatic quantum searching algorithms. Int. J. Quantum Inform. 6, 997. doi:10.1142/S0219749908004249

CrossRef Full Text | Google Scholar

Zheng, Y., Campbell, S., Chiara, G., and Poletti, D. (2015). Cost of transitionless driving and work output. e-print arXiv:1509.01882

Google Scholar

Keywords: quantum computing, quantum information, shortcuts to adiabaticity, superadiabaticity, quantum gates, quantum search

Citation: Coulamy IB, Santos AC, Hen I and Sarandy MS (2016) Energetic Cost of Superadiabatic Quantum Computation. Front. ICT 3:19. doi: 10.3389/fict.2016.00019

Received: 29 March 2016; Accepted: 23 August 2016;
Published: 15 September 2016

Edited by:

Adam Miranowicz, Adam Mickiewicz University in Poznań, Poland

Reviewed by:

Zhao-Ming Wang, Ocean University of China, China
Daniel Matthias Herr, RIKEN, Japan

Copyright: © 2016 Coulamy, Santos, Hen and Sarandy. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Marcelo S. Sarandy, msarandy@id.uff.br