Thermalization of a Lipkin-Meshkov-Glick model coupled to a bosonic bath

We derive a Lindblad master equation that approximates the dynamics of a Lipkin-Meshkov-Glick (LMG) model weakly coupled to a bosonic bath. By studying the time evolution of operators under the adjoint master equation we prove that, for large system sizes, these operators attain their thermal equilibrium expectation values in the long-time limit, and we calculate the rate at which these values are approached. Integrability of the LMG model prevents thermalization in the absence of a bath, and our work provides an explicit proof that the bath indeed restores thermalization. Imposing thermalization on this otherwise non-thermalizing model outlines an avenue towards probing the unconventional thermodynamic properties predicted to occur in ultracold-atom-based realizations of the LMG model.