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.

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.

The use of ultrathin solid foils offers optimal conditions for accelerating protons to high energies from laser–matter interactions. When the target is thin enough that relativistic self-induced transparency sets in, all of the target electrons get heated to high energies by the laser, which maximizes the accelerating electric field and therefore the final ion energy. In this work, we first investigate how ion acceleration by ultraintense femtosecond laser pulses in transparent CH$_2$ solid foils is modified when turning from normal to oblique ($45^\circ$) incidence. Due to stronger electron heating, we find that higher proton energies can be obtained at oblique incidence but in thinner optimum targets. We then show that proton acceleration can be further improved by splitting the laser pulse into two half-pulses focused at opposite incidence angles. An increase by ${\sim }30\,\%$ in the maximum proton energy and by a factor of ${\sim }4$ in the high-energy proton charge is reported compared to the reference case of a single normally incident pulse.

A new optimized quasi-helically symmetric configuration is described that has the desirable properties of improved energetic particle confinement, reduced turbulent transport by three-dimensional shaping and non-resonant divertor capabilities. The configuration presented in this paper is explicitly optimized for quasi-helical symmetry, energetic particle confinement, neoclassical confinement and stability near the axis. Post optimization, the configuration was evaluated for its performance with regard to energetic particle transport, ideal magnetohydrodynamic stability at various values of plasma pressure and ion temperature gradient instability induced turbulent transport. The effects of discrete coils on various confinement figures of merit, including energetic particle confinement, are determined by generating single-filament coils for the configuration. Preliminary divertor analysis shows that coils can be created that do not interfere with expansion of the vessel volume near the regions of outgoing heat flux, thus demonstrating the possibility of operating a non-resonant divertor.

Absorbing boundary conditions or layers are used in simulations to reduce or eliminate wave reflections from the boundary; one of the most widely used absorbing layers is Berenger's perfectly matched layer (PML). In this paper, PML is extended to a compound absorbing layer which has multiple effects of damping and deceleration, and is applied to linear and nonlinear Schrödinger equations. The deceleration extends the time to damp out the modes with higher phase velocities, leading to remarkably reduced total reflection for dispersive waves. By invoking the two effects independently, the flexibility and performance are enhanced. Since this method is based on the WKB formalism, it requires an absorbing layer of a moderate size.

We report measurements of the time-averaged surface floating potential of magnetic and non-magnetic spherical probes (or large dust particles) immersed in a magnetized capacitively coupled discharge. In this study, the size of the spherical probes is taken greater than the Debye length. The surface potential of a spherical probe first increases, i.e. becomes more negative at low magnetic field ($B < 0.05\ \textrm {T}$), attains a maximum value and decreases with further increase of the magnetic field strength ($B > 0.05\ \textrm {T}$). The rate of change of the surface potential in the presence of a $B$-field mainly depends on the background plasma and types of material of the objects. The results show that the surface potential of the magnetic sphere is higher (more negative) compared with the non-magnetic spherical probe. Hence, the smaller magnetic sphere collects more negative charges on its surface than a bigger non-magnetic sphere in a magnetized plasma. The different sized spherical probes have nearly the same surface potential above a threshold magnetic field ($B > 0.03\ \textrm {T}$), implying a smaller role of size dependence on the surface potential of spherical objects. The variation of the surface potential of the spherical probes is understood on the basis of a modification of the collection currents to their surface due to charge confinement and cross-field diffusion in the presence of an external magnetic field.

A disruption database characterizing the current quench of disruptions with ITER-like tungsten divertor has been developed on EAST. It provides a large number of plasma parameters describing the predisruptive plasma, current quench time, eddy current, and mitigation by massive impurity injection, which shows that the current quench time strongly depends on magnetic energy and post-disruption electron temperature. Further, the energy balance and magnetic energy dissipation during the current quench phase has been well analysed. Magnetic energy is also demonstrated to be dissipated mainly by ohmic reheating and inductive coupling, and both of the two channels have great effects on current quench time. Also, massive gas injection is an efficient method to speed up the current quench and increase the fraction of impurity radiation.

We have recently demonstrated that by expanding in small distance from the magnetic axis compared with the major radius, stellarator shapes with low neoclassical transport can be generated efficiently. To extend the utility of this new design approach, here we evaluate measures of magnetohydrodynamic interchange stability within the same expansion. In particular, we evaluate the magnetic well, Mercier's criterion, and resistive interchange stability near a magnetic axis of arbitrary shape. In contrast to previous work on interchange stability near the magnetic axis, which used an expansion of the flux coordinates, here we use the ‘inverse expansion’ in which the flux coordinates are the independent variables. Reduced expressions are presented for the magnetic well and stability criterion in the case of quasisymmetry. The analytic results are shown to agree with calculations from the VMEC equilibrium code. Finally, we show that near the axis, Glasser, Greene and Johnson's stability criterion for resistive modes approximately coincides with Mercier's ideal condition.

The fundamental difference between incompressible ideal magnetohydrodynamics and the dynamics of a non-conducting fluid is that magnetic fields exert a tension force that opposes their bending; magnetic fields behave like elastic strings threading the fluid. It is natural, therefore, to expect that a magnetic field tangled at small length scales should resist a large-scale shear in an elastic way, much as a ball of tangled elastic strings responds elastically to an impulse. Furthermore, a tangled field should support the propagation of ‘magnetoelastic waves’, the isotropic analogue of Alfvén waves on a straight magnetic field. Here, we study magnetoelasticity in the idealised context of an equilibrium tangled field configuration. In contrast to previous treatments, we explicitly account for intermittency of the Maxwell stress, and show that this intermittency necessarily decreases the frequency of magnetoelastic waves in a stable field configuration. We develop a mean-field formalism to describe magnetoelastic behaviour, retaining leading-order corrections due to the coupling of large- and small-scale motions, and solve the initial-value problem for viscous fluids subjected to a large-scale shear, showing that the development of small-scale motions results in anomalous viscous damping of large-scale waves. Finally, we test these analytic predictions using numerical simulations of standing waves on tangled, linear force-free magnetic-field equilibria.

Two-dimensional spatial measurements of magnetic field and plasma parameters have been performed in conventional magnetron DC discharges during the formation of metallic nanoparticles. Correlations between the electron density and temperature distributions, and the magnetic field geometry and strength have been established. A sharp increase of the plasma potential is found on the edge of the last magnetic arch followed by a decrease towards the anode plate and edges. It is shown that the spatial variation of the plasma potential is at the origin of a potential well that can trap negatively charged nanoparticles.

In the present analysis, we study effects of the radiation reaction (RR) on the dynamics of charged particles submitted to the action of localized longitudinal high-frequency carriers travelling at the speed of light. As the wave's crests and troughs keep overtaking particles, dissipative RR forces tend to drag particles alongside the wave in an effort to reduce the relative wave–particle speed. Particles of course never reach the phase velocity of the wave, but are instead driven to an ever-growing velocity, towards the speed of light, while in the wave localization region. We developed a modified average Hamiltonian formalism capable of describing the intricacies of the corresponding dynamics. The modified formalism agrees with simulations and is of particular usefulness in the study of optimum values for the localization length and maximum wave amplitude.

The present work considers the stability of a high-$\beta$, large aspect ratio, circular plasma with diffuse profiles for the safety factor and the angular toroidal frequency (López & Guazzotto, Phys. Plasmas, vol. 24, 032501). An application of the Frieman–Rotenberg formalism results in a system of scalar eigenmode equations whose coupling is retained at the plasma–vacuum transition but is disregarded across the plasma column, which is a standard practice. The solution technique consists of a multidimensional shooting method for the poloidal harmonics; robust initial guesses are constructed by solving the dispersion relation in the static scenario with vanishing magnetic shear. Flow shear appears as a high-$\beta$ toroidal contribution, and we illustrate its destabilizing influence on $n=1$ external kink modes in the presence of ideal and resistive walls. Internal resonances are avoided by means of the selection of appropriate equilibrium parameters. The stabilizing influence of a finite positive average magnetic shear is also exemplified.

The paper presents experimental results from the SMOLA device that is the first facility with a helical mirror section of the magnetic field. This device is built in the Budker Institute of Nuclear Physics for the verification of the helical mirror confinement idea that is the recently introduced technique of the active control of axial losses from a confinement zone. Theory predicts that with rotating plasma, a helical mirror will provide suppression of the axial plasma flow and, simultaneously, density pinching to the axis. Experiments demonstrated that plasma density at the exit from the transport section is suppressed with activation of the helical field, the effect is significant and highly reproducible. The most pronounced effect is observed on the plasma periphery, where the mirror ratio is the highest. The integral suppression ratio reaches 2–2.5 in the discussed experiments. Experimental results are compared with simplified theoretical estimates. The integral suppression ratio matches the simple theoretical estimates even if the transversal diffusion is neglected.

We present a Monte Carlo collisional scheme that models single Compton scattering between leptons and photons in particle-in-cell codes. The numerical implementation of Compton scattering can deal with macro-particles of different weights and conserves momentum and energy in each collision. Our scheme is validated through two benchmarks for which exact analytical solutions exist: the inverse Compton spectra produced by an electron scattering with an isotropic photon gas and the photon–electron gas equilibrium described by the Kompaneets equation. It provides new opportunities for numerical investigation of plasma phenomena where a significant population of high-energy photons is present in the system.

In support of the ICRF experiments planned on the Wendelstein 7-X (W7-X) stellarator, i.e. fast ion generation, wall conditioning, target plasma production and heating, a first experimental study on plasma production has been made in the Uragan-2M (U-2M) stellarator using W7-X-like two-strap antenna. In all the experiments, antenna monopole phasing was used. The W7-X-like antenna operation with launched radiofrequency power of ~100 kW have been performed in helium (p = (4–14) × 10−2 Pa) with the vacuum vessel walls pre-loaded with hydrogen. Production of plasma with a density higher than 1012 cm−3 was observed near the first harmonic of the hydrogen cyclotron frequency. Operation at first hydrogen harmonic is feasible in W7-X future ICRF experiments.

Utilizing the similarity between the spinor representation of the Dirac and the Maxwell equations that has been recognized since the early days of relativistic quantum mechanics, a quantum lattice algorithm (QLA) representation of unitary collision-stream operators of Maxwell's equations is derived for both homogeneous and inhomogeneous media. A second-order accurate 4-spinor scheme is developed and tested successfully for two-dimensional (2-D) propagation of a Gaussian pulse in a uniform medium whereas for normal (1-D) incidence of an electromagnetic Gaussian wave packet onto a dielectric interface requires 8-component spinors because of the coupling between the two electromagnetic polarizations. In particular, the well-known phase change, field amplitudes and profile widths are recovered by the QLA asymptotic profiles without the imposition of electromagnetic boundary conditions at the interface. The QLA simulations yield the time-dependent electromagnetic fields as the wave packet enters and straddles the dielectric boundary. QLA involves unitary interleaved non-commuting collision and streaming operators that can be coded onto a quantum computer: the non-commutation being the very reason why one perturbatively recovers the Maxwell equations.

Statistical properties of the scrape-off layer plasma fluctuations are studied in ohmically heated plasmas in the Alcator C-Mod tokamak. For the first time, plasma fluctuations as well as parameters that describe the fluctuations are compared across measurements from a mirror Langmuir probe (MLP) and from gas-puff imaging (GPI) that sample the same plasma discharge. This comparison is complemented by an analysis of line emission time-series data, synthesized from the MLP electron density and temperature measurements. The fluctuations observed by the MLP and GPI typically display relative fluctuation amplitudes of order unity together with positively skewed and flattened probability density functions. Such data time series are well described by an established stochastic framework that models the data as a superposition of uncorrelated, two-sided exponential pulses. The most important parameter of the process is the intermittency parameter, $\gamma = {\tau _{d}} / {\tau _{w}}$, where ${\tau _{d}}$ denotes the duration time of a single pulse and ${\tau _{w}}$ gives the average waiting time between consecutive pulses. Here we show, using a new deconvolution method, that these parameters can be consistently estimated from different statistics of the data. We also show that the statistical properties of the data sampled by the MLP and GPI diagnostic are very similar. Finally, a synthetic GPI signal using only plasma parameters sampled by the MLP shows qualitatively different fluctuation statistics from the measured GPI signal.

The multiphase simulations are conducted with the kinetic-magnetohydrodynamics hybrid code MEGA to investigate the spatial and the velocity distributions of lost fast ions due to the Alfvén eigenmode (AE) bursts in the Large Helical Device plasmas. It is found that fast ions are lost along the divertor region with helical symmetry both before and during the AE burst except for the promptly lost particles. On the other hand, several peaks are present in the spatial distribution of lost fast ions along the divertor region. These peaks along the divertor region can be attributed to the deviation of the fast-ion orbits from the magnetic surfaces due to the grad-B and the curvature drifts. For comparison with the velocity distribution of lost fast ions measured by the fast-ion loss detector (FILD), the ‘numerical FILD’ which solves the Newton–Lorentz equation is constructed in the MEGA code. The velocity distribution of lost fast ions detected by the numerical FILD during AE burst is in good qualitative agreement with the experimental FILD measurements. During the AE burst, fast ions with high energy (100–180 keV) are detected by the numerical FILD, while co-going fast ions lost to the divertor region are the particles with energy lower than 50 keV.

Here, we propose a thermodynamic model for dusty plasma, where the dust is confined in a small volume within a large plasma background by external fields. In this model, the parameters of dust, e.g. Helmholtz energy, pressure, entropy and enthalpy, etc. can be calculated for given dust density and temperature. The model is solved analytically in the mean field (gaseous) limit and various processes associated with the gaseous phase of dust, e.g. adiabatic/isothermal/constant internal energy expansion/compression, specific heat, free expansion within the plasma background, and the dispersion of novel acoustic waves are studied. Some predictions of the model, e.g. electrostatic pressure of the dust and the isothermal equation of state, have been earlier verified in experiments and numerical simulations. The model is compared with an earlier thermodynamic model of dusty plasma proposed by Hamaguchi and Farouki.