Figure . Population dynamics of the upper polaritonic state (UP), lower polaritonic state (LP), and sum over all dark states (DS) for the HTC model computed with different MQC dynamics methods, including Ehrenfest, GFSH, γ-SQC, and spin-LSC. Different numbers of molecules (Na = 5/10/15) are used in the dynamics simulations. The cavity frequency is ωc = 2.0 eV, the light-matter coupling strength is gc = 68.1 meV, the reorganization energy is λ = 30 meV, the temperature is β = 1000 a.u. and the cutoff frequency during the initial sampling process is ων =24.8 meV. See Model 1 in Table .