arxiv: v2 [quant-ph] 17 Nov 2017

Transcription

1 Dissipation in aiabatic quantum computers: Lessons from an exactly solvable moel arxiv: v2 [quant-ph] 17 Nov 2017 Maximilian Keck, 1 Simone Montangero, 2, 3 Giuseppe E. Santoro, 4, 5 Rosario Fazio, 5, 1 an Davie Rossini 6, 1 1 NEST, Scuola Normale Superiore an Istituto Nanoscienze-CNR, I Pisa, Italy 2 Theoretische Physik, Universität es Saarlanes, D Saarbrücken, Germany 3 Institute for Complex Quantum Systems & Center for Integrate Quantum Science an Technologies, Universität Ulm, D Ulm, Germany 4 SISSA, Via Bonomea 265, I Trieste, Italy 5 ICTP, Straa Costiera 11, I Trieste, Italy 6 Dipartimento i Fisica, Università i Pisa an INFN, Largo Pontecorvo 3, I Pisa, Italy (Date: November 20, 2017) We introuce an stuy the aiabatic ynamics of free-fermion moels subject to a local Linbla bath an in the presence of a time-epenent Hamiltonian. The merit of these moels is that they can be solve exactly, an will help us to stuy the interplay between non-aiabatic transitions an issipation in many-boy quantum systems. After the aiabatic evolution, we evaluate the excess energy (average value of the Hamiltonian) as a measure of the eviation from reaching the target final groun state. We compute the excess energy in a variety of ifferent situations, where the nature of the bath an the Hamiltonian is moifie. We fin a robust evience of the fact that an optimal working time for the quantum annealing protocol emerges as a result of the competition between the non-aiabatic effects an the issipative processes. We compare these results with matrix-prouct-operator simulations of an Ising system an show that the phenomenology we foun applies also for this more realistic case. I. INTRODUCTION The recent experimental avances in the fiel of quantum technologies have rastically enhance our capability to control the quantum coherent ynamics of manyboy systems in a variety of physical systems, ranging from atomic an molecular optics, to trappe ions, an cavity/circuit quantum electroynamics. These progresses have mae real the possibility to experimentally realise quantum simulators [1] as well as the implementation of the first quantum algorithms [2, 3]. Together with the progresses in implementing quantum gates an concatenate them, i.e., by realising stanar circuit computation [4], recently aiabatic quantum computation (AQC) [5] an quantum annealing [6] have receive a tremenous boost thanks to the experiments performe with D-Wave machines [7 10]. The strategy unerlying aiabatic quantum computation [5, 6] is base on the fact that any quantum algorithm can be formulate in terms of ientifying the global minimum (groun state) of a given function (Hamiltonian) over a set of many local minima. On the experimental sie, quantum effects were seen to survive on eight [11, 12], sixteen [13] an even in more than one hunre qubits [14]. Whether these machines hol alreay at present the so calle quantum supremacy or not is still uner ebate [15, 16]. However, it is clearly very important to unerstan the actual moe of operation of these aiabatic computers, in orer to unerstan the limit of their performances an to push it forwar. A key problem in this framework is to unerstan the role of issipation an ecoherence on aiabatic quantum computers. This question amounts to unerstaning the key features that control the aiabatic evolution of a many-boy open quantum system. Let us first state the general problem. Suppose to be able to follow the quantum ynamics of an appropriate time-epenent Hamiltonian: H(t) = [1 f(t)] H in + f(t) H fin, t [t in, t fin ], (1) f(t) being a generic function of time, with f(t in ) = 0 an f(t fin ) = 1. The Hamiltonian H in sets the initial conition as its groun state, while the sought solution to the problem is entaile into the groun state of H fin. If the control time is much larger than the typical inverse gap between the groun state an the first excite state, the system will aiabatically follow its instantaneous groun state ψ 0 (t). Reaching the groun state of H fin by aiabatic evolution is the way AQC works [5, 6]. As long as the evolution is unitary, the only source of errors is ue to excitations generate by non-aiabatic effects in the ynamical evolution. The complexity of the aiabatic algorithms is reflecte in the scaling of the minimum gap with the number of qubits. In general, AQC also requires a special form of f(t) to gain a speeup as compare to classical algorithms, see Ref. [17] for a prominent example of an aiabatic Grover search an Ref. [18] for a review. Moreover, optimal controlle ramps can provie aitional speeups [19]. An alternative protocol, which is more general than Eq. (1), woul incorporate an extra (possibly non-linear) term, such as f(t)[1 f(t)]h E, with H E being properly chosen Hamiltonian. This may have beneficial effects on the minimum gap of the system, an therefore greatly enhance the AQC performance [20]. It is however clear that, especially for long annealing times, another important source of efects is relate to incoherent fluctuations inuce by finite temperature, or more in general by the unavoiable coupling of the system to some external environment. In this case the quantum state of the system will be mixe, escribe by a ensity

2 2 matrix ρ(t), satisfying a ynamical equation of the form: ρ t = i [H(t), ρ] + D[ρ]. (2) The first term on the right-han sie escribes the coherent unitary time evolution, which is rule by a manyboy time-varying Hamiltonian H(t), accoring to the quantum annealing protocol (1). The secon term on the right-han sie of Eq. (2) accounts for the coupling to the environment, an its form will epen on the nature of noise an issipation as well as the form of the coupling between the many-boy system an the external bath. The issipator is a completely positive, trace preserving map, that in general rives the system into a fixe point, that is a steay state or a steay-state manifol (if the fixe point is not unique). As such, it is inucing a ecay in the system towars the steay state. Note that the ifferential form of Eq. (2) generally requires a Markovian bath. Unerstaning the effect of issipation on AQC amounts to quantifying, in some way, how the state of the system eviates from the final groun state, because of the presence of the extra term D[ρ] in Eq. (2). Does the presence of the environment facilitate the reaching of the final groun state or is it etrimental for AQC? It is clear that there cannot be a unique answer to this question: the eviations from the unitary case may epen strongly on the form of D[ρ] in relation to the type of evolution impose by H(t). This variety of possible answers is reflecte in the wie spectrum of cases alreay consiere in the literature. Within the plethora of possible scenarios, it is however important to establish some general trens that may serve as guielines in going eeper in this formiable problem. This type of analysis was first performe [21] in the context of the Kibble-Zurek (KZ) mechanism for efect formation [22, 23]. The KZ mechanism is relate to the fact that, when crossing a gapless critical point of H(t), no matter how slow the variation of the Hamiltonian in time is, the aiabatic theorem is violate an a finite ensity of efects will be prouce. More than thirty years ago, Kibble put forwar a scaling argument aime at preicting the size an the number of such efects [22], while the mechanism yieling the correct scaling was foun later by Zurek [23], roughly iviing the ynamics in either aiabatic or impulsive, accoring to the istance from the critical point. The KZ mechanism has been teste in a variety of quantum toy moels at zero temperature, incluing orere an isorere systems, as well as for crossing isolate or extene critical regions (see, e.g., Ref. [24] for a review). When the annealing velocity is progressively increase, a crossover behaviour sets in between the KZ scaling an the generation of excitations ue to faster quenches, where the ynamics can be escribe by the unerlying classical moel [25]. The proliferation of efects ue to Lanau-Zener transitions is intimately relate to the occurrence of errors in the AQC. While, in the unitary case, the number of efects ecreases on increasing annealing time, the environment will be ominant for long annealing times. In this regime, one expects a efect formation which is almost inepenent on the annealing protocol. This picture was confirme an etaile in Ref. [21], where the scaling in the crossover between the KZ-ominate regime an the environment-ominate regime was also foun. In the presence of spatially correlate noise, aitional intermeiate regimes emerge ue to the comparison of the correlation length of the noise an the correlation length of the system [26]. The role of temperature an external noise was further consiere in the context of AQC in several papers, showing that in some cases it may be beneficial in reaching the target groun state. Work has been one on the comparison between AQC an classical approaches using thermal hopping [27], on the use of quantum iffusion, showing better performance than close-system quantum annealing [28] as well as on the crucial role of noise-inuce thermalisation in AQC that can outperform simulate annealing [29]. Further, the relation of thermally-assiste tunneling to quantum Monte Carlo has been stuie in Ref. [30]. In aition to this, the effect of noise of a thermal environment [31] an ecoherence [32] have been stuie as well. It shoul be kept in min that the aiabatic ynamics of issipative many-boy system is linke to the unerstaning of the Lanau-Zener problem of a two-level system couple to an environment, an extensive stuie problem in many ifferent areas [33 38]. The problem of escribing the aiabatic ynamics of a many-boy open quantum system is a formiable problem an approximations are necessary. It is however of funamental importance to have some non-trivial examples where the outcomes of the analysis are not hampere by any approximation. The aim of this work is to present some simple, yet non-trivial, examples where the aiabatic ynamics can be analyse in full etails. Most importantly, the phenomenology that we will extract is relate to the ynamics of certain spin systems that are very close to the relevant implementations of AQC. This means that the exactly solvable moels that we consier here may be use as a very useful benchmark to test important approximations in more complex cases. The bath we will eal with is Markovian. This means that the issipative term in Eq. (2) can be written in the Linbla form D[ρ] = n ) κ n (L n ρl n 1 2 {ρ, L nl n }, (3) where L n are suitable local Linbla operators that escribe the environment (to be efine later) an κ n are the corresponing couplings, which have to be positive for a Markovian Linbla master equation. The choice of local Linbla operators that we are going to stuy oes not lea to a thermal state in the steay state. We will well in more etails on this point later. We will consier quaratic fermionic moels whose Linbla ynamics can be worke out analytically [39

3 3 41], for ifferent types of local system-bath coupling. The ynamics of this class of moels can be stuie exactly, an it will help us to clarify several features of the interplay between non-aiabatic effects an incoherent transitions ue to the external bath. In orer to unerstan to which extent the results we fin can be applie to more realistic cases, later we will consier a spin-1/2 one-imensional Ising moel. In such case, for the incoherent spin ecay/pumping, the master equation cannot be mappe into local fermionic operators, therefore we will resort to a numerical stuy base on a matrixprouct-operator (MPO) representation of the ensity matrix [42, 43]. As we will iscuss in more etails in the rest of the paper, the overall phenomenology remains unchange, thus reinforcing the fact that the exactly solvable moels introuce here can be very useful benchmarks. We shoul remin that, in general, the thermoynamic properties of low-imensional systems can be strongly affecte by the imensionality: for example, thermal fluctuations wash out quantum fluctuations of finite-temperature systems in one imension, but not anymore in two imensions. Furthermore, methos that are very powerful in one imension, as the Joran-Wigner transformation, are not applicable in higher imensions. However, in the special case of free fermions, the system s thermoynamics is not expecte to change with the imensionality. In all the situations we have aresse, we fin a robust evience of the fact that an optimal working time for the quantum annealing protocol emerges as a result of the competition between the non-aiabatic effects an the issipative processes. There is an optimal time for which the final state is closest to the true final groun state. The scaling of such optimal time (an that of the corresponing generate efects) can be accurately preicte by assuming that the number of efects prouce uring the time evolution is a sum of the two contributions ue to non-aiabaticity an issipation/ecoherence. For larger working times it may happen that, epening on the type of system-bath coupling, an overshooting point sets in, where the ensity of the generate efects is larger than that for an infinitely slow annealing, which woul aiabatically rive the system through the instantaneous steay state. While this kin of behaviour cannot appear in the unitary scenario, where the efects prouction is monotonic non-increasing with the annealing spee, in the system-bath scenario this can emerge even for small systems, being eventually relate to the spectral structure of the Liouvillian an not necessarily to many-boy characteristics. The paper is organise as follows. In Sec. II we efine the Hamiltonian moels, the various issipation schemes, an the annealing problem uner investigation. We then stuy the eparture from the instantaneous groun state, in the presence of issipative processes which may cause ecay or ephasing, both for a translationally invariant free-fermion moel III, an for an Ising spin chain (Sec. IV). We en with a iscussion of our finings an with the concluing remarks in Sec. V. Technical etails on the calculations for the quaratic fermionic moel are provie in the Appenix. II. ADIABATIC DYNAMICS WITH LOCAL DISSIPATION: FROM ISING SYSTEMS TO FERMION CHAINS One of the simplest (an exactly solvable) moels exhibiting a quantum phase transition is the oneimensional Ising chain [44]. This is efine by the Hamiltonian H(t) = J n σ x nσ x n+1 Γ(t) n σ z n, (4) where σn(α α = x, y, z) are the spin-1/2 Pauli operators for the n-th spin of the chain, while J an Γ(t) are respectively the coupling strength between neighbouring spins an the transverse magnetic fiel. We assume perioic bounary conitions, in such a way as to preserve translational invariance of the moel. We will also work in units of = 1 an set the energy scale by fixing J = 1. It is possible to realise quantum annealing in the Ising moel by tuning the parameter Γ(t) in time accoring to a linear ramping, as for example: Γ(t) = t/τ, for t (, 0], (5) where τ is relate to the ramping spee. The choice (5) ensures that, uring the annealing proceure, the system will encounter a critical point in which the groun-state energy gap closes an efects will start to appear. Duality arguments [45] show that the phase transition occurs at Γ c = 1. The system is riven from a paramagnetic phase, where all the spins are aligne along the fiel irection z [i.e., the groun state of the initial Hamiltonian H(t in ) n σz n, since Γ(t in ) = + ], to a oublyegenerate ferromagnetic phase, where the spins are all pointing along the coupling irection x [i.e., the groun state of the final Hamiltonian H(t fin ) = n σx nσn+1, x since Γ(t fin ) = 0]. Several papers have alreay aresse the KZ scaling of efects with τ in the paraigmatic Ising moel, both for the clean system [46, 47] an for the isorere system [48, 49]. Here we are going to a the effect of the coupling to an external environment, moelle through a master equation of the form in Eq. (2). To retain analytic solvability of the full open-system problem [39, 40], we will however start from a mappe version of the Ising chain (4) into a free-fermion moel. The latter can be achieve by employing a Joran- Wigner transformation (JWT), which maps the spin operators in terms of spinless fermions: ( σn = exp iπ c mc m )c n, (6) m<n where σ ± n = 1 2 (σx n ± iσ y n), while c n (c n) enotes the fermionic annihilation (creation) operator on site n,

4 4 obeying the anticommutation relations {c m, c n } = δ m,n, {c m, c n } = 0. The resulting Hamiltonian for an Ising chain of length L is quaratic in such operators an reas: H(t) = [ ] c m A m,n (t) c n (c m B m,n c n + H.c.), (7) m,n where A, B respectively are a symmetric an an antisymmetric L L matrix whose sole non-zero elements are A n,n (t) = Γ(t), A n,n+1 = A n+1,n = J/2, B n,n+1 = B n+1,n = J/2. To enforce perioic bounary conitions, the following matrix elements are non-zero as well: A 1,L = A L,1 = ( 1) NF J/2 an B L,1 = B 1,L = ( 1) NF J/2, where ( 1) NF enotes the parity of the number of fermions N F = n c nc n, which commutes with H. The Hamiltonian (7) can be exactly iagonalize using a Bogolibuov transformation [50, 51]. In the following we will completely relax the requirement on fermion-parity epenent bounary conitions, an we simply assume that anti-perioic bounary conitions for the fermions are always enforce. This assumption is perfectly justifie for a purely coherent evolution, where the fermion parity is conserve an the initial groun state has an even number of fermions. The reason for enforcing the requirement even when consiering the open-system aiabatic ynamics is that, as explaine below, our Linbla operators change the fermion parity, an it woul be impossible to solve the problem by using a parity-epenent bounary conitions. More in etail, we will stuy below the aiabatic ynamics of H(t) uner the action of three ifferent types of memoryless local environments. We moel them in such a way that the Linblaian (3) is a sum of terms that act uniformly (κ n = κ, n) on each site n of the chain: i) L (1) n = c n pumping mechanism, (8) ii) L (2) n = c n ecaying mechanism, (9) iii) L (3) n = c nc n ephasing mechanism. (10) Note that while the fermionic ephasing environment (10) can be irectly mappe into a ephasing environment of spins, since there is irect local mapping c nc n 1 2 (σz n + 1), the same is not true for the ecay an pumping. Inee, when mapping from a spin operator σ n (σ + n ) to a fermionic operator c n (c n), the Joran- Wigner transformation (6) inclues a non-local operator (string). We shoul stress that the naming of these Linbla operators has been chosen with respect to their action on a system with a local, iagonal Hamiltonian, an not with respect to their actual effect on the system that we stuy. The choice of these specific Linbla operators is motivate by the possibility to stuy, in an essentially exact way, the competition between the unitary ynamics an issipative effects, focusing on features that o not epen qualitatively on the form of the coupling to the environment. In the next sections we will first aress analytically the fermionic moel (7) with issipation provie by Linbla terms as in (9)-(10), without making any reference to the mapping with spins. The effect of the nonlocal part of the JWT will be only iscusse in Sec. IV, where we will stuy numerically the Ising moel (4) with a spin-ecay mechanism provie by L (1a) n = σn that generalises Eq. (9) to spins. We restate that this choice of local Linbla operators oes not lea to thermalisation in the steay state, as for this purpose non-local terms woul be require. However they cover a special interest since in several experimental implementations, such as circuit-qed, col-atom settings or trappe ions, this kin of local amping is the relevant one. In orer to quantify the loss of aiabaticity uring the annealing protocol, originating both from the closure of the Hamiltonian gap an from the issipative processes, we are going to stuy the excess energy ε per site, at the en of the annealing. The excess energy at a given time t expresses the ifference between the instantaneous energy uring the annealing, E(t) = Tr [ H(t) ρ(t) ], where ρ(t) is the solution of the master equation (2) at time t, an the groun-state energy E 0 (t) = Tr [ H(t) ψ 0 (t) ψ 0 (t) ] of the instantaneous Hamiltonian system escribe by H(t): { } ε(t) = 1 L Tr[H(t) ρ(t)] ψ 0 (t) H(t) ψ 0 (t). (11) Using the aforementione Bogoliubov transformation, the secon term of Eq. (11) can be compute straightforwarly, while the first term E(t) is non trivial (see the Appenix). For the Ising spin system we will resort to a fully numerical MPO approach. We point out that a relate quantity of interest is the ensity of efects N 1 2L n 1 σx nσn+1, x which, in the case of orere chains an at the en of the annealing, is equivalent to the excess energy ε(0), apart from trivial constants. III. FREE-FERMIONIC SYSTEM We first analyse a fermionic system escribe by the Hamiltonian in Eq. (7), where anti-perioic bounary conitions are impose. A Fourier transform rastically helps in the iagonalisation of the unitary problem, since the ifferent momentum moes ecouple (see App. A). Note that for the sake of simplicity, here we only consier one-imensional systems, but our analysis of fermions can be easily extene to larger imensionalities, since a larger imension will affect calculations only by changing the Brillouin zone. Let us concentrate on the case in which each lattice site is couple to some external bath through a pumping mechanism as in Eq. (8). The master equation (2) uring the annealing protocol can be easily integrate via a straightforwar generalisation of the time-epenent Bogoliubov metho alreay employe by Dziarmaga [47], as etaile in App. B. The crucial point resies in the fact that, as for the Hamiltonian, the issipative part of the Linblaian with L (1) n = c n oes not mix the various

5 5 Excess energy ε(τ) 10 2 κ = 0 κ = κ = κ = κ = κ = κ = κ = κ = 1 εopt τopt Annealing time τ Dissipation strength κ FIG. 1. Final excess energy as a function of the annealing time, for the free-fermion moel (7) couple to an environment which inuces a pumping mechanism, as in Eq. (8): L (1) n = c n. The various ata sets enote ifferent values of the issipative coupling κ, as liste in the legen. Here we simulate the annealing protocol of Eq. (5) for chains of L = 10 3 sites. Black squares enote ata for κ = 0, which obey a power-law behaviour for τ > 1 with the KZ scaling exponent γ = 0.5 (ashe line). moes at ifferent momenta, once a Fourier transform has been employe. As a consequence, the ensity matrix at time t factorizes into ifferent contributions for the various moes: ρ(t) = k ρ k (t). (12) The relevant Hilbert space for each positive momentum k has imension 4, an thus the Liouvillian ynamics can be easily followe insie it. We recall that, for the unitary Schröinger ynamics, a further ecomposition into inepenent 2 2 problems was possible, ue to the aitional conservation of the fermionic parity (which is now violate by the issipative ecaying terms). The excess energy per site ε uring the annealing protocol is thus obtaine via a numerical integration of the linearize Liouville equations for each k moe (B2). For numerical convenience, we restricte the initial point of the annealing proceure (5) to t in = 5τ, an checke that the results are not appreciably affecte by this choice [49]. We stuie systems up to L = 10 3 sites an annealing times up to τ = 10 3 ; a fourth-orer Runge- Kutta integration proceure with time step t = 10 2 has been employe. Fig. 1 shows the behaviour of the excess energy at the en of the annealing protocol, ε(0), for various values of the issipation strength κ, as a function of the annealing time τ. In the absence of issipation (κ = 0), we recover the Kibble-Zurek (KZ) scaling [46, 47] ε(τ) 1/τ γ with γ = 1/2, (13) FIG. 2. Optimal excess energy ε opt (orange squares) an corresponing annealing time τ opt (violet iamons), as a function of the issipation strength κ. Numerical ata (symbols) are obtaine using the same parameters as in Fig. 1, an nicely follow a power-law behaviour (ashe lines) with slope 1/3 an 2/3, respectively. which can be obtaine by the knowlege of the Ising critical exponents associate to the phase transition at Γ c = 1 across which the system is riven. A finite issipation κ > 0 inuces a competition between the KZ mechanism of efect generation ue to the crossing of a gapless point (which is progressively reuce, with increasing annealing time τ), an the prouction of efects generate by the incoherent riving itself. Such competition clearly emerges in Fig. 1 as a non-monotonic behaviour, which generates an optimal working point for the annealing proceure in the presence of issipation. Let us now have a closer look at the non-monotonicity, an focus on the optimal (minimal) value ε opt reache by the excess energy, an on the corresponing annealing time τ opt. Figure 2 isplays how such quantities epen on κ. Our numerical ata nicely agree with a power-law behaviour over more than two ecaes of κ values, such that ε opt κ 1/3 an τ opt κ 2/3. Below we show that this behaviour can be easily preicte by assuming that the KZ prouction of efects is totally inepenent of that generate by the issipation. The above mentione competition is thus explaine in terms of an incoherent summation of the two (inepenent) contributions. A. Scaling of the optimal point We start from the observation that, after the annealing proceure, the final state of the close system can be easily written as a Bogoliubov state where excitations are provie by pairs of quasiparticles with equal an

6 6 opposite momenta [47]: ψ(t fin ) = k>0(α k + β k γ k γ k ) 0. (14) Here 0 inicates the Bogoliubov vacuum corresponing to the final groun state of H(0), α k an β k are complex amplitues, while the momentum k can take L/2 positive values from 0 to π (see App. A for etails). In the issipative case, we will not only have those oubly excite states 1 k, 1 k = γ k γ k 0, but also singly excite states such as 1 k = γ k 0 an 1 k = γ k 0, which represent further sources of efects. Inee, by using the Bogoliubov transformation, we can rotate the Master equation in this frame. This allows us to write own the ynamical equation for 1 k ρ k 1 k. We fin t 1 k ρ k 1 k = κ 0 ρ k 0 f(γ, k), (15) with ( ) 2 Γ cos k + [f(γ, k)] 1 1+Γ2 2Γ cos k = 1+ sin 2 (16) k for the specific choice of L n = c n. In the aiabatic regime where the KZ scaling argument hols an for small issipation, the ensity of efects is much smaller than 1, so that 0 ρ k 0 can be approximate by its initial value 1. Note that, since the ensity of efects N is written as N = k γ k γ k in the Bogoliubov basis, excitations of the form 1 k only contribute to the positive values of k, while excitations ue to coherent ynamics 1 k 1 k contribute to both, k an k. Following this, the incoherent part of ensity of efects can be estimate accoring to N inc = κ L k>0 0 t f[γ(t), k] = 1 κ τ, (17) 2 where the last equality has been obtaine after a change of variables from t to Γ(t) = t/τ, an observing that the summation over k > 0 after the integral over Γ yiels a constant factor L/2. Assuming now that the mechanisms of efect generation ue to KZ an ue to issipation are unrelate [52], we have: N N KZ + N inc = 1 2π 2 τ 1/2 + 1 κ τ. (18) 2 From this expression for the total ensity of efects, the optimal annealing time minimizing the efects prouction can be thus estimate by the conition τ N (τ) τopt = 0. A irect calculation gives ( ) 2/3 1 τ opt = 2π κ 2/3, (19) 2 with a corresponing ensity of efects N opt = N (τ opt ) = 3 ( ) 2/ π κ 1/3. (20) 2 (τ) κ κ = κ = κ = κ = κ = κ = κ = κ = Annealing time τ FIG. 3. Final excess-energy ifference as a function of τ, once rescale by κ. The various ata sets stan for ifferent values of κ, an correspon to those of Fig. 1, where the same colour coe has been use. A straight line inicating the linear scaling in the annealing time τ is shown in black. The preictions given by these equations are in nice agreement with our numerical ata shown in Fig. 2, keeping in min that ε = 2N. To further highlight the role of the issipation uring the annealing proceure, we have also analyse the excess energy at the en of the annealing, after subtracting the corresponing excess energy in the absence of issipation: (τ) = ε(κ, τ) ε(κ = 0, τ). (21) Note that, in orer to properly efine the quantity, we have manifeste in Eq. (21) the κ-epenence of ε. After rescaling such quantity as (τ) (τ)/κ, we observe a fairly goo ata collapse with τ, as plotte in Fig. 3. In aition, our ata obey a linear scaling as a function of the annealing time except for eviations inuce by bigger values of κ (rather than by longer annealing times τ) in the regime where the excess energy is nearly saturate to its maximal value (see also Fig. 1). We have checke that the behaviour of ε(τ) ε(+ ) towars saturation ecays with a power law as τ 1, which is in accorance with Ref. [53]. The observations mae above point towar a substantial inepenence of the role playe by the issipation, with respect to the KZ mechanism. The incoherent coupling to the external bath acts uniformly an irrespective of the aiabaticity conition rule by the groun-state energy gap. B. Interplay between pumping an ecay Here we stuy the interplay between pumping an ecaying mechanism an the question whether the steay

7 7 Excess energy ε(τ) 10 2 κ = 0 κ = κ = κ = κ = κ = κ = κ = κ = Annealing time τ εmax τmax εopt τopt Dissipation strength κ FIG. 4. Same plot as in Fig. 1, but for a free-fermion moel couple to an environment which inuces a ecay mechanism, as in Eq. (9): L (2) n = c n. We observe the same initial tren as for the pumping mechanism, however for longer annealing times we observe an overshooting before the saturation sets in. state of a system subject to both mechanisms is thermal or not. For this, first we focus on the annealing protocol in the presence of a uniform incoherent ecay mechanism = c n. The behaviour of the final excess energy ε(τ) as a function of the annealing time is shown in Fig. 4 for ifferent values of the issipation strength κ. As one can see from the figure, at relatively small annealing times the tren is qualitatively analogous to that obtaine for the incoherent pumping (see Fig. 1). The non-monotonic behaviour of ε(τ) reveals the presence of an optimal working point, where the number of efects is minimal. However for larger times τ we also recognize the appearance of an overshooting point, where the energy efects become larger than those reache for an infinitely slow annealing. Here as well, we have checke that the behaviour of ε(τ) ε(+ ), after such overshooting point, ecays with a power law as τ 1, an again a linear scaling with κ [53]. To better highlight the overshooting behaviour, let us recall that, contrary to the incoherent pumping mechanism, the incoherent ecay will rastically affect the completely fille groun state of the initial Hamiltonian at Γ(t in ) = +, since it woul ten to empty the system an thereby increasing the energy in the system. Consequently in the limit τ, where we can assume the system always to be in the instantaneous steay state, its energy will be E(t) = 2Γ(t) > 0, so it will approach its final value E(t fin ) = 0 from above. Since for 1 τ < we know that the ynamics approximately follows this open aiabatic ynamics, it is reasonable to expect that its instantaneous energy will follow a similar tren, in particular it will approach its final value from above as only, inuce by the Linbla operators L (2) n FIG. 5. Excess energy (orange squares) an corresponing annealing time (violet iamons) as a function of the issipation strength κ, in the presence of an incoherent ecay mechanism, both for the optimal working point (upper panel) an for the overshooting point (lower panel). Numerical ata (symbols) are obtaine using the same parameters as in Fig. 4, an agree with a power-law behaviour (ashe lines) with slopes 1/3 an 0.2/3 (upper), an a constant value as well as a slope of 1 (lower). well, an the corrections ue to finite τ will result in the observe overshooting. In Fig. 5 (upper panel) we analyse the minimum excess energy that is reache at the optimal working point, an the corresponing annealing time. Their behaviour with κ again follows a power law which is similar to the pumping case, as iscusse in Sec. III A. Note that the argument leaing to the scaling preictions for a pumping environment hols as well for a ecaying environment, only the function in Eq. (16) changes. However, this oes not influence the scaling behaviour iscusse here, but only the pre-factors. We also stress that, in the ecay case, the integral involve in this calculation strongly epens on the value 5τ use to replace the initial value of the fiel by, which is reasonable since we have seen that this environment creates efects alreay long before the quantum critical region is reache. As a consequence, the scaling behaviour of ε opt an the corresponing τ behaves in accorance with Eqs. (19)-(20). In the lower panel of Fig. 5 we have repeate a similar analysis for the maximum excess energy at the overshooting point an the corresponing annealing time, as a function of the issipation strength. We observe that such annealing time τ max scales linearly with κ, while the change of the maximum excess energy ε max is relatively small, since it varies by less than 10% over almost two orers of magnitue. For a better unerstaning of the overshooting, in Fig. 6 we show the instantaneous excess energies for ifferent annealing times uring the protocol. For very

8 8 Instantaneous excess energy ε(t) τ = 1 τ = 3 τ = 10 τ = 25 τ = 100 τ = 1000 Instantaneous excess energy ε(t) Γ C L n = c n L n = c n External fiel Γ(t) External fiel Γ(t) 1 0 FIG. 6. Instantaneous excess energy as a function of the external fiel Γ(t), for various annealing times τ an fixe issipation strength κ = 0.1 of the ecaying environment. One observes that for long annealing times (τ = 1000) the system follows the instantaneous steay state, which has an energy 2Γ(t), an at the en it saturates towar ε( ) = 1. Intermeiate times show the same tren but are not following the open aiabatic ynamics as closely, thus resulting in a higher final excess energy ε(τ) > ε( ). For very short annealing times (τ = 1), the influence of the issipation is much smaller an the final excess energy is smaller than ε( ). small annealing times we see that the instantaneous steay-state energy is far away from the actual ynamics an no overshooting takes place. For long annealing times, the excess energy increases hugely in the beginning an then follows the (open) aiabatic ynamics, while the behaviour is similar for intermeiate annealing times, but not as rastic. As a consequence, there is an intermeiate regime where the annealing time τ is big enough such that an overshooting can take place, an the final excess energy ε(τ) will be bigger than in the infinite-time limit ε( ). To unerline the ifference between the two kins of issipation (pumping/ecay), in Fig. 7 we plotte the instantaneous excess energy for the same parameters (τ = 10 3, κ = ), but ifferent type of issipation. We observe that, as state above, in the pumping case the excess energy mostly increases in the last fifth of the protocol, which is close to the quantum critical point. The ecaying scenario shows a completely ifferent behaviour, rather following aiabatically the instantaneous steay state of the system. Now, we turn our attention to the interplay between pumping an ecay an how the overshooting observe for pure ecay is influence. For this we stuy the final excess energy as a function of the ratio between pumping an ecaying, η = κ pump /κ ecay. In Fig. 8 we show the results for values of η ranging from 0 (no pumping), where we observe the biggest overshooting, to 1, where FIG. 7. Instantaneous excess energy uring the annealing protocol as a function of the external fiel Γ(t) for the two ifferent types of issipation. Γ c = 1 locates the critical point where the Ising-like quantum phase transition at zero temperature occurs. Here we fixe κ = 0.1. Excess energy ε(τ) η = 0 η = 1/8 η = 1/4 η = 3/8 η = 1/2 η = 5/8 η = 3/4 η = 7/8 η = Annealing time τ FIG. 8. Final excess energy as a function of the annealing time couple to an environment which inuces both a pumping as well as a ecaying mechanism. The various ata sets enote ifferent values of the ratio between the two η as liste in the legen. Here we simulate the annealing protocol of Eq. (5) for chains of L = 10 3 sites. the overshooting is completely isappeare. Note that the optimal working point oes not change much with varying η since, for the given parameters, the contribution by the pumping mechanism in this regime are far smaller than the one by the ecay. The scaling of the maximum value of the overshooting as a function of η is shown in Fig. 9. An explanation for the iminishing overshooting can be given when looking at the epenence of the instantaneous steay-state energy uring the proto-

9 εmax Ratio of pumping an ecay η Excess energy ε(τ) 10 2 κ = 0 κ = κ = κ = κ = κ = κ = κ = 1 κ = Annealing time τ FIG. 9. Maximum value of the excess energy at the overshooting point uring the annealing protocol, as obtaine from the ata in Fig. 8, as a function of the ratio η. col: if η is smaller than 1, it ecreases from an initial positive value to 0 linearly, such that an overshooting is possible. For η = 1, the instantaneous steay-state energy is constant equal to 0 such that an overshooting ue to aiabatic ynamics is prevente. For η > 1, the instantaneous steay-state energy approaches 0 linearly from below, again preventing an overshooting. Finally, we comment on the issue of thermalisation: A single qubit subjecte to both incoherent pumping (L (1) = c ) an ecay (L (2) = c) processes woul relax to a thermal state whose inverse temperature β is relate to the ratio of the strengths of the two Linbla operators. Since in our case the translational invariant quaratic Hamiltonian H(t) factorizes into many Hamiltonians (each one escribing a moe of pseuomomentum k) whose Hilbert spaces are, each of them, essentially two-imensional, this might raise the question if our system shows thermalisation as well. Inee, the steay state of each of these moes can be approximate by a thermal state with very high fielity (> 98%) for the complete range of physical relevant coupling strengths. However, the corresponing inverse temperature β (k B = 1) of each moe epens on k, an therefore the complete steay state is not well approximate by a thermal state of a single parameter β. C. Dephasing Up to now all the iscussion was base on a systembath coupling scheme which inuces a ecay/pumping mechanism. There is however a complementary effect of ecoherence, where the issipation can generate pure ephasing. This can be easily obtaine through iagonal Linbla terms L (3) n = c nc n (which are proportional to (τ) κ κ = κ = κ = κ = κ = κ = κ = κ = Annealing time τ FIG. 10. Final excess energy ε(τ) (upper panel) an rescale ifference (τ)/κ (lower panel) as a function of the annealing time τ, in the free-fermion moel (7) couple to a ephasing environment L (3) n = c nc n. Here we simulate the annealing protocol of Eq. (5) for chains of L = 501 sites. The other parameters are set as in Fig. 1. the onsite fermionic number operator), as in Eq. (10). As etaile in App. C, espite the translational invariance, in such case the solution to the master equation (2) cannot be trivially written in a tensor structure as that in Eq. (12). As a matter of fact, the Linblaian D[ρ] now transforms into a non-local object, where the ifferent momentum moes are now couple together. Therefore it is more suitable to solve a close set of 4L ifferential linear equations for the relevant two-point correlators [40], see Eq. (C15). By employing a fourth-orer Runge-Kutta integration proceure of those equations, with time step t = 10 2, we were able to reach annealing times up to τ = The main results of our analysis are summarize in Fig. 10, where we plot (upper panel) the excess energy ε(τ) at the en of the annealing, as a function of the

10 ephasing). Finally, we analyse how the final excess energy approaches the τ limit: while for pumping an ecay we observe a behaviour ε(τ) ε( ) τ 1, for ephasing we foun ε(τ) ε( ) exp( τ). εopt τopt Dissipation strength κ FIG. 11. Optimal excess energy (orange squares) an corresponing annealing time (violet iamons) as a function of the issipation strength. Dashe lines enote power-laws with slopes 1/3 an 2/3, respectively for ε opt an for τ opt. Data are taken from Fig. 10, an refer to the free-fermion moel with a ephasing environment. annealing time τ. Comparing these ata with those of Fig. 1, we immeiately recognize a qualitatively analogous tren as for the pumping mechanism. In particular, the non-monotonic behaviour again reveals a competing effect between the KZ mechanism an the incoherent ephasing. Quantitative ifferences are barely visible on the scale of the two figures. We observe a slight worsening of the annealing protocol, for the same value of τ, the excess energy being slightly larger than that of the previous case. As we i previously, we also analyse the excess-energy ifference (τ) rescale by κ (bottom panel). Its scaling with τ is completely analogous to that in Fig. 3, with ata growing linearly with the annealing time, an eventually eviating for sufficiently large values of κ an τ. Finally we recall that the argument of Sec. III A for etermining the scaling of the optimal working point for the annealing protocol as a function of κ hols also in this case. Inee the corresponing ata (with the same power-laws), shown in Fig. 11, are closely similar to those of Fig. 2). Summarizing the results of our analysis on the quantum annealing in a translationally invariant free-fermion moel interacting with a local environment, the emerging scenario for the ifferent types of issipation is the following. For all the three incoherent mechanisms we observe a competition which leas to the onset of an optimal working point for the annealing proceure at a given τ opt rate. On the other sie, for larger values of τ an overshooting point appears only in the presence of a ecay mechanism, ue to the fact that the instantaneous energy approaches the steay-state value ε(τ) = 1 from below (while the opposite happens for the pumping an for the IV. ISING CHAIN Let us now go back to the spin-1/2 language an iscuss the effects of the coupling to an external bath on the quantum annealing of the Ising chain, Eq. (4). We first notice that ephasing can be inuce by a Linbla term L (3a) n = σn, z which is reaily mappe into the local fermionic operator (2c nc n 1), through the JWT of Eq. (6). In such case, one woul thus recover the ephasing mechanism for free fermions (we refer to Sec. III C for etails). On the other han, incoherent pumping/ecay woul be inuce by L (1a) n = σ n + an L (2a) n = σn, respectively; in that case, when mapping into fermions, the appearance of the JW string operator forbis an analytic treatment as the one iscusse previously. Let us thus concentrate on the latter scenario. We employ a numerical metho base on an efficient approximation of the many-boy ensity matrix in terms of a MPO [42, 43]. We expect this to be vali whenever the amount of correlations in the system is sufficiently small to satisfy an area-law scaling for the bipartite entanglement in the operator space. The time evolution is performe by means of the time-evolving block ecimation (TEBD) algorithm, after a Trotter ecomposition of the Liouvillian superoperator in the right-han sie of Eq. (2). In our simulations of the annealing protocol (1) for the Ising moel we consiere systems up to L = 20 sites, using MPOs with a bon link m 250 an aopting a typical Trotter step t = We aopte the same time epenence of the fiel Γ(t) as in Eq. (5), where for practical convenience we starte from t in = 3τ, an verifie that (on the scales of the figures shown below) the results are not affecte by this choice. As etaile below, we foun an emerging physical scenario which is consistent to that previously iscusse in Sec. III, alreay for small sizes L 10. Our numerical results for the annealing in the presence of incoherent ecay, showing the final excess energy as a function of τ an for various issipation strengths κ, are summarize in Fig. 12. Despite the KZ power-law scaling cannot be seen for limite system sizes (not even in the absence of issipation), the non-monotonicity of the various curves for κ 0 clearly emerges as a result of the open-system ynamics. We ascribe this behaviour to the emerging picture escribe in Sec. III, where we iscusse much longer systems of free fermions. Inee, in Fig. 13 we repeate the same analysis for the scaling of the optimal working time τ opt an of the corresponing optimal excess energy ε opt with the issipation strength, fining a similar power-law behaviour. The exponents o agree within 20% of relative ifference. We point out that we were not able to fully resolve the overshooting

11 11 V. CONCLUSIONS Excess energy ε(τ) 10 2 κ = 0 κ = κ = κ = κ = Annealing time τ FIG. 12. Final excess energy as a function of the annealing time, for the Ising chain (4) couple to an environment through Linbla operators inucing a ecay mechanism L (2a) n = σn. The various ata sets enote ifferent values of the issipative coupling κ, as liste in the legen. Fille symbols an continuous lines refer to chains with L = 20 sites, while empty symbols with ashe lines are for L = 10. εopt τopt Dissipation strength κ FIG. 13. Optimal excess energy (orange squares) an corresponing annealing time (violet iamons) as a function of the issipation strength. Numerical ata (symbols) are obtaine using the same parameters as in Fig. 12, for L = 10, an fairly agree with power laws (ashe lines) of slopes 1/3 an 2/3, respectively. behaviour in this case, since it woul require longer annealing times. However this is alreay visible in Fig. 12, for the curve corresponing to κ = 0.1. Moreover, we also checke that the scaling with τ of the final excess-energy ifference (τ) is again linear for sufficiently small values of κ an τ, as for the fermionic moel. We presente an extensive stuy of the aiabatic ynamics of free-fermion moels, being riven across their quantum critical point, within an open-system approach using local Linbla operators. Using the excess energy, we quantifie the eviations from aiabatic ynamics of the groun state, showing a competition between the unitary ynamics following a KZ mechanism an incoherent efect generation ue to issipation. While being local, the stuie environment covers a wie range of possible sources of issipation varying from ecay an pumping, separately or simultaneously, to ephasing an at the same time showing a consistent behaviour for all of them: The competition between the two processes leaing to an optimal working point. This can be moelle by the Ansatz of inepenent processes, which brings to a scaling behaviour that preicts the observe optimal working point in a fairly accurate way. For larger annealing times, we highlighte the possibility to observe an overshooting point, where efects become larger than those reache for an infinitely slow annealing. This effect is intrinsically ue to the coupling with an external bath, which rives the system towar the steay state accoring to the Liouvillian ynamics of the master equation. Furthermore, we stuie the one imensional Ising chain, which is closely relate to the free-fermion moels, by means of a matrix-prouct-operator technique, where we foun the same behaviour for small system sizes as well, suggesting a generic nature of the observe phenomena. Within the framework of free-fermion moels, a generalization to higher imensions is straightforwar. ACKNOWLEDGMENTS We thank S. Barbarino an D. Venturelli for fruitful iscussions. We acknowlege the CINECA awar uner the ISCRA initiative, for the availability of high performance computing resources an support. SM gratefully acknowleges the support of the DFG via a Heisenberg fellowship. This research was partly supporte by the EU projects QUIC an RYSQ, by SFB/TRR21 an by CRP-QSYNC. Appenix A: Unitary ynamics Here we provie technical etails concerning the ynamics of the free-fermion moel in Eq. (7), where perioic bounary conitions are impose. Namely, H(t) = L n=1 { (c n c n+1 + c nc n+1 + H.c.) + 2Γ(t)c nc n }, (A1) with c L+1 = c 1 for the positive parity sector, while c L+1 = c 1 for the negative parity sector. Here we have implicitly set the coupling strength to one.

12 12 The annealing proceure of Eq. (5) in this context has been alreay stuie in Ref. [47]. The approach consists in employing a Fourier transform of the type c n = iπ/4 c k e ikn, L k (A2) where the operators c ( ) k satisfy canonical anticommutation relations for fermions as well, an the inex k takes values (assuming L to be even, without loss of generality) k = ±1 π L, ±3 π L,..., ±(L 1) π L. The resulting Hamiltonian in Fourier space takes the form H = { 2c k c [ ] } k Γ(t) cos k + sin k(c k c k + c kc k ). k (A3) Since H conserves the fermionic parity, the global Hilbert space H can be written as a irect sum over ifferent k subspaces: H = k>0 H k, an inee H = k>0 H k as in Eq. (A3). Each subspace at fixe k > 0 is built from the two states { 0, 1 k, 1 k } or { 1 k, 1 k }, epening on the parity number (even or o, respectively). The groun state is foun in the even parity sector, as one can see from the iagonalization of H k. The Hamiltonian H k at a given fixe time t, as extrapolate from Eq. (A3), can be reaily iagonalize by means of a Bogoliubov transformation [50, 51] c k = u k (t)γ k + v k(t)γ k, (A4) so that the groun state is annihilate by all the quasiparticle operators γ k. The system ynamics once the parameter Γ(t) is varie can thus be foun by employing the time-epenent Bogoliubov metho [54], which makes the Ansatz that the instantaneous system wave function ψ(t) is annihilate by a set of quasiparticle operators γ k(h), in the Heisenberg representation, which are efine though the transformation c k(h) = u k (t) γ k(h) + v k (t) γ k(h). This Ansatz satisfies the Heisenberg equation t c k(h) = i[h(t), c k(h) ], (A5) with the constraint γ k(h) ψ(t) = 0, provie the coefficients u k (t) an v k (t) obey the time-epenent Bogoliubov-e Gennes equations i t u k = 2u k [Γ(t) + cos k] + 2v k sin k, i t v k = +2v k [Γ(t) + cos k] + 2u k sin k. (A6) These equations can be integrate starting from the initial conition Γ( ) = +, an thus mapping them into a Lanau-Zener problem [47]. Appenix B: Fermionic ecay bath Let us now escribe the superimpose action of an environment that inuces ecay in the system, so the Linbla operator on each site n is given by L n = c n. It is important to stress that, when applying the Fourier transform (A2) on the issipative part of the master equation (3), this oes not mix ifferent moes: D[ρ] ( ) D k [ρ] = κ c k ρ c k 1 2 {ρ, c k c k}. (B1) k k The reason resies in the fact that each term contains two fermionic operators c ( ) n, such as for example, n c n ρ c n n 1 L k,k c k ρ c k e i(k k )n, an thus the exponential factor, once summe over n, gives a Kronecker elta δ k,k. As a consequence, the ensity matrix factorizes into ρ(t) = k>0 ρ k (t), an we can ecouple the problem into the same k moes as for the non-issipative case. Notice however that the issipation part violates parity conservation of fermions, therefore here the ifferent subspaces for a given k > 0 are built up from the four states { 0, 1 k, 1 k, 1 k, 1 k }, an not simply from two states. The Hamiltonian in this basis can be explicitly written as H = k>0 H k, where sin k 0 2(Γ + cos k) 0 0 H k = 0 0 2(Γ + cos k) 0, 2 sin k 0 0 4(Γ + cos k) with c k = an c k = As a matter of fact, solving the full quantum ynamics of ρ(t) translates into solving L/2 Linbla equations of imension 4 for ρ k (t). In the vectorize form, they can be written as the following linear ifferential equations with imension 16 (k > 0): t ρ k = κ 2 {i ( 1 H k H k 1 ) [ 1 ( c k c k + c k c ) ( k + c k c k + c k c ) ] k 1 +κ ( c k c k + c k c k ) } ρ k. (B2) In practice, for every linear operator W = m,n W mn m n acting on the four-imensional Hilbert space H k which is spanne by the basis { m } m=1,...,4, we associate a vector in the 16-imensional superoperator space H k H k, which is spanne by the basis { m n } m,n=1,...,4, using the convention W m,n W φ, φ = m + (n 1) L, (B3) with W m,n W mn m n. In this way we have that W 1 W 2 W 3 = (W 1 W3 T ) W 2, where T enotes the transpose operation. A vectorization of the master equation for ρ k (t) using this rule, the fact that H k = Hk T = H k, an that c k = ct k, we finally arrive at Eq. (B2). The excess energy (11) is then reaily obtaine, since ρ(t) = k>0 ρ k (t).