skip to main content


This content will become publicly available on August 1, 2024

Title: Thermal equilibrium of collisional non-neutral plasma in a magnetic dipole trap
This paper discusses thermal equilibrium states of single-species plasmas, such as pure electron plasmas and pure positron plasmas, that are confined in a dipole trap. Thermal equilibrium states for such plasmas are routinely realized in the homogeneous magnetic field of Penning–Malmberg traps. We generalize the theory of these states to include inhomogeneous magnetic dipole fields. The approach to thermal equilibrium takes place in two stages with well separated time scales. On the collision time scale, thermal equilibrium is established along each magnetic field line. On the much longer transport time scale, heat conduction and viscosity bring the plasmas on different flux contours into thermal equilibrium, we call this a state of global thermal equilibrium. We present numerical results for local and global thermal equilibria. These results agree with the analytic predictions for charge collections that are large compared with the Debye length. There is, in principle, no limit to the confinement time of a single-species plasma in a global thermal equilibrium state. Experiments with hot electron–ion plasmas performed in the LDX and RT1 devices give us confidence that, in contrast to a Penning–Malmberg trap, a magnetic dipole field can also confine cold quasi-neutral electron–positron pair plasmas on the time scale of the phenomena of interest. Such pair plasmas are assumed to form in the magnetosphere of neutron stars but have so far not been realized in a laboratory. The creation of an electron–positron pair plasma is the main goal of the APEX collaboration.  more » « less
Award ID(s):
2106332
NSF-PAR ID:
10450436
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Plasma Physics
Volume:
89
Issue:
4
ISSN:
0022-3778
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Inviscid spatial Landau damping is studied experimentally for the case of oscillatory motion of a two-dimensional vortex about its elliptical equilibrium in the presence of an applied strain flow. The experiments are performed using electron plasmas in a Penning–Malmberg trap. They exploit the isomorphism between the two-dimensional Euler equations for an ideal fluid and the drift-Poisson equations for the plasma, where plasma density is the analog of vorticity. Perturbed elliptical vortex states are created using [Formula: see text] strain flows, which are generated by applying voltages to electrodes surrounding the plasma. Measurements of spatial Landau damping (also called critical-layer damping) are in agreement with previous studies in the absence of an applied strain, where the damping is due to a resonance between the local fluid motion and the vortex oscillations. Interestingly, the damping rate does not change significantly over a wide range of applied strain rates. This can be accurately predicted from the initial vorticity profile, even though the resonant frequency is reduced substantially due to the applied strain. For higher amplitude perturbations, nonlinear trapping oscillations also exhibit behavior similar to the strain-free case. In principle, higher-order effects of the applied strain, such as separatrix crossing of peripheral vorticity and interactions with harmonics of the fundamental resonance, are expected to change the damping rate. However, this occurs only for conditions that are not realized in the experiments described here. Vortex-in-cell simulations are used to investigate the possible roles of these effects. 
    more » « less
  2. ABSTRACT

    Black hole and neutron star environments often comprise collisionless plasmas immersed in strong magnetic fields and intense baths of low-frequency radiation. In such conditions, relativistic magnetic reconnection can tap the magnetic field energy, accelerating high-energy particles that rapidly cool by inverse Compton (IC) scattering the dense photon background. At the highest particle energies reached in bright gamma-ray sources, IC scattering can stray into the Klein–Nishina regime. Here, the Comptonized photons exceed pair-production threshold with the radiation background and may thus return their energy to the reconnecting plasma as fresh electron–positron pairs. To reliably characterize observable signatures of such Klein–Nishina reconnection, in this work, we present first-principles particle-in-cell simulations of pair-plasma relativistic reconnection coupled to Klein–Nishina and pair-production physics. The simulations show substantial differences between the observable signatures of Klein–Nishina reconnection and reconnection coupled only to low-energy Thomson IC cooling (without pair production). The latter regime exhibits strong harder-when-brighter behaviour; the former involves a stable spectral shape independent of overall brightness. This spectral stability is reminiscent of flat-spectrum radio quasar (FSRQ) GeV high states, furnishing evidence that Klein–Nishina radiative physics operates in FSRQs. The simulated Klein–Nishina reconnection pair yield spans from low to order-unity and follows an exponential scaling law in a single governing parameter. Pushing this parameter beyond its range studied here might give way to a copious pair-creation regime. Besides FSRQs, we discuss potential applications to accreting black hole X-ray binaries, the M87* magnetosphere, and gamma-ray binaries.

     
    more » « less
  3. Abstract

    We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativistic electrons (and/or positrons) into small-amplitude, kinetic-scale, magnetic-field fluctuations, which pitch-angle scatter the particles and bring the plasma to a near-thermal state of marginal instability. The PIC simulations reveal a nonlinear cyclic evolution of firehose bursts interspersed by periods of stable cooling. We compare the SFHI for electron–positron and electron–ion plasmas. As a byproduct of the growing electron-firehose magnetic-field fluctuations, magnetized ions gain a pressure anisotropy opposite to that of the electrons. If these ions are relativistically hot, we find that they also experience cooling due to collisionless thermal coupling with the electrons, which we argue is mediated by a secondary ion-cyclotron instability. We suggest that the SFHI may be activated in a number of astrophysical scenarios, such as within ejecta from black hole accretion flows and relativistic jets, where the redistribution of energetic electrons from low to high pitch angles may cause transient bursts of radiation.

     
    more » « less
  4. ABSTRACT

    The radiation drag in photon-rich environments of cosmic explosions can seed kinetic instabilities by inducing velocity spreads between relativistically streaming plasma components. Such microturbulence is likely imprinted on the breakout signals of radiation-mediated shocks. However, large-scale, transverse magnetic fields in the deceleration region of the shock transition can suppress the dominant kinetic instabilities by preventing the development of velocity separations between electron–positron pairs and a heavy ion species. We use a 1D five-fluid radiative transfer code to generate self-consistent profiles of the radiation drag force and plasma composition in the deceleration region. For increasing magnetization, our models predict rapidly growing pair multiplicities and a substantial radiative drag developing self-similarly throughout the deceleration region. We extract the critical magnetization parameter σc, determining the limiting magnetic field strength at which a three-species plasma can develop kinetic instabilities before reaching the isotropized downstream. For a relativistic, single ion plasma drifting with γu = 10 in the upstream of a relativistic radiation-mediated shock, we find the threshold σc ≈ 10−7 for the onset of microturbulence. Suppression of plasma instabilities in the case of multi-ion composition would likely require much higher values of σc. Identifying high-energy signatures of microturbulence in shock breakout signals and combining them with the magnetization limits provided in this work will allow a deeper understanding of the magnetic environment of cosmic explosions like supernovae, gamma-ray bursts, and neutron star binary mergers.

     
    more » « less
  5. Abstract

    Statistical and case studies, as well as data‐mining reconstructions suggest that the magnetotail current in the substorm growth phase has a multiscale structure with a thin ion‐scale current sheet embedded into a much thicker sheet. This multiscale structure may be critically important for the tail stability and onset conditions for magnetospheric substorms. The observed thin current sheets are found to be too long to be explained by the models with isotropic plasmas. At the same time, plasma observations reveal only weak field‐aligned anisotropy of the ion species, whereas the anisotropic electron contribution is insufficient to explain the force balance discrepancy. Here we elaborate a self‐consistent equilibrium theory of multiscale current sheets, which differs from conventional isotropic models by weak ion anisotropy outside the sheet and agyrotropy caused by quasi‐adiabatic ion orbits inside the sheet. It is shown that, in spite of weak anisotropy, the current density perturbation may be quite strong and localized on the scale of the figure‐of‐eight ion orbits. The magnetic field, current and plasma density in the limit of weak field‐aligned ion anisotropy and strong current sheet embedding, when the ion scale thin current sheet is nested in a much thicker Harris‐like current sheet, are investigated and presented in an analytical form making it possible to describe the multiscale equilibrium in sharply stretched 2D magnetic field configurations and to use it in kinetic simulations and stability analysis.

     
    more » « less