Quantum annealing must be sufficiently slow to prevent defect formation. We study this problem in the random transverse-field Ising model via the adiabatic theorem to obtain the limiting rate at which defect formation can be avoided. Our approach based on the strong-disorder renormalization group (SDRG) technique views the annealing process as a gradual aggregation of strongly coupled clusters, giving access to the elementary events of error formation. We present a scaling theory of the statistics of these events in general dimensions and confront it to an extensive numerical analysis carried out in one, two, and three dimensions by an efficient implementation of the SDRG method. This approach is also validated by exact diagonalization in one dimension by means of the free-fermion technique. Our results consistently show a logarithmically slow decrease of the error density with the annealing time, characterized by the dimension-dependent critical exponents of the corresponding infinite-disorder fixed point of the model. Furthermore, we provide a map of individual events of potential defect formation in the excitation energy - control parameter plane, which may help in the design of optimal annealing paths.
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
Eres 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.
We analyze mechanisms for universal out-of-equilibrium dynamics near criticality by exploring the effect of randomized quantum resetting (QR) under a finite-time quench across a quantum phase transition. Using the transverse-field Ising chain as a generic model and exploiting its exact solution, QR is found to cause a crossover of the scaling of the topological defect density with the time scale
τ of the quench, from Kibble-Zurek to anti-Kibble-Zurek scaling as
τ increases. This reflects a competition between non-adiabatic quench-driven excitations and QR, giving rise to local minima of the defect densities at optimal annealing times. These times and the corresponding local minima are shown to scale as universal power laws with the rate of QR. Additional results for the scaling of the mean excess energy suggest that a system driven across a quantum critical point exhibits the same scaling behavior under a linear quench with QR as with uncorrelated noise.
Fisher's [Phys. Rev. B
51, 6411 (1995)] real-space renormalization-group (RG) treatment of random transverse-field Ising spin chains is extended into the strongly ordered and strongly disordered Griffiths phases, and asymptotically exact results are obtained. In the noncritical region the asymmetry of the renormalization of the couplings and the transverse fields is related to a nonlinear quantum control parameter \ensuremath{\Delta}, which is a natural measure of the distance from the quantum critical point. \ensuremath{\Delta}, which is found to stay invariant along the RG trajectories, and has been expressed by the initial disorder distributions, stands in the singularity exponents of different physical quantities (magnetization, susceptibility, specific heat, etc.), which are exactly calculated. In this way we have observed a weak-universality scenario: the Griffiths-McCoy singularities do not depend on the form of the disorder, provided the nonlinear quantum control parameter has the same value. The exact scaling function of the magnetization with a small applied magnetic field is calculated, and the critical point magnetization singularity is determined in a simple, direct way.
We investigate defects in the two-dimensional transverse-field Ising ferromagnet on periodic
L×L lattices after quantum annealing from high to vanishing field. With exact numerical solutions for
L≤6, we observe the expected critical Kibble-Zurek (KZ) time scale
∝Lz+1/ν (with
z=1 and
1/ν≈1.59) at the quantum phase transition. We also observe KZ scaling of the ground-state fidelity at the end of the process. The excitations evolve by coarsening dynamics of confined defects, with a time scale
∝L2, and interface fluctuations of system-spanning defects, with life time
∝L3. We build on analogies with classical simulated annealing, where we characterize system-spanning defects in detail and find differences in the dynamic scales of domain walls with winding numbers
W=(1,0)/(0,1) (horizontal/vertical) and
W=(1,1) (diagonal). They decay on time scales
∝L3 (which applies also to system-spanning domains in systems with open boundaries) and
∝L3.4, respectively, when imposed in the ordered phase. As a consequence of
L3.4 exceeding the classical KZ scale
Lz+1/ν=L3.17 the probability of
W=(1,1) domains in SA scales with the KZ exponent even in the final
T=0 state. In QA, also the
W=(1,0)/(0,1) domains are controlled by the KZ time scale
L2.59. The
L3 scale can nevertheless be detected in the excited states, using a method that we develop that should also be applicable in QA experiments.
In transverse-field Ising models, disorder in the couplings gives rise to a drastic reduction of the critical energy gap and, accordingly, to an unfavorable, slower-than-algebraic scaling of the density of defects produced when the system is driven through its quantum critical point. By applying Kibble-Zurek theory and numerical calculations, we demonstrate in the one-dimensional model that the scaling of defect density with annealing time can be made algebraic by balancing the coupling disorder with suitably chosen inhomogeneous driving fields. Depending on the tail of the coupling distribution at zero, balancing can be either perfect, leading to the well-known inverse-square law of the homogeneous system, or partial, still resulting in an algebraic decrease but with a smaller, nonuniversal exponent. We also study defect production during an environment-temperature quench of the open variant of the model in which the system is slowly cooled down to its quantum critical point. According to our scaling and numerical results, uncorrelated disorder induces a logarithmic temporal decrease of the defect density, while balanced disorder leads again to an algebraic decay.