Decoherence induced deformation of the ground state in adiabatic quantum computation
Despite more than a decade of research on adiabatic quantum computation (AQC), its decoherence properties are still poorly understood. Many theoretical works have suggested that AQC is more robust against decoherence, but a quantitative relation between its performance and the qubits’ coherence properties, such as decoherence time, is still lacking. While the thermal excitations are known to be important sources of errors, they are predominantly dependent on temperature but rather insensitive to the qubits’ coherence. Less understood is the role of virtual excitations, which can also reduce the ground state probability even at zero temperature. Here, we introduce normalized ground state fidelity as a measure of the decoherence-induced deformation of the ground state due to virtual transitions. We calculate the normalized fidelity perturbatively at finite temperatures and discuss its relation to the qubits’ relaxation and dephasing times, as well as its projected scaling properties.
Adiabatic quantum computation (AQC) farhi (); gs (), either in its universal form ahar04 (); miz06 (), or in the form of adiabatic quantum optimization exp1 (); exp2 (), or quantum simulations sim (), presents a viable alternative to gate-model quantum computation (GMQC). Although a part of the original motivation for introduction of the AQC gs () was the promise of the increased stability against decoherence due to the energy gap between the ground and excited states, the question of the role of decoherence in AQC remains an open one. This uncertainty makes it important to quantify more precisely the decoherence properties of AQC. A crucial step towards this would be to define a quantitative characteristic of the decoherence strength in AQC, that plays a role similar to the decoherence time for GMQC. However, in the case of AQC, decoherence has qualitatively different, static effect on the qubits, not limiting the operation time of an algorithm Amin09b ().
In AQC, adiabatic evolution of the ground state of a qubit system realizes the solution of a computational problem represented by an appropriately designed Hamiltonian, which is typically written as
where with being the total evolution time. At , one has , , and the system is initialized in the ground state of the initial (driver) Hamiltonian , which usually consists of the uniform superposition of all computational basis states. The energy scales and are then varied monotonically so that at , and . If the evolution is slow enough, an isolated qubit system stays in the ground state with high fidelity throughout the evolution, and at reaches the ground state of the final (problem) Hamiltonian , which provides a solution to a computational problem.
If the qubit system is weakly coupled to a dissipative environment, two effects are expected. First, the low-frequency part of the environmental noise moves the system energy levels relative to each other. This results in a dephasing of the energy eigenstates that eventually suppresses all off-diagonal elements of the qubit density matrix in the energy basis. However, since the population of the ground state is the only important part of the computation and the relative phases of the energy eigenstates do not carry any information, this does not affect AQC. The second effect of the coupling to environment is that it induces thermal transitions between the qubit energy levels pushing the qubit system towards thermal equilibrium at a temperature . In the limit of large , the instantaneous probability to be in the ground state asymptotically approaches the Boltzmann distribution, and so the qubit system loses some of the ground state probability due to thermal occupation of the excited states. Such a thermal loss of probability can be compensated by multiple iterations of an AQC algorithm as long as it does not scale exponentially with the size of the system, i.e., as long as the number of excited states within roughly the energy above the ground state does not grow exponentially.
The preceding arguments provide an intuitive explanation for the predicted robustness of AQC against local environmental noise in the limit of weak coupling Lloyd (); Amin09b (); Amin09 (); Amin08 (); Tiersch (); Sarandy (); Roland (); Childs (). When the strength of the coupling to the environment is increased without changing either the Hamiltonian or the temperature, the qubit system’s Boltzmann distribution is still not directly affected. However, it is known that the decoherence time of the qubits decreases with increased coupling, and strong coupling to the environment eventually makes the qubits completely incoherent, rendering them useless for quantum computation. In GMQC, qubit decoherence leads to computation errors which, without error correction, completely destroy the computation process. This is why the qubits’ quality factor, which is the ratio of the decoherence time and the gate operation time, provides a good measure of the qubit performance in GMQC. It is, however, unclear how an increase in coupling to the environment, or equivalently decrease in qubit quality factor, affects AQC.
In this paper, we look closely at what happens to the eigenstates of the qubit system in AQC when coupling to the environment is non-negligible but still small enough to allow perturbation expansion. We introduce the normalized ground state fidelity, defined as the distance between the open and closed system reduced density matrices normalized to the Boltzmann ground state probability, as a quantitative measure of decoherence-induced deformation of the ground state in AQC, analogous to the decoherence time for GMQC. We calculate the fidelity perturbatively and express it through the same environmental noise correlators that determine the decoherence times in GMQC. Such an equilibrium calculation of the normalized ground state fidelity is accurate in the long limit, but becomes approximate when the evolution is too fast for the system to reach local equilibrium. However, the deviation from the equilibrium distribution is largest when the rate of relaxation between the eigenstates becomes extremely small (e.g., near the end of the evolution in adiabatic quantum optimization). As we shall see, the normalized ground state fidelity is closely related to the relaxation rate and becomes close to 1, independent of the detailed probability distribution, when the relaxation is very slow. Therefore, our calculation of the normalized fidelity can provide a good approximation in all regions as long as the evolution time is not too short.
i.1 Normalized ground state fidelity - definition
We first provide a definition for the normalized ground state fidelity based on the notion of fidelity between two density matrices. To ensure consistent notation throughout this paper, symbols with (without) “” denote quantities related to the coupled (uncoupled) qubit system and environment. We use letters to enumerate the eigenstates and eigenvalues of the qubit system (e.g., , ), letters to enumerate the eigenstates and eigenvalues of the environmental degrees of freedom, and letters to enumerate the eigenstates and eigenvalues of the total system (qubits+environment). The total Hamiltonian is , where and are the environment and interaction Hamiltonians, respectively. In the absence of coupling, , and the eigenstates of the total system are with eigenvalues . When , the new eigenstates are , which typically are entangled superpositions of the unperturbed states . For weak coupling, is very close to and the effect of the environment is thermalization of the qubit system. Once the environment is averaged out, the equilibrium of the total system gives the Boltzmann distribution for the qubits:
where and are the partition functions of the qubit system and the environment.
As the coupling strength increases, the deviation of from grows. In equilibrium, the density matrix of the total system still has the Boltzmann form , where , with being the partition function of the total system. However, the reduced density matrix of the qubit system alone is no longer given by the Boltzmann distribution. The deviation from the Boltzmann form provides a good qualitative measure of how strongly the eigenstates are deformed in comparison to the unperturbed states.
To quantify the loss of fidelity of the ground state due to such deformation of the energy eigenstates of the uncoupled system, it is convenient first to separate this effect from the loss of fidelity due to thermal excitations. This can be done by using normalized ground state fidelity, which we define as the Uhlmann fidelity Uhlmann () between the reduced density matrix and the “ideal” ground state density matrix , normalized to the Boltzmann ground state probability :
where is the equilibrium probability for the qubits to be in the ideal ground state when coupled to the environment. Normalization to the equilibrium Boltzmann ground state probability is natural in the context of this work, since in the calculations presented below we adopt the assumption that the qubit-environment system maintains equilibrium throughout the AQC evolution. In the weak-coupling limit, no deformation of the eigenstates is expected. Then , and Eq. (3) gives . This shows that Eq. (3) indeed correctly separates the effect of the quantum deformation of the ground state, which can be viewed as the result of virtual transitions to the excited states, from the thermal loss of probability. Qualitatively, the effect of the virtual transitions, expressed in , is different from that of the thermal transitions in two important aspects. First, it persists even at , when all the thermal transitions are suppressed. Second, it depends on the strength of coupling to the environment (or decoherence time of the qubits), while thermal equilibrium probabilities only depend on the energy eigenvalues and temperature. Nevertheless, similarly to the thermal transitions, the virtual transitions reduce the occupation probability of the ground state by transferring it to other low energy states. In general, beyond the AQC, deviation of the occupation probabilities of an equilibrium quantum system from the Boltzmann distribution due to non-vanishing strength of coupling to environment has been studied before as introducing corrections to classical thermodynamics (see, e.g., AN (); Nag02 (); Scully03 (); Nori10 (); WLJ () and references therein). In particular, deformation of the ground state by coupling to an environment, described in our case as a suppression of fidelity, is known to lead to several physical effects, e.g., suppression of the persistent current in normal-metal rings PC () or violation of the fluctuation-dissipation theorem for the thermal conductance FDT ().
i.2 Perturbative calculation of fidelity
We calculate the fidelity (3) perturbatively and relate it to measurable parameters of the qubit system and environment. As appropriate for AQC, we assume that the coupling is weak. This allows us to employ the perturbation theory in around the non-interacting state of the qubit system and the environment. To separate the effects of environment from other deviations from the perfect adiabatic evolution of an AQC algorithm, we consider the limit of slow evolution, when the rate of change of the control parameter of the Hamiltonian is small, e.g., smaller than the energy relaxation rate of the qubit system. In this case, non-adiabatic transitions out of the ground state can be neglected, and the occupation probabilities of the excited states (if non-vanishing) correspond to local equilibrium at each moment during the evolution. While this is obviously not the most general case of time evolution in AQC, which can be dominated by non-equilibrium effects, it covers the most important regime of ideal AQC evolution, and is appropriate for the situation of sufficiently slow evolution in the presence of small but finite interaction strength with environment we want to describe. The zeroth-order state of the perturbation theory in this regime is the state with both the uncoupled qubit system and environment in equilibrium at the same temperature , i.e., the total density matrix of the system being the product of the density matrices and , where and are the Boltzmann probabilities.
In general, the reduction of the ground state probability due to finite is caused by two effects. First is the change in the equilibrium probability () as a result of renormalization of the energy eigenvalues (Lamb shifts) of the qubit system. Second effect, conceptually more important for this work, is the probability transfer into and out of the ground state due to renormalization of the qubit system wavefunctions. Explicitly, the probability that defines the normalized fidelity (3) can be expressed as
Introducing interaction-induced corrections to the equilibrium probabilities , where , and wavefunctions: , we can rewrite this expression to the lowest non-vanishing order in as
Using the relation to transform the term in Eq. (5) we obtain
The terms proportional to in Eq. (6) describe the reduction of the ground state probability as a result of virtual transitions between the ground and excited state due to the interaction with the environment.
Next, we calculate and . Quite generally, the interaction Hamiltonian is
where are the Pauli matrices for the th qubit, , and are the corresponding operators of the noise generated by the environment. As usual, the averages of the noise operators vanish, . Then, in the weak coupling regime, the effect of environment is fully characterized by the noise spectral densities:
where is the average over the environmental degrees of freedom. For simplicity, we limit our discussion to the most typical case when the noises with different or are uncorrelated. It is shown in the supplementary information (SI) that the perturbation expansion in in this situation gives
Equation (10) is our main result. The normalized fidelity is well-defined at , when all thermal excitations are suppressed, i.e., for and at . Hence, the values of around , when the denominator in (10) vanishes, do not contribute to the integral. When , the divergences that appear at reflect the fact that environment can also create real thermal excitations of the qubit system. However, the detailed balance relation, , ensures that these divergences cancel each other out and Eq. (10) is well-defined also at (see the SI).
i.3 Normalized fidelity for single qubit
Equation (10) is now applied to specific problems. The first example we consider is a typical individual qubit with the Hamiltonian
coupled as in Eq. (7), but only through , to the environmental noise with spectral density (8). In the usual weak-coupling approximation [see, e.g., Ref. blum, ], the qubit decoherence time is given by
where and are the relaxation and pure dephasing times, given by
We see that the same noise spectral density that defines the relaxation and dephasing rates (13) and (14) of the qubits in the GMQC determines the reduction of the ground-state normalized fidelity in AQC. In this respect, the main difference between the reduction of normalized fidelity and the real-time relaxation and dephasing is that even in the lowest-order perturbation theory, the normalized fidelity is reduced by the whole spectrum of environmental excitations, and not just by limited spectral groups resonant with the qubit energy differences or the low-frequency excitations, as in Eqs. (13) and (14).
To strengthen this comparison, we consider an Ohmic bath characterized by the noise
where, is a dimensionless coefficient and is the cutoff frequency. In this case, the relaxation time is and the normalized fidelity is expressed as
where is the qubit quality factor due to relaxation. Equation (15) gives the following expressions for the factor at low and high temperatures:
Equation (17) relates the normalized ground state fidelity to the qubit quality factor, , as calculated due to relaxation only. This shows that the normalized fidelity can be related more closely to the relaxation () and not dephasing () processes. Adding a low-frequency noise of a realistic magnitude does not change this conclusion (as discussed in more details in the numerical examples below). As expected, a larger leads to a better . Figure 1 shows the factor in Eq. (17) as a function of temperature for different cut-off frequencies . It exhibits the non-monotonic -dependence, and only weak, logarithmic, dependence on , which allows one to estimate without precisely specifying . The factor is maximal around at , which leads to a minimum normalized fidelity . Notice that even a qubit quality factor as low as , which is practically useless for GMQC, leads to %.
i.4 Normalized fidelity for multi-qubit systems
We now consider multi-qubit systems, starting with a system of uncoupled qubits. In this case, the trace in the definition of normalized fidelity (3) can be taken independently over separate qubits, so that the total is the product of fidelities , of the individual qubits: . For instance, a typical starting point of AQC algorithms is to initialize the system in the ground state of the Hamiltonian (20). Then, the state of all qubits is the same and can be characterized by the same normalized fidelity (17). Then,
For independent qubits, scales exponentially with as a result of the exponential scaling of the probability for all qubits to remain in their corresponding ground states. Since is inversely proportional to the noise strength , by decreasing the noise by a factors of, e.g., 10, one can achieve the same with 10 times more qubits.
Next, we focus on how the normalized ground state fidelity behaves in practical AQC systems. We use as an example the D-Wave One quantum annealing processor, as the one installed at the University of Southern California (see Ref. exp2, ). The Hamiltonian implemented by the processor has the form of Eq. (1), with
where and are tunable dimensionless bias and coupling coefficients. The parameters and for this processor are plotted in Fig. 2b. We calculate for a ferromagnetic chain (illustrated in Fig. 2a) with and , otherwise known as a quantum Ising model in a transverse field. Here, the length of the chain is varied from to 16. Although this model is exactly solvable (see, e.g., Ref. is, and references therein), cannot be calculated exactly for practical noise models in which the coupling to environment is dominated by the terms. Hence, we calculate the normalized fidelity numerically. In the limit , the model is known to have a quantum critical point at where . At this point, the chain goes through a quantum phase transition between quantum paramagnetic and ferromagnetic phases. In the ferromagnetic phase, the ground state is doubly degenerate with respect to simultaneous change of signs of all terms. Figure 2c plots several of the lowest energy levels of a 10-qubit chain relative to the ground state energy . In the thermodynamic limit (), the appearance of the doubly-degenerate ground state and the minimum in the energy gap between the ground and the second excited states happen at the quantum critical point. For the 10-qubit chain of Fig. 2, however, these happen at slightly different points than the one defined by .
To calculate for this system, we use a realistic noise model relevant to the D-Wave qubits Harris10 (). In this case, the dominant environmental coupling is to the magnetic flux noise, which couples directly to the qubit computational basis states represented by the operators. The noise spectral density was characterized in the earlier experiments, which were consistent with the noise being a combination of the low-frequency noise and an Ohmic noise at high frequencies env2 (). For calculations of , we take , where is the Ohmic spectral density (16) and . The coefficient appears because the strength of coupling to flux noise depends on the persistent current of the flux qubits which changes as a function of (see SI). Here, is the bias point at which the measurements of and are performed. Based on the experimental data, we use , MHz and . We also assume GHz for the high-frequency cutoff and MHz for the low-frequency cutoff (based on a s evolution time of an algorithm). We found that for these parameters, is dominantly determined by the high-frequency Ohmic noise and not by the noise.
In principle, since the total number of energy levels grows exponentially with , the time required for numerical calculation of also grows exponentially. Fortunately, the value of converges rapidly for a finite number of retained energy states. Here, we keep all energy levels for , and up to 2000 energy levels for larger chains. The normalized ground state fidelity of the 10-qubit chain is plotted as a function of in Fig. 2d. The fidelities of chains with other lengths (and also coupled systems other than chains) are qualitatively the same as the one plotted in Fig. 2d. It is clear from the figure that is minimum close to the critical point . Notice also that the fidelity approaches 1 as , which is the result of commuting with , with only terms and negligible other types of coupling to environment. This again reflects the fact that depends rather on relaxation than dephasing.
Figure 3 shows the numerical results for the normalized ground state fidelity for -qubit chains with to 16 at the critical point. For all chain lengths, is better than 90%. It should be emphasized that these are the minimum fidelities at the quantum critical point . The normalized fidelity at all other points is larger, and near , is very close to 1 as shown in Fig. 2d. We have also plotted in Fig. 3 the normalized ground state fidelity of uncoupled qubits (with the same parameters) at different based on the exponential scaling of Eq. (19). The scaling and magnitude of at large is better for the ferromagnetic chain than for the uncoupled qubits. A plausible reason for this is that the spin-spin interaction introduces additional rigidity into the chain dynamics making it less susceptible to the environmental perturbations, and therefore increasing the fidelity. Unfortunately, it was not possible to pursue numerical calculations beyond 16 qubits, as direct perturbation approximation would break down when strongly deviates from unity. A naive exponential extrapolation of the data points to (representing the worse case) still yields , meaning that the eigenstates could retain their quantum properties without error correction for such a large-size system. As in uncoupled qubits, if one can reduce the noise by a large factor, the size of the chain can be increased by the same factor while keeping unchanged. In addition, other techniques such as dynamical decoupling DDecoupling () or error correction ECorrect () could be employed to enhance the normalized ground state fidelity at large scales.
Finally, we discuss how the normalized ground state fidelity should affect the performance of AQC. First, we notice that the actual equilibrium ground state probability at point is , where is the Boltzmann probability. Therefore, a suppression of the normalized fidelity creates an extra reduction of the ground state probability on top of the thermal reduction. In universal AQC ahar04 (); miz06 (), directly affects the quality of the computation. Indeed, deviations of from 1 mean that the statistics of measurements done on the final state will be different from the one that corresponds to the ideal ground state. For instance, in the case of one qubit with and the Hamiltonian (11), measurement of has a non-vanishing probability of producing the result different from the ground state even at temperatures . However, this effect is absent in the special case when the coupling to environment via commutes with the final Hamiltonian , leading to at the end of evolution, as in the adiabatic quantum optimization discussed above and shown in Fig. 2d. In this case, a small fidelity in the middle of the evolution increases the loss of probability due to the thermal transitions, thereby decreasing the ground state probability even further. Therefore, the probability will be distributed among the low energy states even more than implied by the thermal equilibrium. Part of the probability can be regained later when the gap is larger and is closer to 1. However, since the relaxation time becomes exponentially long near the end of evolution, the majority of the probability that is lost may not be gained back, thus leading to a smaller probability of success. This makes it important to maintain close to unity throughout the evolution. We stress that most treatments of AQC based on the weak coupling master equation, e.g., Childs (); Amin09b (); master_eq (), even with Lamb shift, do not take into account the effect of deformation of the eigenstates that is captured by our calculation of the normalized fidelity.
In summary, we have proposed using normalized ground state fidelity as a quantity for measuring the strength of decoherence effects in AQC. The fidelity plays a role similar to decoherence time in GMQC, but takes into account qualitatively different effects of environment on the ground state relevant to AQC. The fidelity is related to the relaxation processes and is relatively insensitive to the dephasing. Our numerical calculations indicate that a normalized fidelity close to unity can be achieved with a moderate qubit quality factor, even for large numbers of qubits. Normalized ground state fidelity should be a useful measure of the environment related quality of AQC systems in the context of further work on important topics in AQC such as quantum error correction or the threshold theorem.
MHA is grateful to A.J. Berkley, M.W. Johnson, R. Liu, T. Mahon, F. Nori, A. Yu. Smirnov, and B. Wilson for discussions and comments.
- (1) Farhi, E., et al. A Quantum adiabatic evolution algorithm applied to random instances of an NP-complete problem, Science 292, 472 (2001).
- (2) Mizel, A., Mitchell, M.W., and Cohen, M.L. Energy barrier to decoherence, Phys. Rev. A 63, 040302 (2001).
- (3) Aharonov, D., et al. Adiabatic quantum computation is equivalent to standard quantum computation, SIAM J. Comput. 37, 166 (2007).
- (4) Mizel, A., Lidar, D.A., and Mitchell, M. Simple proof of equivalence between adiabatic quantum computation and the circuit model, Phys. Rev. Lett. 99, 070502 (2007).
- (5) Johnson, M.W. et al. Quantum annealing with manufactured spins, Nature 473, 194 (2011).
- (6) Bian, Z., Chudak, F., Macready, W.G., Clark, L., Gaitan F. Experimental determination of Ramsey numbers with quantum annealing, eprint arXiv:1201.1842.
- (7) Biamonte, J.D., Bergholm, V., Whitfield, J.D., Fitzsimons, J., and Aspuru-Guzik, A. Adiabatic quantum simulators, AIP Advances 1, 022126 (2011).
- (8) Amin, M.H.S., Truncik, C.J.S., and Averin, D.V. Role of single qubit decoherence time in adiabatic quantum computation, Phys. Rev. A 80, 022303 (2009).
- (9) Childs, A.M., Farhi, E., and Preskill, J. Robustness of adiabatic quantum computation. Phys. Rev. A 65, 012322 (2001).
- (10) Roland, J. and Cerf, N.J. Noise resistance of adiabatic quantum computation using random matrix theory, Phys. Rev. A 71, 032330 (2005).
- (11) Sarandy, M.S., and Lidar, D.A. Adiabatic quantum computation in open systems, Phys. Rev. Lett. 95, 250503 (2005).
- (12) Tiersch, M., and Schützhold, R. Non-Markovian decoherence in the adiabatic quantum search algorithm, Phys. Rev. A 75, 062313 (2007).
- (13) Amin, M.H.S., Love, P.J., and Truncik, C.J.S. Themally assisted adiabatic quantum computation, Phys. Rev. Lett. 100, 060503 (2008).
- (14) Amin, M.H.S., Averin, D.V., and Nesteroff, J.A. Decoherence in adiabatic quantum computation, Phys. Rev. A 79, 022107 (2009).
- (15) Lloyd, S. Robustness of adiabatic quantum computing, eprint arXiv:0805.2757.
- (16) Uhlmann, A. The transition probability in the state space of a *-algebra, Rep. Math. Phys. 9, 273 (1976).
- (17) Allahverdyan, A.E., and Nieuwenhuizen, T.M. Extraction of work from a single thermal bath in the quantum regime, Phys. Rev. Lett. 85, 1799 (2000).
- (18) Nagaev, K.E., and Büttiker, M. Ground-state energy fluctuations of a system coupled to a bath, Europhys. Lett. 58, 475 (2002).
- (19) Scully, M.O., Zubairy, M.S., Agarwal, G.S., and Walther H. Extracting work from a single heat bath via vanishing quantum coherence, Science 299, 862 (2003).
- (20) Zhang, W., Sun, C.P., and Nori, F. Equivalence condition for the canonical and microcanonical ensembles in coupled spin systems, Phys. Rev. E 82, 041127 (2010).
- (21) Williams, N.S., Le Hur, K., and Jordan, A.N. Effective thermodynamics of strongly coupled qubits, J. Phys. A 44, 385003 (2011).
- (22) Cedraschi, P., Ponomarenko, V.V., and Büttiker, M. Zero-Point fluctuations and the quenching of the persistent current in normal metal rings, Phys. Rev. Lett. 84, 346 (2000).
- (23) Averin, D.V. and Pekola, J.P. Violation of the fluctuation-dissipation theorem in time-dependent mesoscopic heat transport, Phys. Rev. Lett. 104, 220601 (2010).
- (24) K. Blum, Density matrix theory and applications, (Plenum, New York, 1981).
- (25) Harris, R., et al. Experimental demonstration of a robust and scalable flux qubit, Phys. Rev. B 81, 134510 (2010).
- (26) Lanting, T., et al. Probing high-frequency noise with macroscopic resonant tunneling, Phys. Rev. B 83, 180502(R) (2011).
- (27) Dziarmaga J. Dynamics of a quantum phase transition: exact solution of the quantum Ising model, Phys. Rev. Lett. 95, 245701 (2005).
- (28) Lidar, D.A. Towards fault tolerant adiabatic quantum computation, Phys. Rev. Lett. 100, 160506 (2008).
- (29) Jordan, S., Farhi, E., and Shor, P. Error-correcting codes for adiabatic quantum computation, Phys. Rev. A 74, 052322 (2006).
- (30) Albash, T., Boixo, S., Lidar, D.A., and Zanardi P. Quantum adiabatic markovian master equations, New J. Phys. 14, 123016 (2012).
.2 Author contributions
All authors contributed equally to all aspects of this work.
.3 Competing financial interests
The authors declare no competing financial interests in relation to this work.