## Adiabatic tracking of quantum many-body dynamics

Hamed Saberi, 1, 2, ∗ Tomáš Opatrný, 1 Klaus Mølmer, 3 and Adolfo del Campo 4, 5, 6

1 Department of Optics, Faculty of Science, Palacký University, 17. listopadu 12, 77146 Olomouc, Czech Republic

2 Department of Physics and CeOPP, University of Paderborn,

Warburger Straße 100, D-33098 Paderborn, Germany

3 Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark

4 Department of Physics, University of Massachusetts Boston, Boston, MA 02125, USA

5 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA

6 Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA

(Dated: December 5, 2014)

The nonadiabatic dynamics of a many-body system driven through a quantum critical point can be controlled using counterdiabatic driving, where the formation of excitations is suppressed by assisting the dynamics with auxiliary multiple-body nonlocal interactions. We propose an alternative scheme which circumvents practical challenges to realize shortcuts to adiabaticity in mesoscopic systems by tailoring the functional form of the auxiliary counterdiabatic interactions. Adriving scheme resorting in short-range few-body interactions is shown to generate an effectively adiabatic dynamics.

PACS numbers: 03.67.Ac, 64.60.Ht, 05.30.Rt, 37.10.Ty

Introduction. The adiabatic driving of quantum manybody systems is of interest to a wide variety of quantum technologies ranging from quantum simulation to adiabatic quantum computation. However, the implementation of adiabatic driving schemes in the laboratory is often challenging or simply impractical. In recent years a large body of theoretical and experimental progress has been focused on the development of shortcuts to adiabaticity (STA), fast-nonadiabatic protocols that reproduce the preparation of the same final state that would be achieved under slow driving [1]. The feasibility of realizing such shortcuts relies on the control of nonadiabatic excitations away from the ground state manifold of the system of interest. A general technique to achieve this goal is known as counterdiabatic driving (CD) [2-4]. In a nutshell, the adiabatic approximation | Ψ( ) t 〉 to the dynamics generated by a slowly-driven Hamiltonian of interest H λ t 0 [ ( )] , is found to be the exact solution of the so-called counterdiabatic driving Hamiltonian H CD , i.e., i /planckover2pi1 ∂ t | Ψ( ) t 〉 = H CD | Ψ( ) t 〉 , even under fast driving. More precisely, let the instantaneous eigenstates and eigenvalues of H λ t 0 [ ( )] be denoted by {| ε n [ λ t ( )] 〉} and { ε n [ λ t ( )] } , then the CD Hamiltonian can be always written as the sum H CD = H λ t 0 [ ( )] + H aux [ λ t ( )] of the system Hamiltonian H λ t 0 [ ( )] and the auxiliary CD term

$$H _ { \text{aux} } & = i \hbar { \dot { \lambda } } ( t ) \sum _ { n } [ | \partial _ { \lambda } \varepsilon _ { n } \rangle \langle \varepsilon _ { n } | - \langle \varepsilon _ { n } | \partial _ { \lambda } \varepsilon _ { n } \rangle | \varepsilon _ { n } \rangle \langle \varepsilon _ { n } | ] \, \ ( 1 ) \quad \text{pres} \\. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.$$

where the over-dot denotes the time derivative. This auxiliary term suppresses explicitly excitations away from the adiabatic manifold allowing one to design STA.

∗ saberi@optics.upol.cz

Conversely, it vanishes as the adiabatic limit is approached [4, 5]. When applied to lattice systems [5], CD is closely related to the notion of quasiadiabatic continuation [6-8] and is equivalent to transitionless quantum driving [9].

To date, CD has been demonstrated in effective twolevel systems [10-12] but its general implementation is expected to be highly challenging in many-body systems. A remarkable experimental demonstration in a low-dimensional quantum fluid has recently been reported [13], but overall the technique appears to be restricted to self-similar processes [14-16]. For few-particle systems additional progress has targeted the design of experimentally realizable variants of the CD Hamiltonian spanned by a set of implementable control fields [17].

Mimicking adiabatic dynamics is particularly challenging in many-body systems exhibiting a quantum critical point (QCP) λ c as a function of an external parameter λ . In the neighborhood of λ c , the characteristic relaxation time and correlation length exhibit a powerlaw divergence as a function of the reduced parameter /epsilon1 = [ λ -λ c ] /λ c . As a result, the dynamics across QCP is expected to result in the formation of topological defects, signaling the breakdown of adiabatic dynamics [18-20]. Signatures of critical dynamics are still present in finite systems and are masked by finite-size effects as the system size is reduced. Approaches to suppress defect formation in critical systems include nonlinear quenches [21, 22], inhomogeneous driving [23, 24], optimal control [25, 26] and multi-parameter tuning [27], among other examples, see for a brief summary [28]. It has been shown that CD can be used to assist the adiabatic dynamics across a quantum phase transition at the cost of engineering H aux which involves multiple-body non-local counterdiabatic interactions [5]. As a result, implementing H aux in the laboratory remains challeng-

ing despite recent progress in digital quantum simulation [29, 30].

In this Rapid Communication, we propose a practical recipe for suppression on demand of the nonadiabatic transitions in arbitrary many-particle systems. We illustrate the applicability of our method in the context of spin chains. Our scheme is of particular relevance to state preparation in digital quantum simulators, e.g. in trapped ion chains [29-31] and Rydberg gases [32].

Counterdiabatic driving of many-body systems. We consider a quench of a finite transverse Ising chain of N spins with nearest-neighbor interactions under open boundary condition described by

$$H _ { 0 } ( t ) = - B ( t ) \sum _ { j = 1 } ^ { N } \sigma _ { j } ^ { x } + J _ { 0 } \sum _ { j = 1 } ^ { N - 1 } \sigma _ { j } ^ { z } \otimes \sigma _ { j + 1 } ^ { z } \,, \quad ( 2 ) \quad \begin{pmatrix} \tan \alpha \\ \alpha _ { \ell } = \\ \text{mal} \, \mathbb { m } \, \mathbb { m } \end{pmatrix}$$

where B t ( ) is a time-varying external magnetic field, σ , σ x j z j are usual Pauli sigma operators associated with site j , and J 0 denotes uniform magnetic couplings between adjacent spins assumed to be unity throughout the work. In the thermodynamic limit the system exhibits a quantum phase transition as a function of the transverse field, separating a paramagnetic phase ( | B t /J ( ) 0 | &gt; 1 ) from a doubly-degenerate antiferromagnetic one ( | B t /J ( ) 0 | &lt; 1 ).

The quantum simulation of this model and its variants is the goal of current efforts in ion-trap experiments [3336], ultracold atoms [37], and NMR experiments [38], to name just a few examples. The achievable system-size with current technology is still tractable by real-space exact diagonalization (ED), which we employ to derive the instantaneous eigenspectrum {| ε n ( ) t 〉} and the auxiliary CD term via Eq. (1). H aux is found to involve K -

෩

ሺሻ

ॆ

࣬

glyph&lt;c=3,font=/YYLTOZ+Cambria#20Math-KSCms-UHC-H&gt;glyph&lt;c=3,font=/YYLTOZ+Cambria#20Math-KSCms-UHC-H&gt; glyph&lt;c=3,font=/YYLTOZ+Cambria#20Math-KSCms-UHC-H&gt;glyph&lt;c=3,font=/YYLTOZ+Cambria#20Math-KSCms-UHC-H&gt;

࣢

݅

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0524v2-figure-bb213ae1-2d90-4e8b-9471-f8f0b967b1c6.jpg)

݅

݅

Figure 1. (Color online) (a) Schematic of the variational construction of the ansatz of the form Eq. (3) to the 'exact' auxiliary CD term in Eq. (1) for possible experimental implementation in trapped ion chain. ˜ H [ K , R ] aux corresponds to the action of an itinerant 'train' of K -partite Pauli interactions between all possible permutations of K -tuple of spins within the range R and with optimal amplitudes h α ,...,α 1 K i 1 ,...,i K obtained from the contraction pattern depicted in (b).

body interactions ( K = 1 , . . . , N ), in agreement with [5]. We wish to find an alternative CD scheme resorting exclusively in auxiliary control fields associated with a restricted set of operators, assumed to be available in a quantum simulator. The scheme relies on suggesting a new ansatz for the variational construction of the 'exact' auxiliary CD term of the form

$$\begin{array} { c } \mathcal { H } _ { \text{a} } ^ { [ \mathfrak { K }, \mathcal { R } ] } ( t ) = \sum _ { i _ { 1 }, \dots, i _ { \mathfrak { A } } } ^ { \prime } \sum _ { \alpha _ { 1 }, \dots, \alpha _ { \mathfrak { A } } } h _ { i _ { 1 }, \dots, i _ { \mathfrak { A } } } ^ { \alpha _ { 1 }, \dots, \alpha _ { \mathfrak { A } } } ( t ) \ \bigotimes _ { \ell = 1 } ^ { \mathfrak { K } } \sigma _ { i _ { \ell } } ^ { \alpha _ { \ell } } \, \ ( 3 ) \\ \mathfrak { e } \quad \ \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \end{array}$$

where 0 &lt; i | /lscript -i /lscript ′ | ≤ R ∀ /lscript, /lscript ′ is understood in the restricted sum over distinct site indices at maximal distance R = { K -1 , . . . , N -1 } from each other and to be identified as the range of the K -body interaction, α /lscript = { 0 , x, y, z } with σ 0 ≡ 1 , and h α ,...,α 1 K i 1 ,...,i K are the optimal interaction amplitudes to be found via a variational optimization procedure that shall be detailed in the following. This ansatz is general enough to generate an arbitrary unitary. However, we shall see that a ˜ H [ K , R ] aux containing exclusively few-body quasi-local terms suffices to induce an effectively adiabatic dynamics. ˜ H [ K , R ] aux involves nonlocal K -body couplings of spins through the chain as Fig. 1(a) illustrates. Considering the driving of an arbitrary eigenstate | Ψ(0) 〉 = | ε n (0) 〉 prepared at t = 0 , the optimal amplitudes h α ,...,α 1 K i 1 ,...,i K ( ) t that lead to an optimal representation of H aux within the subspace of the ansatz (3) may be obtained by minimizing the quadratic cost function of the form

$$\min _ { h _ { i _ { 1 }, \dots, i _ { \mathcal { A } } } ^ { \alpha _ { 1 }, \dots, \alpha _ { \mathcal { B } } } ( t ) } \| ( H _ { \text{aux} } - \tilde { \mathcal { H } } _ { \text{aux} } ^ { [ \mathcal { R }, \mathcal { R } ] } ) | \Psi _ { \text{a} } ( t ) \rangle \| ^ { 2 } \,, \quad ( 4 )$$

where | Ψ ( ) a t 〉 ≡ | ε n ( ) t 〉 denotes the adiabatic evolution of | Ψ(0) 〉 generated by H 0 ( ) t . The result of such a minimization can generally be written down in the form

$$\sum _ { \{ i _ { \ell } \} } ^ { ^ { \prime } } \sum _ { \{ \alpha _ { \ell } \} } A _ { i _ { 1 } ^ { \prime }, \dots, i _ { \alpha } ^ { \prime }, i _ { 1 } ^ { \prime }, \dots, i _ { \alpha } } ^ { \alpha _ { 1 } ^ { \prime }, \dots, \alpha _ { \beta } ^ { \prime }, \alpha _ { 1 }, \dots, \alpha _ { \beta } } h _ { i _ { 1 }, \dots, i _ { \beta } } ^ { \alpha _ { 1 }, \dots, \alpha _ { \beta } } = C _ { i _ { 1 } ^ { \prime }, \dots, i _ { \alpha } ^ { \prime } } ^ { \alpha _ { 1 } ^ { \prime }, \dots, \alpha _ { \beta } ^ { \prime } } \,, \ ( 5 ) \$$

where

$$& \text{where} \\ & A ^ { \alpha ^ { \prime } _ { 1 }, \dots, \alpha ^ { \prime } _ { \tilde { \lambda } }, \alpha _ { 1 }, \dots, \alpha _ { \tilde { \lambda } } } _ { i ^ { \prime } _ { 1 }, \dots, i ^ { \prime } _ { \tilde { \lambda } }, i ^ { \prime } _ { 1 }, \dots, i _ { \tilde { \lambda } } } = \text{Tr} \left [ \rho ^ { \frac { \tilde { \lambda } } { t } } \{ \bigotimes _ { \ell = 1 } ^ { \tilde { \lambda } } \sigma ^ { \alpha _ { \ell } } _ { i _ { \ell } }, \bigotimes _ { \ell = 1 } ^ { \tilde { \lambda } } \sigma ^ { \alpha ^ { \prime } _ { \ell } } _ { i ^ { \prime } _ { \ell } } \} \right ], \\ & C ^ { \alpha _ { 1 }, \dots, \alpha _ { \tilde { \lambda } } } _ { i _ { 1 }, \dots, i _ { \tilde { \lambda } } } = \text{Tr} \left [ \rho ^ { \frac { \tilde { \lambda } } { t } } \{ H _ { \text{aux} }, \bigotimes _ { \ell = 1 } ^ { \tilde { \lambda } } \sigma ^ { \alpha _ { \ell } } _ { i _ { \ell } } \} \right ] \,, \\ & \text{where the anticommutator} \{ A, B \} = A B + B A \text{ and we}$$

where the anticommutator { A,B } = AB + BA and we have denoted ρ a t = Ψ ( ) | a t 〉〈 Ψ ( ) a t | to stress the flexibility of the formalism in extending to mixed states. At a given time t , the equations above admit a compact form as a single linear tensor equation involving a contraction of the form

$$\text{Ahi} = C \,,$$

with a graphical representation depicted in Fig. 1(b).

The latter tensor equation can be solved by an in situ reshaping of the tensor A of rank 4 K into a regular matrix ¯ A of dimension (4 N ) K × (4 N ) K while combining the set of indices { i /lscript }{ α /lscript } and { i ′ /lscript }{ α ′ /lscript } into single 'superindices' I and I ′ , respectively. This amounts to casting the original tensor equation into a matrix inversion problem involving matrix ¯ A I ,I ′ and reshaped column vectors ¯ h I and ¯ C I ′ given by ¯ = h ¯ A -1 ¯ C . Reshaping back properly the column vector ¯ h so obtained into the original order yields the optimal amplitudes h α ,...,α 1 K i 1 ,...,i K for approximating the exact auxiliary CD term. A costly part of the implementation of the outlined procedure requires an explicit construction of ¯ A and ¯ C which in case of the driving of pure states can be economized on by utilizing the following properties and bringing thereby a significant reduction in CPU time: (i) ¯ A I ,I ′ = ¯ A I,I ′ ; (ii) simplifying the definitions in Eq. (6) as

$$\bar { A } _ { I ^ { \prime }, I } = 2 \text{Re} \{ \langle \Phi _ { I ^ { \prime } } | \Phi _ { I } \rangle \} \, \quad \quad ( 8 \text{a} )$$

$$\dot { \bar { C } } _ { I } = 2 R e \{ \langle \Phi _ { I } | \Phi _ { \text{aux} } \rangle \} \,,$$

where | Φ I 〉 ≡ (⊗ K /lscript =1 σ α /lscript i /lscript ) | Ψ ( ) a t 〉 and | Φ aux 〉 ≡ H aux | Ψ ( ) a t 〉 .

The ansatz (3) is fairly general and its implementation in a quantum simulator can be expected to be complex. In what follows we show that short-range few-body interactions suffice to generate an effectively adiabatic dynamics. For the sake of illustration, we start discussing the case in which ˜ H [ K , R ] aux is restricted to two-body interactions and show that it already suffices to reduce the density of excitations (DoE) by orders of magnitude with respect to the H CD = H 0 case. After discussing it, we shall include higher-order multiple-body interactions.

As an illustrative case, and without loss of generality, we shall henceforth focus on the driving of the ground state of the initial system Hamiltonian, | Ψ(0) 〉 = | ε 0 (0) 〉 . The quality of the adiabatic tracking scheme can be estimated by means of the time-dependent fidelity

$$\mathcal { F } ( t ) \equiv | \langle \Psi _ { \mathfrak { a } } ( t ) | \tilde { \Psi } ( t ) \rangle | ^ { 2 } \,, \quad \quad ( 9 ) \, \underset { \text{sur} } { \mu }$$

defined as the overlap between the adiabatic evolution of the starting state and the instantaneous state of the system denoted by | ˜ ( Ψ ) t 〉 under an evolution dictated by H 0 ( ) t + ˜ H [ K , R ] aux ( ) t while using the units in which /planckover2pi1 = 1 . The DoE n ex can then be computed as the probability of not ending up in the final ground state, i.e., n ex ≡ 1 -F . We consider a quench linear in time of the form B t ( ) = B 0 -vt and study the efficiency of ˜ H [ K , R ] aux in suppressing excitations as a function of the quench rate v for different system sizes following the paradigm of Kibble-Zurek mechanism [18].

Figure 2 summarizes the results with the simplest possible choice of a subclass of the proposed variational ansatz in Eq. (3) of the form

$$\tilde { \mathcal { H } } _ { \text{aux} } ^ { [ 2 ^ { ( y z ) }, \mathcal { R } _ { \max } ] } ( t ) \equiv \sum _ { i _ { 1 }, i _ { 2 } } ^ { ^ { \prime } } h _ { i _ { 1 }, i _ { 2 } } ^ { y, z } ( t ) \, \sigma _ { i _ { 1 } } ^ { y } \otimes \sigma _ { i _ { 2 } } ^ { z } \, \quad ( 1 0 ) \quad \text{$\stackrel { \sim.. } { \text{strat} } \\ \quad \text{$be u$} \\ \quad \text{cons}.$$

2

Figure 2. (Color online) (a) Suppression of DoE following the passage through the QCP induced by the variational ansatz restricted to two-body interactions, as in (10), as a function of the quench rate v in a linear quantum Ising chain for different system sizes. At fast quench rates, the efficiency of the ansatz in suppressing excitations saturates and a further suppression can be achieved by including higher-order multiple-body terms, as (b) depicts.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0524v2-figure-ff01cd23-e381-48bd-b254-9bfa15e22e30.jpg)

motivated by our observation that for K = 2 other permutation of spin components do not contribute to the suppression of excitations. We take the H CD = H 0 case as reference, where n ex increases monotonically until reaching saturation due to finite-size effects [39-41]. It is shown that a variational ansatz restricted to only twobody terms though with an interaction range extending through the whole chain, i.e., R max = N -1 , leads to a successful suppression of n ex by several orders of magnitude depending on the system size. Here the system is evolved in time from an initial value of the quench parameter B &lt; B t 0 ( c ) to a final one deep in the ferromagnetic phase where the value of n ex in the plot is collected. The results therefore suggest a high-fidelity adiabatic passage across the QCP of the model at t c = ( B 0 -1) /v . The residual excitations can be further suppressed by higherorder multiple-body interactions, as shown in Fig. 2(b).

Shortcuts to adiabatic state preparation. We demonstrate that the variational counterdiabatic ansatz can be used to achieve high-fidelity state preparation. We consider for this purpose a time-dependent modulation

F

-

1

0

10

-1

10

-2

10

-3

10

-4

10

-5

10

-6

10

-7

10

-8

10

N

= 8

,

τ

= 2

J , Bf

0

0.2

=

¯

h/J

0

.

= 0 1

0.4

0.6

+ ˜ [2 (

H

yz aux

H

+ ˜ [3

,

,

R

aux

0.8

R

max]

1

s

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0524v2-figure-71ae3e4b-6035-48c9-a218-806352913ab8.jpg)

≡

t/

τ

Figure 3. (Color online) (a) Time-evolution of the state preparation infidelity during a shortcut to the adiabatic driving of the ground state of a quantum Ising chain. The fidelity can be improved significantly by incorporating higher-body interactions at a truncated range of ˜ R = R max -4 . (b) Real-time flow of the interaction amplitudes associated with the implementation of the full-range two-body ansatz as in Eq. (10) during passage through QCP of the model denoted by s c . The color maps in the insets visualize the strength of two-body interactions h y,z i 1 ,i 2 among all pairs of spins at sites ( i 1 , i 2 ) in the chain.

of the magnetic field which enforces H aux to vanish at the beginning and end of the driving scheme, t = 0 , τ . We are thus led to the boundary conditions B (0) = B 0 , B τ ( ) = B f and ˙ B (0) = 0 , ˙ B τ ( ) = 0 which are satisfied by a polynomial quench of the form B s ( ) = B 0 + 3( B f -B 0 ) s 2 -2( B f -B 0 ) s 3 with s ≡ t/τ . The use of the latter quench is further motivated by adiabatic perturbation theory [42]. Figure 3 illustrates the results for preparation of the ground state of the transverse Ising chain in Eq. (2) under such a quench and for various choices of ˜ H [ K , R ] aux . We use the fidelity in Eq. (9) to assess the quality of the preparation procedure. With a full-range two-body interaction, the protocol leaves some room for improvement, as shown in Fig. 3(a), indicating the urge to employ higher-body terms for longer chain lengths. However, for higher-body interactions beyond K = 2 the computational complexity associated with the large number of spins permutations and thereby the dimensions of the matrices to be formed and inverted renders the numerical implementation of the variational procedure intractable. To circumvent such a practical challenge, we suggest to truncate over the range of the K -body ansatz by restricting it to manageable values of ˜ R /lessmuch R max . The latter truncation strategy is further motivated by the nearly tridiagonal structure of h y,z i 1 ,i 2 evident in the inset of Fig. 3(b) which shows remarkably that the dominant contributions to ˜ H [2 ( yz ) , R max ] aux consist of only short-range interactions. Figure 3(a) illustrates the success of the truncation strategy by demonstrating significant improvement in fidelity upon employing shortrange three-body interactions. We point out the number of distinct K -body interactions as the required experimental resources to achieve the maximal-fidelity state preparation within our scheme scales with 1 2 ( 4 K N ! N -K )! which is a polynomial in the size of the system N of the leading order O ( N K ) . The scaling derives from a simple combinatorics corresponding to the total number of possible choices of a K -tuple of spins from N ones in which internal permutation of tuples produces distinct choices due to noncommutativity of the Kronecker product. The prefactors 1 2 and 4 K , moreover, account for the mirror symmetry of the finite open chain under consideration and the multiplicity associated with various spin components of a Pauli sigma operator, i.e., { 0 , x, y, z } , respectively.

Conclusions. We have shown how to engineer an experimentally realizable counterdiabatic control Hamiltonian for the fast driving of many-body spin systems that mimic adiabatic driving. Our approach combines ED with a variational principle to determine the optimal CD scheme with a restricted set of control fields and leads to a suppression of the DoE by several orders of magnitude with respect to the uncontrolled driving dynamics. Although the identification of H aux is a computationally hard problem that poses a challenge to scalability of the method, tests on finite systems are still relevant to currently feasible experiments and the variational approximations to the CD Hamiltonians for these cases may guide future approaches towards large particle numbers. In congruence with recent results in optimal control theory [43-46], our results suggest that the practical implementation of our scheme represents an effort scaling only polynomially with the system size. Our proposal is ideally suited for digital quantum simulation as well as tailoring the nonequilibrium thermodynamics of manybody systems [47]. It further supplements previous adiabatic tracking schemes aimed at accessing highly excited states [48]. The possibility to reduce the level of non-

˜ ]

J

0

sc

(a)

B

0

H

H

H

0

0

0

)

locality in implementing the CD term promised by our approach may facilitate realization of adiabatic quantum computation [49, 50].

We acknowledge stimulating discussions with Daniel Lidar, Marek M. Rams, Alexey Gorshkov, and Xi-Wen Guan. H.S. is grateful to Aarhus University for support

- [1] E. Torrontegui, S. S. Ibáñez, S. Martínez-Garaot, M. Modugno, A. del Campo, D. Guéry-Odelin, A. Ruschhaupt, X. Chen, and J. G. Muga, Adv. At. Mol. Opt. Phys. 62 , 117 (2013).
- [2] M. Demirplak and S. A. Rice, J. Phys. Chem. A 107 , 9937 (2003).
- [3] M. Demirplak and S. A. Rice, J. Phys. Chem. B 109 , 6838 (2005).
- [4] M. Demirplak and S. A. Rice, J. Chem. Phys. 129 , 154111 (2008).
- [5] A. del Campo, M. M. Rams, and W. H. Zurek, Phys. Rev. Lett. 109 , 115703 (2012).
- [6] M. B. Hastings and X.-G. Wen, Phys. Rev. B 72 , 045141 (2005).
- [7] T. J. Osborne, Phys. Rev. A 75 , 032321 (2007).
- [8] K. Van Acoleyen, M. Mariën, and F. Verstraete, Phys. Rev. Lett. 111 , 170501 (2013).
- [9] M. Berry, J. Phys. A 42 , 365303 (2009).
- [10] M. G. Bason, M. Viteau, N. Malossi, P. Huillery, E. Arimondo, D. Ciampini, R. Fazio, V. Giovannetti, R. Mannella, and O. Morsch, Nature Phys. 8 , 147 (2012).
- [11] S. Ibáñez, X. Chen, E. Torrontegui, J. G. Muga, and A. Ruschhaupt, Phys. Rev. Lett. 109 , 100403 (2012).
- [12] J. Zhang, J. H. Shim, I. Niemeyer, T. Taniguchi, T. Teraji, H. Abe, S. Onoda, T. Yamamoto, T. Ohshima, J. Isoya, and D. Suter, Phys. Rev. Lett. 110 , 240501 (2013).
- [13] W. Rohringer, D. Fischer, F. Steiner, I. E. Mazets, J. Schmiedmayer, and M. Trupke, arXiv:1312.5948 .
- [14] C. Jarzynski, Phys. Rev. A 88 , 040101 (2013).
- [15] A. del Campo, Phys. Rev. Lett. 111 , 100502 (2013).
- [16] S. Deffner, C. Jarzynski, and A. del Campo, Phys. Rev. X 4 , 021013 (2014).
- [17] T. Opatrný and K. Mølmer, New J. Phys. 16 , 015025 (2014).
- [18] J. Dziarmaga, Adv. Phys. 59 , 1063 (2010).
- [19] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, Rev. Mod. Phys. 83 , 863 (2011).
- [20] A. del Campo and W. H. Zurek, Int. J. Mod. Phys. A 29 , 1430018 (2014).
- [21] R. Barankov and A. Polkovnikov, Phys. Rev. Lett. 101 , 076801 (2008).
- [22] D. Sen, K. Sengupta, and S. Mondal, Phys. Rev. Lett. 101 , 016806 (2008).
- [23] J. Dziarmaga and M. M. Rams, New J. Phys. 12 , 055007 (2010).
- [24] J. Dziarmaga and M. M. Rams, New J. Phys. 12 , 103002 (2010).
- [25] P. Doria, T. Calarco, and S. Montangero, Phys. Rev. Lett. 106 , 190501 (2011).

and hospitality. This work was financed by the European Social Fund and the state budget of the Czech Republic, project CZ.1.07/2.3.00/30.0041. This research is further supported by the U.S. Department of Energy through the LANL/LDRD Program and a LANL J. Robert Oppenheimer fellowship (AD).

- [26] T. Caneva, T. Calarco, and S. Montangero, Phys. Rev. A 84 , 022326 (2011).
- [27] J. D. Sau and K. Sengupta, arXiv:1402.4805 .
- [28] A. del Campo, T. Kibble, and W. Zurek, J. Phys.: Condens. Matter 25 , 404210 (2013).
- [29] M. Müller, K. Hammerer, Y. L. Zhou, C. F. Roos, and P. Zoller, New J. Phys. 13 , 085007 (2011).
- [30] J. T. Barreiro, M. Müller, P. Schindler, D. Nigg, T. Monz, M. Chwalla, M. Hennrich, C. F. Roos, P. Zoller, and R. Blatt, Nature 470 , 486 (2011).
- [31] J. Casanova, A. Mezzacapo, L. Lamata, and E. Solano, Phys. Rev. Lett. 108 , 190502 (2012).
- [32] H. Weimer, M. Muller, I. Lesanovsky, P. Zoller, and H. P. Buchler, Nature Phys. 6 , 382 (2010).
- [33] A. Friedenauer, H. Schmitz, J. T. Glueckert, D. Porras, and T. Schaetz, Nature Phys. 4 , 757 (2008).
- [34] K. Kim, M.-S. Chang, S. Korenblit, R. Islam, E. E. Edwards, J. K. Freericks, G.-D. Lin, L.-M. Duan, and C. Monroe, Nature 465 , 590 (2010).
- [35] K. Kim, S. Korenblit, R. Islam, E. E. Edwards, M.S. Chang, C. Noh, H. Carmichael, G.-D. Lin, L.-M. Duan, C. C. J. Wang, J. K. Freericks, and C. Monroe, New J. Phys. 13 , 105003 (2011).
- [36] P. Richerme, C. Senko, J. Smith, A. Lee, S. Korenblit, and C. Monroe, Phys. Rev. A 88 , 012334 (2013).
- [37] J. Simon, W. S. Bakr, R. Ma, M. E. Tai, and P. M. P. M. Greiner, Nature 472 , 307 (2011).
- [38] Z. Li, H. Zhou, C. Ju, H. Chen, W. Zheng, D. Lu, X. Rong, C. Duan, X. Peng, and J. Du, Phys. Rev. Lett. 112 , 220501 (2014).
- [39] J. Dziarmaga, Phys. Rev. Lett. 95 , 245701 (2005).
- [40] W. H. Zurek, U. Dorner, and P. Zoller, Phys. Rev. Lett. 95 , 105701 (2005).
- [41] A. Polkovnikov, Phys. Rev. B 72 , 161201 (2005).
- [42] D. A. Lidar, A. T. Rezakhani, and A. Hamma, J. Math. Phys. 50 , 102106 (2009).

[43]

S.

Kallush and

R.

Kosloff,

Phys. Rev. A

83

,

063412 (2011).

- [44] F. Lucas, F. Mintert, and A. Buchleitner, Phys. Rev. A 88 , 032306 (2013).
- [45] S. Lloyd and S. Phys. Rev. Lett. 113 , 010502 (2014).

Montangero,

- [46] N. Wu, A. Nanduri, and H. Rabitz, arXiv:1409.1622 .
- [47] R. Dorner, J. Goold, C. Cormick, M. Paternostro, and V. Vedral, Phys. Rev. Lett. 109 , 160601 (2012).
- [48] M. Moliner and P. Schmitteckert, Phys. Rev. Lett. 111 , 120602 (2013).
- [49] A. T. Rezakhani, W.-J. Kuo, A. Hamma, D. A. Lidar, and P. Zanardi, Phys. Rev. Lett. 103 , 080502 (2009).
- [50] A. T. Rezakhani, A. K. Pimachev, and D. A. Lidar, Phys. Rev. A 82 , 052305 (2010).
Hamed Saberi Tomáš Opatrný Klaus Mølmer Adolfo del Campo
Adiabatic tracking of quantum many-body dynamics
The nonadiabatic dynamics of a many-body system driven through a quantum
critical point can be controlled using counterdiabatic driving, where the
formation of excitations is suppressed by assisting the dynamics with auxiliary
multiple-body nonlocal interactions. We propose an alternative scheme which
circumvents practical challenges to realize shortcuts to adiabaticity in
mesoscopic systems by tailoring the functional form of the auxiliary
counterdiabatic interactions. A driving scheme resorting in few-body
short-range interactions is shown to generate an effectively adiabatic
dynamics.