Statistical mechanics is founded on the assumption that a system can reach thermal equilibrium, regardless of the starting state. Interactions between particles facilitate thermalization, but, can interacting systems always equilibrate regardless of parameter values? The energy spectrum of a system can answer this question and reveal the nature of the underlying phases. However, most experimental techniques only indirectly probe the many-body energy spectrum. Using a chain of nine superconducting qubits, we implement a novel technique for directly resolving the energy levels of interacting photons. We benchmark this method by capturing the intricate energy spectrum predicted for 2D electrons in a magnetic field, the Hofstadter butterfly. By increasing disorder, the spatial extent of energy eigenstates at the edge of the energy band shrink, suggesting the formation of a mobility edge. At strong disorder, the energy levels cease to repel one another and their statistics approaches a Poisson distribution -the hallmark of transition from the thermal to the many-body localized phase. Our work introduces a new many-body spectroscopy technique to study quantum phases of matter.
P. Roushan, C. Neill, J. Tangpanitanon, V.M. Bastidas, A. Megrant, R. Barends, Y. Chen, Z. Chen, B. Chiaro, A. Dunsworth, A. Fowler, B. Foxen, E. Je rey, J. Kelly, E. Lucero, J. Mutus, M. Neeley, C. Quintana, D. Sank, A. Vainsencher, J. Wenner, T. White, H. Neven, D. G. Angelakis, and J. Martinis [under review]