Quantum simulation of the Dicke model in a two-dimensional ion crystal

Quantum many-body systems driven far from equilibrium can exhibit phase transitions, chaos and entanglement, but probing these phenomena in large, nearly closed systems remains challenging. The Dicke model, which couples an ensemble of two-level systems collectively to a single bosonic mode, provides a fundamental framework for exploring these phenomena in light-matter interactions. Here we show that a two-dimensional crystal of approximately 100 trapped ions realizes this model and accesses integrable dynamical phase transitions, signatures of chaotic dynamics and fluctuation-driven spin-phonon excitation growth. Adiabatic elimination of the vibrational mode reveals ferromagnetic and paramagnetic dynamical spin phases, whereas strong spin-phonon coupling produces irregular trajectories and local relaxation consistent with quantum thermalization. Following a quench from an unstable fixed point, correlated excitation growth precedes generalized vacuum Rabi collapses and revivals; calibrated calculations identify two-mode squeezing 2.6 dB below the standard quantum limit, supported by reconstruction from spin measurements. By connecting collective dynamical instabilities to experimentally accessible signatures of quantum correlations, these results establish large ion crystals as a scalable platform for investigating non-equilibrium light-matter dynamics and generating resources for quantum-enhanced sensing.