A new magnetic mirror machine named KAIMIR (KAIST mirror) has been designed and constructed at the Korea Advanced Institute of Science and Technology (KAIST) to study mirror plasma physics and simulate the boundary regions of magnetic fusion plasmas such as in a tokamak. The purpose of this paper is to introduce the characteristics and initial experimental results of KAIMIR. The cylindrical vacuum chamber has a length of 2.48 m and a diameter of 0.5 m and consists of three sub-chambers, namely the source, centre and expander chambers. A magnetic mirror configuration is achieved by electromagnetic coils with a maximum magnetic field strength of 0.4 T at the mirror nozzles and 0.1 T at the centre. The source plasma is generated by a plasma washer gun installed in the source chamber with a pulse forming network system. The typical discharge time is ~12 ms with a ~6 ms (1–7 ms) steady period. Initial results show that the on-axis electron density at the centre is 1019–20 m−3 and the electron temperature is 4–7 eV. Two parameters were varied in this initial phase, the source power and the mirror ratio, which is the ratio of highest to lowest magnetic field strength in the mirror-confined region. We observed that the increase of the electron density was mitigated for a source power above 0.2 MW. It was also found that the electron density increases almost linearly with the mirror ratio. Accordingly, the stored electron energy was also linearly proportional to the mirror ratio, similar to the scaling of the gas dynamic trap.

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Large-amplitude current-driven plasma instabilities, which can transition to the Buneman instability, were observed in one-dimensional simulations to generate high-energy back-streaming ions. We investigate the saturation of multi-dimensional plasma instabilities and its effects on energetic ion formation. Such ions directly impact spacecraft thruster lifetimes and are associated with magnetic reconnection and cosmic ray inception. An Eulerian Vlasov–Poisson solver employing the grid-based direct kinetic method is used to study the growth and saturation of 2D2V collisionless, electrostatic current-driven instabilities spanning two dimensions each in the configuration (D) and velocity (V) spaces supporting ion and electron phase-space transport. Four stages characterise the electric potential evolution in such instabilities: linear modal growth, harmonic growth, accelerated growth via quasi-linear mechanisms alongside nonlinear fill-in and saturated turbulence. Its transition and isotropisation process bears considerable similarities to the development of hydrodynamic turbulence. While a tendency to isotropy is observed in the plasma waves, followed by electron and then ion phase spaces after several ion-acoustic periods, the formation of energetic back-streaming ions is more limited in the 2D2V than in the 1D1V simulations. Plasma waves formed by two-dimensional electrostatic kinetic instabilities can propagate in the direction perpendicular to the net electron drift. Thus, large-amplitude multi-dimensional waves generate high-energy transverse-streaming ions and eventually limit energetic backward-streaming ions along the longitudinal direction. The multi-dimensional study sheds light on interactions between longitudinal and transverse electrostatic plasma instabilities, as well as fundamental characteristics of the inception and sustenance of unmagnetised plasma turbulence.

The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run times and often to set-ups of less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conserving multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison with linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic time scales.

Overtaking collisions of large-amplitude solitons are investigated via fluid simulations for a plasma consisting of cold ions and Boltzmann-distributed electrons. To achieve this, a new fluid simulation code is presented. In addition, a novel approach to construct soliton initial conditions is developed. Using these ideas, initial conditions are combined that allows the simulation of overtaking collisions. It is shown that, in the small-amplitude regime, simulation results agree well with the two-soliton solution obtained from reductive perturbation theory. Interestingly, in the large amplitude regime, both the slow and fast solitons re-emerge after the collision with no significant change, showing that the collisions remain elastic. A comparison between reductive perturbation analysis and the simulations show that the only significant effect of higher order nonlinearities on overtaking collisions is a reduction in the magnitude of the phase shifts of both solitons.

Recent developments in the design of magnetic confinement fusion devices have allowed the construction of exceptionally optimized stellarator configurations. The near-axis expansion in particular has been proven to enable the construction of magnetic configurations with good confinement properties while taking only a fraction of the usual computation time to generate optimized magnetic equilibria. However, not much is known about the overall features of fast-particle orbits computed in such analytical, yet simplified, equilibria when compared with those originating from accurate equilibrium solutions. This work aims to assess and demonstrate the potential of the near-axis expansion to provide accurate information on particle orbits and to compute loss fractions in moderate to high aspect ratios. The configurations used here are all scaled to fusion-relevant parameters and approximate quasi-symmetry to various degrees. This allows us to understand how deviations from quasi-symmetry affect particle orbits and what are their effects on the estimation of the loss fraction. Guiding-centre trajectories of fusion-born alpha particles are traced using gyronimo and SIMPLE codes under the NEAT framework, showing good numerical agreement. Discrepancies between near-axis and magnetohydrodynamic fields have minor effects on passing particles but significant effects on trapped particles, especially in quasi-helically symmetric magnetic fields. Effective expressions were found for estimating orbit widths and passing–trapped separatrix in quasi-symmetric near-axis fields. Loss fractions agree in the prompt losses regime but diverge afterwards.

The pellet ablation mechanism is an interesting subject for plasma fuelling in fusion plasmas. In GAMMA 10/PDX, pellet injection experiments for higher density plasma production are planned to conduct detached plasma experiments in the higher density plasma condition. We measured the pellet ablation cloud by using the two-directional simultaneous photography system in GAMMA 10/PDX. The tomography reconstruction technique was used for considering the pellet trajectory in the plasma and pellet ablation. The three-dimensional pellet trajectory and pellet ablation images in the plasma were clearly obtained for the first time, to the best of our knowledge.

Collisionless charged particle motion and its transport in the two-wire model (TWM) with no axial magnetic fields is investigated numerically. The TWM configuration contains a magnetic X-point, and single particle motions in such a field have two conserved quantities: the total kinetic energy and the base field line value which is a quantity derived from the axial canonical momentum. As gyrating particles travel along the field lines, they may reach near the X-point region where the magnetic moment, the first adiabatic invariant, can be occasionally shifted due to a large gradient of the field. When the magnetic moment becomes large, resulting in a large Larmor radius, particles probabilistically cross the X-point to migrate to the opposite side of the TWM configuration. These phenomena are investigated with single particle simulations. We find that the statistical behaviour of the seemingly chaotic magnetic moment shifts are completely determined by the two aforementioned conserved quantities, and also that there exists a threshold energy, determined by the base field line value, allowing only particles with a higher energy to cross the separatrix and migrate. It is found that the crossing time is distributed exponentially, and that the migration confinement time, which is the average crossing time, is shorter for particles with a base field line closer to the separatrix and a higher energy. We provide an empirical expression, derived with the simulations, for estimating the collisionless migration confinement time.

The excitation conditions of the magnetorotational instability (MRI) are studied for axially unbounded Taylor–Couette (TC) flows of various gap widths between the cylinders. The cylinders are considered as made from both perfect-conducting or insulating material and the conducting fluid with a finite but small magnetic Prandtl number rotates with a quasi-Keplerian velocity profile. The solutions are optimized with respect to the wavenumber and the Reynolds number of the rotation of the inner cylinder. For the axisymmetric modes, we find the critical Lundquist number of the applied axial magnetic field: the lower, the wider the gap between the cylinders. A similar result is obtained for the induced cell structure: the wider the gap, the more spherical the cells are. The marginal rotation rate of the inner cylinder – for a fixed size of the outer cylinder – always possesses a minimum for not too wide and not too narrow gap widths. For perfect-conducting walls the minimum lies at $r_{{\rm in}}\simeq 0.4$, where $r_{{\rm in}}$ is the ratio of the radii of the two rotating cylinders. The lowest magnetic field amplitudes to excite the instability are required for TC flows between perfect-conducting cylinders with gaps corresponding to $r_{{\rm in}}\simeq ~0.2$. For even wider and also for very thin gaps the needed magnetic fields and rotation frequencies are shown to become rather huge. Also the non-axisymmetric modes with $|m|=1$ have been considered. Their excitation generally requires stronger magnetic fields and higher magnetic Reynolds numbers in comparison with those for the axisymmetric modes. If TC experiments with too slow rotation for the applied magnetic fields yield unstable modes of any azimuthal symmetry, such as the currently reported Princeton experiment (Wang et al., Phys. Rev. Lett., vol. 129, 115001), then also other players, including axial boundary effects, than the MRI-typical linear combination of current-free fields and differential rotation should be in the game.

When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter, we derive and implement kinetic sources for the Compton scattering and tritium beta decay runaway electron generation mechanisms in our simulation tool Dream (Hoppe et al., Comput. Phys. Commun., vol. 268, 2021, 108098). To achieve a diverse set of disruption scenarios, Bayesian optimization is used to explore a range of massive material injection densities for deuterium and neon. The cost function is designed to distinguish between successful and unsuccessful disruption mitigation based on the runaway current, current quench time and transported fraction of the heat loss. In the non-activated scenarios, we find that fluid and kinetic disruption simulations can have significantly different runaway electron dynamics, due to an overestimation of the runaway seed by the fluid model. The primary cause of this is that the fluid hot-tail generation model neglects superthermal electron transport losses during the thermal quench. In the activated scenarios, the fluid and kinetic models give similar predictions, which can be explained by the significant influence of the activated sources on the runaway dynamics and the seed.

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul (Phys. Rev. Lett., vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$. An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$. Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. (Phys. Plasma, vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$.

We analysed the excitation of a surface magnetoplasmon wave by the mode conversion of a p-polarized laser beam over a rippled semiconductor (n-type)-free space interface. The pump surface magnetoplasmon wave exerts a ponderomotive force on the free electrons in the semiconductor, imparting a linear oscillatory velocity at the laser modulation frequency to them. This linear oscillatory velocity couples with the modulated electron density to produce a current density, which develops a resonant surface magnetoplasmon wave in the terahertz region. The amplitude of the terahertz surface magnetoplasmon wave can be tuneable with an external magnetic field and the semiconductor's temperature.

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

Laser-driven neutron sources are routinely produced by the interaction of laser-accelerated protons with a converter. They present complementary characteristics to those of conventional accelerator-based neutron sources (e.g. short pulse durations, enabling novel applications like radiography). We present here results from an experiment aimed at performing a global characterization of the neutrons produced using the Titan laser at the Jupiter Laser Facility (Livermore, USA), where protons were accelerated from 23 $\mathrm {\mu }$m thick plastic targets and directed onto a LiF converter to produce neutrons. For this purpose, several diagnostics were used to measure these neutron emissions, such as CR-39, activation foils, time-of-flight detectors and direct measurement of $^7{\rm Be}$ residual activity in the LiF converters. The use of these different, independently operating diagnostics enables comparison of the various measurements performed to provide a robust characterization. These measurements led to a neutron yield of $2.0\times 10^{9}$ neutrons per shot with a modest angular dependence, close to that simulated.

Spontaneous avalanche to plasma begins in the core of an ellipsoidal Rydberg gas of nitric oxide. Ambipolar expansion of NO$^+$ draws energy from avalanche-heated electrons. Then, cycles of long-range resonant electron transfer from Rydberg molecules to ions equalize their relative velocities. This sequence of steps gives rise to a remarkable mechanics of self-assembly, in which the kinetic energy of initially formed hot electrons and ions drives an observed separation of plasma volumes. These dynamics adiabatically sequester energy in a reservoir of mass transport, starting a process that anneals separating volumes to form an apparent glass of strongly coupled ions and electrons. Short-time electron spectroscopy provides experimental evidence for complete ionization. The long lifetime of this system, particularly its stability with respect to recombination and neutral dissociation, suggests that this transformation affords a robust state of arrested relaxation, far from thermal equilibrium. We see this most directly in the excitation spectrum of transitions to states in the initially selected Rydberg series, detected as the long-lived signal that survives a flight time of $500\ \mathrm {\mu }$s to reach an imaging detector. The initial density of electrons produced by prompt Penning ionization, which varies with the selected initial principal quantum number and density of the Rydberg gas, determines a balance between the rising density of ions and the falling density of Rydberg molecules. This Penning-regulated ion-Rydberg molecule balance appears necessary as a critical factor in achieving the long ultracold plasma lifetime to produce spectral features detected after very long delays.

A unique field-reversed configuration (FRC) experiment is presently being assembled at the Plasma Technology Research Institute, KFE. It is a compact small-scale FRC device, which uses a set of radio frequency (RF) antennas to produce an internal E × B that drives the electrons for current-drive, in which E is the electric field and B is the magnetic field. This is very similar to the rotating magnetic field (RMF) current-drive, where the horizontal and vertical antennas are driven 90° out of phase. For this device, the RF antennas are arranged differently than the RMF. The RF antennas, being two separate sets, are positioned inside the vacuum chamber. Each set consists of 8 coils, for a total of 16 coils, where 80~100 kHz sine and cosine waveform currents are applied. One set of coils generates a radial B-field, while the other set provides an E-field in the z-direction. As the phase changes, the E and B fields are switched by these two sets. Nevertheless, E × B propagates in the same θ-direction so that this allows the electrons to rotate around the circumference of the device. The FRC device will test wave heating by launching 2.45 GHz microwaves. Also, passive stabilizers are positioned at each end to provide extra stability while preventing tilt instability. The experiment is expected to produce its first plasma in 2025.

The extended guiding-centre Lagrangian equations of motion are derived by the Lie-transform perturbation method under the assumption of time-dependent and inhomogeneous electric and magnetic fields that satisfy the standard guiding-centre space–time orderings. Polarization effects are introduced into the Lagrangian dynamics by the inclusion of the polarization drift velocity in the guiding-centre velocity and the appearance of finite-Larmor-radius corrections in the guiding-centre Hamiltonian and guiding-centre Poisson bracket.

The two-time energy spectrum of weak magnetohydrodynamic turbulence is found by applying a wave-turbulence closure to the cumulant hierarchy constructed from the dynamical equations. Solutions are facilitated via asymptotic expansions in terms of the small parameter $\varepsilon$, describing the ratio of time scales corresponding to Alfvénic propagation and nonlinear interactions between counter-propagating Alfvén waves. The strength of nonlinearity at a given spatial scale is further quantified by an integration over all possible delta-correlated modes compliant in a given set of three-wave interactions that are associated with energy flux through the said scale. The wave-turbulence closure for the two-time spectrum uncovers a secularity occurring on a time scale of order $\varepsilon ^{-2}$, and the asymptotic expansion for the spectrum is reordered in a manner comparable to the one-time case. It is shown that for the regime of stationary turbulence, the two-time energy spectrum exponentially decays on a lagged time scale $(\varepsilon ^2 \gamma _k^s)^{-1}$ in proportion to the strength of the associated three-wave interactions, characterized by nonlinear decorrelation frequency $\gamma _k^s$. The scaling of the form $k_{\perp } v_0 \chi _0$ exhibited by this frequency is reminiscent of random sweeping by the outer scale with characteristic fluctuation velocity $v_0$ that is modified due to competition with Alfvénic propagation (characterized by $\chi _0$) at the said scale. A brief calculation of frequency broadening of the power spectrum due to nonlinear interactions is also presented.

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC. Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.