Fast digital methods for adiabatic state preparation
We present a quantum algorithm for adiabatic state preparation on a gate-based quantum computer, with complexity polylogarithmic in the inverse error. Our algorithm digitally simulates the adiabatic evolution between two self-adjoint operators H0 and H1, exponentially suppressing the diabatic error by harnessing the theoretical concept of quasi-adiabatic continuation as an algorithmic tool. Given an upper bound α on ∥H0∥ and ∥H1∥ along with the promise that the kth eigenstate ∣ψk(s)⟩ of H(s)≡(1−s)H0+sH1 is separated from the rest of the spectrum by a gap of at least γ>0 for all s∈[0,1], this algorithm implements an operator U such that ∥∣ψk(1)⟩−U∣ψk(s)⟩∥≤ε using O(α2/γ2)polylog(α/γε) queries to block-encodings of H0 and H1. In addition, we develop an algorithm that is applicable only to ground states and requires multiple queries to an oracle that prepares ∣ψ0(0)⟩, but has slightly better scaling in all parameters. We also show that the costs of both algorithms can be further reduced under certain reasonable conditions, such as when ∥H1−H0∥ is small compared to α, or when more information about the gap of H(s) is available. For certain problems, the scaling can even be improved to linear in ∥H1−H0∥/γ up to polylogarithmic factors.