Many-Body Thermodynamics on Quantum Computers via Partition Function Zeros*

Partition functions are ubiquitous in physics: they are important in determining the thermodynamic properties of many-body systems, and in understanding their phase transitions. As shown by Lee and Yang, analytically continuing the partition function to the complex plane allows us to obtains its zeros and thus the entire function. Moreover, the scaling and nature of these zeros can elucidate phase transitions. Here we show how to find partition function zeros on noisy intermediate-scale trapped ion quantum computers in a scalable manner, using the XXZ model as a prototype. We illustrate the transition from XY-like behavior to Ising-like behavior as a function of the anisotropy. While quantum computers cannot yet scale to the thermodynamic limit, our work provides a pathway to do so as hardware improves, allowing the determination of critical phenomena for systems that cannot be solved otherwise.