In all simulations, we use the HTC Hamiltonian for N molecules coupled to the cavity. Each molecule is coupled to an independent phonon bath, discretized by N ν = 20 independent bath modes from the Debye spectral density (see Eq. 9). The bath parameters are sampled using the procedure outlined in Ref. 27 (or more generally, Ref. 40)