Search for: All records

Abstract We report on the first self-consistent multidimensional particle-in-cell numerical simulations of nonhomogeneous pair discharges in polar caps of rotation-powered pulsars. By introducing strong inhomogeneities in the initial plasma distribution in our simulations, we analyze the degree of self-consistently emerging coherence of discharges across magnetic field lines. In 2D, we study discharge evolution for a wide range of physical parameters and boundary conditions corresponding to both the absent and free escape of charged particles from the surface of a neutron star. We also present the results of the first 3D simulations of discharges in a polar cap with a distribution of the global magnetospheric current appropriate for a pulsar with 60° inclination angle. For all parameters, we find the coherence scale of pair discharges across magnetic field lines to be of the order of the gap height. We also demonstrate that the popular “spark” model of pair discharges is incompatible with the universally adopted force-free magnetosphere model: intermittent discharges fill the entire zone of the polar cap that allows pair cascades, leaving no space for discharge-free regions. Our findings disprove the key assumption of the spark model about the existence of isolated distinct discharge columns. 

Abstract We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫σ, and form a hard power-law energy distribution,f_i∝γ⁻¹. This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to eitherσin the strong cooling regime or the radiation burnoff limit,γ_syn, when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit,ε_c≈ 16 MeV, is only efficiently produced in the strong cooling regime,γ_syn<σ. In this regime, we find that the spectral cutoff scales asε_cut≈ε_c(σ/γ_syn) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula. 

Abstract Active galactic nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region near the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares near the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations, we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares.