Quantum simulation of thermalization dynamics of a nonuniform Dicke model

Previous experimental realizations of Dicke model in atomic or ionic systems are based on global observables assuming uniform spin-boson coupling, while inevitable experimental nonuniformity on the one hand requires site-resolved measurement of spin states, and on the other hand provides potential quantum advantage on the simulation of multi-spin distributions. Here we report the quantum simulation of a nonuniform Dicke-like model in a two-dimensional (2D) crystal of up to 200 ions. We explicitly demonstrate the sensitivity of few-spin observables and multi-spin distributions to the spatial inhomogeneity of the model, and examine the thermalization dynamics of the nonuniform model by measuring the subsystem entropies of selected ion groups. Our work enables the study of Dicke-like models beyond the symmetric subspace, paving the way toward understanding the role of disorder in its thermalization and quantum chaos behavior.