The heat conductivity of a plasma is usually much higher along the magnetic field than across it, and, as a result, the presence of a magnetic island can significantly affect the temperature profile in its vicinity. Radiation energy losses, which depend sensitively on temperature, are thus strongly affected by magnetic islands. This phenomenon is explored in a simple mathematical setting, and it is shown that the presence of a magnetic island greatly enhances a plasma's capacity to radiate energy. In the limit of highly anisotropic heat conductivity, the steady-state heat conduction equation can be reduced to an ordinary differential equation. Although this equation operates in one dimension, the topology is not that of the real line, but corresponds to a rod with a cooling fin. As parameters such as the incoming heat flux or the radiation amplitude are varied, the radiation has a tendency to linger around the island, in particular in the region of the separatrix, and the total radiated energy is then significantly increased. The island acts as a ‘cooling fin’ to the plasma. Furthermore, the solutions exhibit bifurcations, where the location of the radiation zone suddenly changes.

Electron cyclotron resonance ion thrusters (ECRITs) have the potential to be used for space gravitational wave detection due to their wide thrust range. However, an unclear understanding of dynamic processes of ECRITs with strongly coupled multi-operating parameters limits further improvements on thrust noise and response velocity by feedback control systems. An integrative mathematical model considering the non-Maxwell electron energy distribution function for ECRITs is validated by experiments and used to study the influence of operating parameters on the dynamic processes of thrusters, which provides a new simplified grid model. Simulation results show the response processes with microwave (MW) power can be divided into two stages. The characteristic times of the first and second stages are respectively several microseconds and 10 ms, which are respectively dominated by plasma motion and the volume effect. The overshoot of screen grid (SG) current decreases, while its response time remains unchanged when the response time of MW power is prolonged. The response time of SG current with a step increase of flow rate is approximately 10 ms, consistent with the volume effect. The SG current decreases with rise of flow rate for high flow rate operations due to the small increment of ion density limited by low electron temperature, the decrease of ion Bohm velocity and reduction of sheath extraction area. The influence of grid voltage on the dynamic process of the SG current depends on variation ranges of extraction capabilities. When variations of sheath extraction area are limited, the response time is 5 μs, consistent with plasma response time. It is prolonged to 0.5 ms if sheath extraction area variations are large because they cause obvious variations of plasma parameters in the discharge chamber. These dynamic results can not only facilitate designing feedback controllers of micro-propulsion systems for high-precision space missions, but also provide guidance for ion sources to generate highly stable or rapid-response ion beam.

MagNetUS is a network of scientists and research groups that coordinates and advocates for fundamental magnetized plasma research in the USA. Its primary goal is to bring together a broad community of researchers and the experimental and numerical tools they use in order to facilitate the sharing of ideas, resources and common tasks. Discussed here are the motivation and goals for this network and details of its formation, history and structure. An overview of associated experimental facilities and numerical projects is provided, along with examples of scientific topics investigated therein. Finally, a vision for the future of the organization is given.

Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared with unoptimised configurations. Although stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratios and rotational transform values for which stellarators with excellent quasisymmetry on the boundary can be obtained. A large number of Fourier harmonics is included in the boundary representation, which is made computationally tractable by the use of adjoint methods to enable fast gradient-based optimisation and by the direct optimisation of vacuum magnetic fields, which converge more robustly compared with solutions from magnetohydrostatics. Several novel configurations are presented, including stellarators with record levels of quasisymmetry on a surface, three field period quasiaxisymmetric stellarators with substantial magnetic shear, and compact quasisymmetric stellarators at low aspect ratios similar to tokamaks.

A first-order model is derived for quasisymmetric stellarators where the vacuum field due to coils is dominant, but plasma-current-induced terms are not negligible and can contribute to magnetic differential equations, with $\beta$ of the order of the ratio induced to vacuum fields. Under these assumptions, it is proven that the aspect ratio must be large and a simple expression can be obtained for the lowest-order vacuum field. The first-order correction, which involves both vacuum and current-driven fields, is governed by a Grad–Shafranov equation and the requirement that flux surfaces exist. These two equations are not always consistent, and so this model is generally overconstrained, but special solutions exist that satisfy both equations simultaneously. One family of such solutions is the set of first-order near-axis solutions. Thus, the first-order near-axis model is a subset of the model presented here. Several other solutions outside the scope of the near-axis model are also found. A case study comparing one such solution to a VMEC-generated solution shows good agreement.

Physics of the plasma rotation driven by biasing in linear traps is analysed for two limiting cases. The first, relevant for traps with low effective viscosity, considers the line-tying effects to be responsible for the drive as well as for the dissipation of the angular momentum. Meanwhile, in long and thin traps with collisional plasma or developed turbulence, the radial transport of the angular momentum becomes its primary loss channel. The momentum flux goes into the scrape-off layer, which makes conditions there partially responsible for the achievable rotation limits.

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.

New plasma applications in advanced fields, such as fundamental nuclear fusion research, require high-density, large-diameter plasma with a strong and non-uniform magnetic field. A helicon plasma (HP) source using a flat-type antenna is expected to be one of the promising methods for such applications. In this study, we developed an HP source with a two-turn flat-loop antenna connected to a 30 kW radio frequency power supply in the Compact Test Plasma device. In the argon plasma generation experiment with various magnetic fields, HP generation was observed for the first time in this device. The electron density was calculated from the dispersion relation with the magnetic field strength at 45 cm from the antenna surface, assuming a fundamental radial mode and an azimuthal mode of $m=0$. The electron density expected from the experimental result was approximately in the same range as the calculation result by a factor of 2.3 to 3.5. In addition, the magnetic field strength and shape around the antenna are important factors in the plasma properties. This plasma source has been installed in the pilot GAMMA PDX-SC, which is under development for nuclear plasma research, and it contributes to the study of the HP generation process.

The paper presents experimental results from the SMOLA device, which was built in the Budker Institute of Nuclear Physics for the verification of the helical mirror confinement idea. This concept involves active control of axial losses from the confinement zone in an open magnetic trap through the use of multiple mirrors that move in the plasma frame of reference. The discussed experiments focused on determining the cumulative effect of a helical mirror system in combination with a short segment of a stronger magnetic field. Combination of these two methods of axial flow suppression results in higher efficiency compared with each method individually. Different combinations of the mirrors were tested. The most effective flow suppression was observed if the short mirror was placed between the confinement region and the helical mirror. In this configuration, an effective mirror ratio of $R_{{\rm eff}} = 32.6\pm 7.8$ was achieved, along with a more than three-fold increase in plasma density within the confinement region. The possibility of a cumulative effect of different types of magnetic mirrors offers a way to improve the confinement performance of the reactor-grade mirror confinement devices.

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.

Many stellarator coil design problems are plagued by multiple minima, where the locally optimal coil sets can sometimes vary substantially in performance. As a result, solving a coil design problem a single time with a local optimization algorithm is usually insufficient and better optima likely do exist. To address this problem, we propose a global optimization algorithm for the design of stellarator coils and outline how to apply box constraints to the physical positions of the coils. The algorithm has a global exploration phase that searches for interesting regions of design space and is followed by three local optimization algorithms that search in these interesting regions (a ‘global-to-local’ approach). The first local algorithm (phase I), following the globalization phase, is based on near-axis expansions and finds stellarator coils that optimize for quasisymmetry in the neighbourhood of a magnetic axis. The second local algorithm (phase II) takes these coil sets and optimizes them for nested flux surfaces and quasisymmetry on a toroidal volume. The final local algorithm (phase III) polishes these configurations for an accurate approximation of quasisymmetry. Using our global algorithm, we study the trade-off between coil length, aspect ratio, rotational transform and quality of quasi-axisymmetry. The database of stellarators, which comprises approximately 200 000 coil sets, is available online and is called QUASR, for ‘quasi-symmetric stellarator repository’.

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

This paper is based on a report at the 2nd International Fusion Plasma Conference & 13th International Conference on Open Magnetic Systems for Plasma Confinement (iFPC & OS 2023), August 21–25, 2023, Busan, Korea and provides a brief overview of the status of work at the Budker Institute on the study of hot plasma confinement in open-type magnetic traps with a linear axisymmetric configuration. The main attention is paid to key problems: magnetohydrodynamics (MHD) stability in regimes with extremely high relative pressure, longitudinal electronic thermal conductivity, stability with respect to the development of kinetic modes and transverse transport. This paper provides an overview of the methods we are developing to address these problems, the experimental and theoretical results achieved and plans for future development. The last section of the article provides brief information about the preliminary design of the gas-dynamic multiple-mirror trap device, the development of which has been completed.

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.

Influence of the local-ionization-induced neutral depletion on the thrust imparted by the magnetic nozzle plasma thruster is discussed by simply considering reduction of the neutral density due to the ionization in the thruster model combining the global source model and the one-dimensional magnetic nozzle model. When increasing the rf power, it is shown that the increase rate of the plasma density is reduced, while the electron temperature continues to increase due to a decrease in the neutral density. Since the major components of the thrust are originated from the electron pressures in the source and in the magnetic nozzle, the increase in the electron temperature contributes to the increase in the thrust in addition to the gradual density increase by the rf power. The model qualitatively predicts the reduction of the thruster efficiency by the neutral depletion for the high-power condition, compared with the constant neutral density model.

We investigate the effects of different guard geometries on the currents to needle-type Langmuir probes. The results are based on particle-in-cell numerical simulations. We show that if the guard length is less than 6–8 Debye lengths there can be a significant effect on the currents to the probe. A guard radius should not be larger than the Debye length, otherwise it can also significantly affect the currents. However, since guard radii are often close to the probe radius, the second condition is usually satisfied.

Radio frequency quadrupole coolers (RFQCs) are very suitable to cool ion beams with moderate energy spread, typically ions of exotic nuclear species (like $^{132}$Sn$^{1+}$) as in the Selective Production of Exotic Species project at the Laboratori Nazionali di Legnaro, whose ion source supplies 40 keV ions. Beam dynamics includes ion–gas collisions (with a balance of cooling and diffusion effects), acceleration and deceleration and radiofrequency confinement, which can be supplemented by static magnetic field effects. Insertion of a prototype RFQC in the solenoid of the Eltrap machine is also discussed here, with innovations in beam extraction and in modelling, now also based on stochastic equations. Practical consideration on gas pumping and voltage distribution are also included. Typical limits of RFQC are discussed, with special attention to the extracted beam root mean square emittance, which is shown to strongly depend not only on cooler parameters, but also on extraction optics.

In this paper, the modelling and manufacturing of a two-stage supersonic gas jet nozzle enabling the formation of adaptive plasma concentration profiles for injection and acceleration of electrons using few-cycle laser beams are presented. The stages are modelled using the rhoSimpleFoam algorithm of the OpenFOAM computational fluid dynamics software. The first 200–300 ${\rm \mu}$m diameter nozzle stage is dedicated to 1 % N2 + He gas jet formation and electron injection. By varying the pressure between the first and second stages of the injectors, the electron injection location could be adjusted, and the maximum acceleration distance could be ensured. By changing the concentration of the nitrogen in the gas mixture, the charge of the accelerated electrons could be controlled. The second nozzle stage is designed for acceleration in fully ionised He or hydrogen gas and forms the optimal plasma concentration for bubble formation depending on the laser pulse energy, duration and focused beam diameter. In order to reduce the diameter of the plasma profile formed by the first nozzle and the concentration drop gap between the two nozzles, a one-side straight section was introduced in the first nozzle. The shock wave reflected from the straight section of the wall propagates parallel to the shock wave of the intersecting supersonic jets and ensures a minimal gap between the jets. The second-stage longitudinal plasma concentration profile could have an increasing gas density gradient to compensate for dephasing between the electron bunch and the plasma wave due to wave shortening with increasing plasma concentration.

This paper explores the feasibility of a break-even-class mirror referred to as BEAM (break-even axisymmetric mirror): a neutral-beam-heated simple mirror capable of thermonuclear-grade parameters and $Q\sim 1$ conditions. Compared with earlier mirror experiments in the 1980s, BEAM would have: higher-energy neutral beams, a larger and denser plasma at higher magnetic field, both an edge and a core and capabilities to address both magnetohydrodynamic and kinetic stability of the simple mirror in higher-temperature plasmas. Axisymmetry and high-field magnets make this possible at a modest scale enabling a short development time and lower capital cost. Such a $Q\sim 1$ configuration will be useful as a fusion technology development platform, in which tritium handling, materials and blankets can be tested in a real fusion environment, and as a base for development of higher-$Q$ mirrors.

Our earth is immersed in the near-earth space plasma environment, which plays a vital role in protecting our planet against the solar-wind impact and influencing space activities. It is significant to investigate the physical processes dominating the environment, for deepening our scientific understanding of it and improving the ability to forecast the space weather. As a crucial part of the National Major Scientific and Technological Infrastructure–Space Environment Simulation Research Infrastructure (SESRI) in Harbin, the Space Plasma Environment Research Facility (SPERF) builds a system to replicate the near-earth space plasma environment in the laboratory. The system aims to simulate the three-dimensional (3-D) structure and processes of the terrestrial magnetosphere for the first time in the world, providing a unique platform to reveal the physics of the 3-D asymmetric magnetic reconnection relevant to the earth's magnetopause, wave–particle interaction in the earth's radiation belt, particles’ dynamics during the geomagnetic storm, etc. The paper will present the engineering design and construction of the near-earth space plasma simulation system of the SPERF, with a focus on the critical technologies that have been resolved to achieve the scientific goals. Meanwhile, the possible physical issues that can be studied based on the apparatus are sketched briefly. The earth-based system is of great value in understanding the space plasma environment and supporting space exploration.