where σ ε is the width of the Gaussian centered at ε. This introduces static inhomogeneity within the molecular excitation energies. The sampling is performed independently for each trajectory. Following the previous work, we also study the effect of orientation disorder by sampling the molecular orientation (with respect to the cavity mode) from Gaussian distribution,