The global coupling of few-level quantum systems (``spins'') to a discrete set of bosonic modes is a key ingredient for many applications in quantum science, including large-scale entanglement generation, quantum simulation of the dynamics of long-range interacting spin models, and hybrid platforms for force and spin sensing. We present a general numerical framework for treating the out-of-equilibrium dynamics of such models based on matrix product states. Our approach applies for generic spin-boson systems: it treats any spatial and operator dependence of the two-body spin-boson coupling and places no restrictions on relative energy scales. We show that the full counting statistics of collective spin measurements and infidelity of quantum simulation due to spin-boson entanglement, both of which are difficult to obtain by other techniques, are readily calculable in our approach. We benchmark our method using a recently developed exact solution for a particular spin-boson coupling relevant to trapped ion quantum simulators. Finally, we show how decoherence can be incorporated within our framework using the method of quantum trajectories, and study the dynamics of an open-system spin-boson model with spatially nonuniform spin-boson coupling relevant for trapped atomic ion crystals in the presence of molecular ion impurities.
Using a renormalization group approach, we solve the time evolution of random Ising spin chains with generic interactions starting from initial states of arbitrary energy. As a function of the Hamiltonian parameters, the system is tuned through a dynamical transition, similar to the ground-state critical point, at which the local spin correlations establish true long-range temporal order. In the state with a dominant transverse field, a spin that starts in an up state loses its orientation with time, while in the ``ordered'' state it never does. As in ground-state quantum phase transitions, the dynamical transition has unique signatures in the entanglement properties of the system. When the system is initialized in a product state, the entanglement entropy grows as
log(t) in the two ``phases,'' while at the critical point it grows as ${\mathrm{log}}^{\ensuremath{\alpha}}(t)$, with $\ensuremath{\alpha}$ a universal number. This universal entanglement growth requires generic (``integrability breaking'') interactions to be added to the pure transverse field Ising model.
Generically a quantum system initialized in an unentangled state will, under unitary dynamics, rapidly become entangled, a process closely related to information transport and to thermalization. Disorder can suppress the growth of entanglement and result in memory of initial conditions. In non-interacting systems this arises from localization of single-particle states, the occupancy of which is fixed by the initial condition. In interacting systems similar localized conserved quantities persist, but with the added feature that they are coupled, resulting in entanglement growth which is distinct from both non-interacting localized systems and from generic ergodic systems. The Fermi-Hubbard model has two degrees of freedom per site -- charge and spin -- and disorder may be present in both of these. We study the growth of entanglement in two scenarios -- disorder in charge equal and unequal to that in spin, and determine the distinct contributions of charge and spin degrees of freedom by expanding the Hamiltonian in terms of a set of optimally localized conserved quantities with separate charge and spin character. We find that coupling between charge and spin is significantly weaker than charge-charge and spin-spin coupling. While this decoupling is present in all our results, it is only apparent when the strength of the disorder in the two sectors is different such that there is a separation between the characteristic timescales of the contributions to entanglement made by charge and by spin.