2025-04-03 | | Total: 3
An Improved Boris algorithm for simulating charged particle motion in electromagnetic fields has been derived. It addresses the issue of inaccurate fast-scale cyclotron phase calculations in the original Boris algorithm, while still maintaining its advantage in simulating slow-scale guiding center motion, hence achieves a balance between low and high-frequency dynamics, overcoming the limitation of traditional second-order volume-preserving algorithms(VPAs) that are constrained to a single characteristic frequency. Test particle simulations indicate that, in most cases, the improved Boris algorithm achieves an exceedingly higher accuracy than conventional VPAs when simulating cases involving various frequencies of electric field within a typical Tokamak magnetic field, highlighting its superior efficacy in handling problems over a large range of characteristic frequencies.
As powerful gamma-ray engines, blazars -- relativistic plasma jets launched toward Earth from active galactic nuclei -- are manifestly high-energy particle accelerators. Yet, exactly how these jets accelerate particles as well as what they are made of both remain largely mysterious. In this work, we argue that these issues may be linked through the gamma-ray emission for which blazars are renowned. Namely, high-energy photons produced at sites of intense particle acceleration could be absorbed by soft radiation within the jet, enriching it with electron-positron pairs. We explore this possibility in the specific context of particle acceleration by magnetized radiative relativistic turbulence. Using a combination of theory, particle-in-cell simulations, and Fokker-Planck modeling, we identify and characterize a novel pair-production-mediated equilibration mechanism in such turbulence. Initially, turbulent energy injection outpaces radiative cooling, leading to runaway particle acceleration and gamma-ray radiation. Then, gamma-ray absorption begets copious newborn pairs, slowing subsequent particle acceleration. This eventually brings particle acceleration into balance with radiative cooling and shuts down pair production: a pair-enriched final equilibrium. We estimate that this process could significantly load jets of flat-spectrum radio quasars with fresh pairs. These results represent an important connection between particle acceleration and plasma composition in blazar jets.
Black hole (BH) accretion disks are often coupled to ultramagnetized and tenuous plasma coronae close to their central BHs. The coronal magnetic field can exchange energy between the disk and the BH, power X-ray emission, and lead to jetted outflows. Up until now, the coronal physics of BH accretion has only been studied using fluid modeling. We construct the first model of a BH feeding on a zero-net-flux accretion disk corona based on kinetic plasma physics. This allows us to self-consistently capture how collisionless relativistic magnetic reconnection regulates the coronal dynamics. We present global, axisymmetric, general relativistic particle-in-cell simulation of a BH coupled, via a series of magnetic loops, to a razor-thin accretion disk. We target the jet-launching regime where the loops are much larger than the BH. We ray-trace high-energy synchrotron lightcurves and track the flow of Poynting flux through the system, including along specific field-line bundles. Reconnection on field lines coupling the BH to the disk dominates the synchrotron output, regulates the flux threading the BH, and ultimately untethers magnetic loops from the disk, ejecting them via a magnetically striped Blandford-Znajek jet. The jet is initially Poynting-dominated, but reconnection operates at all radii, depleting the Poynting power logarithmically in radius. Coronal emission and jet launch are linked through reconnection in our model. This link might explain coincident X-ray flaring and radio-jet ejections observed during hard-to-soft X-ray binary state transitions. It also suggests that striped jet launch could be heralded by a bright coronal counterpart. Our synchrotron signatures resemble variability observed from the peculiar changing-look AGN, 1ES 1927+654, and from Sagittarius A*, hinting that processes similar to our model may be at work in these contexts.