Where and under what conditions the transfer of energy between electromagnetic fields and particles takes place in the solar wind remains an open question. We investigate the conditions that promote the growth of kinetic instabilities predicted by linear theory to infer how turbulence and temperature-anisotropy-driven instabilities are interrelated. Using a large dataset from Solar Orbiter, we introduce the radial rate of strain, a novel measure computed from single-spacecraft data, which we interpret as a proxy for the double-adiabatic strain rate. The solar wind exhibits high absolute values of the radial rate of strain at locations with large temperature anisotropy. We measure the kurtosis and skewness of the radial rate of strain from the statistical moments to show that it is non-Gaussian for unstable intervals and increasingly intermittent at smaller scales with a power-law scaling. We conclude that the velocity field fluctuations in the solar wind contribute to the presence of temperature anisotropy sufficient to create potentially unstable conditions.

The propagation of electromagnetic waves in a linearly varying index of refraction is a fundamental problem in wave physics, being relevant in fusion science for describing certain wave-based heating and diagnostic schemes. Here, an exact solution is obtained for a given incoming wavefield specified on the boundary transverse to the direction of inhomogeneity by performing a spectral, rather than asymptotic, matching. Two case studies are then presented: a Gaussian beam at oblique incidence and a speckled wavefield at normal incidence. For the Gaussian beam, it is shown that when the waist $W$ is sufficiently large, oblique incidence manifests simply as rigid translation and focal shift of the corresponding diffraction pattern at normal incidence. The destruction of the hyperbolic umbilic caustic (corresponding to a critically focused beam) as $W$ is reduced is then demonstrated. The caustic disappears once $W \lesssim \delta _a \sqrt {L}$ ($L$ being the medium length scale normalized by the Airy skin depth $\delta _a$), at which point the wave behaviour is increasingly described by Airy functions, but experiences less focusing as a result. To maximize the intensity of a launched Gaussian beam at a turning point, one should therefore minimize the imaginary part of the launched complex beam parameter while having the real part satisfy critical focusing. For the speckled wavefield, it is shown that the transverse speckle pattern only couples to the Airy longitudinal pattern when the coupling parameter $\eta = \sqrt {L}/f_{\#}$ is large, with $f_{\#}$ being the f-number of the launching aperture. When $\eta \ll 1$, a reduced description of the total wavefield can be obtained by simply multiplying the incoming speckle pattern with the Airy swelling.

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

The response tensor is derived for a relativistically streaming, strongly magnetized, one-dimensional Jüttner distribution of electrons and positrons, referred to as a pulsar plasma. This is used to produce a general treatment of wave dispersion in a pulsar plasma. Specifically, relativistic streaming, the spread in Lorentz factors in a pulsar rest frame and cyclotron resonances are taken into account. Approximations to the response tensor are derived by making approximations to relativistic plasma dispersion functions appearing in the general form of the response tensor. The cold-plasma limit, the highly relativistic limit and limits related to cyclotron resonances are considered. The theory developed in this paper has applications to generalized Faraday rotation in pulsars and magnetars.

We present a framework for analysing plasma flow in a rotating mirror. By making a series of physical assumptions, we reduce the magnetohydrodynamic (MHD) equations in a three-dimensional cylindrical system to a one-dimensional system in a shallow, cuboidal channel within a transverse magnetic field, similar to the Hartmann flow in ducts. We then solve the system both numerically and analytically for a range of values of the Hartmann number and calculate the dependence of the plasma flow speed on the thickness of the insulating end cap. We observe that the mean flow overshoots and decelerates before achieving a steady-state value, a phenomenon that the analytical model cannot capture. This overshoot is directly proportional to the thickness of the insulating end cap and the external electric field, with a weak dependence on the external magnetic field. Our simplified model can act as a benchmark for future simulations of the supersonic mirror device CMFX (centrifugal magnetic fusion experiment), which will employ more sophisticated physics and realistic magnetic field geometries.

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, $\beta \gg 1$, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-$\beta$, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (${\sim }5\,\%$) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

We draw attention to an interesting possibility in the design and operation of stellarator fusion reactors, which has hitherto been considered unrealistic under burning-plasma conditions. Thanks to recent advances in stellarator optimisation theory, it appears possible to create a positive (outward-pointing) radial electric field in the plasma core by carefully tailoring the geometry of the magnetic field. This electric field is likely to expel highly charged impurities from the centre of the plasma through neoclassical transport and thus eliminate, or at least mitigate, a long-standing problem in stellarator physics. Further out, the electric field is expected to suddenly change sign from positive to negative, thus creating a region of strongly sheared flow, which could locally suppress turbulent transport and enhance overall energy confinement.

In magnetic confinement fusion devices close to axisymmetry, such as tokamaks, a key element is the winding profile of the magnetic field lines, or its inverse, the safety profile $q=q_{\boldsymbol {B}}$. A corresponding profile, $q_{\boldsymbol {J}}$, can be defined for the current density field lines. Ampère's law relates any mode of current perturbation $\delta \boldsymbol {J}_{m,n}$ with a mode of magnetic perturbation $\delta \boldsymbol {B}_{m,n}$. It is shown that the knowledge of the pair $(q_{\boldsymbol {B}},q_{\boldsymbol {J}})$ allows us then to characterize the resonant, or non-resonant, nature of the modes for both the magnetic and current density field lines. The expression of $q_{\boldsymbol {J}}$ in the flux coordinate is derived. Including this calculation in real-time Grad–Shafranov equilibrium reconstruction codes would yield a comprehensive view of the magnetics. The monitoring of the pair $(q_{\boldsymbol {B}},q_{\boldsymbol {J}})$ would then allow us to investigate the role played by the resonant modes for the current density, that are current filamentary modes, in the plasma small-scale turbulence. By driving the magnetic and current density profiles apart so that the images of $q_{\boldsymbol {B}}$ and $q_{\boldsymbol {J}}$ are disjoint, these filamentary modes would not impact the magnetic field topology, being not associated with magnetic islands but with non-resonant magnetic modes. It remains to be explored to what extent such a configuration, where the spectrum of tiny current density filaments produces a spectrum of magnetic modes that has practically no effect on heat transport, is beneficial.