Matthew Hagan, Nathan Wiebe (Feb 06 2025).

Abstract: In this work we investigate a model of thermalization wherein a single ancillary qubit randomly interacts with the system to be thermalized. This not only sheds light on the emergence of Gibbs states in nature, but also provides a routine for preparing arbitrary thermal states on a digital quantum computer. For desired $\beta$ and random interaction $G$ the routine boils down to time independent Hamiltonian simulation and is represented by the channel $\Phi : \rho \mapsto \mathbb{E}_G {\rm Tr}_{\rm Env} \left[ e^{-i(H + \alpha G)t} \left(\rho \otimes \frac{e^{-\beta H_E}}{\mathcal{Z}}\right) e^{i (H + \alpha G)t} \right]$. We rigorously prove that these dynamics reduce to a Markov chain process in the weak-coupling regime with the thermal state as the approximate fixed point. We upper bound the total simulation time required in terms of the Markov chain spectral gap $\lambda_\star$, which we compute exactly in the ground state limit. These results are independent of any eigenvalue knowledge of the system, but we are further able to show that with knowledge of eigenvalue differences $\lambda_S(i) - \lambda_S(j)$, then the total simulation time is dramatically reduced. The ratio of the complete ignorance simulation cost to the perfect knowledge simulation cost scales as $\widetilde{O} \left({\frac{\|{H_S}\|^7}{\delta_{\rm min}^7 \epsilon^{3.5} \lambda_\star(\beta)^{3.5}}}\right)$, where $\delta_{\min}$ is related to the eigenvalue differences of the system. Additionally, we provide more specific results for single qubit and harmonic oscillator systems as well as numeric experiments with hydrogen chains. In addition to the algorithmic merits, these results can be viewed as broad extensions of the Repeated Interactions model to generic Hamiltonians with unknown interactions, giving a complete picture of the thermalization process for quantum systems.

Arxiv: https://arxiv.org/abs/2502.03410