Preparing thermal states of frustrated quantum spin systems using 139 qubits

Finite-temperature properties of strongly correlated quantum matter are central to condensed matter, chemistry, and high-energy physics, yet are often inaccessible to classical methods such as quantum Monte Carlo (QMC). Here, we investigate dissipative thermal state preparation of frustrated spin systems using digital quantum computers. We focus on two paradigmatic models on the kagome lattice: the antiferromagnetic Heisenberg model (AFHM), whose finite-temperature properties are inaccessible to QMC due to a severe sign problem, and the antiferromagnetic Ising model (AFIM), which serves as a sign-problem-free benchmark. Using IBM quantum processors, we prepare approximate thermal states of the AFIM on kagome lattices with up to $79$ spins coupled to $60$ environment qubits. We observe the emergence of a robust steady state with an adjustable effective temperature that persists in circuits with over 1000 layers of two-qubit gates. We further study the scalability of the dissipative protocol through classical statevector simulations of the AFIM and AFHM. On lattices with up to 24 sites, we find that the circuit depth to reach thermal equilibrium is independent of system size and grows at most linearly with inverse temperature. These results establish engineered dissipation as a promising approach to finite-temperature quantum simulation of frustrated matter, and point toward regimes where quantum devices may outperform classical methods.