Many-body phases in Jaynes-Cummings-Hubbard arrays

Disorder in one-dimensional (1D) many-body systems facilitates abundant phases such as many-body localization (MBL) and thermalization. However, it remains unclear regarding their existence and behavior within hybrid quantum systems. Here, based on a simple bosonic-spin hybrid model, known as the Jaynes-Cummings-Hubbard (JCH) array, we investigate the effect of disorder compared to the phenomena in the clean system with the variation of atom-photon coupling strength. By using the level-spacing ratio, entanglement entropy, and the properties of observable diagonal and off-diagonal matrix elements, we find that strong disorder results in the appearance of a MBL phase in the JCH model that strongly violates the eigenstate thermalization hypothesis (ETH), while a conditional prethermal behavior can exist in the weak disorder regime. The conditional prethermal dynamics is based on the choice of initial product states. This work systematically reveals abundant many-body phases in the 1D JCH model and clarifies the discrepancies in the thermalization properties of systems with and without disorder.