In the presence of strong magnetic fields the electronic bandstructure of graphene drastically changes. The Dirac cone collapses into discrete non-equidistant Landau levels, which can be externally tuned by changing the magnetic field. In contrast to conventional materials, specific Landau levels are selectively addressable using circularly polarized light. Exploiting these unique properties, we propose the design of a tunable laser operating in the technologically promising terahertz spectral range. To uncover the many-particle physics behind the emission of light, we perform a fully quantum mechanical investigation of the non-equilibrium dynamics of electrons, phonons, and photons in optically pumped Landau-quantized graphene embedded into an optical cavity. The gained microscopic insights allow us to predict optimal experimental conditions to realize a technologically promising terahertz laser.The terahertz (THz) regime of the electromagnetic spectrum can be exploited in a wide range of applications including medical diagnostics, atmosphere and space science as well as security and information technology [1][2][3][4]. Although THz research has progressed significantly in the last 20 years, the transition from laboratory demonstration to practical environment has occurred slowly and only for some niche applications. The largest challenge is the lack of adequate, tunable THz radiation sources. In 1986, H. Aoki proposed to design Landau level (LL) lasers exploiting the discreteness of LLs in two-dimensional electron gases [5]. Here, the energetic LL spacing and thus the possible laser frequency can be externally tuned through the magnetic field. The challenge for the realization of such a laser is to obtain a stable population inversion, i.e. a larger carrier occupation within an energetically higher LL. Since conventional semiconductors exhibit an equidistant spectrum of LLs, strong Coulomb scattering acts in favor of an equilibrium Fermi-Dirac distribution and strongly counteracts the build-up of a population inversion. In contrast, graphene as a two-dimensional material with a linear dispersion exhibits a non-equidistant LL separation offering entirely different conditions for many-particle processes [6][7][8]. Exploiting these remarkable properties of Landauquantized graphene, we propose an experimentally accessible scenario to achieve continuous wave lasing with tunable frequencies in the technologically promising terahertz spectral regime. The non-equidistant arrangement of energy levels [9,10] combined with selection rules for circularly polarized light, allow to selectively address individual inter-LL * samuel.brem@physik.tu-berlin.de Figure 1. Laser scheme. Sketch of the energetically lowest Landau levels in graphene embedded in an optical cavity. (A) The linearly polarized optical pump field induces transitions l = −3 → +2 and l = −2 → +3 (yellow arrows), which results in a population inversion between l = +1 and l = +2 in the conduction band and l = −1 and l = −2 in the valance band.(B) The emission o...