Adiabatic quantum annealing dynamics in a random Ising chain

Alexander I. Nesterov, G. P. Berman, Juan Carlos Beas Zepeda +1 more

Davide Venturelli, Vadim Smelyanskiy, Alejandro Perdomo‐Ortiz +3 more

Andrew D. King, Sei Suzuki, Jack Raymond +21 more

Sei Suzuki, Hiroki Oshiyama, Naokazu Shibata

We study quantum annealing of pure and random transverse Ising chains that are coupled to boson baths, using a novel numerical method based on the combination of the quasi-adiabatic propagation path integral (QUAPI) and the matrix product state (MPS) formalisms. We present numerical results on systems up to 64 spins and conclude that the baths disturb quantum annealing of pure and random Ising chains. The numerical method to compute the reduced density matrix is described in detail.

Óscar Promio Muñoz

Tommaso Caneva, Rosario Fazio, Giuseppe E. Santoro

We present here our study of the adiabatic quantum dynamics of a random Ising chain across its quantum critical point. The model investigated is an Ising chain in a transverse field with disorder present both in the exchange coupling and in the transverse field. The transverse field term is proportional to a function $\ensuremath{\Gamma}(t)$ which, as in the Kibble-Zurek mechanism, is linearly reduced to zero in time with a rate ${\ensuremath{\tau}}^{\ensuremath{-}1}$, $\ensuremath{\Gamma}(t)=\ensuremath{-}t∕\ensuremath{\tau}$, starting at $t=\ensuremath{-}\ensuremath{\infty}$ from the quantum disordered phase $(\ensuremath{\Gamma}=\ensuremath{\infty})$ and ending at $t=0$ in the classical ferromagnetic phase $(\ensuremath{\Gamma}=0)$. We first analyze the distribution of the gaps, occurring at the critical point ${\ensuremath{\Gamma}}_{c}=1$, which are relevant for breaking the adiabaticity of the dynamics. We then present extensive numerical simulations for the residual energy ${E}_{\mathrm{res}}$ and density of defects ${\ensuremath{\rho}}_{k}$ at the end of the annealing, as a function of the annealing inverse rate $\ensuremath{\tau}$. Both the average ${E}_{\mathrm{res}}(\ensuremath{\tau})$ and ${\ensuremath{\rho}}_{k}(\ensuremath{\tau})$ are found to behave logarithmically for large $\ensuremath{\tau}$, but with different exponents, $[{E}_{\mathrm{res}}(\ensuremath{\tau})∕L{]}_{\mathrm{av}}\ensuremath{\sim}1∕{\mathrm{ln}}^{\ensuremath{\zeta}}(\ensuremath{\tau})$ with $\ensuremath{\zeta}\ensuremath{\approx}3.4$, and $[{\ensuremath{\rho}}_{k}(\ensuremath{\tau}){]}_{\mathrm{av}}\ensuremath{\sim}1∕{\mathrm{ln}}^{2}(\ensuremath{\tau})$. We propose a mechanism for $1∕{\mathrm{ln}}^{2}\ensuremath{\tau}$ behavior of $[{\ensuremath{\rho}}_{k}{]}_{\mathrm{av}}$ based on the Landau-Zener tunneling theory and on a Fisher's-type real-space renormalization group analysis of the relevant gaps. The model proposed shows therefore a paradigmatic example of how an adiabatic quantum computation can become very slow when disorder is at play, even in the absence of any source of frustration.