Deep Electric Field Predictions by Drift-Reduced Braginskii Theory with Plasma-Neutral Interactions Based on Experimental Images of Boundary Turbulence.

We present two-dimensional turbulent electric field calculations via physics-informed deep learning consistent with (i) drift-reduced Braginskii theory under the framework of an axisymmetric fusion plasma with purely toroidal field and (ii) experimental estimates of the fluctuating electron density and temperature on open field lines obtained from analysis of gas puff imaging of a discharge on the Alcator C-Mod tokamak. The inclusion of effects from the locally puffed atomic helium on particle and energy sources within the reduced plasma turbulence model is found to strengthen correlations between the electric field and electron pressure. The neutrals are also directly associated with broadening the distribution of turbulent field amplitudes and increasing E×B shearing rates. This demonstrates a novel approach in plasma experiments by solving for nonlinear dynamics consistent with partial differential equations and data without encoding explicit boundary nor initial conditions.

Explaining Cold-Pulse Dynamics in Tokamak Plasmas Using Local Turbulent Transport Models.

A long-standing enigma in plasma transport has been resolved by modeling of cold-pulse experiments conducted on the Alcator C-Mod tokamak. Controlled edge cooling of fusion plasmas triggers core electron heating on time scales faster than an energy confinement time, which has long been interpreted as strong evidence of nonlocal transport. This Letter shows that the steady-state profiles, the cold-pulse rise time, and disappearance at higher density as measured in these experiments are successfully captured by a recent local quasilinear turbulent transport model, demonstrating that the existence of nonlocal transport phenomena is not necessary for explaining the behavior and time scales of cold-pulse experiments in tokamak plasmas.

Observation of Efficient Lower Hybrid Current Drive at High Density in Diverted Plasmas on the Alcator C-Mod Tokamak.

Efficient lower hybrid current drive (LHCD) is demonstrated at densities up to n[over ¯]_{e}≈1.5×10^{20} m^{-3} in diverted plasmas on the Alcator C-Mod tokamak by operating at increased plasma current-and therefore reduced Greenwald density fraction. This density exceeds the nominal "LH density limit" at n[over ¯]_{e}≈1.0×10^{20} m^{-3} reported previously, above which an anomalous loss of current drive efficiency was observed. The recovery of current drive efficiency to a level consistent with engineering scalings is correlated with a reduction in density shoulders and turbulence levels in the far scrape-off layer. Concurrently, rf wave interaction with the edge and/or scrape-off-layer plasma is reduced, as indicated by a minimal broadening of the wave frequency spectrum measured at the plasma edge. These results have important implications for sustaining steady-state tokamak operation and indicate a pathway forward for implementing efficient LHCD in a reactor.

Turbulence Nonlinearities Shed Light on Geometric Asymmetry in Tokamak Confinement Transitions.

A comprehensive study of fully frequency-resolved nonlinear kinetic energy transfer has been performed for the first time in a diverted tokamak, providing new insight into the parametric dependences of edge turbulence transitions. Measurements using gas puff imaging in the turbulent L-mode state illuminate the source of the long known but as yet unexplained "favorable-unfavorable" geometric asymmetry of the power threshold for transition to the turbulence-suppressed H mode. Results from the recently discovered I mode point to a competition between zonal flow (ZF) and geodesic-acoustic modes (GAM) for turbulent energy, while showing new evidence that the I-to-H transition is still dominated by ZFs. The availability of nonlinear drive for the GAM against net heat flux through the edge corresponds very well to empirical scalings found experimentally for accessing the I mode.

Fast Low-to-High Confinement Mode Bifurcation Dynamics in a Tokamak Edge Plasma Gyrokinetic Simulation.

Transport barrier formation and its relation to sheared flows in fluids and plasmas are of fundamental interest in various natural and laboratory observations and of critical importance in achieving an economical energy production in a magnetic fusion device. Here we report the first observation of an edge transport barrier formation event in an electrostatic gyrokinetic simulation carried out in a realistic diverted tokamak edge geometry under strong forcing by a high rate of heat deposition. The results show that turbulent Reynolds-stress-driven sheared E×B flows act in concert with neoclassical orbit loss to quench turbulent transport and form a transport barrier just inside the last closed magnetic flux surface.

Observation of edge instability limiting the pedestal growth in tokamak plasmas.

With fusion device performance hinging on the edge pedestal pressure, it is imperative to experimentally understand the physical mechanism dictating the pedestal characteristics and to validate and improve pedestal predictive models. This Letter reports direct evidence of density and magnetic fluctuations showing the stiff onset of an edge instability leading to the saturation of the pedestal on the Alcator C-Mod tokamak. Edge stability analyses indicate that the pedestal is unstable to both ballooning mode and kinetic ballooning mode in agreement with observations.

Upgrade of the Lyman-alpha diagnostic system on DIII-D for main chamber edge neutral studies.

The LLAMA (Lyman Alpha Measurement Apparatus) pinhole camera diagnostic had previously been deployed on DIII-D to measure radial profiles of the Lyman-α (Ly-α) deuterium neutral line brightness across the plasma boundary in the lower chamber to infer neutral deuterium density and ionization rate profiles. This system has recently been upgraded with a new diagnostic head, named ALPACA, that also encloses two pinhole cameras and duplicates the LLAMA views in the upper chamber. Similar to LLAMA, ALPACA provides two times 20 lines of sight, viewing the plasma edge on the inboard and outboard sides with a radial resolution of ∼2.5 cm (FWHM) and an effective time resolution of ∼1 ms that allows for the investigation of inter-ELM dynamics. The extended Ly-α system provides better coverage to study neutrals in experiments with various plasma shapes utilizing both the upper and lower divertors. Furthermore, post-campaign calibration of the LLAMA diagnostic has successfully been demonstrated for the first time. This was facilitated by various upgrades to the calibration set-up and detailed measurements of the emissivity distribution of the Ly-α calibration source using a pinhole collimator. It was found that the sensitivity of the inboard LLAMA pinhole camera was reduced by a factor of 2.0 ± 0.2 over the course of six months of plasma operation in 2021. The upgraded Ly-α system, equipped with improved absolute calibration, will provide key input for neutral fueling and pedestal particle transport studies and for 2D edge transport code validation on the DIII-D tokamak.

Overview of the early campaign diagnostics for the SPARC tokamak (invited).

The SPARC tokamak is a high-field, Bt0 ∼12 T, medium-sized, R0 = 1.85 m, tokamak that is presently under construction in Devens, MA, led by Commonwealth Fusion Systems. It will be used to de-risk the high-field tokamak path to a fusion power plant and demonstrate the commercial viability of fusion energy. SPARC's first campaign plan is to achieve Qfus > 1 using an ICRF-heated, <10 MW, high current, Ip ∼ 8.5 MA, L-mode fueled by D-T gas injection, and its second campaign will investigate H-mode operations in D-D. To facilitate plasma control and scientific learning, a targeted set of ∼50 plasma diagnostics are being designed and built for operation during these campaigns. While nearly all diagnostics are based on established techniques, the pace of deployment, relative to the first plasma, and the harshness of the thermal, electromagnetic, and radiation environment are unprecedented for medium-sized tokamaks. An overview of the SPARC diagnostic set is given, providing context to further details communicated by the SPARC team in companion publications that are system-specific. The system engineering philosophy for SPARC diagnostics is outlined, and the design and engineering verification process for components inside and outside the primary vacuum boundary are described. Diagnostics are mounted directly to the vacuum vessel as well as housed within a series of eight midplane and 24 off-midplane replaceable port plugs. With limited exceptions, signal conditioning, digitization electronics and cameras as well as lasers and microwave sources are localized to a series of five Diagnostic Lab spaces, totaling ∼350 m2, located >15 m from the center of the tokamak, on the other side of a 2.4 m concrete shielding wall. A series of 31 large-scale penetrations have been included in the SPARC Tokamak Hall to facilitate integration of early campaign diagnostics and to provide upgradability.

Formation and stability of impurity "snakes" in tokamak plasmas.

New observations of the formation and dynamics of long-lived impurity-induced helical "snake" modes in tokamak plasmas have recently been carried out on Alcator C-Mod. The snakes form as an asymmetry in the impurity ion density that undergoes a seamless transition from a small helically displaced density to a large crescent-shaped helical structure inside q<1, with a regularly sawtoothing core. The observations show that the conditions for the formation and persistence of a snake cannot be explained by plasma pressure alone. Instead, many features arise naturally from nonlinear interactions in a 3D MHD model that separately evolves the plasma density and temperature.

Effects of magnetic shear on toroidal rotation in tokamak plasmas with lower hybrid current drive.

Application of lower hybrid (LH) current drive in tokamak plasmas can induce both co- and countercurrent directed changes in toroidal rotation, depending on the core q profile. For discharges with q(0) <1, rotation increments in the countercurrent direction are observed. If the LH-driven current is sufficient to suppress sawteeth and increase q(0) above unity, the core toroidal rotation change is in the cocurrent direction. This change in sign of the rotation increment is consistent with a change in sign of the residual stress (the divergence of which constitutes an intrinsic torque that drives the flow) through its dependence on magnetic shear.

Robust identification of multiple-input single-output system response for efficient pickup noise removal from tokamak diagnostics.

Electromagnetic pickup noise in the tokamak environment imposes an imminent challenge for measuring weak diagnostic photocurrents in the nA range. The diagnostic signal can be contaminated by an unknown mixture of crosstalk signals from coils powered by currents in the kA range. To address this issue, an algorithm for robust identification of linear multi-input single-output (MISO) systems has been developed. The MISO model describes the dynamic relationship between measured signals from power sources and observed signals in the diagnostic and allows for a precise subtraction of the noise component. The proposed method was tested on experimental diagnostic data from the DIII-D tokamak, and it has reduced noise by up to 20 dB in the 1-20 kHz range.

Deep modeling of plasma and neutral fluctuations from gas puff turbulence imaging.

The role of turbulence in setting boundary plasma conditions is presently a key uncertainty in projecting to fusion energy reactors. To robustly diagnose edge turbulence, we develop and demonstrate a technique to translate brightness measurements of HeI line radiation into local plasma fluctuations via a novel integrated deep learning framework that combines neutral transport physics and collisional radiative theory for the 33D - 23P transition in atomic helium with unbounded correlation constraints between the electron density and temperature. The tenets for experimental validity are reviewed, illustrating that this turbulence analysis for ionized gases is transferable to both magnetized and unmagnetized environments with arbitrary geometries. Based on fast camera data on the Alcator C-Mod tokamak, we present the first two-dimensional time-dependent experimental measurements of the turbulent electron density, electron temperature, and neutral density, revealing shadowing effects in a fusion plasma using a single spectral line.

Rotation reversal bifurcation and energy confinement saturation in tokamak Ohmic L-mode plasmas.

Direction reversals of intrinsic toroidal rotation have been observed in diverted Alcator C-Mod Ohmic L-mode plasmas following electron density ramps. For low density discharges, the core rotation is directed cocurrent, and reverses to countercurrent following an increase in the density above a certain threshold. Such reversals occur together with a decrease in density fluctuations with 2 cm(-1)≤k(θ)≤11 cm(-1) and frequencies above 70 kHz. There is a strong correlation between the reversal density and the density at which the Ohmic L-mode energy confinement changes from the linear to the saturated regime.

Edge temperature gradient as intrinsic rotation drive in Alcator C-Mod tokamak plasmas.

Intrinsic rotation has been observed in I-mode plasmas from the C-Mod tokamak, and is found to be similar to that in H mode, both in its edge origin and in the scaling with global pressure. Since both plasmas have similar edge ∇T, but completely different edge ∇n, it may be concluded that the drive of the intrinsic rotation is the edge ∇T rather than ∇P. Evidence suggests that the connection between gradients and rotation is the residual stress, and a scaling for the rotation from conversion of free energy to macroscopic flow is calculated.

Uncovering turbulent plasma dynamics via deep learning from partial observations.

One of the most intensely studied aspects of magnetic confinement fusion is edge plasma turbulence which is critical to reactor performance and operation. Drift-reduced Braginskii two-fluid theory has for decades been widely applied to model boundary plasmas with varying success. Towards better understanding edge turbulence in both theory and experiment, we demonstrate that a physics-informed deep learning framework constrained by partial differential equations can accurately learn turbulent fields consistent with the two-fluid theory from partial observations of electron pressure which is not otherwise possible using conventional equilibrium models. This technique presents a paradigm for the advanced design of plasma diagnostics and validation of magnetized plasma turbulence theories in challenging thermonuclear environments.

Feasibility study for a high-k temperature fluctuation diagnostic based on soft x-ray imaging.

A new pseudolocal tomography algorithm is developed for soft X-ray(SXR) imaging measurements of the turbulent electron temperature fluctuations (δ Te) in tokamaks and stellarators. The algorithm overcomes the constraints of limited viewing ports on the vessel wall (viewing angle) and limited number of lines of sight (LOS). This is accomplished by increasing the number of LOS locally in a region of interest. Numerical modeling demonstrates that the wavenumber spectrum of the turbulence can be reliably reconstructed, with an acceptable number of viewing angles and LOS and suitable low SNR detectors. We conclude that a SXR imaging diagnostic for measurements of turbulent δ Te using a pseudolocal reconstruction algorithm is feasible.

A 1D Lyman-alpha profile camera for plasma edge neutral studies on the DIII-D tokamak.

A one dimensional, absolutely calibrated pinhole camera system was installed on the DIII-D tokamak to measure edge Lyman-alpha (Ly-α) emission from hydrogen isotopes, which can be used to infer neutral density and ionization rate profiles. The system is composed of two cameras, each providing a toroidal fan of 20 lines of sight, viewing the plasma edge on the inboard and outboard side of DIII-D. The cameras' views lie in a horizontal plane 77 cm below the midplane. At its tangency radius, each channel provides a radial resolution of ∼2 cm full width at half maximum (FWHM) with a total coverage of 22 cm. Each camera consists of a rectangular pinhole, Ly-α reflective mirror, narrow-band Ly-α transmission filter, and a 20 channel AXUV photodetector. The combined mirror and transmission filter have a FWHM of 5 nm, centered near the Ly-α wavelength of 121.6 nm and is capable of rejecting significant, parasitic carbon-III (C-III) emission from intrinsic plasma impurities. To provide a high spatial resolution measurement in a compact footprint, the camera utilizes advanced engineering and manufacturing techniques including 3D printing, high stability mirror mounts, and a novel alignment procedure. Absolutely calibrated, spatially resolved Ly-α brightness measurements utilize a bright, isolated line with low parasitic surface reflections and enable quantitative comparison to modeling to study divertor neutral leakage, main chamber fueling, and radial particle transport.

Absolute calibration of the Lyman-α measurement apparatus at DIII-D.

The LLAMA (Lyman-Alpha Measurement Apparatus) diagnostic was recently installed on the DIII-D tokamak [Rosenthal et al., Rev. Sci. Instrum. (submitted) (2020)]. LLAMA is a pinhole camera system with a narrow band Bragg mirror, a bandpass interference filter, and an absolute extreme ultraviolet photodiode detector array, which measures the Ly-α brightness in the toroidal direction on the inboard, high field side (HFS) and outboard, low field side (LFS). This contribution presents a setup and a procedure for an absolute calibration near the Ly-α line at 121.6 nm. The LLAMA in-vacuum components are designed as a compact, transferable setup that can be mounted in an ex situ vacuum enclosure that is equipped with an absolutely calibrated Ly-α source. The spectral purity and stability of the Ly-α source are characterized using a vacuum ultraviolet spectrometer, while the Ly-α source brightness is measured by a NIST-calibrated photodiode. The non-uniform nature of the Ly-α source emission was overcome by performing a calibration procedure that scans the Ly-α source position and employs a numerical optimization to determine the emission pattern. Nominal and measured calibration factors are determined and compared, showing agreement within their uncertainties. A first conversion of the measured signal obtained from DIII-D indicates that the Ly-α brightness on the HFS and LFS is on the order of 1020 Ph sr-1 m-2 s-1. The established calibration setup and procedure will be regularly used to re-calibrate the LLAMA during DIII-D vents to monitor possible degradation of optical components and detectors.

Mineral separation protocol for accurate and precise rhenium-osmium (Re-Os) geochronology and sulphur isotope composition of individual sulphide species.

A temporal framework for mineral deposits is essential when addressing the history of their formation and conceptualizing genetic models of their origin. This knowledge is critical to understand how crust-forming processes are related to metal accumulations at specific time and conditions of Earth evolution. To this end, high-precision absolute geochronology utilising the rhenium-osmium (Re-Os) radiometric system in specific sulphide minerals is becoming a method of choice. Here, we present a procedure to obtain mineral separates of individual sulphide species that may coexist within specific mineralized horizons in ore deposits. This protocol is based on preliminary petrographic and paragenetic investigations of sulphide and gangue minerals using reflected and transmitted light microscopy. Our approach emphasizes the key role of a stepwise use of a Frantz isodynamic separator to produce mineral separates of individual sulphide species that are subsequently processed for Re-Os and sulphur isotope geochemistry.•Detailed method and its graphical illustration modified from an original procedure introduced by [1], [2].•Quality control and validation of monophasic mineral separates made by microscopic investigations and qualitative analysis of aliquots embedded in epoxy mounts.•The present method, which contributed to the successful results presented in the co-publication by Saintilan et al. (2020), demonstrates why other studies reporting Re-Os isotope data for mixtures of sulphide minerals should be considered with caution.

Temperature gradient scale length measurement: A high accuracy application of electron cyclotron emission without calibration.

Calibration is a crucial procedure in electron temperature (Te) inference from a typical electron cyclotron emission (ECE) diagnostic on tokamaks. Although the calibration provides an important multiplying factor for an individual ECE channel, the parameter ΔTe/Te is independent of any calibration. Since an ECE channel measures the cyclotron emission for a particular flux surface, a non-perturbing change in toroidal magnetic field changes the view of that channel. Hence the calibration-free parameter is a measure of Te gradient. BT-jog technique is presented here which employs the parameter and the raw ECE signals for direct measurement of electron temperature gradient scale length.