Turbulent space and astrophysical plasmas exhibit a complex dynamics, which involves nonlinear coupling across different temporal and spatial scales. There is growing evidence that impulsive events, such as magnetic reconnection instabilities, lead to a spatially localized enhancement of energy dissipation, thus speeding up the energy transfer at small scales. Capturing such a diverse dynamics is challenging. Here, we employ the Multidimensional Iterative Filtering (MIF) method, a novel technique for the analysis of non-stationary multidimensional signals. Unlike other traditional methods (e.g. based on Fourier or wavelet decomposition), MIF does not require any previous assumption on the functional form of the signal to be identified. Using MIF, we carry out a multiscale analysis of Hall-magnetohydrodynamic (HMHD) and hybrid particle-in-cell (HPIC) numerical simulations of decaying plasma turbulence. The results assess the ability of MIF to spatially identify and separate the different scales (the MHD inertial range, the sub-ion kinetic and the dissipation scales) of the plasma dynamics. Furthermore, MIF decomposition allows localized current structures to be detected and their contribution to the statistical and spectral properties of turbulence to be characterized. Overall, MIF arises as a very promising technique for the study of turbulent plasma environments.

Magnetic field quantization is an important issue for degenerate environments such as neutron stars, radio pulsars and magnetars etc., due to the fact that these stars have a magnetic field higher than the quantum critical field strength of the order of $4.4\times 10^{13}~\text{G}$, accordingly, the cyclotron energy may be equal to or even much more than the Fermi energy of degenerate particles. We shall formulate here the exotic physics of strongly magnetized neutron stars, known as pulsars, specifically focusing on the outcomes of the quantized magnetic pressure. In this scenario, while following the modified quantum hydrodynamic model, we shall investigate both linear and nonlinear fast magnetosonic waves in a strongly magnetized, weakly ionized degenerate plasma consisting of neutrons and an electron–ion plasma in the atmosphere of a pulsar. Here, linear analysis depicts that sufficiently long, fast magnetosonic waves may exist in a weakly dispersive pulsar having finite phase speed at cutoff. To investigate one-dimensional nonlinear fast magnetosonic waves, a neutron density expression as a function of both the electron magnetic and neutron degenerate pressures, is derived with the aid of Riemann’s wave solution. Consequently, a modified Korteweg–de Vries equation is derived, having a rarefractive solitary wave solution. It is found that the basic properties such as amplitude, width and phase speed of the fast magnetoacoustic waves are significantly altered by the electron magnetic and the neutron degenerate pressures. The results of this theoretical investigation may be useful for understanding the formation and features of the solitary structures in astrophysical compact objects such as pulsars, magnetars and white dwarfs etc.

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in a magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry–Horton system has been proposed by St-Onge (J. Plasma Phys., vol. 83, 2017, 905830504) as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U'' \doteq {\rm d}^2 U/{\rm d} x^2 $. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $| {\rm d} U / {\rm d} x |$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry–Horton system.

A distribution of the electric potentials around a charged absorbing particle in a drifting weakly ionised collisional plasma with negative ions is calculated in the linear hydrodynamic approach. Coulomb-like asymptote of the electric potential around the absorbing particle deforms under the action of the negative ions’ flow and exhibits a valley profile along the flow behind the particle. The presence of the flowing negative ions can be conducive to string formation in the dust structures at relatively large pressures.

The isotope dependence of plasma transport will have a significant impact on the performance of future D-T experiments in JET and ITER and eventually on the fusion gain and economics of future reactors. In preparation for future D-T operation on JET, dedicated experiments and comprehensive transport analyses were performed in H, D and H-D mixed plasmas. The analysis of the data has demonstrated an unexpectedly strong and favourable dependence of the global confinement of energy, momentum and particles in ELMy H-mode plasmas on the atomic mass of the main ion species, the energy confinement time scaling as ${\tau _E}\sim {A^{0.5}}$ (Maggi et al., Plasma Phys. Control. Fusion, vol. 60, 2018, 014045; JET Team, Nucl. Fusion, vol. 39, 1999, pp. 1227–1244), i.e. opposite to the expectations based only on local gyro-Bohm (GB) scaling, ${\tau _E}\sim {A^{ - 0.5}}$, and stronger than in the commonly used H-mode scaling for the energy confinement (Saibene et al., Nucl. Fusion, vol. 39, 1999, 1133; ITER Physics Basis, Nucl. Fusion, vol. 39, 1999, 2175). The scaling of momentum transport and particle confinement with isotope mass is very similar to that of energy transport. Nonlinear local GENE gyrokinetic analysis shows that the observed anti-GB heat flux is accounted for if collisions, E × B shear and plasma dilution with low-Z impurities (9Be) are included in the analysis (E and B are, respectively the electric and magnetic fields). For L-mode plasmas a weaker positive isotope scaling ${\tau _E}\sim {A^{0.14}}$ has been found in JET (Maggi et al., Plasma Phys. Control. Fusion, vol. 60, 2018, 014045), similar to ITER97-L scaling (Kaye et al., Nucl. Fusion, vol. 37, 1997, 1303). Flux-driven quasi-linear gyrofluid calculations using JETTO-TGLF in L-mode show that local GB scaling is not followed when stiff transport (as is generally the case for ion temperature gradient modes) is combined with an imposed boundary condition taken from the experiment, in this case predicting no isotope dependence. A dimensionless identity plasma pair in hydrogen and deuterium L-mode plasmas has demonstrated scale invariance, confirming that core transport physics is governed, as expected, by the 4 dimensionless parameters ρ*, ν*, β, q (normalised ion Larmor radius, collisionality, plasma pressure and safety factor) consistently with global quasi-linear gyrokinetic TGLF calculations (Maggi et al., Nucl. Fusion, vol. 59, 2019, 076028). We compare findings in JET with those in different devices and discuss the possible reasons for the different isotope scalings reported from different devices. The diversity of observations suggests that the differences may result not only from differences affecting the core, e.g. heating schemes, but are to a large part due to differences in device-specific edge and wall conditions, pointing to the importance of better understanding and controlling pedestal and edge processes.

This paper examines the velocity distribution function and cyclotron resonance conditions for a beam of electrons moving in a magnetic field which gradually changes with time. A spatial gradient of magnetic field is known to result in an unstable horseshoe distribution of electrons. The field gradient in time adds additional effects due to an induced electric field. The resultant anisotropic velocity distribution function, which we call a Luvdisk distribution, has some distinctive properties when compared to the horseshoe. Fitting the cyclotron resonance condition circle shows that the frequency of the resultant emission is under the local cyclotron frequency. While the spatial gradient results in the emission coming almost perpendicularly to the field, the direction of the radiation under a time-changing field has more variability. The Luvdisk distribution also arises when the magnetic field has a gradient both in space and time. The beam can be unstable if those gradients are added or subtracted from each other (if the gradients are of equal or different sign), which occurs even when the total change of magnetic field is negative. While the frequency of the emission is related to the final magnetic field value, its direction is indicative of the field’s history which produced the instability.

In this paper, we propose an efficient diagnostic technique for determining spatially resolved measurements of the ion density ratio in a magnetized two-ion species plasma. Shear Alfvén waves were injected into a mixed helium–neon plasma using a magnetic loop antenna, for frequencies spanning the ion cyclotron regime. Two distinct propagation bands are observed, bounded by $\omega < \varOmega _\textrm {Ne}$ and $\omega _{ii} < \omega < \varOmega _\textrm {He}$, where $\omega _{ii}$ is the ion–ion hybrid cutoff frequency and $\varOmega _\textrm {He}$ and $\varOmega _\textrm {Ne}$ are the helium and neon cyclotron frequencies, respectively. A theoretical analysis of the cutoff frequency was performed and shows it to be largely unaffected by kinetic electron effects and collisionality, although it can deviate significantly from $\omega _{{ii}}$ in the presence of warm ions due to ion finite Larmor radius effects. A new diagnostic technique and accompanying algorithm was developed in which the measured parallel wavenumber $k_\parallel$ is numerically fit to the predicted inertial Alfvén wave dispersion in order to resolve the local ion density ratio. A major advantage of this algorithm is that it only requires a measurement of $k_\parallel$ and the background magnetic field in order to be employed. This diagnostic was tested on the Large Plasma Device at UCLA and was successful in yielding radially localized measurements of the ion density ratio.

Radiative-magnetohydrodynamic (RMHD) equations along with a full set of Maxwell's equations are followed to formulate the charged surface waves at the interface of an incompressible, radiative, magnetized dusty plasma and vacuum, while assuming that the characteristic wave frequency is much smaller than the ion gyrofrequency, having an equilibrium background state. It is found that the separation of charges on the surface is followed by thermal motion, which further leads to a negative pressure gradient normal to the surface, hence the plasma–vacuum interface is under tension due to two different types of oppositely directed pressures. The dusty plasma RMHD set of equations admits a linear dispersion relation of surface Jeans instability of an incompressible dusty plasma, which exhibits a strong coupling between the electron surface charge and dust surface mass densities and we conclude that the surface densities of both electrons and dust as well as the dust inertia play major roles in the gravitational collapse of the surface of astrophysical objects such as stars, galaxies etc. Further, the growth rate of radiative surface waves is found to be function of both the temperature inhomogeneity, appearing due to thermal radiation heat flux, as well as the thermal radiation pressure. The present findings of charged surface waves may seek application at the astroscales.

Although the interaction of a flat foil with currently available laser intensities is now considered a routine process, during the last decade, emphasis has been given to targets with complex geometries aiming at increasing the ion energy. This work presents a target geometry where two symmetric side holes and a central hole are drilled into the foil. A study of the various side-hole and central-hole length combinations is performed with two-dimensional particle-in-cell simulations for polyethylene targets and a laser intensity of $5.2\times 10^{21}~\text{W}~\text{cm}^{-2}$. The holed targets show a remarkable increase of the conversion efficiency, which corresponds to a different target configuration for electrons, protons and carbon ions. Furthermore, diffraction of the laser pulse leads to a directional high energy electron beam, with a temperature of ${\sim}40~\text{MeV}$, or seven times higher than in the case of a flat foil. The higher conversion efficiency consequently leads to a significant enhancement of the maximum proton energy from holed targets.

A new concept for the generation of thrust for space propulsion is introduced. Energetic thrust is generated in the form of plasmoids (confined plasma in closed magnetic loops) when magnetic helicity (linked magnetic field lines) is injected into an annular channel. Using a novel configuration of static electric and magnetic fields, the concept utilizes a current-sheet instability to spontaneously and continuously create plasmoids via magnetic reconnection. The generated low-temperature plasma is simulated in a global annular geometry using the extended magnetohydrodynamic model. Because the system-size plasmoid is an Alfvenic outflow from the reconnection site, its thrust is proportional to the square of the magnetic field strength and does not ideally depend on the mass of the ion species of the plasma. Exhaust velocities in the range of 20 to $500\ \textrm {km}\ \textrm {s}^{-1}$, controllable by the coil currents, are observed in the simulations.

We consider the tilting instability of a magnetically confined spheromak using three-dimensional magnetohydrodynamic and relativistic particle-in-cell calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter-alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially, the spheromak rotates – tilts – trying to lower its magnetic potential energy. As a result, a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on a few Alfvén time scales. In the case of a higher-order (second-order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfvén waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetar giant flares, fast development of tilting instabilities and no stabilization of the higher-order spheromaks make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds.

The results of flux-driven, two-fluid simulations in single-null configurations are used to investigate the processes determining the turbulent transport in the tokamak edge. Three turbulent transport regimes are identified: (i) a developed transport regime with turbulence driven by an interchange instability, which shares a number of features with the standard L-mode of tokamak operation; (ii) a suppressed transport regime, characterized by a higher value of the energy confinement time, low-amplitude relative fluctuations driven by a Kelvin–Helmholtz instability, a strong $\boldsymbol {E}\times \boldsymbol {B}$ sheared flow and the formation of a transport barrier, which recalls the H-mode; and (iii) a degraded confinement regime, characterized by a catastrophically large interchange-driven turbulent transport, which recalls the crossing of the Greenwald density limit. We derive an analytical expression of the pressure gradient length in the three regimes. The transition from the developed transport regime to the suppressed transport regime is obtained by increasing the heat source or decreasing the collisionality and vice versa for the transition from the developed transport regime to the degraded confinement regime. An analytical expression of the power threshold to access the suppressed transport regime, linked to the power threshold for H-mode access, as well as the maximum density achievable before entering the degraded confinement regime, related to the Greenwald density, are also derived. The experimental dependencies of the power threshold for H-mode access on density, tokamak major radius and isotope mass are retrieved. The analytical estimate of the density limit contains the correct dependence on the plasma current and on the tokamak minor radius.

The relationship between nonlinear large-scale dynamo action and the generation and transport of magnetic helicity is investigated at moderate values of the magnetic Reynolds number ($Rm$). The model consists of a helically forced, sheared flow in a Cartesian domain. The boundary conditions are periodic in the horizontal and impenetrable for the vertical. The magnetic field is required to be vertical at the upper and lower boundaries. There are two consequences of this choice; one is that the magnetic helicity is not gauge invariant, the second is that fluxes of magnetic helicity are allowed in and out of the domain. We select the winding gauge, define all the contributions to the evolution of the helicity in this gauge and measure these contributions for various solutions of the dynamo equations. We vary $Rm$ and the shear strength, and find a rich landscape of dynamo solutions including travelling waves, pulsating waves and non-wave-like solutions. We find that, at the $Rm$ considered, the main contribution to the growth of magnetic helicity comes from processes throughout the volume of the fluid and that boundary terms respond by limiting the growth. We find that, in this magnetic Reynolds number regime, helicity conservation is not a strong constraint on large-scale dynamo action. We speculate on what may happen at higher $Rm$.

When two collisionless plasma shells collide, they interpenetrate and the overlapping region may turn Weibel unstable for some values of the collision parameters. This instability grows magnetic filaments which, at saturation, have to block the incoming flow if a Weibel shock is to form. In a recent paper (Bret, J. Plasma Phys., vol. 82, 2016b, 905820403), it was found by implementing a toy model for the incoming particle trajectories in the filaments, that a strong enough external magnetic field $\unicode[STIX]{x1D63D}_{0}$ can prevent the filaments blocking the flow if it is aligned with them. Denoting by $B_{f}$ the peak value of the field in the magnetic filaments, all test particles stream through them if $\unicode[STIX]{x1D6FC}=B_{0}/B_{f}>1/2$. Here, this result is extended to the case of an oblique external field $B_{0}$ making an angle $\unicode[STIX]{x1D703}$ with the flow. The result, numerically found, is simply $\unicode[STIX]{x1D6FC}>\unicode[STIX]{x1D705}(\unicode[STIX]{x1D703})/\cos \unicode[STIX]{x1D703}$, where $\unicode[STIX]{x1D705}(\unicode[STIX]{x1D703})$ is of order unity. Noteworthily, test particles exhibit chaotic trajectories.

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number (${Pm}\gtrsim {1}$) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. In all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

The possibility of performing electron density and temperature measurements in a high power helicon plasma is a crucial issue in the framework of the AWAKE (Advanced WAKefield Experiment) project, which demonstrates acceleration of particles using $\text{GeV}~\text{m}^{-1}$ electric fields in plasmas. For AWAKE, a helicon is currently envisaged as a candidate plasma source due to its capability for low electron and ion temperature, high electron density and production of an elongated plasma column. A plasma diagnostic to accurately determine the electron density in AWAKE regimes would be a valuable supporting tool. A demonstration Thomson scattering (TS) diagnostic was installed and successfully tested on the resonant antenna ion device (RAID) at the Swiss Plasma Center of Ecole Polytechnique Fédérale de Lausanne. RAID produces a helicon plasma column with characteristics similar to those of the AWAKE helicon source, and is therefore an optimal testbed for application to the AWAKE device. The spectrometer employed in RAID is based on polychromators which collect the light scattered by plasma electrons in spectrally filtered wavelength regions. Results from TS on RAID demonstrate conditions of electron density and temperature respectively of $n_{e}=1.10\,(\pm 0.19)\times 10^{19}~\text{m}^{-3}$ and $T_{e}=2.3\,(\pm 0.6)~\text{eV}$ in a steady-state discharge in an Ar plasma with 5 kW of RF power. If the same polychromator system is used for AWAKE, where the electron density attained is $2\times 10^{20}~\text{m}^{-3}$, the contribution to measurement error due to coherent scattering is ${\sim}2.5\,\%$. Presented here are details of the TS diagnostic and the first tests in RAID, and the expectations for the system when employed on the AWAKE device.

Self-interaction is the process by which a microinstability eigenmode that is extended along the direction parallel to the magnetic field interacts non-linearly with itself. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics and is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near low-order mode rational surfaces (Weikl et al. Phys. Plasmas, vol. 25, 2018, 072305). We find that self-interaction, in fact, plays a very significant role in also generating fluctuating zonal flows, which is critical to regulating turbulent transport throughout the radial extent. Unlike the usual picture of zonal flow drive in which microinstability eigenmodes coherently amplify the flow via modulational instabilities, the self-interaction drive of zonal flows from these eigenmodes are uncorrelated with each other. It is shown that the associated shearing rate of the fluctuating zonal flows therefore reduces as more toroidal modes are resolved in the simulation. In simulations accounting for the full toroidal domain, such an increase in the density of toroidal modes corresponds to an increase in the toroidal system size, leading to a finite system size effect that is distinct from the well-known profile shearing effect.

It is demonstrated that finite-pressure, approximately quasi-axisymmetric stellarator equilibria can be directly constructed (without numerical optimization) via perturbations of given axisymmetric equilibria. The size of such perturbations is measured in two ways, via the fractional external rotation and, alternatively, via the relative magnetic field strength, i.e. the average size of the perturbed magnetic field, divided by the unperturbed field strength. It is found that significant fractional external rotational transform can be generated by quasi-axisymmetric perturbations, with a similar value of the relative field strength, despite the fact that the former scales more weakly with the perturbation size. High mode number perturbations are identified as a candidate for generating such transform with local current distributions. Implications for the development of a general non-perturbative solver for optimal stellarator equilibria are discussed.