Coupled carrier–phonon nonequilibrium dynamics in terahertz quantum cascade lasers: a Monte Carlo analysis