The design of an optical diagnostic system to localize the radiation front in a mid-leg pumped divertor configuration on DIII-D is presented. Divertor detachment is a key mechanism for handling power exhaust in tokamaks, and mid-leg pumping offers a promising approach to radiatively dissipate power while maintaining acceptable core performance. To predict the location of the radiation front and design a spectroscopic diagnostic to measure its position, a database of SOLPS-ITER simulations across a range of input powers and gas puffing rates representative of DIII-D operation was employed. These simulations provide self-consistent plasma backgrounds for a Cherab–Raysect synthetic diagnostic framework, which incorporates detailed tokamak geometry and physically accurate ray-tracing. Synthetic simulations of optical sightlines and viewing cone geometry were used to calculate line-integrated emission from the plasma, particularly of the C III 465 nm line, to serve as a proxy for the Te≈7–10 eV temperature region associated with the onset and evolution of divertor detachment. Neutral deuterium emission is also evaluated for comparison. The synthetic diagnostic is used to assess the performance of existing DIII-D optical systems, including filterscopes and the Multichord Divertor Spectrometer, and to optimize line-of-sight placement within mechanical and installation constraints. The results provide quantitative guidance for diagnostic implementation on DIII-D and demonstrate the advantage of integrated synthetic diagnostics for divertor design studies and future advanced divertor concepts.

A series of L-mode plasma discharges was performed in the DIII-D tokamak to assess the impact of outer strike point (OSP) position and toroidal magnetic field direction on erosion and core contamination potential of the recently-installed, tungsten-coated Small Angle Slot (SAS-VW) divertor. In one discharge, in-slot emission spectroscopy measured an < 48 % increase in the W gross erosion rate when the OSP was moved 3 cm outwards, away from the V-shaped vertex of the slot divertor. However, the effective W yield (erosion rate divided by the incident D flux) was, overall, insensitive to changes in OSP location. Consistently low estimates of the effective W yield based on measurements taken a few cm outwards from the vertex suggest potentially significant C surface contamination. No W emission signal was detected when orienting the toroidal magnetic field such that the ion Bx del B drift direction is pointed away from the X-point. However, measurements of W content in the plasma core for both toroidal magnetic field directions suggest the presence of additional, unmeasured sources of erosion. The difference in the measured core W density with OSP position is much greater than the difference in the measured erosion rates, which may suggest that the leakage of eroded impurities out of the divertor is governed primarily through the parallel ion temperature gradient and friction forces.

At the high ne proposed for high-field fusion reactors, it is uncertain whether ionization, as opposed to plasma transport, will be most influential in determining ne at the pedestal and separatrix. A database of Alcator C-Mod discharges is analyzed to evaluate the impact of source modification via cryopumping. The database contains similarly-shaped H-modes at fixed IP=0.8 MA and Bt=5.4 T, spanning a large range in Pnet and ionization. Measurements from an edge Thomson scattering system are combined with those from a midplane-viewing Lyα camera to evaluate changes to ne and Te in response to changes to ionization rates, Sion. nesep and Teped are found to be most sensitive to changes to Sionsep, as opposed to neped and Tesep. Dimensionless quantities, namely αMHD and ν∗, are found to regulate attainable pedestal values. Select discharges at different values of Pnet and in different pumping configurations are analyzed further using SOLPS-ITER. It is determined that changes to plasma transport coefficients are required to self-consistently model both plasma and neutral edge dynamics. Pumping is found to modify the poloidal distribution of atomic neutral density, n0, along the separatrix, increasing n0 at the active X-point. Opaqueness to neutrals from high ne in the divertor is found to play a role in mediating neutral penetration lengths and hence, the poloidal distribution of neutrals along the separatrix. Pumped discharges thus require a larger particle diffusion coefficient than that inferred purely from 1D experimental profiles at the outer midplane.

Experimental measurements of plasma and neutral profiles across the pedestal are used in conjunction with 2D edge modeling to examine pedestal stiffness in Alcator C-Mod H-mode plasmas. Enhanced D alpha experiments on Alcator C-Mod observed pedestal degradation and loss in confinement below a critical value of net power crossing the separatrix, Pnet=Pnetcrit approximate to 2.3 MW, in the absence of any external fueling. New analysis of ionization and particle flux profiles reveal saturation of the pedestal electron density, neped, despite continuous increases in ionization throughout the pedestal, inversely related to Pnet. A limit to the pedestal del ne emerges as the particle flux, Gamma D, continues to grow, implying increases in the effective particle diffusivity, Deff. This is well-correlated with the separatrix collisionality, nu sep & lowast; and a turbulence control parameter, alpha t, implying a possible transition in type of turbulence. The transition is well correlated with the experimentally observed value of Pnetcrit. SOLPS-ITER modeling is performed for select discharges from the power scan, constrained with experimental electron and neutral densities, measured at the outer midplane. The modeling confirms general growth in Deff, consistent with experimental findings, and additionally suggests even larger growth in chi e at the same Pnetcrit.

This paper presents experimental measurements of the location of the impurity flow stagnation point in the scrape-off-layer (SOL) of a tokamak plasma. Coherence imaging of carbon-2+ emission (465 nm) is used to track the main-chamber impurity velocity of DIII-D L-mode plasmas with B×∇B out of the divertor. The C2+ flow stagnates near the top or crown of the plasma when an open divertor (no baffling) is used. In contrast, with matched conditions and using a divertor with baffling, the C2+ flow stagnates near the outer divertor leg (X-point). The C2+ poloidal emission is hollow, peaking near the divertor legs, in the open configuration. In contrast, in the closed configuration, the C2+ emission is flat through most of the main-chamber SOL. Changing divertor dissipation from attached to detached conditions had only a minor effect on the main-chamber midplane impurity velocity. Numerical simulations using the multi-fluid edge transport code UEDGE including cross-field drifts show qualitative agreement with the open divertor experimental result.

Current profile reconstructions are obtained for high current ( I p similar or equal to 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.

Experiments performed during strongly-shaped high-power diverted negative triangularity (NT) experiments in DIII-D achieved detached divertor conditions and a transient-free edge, showcasing the potential for application of NT to a core-edge integrated reactor-like scenario and providing the first characterization of the parametric dependencies for detachment onset. Detached divertor conditions will be required in future devices to mitigate divertor heat fluxes. Access to dissipative divertor conditions was investigated via an increase in upstream density. Detachment onset at the outer strike point was achieved with H-mode level confinement H98-y2 similar to 1 and reactor-relevant normalized pressures beta N similar to 2 . Confinement degradation was observed with deeper detachment, associated with the loss of an electron temperature pedestal. Differences in geometry, radial transport, impact of cross field drifts are discussed to explain differences in access to detachment in NT discharges. Higher normalized densities, with respect to equivalent discharges in positive triangularity, were necessary to achieve detachment, partially explained by the shorter parallel connection length to the targets. The effect of cross-field particle drifts (ExB, B x del B) on access to detachment was demonstrated by the lower upstream density needed to access detachment with ion B x del B drift directed outside of the active divertor (Greenwald fraction fGw similar to 0.9-1.0 vs fGw similar to 1.3). The upstream density at detachment onset was observed to increase linearly with plasma current with ion B x del B drift into the divertor, consistent with the observed narrowing of the scrape-off layer heat flux width lambda q . Edge fluid simulations capture separatrix densities needed to achieve detachment in NT plasma and their dependence on drift direction. The ability to reproduce detachment dynamics in NT plasma increases the confidence in future design studies for NT divertors.

In ITER, pellets are calculated to require more than 8 times the mass than currently planned to reliably trigger edge-localized modes (ELMs). Unmitigated heat flux impulses from ELMs are intolerable in ITER at full power and current. Therefore, ITER operation relies on multiple approaches to control ELM heat fluxes. One method is pellet ELM pacing to instigate small rapid ELMs with low heat flux. Predicting the performance of pellet pacing is critical for ITER, which is expected to operate in a regime with a low-collisionality, peeling-limited pedestal. However, to trigger ELMs the local pressure increase in the expanding pellet cloud pushes the equilibrium over the ballooning stability limit. In this work, linear and nonlinear M3D-C1 simulations are used to predict pellet mass thresholds in DIII-D discharges and ITER scenarios with peeling-limited pedestals. It is found that the distance of the equilibrium's operational point from the ballooning branch of the pedestal stability boundary strongly changes thresholds. Linear M3D-C1 simulations find a strong dependence of the pellet mass threshold on the poloidal injection location for ITER's 15 MA, Q = 10 scenario. The required pellet mass at the planned injection locations is 8 to 17 times larger than currently considered. However, such linear simulations do not include pellet ablation physics or time evolution of density and temperature. A new scheme of 2D nonlinear simulations, coupled with linear stability analysis at various steps throughout the nonlinear time evolution, was developed to include such physics and improve on the linear results. These new nonlinear-to-linear simulations confirm previous findings. This result suggests that pellet ELM triggering in ITER could require pellets much larger than those currently planned, which makes ELM-pacing operationally challenging. On the other hand, fueling pellets injected from the high-field side will likely not unintentionally trigger ELMs in an otherwise ELM-stable plasma.

The DIII-D tokamak has elucidated crucial physics and developed projectable solutions for ITER and fusion power plants in the key areas of core performance, boundary heat and particle transport, and integrated scenario operation, with closing the core-edge integration knowledge gap being the overarching mission. New experimental validation of high-fidelity, multi-channel, non-linear gyrokinetic turbulent transport models for ITER provides strong confidence it will achieve Q 10 operation. Experiments identify options for easing H-mode access in hydrogen, and give new insight into the isotopic dependence of transport and confinement. Analysis of 2,1 islands in unoptimized low-torque IBS demonstration discharges suggests their onset time occurs randomly in the constant beta phase, most often triggered by non-linear 3-wave coupling, thus identifying an NTM seeding mechanism to avoid. Pure deuterium SPI for disruption mitigation is shown to provide favorable slow cooling, but poor core assimilation, suggesting paths for improved SPI on ITER. At the boundary, measured neutral density and ionization source fluxes are strongly poloidally asymmetric, implying a 2D treatment is needed to model pedestal fuelling. Detailed measurements of pedestal and SOL quantities and impurity charge state radiation in detached divertors has validated edge fluid modelling and new self-consistent 'pedestal-to-divertor' integrated modeling that can be used to optimize reactors. New feedback adaptive ELM control minimizes confinement reduction, and RMP ELM suppression with sustained high core performance was obtained for the first time with the outer strike point in a W-coated, compact and unpumped small-angle slot divertor. Advances have been made in integrated operational scenarios for ITER and power plants. Wide pedestal intrinsically ELM-free QH-modes are produced with more reactor-relevant conditions, Low torque IBS with W-equivalent radiators can exhibit predator-prey oscillations in T-e and radiation which need control. High-beta(P) scenarios with q(min) > 2, q(95)-7.9, beta(N) > 4, beta(T)-3.3% and H-98y2 > 1.5 are sustained with high density ((n) over bar = 7E19 m(-3), f(G)-1) for 6 tau(E), improving confidence in steady-state tokamak reactors. Diverted NT plasmas achieve high core performance with a non-ELMing edge, offering a possible highly attractive core-edge integration solution for reactors.

Tokamak plasmas with strong negative triangularity (NT) shaping typically exhibit fundamentally different edge behavior than conventional L-mode or H-mode plasmas. On DIII-D, every plasma with sufficiently negative triangularity (δ < δcrit ≃ −0.12) is found to be inherently free of edge localized modes (ELMs), even at injected powers well above the predicted L-H power threshold. It is also possible to access an ELM-free state at weaker average triangularities provided that at least one of the two x-points is still sufficiently negative. Access to the ELM-free NT scenario is found to coincide with the closure of the second stability region for infinite-n ballooning modes, suggesting that ballooning stability may play a role in limiting the accessible pressure gradient in NT plasmas. Despite this, NT plasmas are able to support small pedestals and are typically characterized by an enhancement of edge pressure gradients beyond those found in traditional L-mode plasmas. Further, the pressure gradient inside of this small pedestal is unusually steep, allowing access to high core performance that is competitive with other ELM-free regimes previously achieved on DIII-D. Since ELM-free operation in NT is linked directly to the magnetic geometry, NT fusion pilot plants are predicted to maintain advantageous edge conditions even in burning plasma regimes, potentially eliminating reactor core-integration issues caused by ELMs.

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and similar to 25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

On December 5, 2022, an indirect drive fusion implosion on the National Ignition Facility (NIF) achieved a target gain G_{target} of 1.5. This is the first laboratory demonstration of exceeding "scientific breakeven" (or G_{target}>1) where 2.05 MJ of 351 nm laser light produced 3.1 MJ of total fusion yield, a result which significantly exceeds the Lawson criterion for fusion ignition as reported in a previous NIF implosion [H. Abu-Shawareb et al. (Indirect Drive ICF Collaboration), Phys. Rev. Lett. 129, 075001 (2022)PRLTAO0031-900710.1103/PhysRevLett.129.075001]. This achievement is the culmination of more than five decades of research and gives proof that laboratory fusion, based on fundamental physics principles, is possible. This Letter reports on the target, laser, design, and experimental advancements that led to this result.

The path of tokamak fusion and International thermonuclear experimental reactor (ITER) is maintaining high-performance plasma to produce sufficient fusion power. This effort is hindered by the transient energy burst arising from the instabilities at the boundary of plasmas. Conventional 3D magnetic perturbations used to suppress these instabilities often degrade fusion performance and increase the risk of other instabilities. This study presents an innovative 3D field optimization approach that leverages machine learning and real-time adaptability to overcome these challenges. Implemented in the DIII-D and KSTAR tokamaks, this method has consistently achieved reactor-relevant core confinement and the highest fusion performance without triggering damaging bursts. This is enabled by advances in the physics understanding of self-organized transport in the plasma edge and machine learning techniques to optimize the 3D field spectrum. The success of automated, real-time adaptive control of such complex systems paves the way for maximizing fusion efficiency in ITER and beyond while minimizing damage to device components.

The filterscope diagnostic on DIII-D utilizes photomultiplier tubes to measure visible light emission from the plasma. The system has undergone a substantial upgrade since previous attempts to cross-calibrate the filterscope with other spectroscopic diagnostics were unsuccessful. The optics now utilize a dichroic mirror to initially split the light at nearly 99% transmission or reflectance for light below or above 550 nm. This allows the system to measure Dα emission without degrading visible light emission from the plasma for wavelengths below 550 nm (to measure Dβ, Dγ, W–I, C-III, etc.). Additional optimization of the optical components and calibration techniques reduce the error in the signal up to 10% in some channels compared to previous methods. Cross-calibration measurements with two other high resolution spectroscopic diagnostics now show excellent agreement for the first time. This expands the capabilities of the filterscope system allowing measurement of divertor detachment, emission profiles, edge-localized mode behavior, and plasma–wall interactions. It also enables direct comparisons against calculations from boundary plasma simulations. These were not possible before.

Edge-plasma simulations of a baffled, long-legged divertor in DIII-D, performed using the multi-fluid code UEDGE, indicate that the position of the detachment front is constrained to the location of the pump duct along the low-field side (LFS) baffle. Simulations including magnetic and E×B drifts were performed for 12.5 MW deuterium plasmas including intrinsic carbon and seeded neon to assess the optimal location of the LFS divertor pump to create a stable detachment front between the target and the X-point. The radiation front position in the simulations, taken to be indicative of the detachment front, can be controlled between the pump and X-point in the favorable magnetic field direction for H-mode access by moving the pump duct location upstream of the target along the LFS baffle. In the unfavorable magnetic field direction, the radial Eθ×B drift flows are directed towards the pumping surface, efficiently removing the injected deuterium gas and limiting the sensitivity of the radiation front location to the gas injection rate. The role of pumping rate and drift direction on the pumping efficiency are also found to affect the divertor plasma conditions and detachment front location in UEDGE simulations.

Many existing and future tokamaks with tungsten divertors operate, or will operate, with low- Z impurity seeding, but the direct effect of these seeded impurities on tungsten Scrape-off-Layer (SOL) transport has not been explored in detail. This paper reports on a DIII-D experiment designed to test how tungsten divertor leakage from the Small-Angle Slot V-Shaped, tungsten-coated divertor is impacted by neon seeding at a variety of injection rates and poloidal injection locations. Measurements from the experiment show an inverse relationship between the neon injection rate and the tungsten core penetration factor. Interpretive modeling is performed with a combination of the SOLPS-ITER and DIVIMP codes to assess the underlying tungsten behavior. The modeling results show that the reduction in tungsten divertor leakage is driven by both an increase in the divertor collisionality as well as a reduction in the ion temperature gradient near the divertor target. Collisions between low- Z impurities and tungsten impurities are found to have a significant impact on the tungsten SOL transport, such that ignoring the low- Z impurity collisional effects on the tungsten transport can result in an overestimate of the divertor leakage by an order-of-magnitude. Given the importance of these localized interactions, neon seeding from the closed, slot-like divertor has a clear advantage in being able to reduce tungsten divertor leakage without the high levels of neon core contamination that occur when seeding from other poloidal locations.

It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E x B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

UEDGE-CRUMPET simulations indicate the impact of the molecular hydrogenic isotopologue effect under high-recycling LFS divertor conditions in DIII-D to be negligible for electron density and temperature profiles at the LFS target plate. A 30% decrease in molecular content, accompanied by a 10% increase in atomic content, is predicted for deuterium compared to hydrogen. The predicted isotopologue effect on the radiative power balance, validated with calibrated spectroscopy, is found to be small despite a 20% increase in LFS divertor molecular band emission for deuterium compared to hydrogen. The predictions and measurements show a negligible contribution of molecularly-induced atomic and direct molecular emission to the total radiative power balance under high-recycling conditions, consistent with previous EDGE2D-EIRENE investigations. The UEDGE-CRUMPET simulations were performed using effective hydrogen and deuterium rates considering molecular breakup and excitation processes for H2 and D2, calculated by the CRUMPET collisional-radiative model.

An experiment was performed in the DIII-D lower divertor to validate numerical SOL tungsten (W) impurity erosion and redeposition simulations against experimental data. The net and gross erosion of W were calculated as a function of the voltage (or bias) applied to the exposed material. Five samples were inserted into the DIII-D lower divertor using the Divertor Material Evaluation System (DiMES) manipulator and exposed to constant L-mode attached plasma conditions. Each sample was partially coated with W. During plasma shots, samples were biased with respect to the machine vessel ground, ranging from −60 V to 25 V. The ERO2.0 code was used to numerically simulate the experiment aiming to compare the numerical results with experimental measures. A good agreement is found between estimated and measured tungsten erosion at least for negative biases.