Electron Temperature Research Articles

Detection of alpha heating in JET-ILW DT plasmas by a study of the electron temperature response to ICRH modulation

In the JET DTE2 campaign a new method was successfully tested to detect the heating of bulk electrons by α-particles, using the dynamic response of the electron temperature T e to the modulation of ion cyclotron resonance heating (ICRH). A fundamental deuterium (D) ICRH scheme was applied to a tritium-rich hybrid plasma with D-neutral beam injection (NBI). The modulation of the ion temperature T i and of the ICRH accelerated deuterons leads to modulated α-heating with a large delay with respect to other modulated electron heating terms. A significant phase delay of ∼40° is measured between central T e and T i, which can only be explained by α-particle heating. Integrated modelling using different models for ICRH absorption and ICRH/NBI interaction reproduces the effect qualitatively. Best agreement with experiment is obtained with the European Transport Solver/Heating and Current Drive workflow.

Tokamak divertor plasma emulation with machine learning

Future tokamak devices that aim to create conditions relevant to power plant operations must consider strategies for mitigating damage to plasma facing components in the divertor. One of the goals of MAST-U tokamak operations is to inform these considerations by researching advanced divertor configurations that aid stable plasma detachment. Machine design, scenario planning and detachment control would all greatly benefit from tools that enable rapid calculation of scenario-relevant quantities given some input parameters. This paper presents a method for generating large, simulated scrape-off layer data sets, which was applied to generate a data set of steady-state Hermes-3 simulations of the MAST-U tokamak. A machine learning model was constructed using a Bayesian approach to hyperparameter optimisation to predict diagnosable output quantities given control-relevant input features. The resulting best-performing model, which is based on a feedforward neural network, achieves high accuracy when predicting electron temperature at the divertor target and carbon impurity radiation front position and runs in around 1 ms in inference mode. Techniques for interpreting the predictions made by the model were applied, and a high-resolution parameter scan of upstream conditions was performed to demonstrate the utility of rapidly generating accurate predictions using the emulator. This work represents a step forward in the design of machine learning-driven emulators of tokamak exhaust simulation codes in operational modes relevant to divertor detachment control and plasma scenario design.

Empirical line ratios for optical emission spectroscopy in low-temperature, low-density, magnetised Ar plasmas

Spectroscopic evaluation of electron temperature and density in low-temperature, low-density, magnetized plasmas can be difficult, but necessary in situations where chemical erosion and physical sputtering prevent the use of other diagnostics, such as Langmuir probes (LP). Further, in such cases, because of the low densities and temperatures, the vessel and environment involved, theoretical line ratios derived from Collisional-Radiative models may not be easily applicable. This is the case, for example, of low-temperature (<15 eV), low-density (<1011 cm−3), magnetised plasma used for plasma-material interaction studies where chemical erosion and physical sputtering can be significant. The aim of the present work is to define an empirical line ratio (ELR) method derived from an extensive calibration campaign with the two diagnostics, using LP measurements as a reference. The ELR method is useful to permit the use of optical emission spectroscopy independent of LP in conditions that are critical for the latter, resulting in an effective instrument for the evaluation of plasma parameters. Further, the use of two different lines of sight and the influence of the magnetic field intensity on the measurements are also discussed. The proposed ELR method is demonstrated here for pure Ar linear plasmas and is in principle applicable also to other similar cases.

Prediction of transport in the JET DTE2 discharges with TGLF and NEO models using the TGYRO transport code

The JET Deuterium-Tritium-Experiment Campaign 2 (DTE2) has demonstrated the highest-ever fusion energy production. To forecast the transport dynamics within these discharges, the TGLF and NEO models within the TGYRO transport code were employed. A critical development in this study is the new quasilinear transport model, TGLF-SAT2, specifically designed to resolve discrepancies identified in JET deuterium discharges. This model accurately describes the saturated three-dimensional (3D) fluctuation spectrum, aligning closely with a database of nonlinear CGYRO turbulence simulations, thereby enhancing the predictive accuracy of TGYRO simulations. In validating against the JET DTE2 discharges across two primary operating scenarios, TGYRO effectively predicted the temperature profiles within a broad radial window (ρ ∼ 0.2–0.85), though with minor ion temperature discrepancies near the core. However, a consistent underprediction of electron density profiles by 20% across the simulation domain was noted, indicating areas for future refinement. To achieve a self-consistent steady-state solution based on the JET DTE2 discharges, an integrated modeling workflow TGYRO-STEP within the OMFIT framework was introduced. This workflow iterates among the core transport, the pedestal pressure and the MHD equilibrium, ultimately yielding a converged solution that significantly reduces dependence on experimental boundary conditions for temperature and density profiles. The integrated simulation results show negligible differences in electron density and temperature profiles compared to standalone TGYRO modeling, while the ion temperature profile is lower due to the updated boundary condition in TGYRO-STEP. The application of the TGYRO-STEP workflow to JET DTE2 discharges serves as a crucial test to validate its robustness and highlights its limitations, providing valuable insights for its potential future application in ITER and Fusion Power Plant deuterium and tritium prediction modeling.

Effects of prepulse on hot electron emission from mesoscopic particles

In this paper we make an attempt to study the effect of a nanosecond prepulse on electron acceleration from microscopic solid particles. While standard scaling laws predict temperatures of ∼50 keV, our investigations showed an anomalous enhancement ( ∼200 keV) in electron temperatures from microscopic solids subjected to intense laser fields. Conventional isotropic expansion of the plasma target as a result of pre-pulse interaction is unable to explain the enhancement, prompting an experimental inquiry into the expansion form. Tentacle-like structures emerging from the spherical plasma are detected as an outcome of the pre-pulse interaction, with the region between the tentacles being targeted by the main laser pulse. We resort to PIC simulations to understand and predict the exact plasma condition that exists in this region. In this paper, we study a series of approximate pre-plasma conditions, particle shapes, and electron densities, examining them for their role in electron acceleration from this extended over-critical plasma target.

ICRH-related impurity source and control across experiments in H, D, T plasmas at JET-ILW

The experimental and theoretical analysis were focused on experiments conducted to assess the effect of plasma isotopes, protium (H), deuterium (D), and tritium (T) on ion cyclotron resonance heating (ICRH) related plasma wall interactions. Comparison of L-mode discharges with N = 1 3He and N = 1 H minority ICRH heating scenarios were done for different isotopes. For the selected pulses, the behaviour of high-Z, mid-Z and low-Z intrinsic impurity and radiated power behaviour was investigated based on data from VUV, visible spectroscopy, and bolometry diagnostic at Joint European Torus. It was found that for N = 1 3He scenario during radiofrequency antennas operation, core W, Ni content, Be source and the radiated power are higher for π/2 in comparison to dipole antenna phasing. Lowest core Ni, W content and radiated power is clearly observed for H plasmas in comparison to D and T, where for this ICRH scenario behaviour was similar. However, lower Be photon flux is observed for T in comparison to D plasmas. Be sputtering by He particles is responsible for such an effect. Additionally, several computer simulations were conducted using the COREDIV code. The difference in the electron temperature was due to the difference in the isotope masses. Increased temperature in the central plasma in the case of T plasmas leads to higher radiation in the central plasma in comparison to H plasmas. As a result, the power across separatrix is lower and the temperature on the divertor plate decreases with the increase of the isotope mass. At these temperatures on the divertor plate, W is not sputtered by the main plasma ions H, D and T and by He. For the N = 1 H ICRH scenario clear difference between D and T plasma was observed with higher metallic impurity content for T plasma in comparison to D. Impurity content in the plasmas is found to be sensitive to the power balance between the antenna straps. Its minimum is observed for the maximum of P cen/P tot.

Plasma properties in the vicinity of the last closed flux surface in hydrogen and helium fusion plasma discharges

The origin of the bi-Maxwellian electron energy distribution function (EEDF) observed in the scrape-off layer (SOL) of tokamak plasmas by means of Langmuir probes is still under discussion. It has been assumed that the ionization of hydrogen and deuterium neutrals by thermal electrons penetrating the SOL from the bulk plasma is the main reason for the appearance of a second Maxwellian. To validate this assumption, radial measurements of the electron temperatures and densities, or the plasma properties in helium plasmas in the GOLEM tokamak and the TJ-II stellarator were performed. The radial profiles of the low-temperature electron group densities follow the trend of the calculated radial profiles of the electron sources arising from the ionization of neutrals in both deuterium and helium plasmas in TJ-II. The difference in the radial location where the bi-Maxwellian EEDF appears can be explained by the difference in the rate coefficients for ionization of deuterium and helium. The results of probe measurements in GOLEM and the WEST tokamak divertor, at one radial location in the SOL, are compatible with the hypothesis concerning the ionization of neutral atoms and the type of the EEDF.

Noise reduction of stochastic density functional theory for metals.

Density Functional Theory (DFT) has become a cornerstone in the modeling of metals. However, accurately simulating metals, particularly under extreme conditions, presents two significant challenges. First, simulating complex metallic systems at low electron temperatures is difficult due to their highly delocalized density matrix. Second, modeling metallic warm-dense materials at very high electron temperatures is challenging because it requires the computation of a large number of partially occupied orbitals. This study demonstrates that both challenges can be effectively addressed using the latest advances in linear-scaling stochastic DFT methodologies. Despite the inherent introduction of noise into all computed properties by stochastic DFT, this research evaluates the efficacy of various noise reduction techniques under different thermal conditions. Our observations indicate that the effectiveness of noise reduction strategies varies significantly with the electron temperature. Furthermore, we provide evidence that the computational cost of stochastic DFT methods scales linearly with system size for metal systems, regardless of the electron temperature regime.

Effect of electron and ion mobility on edge biasing in tokamak plasmas

We present an improved model for the study of edge biasing in a tokamak plasma that incorporates electron and ion mobility contributions. The non-ambipolar nature of the drifts due to the electron/ion mobility terms influences the space charge separation due to edge biasing and affects plasma dynamics in the edge and SOL regions in a significant manner. In contrast to earlier studies, the present model enables simulation studies at higher biasing voltages. The inclusion of mobility enhances/decreases the effect of negative/positive biasing. The radial profiles of plasma density, electron temperature, radial electric field, and its shear for positive as well as negative biasing are investigated as a function of mobility.

Enhancing predictive capabilities in fusion burning plasmas through surrogate-based optimization in core transport solvers

This work presents the PORTALS framework (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), which leverages surrogate modeling and optimization techniques to enable the prediction of core plasma profiles and performance with nonlinear gyrokinetic simulations at significantly reduced cost, with no loss of accuracy. The efficiency of PORTALS is benchmarked against standard methods, and its full potential is demonstrated on a unique, simultaneous 5-channel (electron temperature, ion temperature, electron density, impurity density and angular rotation) prediction of steady-state profiles in a DIII-D ITER Similar Shape plasma with GPU-accelerated, nonlinear CGYRO (Candy et al 2016 J. Comput. Phys. 324 73–93). This paper also provides general guidelines for accurate performance predictions in burning plasmas and the impact of transport modeling in fusion pilot plants studies.

Experimental study on the discharge characteristics of air rotating gliding arc

Abstract In this study, the discharge characteristics of air rotating gliding arc are investigated by the synchronous measurements of digital oscilloscope and high-speed camera, and emission spectrum. The discharge evolution in one complete motion cycle exhibits “breakdown-elongation-extinction” process accompanied by the jump of arc root and back-breakdown phenomenon. The discharge evolves from the unstable breakdown mode (U-B), to the transition mode and finally to the stable gliding mode (S-G) by increasing the input voltage or decreasing the tangential and axial gas flow rates. The U-B mode at the input voltage of 120 V is featured by the large reduced electric field and high electron temperature of 1.90 eV, but the arc length and existence time are very short. The S-G mode at the input voltage of 270 V has relatively low breakdown frequency of 0.33 kHz and average breakdown current of 1.31 A, implying that the arc steadily glides and rotates along the spiral electrode. The average electron temperature is 0.64 eV in S-G mode, while the arc length and existence time are longer. The rotational and vibrational temperatures of N2 state are respectively measured to 2200 K and 4400 K in U-B mode, and in S-G mode are 2600 K and 4820 K. From the evolution of emission intensities of measured excited species, it is found that the NOγ band emission intensity generally rises from U-B mode to S-G mode since the gas temperature and arc existence time rise, indicating that S-G mode may be beneficial for the vibrationally-promoted Zeldovich reactions. This study could deepen the understanding of arc characteristics in air rotating gliding arc for selecting a suitable mode to achieve better plasma performance in practical applications.

RFX-mod2 diagnostic capability enhancements for the exploration of multi-magnetic-configurations

The RFX-mod2 device, the upgraded version of the previous RFX-mod with a modified magnetic boundary, is presently under realization and will start to be operated in 2025. Significant upgrades of the diagnostic capabilities have been proposed and are under development. These include a largely increased number of in-vessel magnetic and electrostatic sensors, a new fast reciprocating manipulator for the exploration of the edge plasma in a wide range of experimental conditions, the improved Thomson scattering and soft x-ray diagnostics system for a detailed determination of the behavior of the electron temperature profile, new dedicated systems for the space and time resolved analysis of x-ray spectra and neutron rate, a reflectometric diagnostic for real-time determination of plasma position, two diagnostics devoted to the imaging of light impurities and influxes behavior along with arrays of halo current sensors. These diagnostic upgrades will be accompanied by a significant effort to improve the control of the electron density and of the impurity influxes by means of proper treatment of plasma facing components with in-vessel fixed electrodes distributed over the first wall. The described advancements will allow a deeper understanding of physics phenomena in the wide variety of magnetic configurations, including the tokamak, the reversed-field pinch and the Ultra-low q, which can be produced in RFX-mod2 thanks to its flexibility and unique MHD control capabilities.

Predictive modeling of NSTX discharges with the updated multi-mode anomalous transport module

The objective of this study is twofold: firstly, to demonstrate the consistency between the anomalous transport results produced by updated Multi-Mode Model (MMM) version 9.0.4 and those obtained through gyrokinetic simulations; and secondly, to showcase MMM’s ability to predict electron and ion temperature profiles in low aspect ratio, high beta NSTX discharges. MMM encompasses a range of transport mechanisms driven by electron and ion temperature gradients, trapped electrons, kinetic ballooning, peeling, microtearing, and drift resistive inertial ballooning modes. These modes within MMM are being verified through corresponding gyrokinetic results. The modes that potentially contribute to ion thermal transport are stable in MMM, aligning with both experimental data and findings from linear CGYRO simulations. The isotope effects on these modes are also studied and higher mass is found to be stabilizing, consistent with the experimental trend. The electron thermal power across the flux surface is computed within MMM and compared to experimental measurements and nonlinear CGYRO simulation results. Specifically, the electron temperature gradient modes (ETGM) within MMM account for 2.0 MW of thermal power, consistent with experimental findings. It is noteworthy that the ETGM model requires approximately 5.0 ms of computation time on a standard desktop, while nonlinear CGYRO simulations necessitate 8.0 h on 8 K cores. MMM proves to be highly computationally efficient, a crucial attribute for various applications, including real-time control, tokamak scenario optimization, and uncertainty quantification of experimental data.

Inter-ELM pedestal turbulence dynamics dependence on q95 and temperature gradient

Abstract A series of dedicated experiments from DIII-D tokamak provide spatially and temporally resolved measurements of electron density and temperature, and multiscale and multichannel fluctuations over a wide range of conditions. Measurements of long wavelength density fluctuations in the type-I ELMing H-mode pedestal routinely revealed a coexistence of multiple instabilities that exhibit dramatic different dynamic behaviors as q95 and temperature gradient are varied, apparently responsible for limiting pedestal temperature profiles. Two distinct frequency bands of density fluctuations are modulated with ELM cycle with frequency above 200kHz propagating in the electron diamagnetic direction in the lab frame (electron mode) and below 200kHz propagating in the ion diamagnetic direction (ion mode). The electron mode amplitude peaks near the electron temperature gradient region and increases with q95 which seems to be correlated with the increased χe at higher q95, similar to the characteristics expected for Micro-tearing Mode (MTM). At higher q95, during the inter-ELM period, the ion mode decays at later phase of the ELM cycle. Consistently, the poloidal correlation length of the ion mode is also found to reduce which suggests the possible E×B flow shear suppression of ion mode at later phase of the ELM cycle as the Er well recovers. In contrast, the electron mode grows during the ELM cycle and reaches saturation at around 50-60% of ELM period. Linear gyrokinetic simulations find the MTMs as the most unstable mode in the pedestal electron temperature gradient region. The higher q95 and lower magnetic shear destabilize MTMs. These observations provide key insights of the underlying physics of multifield properties and rich dataset of experimental ‘fingerprints’ that enable new tests of the theoretical pedestal models and lead to developing a predictive model for pedestal formation on ITER and future burning plasma experiments.

OES diagnostics of atmospheric pressure argon plasma: Electron temperature and density assessment through visible bremsstrahlung inversion method and collisional-radiative model

This study determined the electron temperature and density in atmospheric pressure argon plasma using optical emission spectroscopy. The analysis combined continuum and line spectral data. Visible bremsstrahlung inversion was used to derive a partial electron energy probability function (EEPF) from the continuum spectrum. Subsequently, electron temperature was estimated assuming a two-temperature distribution based on the derived EEPF. Electron density was obtained by fitting a collisional-radiative (CR) model to the line spectrum, incorporating the obtained EEPF instead of assuming a Druyvesteynian EEPF. Comparative analysis revealed that the electron densities determined using the approach were approximately one order of magnitude lower than those derived from the CR model with the Druyvesteynian EEPF. However, they exhibited strong agreement with the results obtained by the CR model using a two-temperature distribution. This approach demonstrated favorable performance in reproducing both continuum and line spectra, revealing its high reliability and accuracy for atmospheric pressure argon plasma diagnosis.

Space-resolved electron density and temperature evaluation by x-ray pinhole camera method in an ECR plasma

X-ray emission characterization provides valuable insights about electron cyclotron resonance (ECR) plasmas. In principle, space-resolved spectroscopic techniques can be used to reveal spatial distributions of electron density and temperature. In the PANDORA (Plasma for Astrophysics, Nuclear Decay Observation, and Radiation for Archaeometry) project framework, and within the collaboration between the Atomki and INFN-LNS laboratories, we developed a high-resolution full-field x-ray pinhole setup. This setup incorporates advanced analysis techniques for single photon counted imaging in high dynamical range mode, enabling x-ray imaging and space-resolved spectroscopy at high spatial and energy resolution (560 μm and 242 eV @ 8.1 keV, respectively). Here, we introduce an innovative technique for quantitatively evaluating the local electron density and temperature of plasma, as the first application of such a method in an ECR setup. Specifically, we examine an argon plasma heated by 200 W microwave power at 14 GHz. Our analysis includes a retrospective comparison with past x-ray data collected from other ECR ion source setups. Our findings clearly reveal the formation of a plasmoid–halo structure within the plasma chamber, characterized by a dense and hot plasma almost totally enclosed inside the ECR magnetic iso-surface (the plasmoid). This plasmoid exhibits nearly uniform distribution of electron density and temperature, with only gentle gradients of both the parameters toward its edges. Inside the halo, x-ray emission is minimal or even negligible. Notably, cusp structures correspond to magnetic branches where deconfined electrons impinge upon the plasma chamber walls and endplates. The average values of temperature and density measured inside the plasmoid are 12.44±1.84 keV and (1.66±0.15)×1017 m−3, respectively.

Elevated subgingival temperature infers high bacterial pathogen counts in severe periodontitis.

Periodontal inflammation may be assessed by bleeding on probing and subgingival temperature. This pilot study evaluated the intrapatient relationship between subgingival temperature and selected bacterial groups/species in deep periodontal pockets with bleeding on probing. In each of eightadults, an electronic temperature probe identified three"hot" pockets with elevated subgingival temperature and three"cool" pockets with normal subgingival temperature among premolars/molars with 6‒10 mm probing depths and bleeding on probing. Microbial samples collected separately from the hot and cool periodontal pockets were cultured for selected periodontal pathogens. Hot compared to cool periodontal pockets revealed significantly higher absolute and normalized subgingival temperatures and yielded higher mean proportions of Porphyromonas gingivalis (10.2% for hot vs. 2.5% for cool, p = 0.030) and total red/orange complex periodontal pathogens (48.0% for hot vs. 24.6% for cool, p = 0.012). Hot versus cool deep periodontal pockets harbored significantly higher levels of major periodontal pathogens. Subgingival temperature measurements may potentially be useful to assess risk of periodontitis progression and the efficacy of periodontal therapy.

Revisit of electron temperature effect on stimulated Brillouin scattering in homogenous plasma

Stimulated Brillouin scattering (SBS) has complex dependence on the plasma electron temperature via the Landau damping, particle trapping, and consequent nonlinear frequency shift. It is found from our numerical simulation that the SBS reflectivity in its saturation stage tends to increase with the plasma electron temperature within a certain range, although the linear growth rate of SBS normally reduces with the increasing electron temperature. This is because the phase velocity of an ion acoustic wave (IAW) increases with the electron temperature, which tends to reduce the Landau damping of the IAWs and hence reduce ion trapping. In the nonlinear saturation stage, ion trapping will modify the ion distribution function and induce a negative frequency shift in the IAW. This nonlinear frequency shift will break the three-wave coupling, thereby causing saturation of the SBS. With further increase in the electron temperature, however, electron trapping will dominate over ion trapping, which induces a positive frequency shift in the IAW and can lead to the SBS saturation as well. As a result, the SBS reflectivity first increases and then decreases with increase in the electron temperature. At around the peak of the SBS reflectivity, the positive frequency shift of IAW induced by electron trapping roughly offsets the negative frequency shift induced by ion trapping.

Numerical study of plasmas start-up by electron cyclotron waves in NCST spherical tokamak and CN-H1 stellarator

According to the physics of tokamak start-up, this study constructs a zero-dimensional (0D) model applicable to electron cyclotron (EC) wave assisted start-up in NCST spherical torus (spherical tokamak) and CN-H1 stellarators. Using the constructed 0D model, the results obtained in this study under the same conditions are compared and validated against reference results for pure hydrogen plasma start-up in tokamak. The results are in good agreement, especially regarding electron temperature, ion temperature and plasma current. In the presence of finite Ohmic electric field in the spherical tokamak, a study on the EC wave assisted start-up of the NCST plasma at frequency of 28 GHz is conducted. The impact of the vertical magnetic field B v on EC wave assisted start-up, the relationship between EC wave injection power P inj, Ohmic electric field E, and initial hydrogen atom density n H0 are explored separately. It is found that under conditions of Ohmic electric field lower than ITER (~ 0.3 V m−1), EC wave can expand the operational space to achieve better plasma parameters. Simulating the process of 28 GHz EC wave start-up in the CN-H1 stellarator plasma, the plasma current in the zero-dimensional model is replaced with the current in the poloidal coil of the stellarator. Plasma start-up can be successfully achieved at injection powers in the hundreds of kilowatts range, resulting in electron densities on the order of 1017–1018 m–3.

Kilonova evolution: The rapid emergence of spectral features

Kilonovae (KNe) are one of the fastest types of optical transients known, cooling rapidly in the first few days following their neutron-star merger origin. We show here that KN spectral features go through rapid recombination transitions, with features due to elements in the new ionisation state emerging quickly. Due to time-delay effects of the rapidly expanding KN, a `wave' of these new features passing though the ejecta should be a detectable phenomenon. In particular, isolated line features will emerge as blueshifted absorption features first, gradually evolving into P\,Cygni features and then pure emission features. In this analysis, we present the evolution of individual exposures of the KN AT2017gfo observed with VLT/X-shooter, which together comprise X-shooter's first epoch spectrum (1.43\,days post-merger). The spectra of these `sub-epochs' show a significant evolution across the roughly one hour of observations, including a decrease in the blackbody temperature and photospheric velocity. The early cooling is even more rapid than that inferred from later photospheric epochs and suggests that a fixed power-law relation between the temperature and time does not describe the data. The cooling constrains the recombination wave, where a interpretation of the AT2017gfo $ feature predicts both a specific timing for the feature emergence and its early spectral shape, including the very weak emission component observed at about 1.43\,days. This empirically indicates a strong correspondence between the radiation temperature of the blackbody and the ejecta's electron temperature. Furthermore, this reverberation analysis suggests that temporal modelling is important for interpreting individual spectra and that higher-cadence spectral series, especially when concentrated at specific times, can provide strong constraints on KN line identifications and the ejecta physics. Given the use of such short-timescale information, we lay out improved observing strategies for future KN monitoring.