The high-power helicon wave system in the DIII-D tokamak introduces new plasma–material interaction (PMI) challenges due to rectified RF sheath potentials forming near antenna structures and surrounding tiles. Using the STRIPE modeling framework-which integrates SOLPS-ITER, COMSOL, RustBCA, and GITR/GITRm-we simulate carbon erosion, re-deposition, and global impurity transport in two H-mode discharges with varying antenna–plasma gaps and RF powers. COMSOL predicts rectified sheath potentials of 1-5 kV, localized near the bottom of the antenna where magnetic field lines intersect at grazing angles. Erosion is dominated by carbon self-sputtering, with RF-accelerated D+ ions contributing up to 1

Experiments on DIII-D have demonstrated a density-confinement synergy that enables sustainment of high performance in a previously unattained parameter regime of simultaneous very high energy confinement quality (H _98y2 ⩾ 1.5), very high line-average density Greenwald fraction ( ƒ _Gr = πa ^2 < n >/ I _P ⩾ 1.4), and high toroidal beta ( ${\beta _{\text{T}}}$ ⩾ 3%). Tokamak operation in this regime is essential for a compact steady-state FPP, as well as for Q = 10 with 500 MW of fusion power in ITER at I _P ≪ 15 MA. These experiments leveraged the knowledge that, in the high-poloidal-beta ( ${\beta _{\text{P}}}$ ) regime, density gradients can enhance turbulence stabilization caused by high ${\alpha _{{\text{MHD}}}}$ ( ${\alpha _{{\text{MHD}}}}\sim {\text{d}}{\beta _{\text{P}}}/{\text{d}}r$ ). This was described by theoretical predictions and gyrokinetic transport simulations (Kotschenreuther et al 2024 Nucl. Fusion 64 076033), and later confirmed by experiments on DIII-D (Ding et al 2024 Nature 629 555). To increase both ${\beta _{\text{P}}}$ and ${\beta _{\text{T}}}$ , the new experiments increased the ideal-wall stability ${\beta _{\text{N}}}$ -limit by using a smaller plasma-outer wall distance and higher triangularity in the plasma cross section (top/bottom average δ ∼ 0.9), enabled by the recent ‘shape & volume rise’ (SVR) modification to the DIII-D divertor. The higher triangularity also contributed to achieving higher ƒ _Gr by enabling higher pedestal density. At high density, the pedestal is ballooning limited and exhibits small and frequent ELMs, while the divertor is near detachment even without any impurity seeding. Highplasma performance was attained and sustained reproducibly, with the eventual terminations brought about by an MHD mode destabilized as the current profile slowly continued to evolve. A path to stationary fully noninductive operation might include ECH injection to reduce both core impurity accumulation and the electron collisionality, thus increasing the bootstrap current. These experiments provide the first experimental demonstration of the ƒ _Gr , H _98 _y _2 , and ${\beta _{\text{T}}}$ values required simultaneously for ITER Q = 10 at I _P < 10 MA, pointing to practical ways to improve the energy confinement in a fusion reactor.

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high- β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations ( $\delta B_{\parallel}$ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

Helicon waves (also known as whistler waves) satisfying the normal wave-particle cyclotron resonance are observed to limit the growth and maximum energy of relativistic electrons (REs) in low-density Ohmic DIII-D tokamak plasmas. Following the application of helicon waves, pitch-angle scattering of high-energy REs causes an increase in both synchrotron and electron-cyclotron emissions. The hard x-ray emission, a proxy for the RE population, ceases to grow. Energy-resolved hard x-ray measurements also show a striking decrease in the number of high-energy REs (above the resonance at approximately 8 MeV) to below the noise floor and an increase in low-energy (∼4  MeV) REs. This occurs despite the toroidal electric field remaining high enough to drive exponential RE growth in the absence of helicon waves. These results open new directions for limiting the maximum energy of RE populations in laboratory and fusion plasmas.

A new frequency-stepped Doppler backscattering (DBS) system has been integrated into a real-time steerable electron cyclotron heating launcher system to simultaneously probe local background turbulence (f < 10 MHz) and high-frequency (20-550 MHz) density fluctuations in the DIII-D tokamak. The launcher allows for 2D steering (horizontally and vertically) over wide angular ranges to optimize probe location and wavenumber response. The vertical steering can be optimized during a discharge in real time. The new DBS system employs a programmable frequency synthesizer with adjustable dwell time as a source to launch either O or X-mode polarized millimeter waves. This system can step in real-time over the entire E-band frequency range (60-90 GHz). This combination of capabilities allows for the diagnosis of the complex internal spatial structure of high power (>200 kW) helicon waves (476 MHz) injected from an external antenna during helicon current drive experiments in DIII-D. Broadband density fluctuations around the helicon frequency are observed during real-time scans of measurement location and wavenumber during these experiments. Analysis indicates that these broadband high-frequency fluctuations are a result of backscattering of the DBS millimeter-wave probe beam from plasma turbulence modulated by the helicon wave. It is observed that background turbulence is effectively locally "tagged" with the helicon wave electric field, forming images of the turbulent spectrum in the overall density fluctuation spectrum that appear as high-frequency sidebands of the turbulence. These observations of background turbulence and high-frequency fluctuations open up the possibility of monitoring local helicon wave amplitude by comparing the high-frequency signal amplitude to the simultaneously measured background turbulence. In combination with the real-time measurement location and wavenumber scanning capabilities (offered by real-time frequency-stepping and steering), this allows rapid determination of the spatial distribution of the helicon wave power during steady-state plasma operation. In the long term, such measurements may be used to validate predictive modeling (GENRAY [Smirnov and Harvey, Bull. Am. Phys. Soc. 40, 1837 (1995)] or AORSA [Lau et al., Nucl. Fusion 58 066004 (2018)]) of helicon current drive in DIII-D plasmas.

A set of experiments were conducted on the LArge Plasma Device (LAPD) at UCLA to test the operational principles of a traveling wave antenna of the comb-line type. This antenna was designed to launch helicon waves (fast waves in the lower hybrid range of frequencies) on DIII-D. With the order-of-magnitude lower static magnetic field on LAPD, the antenna excites waves in a different regime. Whenever fast waves can propagate in LAPD, slow waves are also supported by the plasma so it is necessary to distinguish between the two cold-plasma branches in evaluating the effectiveness of the launcher. The results show that the launcher couples well to fast waves when the plasma supports fast-wave propagation; control of the principal imposed parallel wavenumber can be achieved through varying the launch frequency on the antenna within its bandwidth of operation; and that the launched waves exhibit strong directionality. We also investigate the role of the plasma profile and wave mode on the loading characteristics. Additionally, a comparison with full-wave modeling of the propagating waves is shown using both a cold-plasma model in COMSOL and a hot-plasma model in RFPisa, which obtain similar results in the present regime.

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker-Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k phi phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

Ion cyclotron emission (ICE) is widely observed from spatially localised minority energetic ion populations in toroidal magnetically confined fusion (MCF) plasmas, both tokamaks and stellarators. Its spectral structure is typically regular with narrow suprathermal peaks, whose frequency separation matches a local energetic ion cyclotron frequency. Here we report the first computational (fully nonlinear self-consistent kinetic particle-in-cell code) and analytical (linear magnetoacoustic cyclotron instability (MCI)) studies of ICE observations from cylindrical plasmas contained in the Large Plasma Device (LAPD). Because LAPD is cylindrical, the plasma physics giving rise to the observed ICE spectrum necessarily excludes toroidal effects. Our approach, previously successful for toroidal plasma ICE, assumes slab geometry and hence is well adapted to LAPD. ICE from LAPD is strongly electrostatic, as distinct from electromagnetic, and is driven by 15 keV beam ions for which the ratio of perpendicular speed to the local Alfven speed, v perpendicular to/vA, is 0.15, lower than in MCF plasmas from which beam-driven ICE has previously been observed. Our results are in good agreement with these observations. There is congruence between simulated ICE spectra, obtained in the saturated nonlinear regime of our computations, and observed ICE spectra; the underlying physics is essentially the same as in toroidal plasmas; and there is alignment with linear analytical theory where appropriate. The present work establishes a mapping from the cylindrical LAPD ICE observations to toroidal MCF ICE observations. The LAPD spectra are instances of beam-driven MCI-type ICE in its sub-Alfvenic, predominantly electrostatic manifestation, which has precedents in MCF stretching back to the 1990s. An interesting corollary is that, for many purposes, ICE in toroidal MCF plasmas 'might as well' be occurring in a cylinder.

A set of experiments using controlled, skin depth-sized plasma pressure filaments in close proximity have been carried out in a large linear magnetized plasma device. Two- and three-filament configurations have been used to determine the scale of cross field nonlinear interaction. When the filaments are separated by a distance of approximately five times the size of a single filament or less, a significant transfer of charge and energy occurs, leading to the generation of inter-filament electric fields. This has the effect of rotating the filaments and influencing the merging dynamics. Nonlinear gyrokinetic simulations using seeded filaments confirm the presence of unstable drift-Alfvén modes driven by the steep electron temperature gradient. When the filaments are within a few collisionless electron skin depths (separations twice the size of a single filament), the unstable perturbations drive the convective mixing of the density and temperature and rearrange the gradients such that they maximize in the region surrounding the filament bundle.

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.

The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

Helicon waves are thought to be promising in various tokamaks, such as DIII-D, because they can penetrate reactor-grade high-density cores and drive the off-axis current with higher efficiency. In the frequency regime similar to 476 MHz, both slow electrostatic and fast electromagnetic helicon waves can coexist in DIII-D. If the antenna parasitically excites the slow mode, these waves can propagate along the magnetic field line into the scrape-off layer (SOL). Although the importance of the misalignment of the Faraday screen and the electron density in the SOL on the excitation and propagation of slow modes is well known, the conditions for minimizing slow mode excitation have yet to be optimized. Using the Petra-M simulation code in the 2D domain, we analyze the effects of the misalignment of the antenna in the poloidal direction, the misalignment of the Faraday screen in the toroidal direction, and the density in front of the antenna on slow mode generation. Our results suggest that the misalignment of the Faraday screen is a critical factor in reducing the slow mode and that the misalignment angle should be below similar to 5 degrees to minimize the slow wave excitation. When the electron density is higher than 3.5x10(18 )m(-3) in the SOL, the generation of the slow mode from the antenna is minimized and unaffected by the misalignment of the Faraday screen.

More than 0.6 MW of rf power at 476 MHz has been coupled to DIII-D plasmas by launching helicon (whistler) waves with a traveling-wave antenna (comb-line) in the fast-wave polarization (Van Compernolle et al 2021 Nucl. Fusion 61 116034) which resulted in the observation of electron heating of the core plasma with single-pass absorption based on ray-tracing in L-mode discharges. The coupling performance of the 1.5 m wide 30-element comb-line traveling-wave antenna has been consistent with expectations based on the 2015-2016 experiments on DIII-D with a low-power 12-element prototype (Pinsker et al 2018 Nucl. Fusion 58 106007). The conditioning process that was necessary to carry out high-power experiments is discussed; rf-specific impurities have not been observed. Parametric decay instabilities have been observed and are being investigated as a potential edge absorption mechanism (Porkolab et al 2023 AIP Conf. Proc. 2984 070004).

A high-level design study for a new experimental tokamak shows that advances in fusion science and engineering can be leveraged to narrow the gaps in energy confinement and exhaust power handling that remain between present devices and a future fusion pilot plant (FPP). This potential new U.S. facility, an Exhaust and Confinement Integration Tokamak Experiment (EXCITE), will access an operational space close to the projected FPP performance regime via a compact, high-field, high-power-density approach that utilizes advanced tokamak scenarios and high-temperature superconductor magnets. Full-device optimization via system code calculations, physics-based core-edge modeling, plasma control simulations, and finite element structural and thermal analysis has converged on a BT=6 T, IP=5 MA, R0=1.5 m, A=3, D-D tokamak with strong plasma shaping, long-legged divertors, and 50 MW of auxiliary power. Such a device will match several absolute FPP parameters: plasma pressure, exhaust heat flux, and toroidal magnetic field. It will also narrow or close the gap in key dimensionless parameters: toroidal beta, bootstrap fraction, collisionality, and edge neutral opacity. Integrated neutron shielding preserves personnel access by limiting nuclear activation and maximizes experimental run time by reducing site radiation. In addition to design study results and optimization details, parameter sensitivities and uncertainties are also discussed.

A novel quadrature Doppler Backscattering (DBS) system has been developed and optimized for the E-band (60-90 GHz) frequency range using either O-mode or X-mode polarization in DIII-D plasmas. In general, DBS measures the amplitude of density fluctuations and their velocity in the lab frame. The system can simultaneously monitor both low-frequency turbulence (f < 10 MHz) and radiofrequency plasma density fluctuations over a selectable frequency range (20-500 MHz). Detection of high-frequency fluctuations has been demonstrated for low harmonics of the ion cyclotron frequency (e.g., 2fci ∼ 23 MHz) and externally driven high-frequency helicon waves (f = 476 MHz) using an adjustable frequency down conversion system. Importantly, this extends the application of DBS to a high-frequency spectral domain while maintaining important turbulence and flow measurement capabilities. This unique system has low phase noise, good temporal resolution (sub-millisecond), and excellent wavenumber coverage (kθ ∼ 1-20 cm-1 and kr ≲ 30 cm-1). As a demonstration, localized internal DIII-D plasma measurements are presented from turbulence (f ≤ 5 MHz), Alfvenic waves (f ∼ 6.5 MHz), ion cyclotron waves (f ≥ 20 MHz), as well as fluctuations around 476 MHz driven by an external high-power 476 MHz helicon wave antenna. In the future, helicon measurements will be used to validate GENRAY and AORSA modeling tools for prediction of helicon wave propagation, absorption, and current drive location for the newly installed helicon current drive system on DIII-D.

In this paper, we present an overview of the development of two high power traveling wave antennas (TWAs) termed the “comb-line” and the “finline” for launching fast and slow waves in plasmas. The comb-line antenna which has been deployed at the DIII-D National Fusion Facility tokamak can couple megawatt (MW) level RF power into the tokamak plasma with $\sim$ 2% reflected power and $\sim$ 2% dissipated power per array element at its design frequency, 476 MHz. A prototype of the finline antenna which has been designed to launch lower hybrid slow waves at 4.6 GHz, has reflection coefficient below 10% with 1.8 dB insertion loss.

The ponderomotive force has previously been identified as a possible driver of observed density modifications close to radio frequency actuators during operation. This nonlinear force redistributes density in regions of gradients in the magnitude of an oscillating electric field and describes the influence of the fast time scale RF wave dynamics on slow time scale plasma transport. Depletion of the saturation current (a proxy for the density) measured at the Large Plasma Device (LAPD) was 30–35% during ion cyclotron range of frequencies operation. A coupled 1D plasma transport and cold plasma frequency domain wave solver was developed to self-consistently describe ponderomotive effects and was used to compare with results obtained from the LAPD experiment. The scaled current density driver for the wave model yielded an RF B field in close agreement with two components of the experimental data. However, the 1D parallel model did not accurately reproduce the amplitude or spatial distribution observed in experimental measurements of By. Within the limitations of the 1D model, initial simulation results showed that the ponderomotive force depleted up to 8% for high power (1 MW) and around 1% for the experimental power of 120 kW. This could suggest that the ponderomotive force is not the main driver of density modification for the LAPD experiments presented in this paper. Higher fidelity tools of at least 2D will be required to give a more realistic description of the RF E fields and the effect of the ponderomotive force on the LAPD.