Solving boundary time crystals via the superspin method
Boundary time crystals have been extensively studied through numerical and mean-field methods. However, a comprehensive analytical approach in terms of the Liouvillian eigenvalues, without semiclassical approximations, has yet to be developed. In this work, we analyse the Liouvillian spectrum of dissipative spin models, within a perturbative framework in the weak-dissipation limit. Introducing the superspin method, we compute the eigenvalues to first order in perturbation theory, providing a direct and transparent explanation for the emergence of the time crystal phase. We analytically demonstrate how spontaneous symmetry breaking occurs, leading to persistent oscillations. Our method can be used as a tool to unequivocally identify new dissipative collective spin models that support a boundary time crystal phase. We demonstrate this by applying the technique to four distinct, collective spin Liouvillians, which experience some dissipation. The first of these is the paradigmatic model studied in [1]. We demonstrate how the boundary time crystal phase is a general feature of models that yield a solution solely in terms of a quantity that we term the superspin. Furthermore, we analyse the stability of this phase and demonstrate how, in some cases, it is stable against increasing dissipation.