Andréas Sundström

In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (MHD) coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ∼ 1 M A compared to ∼ 8.7 M A of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r / a < 0.3 and with kinetic energies ∼ 0.2 – 15 M e V , contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ∼ 18 m 2 s − 1 , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field’s effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.

We present a scheme for amplifying an extreme-ultraviolet (XUV) seed isolated attosecond pulse via stimulated Raman scattering of a pulse-train pump. At sufficient seed and pump intensity, the amplification is nonlinear, and the amplitude of the seed pulse can reach that of the pump, one order of magnitude higher than the initial seed amplitude. In the linear amplification regime, we find that the spectral signature of the pump pulse train is imprinted on the spectrum of the amplified seed pulse. Since the spectral signature is imprinted with its frequency downshifted by the plasma frequency, it is possible to deduce the electron density in the region of interaction. This region can be of micrometer length scale longitudinally. By varying the delay between the seed and the pump, this scheme provides a local electron-density measurement inside solid-density plasmas that cannot be probed with optical frequencies, with micrometer resolution.

The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.

Extreme-ultraviolet pulses can propagate through ionised solid-density targets, unlike optical pulses and, thus, have the potential to probe the interior of such plasmas on sub-femtosecond timescales. We present a synthetic diagnostic method for solid-density laser-generated plasmas based on the dispersion of an extreme-ultraviolet attosecond probe pulse, in a pump–probe scheme. We demonstrate the theoretical feasibility of this approach through calculating the dispersion of an extreme-ultraviolet probe pulse propagating through a laser-generated plasma. The plasma dynamics is calculated using a particle-in-cell simulation, whereas the dispersion of the probe is calculated with an external pseudo-spectral wave solver, allowing for high accuracy when calculating the dispersion. The application of this method is illustrated on thin-film plastic and aluminium targets irradiated by a high-intensity pump pulse. By comparing the dispersion of the probe pulse at different delays relative to the pump pulse, it is possible to follow the evolution of the plasma as it disintegrates. The high-frequency end of the dispersion provides information on the line-integrated electron density, whereas lower frequencies are more affected by the highest density encountered along the path of the probe. In addition, the presence of thin-film interference could be used to study the evolution of the plasma surface.

The potential formation of multi-mega-ampere beams of relativistic ‘runaway’ electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak’s development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium–tritium plasmas.

We perform fully kinetic simulations of flows known to produce dynamo in magnetohydrodynamics (MHD), considering scenarios with low Reynolds number and high magnetic Prandtl number, relevant for galaxy cluster scale fluctuation dynamos. We find that Landau damping on the electrons leads to a rapid decay of magnetic perturbations, impeding the dynamo. This collisionless damping process operates on spatial scales where electrons are nonmagnetized, reducing the range of scales where the magnetic field grows in high magnetic Prandtl number fluctuation dynamos. When electrons are not magnetized down to the resistive scale, the magnetic energy spectrum is expected to be limited by the scale corresponding to magnetic Landau damping or, if smaller, the electron gyroradius scale, instead of the resistive scale. In simulations we thus observe decaying magnetic fields where resistive MHD would predict a dynamo.

We numerically investigate the impact of Coulomb collisions on the ion dynamics in high-Z, solid density caesium hydride and copper targets, irradiated by high-intensity ( I ≈ 2 − 5 × 10 20   W   c m − 2 ), ultrashort (∼10 fs), circularly polarized laser pulses, using particle-in-cell simulations. Collisions significantly enhance electron heating, thereby strongly increasing the speed of a shock wave launched in the laser-plasma interaction. In the caesium hydride target, collisions between the two ion species heat the protons to ∼100−1000 eV temperatures. However, in contrast to previous work (A E Turrell et al 2015 Nat. Commun. 6 8905), this process happens in the upstream only, due to nearly total proton reflection. This difference is ascribed to distinct models used to treat collisions in dense/cold plasmas. In the case of a copper target, ion reflection can start as a self-amplifying process, bootstrapping itself. Afterwards, collisions between the reflected and upstream ions heat these two populations significantly. When increasing the pulse duration to 60 fs, the shock front more clearly decouples from the laser piston, and so can be studied without direct interference from the laser. The shock wave formed at early times exhibits properties typical of both hydrodynamic and electrostatic shocks, including ion reflection. At late times, the shock is seen to evolve into a hydrodynamic blast wave.

In strictly collisionless electrostatic shocks, the ion distribution function can develop discontinuities along phase-space separatrices, due to partial reflection of the ion population. In this paper, we depart from the strictly collisionless regime and present a semi-analytical model for weakly collisional kinetic shocks. The model is used to study the effect of small but finite collisionalities on electrostatic shocks, and they are found to smooth out these discontinuities into growing boundary layers. More importantly, ions diffuse into and accumulate in the previously empty regions of phase space, and, by upsetting the charge balance, lead to growing downstream oscillations of the electrostatic potential. We find that the collisional age of the shock is the more relevant measure of the collisional effects than the collisionality, where the former can become significant during the lifetime of the shock, even for weak collisionalities.

We describe a low-cost yet experimentally challenging method to measure the acceleration of gravity, g, using a liquid in a rotating bowl and a laser pointer. The idea underpinning this novel method is that the rotating liquid surface will form a parabolic reflector which will focus light into a unique focal point. By measuring the height of the focal point, g could be determined to m s⁻². We discuss the pedagogical merits of this method compared to more traditional methods for measuring g, and how it can be implemented as an experimental problem at different educational levels.