An official website of the United States government
Abstract Using 3D particle-in-cell simulation, we characterize energy conversion, as a function of guide magnetic field, in a thin current sheet in semirelativistic plasma, with relativistic electrons and subrelativistic protons. There, magnetic reconnection, the drift-kink instability (DKI), and the flux-rope kink instability all compete and interact in their nonlinear stages to convert magnetic energy to plasma energy. We compare fully 3D simulations with 2D in two different planes to isolate reconnection and DKI effects. In zero guide field, these processes yield distinct energy conversion signatures: ions gain more energy than electrons in 2Dxy(reconnection), while the opposite is true in 2Dyz(DKI), and the 3D result falls in between. The flux-rope instability, which occurs only in 3D, allows more magnetic energy to be released than in 2D, but the rate of energy conversion in 3D tends to be lower. Increasing the guide magnetic field strongly suppresses DKI, and in all cases slows and reduces the overall amount of energy conversion; it also favors electron energization through a process by which energy is first stored in the motional electric field of flux ropes before energizing particles. Understanding the evolution of the energy partition thus provides insight into the role of various plasma processes, and is important for modeling radiation from astrophysical sources such as accreting black holes and their jets.
more »
« less
Classical trajectory Monte Carlo simulation of plasma fueling using magnetic plasma expulsion
The possibility of fueling a magnetically confined plasma using particle sources located inside of the plasma is studied by computer simulation. Magnetic plasma expulsion [R. E. Phillips and C. A. Ordonez, Phys. Plasmas 25, 012508 (2018)] would serve to keep the magnetically confined plasma away from the particle sources without adversely affecting plasma confinement. The simulations show how charged particles can be injected into a plasma by using particle sources located directly between two current-carrying wires that create a magnetic expulsion field. Plasma fueling with the average energy of injected particles greater than the average energy of plasma particles may serve for heating the plasma. Also, plasma fueling with positive and negative particles injected at different rates may serve for changing the neutrality of the plasma. Conditions for plasma fueling are investigated using a classical trajectory Monte Carlo simulation. Two types of particle sources are considered, and the fraction of emitted particles that reach (and fuel) the magnetically confined plasma is evaluated for each.
more »
« less
- Award ID(s):
- 1803047
- PAR ID:
- 10597188
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- AIP Advances
- Volume:
- 9
- Issue:
- 7
- ISSN:
- 2158-3226
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT High-energy astrophysical systems frequently contain collision-less relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization (0.04 ≲ σ ≲ 11). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of non-thermal energetic particles (extending the tail by a factor of ∼2). This indicates that non-thermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g. second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, non-thermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*).more » « less
-
Abstract Relativistic magnetic turbulence has been proposed as a process for producing nonthermal particles in high-energy astrophysics. The particle energization may be contributed by both magnetic reconnection and turbulent fluctuations, but their interplay is poorly understood. It has been suggested that during magnetic reconnection the parallel electric field dominates the particle acceleration up to the lower bound of the power-law particle spectrum, but recent studies show that electric fields perpendicular to the magnetic field can play an important, if not dominant role. In this study, we carry out two-dimensional fully kinetic particle-in-cell simulations of magnetically dominated decaying turbulence in a relativistic pair plasma. For a fixed magnetization parameterσ0 = 20, we find that the injection energyεinjconverges with increasing domain size toεinj ≃ 10mec2. In contrast, the power-law index, the cut-off energy, and the power-law extent increase steadily with domain size. We trace a large number of particles and evaluate the contributions of the work done by the parallel (W∥) and perpendicular (W⊥) electric fields during both the injection phase and the postinjection phase. We find that during the injection phase, theW⊥contribution increases with domain size, suggesting that it may eventually dominate injection for a sufficiently large domain. In contrast, on average, both components contribute equally during the postinjection phase, insensitive to the domain size. For high energy (ε ≫ εinj) particles,W⊥dominates the subsequent energization. These findings may improve our understanding of nonthermal particles and their emissions in astrophysical plasmas.more » « less
-
Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.more » « less
-
Comptonization by reconnection plasmoids in black hole coronae I: Magnetically dominated pair plasmaABSTRACT We perform 2D particle-in-cell simulations of reconnection in magnetically dominated electron–positron plasmas subject to strong Compton cooling. We vary the magnetization σ ≫ 1, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting black holes, which produce hard X-rays through Comptonization of seed soft photons. We find that the particle energy spectrum is dominated by a peak at mildly relativistic energies, which results from bulk motions of cooled plasmoids. The peak has a quasi-Maxwellian shape with an effective temperature of ∼100 keV, which depends only weakly on the flow magnetization and the strength of radiative cooling. The mean bulk energy of the reconnected plasma is roughly independent of σ, whereas the variance is larger for higher magnetizations. The spectra also display a high-energy tail, which receives ∼25 per cent of the dissipated reconnection power for σ = 10 and ∼40 per cent for σ = 40. We complement our particle-in-cell studies with a Monte Carlo simulation of the transfer of seed soft photons through the reconnection layer, and find the escaping X-ray spectrum. The simulation demonstrates that Comptonization is dominated by the bulk motions in the chain of Compton-cooled plasmoids and, for σ ∼ 10, yields a spectrum consistent with the typical hard state of accreting black holes.more » « less
