Search for: All records

Abstract We use 3D fully kinetic simulations to study different turbulence modes and turbulence anisotropy of relativistic turbulence in magnetically dominated collisionless plasmas. We extend the method developed by J. Cho & A. Lazarian for decomposing nonrelativistic magnetohydrodynamic (MHD) turbulence into Alfvén, fast, and slow modes to the regime of collisionless plasmas. We find that Alfvén and slow modes are anisotropic, following the P. Goldreich & H. Sridhar scaling, while fast modes are isotropic. We observe a larger kinetic energy fraction of fast modes compared to that found in the nonrelativistic MHD turbulence, suggesting a stronger coupling of Alfvén and fast modes in relativistic magnetized turbulence in collisionless plasmas. We further examine the dynamic alignment and find a weaker scale dependence of the alignment angle than previously proposed. The dominant thermal fluctuations in the kinetic range can cause flattening of the turbulent velocity structure function and weakening of the turbulence anisotropy and dynamic alignment near the kinetic scales. 

Abstract The TeV neutrino signal observed by IceCube from the active galactic nucleus NGC 1068 can probe its innermost coronal regions. If these neutrinos originate from hadrons accelerated within a magnetized turbulent corona, their intensity and spectrum depend on the turbulent magnetic field strength and turbulence coherence scale. The gamma rays accompanying neutrino production are absorbed in this optically thick environment, in a way that depends sensitively on the size of the corona. By a joint fit of the IceCube and Fermi-LAT observations, we translate the multimessenger signal from NGC 1068 and the tentative signal from NGC 7469 into quantitative constraints on coronal properties. NGC 1068, with a significant TeV neutrino excess, favors a compact, strongly magnetized corona with a large turbulence coherence length relative to the coronal size. NGC 7469, with two ∼100 TeV neutrino events, points instead to a somewhat larger corona with much smaller coherence length and high magnetization, but a very small fraction of energy in nonthermal protons. We obtain the diffuse flux from a population of Seyfert galaxies identical to either NGC 1068 or NGC 7469. Finally, we consider a third scenario, motivated by the spectral break observed in the diffuse neutrino flux at tens of TeV, with coronal properties intermediate between the two point-source-inspired models. To enable detailed comparisons with the IceCube and electromagnetic observations, we release our model predictions in a GitHub repository. 

Abstract Coherent structures created through turbulent cascades play a key role in energy dissipation and particle acceleration. In this work, we investigate both current and vorticity sheets in 3D particle-in-cell simulations of decaying relativistic turbulence in pair plasma by training a self-organizing map to recognize these structures. We subsequently carry out an extensive statistical analysis to reveal their geometric and structural properties. This analysis is systematically applied across a range of magnetizations (σ) and fluctuating-to-mean magnetic field strengths (δB₀/B₀) to assess how these parameters influence the resulting structures. We find that the structures’ geometric properties form power-law distributions in their probability density functions, with the exception of the structure width, which generally exhibits an exponential distribution peaking around two electron skin depths. The measurements show a weak dependence onσbut a strong dependence onδB₀/B₀. Finally, we investigate the spatial relationship between current sheets and vorticity sheets. We find that most current sheets are directly associated with at least one vorticity sheet neighbor and are often situated between two vorticity sheets. These findings provide a detailed statistical framework for understanding the formation and organization of coherent structures in relativistic magnetized turbulence, allowing for their incorporation into updated theoretical models for structure-based energy dissipation and particle acceleration processes crucial for interpreting high-energy astrophysical observations. 

Abstract We investigate the onset of driven collisionless reconnection and plasmoid formation in a magnetically dominated pair plasma, using 2D particle-in-cell simulations. Two force-free flux tubes of radiusRare initially pushed together with a prescribed velocity, forming a current sheet whose width shrinks until reconnection sets in. Even in our largest simulation withR ≈ 1600 plasma skin depths, the sheet thickness at reconnection onset is comparable to the skin depth. Plasmoid chains develop when the sheet length-to-width aspect ratioA ≳ 30. In the strongly magnetized limit, the onset of reconnection occurs in roughly 2–6 light-crossing times, depending on the imposed driving timescale, which controls the duration of the thinning phase. In the subsequent nonlinear merging phase, the evolution becomes effectively independent of the initially imposed velocity, leading to magnetic-energy dissipation consistent with a normalized reconnection rate ∼0.1. Our results have important implications for explosive release of magnetic energy in magnetospheres of astrophysical compact objects and their surroundings. 

Abstract We employ first-principles, fully kinetic particle-in-cell simulations to investigate magnetic field-line curvature in magnetically dominated turbulent plasmas and its role in particle acceleration through curvature-drift motion along the motional electric field. By varying the fluctuation-to-mean-magnetic-field ratioδB₀/B₀, we examine curvatureκstatistics and their connection to particle acceleration. The curvature probability densities display broad power-law wings, scaling linearly inκbelow the peak and developing hard high-κtails forδB₀/B₀ ≳ 1. As the mean field strengthens, the high-κtails steepen, and large-curvature events are suppressed whenδB₀/B₀ ≪ 1. The probability density functions of magnetic field-line contraction,v_E · κ, withv_Ethe field-line velocity, develop power-law tails well described by a symmetric Pareto distribution, characteristic of stochastic and intermittent energy exchanges, with the tails becoming harder asδB₀/B₀increases. Our guiding-center analysis shows that curvature-drift acceleration accounts for a substantial fraction of the energization via the motional electric field, and that it strengthens with increasingδB₀/B₀. For well-magnetized particles, curvature-drift acceleration typically exceeds∇Bdrift, polarization drift, and betatron contributions. These results identify curvature-drift acceleration as a principal pathway through which magnetized turbulence transfers energy to nonthermal particles in astrophysical plasmas. 

Abstract Active galactic nuclei (AGN) can accelerate protons to energies of ∼10–100 TeV, with secondary production of high-energy neutrinos. If the acceleration is driven by magnetized turbulence, the main properties of the resulting proton and neutrino spectra can be deduced based on insights from particle-in-cell simulations of magnetized turbulence. We have previously shown that these properties are consistent with the TeV neutrino signal observed from the nearby active galaxy NGC 1068. In this work, we extend this result to a population study. We show that the produced neutrino flux depends mainly on the energetics of the corona—the relative fractions of X-ray, magnetic, and nonthermal proton energy—and on the spectral energy distribution of the AGN. We find that coronae with similar properties can explain neutrinos from the candidate AGN for which IceCube has reported an excess, albeit less significant than NGC 1068. Building on this framework, we show how the neutrino signal evolves with AGN luminosity, and use this AGN sequence to predict the diffuse neutrino flux from the extragalactic population, showing that it can account for the diffuse neutrino signal observed by IceCube in the ∼1–100 TeV energy range. 

Abstract A ubiquitous feature of accreting black hole systems is their hard X-ray emission which is thought to be produced through Comptonization of soft photons by electrons and positrons in the vicinity of the black hole, in a region with optical depth of order unity. The origin and composition of this Comptonizing region, known as the corona, is a matter open for debate. In this paper we investigate the role of relativistic protons accelerated in black-hole magnetospheric current sheets for the pair enrichment and neutrino emission of AGN coronae. Our model has two free parameters, namely the proton plasma magnetizationσ_p, which controls the peak energy of the neutrino spectrum, and the Eddington ratio λ_X,Edd(defined as the ratio between X-ray luminosityL_Xand Eddington luminosityL_Edd), which controls the amount of energy transferred to secondary particles. For sources with λ_X,Edd≳ λ_Edd,crit(where λ_Edd,crit∼ 10^-1forσ_p= 10⁵or ∼ 10^-2forσ_p= 10⁷), proton-photon interactions andγγannihilation produce enough secondary pairs to achieve Thomson optical depthsτ_T∼ 0.1-10. In the opposite case of λ_X,Edd≲ λ_Edd,crit, the coronal pairs cannot originate only from hadronic interactions. Additionally, we find that the neutrino luminosity scales asL²_X/L_Eddfor λ_X,Edd≲ λ_Edd,crit, while it is proportional toL_Xfor higher λ_X,Eddvalues. We apply our model to four Seyfert galaxies, including NGC 1068, and discuss our results in light of recent IceCube observations. 

Abstract The cores of active galactic nuclei are potential accelerators of 10–100 TeV cosmic rays, in turn producing high-energy neutrinos. This picture was confirmed by the compelling evidence of a TeV neutrino signal from the nearby active galaxy NGC 1068, leaving open the question of what is the site and mechanism of cosmic-ray acceleration. One candidate is the magnetized turbulence surrounding the central supermassive black hole. Recent particle-in-cell simulations of magnetized turbulence indicate that stochastic cosmic-ray acceleration is nonresonant, in contrast to the assumptions of previous studies. We show that this has important consequences on a self-consistent theory of neutrino production in the corona, leading to a more rapid cosmic-ray acceleration than previously considered. The turbulent magnetic-field fluctuations needed to explain the neutrino signal are consistent with a magnetically powered corona. We find that strong turbulence, with turbulent magnetic energy density higher than 1% of the rest-mass energy density, naturally explains the normalization of the IceCube neutrino flux, in addition to the neutrino spectral shape. Only a fraction of the protons in the corona, which can be directly inferred from the neutrino signal, are accelerated to high energies. Thus, in this framework, the neutrino signal from NGC 1068 provides a testbed for particle acceleration in magnetized turbulence. 

Abstract Ultra-high-energy cosmic rays (UHECRs), particles characterized by energies exceeding 10¹⁸eV, are generally believed to be accelerated electromagnetically in high-energy astrophysical sources. One promising mechanism of UHECR acceleration is magnetized turbulence. We demonstrate from first principles, using fully kinetic particle-in-cell simulations, that magnetically dominated turbulence accelerates particles on a short timescale, producing a power-law energy distribution with a rigidity-dependent, sharply defined cutoff well approximated by the form $f_{\mathrm{cut}}\left(E,E_{\mathrm{cut}}\right)=\mathrm{sech}\left({\left(E/E_{\mathrm{cut}}\right)}^2\right)$. Particle escape from the turbulent accelerating region is energy dependent, witht_esc∝E^−δandδ∼ 1/3. The resulting particle flux from the accelerator follows $\mathit{dN}/\mathit{dEdt}\propto E^{-s}\mathrm{sech}\left({\left(E/E_{\mathrm{cut}}\right)}^2\right)$, withs∼ 2.1. We fit the Pierre Auger Observatory’s spectrum and composition measurements, taking into account particle interactions between acceleration and detection, and show that the turbulence-associated energy cutoff is well supported by the data, with the best-fitting spectral index being $s={2.1}_{-0.13}^{+0.06}$. Our first-principles results indicate that particle acceleration by magnetically dominated turbulence may constitute the physical mechanism responsible for UHECR acceleration. 

Context.IceCube observations point to active galactic nuclei (AGNs) as promising contributors to the observed astrophysical neutrino flux. Close to the central black hole, protons can be accelerated through magnetic reconnection to very high energies and subsequently interact with abundant X-ray photons in the source, leading to neutrino production. Aims.We investigate whether the diffuse neutrino flux observed by IceCube can originate, via proton acceleration, in reconnection-powered coronae of non-jetted AGNs. Methods.We created a library of neutrino spectral templates over a large grid of values for the three key model parameters: the proton plasma magnetization of the coronaσ_p, the X-ray coronal luminosity, and the black hole mass. Synchrotron cooling of pions and muons plays a significant role due to the large coronal magnetic fields. To infer the diffuse neutrino flux, we coupled the single-source model with a mock AGN catalog consistent with the observed X-ray and mid-infrared AGN samples at redshiftsz = 0 − 4. Results.The coronal emission satisfactorily explains the most recent IceCube measurements of the diffuse neutrino flux up to energies of ∼1 PeV, provided that ∼10% of the AGN coronae haveσ_p ∼ 10⁵, while the rest are distributed over a range of lower magnetizations. Coronal emission is suppressed at higher energies by pion and muon cooling so that another population is required, with jetted AGNs being strong candidates.