Simulation studies of non-neutral plasma equilibria in an electrostatic trap with a magnetic mirror

Abstract

The equilibrium of an infinitely long, strongly magnetized, non-neutral plasma confined in a Penning-Malmberg trap with an additional mirror coil has been solved analytically [J. Fajans, Phys. Plasmas 10, 1209 (2003)] and shown to exhibit unusual features. Particles not only reflect near the mirror in the low field region, but also may be weakly trapped in part of the high field region. The plasma satisfies a Boltzmann distribution along field lines; however, the density and the potential vary along field lines. Some other simplifying assumptions were employed in order to analytically characterize the equilibrium; for example the interface region between the low and high field regions was not considered. The earlier results are confirmed in the present study, where two-dimensional particle-in-cell (PIC) simulations are performed with the Warp code in a more realistic configuration with an arbitrary (but physical) density profile, realistic trap geometry and magnetic field. A range of temperatures and radial plasma sizes are considered. Particle tracking is used to identify populations of trapped and untrapped particles. The present study also shows that it is possible to obtain local equilibria of non-neutral plasmas using a collisionless PIC code, by a scheme that uses the inherent numerical collisionality as a proxy for physical collisions.