We present a two-quantum well THz intersubband laser operating up to 192 K. The structure has been optimized with a non-equilibrium Green's function model. The result of this optimization was confirmed experimentally by growing, processing and measuring a number of proposed designs. At high temperature (T > 200 K), the simulations indicate that lasing fails due to a combination of electron-electron scattering, thermal backfilling, and, most importantly, re-absorption coming from broadened states.Terahertz quantum cascade lasers (QCLs) 1 are interesting candidates for a wide variety of potential applications 2,3 . However, to date, their operation is limited to ∼200 K 4 and the necessity of cryogenic cooling hinders a widespread use of these devices. In the last decade, significant scientific effort has been directed towards identifying the main temperaturedegrading mechanisms 5-8 , as well as finding optimized QCL designs 9-14 . The degrading mechanisms include thermal backfilling 3,15 , thermally activated LO phonon emission [6][7][8]16 , increased broadening [17][18][19] , and carrier leakage into continuum states 20 . When numerically optimizing a design, it is important to take all of these effects into consideration, in order to ensure a close correspondence between the model and the real device. Combined with the fact that the optimization parameters are typically trade-offs for one another, the task is very complex. Here, typically simpler rate equation or density matrix models are used in order to more quickly sweep the parameter space 21-23 , while more advanced models, such as non-equilibrium Green's functions (NEGF) or Monte-Carlo, are used to validate and analyze the final designs 13,[24][25][26] . In contrast, in this work we will employ an advanced model directly at the optimization stage. Specifically, we shall use a NEGF model 27 , capable of accurately simulating experimental devices 13,26,28 and including the most general treatment of scattering, from all relevant processes.The goal of the optimization is to achieve the highest possible operating temperature. Thus, the gain of the active medium should be maximized at high lattice temperature, and simultaneously the external losses minimized. The key figures for gain are inversion, oscillator strength, and line width 29 . These are mainly controlled by the doping density, the energy difference E ex between the lower laser level ll and the extractor state e, and the width of the two barriers: the laser and injection barriers. Population inversion increases with doping, although a too high level promotes detrimental effects, such as electron-electron scattering. E ex , which is chosen to be close to the LO phonon resonance E LO in order to have a short ll lifetime, and the laser frequency ω are mainly determined by the well widths. The laser barrier width determines the oscillator strength, which at the same time affects inversion; a more vertical transition