We perform experiments in the Large Plasma Device at the University of California, Los Angeles, studying how different end-electrode biasing schemes modify the radial potential profile in the machine. We impose biasing profiles of different polarities and gradient signs on a set of five concentric electrodes placed 12 m downstream from the plasma source. We find that imposing concave-down profiles (negative potential radial gradient) on the electrodes creates radial potential profiles halfway up the plasma column that are comparable to those imposed on the electrodes and a few electron temperature in height, regardless of the biasing polarity. On the other hand, imposing concave-up profiles (positive potential radial gradient) leads to non-monotonic radial potential profiles. This observation can be explained by the current drawn through the electrodes and the parallel plasma resistivity, highlighting their important role in controlling the rotation of plasma. Concave-down plasma potential profiles, obtained by drawing electrons on the axis, are predicted to drive azimuthal drift velocities that can approach significant fractions of the ion sound speed in the central region of the plasma column.

During a space mission, switching to an electric propulsion system from chemical propulsion, once the spacecraft is out of the Earth's gravity, significantly reduces the mission's overall cost. In electric propulsion, the Hall thruster and gridded ion thruster are established technologies. These thrusters compromise mission longevity due to continuous erosion of the device electrode material. To overcome this issue, an electrode less expanding magnetic field plasma thruster or helicon plasma thruster (HPT), was proposed and research is on going worldwide. The HPT shows scaling of thrust with input power while Hall thrusters and ion thrusters do not. Typically, an inert xenon gas is used as a fuel in HPT devices due to a low ionization potential and non-hazardous nature. Xenon is not easily available in nature and during a space mission it needs to be stored in high pressure tanks. Recently, iodine has been proposed as an alternate to xenon as it is easily available and does not have any storage issues. In most of the numerical simulations, argon is used as a fuel gas to reduce the simulation cost. Using a 1D3V particle-in-cell Monte Carlo collision code, we present here a net thrust generation for different fuel gases such as argon, xenon and iodine. We compare plasma flow rates and directed ion beam velocity for different fuel gases having identical inputs. Thrust and plasma flow are investigated for different magnetic field gradients in the plasma expansion region for unidirectional and bidirectional HPT and is reported here. Using iodine fuel, a significant increase in net thrust is obtained for higher magnetic field divergence for identical simulation input parameters while comparing with xenon fueled cases.

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.

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.

The main emphasis of the Laser Ion Generation, Handling and Transport (LIGHT) beamline at GSI Helmholtzzentrum für Schwerionenforschung GmbH are phase-space manipulations of laser-generated ion beams. In recent years, the LIGHT collaboration has successfully generated and focused intense proton bunches with an energy of 8 MeV and a temporal duration shorter than 1 ns (FWHM). An interesting area of application that exploits the short ion bunch properties of LIGHT is the study of ion-stopping power in plasmas, a key process in inertial confinement fusion for understanding energy deposition in dense plasmas. The most challenging regime is found when the projectile velocity closely approaches the thermal plasma electron velocity ($v_{i}\approx v_{e,\text {th}}$), for which existing theories show high discrepancies. Since conclusive experimental data are scarce in this regime, we plan to conduct experiments on laser-generated plasma probed with ions generated with LIGHT at a higher temporal resolution than previously achievable. The high temporal resolution is important because the parameters of laser-generated plasmas are changing on the nanosecond time scale. To meet this goal, our recent studies have dealt with ions of lower kinetic energies. In 2021, laser accelerated carbon ions were transported with two solenoids and focused temporally with LIGHT's radio frequency cavity. A bunch length of 1.2 ns (FWHM) at an energy of 0.6 MeV u$^{-1}$ was achieved. In 2022, protons with an energy of 0.6 MeV were transported and temporally compressed to a bunch length of 0.8 ns. The proton beam was used to measure the energy loss in a cold foil. Both the ion and proton beams will also be employed for energy loss measurements in a plasma target.

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.

The ASACUSA (atomic spectroscopy and collisions using slow antiprotons) Cusp experiment requires the production of dense positron plasmas with a high repetition rate to produce a beam of antihydrogen. In this work, details of the positron production apparatus used for the first observation of the antihydrogen beam, and subsequent measurements, are described in detail. This apparatus replaced the previous compact trap design resulting in an improvement in the positron accumulation rate by a factor of $52\pm 3$.

A high-intensity laser pulse propagating through a gas target disturbs the uniform plasma distribution. Plasma density structures, created by high-order Bessel–Gauss beams for guiding the accelerating Gaussian beam and laser wakefield acceleration of electrons, are analysed using Wake-T and Fourier–Bessel particle-in-cell (FBPIC) simulation tools. The use of Bessel–Gauss doughnut beams increases the acceleration distance and energy of accelerated electrons up to 2.3 times at a 2 mm distance relative to the Gaussian beam of the same intensity.

A novel curved surface nanowire target is proposed to improve the cutoff energy of accelerated protons via target normal sheath acceleration. The interaction of a laser of intensity $1.37\times 10^{20}\ {\rm W}\ {\rm cm}^{-2}$ with a curved surface nanowire target is studied by two-dimensional particle-in-cell simulations. The numerical results indicate that the sheath electric field at the target rear side is significantly enhanced by this simple target design, compared with using the planar nanowire target. The transverse motion of hot electrons is effectively confined and the energy density of electrons is naturally increased. A series of simulations with various target parameters is carried out to investigate the performance of this novel target. This tailored target may provide implications for generating high-quality proton beams in experiments.

Plasma jets are widely investigated both in the laboratory and in nature. Astrophysical objects such as black holes, active galactic nuclei and young stellar objects commonly emit plasma jets in various forms. With the availability of data from plasma jet experiments resembling astrophysical plasma jets, classification of such data would potentially aid in not only investigating the underlying physics of the experiments but also the study of astrophysical jets. In this work we use deep learning to process all of the laboratory plasma images from the Caltech Spheromak Experiment spanning two decades. We found that cosine similarity can aid in feature selection, classify images through comparison of feature vector direction and be used as a loss function for the training of AlexNet for plasma image classification. We also develop a simple vector direction comparison algorithm for binary and multi-class classification. Using our algorithm we demonstrate 93 % accurate binary classification to distinguish unstable columns from stable columns and 92 % accurate five-way classification of a small, labelled data set which includes three classes corresponding to varying levels of kink instability.

This paper focuses on utilizing several different optical diagnostics to experimentally characterize a pure helium atmospheric pressure plasma jet. Axial electric field measurements were carried out along the plasma plume through the use of a non-perturbing method based on polarization-dependent Stark spectroscopy of the helium $492.2$ nm line. The electric field is shown to increase with distance along the plume length, reaching values as high as $24.5$ kV cm$^{-1}$. The rate of increase of the electric field is dependent on the helium gas flow rate, with lower gas flows rising quicker with distance in comparison with larger flow rates, with the typical values remaining within the same range. This sensitivity is linked to gas mixing between the helium and surrounding ambient air. Schlieren imaging of the gas flow is included to support this. The addition of a target is shown to further increase the measured electric field in close range to the target, with the magnitude of this increase being strongly dependent on the distance between the tube exit and target. The relative increase in the electric field is shown to be on average greater for a conducting target of water in comparison with plastic. A minimal equipment optical configuration, which is here referred to as fast two-dimensional monochromatic imaging, is introduced as an approach for estimating excited state densities within the plasma. Densities of the upper helium states for transitions, $1s3s$ $^{3}S_{1}$ $\rightarrow$ $1s2p$ $^{3}P^{0}_{0,1,2}$ at $706.5$ nm and $1s3s$ $^{1}S_{0}$ $\rightarrow$ $1s2p$ $^{1}P^{0}_{1}$ at $728.1$ nm, were estimated using this approach and found to be of the order of $10^{10}$–$10^{11}$ cm$^{-3}$.

We report a generation of energetic protons by the interaction of a high-energy electron driving beam with an underdense plasma slab. After an interaction period of approximately 4000 fs, a proton beam with maximum energy greater than 250 MeV can be achieved by applying a driving beam with energy 1.0 GeV to a 200 $\mathrm {\mu }$m plasma slab. Our two-dimensional particle-in-cell simulations also show that the proton acceleration process can be divided into two stages. In the first stage, a strong positive longitudinal electric field appears near the rear boundary of the plasma slab after the driving beam has passed through it. This acceleration process is similar to the target normal sheath acceleration scheme by the interaction between intense pulsed laser with overdense plasma targets. In the second stage, the accelerated protons experience a long-range acceleration process with a two-stream instability between the high-energy driving beam and the proton beam. Further analyses show that this accelerated proton beam is equipped with the property of good collimation and high energy. This scheme presents a new way for proton or ion acceleration on some special occasions.

In the construction of a stellarator, the manufacturing and assembling of the coil system is a dominant cost. These coils need to satisfy strict engineering tolerances, and if those are not met the project could be cancelled as in the case of the National Compact Stellarator Experiment (NCSX) project (R.L. Orbach, 2008, https://ncsx.pppl.gov/DOE_NCSX_052208.pdf). Therefore, our goal is to find coil configurations that increase construction tolerances without compromising the performance of the magnetic field. In this paper, we develop a gradient-based stochastic optimization model which seeks robust stellarator coil configurations in high dimensions. In particular, we design a two-step method: first, we perform an approximate global search by a sample efficient trust-region Bayesian optimization; second, we refine the minima found in step one with a stochastic local optimizer. To this end, we introduce two stochastic local optimizers: BFGS applied to the sample average approximation; and Adam, equipped with a control variate for variance reduction. Numerical simulations performed on a W7-X-like coil configuration demonstrate that our global optimization approach finds a variety of promising local solutions at less than $0.1\,\%$ of the cost of previous work, which considered solely local stochastic optimization.

A spacecraft during mission typically switches from chemical propulsion to electric propulsion once it lifts out of the Earth's gravity as the thrust requirement to drive it reduces substantially. Consequently, electric propulsion technology is commonly used for deep space mission, satellite orbit keeping and orbit correction. In the last few decades, the amount of man-made space junk (space debris) has increased enormously and has become a potential danger for space stations, space shuttles and other live satellites. A bi-directional plasma thruster, mounted on a satellite, can be used to remove space debris during satellite operation (Takahashi et al., Sci. Rep., vol. 8, 2018, p. 14417). A directed ion beam ejected from a plasma thruster imparts a net force on space debris to decelerate and facilitates manoeuvring and re-entry of space debris into the Earth's atmosphere where it can burn out. We present a detailed 1D3V PIC-MCC (particle in cell-Monte Carlo collision) simulation of a bi- directional plasma thruster. To this end, a PIC-MCC solver which resolves thruster axial direction and all three velocity dimensions is used to study a magnetic nozzle plasma thruster with both ends open (bi-directional plasma thruster). We show that such a bi-directional plasma thruster can be used to accelerate–decelerate a live satellite and also to remove space debris by altering the magnetic field spatial profile in the plasma expansion region. A detailed study is presented.

Particle acceleration via magnetic reconnection is a fundamental process in astrophysical plasmas. Experimental architectures are able to confirm a wide variety of particle dynamics following the two-dimensional Sweet–Parker model, but are limited in their reproduction of the fan-spine magnetic field topology about three-dimensional (3-D) null points. Specifically, there is not yet an experiment featuring driven 3-D torsional magnetic reconnection. To move in this direction, this paper expands on recent work toward the design of an experimental infrastructure for inducing 3-D torsional fan reconnection by predicting feasible particle acceleration profiles. Solutions to the steady-state, kinematic, resistive magnetohydrodynamic equations are used to numerically calculate particle trajectories from a localized resistivity profile using well-understood laboratory plasma parameters. We confine a thin, 10 eV helium sheath following the snowplough model into the region of this localized resistivity and find that it is accelerated to energies of ${\approx }2$ keV. This sheath is rapidly accelerated and focused along the spine axis propagating a few centimetres from the reconnection region. These dynamics suggest a novel architecture that may hold promise for future experiments studying solar coronal particle acceleration and for technology applications such as spacecraft propulsion.

The possibility of controlling the electrostatic field distribution in plasma has yielded wide prospects for modern technologies. As a magnetic field primarily allows for creating electric fields in plasma, it serves as an additional obstacle for the current flow through a medium. In the present paper, an axially symmetric system is considered in which the magnetic field is directed along the axis and concentric electrodes are located at the ends. The electrodes are negatively biased. A model which solves the problem of the radial distribution of the plasma potential inside the cylindrical plasma column supported by the end electrodes is proposed. The most commonly encountered configurations of the electrical connection for the end electrodes are considered, and the particular solutions to the problem of the radial distribution are presented. The contribution of ions and electrons to the transverse conductivity is evaluated in detail. The influence of a thermionic element on the radial profile of the plasma potential is considered. To verify the proposed model, an experimental study of the reflex discharge is carried out with both cold electrodes and a thermionic element on the axis. A comparison of the computational model results with experimental data is given. The presented model makes it possible to solve the problem concerning the plasma potential distribution in the case of an arbitrary number of end electrodes, and also to take into account the inhomogeneity of the distribution of plasma density, neutral gas pressure and electron temperature along the radius.

The purpose of the present paper is to investigate the generation of soft X-ray emission from an anharmonic collisional nanoplasma by a laser–nanocluster interaction. The electric field of the laser beam interacts with the nanocluster and leads to ionization of the cluster atoms, which then produces a nanoplasma. Because of the nonlinear restoring force in an anharmonic nanoplasma, the fluctuations and heating rate of, as well as the power radiated by, the electrons in the nanocluster plasma will be notably different from those arising from a linear restoring force. By comparing the nonlinear restoring force state (which arises from an anharmonic cluster) with that of the linear restoring force (in harmonic clusters), the cluster temperature specifically changes at the resonant frequency relative to the linear restoring force, while the variation of the anharmonic cluster radius is almost identical to that of the harmonic cluster radius. In addition, it is revealed that a sharp peak of X-ray emission arises after some picoseconds in deuterium, helium, neon and argon clusters.

In this article, the electron trapping and acceleration in the wake field driven by an ultrarelativistic hollow electron beam is studied. When the hollow driver injects into plasma, there is a doughnut-shaped electron bubble formed because of the existence of a special ‘backflow’ beam in the centre of the electron bubble. At the same time, there is a transverse convergence of the hollow driver, which leads to the weakening of the backflow beam. This results in a local electron density transition at the rear of the bubble. During this process, there is an expansion of the longitudinal electron bubble size, and a bunch of background electrons is trapped by the wake field at the rear of the bubble. The tracks for the trapped electrons show that there are two sources: one is from the bubble sheath and the other is from the unique backflow beam. In the particle-in-cell simulation where the driving beam has initial energy of $1.0$ GeV per particle, the trapped beam can be accelerated to energy of more than $1.5$ GeV per particle and the corresponding transformer ratio is $1.5$. With the increase of driving beam energy up to $40.0$ GeV, a transformer ratio of $1.4$ still can be achieved. By adjusting the hollow beam density, it is possible to control the trapped beam charge value and beam quality, such as its energy spread and transverse emittance.

Following up on earlier work which demonstrated an improved numerical stellarator coil design optimization performance by the use of stochastic optimization (Lobsien et al., Nucl. Fusion, vol. 58 (10), 2018, 106013), it is demonstrated here that significant further improvements can be made – lower field errors and improved robustness – for a Wendelstein 7-X test case. This is done by increasing the sample size and applying fully three-dimensional perturbations, but most importantly, by changing the design sequence in which the optimization targets are applied: optimization for field error is conducted first, with coil shape penalties only added to the objective function at a later step in the design process. A robust, feasible coil configuration with a local maximum field error of 3.66 % and an average field error of 0.95 % is achieved here, as compared to a maximum local field error of 6.08 % and average field error of 1.56 % found in our earlier work. These new results are compared to those found without stochastic optimization using the FOCUS and ONSET suites. The relationship between local minima in the optimization space and coil shape penalties is also discussed.

This paper investigates the motion of three-dimensional ideal magnetohydrodynamics with incompressible flows. The governing equation is performed at steady state, with the magnetic field parallel to the plasma flow. The equations of stationary equilibrium are derived and described mathematically in Cartesian space. Two approaches for derivation of general three-dimensional solutions for Alfvénic and non-Alfvénic flows at constant and variable fluid densities are constructed. The general vector and scalar potentials of the velocity field are used to derive general formulas of general three-dimensional solutions for Alfvénic and non-Alfvénic flows. To verify the general results we have obtained, some examples are presented. An application that may be of interest for coronal loops and solar prominences is presented.