Superadiabatic Control of Quantum Operations
Abstract
Adiabatic pulses are used extensively to enable robust control of quantum operations. We introduce a new approach to adiabatic control that uses the superadiabatic quality or factor as a performance metric to design robust, high fidelity pulses. This approach permits the systematic design of quantum control schemes to maximize the adiabaticity of a unitary operation in a particular time interval given the available control resources. The interplay between adiabaticity, fidelity and robustness of the resulting pulses is examined for the case of singlequbit inversion, and superadiabatic pulses are demonstrated to have improved robustness to control errors. A numerical search strategy is developed to find a broader class of adiabatic operations, including multiqubit adiabatic unitaries. We illustrate the utility of this search strategy by designing control waveforms that adiabatically implement a twoqubit entangling gate for a model NMR system.
.1 Introduction
Speed and robustness are two essential characteristics of quantum control schemes that can often seem to be at odds with one another. On one hand, fast diabatic gates designed using optimal control techniques can approach the quantum speed limit (QSL) and minimize errors due to decoherence Caneva2009 (). However, these pulses are often sensitive to variations in the experimental control parameters and to uncertainties in the system Hamiltonian. Additionally, the pulse shapes produced by these techniques are typically not smooth and frequently push the hardware limits of the system, requiring careful tuning and calibration to ensure high fidelity Borneman2012 (); Martinis2014 (). On the other hand, smoothly varying gates can be made robust to control errors and are typically easier to implement due to the simpler hardware requirements. In particular, the utility of adiabatic gates that rely on the wellknown adiabatic theorem Messiah (1965) has been demonstrated for a variety of control tasks for quantum information processing Recati2002 (); Duan2003 (); Fuchs2009 (); Chen2014 (). The transitionfree driving of a quantum system enabled by adiabatic gates is particularly important if some excited states of the system are more susceptible to decoherence. Hybrid approaches that combine both diabatic and adiabatic control have also been explored Bacon2009 (); Hen2015 (); Chasseur2015 ().
Amplitude and frequencymodulated “adiabatic pulses” have long been used in nuclear magnetic resonance (NMR) to efficiently invert nuclear spin states Silver et al. (1984); Baum et al. (1985); Garwood and DelaBarre (2001) and provide robustness against inhomogeneities in both the static and radiofrequency (RF) magnetic fields, finding applications in both high resolution NMR spectroscopy and in vivo magnetic resonance imaging Tannús and Garwood (1997). Similar schemes have been used to optically control population transfers in atomic and molecular gases Allen and Eberly (1987); Gaubatz et al. (1990).
Finite time operations can only approximately satisfy the adiabatic condition, an issue that becomes critical in the context of adiabatic quantum computation Farhi2001 (). Long control pulses are also susceptible to decoherence introduced by interactions with unwanted environmental degrees of freedom. This raises an important question: what is the minimal time required to perform a highfidelity adiabatic transition? Counterdiabatic driving strategies — called shortcuts to adiabaticity (STA) — enable transitionless driving on much shorter timescales Torrontegui et al. (2012), even approaching the QSL Santos2015 (). These techniques have found applications in quantum state engineering Choi et al. (2011, 2012), quantum computing Sarandy et al. (2011); Herrera et al. (2014), manybody physics Saberi et al. (2014) and quantum simulations del Campo (2011); del Campo et al. (2012), and have been shown to have robustness against control parameter variations Bason2011 (). The DRAG pulses used in superconducting qubit implementations share many of these features as well Motzoi2009 (); Gambetta2011 (); Chow2010 (). One challenge to implementing counterdiabatic driving strategies, particularly for systems of more than onequbit, is that it may be difficult to generate the necessary counterdiabatic driving terms to ensure transitionfree evolution using the available experimental controls.
Here, we introduce a new approach to adiabatic control, based on Berry’s “superadiabatic” formalism Berry (1987), that enables the systematic design of quantum control schemes to maximize the adiabaticity of a unitary operation in a particular time interval, given the available controls. We explicitly use the superadiabatic quality or factor as a performance metric to optimize the available quantum control parameters. The idea of a superadiabatic factor was introduced by Deschamps et al. to explain the unexpectedly high fidelity of certain adiabatic pulses used in NMR Deschamps et al. (2008). We show that maximizing superadiabatic factors improves the performance of standard onequbit inversion pulses used in NMR and introduce a numerical search strategy to find a broader class of adiabatic unitaries when analytical solutions are not available. We numerically examine the interplay between adiabaticity, fidelity and robustness of the resulting pulses and show that superadiabatic pulses also improve robustness. Finally, we show how the search technique can be used to create control waveforms that adiabatically implement a twoqubit entangling gate. While we explore these ideas in the context of NMR experiments, the ideas are broadly applicable to other modalities.
.2 Superadiabatic Factors
Consider a timedependent Hamiltonian with instantaneous eigenbasis at time . Transforming to an interaction frame under the unitary operator that instantaneously diagonalizes the Hamiltonian yields an interaction frame Hamiltonian of the form , where is diagonal and is a nondiagonal correction term (called an inertial term) arising from the time dependence of the Hamiltonian. A transition has typically been considered adiabatic if or for the duration of the transition, where
(1) 
The “adiabatic factor” is then defined as
(2) 
For finitetime processes, is nonzero and the transformed Hamiltonian is nondiagonal. In many STA approaches, a counterdiabatic driving term is introduced to explicitly cancel this nondiagonal inertial term Torrontegui et al. (2012). Note that this is only possible if such an effective Hamiltonian can be generated with the available controls.
The above procedure for diagonalizing the instantaneous Hamiltonian can be applied to the transformed Hamiltonian , yielding a new Hamiltonian . Repeated indefinitely, this iterative procedure yields a countably infinite family of transformed Hamiltonians. Consider, for example, the Hamiltonian . If the set forms the instantaneous eigenbasis of , the unitary operator diagonalizes . In the interaction picture in which is instantaneously diagonalized, the Hamiltonian takes the form , where and . By direct analogy with Eqs. (1) and (2), the adiabatic factor in frame takes the form
(3) 
Counterdiabatic driving STA strategies can also be derived for superadiabatic interaction frames Ibáñez et al. (2012).
Deschamps et al. suggested that in a superadiabatic transformation, if the system starts out in one of the eigenstates of , it will evolve adiabatically to the target state in one of the superadiabatic frames as long as
(4) 
where is defined to be the superadiabatic factor Deschamps et al. (2008).
Scaling of
The metric shows two important features:
1. If is a timedependent Hamiltonian and , then for .
Proof: Let be the eigenvectors of . Then are eigenvectors of , and hence
from which we have that
and
Therefore
2. If where , then
for .
Proof: Let . Then is the unitary that diagonalizes , , and
Therefore
.3 Analytical NMR Inversion Pulses
To demonstrate the utility of the superadiabatic formalism, we examine the wellknown tanh/tan adiabatic inversion pulse, one of a family of single spin1/2 adiabatic inversion pulses used in NMR Hwang et al. (1998); Tannús and Garwood (1997). For this system, the Hamiltonian during the pulse in a reference frame rotating at the nuclear spin Larmor frequency () takes the form:
(5) 
where is the resonance offset, encodes the frequency and phase of the pulse, , is the nuclear gyromagnetic ratio, is the amplitude of the applied RF field, and and are Pauli spin operators. Here and throughout this paper, has been set to . The goal of the pulse is to invert the state from to .
For a tanh/tan pulse of length , the first half of the pulse () can be described by Hwang et al. (1998):
(6) 
and
(7) 
where corresponds to the maximum RF field strength, and , , and are parameters that can be optimized for a particular system. For the second half of the pulse (), and .
In the simulations here, the maximum RF amplitude was set at krad/sec (12.7 kHz), a typical value for a liquidstate NMR spectrometer. This corresponds to a minimum gate time of 39.27 s for a rectangular inversion pulse. The remaining three parameters (, , ) were numerically optimized using bruteforce search to generate pulses that either (a) maximized the traditional adiabatic factor , or (b) maximized the superadiabatic factor . Since for the pulse lengths examined, was calculated by computing the maximum value of the first ten factors, using the analytical forms derived for these pulses by Deschamps et al Deschamps et al. (2008). The optimal pulse parameters are shown in Table 1.
pulse  (rad/sec)  

6.9  16.1  
(120 s)  65.8  49.2  
(50 s)  36.3  41.6 
The optimization was first performed for pulse length s, about 3 times longer than the hardpulse time. Figure 1(a) compares the values of (where is defined by Eq. (3)) for the two pulses at this pulse length. For both optimized pulses, initially increases with until it reaches a peak value, which is the superadiabatic factor, ; in this case, for the optimized pulse, , while for the optimized pulse, . As Figure 1(a) shows, for , begins to decrease or “diverge,” a phenomenon that has been attributed to the finite time of the transition Berry (1987); Deschamps et al. (2008).
The overall fidelity of the pulse was characterized by the overlap . Figure 1(b) compares the performance, using the infidelity (), of the two optimized tanh/tan pulse shapes as their duration was changed from 0 to 250 s, demonstrating the improvement in fidelity provided by the superadiabatic pulse for pulse lengths s. Note the oscillations in the optimized pulse that occasionally give very high fidelity at certain times.
As a visual representation of the adiabatic dynamics, Figure 1(c) and (d) show how the Bloch vector corresponding to the state tracks the Hamiltonian of the optimized pulses on the Bloch sphere for the and the optimized pulses respectively. The timedependent Hamiltonian can also be represented as a vector on the Bloch sphere . Since for most values of , we plot the projection of onto the Bloch sphere instead of itself. The instantaneous deviation between and can be quantified in any superadiabatic frame by calculating the angle between and in that frame:
(8) 
Figures 1(e) and (f) show (dashed) and (solid) for the optimized pulse () and the optimized pulse () respectively. For the optimized pulse, and are on the same order of magnitude, accounting for the quantum state’s failure to reach the target state at this pulse length. For the optimized pulse, on the other hand, is negligible compared to , suggesting that the state is locked to the superadiabatic Hamiltonian, , but not to . For = 120 s, the infidelity of the pulse is seen to be quite large, which is reflected in both Figures 1(b) and (c).
We next examine the more general problem of engineering an optimally adiabatic pulse for a given pulse length . As shown earlier scales linearly with the length of the pulse if the pulse shape is held fixed, so a pulse shape that is optimal for a particular pulse length will remain optimal for all pulse lengths. Importantly, this property does not hold for higherorder factors, and hence a optimal pulse at one pulse length may not be optimal at a different pulse length, suggesting that a separate optimization needs to be performed for each pulse length of interest.
Figures 2(a) and (b) compare the performance of two optimized tanh/tan inversion pulses that were optimized for inversion times of 50 s and 120 s to the original optimized pulse. The pulse designed for 50 s is seen to perform better at shorter pulse lengths near 50 s (in terms of both fidelity and superadiabatic factor), while the pulse optimized for 120 s performs better at longer times, with the behavior appearing to switch around 77 s. The dotted vertical vertical line in Figure 2(a) indicates the duration of a “hard” rectangular pulse using the maximum available RF field of 80 krad/s, and the dashed line shows the fidelity achieved with this pulse. The fidelities of the three adiabatic pulses approach that of the ideal hard pulse at short times, but never exceed it. However the adiabaticity of the pulses is seen to rapidly fall as the pulse durations are reduced. For these single qubit inversion pulses, we found that preserved the desired robustness properties for the adiabatic pulses.
.4 Generalized Numerical Search Scheme
In the discussion above we considered the optimization of for singlespin pulses of a specific analytical form. In order to consider other unitaries, and to provide an optimization scheme that can be readily extended to higherdimensional spaces where closedform expressions for are generally not available, we have designed an evolutionary search strategy that iterates on an initial guess pulse to produce numerically optimized pulse shapes that maximize for any frame of interest . It should be noted that numerical optimization techniques have previously been used both to find the optimal pulse parameters of standard adiabatic NMR pulse shapes as well as to optimize arbitrary pulse shapes that maximize Baum et al. (1985); Silver et al. (1984); Garwood and DelaBarre (2001). The algorithm described below is similar to other derivativefree pulseshaping methods that have been used in the past Warren ().
We assume our Hamiltonian has the form , where is the timeindependent part of the Hamiltonian and are the control parameters corresponding to the Hermitian control operators . Let be the number of time steps used to define the pulse. An initial guess pulse (set of ) is chosen that satisfies the necessary boundary conditions at and to ensure that the initial and final states are eigenstates of and respectively. The key steps in our method are outlined here (see Appendix for additional details). (i) The parameters of the guess pulse are perturbed in a time interval and recalculated by numerically diagonalizing to find all the and as outlined above. Perturbations that improve are preserved and used to update the guess pulse. (ii) The center of the perturbation (), the size of the perturbed region () and the amplitude of the perturbation are all cycled systematically during the search as is maximized.
It is important to note that the evolutionary search does not guarantee convergence to a globally optimal pulse shape. As with many numerical search strategies, it is possible for the algorithm to get trapped in a local optimum. This may present a particular challenge as the size and complexity of the search space increases.
In Figure 3, this search technique has been applied to the case of onespin inversion. The chosen guess pulse consists of a linear ramp with arbitrarily chosen slope for the RF frequency offset and a parabola for the RF amplitude with zeros at the endpoints and a maximum value of at . The evolutionary algorithm was first applied to the guess pulse to maximize . The fidelity profile of the resulting pulse is plotted as a dasheddotted line in Figure 1, showing considerable improvement over the guess pulse fidelity. This optimized pulse was then used as the starting point for a second round of optimization, this time maximizing at a pulse length 50 s. The fidelity of the resulting pulse is also plotted in Figure 3 as a solid line. For comparison, the 50 s optimized tanh/tan pulse shown in Figure 2 is reproduced here as a dashed line. For pulse lengths around 50 s, the numerically optimized pulse outperforms the optimal tanh/tan pulse, demonstrating the potential benefits of numerical pulsefinding.
Optimal control techniques have also been used to maximize adiabaticity Rosenfeld1996 (); Brif2014 (). Previous work using optimal control techniques to find adiabatic pulses used the integral of as a global metric to search for pulses Rosenfeld1996 (). It should be possible to adapt such a technique to maximize the integral of the superadiabatic curve as well, which may enable the use of gradient based methods. However, it is uncertain whether maximizing the integral of will preserve transitionfree steering of the system at all times.
.5 Robustness against inhomogeneity
One of the principal benefits of adiabatic pulses is robustness against inhomogeneity in both the and terms of the Hamiltonian described by Eq. 4. Such robustness is desirable both for ensemble experiments in which there is a distribution of Hamiltonians (of both the system and control Hamiltonians either in space or in time), or if there is uncertainty in the Hamiltonian parameters. We consider here the performance of the  and optimized pulses discussed above when they are subjected to variations in both the frequency offset and the amplitude . Consider a onespin pulse described by the vector . We examine two distinct cases: (i) The RF amplitude is held fixed and a frequency offset term is added to , yielding the modified pulse
(ii) The original frequency offset is preserved, but the RF amplitude is multiplied by a scale factor , yielding the modified pulse
In Figure 4, the robustness of several of the pulses discussed above is examined for offsets krad/s and scalings . The dotted line corresponds to a hard pulse with an 80 krad/sec RF field and is the least robust of the pulses shown. The robustness of the optimal tanh/tan pulse is examined at pulse length 46 s, which is the shortest pulse length at which the pulse achieves nearly perfect fidelity (see Fig. 1(b) and 2(b)). The fidelity curves coincide for krad/s and , suggesting that, for this pulse length, the optimized pulse behaves like a hard pulse and confers little advantage in terms of robustness. The figure also compares the optimized numerical pulse with the pulse length also set at 46 s. Though the numerical pulse performs worse under ideal conditions, it achieves a higher fidelity than both the hard pulse and the optimal pulse for krad/s and . The advantages of the adiabatic pulses are more pronounced for longer pulse lengths. The fidelities of the optimized and 120 s optimized pulses at pulse length 120 s are plotted as dashed lines. While both exhibit robustness for a wide range of offsets, the superadiabatic tanh/tan pulse outperforms the optimized pulse for all and all . Furthermore, the superadiabatic pulse achieves nearly perfect fidelity for . This suggests that superadiabatic pulses offer an advantage not only in fidelity as a function of pulse length, as shown in Figures 1–3, but also in robustness against variations in the Hamiltonian parameters.
.6 Multiple qubits
Our approach can be extended, in principle, to a larger number of qubits. However, since it requires the diagonalization of the instantaneous Hamiltonian to optimize the trajectory, it is not a scalable approach, a property it shares with most optimal control schemes. We consider a twoqubit system whose Hamiltonian is given by
(9) 
where is the 2by2 identity operator, and are the qubit controls for qubits and respectively, and is a fixed coupling constant in units of Hz. This Hamiltonian arises in liquidstate NMR experiments and has also been implemented with superconducting qubits Plantenberg2007 (). Here we demonstrate the use of our numerical strategy to adiabatically evolve a nonentangled pure state to the maximally entangled Bell state without controlling . In some systems can also be a timedependent control when it can be experimentally varied Chen2014 ().
To design an adiabatic transition between and , we must first identify an initial Hamiltonian with eigenstate and a final Hamiltonian with eigenstate . Importantly, for the adiabatic theorem to hold, the two eigenstates must be nondegenerate and the order of the eigenstates must be preserved. Setting and requiring that and , can be written in matrix form as:
If and , the initial state is the eigenvector of corresponding to the secondlargest eigenvalue .
The condition on can be satisfied by setting and further requiring that and . Again, is chosen so that . In matrix form, with these conditions applied, becomes:
and the normalized eigenvector of with the secondlargest eigenvalue is the Bell state .
Simulating a liquid state NMR experiment, we used a fixed value of 209.4 Hz for the coupling, corresponding to the measured protoncarbon coupling in a carbon13 labeled chloroform sample. For the initial guess pulse the RF amplitudes (carbon) and (proton) were chosen to vary linearly from krad/s at time to krad/s ( kHz) at time . The resonance offsets and were also chosen to be linear, with . The values krad/s and krad/s were chosen to maximize (See Appendix for details). The search algorithm was then used to iterate on this initial guess to find a control sequence that maximizes . The algorithm was carried out at an arbitrary pulse length since scales linearly with the length of the pulse.
The shape, performance, and fidelity of the resulting pulse depend on how long the algorithm is allowed to iterate on the initial guess pulse. Here, a round of optimization is taken to be the number of times that each point in the pulse serves as a center of perturbation. Figure 5 shows the infidelity of the guess pulse and two optimized pulses (following one and three rounds of optimization) as the length of the pulse is varied, while the inset plots of both the guess pulse and optimized pulses, showing improvement in . This improvement in adiabaticity is matched by an improvement in fidelity, with the optimized pulse outperforming the guess pulse for many of the depicted pulse lengths. It is interesting to note that the QSL for a nonadiabatic gate in this twoqubit system is on the order of 1 ms with the same control resources (indicated by the dotted vertical line), which is significantly shorter than the high fidelity adiabatic pulses obtained here. The nonadiabatic entangling gate consists of pulses on both spins, followed by a delay , which is then followed by a pulse on the protons. The dashed line shows the drop in the fidelity of the nonadiabatic gate as the delay is reduced below . It is the small size of the coupling that necessitates long adiabatic gates in this case.
Figure 6 shows how the instantaneous eigenvalues of the system change during the evolution of the final pulses (optimized 3 times) with the values of and above. The second largest eigenvalue, corresponding to the transition under consideration, is plotted as a solid line. The figure confirms that the eigenvalues remain nondegenerate during the entire gate, with the size of the minimum energy gap set by the strength of the coupling. We plan to explore superadiabatic control of multiqubit systems in more detail in future work.
.7 Conclusions
We have introduced a new approach to transitionfree driving of quantum systems. This approach uses the superadiabatic factor as a performance metric to design robust, high fidelity pulses that maximize the adiabaticity of the quantum operation in a particular interval, given the available experimental controls. These smoothlyvarying superadiabatic pulses are also easier to implement due to their typically simpler hardware requirements.
For the case of single qubit inversion pulses, we found that optimizing instead of improved both fidelity and robustness over a wide range of pulse lengths. At shorter pulse lengths a tradeoff was observed between fidelity and robustness, in which pulses that perform at high fidelity near the quantum speed limit tended to be less robust against inhomogeneity in the control parameters. We also introduced a simple numerical search strategy to implement a broader class of adiabatic operations, including multiqubit adiabatic unitaries, and designed an adiabatic control sequence to implement a twoqubit entangling gate. Our investigation highlights the generality of the factor formalism, which can readily be extended to even larger Hilbert spaces or to systems characterized by entirely different Hamiltonians.
The proposed method promises to offer another option in the toolbox of quantum control techniques. Ultimately, it would be useful to systematically characterize the landscape of control in terms of potential tradeoffs between desirable features such as speed, robustness against control errors and adiabaticity (or transition free driving). This would allow experimentalists to tailor their control strategy to the specific experimental constraints in their setups.
.8 Acknowledgements
Jonathan Vandermause acknowledges support of a James O. Freedman Presidential Scholarship and a research grant from the Paul K. Richter and Evelyn E. Cook Richter Memorial Fund at Dartmouth College. This research is based in part upon work supported by the National Science Foundation under CHE1410504.
Appendix A Appendix
a.1 Search Algorithm Details
Figure 7 outlines the steps of the algorithm. Assume that the control pulses of length are divided into equal intervals such that . The control waveform is parameterized by , .

Choose the center of the perturbation () of the initial curve .

Choose the radius of perturbation . For each center of perturbation , the radii were allowed to vary from (alters the entire curve) to (smallest local perturbation).

Introduce a parabolic perturbation centered at with radius : For every point , is changed to such that
where is a constant that controls the size of the perturbation. For any given combination of and , we perturb the curve ten times (chosen arbitrarily), in each case choosing to be a random value in the interval . The figure shows 4 such perturbations.

The perturbation that maximizes the chosen adiabatic factor is preserved. If, for any given radius of perturbation, none of the ten perturbations improved the adiabaticity, we return to step (b), this time choosing a smaller perturbation radius.
When a perturbation that improves adiabaticity is found, the fourstep procedure is repeated for a new center (mod N). If the algorithm does not find an improvement for any of the radii between and , the center is changed. Finally, since a pulse consists of two functions, and , the algorithm toggled between the two: was perturbed at center , and before perturbing again with , was perturbed at center .
a.2 Initial guess pulses for twoqubit control
Figure 8 motivates the choice of (carbon, system ) and (hydrogen, system ) mentioned in the text. Each point in the 2D grid represents a different linear guess pulse. The colorbar indicates the value of – the lighter the color of the grid point, the higher is for the corresponding pulse. The highest values of occur when , which is forbidden by the adiabatic theorem since it leads to degeneracy in the eigenvalues. Instead, the point indicated by the black arrow was chosen, with krad/s and krad/s. Figure 9 plots the shape of this linear guess pulse as well as the optimized numerical pulse.
References
 T. Caneva, M. Murphy, T. Calarco, R. Fazio, S. Montangero, V. Giovannetti, G. E. Santoro, Phys. Rev. Lett., 103, 240501 (2009).
 T. W. Borneman and D. G. Cory, J. Magn. Reson., 225, 120–129 (2012).
 J. M. Martinis and M. R. Geller, Phys. Rev. A, 90, 022307 (2014).
 A. Messiah, Quantum Mechanics (NorthHolland, Amsterdam, 1965).
 A. Recati, T. Calarco, P. Zanardi, J. I. Cirac, and P. Zoller, Phys. Rev. A, 66, 032309 (2002).
 L.M. Duan and H. J. Kimble, Phys. Rev. Lett., 90, 253601 (2003).
 G.D. Fuchs, V.V. Dobrovitski, D.M. Toyli, F.J. Heremans, D.D. Awschalom, Science, 326, 1520–1522 (2009).
 Y. Chen, C. Neill, P. Roushan, N. Leung, M. Fang, R. Barends, J. Kelly, B. Campbell, Z. Chen, B. Chiaro, A. Dunsworth, E. Jeffrey, A. Megrant, J. Y. Mutus, P. J. J. O’Malley, C. M. Quintana, D. Sank, A. Vainsencher, J. Wenner, T. C. White, M. R. Geller, A. N. Cleland, and J. M. Martinis, Phys. Rev. Lett., 113, 220502 (2014).
 D. Bacon and S.T. Flammia, Phys. Rev. Lett., 103, 120504 (2009).
 I. Hen, Phys. Rev. A, 91, 022309 (2015).
 T. Chasseur, L. S. Theis, Y. R. Sanders, D. J. Egger, F. K. Wilhelm, Phys. Rev. A, 91, 043421 (2015).
 M. S. Silver, R. I. Joseph, and D. I. Hoult, Journal of Magnetic Resonance 59, 347 (1984).
 J. Baum, R. Tycko, and A. Pines, Phs. Rev. A 32, 6 (1985).
 M. Garwood and L. DelaBarre, Adv. Magn. Reson. 153–177, 155 (2001).
 A. Tannús and M. Garwood, NMR in Biomed. 10, 423 (1997).
 L. Allen and J. H. Eberly, Optical Resonance and TwoLevel Atoms (Dover, New York, 1987).
 U. Gaubatz, P. Rudecki, S. Schiemann, and K. Bergmann, J. Chem. Phys. 92, 5363 (1990).
 E. Farhi, J. Goldstone, S. Gutmann and M. Sipser, arXiv.org:quantph/0001106 (2001).
 E. Torrontegui, S. Ibáñez, S. MartínezGaraot, M. Modugno, A. del Campo, D. GuéryOdelin, A. Ruschhaupt, X. Chen, and J. Muga, Advances In Atomic, Molecular, and Optical Physics 62, 117–169 (2013) .
 A.C. Santos and M. S. Sarandy, Scientific Reports, 5, 15775 (2015).
 S. Choi, R. Onofrio, and B. Sundaram, Phys. Rev. A 84, 051601(R) (2011).
 S. Choi, R. Onofrio, and B. Sundaram, Phys. Rev. A 86, 043436 (2012).
 M. S. Sarandy, E. I. Duzzioni, and R. M. Serra, Phys. Lett. A 375, 3343 (2011).
 M. Herrera, M. S. Sarandy, E. I. Duzzioni, and R. M. Serra, Phys. Rev. A 89, 022323 (2014).
 H. Saberi, T. Opatrný, K. M. lmer, and A. del Campo, Phys. Rev. A 90, 060301(R) (2014).
 A. del Campo, Phys. Rev. A 84, 031606(R) (2011).
 A. del Campo, M. M. Rams, and W. H. Zurek, Phys. Rev. Lett. 109, 115703 (2012).
 M. G. Bason, M. Viteau, N. Malossi, P. Huillery, E. Arimondo, D. Ciampini, R. Fazio, V. Giovannetti, R Mannella, O. Morsch, Nature Physics, 8, 147, (2012).
 F. Motzoi, J.M. Gambetta, P. Rebentrost, F.K. Wilhelm, Phys. Rev. Lett., 103, 110501 (2009).
 J.M. Gambetta, F. Motzoi, S.T. Merkel, F.K. Wilhelm, Phys. Rev. A, 83, 012308, (2011).
 J.M. Chow, L. DiCarlo, J.M. Gambetta, F. Motzoi, L. Frunzio, S.M. Girvin, R.J. Schoelkopf, Phys. Rev. A, 82, 040305(R), (2010).
 M. Berry, Proc. R. Soc. Lond. A 414, 31 (1987).
 M. Deschamps, G. Kervern, D. Massiot, G. Pintacuda, L. Emsley, and P. J. Grandinetti, J. Chem. Phys. 129, 204110 (2008).
 S. Ibáñez, X. Chen, E. Torrontegui, J. Muga, and A. Ruschhaupt, Phys. Rev. Lett. 109, 100403 (2012).
 T.L. Hwang, P. C. van Zijl, and M. Garwood, J. Magn. Reson. 133, 200 (1998).
 W.S. Warren and M.S. Silver, Adv. Magn. Reson., 12, 247, (1988).
 D. Rosenfeld, Y. Zur, Magn. Reson. Med., 36, 401–409 (1996).
 C. Brif, M.D. Grace, M. Sarovar, K.C. Young, New J. Phys., 16, 065013 (2014).
 J.H. Plantenberg, P.C. de Groot, C.J.P.M. Harmans and J.E. Mooij, Nature, 447, 836 (2007).