# Robust dynamical decoupling

## Abstract

Quantum computers, which process information encoded in quantum mechanical systems, hold the potential to solve some of the hardest computational problems. A substantial obstacle for the further development of quantum computers is the fact that the life time of quantum information is usually too short to allow practical computation. A promising method to increase the life time, known as dynamical decoupling, consists of applying a periodic series of inversion pulses to the quantum bits. In the present review, we give an overview of this technique and compare different pulse sequences proposed earlier. We show that pulse imperfections, which are always present in experimental implementations, limit the performance of dynamical decoupling. The loss of coherence due to the accumulation of pulse errors can even exceed the perturbation from the environment. This effect can be largely eliminated by a judicious design of pulses and sequences. The corresponding sequences are largely immune to pulse imperfections and provide an increase of the coherence time of the system by several orders of magnitude.

###### pacs:

03.65.Yz,03.67.Pp,76.60.-k ,76.60.Lz## I Introduction

During the last decade, it was shown that quantum mechanical systems have the potential for processing information more efficiently than classical systems Shor (1994); DiVincenzo (1995); Nielsen and Chuang (2000); Bennett and DiVincenzo (2000). However, it remains difficult to realize this potential because quantum systems are extremely sensitive to perturbations. These perturbations come from external degrees of freedom or from the finite precision with which the systems can be realized and controlled by external fields. This loss of quantum information to the environment is called decoherence Zurek (2003). Different results show that a state is more sensitive to decoherence as the number of qubits increases Krojanski and Suter (2004, 2006a); Krojanski and Suter (2006b); Cho et al. (2006); Lovric et al. (2007); Sánchez et al. (2007); Doronin et al. (2011); Zobov and Lundin (2011). This is also manifested by the impossibility to time reverse a quantum evolution when an initially localized excitation spreads over a system Rhim et al. (1970); Zhang et al. (1992); G. Usaj and Levstein (1998); Pastawski et al. (2000). As the information is distributed over a increasing number of qubits, the evolution becomes more sensitive to the perturbation Pastawski et al. (2000); Jalabert and Pastawski (2001). In a similar vein, this sensitivity with the systems size limits the distance over which one can transfer information or analogously limits the number of qubits that one can control reliably Álvarez and Suter (2010, 2011). This is manifested as a localization effect for the quantum information Anderson (1958); Pomeransky and Shepelyansky (2004); Chiara et al. (2005); Burrell and Osborne (2007); Keating et al. (2007); Apollaro and Plastina (2007); Allcock and Linden (2009); Álvarez and Suter (2010, 2011). In order to overcome these limitations for allowing quantum information processing with large number of qubits, methods for reducing the decoherence effects have to be developed.

One can tackle the problem by correcting the errors generated by the perturbations, but this is only possible if the perturbation is small enough to keep the quantumness of the system Knill et al. (1998); Preskill (1998); Knill (2005). Therefore, one needs first to reduce the perturbation effects. The pioneering strategies for reducing decoherence were introduced in the Nuclear Magnetic Resonance (NMR) community, in particular by Erwin Hahn who showed that inverting a spin-1/2 system (a qubit) corresponds to an effective change of the sign of the perturbation Hamiltonian and therefore generates a time reversal of the corresponding evolution Hahn (1950). This leads to the formation of an echo that later was formalized as a Loschmidt echo Jalabert and Pastawski (2001). These manipulations were extended to the so-called decoupling methods Carr and Purcell (1954); Meiboom and Gill (1958); Waugh et al. (1968a); Maudsley (1986); Gullion et al. (1990), which disconnect effectively the environment.

In the context of this review, we describe the environment as a spin-bath, without loss of generality. Considering spins 1/2 as qubits, two different types of decoupling methods can be distinguished. In the first case, the qubit system is well distinguished from the environment. Its energy level splitting differs significantly from that of the bath. As a result, the coupling between them is much smaller than the difference of their energy level splittings. The interaction can then always be approximated by an Ising-type () interaction, which causes dephasing of the system qubit but no qubit flips. The decoupling methods required for these cases are called heteronuclear decoupling within the NMR community and they can involve manipulation only on the spin-system Hahn (1950); Carr and Purcell (1954); Meiboom and Gill (1958); Maudsley (1986); Gullion et al. (1990) or only at the environment J.S and Waugh (1982).

In the second case, the system and the environment have similar energy level splittings. This is the case of a homonuclear system where the general form of the coupling must be retained and it can induce flips of the system qubit as well as dephasing. In this case, decoupling will generally affect the complete system plus “bath” Waugh et al. (1968b); Mansfield et al. (1973); Rhim et al. (1973, 1974); Burum and Rhim (1979); Bur (1981).

In this review we focus on decoupling the system from the environment by applying control pulses only to the system. During the last years this technique has gathered a lot of interest because it requires relatively modest resources: it requires no overhead of information encoding, measurements or feedback. The method is known as dynamical decoupling (DD). Since its initial introduction Viola et al. (1999), a lot of effort has been invested to find sequences with improved error suppression Viola and Knill (2003); Khodjasteh and Lidar (2005); Uhrig (2007); Gordon et al. (2008); Biercuk et al. (2009); Uhrig (2009); Yang et al. (2010); Clausen et al. (2010); West et al. (2010a); Souza et al. (2011a). The optimal design of DD sequences depends of the different sources of errors that have to be eliminated. The first to be considered is the nature of the SE interaction, i.e., a pure dephasing when only the -spin operator is present, a pure spin-flip interaction ( and/or are present) or a general interaction with , and . Sequences like Carr-Purcell (CP) and Carr-Purcell-Meiboom-Gill (CPMG), which use rotations around a single axis, are useful when only two operators of , and are present, but in order to fight against a general interaction, pulses along different spatial direction have to be applied. The shortest sequence that fulfills this condition is the XY-4 sequence Maudsley (1986); Viola et al. (1999). An actual implementation must take into account, in addition to the above issues, the effect of pulse imperfections Viola and Knill (2003); Álvarez et al. (2010); Ryan et al. (2010); Souza et al. (2011a); Wang and Dobrovitski (2011); Xiao et al. (2011). Fighting the effect of pulse errors was in fact the original motivation for the development of the XY-4 sequence Maudsley (1986). Another experimental consideration is the amount of power deposited in the system, which often must be kept small to avoid heating effect or damage to the sample.

DD technique is becoming an important tool for quantum information processing Morton et al. (2006, 2008); Biercuk et al. (2009); Du et al. (2009); Álvarez et al. (2010); G. de Lange, et al. (2010); Ryan et al. (2010); Barthel et al. (2010); Bluhm et al. (2010); Ajoy et al. (2011); Naydenov et al. (2011); Souza et al. (2011a, b) as well as in spectroscopy Jenista et al. (2009); Bylander et al. (2011); Almog et al. (2011); Alvarez and Suter (2011) and imaging Taylor et al. (2008); Meriles et al. (2010); Hall et al. (2010); de Lange et al. (2011). In most cases, the goal is to preserve a given input state, but it may also be combined with gate operations Khodjasteh and Viola (2009); Khodjasteh et al. (2010); West et al. (2010b); Ng et al. (2011a); Boixo and Somma (2008). In many cases, experimental results show that one of the main limitations for improving the decoupling efficiency are the non-ideal properties of the decoupling pulses Álvarez et al. (2010); Ryan et al. (2010); Souza et al. (2011a, b). Decoherence effects due to the environment can in principle be reduced by shortening the cycle time . However, in reality the pulse lengths are finite and the DD cycle time is limited by the available hardware. It is therefore important to find the best finite time sequences for decoupling Viola and Knill (2003); Khodjasteh and Lidar (2007); Hodgson et al. (2010); Uhrig and Lidar (2010); Álvarez et al. (2010); Souza et al. (2011a, b); Ng et al. (2011b). Furthermore, the pulses do not only have finite lengths, they also do not implement perfect rotations. Thus, increasing the number of pulses can result in a large overall error that destroys the qubit coherence instead of preserving it. As a result, the performance of the decoupling sequence may have an optimum at a finite cycle time Álvarez et al. (2010); Souza et al. (2011a, b). In this review, we summarize the different DD strategies for fighting the effect of pulse imperfections and we show that they must be considered in the designing of useful DD sequences.

The paper is structured as follows: In section II we give some basics of dynamical decoupling, in section III we introduce the effects of pulse imperfections and in section IV we describe the strategies to fight against the imperfections effects. In the last section we draw some conclusions and perspectives.

## Ii Basics of dynamical decoupling

### ii.1 The system

We consider a single qubit as the system that is coupled to the bath. In a resonantly rotating frame of reference Abragam (1961), the free evolution Hamiltonian is

(1) |

where is the environment Hamiltonian and

(2) |

is a general interaction between the system and the environment. are operators of the environment and the SE coupling strength. The index runs over all modes of the environment. Dephasing is due to an interaction that affects the component of the spin-system operator, and spin-flips are done trough the and/or operators. A heteronuclear spin-spin interaction involves a pure dephasing interaction. This type of interaction that is naturally encountered when the system can be distinguished from the environment can be found in a wide range of solid-state spin systems, as for example nuclear spin systems in NMR Carr and Purcell (1954); Meiboom and Gill (1958); Álvarez et al. (2010); Ajoy et al. (2011), electron spins in diamonds Ryan et al. (2010), electron spins in quantum dots Hanson et al. (2007), donors in silicon Kane (1998), etc. In other cases when the system and environment have similar energies, the SE interaction can include terms along the , and axis.

### ii.2 DD sequences with a single rotation axis

DD is achieved by iteratively applying to the system a series of stroboscopic control pulses in cycles of period Viola et al. (1999). Over that period, the time-averaged SE interaction can be described by an averaged or effective Hamiltonian Haeberlen (1976). The goal of DD is the elimination of the effective SE interaction. This can be seen by looking at Hahn’s pioneering spin-echo experiment Hahn (1950) (See Fig. 1b). It is based on the application of a -pulse to the spin system at a time after the spins were left to evolve in the magnetic field. This pulse effectively changes the sign of the system-environment (SE) interaction, in this case the Zeeman interaction with the magnetic field. Letting the system to evolve for a refocusing period or time reversed evolution during the same duration generates the echo. If the magnetic field is static, the dynamics is completely reversed and the initial state of the spin recovered. This is observed as the echo. However if the magnetic field fluctuates, its effect cannot be reversed completely. Thus, the echo amplitude decays as a function of the refocusing time Hahn (1950); Carr and Purcell (1954). This decay contains information about the time-dependence of the environment.

To reduce the decay rate of the echo due to a time-dependent environment, Carr and Purcell introduced a variant of the Hahn-echo sequence, where the single -pulse is replaced by a series of pulses separated by intervals of duration Carr and Purcell (1954). This CP sequence reduces the changes induced by the environment if the pulse intervals are shorter than the correlation time of the environment. However, as the number of pulses increases, pulse errors tend to accumulate. Their combined effect can destroy the state of the system, rather than preserving it against the effect of the environment. This was noticed by Meiboom and Gill who proposed a modification of the CP sequence for compensating pulse errors, the CPMG sequence Meiboom and Gill (1958).

CP and CPMG sequences are useful only when two of the spin operators , and are affected by the environment. They can be written as

(3) |

where is a free evolution operator and is a pulse around the axis (and the analogous for ). The difference between the CP and CPMG sequence is the orientation of the rotation axis with respect to the initial condition. For applications in quantum information processing, this distinction is not relevant, since all operations have to be independent of the initial condition. CP or CPMG is the shortest sequence of pulses for decoupling a SE interaction that affects only two of the spin components , and Viola et al. (1999).

Usually the average Hamiltonian generated by DD sequences can be described by a series expansion, such as the Magnus expansionMagnus (1954). All the higher-order terms in this expansion describe imperfections, which reduce the fidelity of the sequence and should be eliminated. Improving the DD performance is therefore related to reducing the contribution of higher order terms. This is closely related to efforts for developing better decoupling sequences for NMR Waugh et al. (1968a). For quantum information processing (QIP), this lead to the design of sequences that make DD more effective, such as concatenated dynamical decoupling Khodjasteh and Lidar (2005, 2007). An important innovation was due to G. Uhrig Uhrig (2007), who proposed a sequence with non-equidistant pulse spacings, while all the standard sequences like CPMG are based on equidistant pulses. The UDD sequence is defined by

(4) |

where the delays are determined by the positions

(5) |

of the pulses with and . The lowest nontrivial order is equal to the CPMG sequence, UDD CPMG.

The effect of the non equidistant pulses can be discussed in the context of filter theory: DD can be considered as an environmental noise filter, where the distribution of pulses generates different filter shapes as a function of the frequencies. The overlap of this filter function with the spectral distribution of the environmental noise determines the decoherence rate Cywinski et al. (2008). Analogously, the filter shapes can be connected to diffraction patterns induced by interferences in the time domain Ajoy et al. (2011). The Uhrig DD sequence was shown theoretically to be the optimal sequence for reducing low frequency noise Uhrig (2007); Cywinski et al. (2008); Uhrig (2008). This prediction was confirmed experimentally Biercuk et al. (2009); Du et al. (2009); Jenista et al. (2009). However, it appears that non-equidistant sequences perform better only for particular noise spectral densities that increase for higher frequencies and have a strong cut off. In usual case, where the spectral density decrease smoothly with the frequency, as usually happens with NMR spin baths, equidistant sequences were predicted Cywinski et al. (2008); Uhrig (2008); Pasini and Uhrig (2010) and demonstrated Biercuk et al. (2009); Álvarez et al. (2010); G. de Lange, et al. (2010); Ryan et al. (2010); Barthel et al. (2010); Ajoy et al. (2011)to be the best option Ajoy et al. (2011).

This filter function description can be traced back to previous NMR approaches Garroway (1977) and to work on universal dynamical control Kofman and Kurizki (2001). Choosing the times for the pulses leads to a variety of sequences that can be optimized according to the spectral density of the bath Kofman and Kurizki (2001); Gordon et al. (2008); Biercuk et al. (2009); Uys et al. (2009); Clausen et al. (2010); Pasini and Uhrig (2010); Pan et al. (2010).

If the SE interaction is a pure dephasing one, it is sufficient to apply pulses in one direction. However because every experimental setup has finite precision, pulse errors create an effective general SE interaction Khodjasteh and Lidar (2007); Álvarez et al. (2010); Souza et al. (2011b). In this case, it was shown that sequences that are designed for general SE interaction perform better than 1D sequences Álvarez et al. (2010); Ryan et al. (2010); Souza et al. (2011a).

### ii.3 DD sequences with multiple rotation axis

If the system-environment interaction includes all three components of the system spin operator, decoupling can only be achieved if the sequence includes rotations around at least two different axes. The first decoupling sequence of this type is the XY-4 sequence, which alternates rotations around the x- and y-axes (see Fig. 1 d, ). This sequence was initially used to eliminate the effect of pulse errors in the CP and CPMG sequences Maudsley (1986). It is also the shortest sequence for DD for general SE interactions Viola et al. (1999). In quantum information processing, where we consider the initial state to be unknown, the CP and CPMG sequences, which correspond to a train of identical -pulses, are identical Álvarez et al. (2010). However, their effect on the quantum state depends on the (generally unknown) initial condition: If the initial condition is oriented along the rotation axis of the pulses, flip angle errors of the first pulse are refocused by the second pulse. However, for components perpendicular to the rotation axis, the pulse errors of all pulses add and cause rapid decay of the coherence, even in the absence of system-environment interactions Maudsley (1986); Álvarez et al. (2010). This motivated the development of the XY-4 sequence. In addition, pulse imperfections convert an ising-type system environment interaction into an effective general SE interaction Khodjasteh and Lidar (2007); Álvarez et al. (2010); Souza et al. (2011b), which is not eliminated by the CP/CPMG sequence, but is partially eliminated by the XY-4 sequence. In the QIP community, the XY-4 sequence is usually referred to as periodic dynamical decoupling (PDD).

The XY-4 sequence is also the building block for concatenated DD (CDD) sequences that improve the decoupling efficiency Khodjasteh and Lidar (2005, 2007). The CDD scheme recursively concatenates lower order sequences to increase the decoupling power. The CDD evolution operator for its original version for a recursion order of is given by

(6) |

where and =XY-4. Figure 1 shows a general scheme for this process. Each level of concatenation reduces the norm of the first non-vanishing order term of the Magnus expansion of the previous level, provided that the norm was small enough to begin with Khodjasteh and Lidar (2005, 2007). This reduction comes at the expense of an extension of the cycle time by a factor of four. If the delays between the pulses are allowed to be non-equidistant like in UDD, it becomes possible to create hybrid sequences, such as CUDD Uhrig (2009) and QDD West et al. (2010a); Quiroz and Lidar (2011); Kuo and Lidar (2011).

## Iii Effects of imperfections

Since the precision of any real operation is finite, the control fields used for decoupling introduce errors. Depending on the sequence, these errors can accumulate. If the number of pulses is large and the sequence is not properly designed, the accumulated pulse errors can lead to severe loss of coherence than the effect of the environment. Designing effective decoupling sequences that suppress environmental effects without degrading the system, even if the control fields have errors, requires a careful analysis of the relevant errors and appropriate strategies for combining rotations in such a way that the errors cancel rather than accumulate.

One example of a non-ideal control pulse is its finite duration, which implies a minimum achievable cycle time. The effects introduced by finite pulse lengths have been considered in different theoretical works Viola and Knill (2003); Khodjasteh and Lidar (2007); Hodgson et al. (2010). These works predict that high order CDD or UDD sequences in general lose their advantages when the delays between pulses or pulse length are strongly constrained. While the limitation on the cycle reduces the maximal achievable DD performance, pulse errors can be even more destructive. In most cases, the dominant cause of errors is a deviation between the ideal and the actual amplitude the of control fields. The result of this amplitude error is that the rotation angle experimentally implemented deviates from , typically by a few percent. Another important error occurs when the control field is not applied in resonance with the transition of the qubit. This off-resonant effect produces a rotation in which the flip angle and the rotation axis deviate from the ideal ones.

An example of the destructive effects of pulse imperfection is illustrated in the left panel of figure 2. Here, we measured the magnetization decay of the C nuclear spins during two different DD sequences. The sample used for this experiment was polycrystalline adamantane Álvarez et al. (2010); Ajoy et al. (2011); Souza et al. (2011a, b). The dephasing of the nuclear spins originates from the interaction with an environment consisting of H nuclear spins and can be considered as a pure dephasing process. The first sequence considered in the figure is CPMG. In this case we can observe that the decay of the magnetization is imperceptibly slow when the system is initially oriented parallel to the rotation axis of the pulse (longitudinal state). As we discuss below, this is an indication that the pulse errors have no effect on this initial state. In contrast, for a transverse initial state, the errors of the individual pulses accumulate and lead to a rapid decay, as shown in Fig. 2. A similar behavior is found for the UDD sequence Álvarez et al. (2010); Ajoy et al. (2011).

The second DD sequence considered in Fig. 2 is the XY-4 sequence, which consists of pulses applied along the x and y axes. The alternating phases of the pulses results in a partial cancellation of pulse errors, independent of the initial condition Maudsley (1986); Gullion et al. (1990). As a result, the performance of this sequence is much more symmetric with respect to the initial state in the xy-plane and the average decay times are significantly longer Álvarez et al. (2010); Souza et al. (2011a).

In the context of quantum information processing, it is important that the performance of gate operations be independent of the initial conditions (which typically are unknown). A common choice for quantifying the performance of a general quantum operation is then the fidelity Wang et al. (2008):

(7) |

Here, is the target propagator for the process and the actual propagator. For the present situation, where the goal is a quantum memory, the target propagator is the identity operation .

We can not assume that the actual propagators are unitary. We therefore write the process as

(8) |

where and are the density matrices at the beginning and end of the process. The operators must form a basis. For the present case, we choose them as . The ideal and actual processes can therefore be quantified by the matrix elements . For the target evolution, the -matrix is

(9) |

The matrix elements for the actual process are determined experimentally by quantum process tomography Chuang and Nielsen (1997); Nielsen and Chuang (2000). We use them to calculate the process fidelity from eq. (7). In figure 2, we compare the real part of the experimental matrices for the two sequences. After one cycle, the process matrices for both sequences represent a quantum operation that is close to the identity operation. The fidelity between the experimental matrices and the ideal matrix (9) is and for the CPMG and XY-4 cycles, respectively. Measured over 40 cycles, the process matrices of the two sequences differ significantly from the identity operation but also from each order. For the XY-4 sequence, the non-vanishing elements are and , while for CPMG the non-vanishing elements are and . The two matrices represent therefore qualitatively different processes. The XY-4 sequence destroys all transversal magnetization, in this case the resulting density matrix is

(10) |

This corresponds to a projection of the density operator onto the -axis, i.e. to a complete dephasing of the transverse components in the -plane.

The CPMG sequence, conversely, projects the density operator onto the -axis:

(11) |

The CPMG is a spin lock Santyr et al. (1988) sequence that allows retains magnetization in the direction but destroys all components perpendicular to it. Since a real experimental implementation always generates a distribution of control field amplitudes, spins at different positions precess with different rates around the direction of the rf field. As a result, the perpendicular components become completely randomized after a sufficiently large overall flip angle as shown in the left panel of Fig. 2.

The performance of experimentally accessible DD sequences is limited by the pulse errors Álvarez et al. (2010); Ryan et al. (2010); Souza et al. (2011a, b). In many situations, the dominating error contributions are flip angle and offset errors. In the next section we show different strategies to make decoupling sequences robust against such errors.

## Iv Robust decoupling sequences

We have to make decoupling insensitive to pulse imperfections. Different possibilities for generating high-fidelity sequences have been proposed in the context of quantum information processing Khodjasteh and Lidar (2005, 2007); Souza et al. (2011a). Here, we discuss two possible approaches: first we show that it is possible to replace individual refocusing pulses by compensated pulses that implement very precise inversions, and then we discuss sequences that are inherently robust, i.e. insensitive to the imperfections of the individual pulses.

### iv.1 Robust Pulses

The simplest approach to make a sequence robust is by replacing every standard pulse by a robust composite pulse Levitt (a). The composite pulses are sequences of consecutive pulses designed to be robust against various classes of imperfections generating therefore rotations that are close to the ideal rotation even in the presence errors. Particularly useful for quantum information applications are those composite pulses that produce compensated rotations for any initial condition, denominated in the NMR literature as class-A pulses Levitt (a).

Recent experiments have successfully used composite pulses to demonstrate the resulting increase of the performance of different DD sequences Ryan et al. (2010); Souza et al. (2011a). These works have implemented the composite pulse defined as:

(12) |

which is equivalent to a robust rotation around the axis defined by followed by a rotation around the axis Tycko et al. (1985). For cyclic sequences, which always consist of even numbers of rotations, the effect of the additional rotation vanishes if the flip angle errors are sufficiently small.

A comparison between standard sequences (not using robust pulses) against sequences with robust pulses has been reported in Souza et al. (2011a). It was observed that robust pulses improve the performance at high duty cycles. However, for low duty cycles, standard sequences perform better. This is due to the shorter cycle time of the standard sequence if constant duty cycles are compared. Thus, if the objective is only to preserve a quantum state, the best performance is obtained at high duty cycles, using robust pulses. The best sequences that are suitable for parallel application of quantum gate operations are the self-correcting sequences discussed below.

The composite pulses are usually designed to correct flip angle errors and offset errors. For compensating the effects introduced by the finite length of pulses, some theoretical works have proposed that a finite pulse could be approximated as an instantaneous one by using an appropriate shaped pulse Pryadko and Quiroz (2009); Uhrig and Pasini (2010); Pasini et al. (2011). These works has provided analytical and numerical evidence that a careful shape design can strongly affects the performance of decoupling sequences.

### iv.2 Self-correcting sequences

An alternative to the use of composite pulses consists in making the decoupling sequences fault-tolerant without compensating the error of each pulse, but by designing them in such a way that the error introduced by one pulse is compensated by the other pulses of the cycle Souza et al. (2011a). A straightforward strategy for designing improved sequences consists in concatenating one basic building block cycle into a longer and robust cycle.

The XY-4 cycle is often chosen as the building block for constructing self-correcting sequences. This cycle is the shortest DD sequence that cancels the zero-order average Hamiltonian for a general SE interaction and also has the advantage of being partially robust to pulse imperfections.

In the spectroscopy and quantum computing communities, two versions of the XY-4 sequence are used Souza et al. (2011b). The basic cycle originally introduced in NMR shows reflection symmetry with respect to the center of the cycle. In contrast to that, the sequence used in the quantum information community is time asymmetric. One consequence of this small difference is that in the symmetric version, the echoes are formed in the center of the windows between any two pulses, while in the case of asymmetric cycles, the echoes coincide with every second pulse, as shown in figure 3. The separation in time between the echoes is therefore twice as long in this case. If the environment is not static, the larger separation of the echoes leads to a faster decay of the echo amplitude Souza et al. (2011b).

If a robust pulse only contains rotations, as in the case of (12), it is also possible to convert such a composite pulse into a decoupling cycle by inserting delays between the individual rotations. This approach has been used in Souza et al. (2011a) to build an self-correcting sequence. The basic cycle is then

(13) |

The self-correcting sequence is created by combining 5-pulse blocks shifted in phase by , such as [KDD - KDD], where the lower index gives the overall phase of the block. The cyclic repetition of these 20 pulses is referred to as the KDD sequence Souza et al. (2011a).

In Fig. 4, we show how strongly errors in the flip angles of individual pulses affect the fidelity of the pulse sequence. Neglecting the effect of the environment, we calculate the fidelity after the application of 20 pulses as a function of the flip angle error. The figure compares the fidelities for the CPMG, XY-4 and KDD cycles. For the CPMG sequence, the fidelity drops to <95% if the flip angle error exceeds 2 %. For the XY-4 sequence, this bandwidth increases to 10 % and for KDD to 30 %. KDD and XY-4 are obviously much less susceptible to pulse imperfections than CPMG. The low fidelities observed for CPMG is experimentally manifested by the fast decay of the transverse components, such as in Fig. 2.

### iv.3 Combining basic cycles

Every decoupling sequence contains unwanted terms in the average Hamiltonian. They can be reduced by combining different versions of the basic cycles in such a way that some of the error terms cancel. Two different versions of this procedure have been used: The basic cycles can be applied subsequently Gullion et al. (1990) or one cycle can be inserted into the delays of another cycle Khodjasteh and Lidar (2005, 2007). The first approach was introduced in NMR, e.g. for designing high-performance homonuclear decoupling sequences Waugh et al. (1968b); Mansfield et al. (1973); Rhim et al. (1973, 1974); Burum and Rhim (1979); Bur (1981) or in high-resolution heteronuclear decoupling Hahn (1950); Carr and Purcell (1954); Meiboom and Gill (1958); Maudsley (1986); Gullion et al. (1990); J.S and Waugh (1982). Examples of DD sequences that can be constructed using this approach are XY-8 and XY-16 sequences Gullion et al. (1990). Here the XY-8 is created combining a XY-4 cycle with its time-reversed image while XY-16 is created combining the XY-8 with its phase-shifted copy.

The second approach is the concatenation scheme proposed by Khodjasteh and Lidar Khodjasteh and Lidar (2005, 2007). It generates the CDD sequence of order by inserting CDD cycles into the delays of the XY-4 sequence (see figure 1). Ideally, each level of concatenation improves the decoupling performance and the tolerance to pulse imperfections; in practice, higher order sequences do not always perform better. It has been theoretically predicted Viola and Knill (2003); Khodjasteh and Lidar (2007); Hodgson et al. (2010) and later observed experimentally that the finite duration of the pulses and constrained delays between pulses result in the existence of optimal levels of concatenation Álvarez et al. (2010); Souza et al. (2011a).

In figure 5, we demonstrate how the well-designed combination of basic cycles can lead to extended cycles with better error compensation. Here, we consider as the leading experimental imperfections deviations of the amplitude and frequency of the pulse. Neglecting the effect of the environment, we calculate the fidelity after applying 1680 pulses to the system as a function of the two error parameters. Each panel contains the color-coded fidelity for one of six different decoupling sequences. In the top panels we can clearly see the improvement in the error tolerance due to the CDD scheme of concatenation. Panels 1, 4 and 5 show the same result for the sequential concatenation scheme, where only two cycles are combined at each step: concatenation of the XY-4 cycle with its time-inverted and phase-shifted copies forms the XY-8 and XY-16 sequences. The 16-pulse XY-16 cycle is significantly more robust than the 84-pulse CDD cycle. The best performance is achieved by the KDD sequence, whose cycle consists of 20 pulses .

In Fig. 6, we compare the experimental performance of different self-correcting sequences. The performance of the CDD sequences always saturates or decreases with increasing duty cycle under the current experimental conditions Souza et al. (2011a). However, instead of saturating, the relaxation time for the KDD sequence continues to increase, as in the case of sequences with robust pulses Souza et al. (2011b). The KDD sequence combines the useful properties of robust sequences with those of sequences of robust pulses and can thus be used for both quantum computing and state preservation.

### iv.4 Time reversal symmetry

The symmetry of the basic building blocks has a key role in determining the performance of the concatenated higher order sequences. Two sequences constructed according to the same rules from a basic block have different propagators if the basic blocks are symmetric or not Souza et al. (2011a, b). If we concatenate four XY-4 cycles to the XY-16 sequence, e.g., we obtain new cycles, which are time-symmetric, independent of which version of the XY-4 sequence was used for the building blocks. Although all the odd order terms vanish in the average Hamiltonians of both versions, the even order terms of the sequences that are built from asymmetric blocks contain additional unwanted terms Souza et al. (2011b). The different behavior of sequences consisting of symmetric vs. asymmetric blocks is illustrated in figure 7. If we start from the symmetric form of XY-4, the resulting XY-16 sequence shows much better performance than the sequence using the asymmetric XY-4 as the building block. Analogous results were obtained for the two versions of the XY-8 sequence Souza et al. (2011b).

Earlier experiments showed two different contributions to the overall fidelity loss Souza et al. (2011b): A precession around the -axis, which can be attributed to the combined effect of flip-angle errors and an overall reduction of the amplitude, which results from the system-environment interaction. The combination of precession and reduction of amplitude is illustrated in right panel of Fig. 7. In this figure we show the xy-components of the magnetization at different times during the XY-16 sequence. If the XY-16 sequence is built by the asymmetric form of XY-4 a distinct precession around the -axis is observed. This causes a deviation from the desired evolution and reduces therefore the fidelity of the process. However, for the sequence consisting of symmetric blocks, the precession is negligible. These results suggest that pulse errors are better compensated by concatenating symmetric building blocks.

The same concept can also be applied to CDD sequences. The conventional concatenation scheme of eq. (6) uses asymmetric building blocks and is not compatible with the symmetric version of XY-4. A new concatenation scheme was therefore proposed in Souza et al. (2011a, b). In this scheme, the symmetrized version of CDD is constructed as

(14) |

In Fig. 8 we compare the process fidelities for the two versions of the CDD sequence. As in the case of XY sequences, clearly, the symmetrized version, CDD(S), shows a significantly improved performance, compared to the standard CDD(A) version. In Ref. Souza et al. (2011b), it was experimentally observed that the performance of all DD sequences based on symmetric building blocks is better or equal to that of sequences using non-symmetric building blocks. This behavior is consistent with general arguments based on average Hamiltonian theory Levitt (2008, b).

## V Conclusion and perspectives

Dynamical decoupling is becoming a standard technique for preserving the coherence of quantum mechanical systems, which does not need control over the environmental degrees of freedom. The technique aims to reduce decoherence rates by attenuating the system-environment interaction with a periodic sequence of pulses applied to the qubits. The pioneering strategies for decoupling were introduced in the context of NMR spectroscopy Hahn (1950). Since then, many different decoupling sequences have been developed in the context of NMR Carr and Purcell (1954); Meiboom and Gill (1958); Maudsley (1986); Gullion et al. (1990) or quantum information processing Viola et al. (1999); Khodjasteh and Lidar (2005, 2007); Uhrig (2007, 2009); West et al. (2010a); Quiroz and Lidar (2011); Kuo and Lidar (2011); Souza et al. (2011a).

Generally, we can divide the DD sequences in two groups: (i) Sequences that involve pulses in a single spatial direction and (ii) sequences that contain pulses in different spatial directions. The type (i) sequences are strongly sensitive to pulse errors and are only capable of suppressing the effects of a purely dephasing environment or pure spin-flip interaction. Examples of such DD sequences are CPMG and UDD. The second group can suppress a general SE interaction and usually exhibits better tolerance to experimental imperfections. Example of such sequences are the XY-family (XY-4, XY-8 and XY-16), the CDD sequences and the KDD sequence.

Recent experiments have successfully implemented DD methods and demonstrated the resulting increase of the coherence times by several orders of magnitude Álvarez et al. (2010); Ryan et al. (2010); Souza et al. (2011a, b). These works also showed that the main limitation to the reduction of the decay rates are the imperfections of the pulse. Two approaches have been used to correct this. The first approach replaces the inversion pulses by robust composite pulses Levitt (a), which generate rotations that are close to the target value even in the presence of pulse errors. In this case, the pulses are corrected individually. The second approach consists in designing the decoupling sequences in such a way that the error introduced by one pulse is compensated by the other pulses, without compensating the error of each pulse individually. The properties of basic decoupling cycles can be further improved by concatenating basic cycles into longer and more robust cycles. The concatenation can be made either by combining symmetry-related copies of a basic cycle subsequently Gullion et al. (1990) (resulting in the XY-8 and XY-16 sequences) or by inserting the basic cycle into the delays of another cycle Khodjasteh and Lidar (2005, 2007) (CDD sequences).

The time reversal symmetry of the basic building blocks is a useful criterion for minimizing error contributions. It has been demonstrated that the sequences built from symmetric building blocks often perform better and never worse than sequences built from non-symmetric blocks Souza et al. (2011a, b). This is a significant advantage, considering that the complexity of the sequences based on symmetric or asymmetric blocks are identical.

Earlier experiments Souza et al. (2011a) showed that the best sequences that are suitable for parallel application of quantum gate operations are the symmetric self-correcting sequences. However, as the delay between pulses decreases, sequences with robust pulses perform better. Thus, if the objective is only to preserve a quantum state, the best performance is achieved by using robust pulses to correct pulse errors. On the other hand, the KDD sequence introduced in Souza et al. (2011a) combines the useful properties of self-correcting sequence with those of robust pulses and can thus be used for both quantum memory and quantum computing.

During the last years, many advances have been achieved. However, for the application of the technique in real quantum devices, further studies will certainly be required. So far, most work has focused on single qubit systems. In the future, more experimental tests will be needed with multi-qubit systems. In the field of quantum computation, another important development may result from the combination of dynamical decoupling sequences with those techniques used to implement robust quantum gates Knill et al. (1998); Knill (2005); Bravyi and Kitaev (2005); Souza et al. (2011c). Since DD does not require auxiliary qubits or measurements, it can be used as an economical alternative to complement quantum error correction. Some theoretical works Khodjasteh and Viola (2009); Khodjasteh et al. (2010); West et al. (2010b); Ng et al. (2011a) proposed methods for combining the two methods but no experimental test was carried out to date. Future research on dynamical decoupling will also focus on applications outside of quantum information processing. Recent experiments have applied DD pulse sequences, for example, to probe the noise spectrum directly Bylander et al. (2011); Almog et al. (2011); Alvarez and Suter (2011) and detect weak magnetic fields Taylor et al. (2008); Meriles et al. (2010); Hall et al. (2010); de Lange et al. (2011).

###### Acknowledgements.

We acknowledge useful discussions with Daniel Lidar and Gregory Quiroz. This work is supported by the DFG through Su 192/24-1.### References

- P. W. Shor, in Proceedings of the 35th Annual Symposium on the Foundations of Computer Science, edited by S. Goldwasser (IEEE Computer Society Press, Los Alamitos, CA, 1994), p. 124.
- D. P. DiVincenzo, Science 270, 255 (1995).
- M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, Cambridge, 2000).
- C. H. Bennett and D. P. DiVincenzo, Nature 404, 247 (2000).
- W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003).
- H. G. Krojanski and D. Suter, Phys. Rev. Lett. 93, 090501 (2004).
- H. G. Krojanski and D. Suter, Phys. Rev. Lett. 97, 150503 (2006a).
- H. G. Krojanski and D. Suter, Phys. Rev. A 74, 062319 (2006b).
- H. Cho, P. Cappellaro, D. G. Cory, and C. Ramanathan, Phys. Rev. B 74, 224434 (2006).
- M. Lovric, H. Krojanski, and D. Suter, Phys. Rev. A 75, 042305 (2007).
- C. M. Sánchez, H. M. Pastawski, and P. R. Levstein, Physica B 398, 472 (2007).
- S. I. Doronin, E. B. Fel’dman, and A. I. Zenchuk, J. Chem. Phys. 134, 034102 (2011).
- V. E. Zobov and A. A. Lundin, J. Exp. Theor. Phys. 112, 451 (2011).
- W. K. Rhim, A. Pines, and J. S. Waugh, Phys. Rev. Lett. 25, 218 (1970).
- S. Zhang, B. H. Meier, and R. R. Ernst, Phys. Rev. Lett. 69, 2149 (1992).
- H. P. G. Usaj and P. Levstein, Mol. Phys. 95, 1229 (1998).
- H. M. Pastawski, P. R. Levstein, G. Usaj, J. Raya, and J. A. Hirschinger, Physica A 283, 166 (2000).
- R. A. Jalabert and H. M. Pastawski, Phys. Rev. Lett. 86, 2490 (2001).
- G. A. Álvarez and D. Suter, Phys. Rev. Lett. 104, 230403 (2010).
- G. A. Álvarez and D. Suter, Phys. Rev. A 84, 012320 (2011).
- P. Anderson, Phys. Rev. 109, 1492 (1958).
- A. Pomeransky and D. Shepelyansky, Phys. Rev. A 69, 014302 (2004).
- G. D. Chiara, D. Rossini, S. Montangero, and R. Fazio, Phys. Rev. A 72, 012323 (2005).
- C. K. Burrell and T. J. Osborne, Phys. Rev. Lett. 99, 167201 (2007).
- J. Keating, N. Linden, J. Matthews, and A. Winter, Phys. Rev. A 76, 012315 (2007).
- T. Apollaro and F. Plastina, Open Systems & Information Dynamics 14, 41 (2007).
- J. Allcock and N. Linden, Phys. Rev. Lett. 102, 110501 (2009).
- E. Knill, R. Laflamme, and W. H. Zurek, Science 279, 5349 (1998).
- J. Preskill, Proc. R. Soc. Lond. A 454, 385 (1998).
- E. Knill, Nature 434, 39 (2005).
- E. L. Hahn, Phys. Rev. 4, 580 (1950).
- H. Y. Carr and E. M. Purcell, Phys. Rev. 94, 630 (1954).
- S. Meiboom and D. Gill, Rev. Sci. Instrum. 29, 688 (1958).
- J. S. Waugh, C. H. Wang, L. M. Huber, and R. L. Vold, J. Chem. Phys. 48, 662 (1968a).
- A. A. Maudsley, J. Magn. Reson. 69, 488 (1986).
- T. Gullion, D. B. Baker, and M. S. Conradi, J. Magn. Reson. 89, 479 (1990).
- J.S and Waugh, J. Magn. Reson. (1969) 50, 30 (1982), ISSN 0022-2364.
- J. S. Waugh, L. M. Huber, and U. Haeberlen, Phys. Rev. Lett. 20, 180 (1968b).
- P. Mansfield, M. J. Orchard, D. C. Stalker, and K. H. B. Richards, Phys. Rev. B 7, 90 (1973).
- W. K. Rhim, D. D. Elleman, and R. W. Vaughan, J. Chem. Phys. 59, 3740 (1973), ISSN 00219606.
- W. K. Rhim, D. D. Elleman, L. B. Schreiber, and R. W. Vaughan, J. Chem. Phys. 60, 4595 (1974), ISSN 00219606.
- D. P. Burum and W. K. Rhim, J. Chem. Phys. 71, 944 (1979), ISSN 00219606.
- J. Magn. Reson. (1969) 44, 173 (1981), ISSN 0022-2364.
- L. Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 82, 2417 (1999).
- L. Viola and E. Knill, Phys. Rev. Lett. 90, 037901 (2003).
- K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 180501 (2005).
- G. S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007).
- G. Gordon, G. Kurizki, and D. A. Lidar, Phys. Rev. Lett. 101, 010403 (2008).
- M. J. Biercuk, H. Uys, A. P. VanDevender, N. Shiga, W. M. Itano, and J. J. Bollinger, Nature 458, 996 (2009).
- G. S. Uhrig, Phys. Rev. Lett. 102, 120502 (2009).
- W. Yang, Z. Wang, and R. Liu, arXiv: 1007.0623 (2010).
- J. Clausen, G. Bensky, and G. Kurizki, Phys. Rev. Lett. 104, 040401 (2010).
- J. R. West, B. H. Fong, and D. A. Lidar, Phys. Rev. Lett. 104, 130501 (2010a).
- A. M. Souza, G. A. Álvarez, and D. Suter, Phys. Rev. Lett. 106, 240501 (2011a).
- G. A. Álvarez, A. Ajoy, X. Peng, and D. Suter, Phys. Rev. A 82, 042306 (2010).
- C. A. Ryan, J. S. Hodges, and D. G. Cory, Phys. Rev. Lett. 105, 200402 (2010).
- Z. Wang and V. V. Dobrovitski, J. Phys. B: At., Mol. Opt. Phys. 44, 154004 (2011).
- Z. Xiao, L. He, and W.-g. Wang, Phys. Rev. A 83, 032322 (2011).
- J. J. L. Morton, A. M. Tyryshkin, A. Ardavan, S. C. Benjamin, K. Porfyrakis, S. A. Lyon, and G. A. D. Briggs, Nature Phys. 2, 40 (2006).
- J. J. L. Morton, A. M. Tyryshkin, R. M. Brown, S. Shankar, B. W. Lovett, A. Ardavan, T. Schenkel, E. E. Haller, J. W. Ager, and S. A. Lyon, Nature 455, 1085 (2008).
- J. Du, X. Rong, N. Zhao, Y. Wang, J. Yang, and R. B. Liu, Nature 461, 1265 (2009).
- G. de Lange, et al., Science 330, 60 (2010).
- C. Barthel, J. Medford, C. M. Marcus, M. P. Hanson, and A. C. Gossard, Phys. Rev. Lett. 105, 266808 (2010).
- H. Bluhm, S. Foletti, I. Neder, M. Rudner, D. Mahalu, V. Umansky, and A. Yacoby, arXiv:1005.2995 (2010).
- A. Ajoy, G. A. Álvarez, and D. Suter, Phys. Rev. A 83, 032303 (2011).
- B. Naydenov, F. Dolde, L. T. Hall, C. Shin, H. Fedder, L. C. L. Hollenberg, F. Jelezko, and J. Wrachtrup, Phys. Rev. B 83, 081201 (2011).
- A. M. Souza, G. A. Álvarez, and D. Suter, arXiv:1110.1011 (2011b).
- E. R. Jenista, A. M. Stokes, R. T. Branca, and W. S. Warren, J. Chem. Phys. 131, 204510 (2009).
- J. Bylander, S. Gustavsson, F. Yan, F. Yoshihara, K. Harrabi, G. Fitch, D. G. Cory, Y. Nakamura, J.-S. Tsai, and W. D. Oliver, Nat. Phys. 7, 565 (2011).
- I. Almog, Y. Sagi, G. Gordon, G. Bensky, G. Kurizki, and N. Davidson, J. Phys. B: At. Mol. Opt. Phys. 44, 154006 (2011).
- G. A. Alvarez and D. Suter, arXiv:1106.3463v2 (2011).
- J. M. Taylor, P. Cappellaro, L. Childress, L. Jiang, D. Budker, P. R. Hemmer, A. Yacoby, R. Walsworth, and M. D. Lukin, Nat. Phys. 4, 810 (2008).
- C. A. Meriles, L. Jiang, G. Goldstein, J. S. Hodges, J. Maze, M. D. Lukin, and P. Cappellaro, J. Chem. Phys. 133, 124105 (2010).
- L. T. Hall, C. D. Hill, J. H. Cole, and L. C. L. Hollenberg, Phys. Rev. B 82, 045208 (2010).
- G. de Lange, D. Ristè, V. V. Dobrovitski, and R. Hanson, Phys. Rev. Lett. 106, 080802 (2011).
- K. Khodjasteh and L. Viola, Phys. Rev. Lett. 102, 080501 (2009).
- K. Khodjasteh, D. A. Lidar, and L. Viola, Phys. Rev. Lett. 104, 090501 (2010).
- J. R. West, D. Lidar, B. H. Fong, and M. F. Gyure, Phys. Rev. Lett. 105, 230503 (2010b).
- H. K. Ng, D. A. Lidar, and J. Preskill, Phys. Rev. A 84, 012305 (2011a).
- S. Boixo and R. D. Somma, Phys. Rev. A 77, 052320 (2008).
- K. Khodjasteh and D. A. Lidar, Phys. Rev. A 75, 062310 (2007).
- T. E. Hodgson, L. Viola, and I. D’Amico, Phys. Rev. A 81, 062321 (2010).
- G. S. Uhrig and D. A. Lidar, Phys. Rev. A 82, 012301 (2010).
- H. K. Ng, D. A. Lidar, and J. Preskill, Phys. Rev. A 84, 012305 (2011b).
- A. Abragam, Principles of Nuclear Magnetism (Oxford University Press, London, 1961).
- R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2007).
- B. E. Kane, Nature 393, 133 (1998).
- U. Haeberlen, High Resolution NMR in Solids: Selective Averaging (Academic Press, New York, 1976).
- W. Magnus, Commun. Pure Appl. Math. 7, 649 (1954).
- L. Cywinski, R. M. Lutchyn, C. P. Nave, and S. DasSarma, Phys. Rev. B 77, 174509 (2008).
- G. S. Uhrig, New J. Phys. 10, 083024 (2008).
- S. Pasini and G. S. Uhrig, Phys. Rev. A 81, 012309 (2010).
- A. N. Garroway, J. Magn. Reson. 28, 365 (1977).
- A. G. Kofman and G. Kurizki, Phys. Rev. Lett. 87, 270405 (2001).
- H. Uys, M. J. Biercuk, and J. J. Bollinger, Phys. Rev. Lett. 103, 040501 (2009).
- Y. Pan, Z. Xi, and W. Cui, Phys. Rev. A 81, 022309 (2010).
- G. Quiroz and D. A. Lidar, Phys. Rev. A 84, 042328 (2011).
- W.-J. Kuo and D. A. Lidar, Phys. Rev. A 84, 042329 (2011).
- X. Wang, C.-S. Yu, and X. Yi, Phys, Lett. A 373, 58 (2008).
- I. L. Chuang and M. A. Nielsen, J. Mod. Opt. 44, 2455 (1997).
- G. E. Santyr, R. Henkelman, and M. J. Bronskill, J. Magn. Reson. (1969) 79, 28 (1988), ISSN 0022-2364.
- M. H. Levitt, Composite pulses, in Encyclopedia of NMR, Eds. D. M. Grant and R. K. Harris (Wiley, 1996).
- R. Tycko, A. Pines, and J. Guckenheimer, J. Chem. Phys. 83, 2775 (1985).
- L. P. Pryadko and G. Quiroz, Phys. Rev. A 80, 042317 (2009).
- G. S. Uhrig and S. Pasini, New J. Phys. 12, 045001 (2010).
- S. Pasini, P. Karbach, and G. S. Uhrig, EPL (Europhysics Letters) 96, 10003 (2011).
- M. H. Levitt, J. Chem. Phys. 128, 052205 (2008).
- M. H. Levitt, Symmetry-based pulse sequences in magic-angle spinning solid-state nmr, in Encyclopedia of NMR (Wiley, 2002).
- S. Bravyi and A. Kitaev, Phys. Rev. A 71, 022316 (2005).
- A. M. Souza, J. Zhang, and R. Laflamme, Nat. Commun. 2, 169 (2011c).