SEARCH · Engineering Papers

Results for “electron density profile”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Oxygen and hydrogen ion densities above Millstone Hill.

Measurements of the vertical flux of oxygen ions, when combined with simultaneous measurements of electron density and of electron and ion temperatures, present a unique opportunity to examine conditions in the topside ionosphere. The measured fluxes and densities may be used directly to evaluate terms in the O+ continuity equation without requiring any assumptions to be made about neutral winds, electric fields, or the ambipolar diffusion coefficient. From observations at Millstone Hill we have, in this way, derived the rate of loss of O+ in the charge-exchange reaction with hydrogen, the rate of photoionization of atomic oxygen, and the rate of loss of O+ in reactions with N2 and O2. In combination with laboratory and theoretical results, these rates can be interpreted to yield number densities of the corresponding neutral species. When the measured fluxes are examined for consistency with measured density gradients and temperature profiles, we find clear evidence of wind-induced or electrodynamic vertical drifts larger than 60 m/sec.

Schunk, R. W.↗

DOE OSTI2020

MHD stability and disruptions in the SPARC tokamak

SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Inviscid damping of an elliptical vortex subject to an external strain flow

Inviscid spatial Landau damping is studied experimentally for the case of oscillatory motion of a two-dimensional vortex about its elliptical equilibrium in the presence of an applied strain flow. Here, the experiments are performed using electron plasmas in a Penning–Malmberg trap. They exploit the isomorphism between the two-dimensional Euler equations for an ideal fluid and the drift-Poisson equations for the plasma, where plasma density is the analog of vorticity. Perturbed elliptical vortex states are created using E x B strain flows, which are generated by applying voltages to electrodes surrounding the plasma. Measurements of spatial Landau damping (also called critical-layer damping) are in agreement with previous studies in the absence of an applied strain, where the damping is due to a resonance between the local fluid motion and the vortex oscillations. Interestingly, the damping rate does not change significantly over a wide range of applied strain rates. This can be accurately predicted from the initial vorticity profile, even though the resonant frequency is reduced substantially due to the applied strain. For higher amplitude perturbations, nonlinear trapping oscillations also exhibit behavior similar to the strain-free case. In principle, higher-order effects of the applied strain, such as separatrix crossing of peripheral vorticity and interactions with harmonics of the fundamental resonance, are expected to change the damping rate. However, this occurs only for conditions that are not realized in the experiments described here. Vortex-in-cell simulations are used to investigate the possible roles of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Isotopic mass effects of tritium-fueled high-performance TFTR supershots

An increase in total stored energy correlated with the addition of tritium fuel was observed in supershots during the TFTR DT campaign. This supershot regime had strikingly high, centrally peaked ion and electron temperatures, and the largest neutron emission rates observed in TFTR. This paper presents a study of the causes of this increase in stored energy in supershots. Twenty-six supershots have been recently reanalyzed with the TRANSP plasma analysis code. Early TRANSP simulations did not accurately match the measured magnitude and time evolutions of the neutron emission rates. This mismatch is attributed to neglecting apparent increases of trace amounts of heavy impurities during neutral beam injection. The new TRANSP runs were tuned to accommodate this and match the measured global neutron emission rates. These new runs also had improved fidelity in predicting the time histories and radial dependencies of measured DT neutron emission rates. That in turn adds confidence in the simulated thermal deuterium and tritium density profiles that are needed for calculating the average hydrogenic atomic mass profiles. Six subsets of these supershots had well matched toroidal field B tor , plasma current I p , flux geometry, and total injected neutral beam power. The mix of D and T beam ions was varied for different discharges. The magnitude of the increase of the thermal ion energy W i with added tritium was relatively small, and the total W tot increase is dominated by the increase in fast beam ions with T. Analyses at times before the occurrences of deleterious MHD instabilities yielded scaling of W tot with the volume-average isotopic mass consistent with previous publications. The relative fraction of fast energy ions is expected to be small in practical tokamak reactors. Thus the increase in stored energy W tot observed in TFTR supershots does not appear likely to be significantly helpful for producing useful fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1991

Physics of heavy ions (1989-1990)

The results from studies on polar wind ion heating due to kinetic ion beam instabilities and the effects of such ion heating on the outflow of O(+) in the polar wind are presented and discussed. First, the linear instabilities associated with an O(+) and H(+) polar wind plasma in the presence of O(+) and H(+) beams for a range of O(+)/H(+) beam densities, T(sub e)/T(sub i), and ion beam speeds were examined. Then, nonlinear heating of the polar wind ions was studied, using numerical simulations. The O(+) and H(+) polar wind ions were modeled by isotropic Maxwellian distributions, and the electrons, O(+) beams, and H(+) beams were modeled by drifting Maxwellian distributions. The effects of the kinetic ion heating on the outflow of the polar wind ions were examined from the ionosphere, using a time-dependent hydrodynamic model. A numerical code to solve the O(+) and H(+) continuity and momentum equations in a flux tube from ionospheric to magnetospheric altitudes were developed. The effects of ion heating were included by allowing for the altitudinal variation of the ion temperatures in the momentum equation. The ion temperature profiles were specified based on the ion heating characteristics found from previous kinetic simulations. It was assumed that heating occurred above 1500 km and increased to a saturated value of temperature that was obtained directly from the kinetic simulation study. The characteristics of the dynamical polar wind without ion heating were studied, and a flux tube on closed field lines that suddenly became open at t = 0 was simulated. Then, the effects of ion heating were included. To gain some physical insight, two limiting cases were considered: preferential H(+) heating and preferential O(+) heating. How O(+) heating can lead to enhanced polar wind O(+) fluxes in the polar magnetosphere is shown.

Ashour-Abdalla, M.↗

DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗

DOE OSTI2023

Towards fast, accurate predictions of RF simulations via data-driven modeling: Forward and lateral models

Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. More accurate simulations with fast electron diffusion are even slower, requiring multiple hours of run time with parallel processing. The machine learning models use a database of 16,000+ GEN-RAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods implemented in πScope ensure that the database covers the range of 9 input parameters (n e0 , T e0 , I p , B t , R 0 , n ∥︀ , Z e f f , V loop , P LHCD ) with sufficient density in all regions of parameter space. The surrogate models reduce the computation time from minutes-hours to ms with high accuracy across the input parameter space. Data-driven surrogate models also allow for solving inverse and “lateral” problems. A surrogate model for the inverse problem maps from a desired current drive or power deposition profile to a set of input parameters that would result in such a profile, while a surrogate model for the lateral problem maps from a measured experimental quantity such as hard x-ray emission to a current drive or power deposition profile. In conclusion, the πScope database creation workflow is flexible and applicable to other RF simulation codes such as TORIC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

High-performance plasmas after pellet injections in Wendelstein 7-X

A significant improvement of plasma parameters in the optimized stellarator W7-X is found after injections of frozen hydrogen pellets. The ion temperature in the post-pellet phase exceeds 3 keV with 5 MW of electron heating and the global energy confinement time surpasses the empirical ISS04-scaling. The plasma parameters realized in such experiments are significantly above those in comparable gas-fuelled discharges. Here, we present details of these pellet experiments and discuss the main plasma properties during the enhanced confinement phases. Local power balance is applied to show that the heat transport in post-pellet phases is close to the neoclassical level for the ion channel and is about a factor of two above that level for the combined losses. In comparable gas-fuelled discharges, the heat transport is by about ten times larger than the neoclassical level, and thus is largely anomalous. It is further observed that the improvement in the transport is related to the peaked density profiles that lead to a stabilization of the ion-scale turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

Stellar model chromospheres. XIII - M dwarf stars

Single-component, homogeneous model chromospheres that are consistent with high-resolution profiles of the Ca II K line calibrated in surface flux units for three dMe and 2 dM stars observed at quiescent times are constructed. The models reveal several systematic trends. Large values of the ratio of T(min) to T(eff) are derived, indicating a large amount of nonradiative heating present in the upper photospheres of M dwarf stars. It is also found that the lower chromospheric temperature gradient is similar for all the M dwarf stars. Since for the models here the chromospheric K line emission strength is most sensitive to the total amount of chromospheric material present within the approximate temperature range T(min)-6000 K, increasing the emission strength is not simply due to increasing chromospheric temperature gradients. It is also found that both the electron density and electron temperature at one thermalization length in the K line below the top of the chromospheres are greater in the dMe stars than in the dM stars. The M dwarf models here have microturbulent velocities between 1 and 2 km/sec, which are much smaller than for solar chromosphere models.

Giampapa, M. S.↗

NASA NTRS2000

2D Quantum Mechanical Study of Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25, 50 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. Surprisingly, the self-consistent potential profile shows lower injection barrier in the channel in quantum case. These results are qualitatively consistent with ID Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗

DOE OSTI2023

Role of Simple Spatial Gradient in Reinforcing the Accuracy of Temperature Determination of HED Plasma via Spectral Line-Area Ratios

We report on the simulation of temperature gradients in tamped NaFMgO target-foil plasma, heated and backlit by z-pinch dynamic hohlraum radiation. Our approach compares the spectroscopic output of a collisional-radiative model (prismspect) with soft X-ray absorption spectra collected on Sandia National Laboratories’ (SNL) Z Pulsed Power Facility. The pattern of minimum χ2 is seen to agree with an efficient, three-parameter model. Results show that a negligible gradient in electron temperature Te is consistent with experimental data, justifying the assumptions of previous work. The predicted sensitivity of line spectra to the gradient-aligned profile of Te is documented for each spectral feature, so that the line-area ratio between a pair of spectral features may be assessed as a proxy for the existence and quantification of such gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Isotope mass scaling and transport comparison between JET Deuterium and Tritium L-mode plasmas

The dimensionless isotope mass scaling experiment between pure Deuterium and pure Tritium plasmas with matched $\rho^*$, $\nu^*$ , β n , q and $T_e/T_i$ has been achieved in JET L-mode with dominant electron heating (NBI+ohmic) conditions. 28% higher scaled energy confinement time $B_t\tau_{E,\mathrm{th}}/A$ is found in favour of the Tritium plasma. This can be cast in the form of the dimensionless energy confinement scaling law as $\Omega_i \tau_{E,\mathrm{th}} \sim A^{0.48 \pm 0.16}$. This significant isotope mass scaling is consequently seen in the scaled one-fluid heat diffusion coefficient $A \chi_{\mathrm{eff}}/B_t$ which is around 50% lower in the Tritium plasma throughout the whole plasma radius. The isotope mass dependence in the particle transport channel is negligible, supported also by the perturbative particle transport analysis with gas puff modulation. The comparison of the edge particle fuelling or ionisation profiles from the EDGE2D-EIRENE simulations show that the absolute density differences that are necessary for the dimensionless match in the confined plasma dominate over any isotope mass dependencies of particle fuelling and ionization profiles at the plasma edge. Local GENE simulation results indicate a mild anti-gyroBohm effect at $\rho_\mathrm{tor}$ = 0.6 and thereby a small isotope mass dependence in favour of Tritium on heat transport and a negligible effect on particle transport. A significant fraction of the isotope scaling and reduced heat transport observed in the Tritium plasma is not captured in the GENE and ASTRA-TGLF-SAT2 simulations by simply changing the isotope mass for the same input profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗

DOE OSTI2021

Ion temperature anisotropy in the tokamak scrape-off layer

Understanding tokamak exhaust-power heat loads on divertor plates depends critically on having a realistic model of the scrape-off layer (SOL) plasma. The Braginskii fluid model is often solved to understand the SOL plasma behavior. This model is based on the collisional limit for transport along the magnetic field B → . The ions and electron gyrofrequencies are assumed to be much larger than the Coulomb collision frequencies, which are nonetheless, sufficiently large to yield common parallel and perpendicular temperatures for each species, i.e. the temperatures are assumed to be isotropic. Additionally, in certain circumstances such as encountered for the tokamak H-mode, the ion temperature can be quite anisotropic. In this work, the anisotropy effects are implemented in the two-dimensional (2D) transport code UEDGE. Various geometries (1D slab, 2D slab and a toroidal tokamak geometry) are used to study the 2D structure of ion temperature anisotropy and its effects on plasma transport in detail. Results show that the effects of ion temperature anisotropy on the plasma parallel transport are substantial near the magnetic X-point, which leads to different steady state density profiles in the divertor regions. The extra mirror force introduced by ion temperature anisotropy can be one of the main forces contributing to the plasma flow in the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Modeling and Characterization of Surface Discharges in Insulating Material for Spacers: Electrode Shape, Discharge Mode, and Revision of the Creepage Concept

In the design of MV AC and DC spacers, the predominant factors are surface and interface conditions. Design is generally carried out on specifications and standards which are based on long-term experience and lab testing. However, the diffusion of power electronics with a trend to increase electric field, switching frequency, and rise time to achieve higher power density calls for an innovative, global approach to optimized insulation system design. A new methodology, based on field simulation, discharge modeling, and partial discharge inception measurements, called the three-leg approach, can form the basis to optimize insulation design for any type of supply voltage waveform. This paper focuses on the influence of the type of electrode on the inception and phenomenology of surface discharges and, as a consequence, on the interpretation of the results used for application of the three-leg approach. It is demonstrated that a typical electrode system used for insulating material testing can generate both gas and surface discharges at the triple point, when the electrodes have a smooth profile that is used to avoid corona or flashover. Hence, testing partial discharge may not provide a straightforward indication of the surface discharge inception and, thus, be partially misleading for insulation design. Another takeover is that such analysis must benefit from PD testing tools endowed with analytics able to provide automatic identification of the type of defect generating PD, i.e., internal, surface, and corona, since design and remedy actions can be taken, and adequate insulating materials developed, only knowing the type of source generating PD. Hence, testing partial discharge may not provide a straightforward indication of surface discharge inception and, thus, be partially misleading for insulation design. In addition to the importance of the three-leg approach to favor reliable and optimized design of insulation systems, there is a clear need to have a PD testing tool endowed with analytics. It should preferably be able to provide automatic identification of the type of defect generating PD, i.e., internal, surface, and corona.

36 MATERIALS SCIENCE↗

NASA NTRS1994

The downward flux of O(+) over the nightside of Venus

We have constructed a map of the downward flux of O(+) over the nightside of Venus at high and low solar activities through a combination of modeling and analysis of Pioneer Venus ion mass spectrometer data. O(+) density profiles were obtained for almost 40 inbound or outbound segments of orbits from the first 2 years of the mission and about 40 more in the recent reentry phase of the mission. We have determined the nearly linear relationship between the 0(+) maximum density and the downward O(+) flux for several solar zenith angles and local times by constructing models of the nightside ionosphere of Venus for a range of downward ion fluxes at the upper boundaries. We find that the largest downward fluxes occur near the terminators, and the fluxes fall off sharply toward the antisolar point. Although the standard deviations in the data are large, there is a suggestion of a local maximum near 155 deg, and the location of this maximum correlates fairly well with structure in the peak electron density as a function of solar zenith angle reported for the Pioneer Venus radio occulation experiment. The average downward ion flux is inferred to be about 1.7 x 10(exp 8)/sq cm/sec over the nightside hemisphere at solar maximum and a factor of 7 less at moderately low solar activity.

Brannon, J. J.↗

DOE OSTI2024

Quasilinear modelling of collisional trapped electron modes

The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$. One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.

fusion plasma↗

DOE OSTI2024

Analyzing the effect of doping concentration in split-well resonant-phonon terahertz quantum cascade lasers

The effect of doping concentration on the temperature performance of the novel split-well resonant-phonon (SWRP) terahertz quantum-cascade laser (THz QCL) scheme supporting a clean 4-level system design was analyzed using non-equilibrium Green’s functions (NEGF) calculations. Experimental research showed that increasing the doping concentration in these designs led to better results compared to the split-well direct-phonon (SWDP) design, which has a larger overlap between its active laser states and the doping profile. However, further improvement in the temperature performance was expected, which led us to assume there was an increased gain and line broadening when increasing the doping concentration despite the reduced overlap between the doped region and the active laser states. Through simulations based on NEGF calculations we were able to study the contribution of the different scattering mechanisms on the performance of these devices. We concluded that the main mechanism affecting the lasers’ temperature performance is electron-electron (e-e) scattering, which largely contributes to gain and line broadening. Interestingly, this scattering mechanism is independent of the doping location, making efforts to reduce overlap between the doped region and the active laser states less effective. Optimization of the e-e scattering thus could be reached only by fine tuning of the doping density in the devices. By uncovering the subtle relationship between doping density and e-e scattering strength, our study not only provides a comprehensive understanding of the underlying physics but also offers a strategic pathway for overcoming current limitations. This work is significant not only for its implications on specific devices but also for its potential to drive advancements in the entire THz QCL field, demonstrating the crucial role of e-e scattering in limiting temperature performance and providing essential knowledge for pushing THz QCLs to new temperature heights.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗