Search NASA⌕ Search

SEARCH · Search NASA

Results for “collisionless”

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

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

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A magnetic analog of pressure–strain interaction

We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection.

Barbhuiya, M. Hasan (ORCID:0000000163301650)↗
DOE OSTI2025

Microtearing stability and turbulence in the pedestal: Linear gyrokinetics, reduced models, and nonlinear turbulent transport

Microtearing modes can play a crucial role in electron heat transport in tokamak plasmas, affecting both energy confinement and overall performance. This study investigates microtearing modes (MTM) stability and turbulence in a JET pedestal through gyrokinetic simulations using the Gene code, complemented by a reduced eigenvalue model. The focus is on how MTM properties depend on key plasma parameters, including collisionality and plasma beta β—the ratio of plasma pressure to magnetic pressure—the normalized toroidal wavenumber k y ρ s ⁠, where ρ s denotes the ion sound gyroradius (typically a few millimeters in edge plasmas) and isotope mass. Collisionality enhances MT growth rates, while increasing β leads to a shift from MTMs to kinetic-ballooning modes, typically for k y ρ s ⁠, where ρ s ≲ 0.2⁠. A purely collisionless branch of MTMs persists at low k y ρ s ⁠, where ρ s with distinctive properties including non-negligible particle flux and ion thermal transport. Isotope mass scans reveal modest reduction of MTM growth rates as ion mass decreases. Nonlinear simulations produce experimentally relevant transport levels. Numerical experiments turning off zonal flows and fields identify the critical role of zonal flows and zonal fields in regulating MTM turbulence. Their removal leads to a significant increase in electron heat flux. These findings provide new insight into MTM-driven transport and its impact on tokamak confinement and lay a foundation for reduced modeling and predictive capabilities.

Electrostatics↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Understanding cold electron impact on parallel-propagating whistler chorus waves via moment-based quasilinear theory

Earth's magnetosphere hosts a wide range of collisionless particle populations that interact through various wave-particle processes. Among these, cold electrons, with energies below 100 eV, often dominate the plasma density but remain poorly characterized due to measurement challenges such as spacecraft charging and photoelectron contamination. Understanding the contribution of these cold populations to wave–particle interaction is of significant interest. Recent kinetic simulations identified a secondary drift-driven instability, in which parallel-propagating whistler-mode chorus waves excite oblique electrostatic whistler waves near the resonance cone and Bernstein-mode turbulence. These secondary modes enable a new channel of energy transfer from the parallel-propagating whistler wave to the cold electrons. In this work, we develop a moment-based quasilinear theory of the secondary instabilities to quantify such energy exchange. Our results show that these secondary instabilities persist for a wide range of parameters and, in many cases, lead to nearly complete damping of the primary wave. Such secondary instability might limit the amplitude of parallel-propagating whistler waves in Earth's magnetosphere and might explain why high-amplitude oblique whistler or electron Bernstein waves are rarely observed simultaneously with high-amplitude field-aligned whistler waves in the inner magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electron heating in kinetic-Alfvén-wave turbulence

We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). In conclusion, the use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fast ion confinement in the presence of core magnetic islands in Wendelstein 7-X

The effect of magnetic islands in the core region of Wendelstein 7-X (W7-X) on fast ion confinement is explored through simulations with the BEAMS3D code. A magnetic configuration where the n/m = 5/5 island chain is shifted to r/a ~ 0.7 allows the exploration of core island physics in W7-X. The control coil system on W7-X allows the tuning of the island size either increasing the island width or decreasing it. A coupling of the BEAMS3D code to the FIELDLINES code provides a versatile mechanism for incorporating magnetic islands and stochastic regions into the BEAMS3D code. Collisionless simulations suggest that the presence of core islands degrade the confinement of passing particles in the region of the island chain. Full neutral beam simulations of W7-X show a similar behavior with confinement decreasing as the island width is increased. Comparisons between a vacuum magnetic field and low beta HINT2 simulation are made showing similar fast ion behavior. Measurements of lost fast ions in W7-X confirm this trend with the control coil suppressed island configuration showing lower losses than that with no control coils applied. Simulations of fast ion wall loads are performed suggesting no drastic change in loss pattern and a slight reduction in losses with minimized islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Sheath constraints on turbulent magnetised plasmas

Abstract A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale λ D . In magnetised plasmas where the magnetic field is inclined at an oblique angle α with the target, some of the sheath potential variation occurs also on the ion sound gyroradius length scale ρ S cos ⁡ α , caused by finite ion gyro-orbit distortion and losses. We consider a collisionless and steady-state magnetised plasma sheath whose thickness l ms ∼ max ( λ D , ρ S cos ⁡ α ) is smaller than the characteristic length scale L of spatial fluctuations in the bulk plasma, such that the limit l ms / L → 0 is appropriate. Spatial structures are assumed to be magnetic field-aligned. In the case of small magnetic field angle α ∼ δ ≡ ρ S / L ≪ 1 , electric fields tangential to the target transport ions towards the target via E × B drifts at a rate comparable to the one from parallel streaming. A generalised form of the kinetic Bohm–Chodura criterion at the sheath entrance is derived by requiring that the sheath electric field have a monotonic spatial decay far from the target. The criterion depends on tangential gradients of potential and ion distribution function, with additional nontrivial conditions.

kinetic Bohm–Chodura condition↗
DOE OSTI2021

Effects of emitting surfaces and trapped ions on the sheath physics and current flow in multidimensional plasma systems

Recent one-dimensional simulations of planar sheaths with strong electron emission have shown that trapping of charge-exchange ions causes transitions from space-charge limited (SCL) to inverse sheaths. However, multidimensional emitting sheath phenomena with collisions remained unexplored, due in part to high computational cost. As such, we developed a novel continuum kinetic code to study the sheath physics, current flow and potential distributions in two-dimensional unmagnetized configurations with emitting surfaces. For small negatively biased thermionic cathodes in a plasma, the cathode sheath can exist in an equilibrium SCL state. The SCL sheath carries an immense density of trapped ions, neutralized by thermoelectrons, within the potential well of the virtual cathode. For further increases of emitted flux, the trapped ion cloud expands in space. The trapped ion space charge causes an increase of thermionic current far beyond the saturation limit predicted by conventional collisionless SCL sheath models without ion trapping. For sufficiently strong emission, the trapped ion cloud consumes the entire 2D plasma domain, forming a mode with globally confined ions and an inverse sheath at the cathode. In situations where the emitted flux is fixed and the bias is swept (e.g. emissive probe), the trapped ions cause a large thermionic current to escape for all biases below the plasma potential. Strong suppression of the thermionic emission, required for the probe to float, only occurs when the probe is above the plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

Abstract It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage V L on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage V L at 2 MHz varies from 0 to 700 V while the HF voltage V H at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. V L = 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as V L increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher V L , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔEbecomes wider with the increase of V L . The IEDF at the radial center of the powered electrode exhibits a broader ΔEthan that at the edge. For the IADF, an increased V L results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Physics↗
DOE OSTI2023

Superfluid dark matter around black holes

Superfluid dark matter, consisting of self-interacting light particles that thermalize and condense to form a superfluid in galaxies, provides a novel theory that matches the success of the standard ΛCDM model on cosmological scales while simultaneously offering a rich phenomenology on galactic scales. Within galaxies, the dark matter density profile consists of a nearly homogeneous superfluid core surrounded by an isothermal envelope. In this work we compute the density profile of superfluid dark matter around supermassive black holes at the center of galaxies. Here, in this work, we show that, depending on the fluid equation of state, the dark matter profile presents distinct power-law behaviors, which can be used to distinguish it from the standard results for collisionless dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamical friction in dark matter superfluids: The evolution of black hole binaries

Abstract The theory of superfluid dark matter is characterized by self-interacting sub-eV particles that thermalize and condense to form a superfluid core in galaxies. Massive black holes at the center of galaxies, however, modify the dark matter distribution and result in a density enhancement in their vicinity known as dark matter spikes. The presence of these spikes affects the evolution of binary systems by modifying their gravitational wave emission and inducing dynamical friction effects on the orbiting bodies. In this work, we assess the role of dynamical friction for bodies moving through a superfluid core enhanced by a central massive black hole. As a first step, we compute the dynamical friction force experienced by bodies moving in a circular orbit. Then, we estimate the gravitational wave dephasing of the binary, showing that the effect of the superfluid drag force is beyond the reach of space-based experiments like LISA, contrarily to collisionless dark matter, therefore providing an opportunity to distinguish these dark matter models.

Astronomy & Astrophysics↗
DOE OSTI2025

SASHIMI-SIDM: semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales

We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that SASHIMI-SIDM agrees with SIDM subhalo populations from high-resolution cosmological zoom-in simulations in resolved regimes. Crucially, we predict that the fraction of core-collapsed subhalos peaks at a mass scale determined by the input SIDM cross section and decreases toward higher halo masses, consistent with the predictions of gravothermal models and cosmological simulations. For the first time, we also show that the core-collapsed fraction decreases toward lower halo masses. While the dependence of the collapse time on mass and concentration implies such behaviour, our semi-analytical approach allows us to quantify and illustrate this trend clearly across the full mass spectrum of subhalos, including for subhalo masses below the resolution limit of any current cosmological SIDM simulation. As a proof of principle, we apply SASHIMI-SIDM to predict the boost to the local dark matter density and annihilation rate from core-collapsed SIDM subhalos, which can be enhanced relative to CDM by an order of magnitude for viable SIDM models. Thus, SASHIMI-SIDM provides an efficient and reliable tool for scanning SIDM parameter space and testing it with astrophysical observations. The code is publicly available at https://github.com/shinichiroando/sashimi-si.

cosmological simulations↗
DOE OSTI2025

Abundance and properties of dark radiation from the cosmic microwave background

We study the cosmological signatures of new light relics that are collisionless like standard neutrinos or are strongly interacting. We provide a simple and succinct rephrasing of their physical effects in the cosmic microwave background, as well as the resulting parameter degeneracies with other cosmological parameters, in terms of the total radiation abundance and the fraction thereof that freely streams. In these more general terms, interacting and noninteracting light relics are differentiated by their respective decrease and increase of the free-streaming fraction, and, moreover, the scale-dependent interplay thereof with a common, correlated reduction of the fraction of matter in baryons. We then derive updated constraints on various dark-radiation scenarios with the latest cosmological observations, employing this language to identify the physical origin of the impact of each dataset. The “PR4” reanalyses of Planck CMB data prefer a larger primordial helium yield and therefore also slightly more radiation than the 2018 analysis; we investigate the differences between the two releases that drives these shifts. Smaller free-streaming fractions are disfavored by the excess lensing of the CMB measured in lensing reconstruction data from Planck and the Atacama Cosmology Telescope. On the other hand, baryon acoustic oscillation measurements from the Dark Energy Spectroscopic Instrument drive marginal detections of new, strongly interacting light relics due to that data's preference for lower matter fractions. Finally, we forecast measurements from the CMB-S4 experiment.

cosmological parameters from CMBR↗
DOE OSTI2025

The Atacama Cosmology Telescope: DR6 constraints on extended cosmological models

We use new cosmic microwave background (CMB) primary temperature and polarization anisotropy measurements from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) to test foundational assumptions of the standard cosmological model, ΛCDM, and set constraints on extensions to it. We derive constraints from the ACT DR6 power spectra alone, as well as in combination with legacy data from the Planck mission. To break geometric degeneracies, we include ACT and Planck CMB lensing data and baryon acoustic oscillation data from DESI Year-1. To test the dependence of our results on non-ACT data, we also explore combinations replacing Planck with WMAP and DESI with BOSS, and further add supernovae measurements from Pantheon+ for models that affect the late-time expansion history. We verify the near-scale-invariance (running of the spectral index dn s /d ln k = 0.0062 ± 0.0052) and adiabaticity of the primordial perturbations. Neutrino properties are consistent with Standard Model predictions: we find no evidence for new light, relativistic species that are free-streaming (N eff = 2.86 ± 0.13, which combined with astrophysical measurements of primordial helium and deuterium abundances becomes N eff = 2.89 ± 0.11), for non-zero neutrino masses (∑m ν < 0.089 eV at 95% CL), or for neutrino self-interactions. We also find no evidence for self-interacting dark radiation (N idr < 0.134), or for early-universe variation of fundamental constants, including the fine-structure constant (α EM /α EM,0 = 1.0043 ± 0.0017) and the electron mass (m e /m e,0 = 1.0063 ± 0.0056). Our data are consistent with standard big bang nucleosynthesis (we find Y p = 0.2312 ± 0.0092), the COBE/FIRAS-inferred CMB temperature (we find T CMB = 2.698 ± 0.016 K), a dark matter component that is collisionless and with only a small fraction allowed as axion-like particles, a cosmological constant (w = -0.986 ± 0.025), and the late-time growth rate predicted by general relativity (γ = 0.663 ± 0.052). We find no statistically significant preference for a departure from the baseline ΛCDM model. In fits to models invoking early dark energy, primordial magnetic fields, or an arbitrary modified recombination history, we find H 0 = 69.9 +0.8 -1.5 , 69.1 ± 0.5, or 69.6 ± 1.0 km/s/Mpc, respectively; using BOSS instead of DESI BAO data reduces the central values of these constraints by 1–1.5 km/s/Mpc while only slightly increasing the error bars. In general, models introduced to increase the Hubble constant or to decrease the amplitude of density fluctuations inferred from the primary CMB are not favored over ΛCDM by our data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

New heat flux model for non-axisymmetric divertor infrared structures

A convective heat flux model for perturbed plasmas, based on guiding center ion drift in vacuum elds (A. Wingen, et al., Phys. Plasmas 21 (2014) 012509), has been updated. The old model only considered ion heat flux, while here also electron heat flux is included. The updated model predicts divertor heat flux distributions in non-axisymmetric (3D) plasmas with applied Resonant Magnetic Perturbation (RMP) fields, and includes electric scalar potentials. It is found that a radial electric field in the near Scrape-o Layer (SOL) can considerably shift the footprints toroidally, leading to a smearing out eff ect of the incident heat flux, while a simple model for sheath potential has little impact on footprints. Various approaches to model electron heat flux are studied. A convective electron model, based on collisionless free streaming, is found to yield the best agreement with measurements, while a conductive model requires a at temperature gradient inside lobes to yield acceptable peak heat ux values. A heuristic heat flux layer approach, based on a fixed layer width also requires a limited heat flux inside the last closed flux surface (LCFS); by selecting various locations of the LCFS, the results of the conductive or convective model can be recovered respectively. The sum of ion and electron heat fluxes, both obtained by the convective model, is compared to experimental data for multiple time slices in DIII-D. Strike point splitting is observed with peak heat fluxes and layer widths that compare well to infrared camera (IR) measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗