Publications by Year: 2024

2024
Rodrigues X, Karl M, Padovani P, Giommi P, Paiano S, Falomo R, Petropoulou M, Oikonomou F. The Spectra of IceCube Neutrino (SIN) candidate sources: V. Modeling and interpretation of multiwavelength and neutrino data. [Internet]. 2024;689:A147. WebsiteAbstract
Context. A correlation has been reported between the arrival directions of high-energy IceCube events and γ-ray blazars classified as intermediate- and high-synchrotron-peaked BL Lacs. Subsequent studies have investigated the optical properties of these sources, compiled and analyzed public multiwavelength data, and constrained their individual neutrino emission based on public IceCube point-source data. Aims. We provide a theoretical interpretation of public multiwavelength and neutrino point source data for the 32 BL Lac objects in the sample previously associated with an IceCube alert event. We combined the individual source results to draw conclusions regarding the multimesssenger properties of the sample and the required power in relativistic protons. Methods. We performed particle interaction modeling using open-source numerical simulation software. We constrained the model parameters using a novel and unique approach that simultaneously describes the host galaxy contribution, the observed synchrotron peak properties, the average multiwavelength fluxes, and, where possible, the IceCube point source constraints. Results. We show that a single-zone leptohadronic model can describe the multiwavelength broadband fluxes from all 32 IceCube candidates. In some cases, the model suggests that hadronic emission may contribute a considerable fraction of the γ-ray flux. The required power in relativistic protons ranges from a few percent to a factor of ten of the Eddington luminosity, which is energetically less demanding compared to other leptohadronic blazar models in recent literature. The model can describe the 68% confidence level IceCube flux for a large fraction of the masquerading BL Lacs in the sample, including TXS 0506+056; whereas, for true BL Lacs, the model predicts a low neutrino flux in the IceCube sensitivity range. Physically, this distinction is due to the presence of photons from broad line emission in masquerading BL Lacs, which increase the efficiency of hadronic interactions. The predicted neutrino flux peaks between a few petaelectronvolt and 100 PeV and scales positively with the flux in the gigaelectronvolt, megaelectronvolt, X-ray, and optical bands. Based on these results, we provide a list of the brightest neutrino emitters, which can be used for future searches targeting the 10–100 PeV regime.
Marcotulli L, Ajello M, Böttcher M, Coppi P, Costamante L, Di Gesu L, Errando M, García JA, Gokus A, Liodakis I, et al. The high energy X-ray probe (HEX-P): the most powerful jets through the lens of a superb X-ray eye. [Internet]. 2024;11:1290057. WebsiteAbstract
A fraction of the active supermassive black holes at the centers of galaxies in our Universe are capable of launching extreme kiloparsec-long relativistic jets. These jets are known multiband (radio to γ-ray) and multimessenger (neutrino) emitters, and some of them have been monitored over decades at all accessible wavelengths. However, several open questions remain unanswered about the processes powering these highly energetic phenomena. These jets intrinsically produce soft-to-hard X-ray emission that extends from E>0.1keV up to E>100keV, and simultaneous broadband X-ray coverage, combined with excellent timing and imaging capabilities, is required to uncover the physics of jets. Indeed, truly simultaneous soft-to-hard X-ray coverage, in synergy with current and upcoming high-energy facilities (such as IXPE, COSI, CTAO, etc.) and neutrino detectors (e.g., IceCube), would enable us to disentangle the particle population responsible for the high-energy radiation from these jets. A sensitive hard X-ray survey (F20−80keV<10−15ergcm−2s−1) could unveil the bulk of their population in the early Universe. Acceleration and radiative processes responsible for the majority of their X-ray emission would be pinned down by microsecond timing capabilities at both soft and hard X-rays. Furthermore, imaging jet structures for the first time in the hard X-ray regime could unravel the origin of their high-energy emission. The proposed Probe-class mission concept High Energy X-ray Probe (HEX-P) combines all these required capabilities, making it the crucial next-generation X-ray telescope in the multi-messenger, time-domain era. HEX-P will be the ideal mission to unravel the science behind the most powerful accelerators in the Universe.
Karavola D, Petropoulou M. A closer look at the electromagnetic signatures of Bethe-Heitler pair production process in blazars. [Internet]. 2024;2024:006. WebsiteAbstract
The "twin birth" of a positron and an electron by a photon in the presence of a nucleus, known as Bethe-Heitler pair production, is a key process in astroparticle physics. The Bethe-Heitler process offers a way of channeling energy stored in a population of relativistic protons (or nuclei) to relativistic pairs with extended distributions. Contrary to accelerated leptons, whose maximum energy is limited by radiative losses, the maximal energy of pairs is determined by the kinematics of the process and can be as high as the parent proton energy. We take a closer look at the features of the injected pair distribution, and provide a novel empirical function that describes the spectrum of pairs produced by interactions of single-energy protons with single-energy photons. The function is the kernel of the Bethe-Heitler pair production spectrum that replaces a double numerical integration involving the complex differential cross section of the process, and can be easily implemented in numerical codes. We further examine under which conditions Bethe-Heitler pairs produced in blazar jets can emit γ-ray photons via synchrotron radiation, thus providing an alternative to the inverse Compton scattering process for high-energy emission in jetted active galactic nuclei. For this purpose, we create 36 numerical models using the code ATHEνA optimized so that the Bethe-Heitler synchrotron emission dominates their γ-ray emission. After taking into consideration the broadband spectral characteristics of the source, the jet energetics, and the properties of radiation fields present in the blazar environment, we conclude that γ-rays in high-synchrotron-peaked blazars are unlikely to be produced by Bethe-Heitler pairs, because the emitting region is found to be opaque in photon-photon pair production at photon energies ≳ 10 GeV. On the contrary, γ-ray spectra of low-synchrotron-peaked blazars may arise from Bethe-Heitler pairs in regions of the jet with typical transverse size ∼ 1015 – 1016 cm and co-moving magnetic field 50 – 500 G. For such cases, an external thermal target photon field with temperatures T ∼ 4 · 102– 6 · 103 K is needed. The latter values could point to the dusty torus of the AGN. Interestingly, a Bethe-Heitler-dominated high-energy component is mostly found in models of intermediate-synchrotron peaked blazars, for a wide range of magnetic fields and source radii.
Petropoulou M, Mastichiadis A, Vasilopoulos G, Paneque D, Becerra González J, Zanias F. TeV pion bumps in the gamma-ray spectra of flaring blazars. [Internet]. 2024;685:A110. WebsiteAbstract
Context. Very high-energy (VHE, E > 100 GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Aims: Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions (π0 bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. Methods: We explored the conditions needed for the emergence of narrow π0 bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies ∼100 keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources. Results: We find that synchrotron-dominated flares with peak energies ≳100 keV can be ideal periods for the search of π0 bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct π0 bump centered at 3 TeV. Spectral fitting of simulated SSC+π0 spectra for the Cherenkov Telescope Array (CTA) show that a π0 bump could be detected at a 5σ significance level with a 30-min exposure. Conclusions: A harder VHE γ-ray spectrum than the usual SSC prediction or, more occasionally, a distinct narrow bump at VHE energies during hard X-ray flares, can be suggestive of a relativistic hadronic component in blazar jets that otherwise would remain hidden. The production of narrow features or spectral hardenings due to π0 decay in the VHE spectra of blazars is testable with the advent of CTA.
Stathopoulos SI, Petropoulou M, Vasilopoulos G, Mastichiadis A. LeHaMoC: A versatile time-dependent lepto-hadronic modeling code for high-energy astrophysical sources. [Internet]. 2024;683:A225. WebsiteAbstract
Context. Recent associations of high-energy neutrinos with active galactic nuclei (AGN) have revived the interest in leptohadronic models of radiation from astrophysical sources. The rapid increase in the amount of acquired multi-messenger data will require fast numerical models that may be applied to large source samples. Aims: We develop a time-dependent leptohadronic code, LeHaMoC, that offers several notable benefits compared to other existing codes, such as versatility and speed. Methods: LeHaMoC solves the Fokker-Planck equations of photons and relativistic particles (i.e. electrons, positrons, protons, and neutrinos) produced in a homogeneous magnetized source that may also be expanding. The code utilizes a fully implicit difference scheme that allows fast computation of steady-state and dynamically evolving physical problems. Results: We first present test cases where we compare the numerical results obtained with LeHaMoC against exact analytical solutions and numerical results computed with ATHEvA, a well-tested code of similar philosophy but a different numerical implementation. We find a good agreement (within 10-30%) with the numerical results obtained with ATHEvA without evidence of systematic differences. We then demonstrate the capabilities of the code through illustrative examples. First, we fit the spectral energy distribution from a jetted AGN in the context of a synchrotron-self Compton model and a proton-synchrotron model using Bayesian inference. Second, we compute the high-energy neutrino signal and the electromagnetic cascade induced by hadronic interactions in the corona of NGC 1068. Conclusions: LeHaMoC is easily customized to model a variety of high-energy astrophysical sources and has the potential to become a widely utilized tool in multi-messenger astrophysics. Instructions for downloading the code, accessing online documentation, and reproducing applications presented in this paper can be found at https://github.com/mariapetro/LeHaMoC Github repository.
Tzavellas A, Vasilopoulos G, Petropoulou M, Mastichiadis A, Stathopoulos SI. Application of neural networks to synchro-Compton blazar emission models. [Internet]. 2024;683:A185. WebsiteAbstract
Context. Jets from supermassive black holes at the centers of active galaxies are the most powerful and persistent sources of electromagnetic radiation in the Universe. To infer the physical conditions in the otherwise out-of-reach regions of extragalactic jets, we usually rely on fitting their spectral energy distributions (SEDs). The calculation of radiative models for the jet's non-thermal emission usually relies on numerical solvers of coupled partial differential equations. Aims: In this work, we use machine learning to tackle the problem of high computational complexity to significantly reduce the SED model evaluation time, which is necessary for SED fittings carried out with Bayesian inference methods. Methods: We computed the SEDs based on the synchrotron self-Compton model for blazar emission using the radiation code ATHEvA. We used them to train neural networks (NNs) to explore whether they can replace the original code, which is computationally expensive. Results: We find that a NN with gated recurrent unit neurons (GRUN) can effectively replace the ATHEvA leptonic code for this application, while it can be efficiently coupled with Markov chain Monte Carlo (MCMC) and nested sampling algorithms for fitting purposes. We demonstrate this approach through an application to simulated data sets, as well as a subsequent application to observational data. Conclusions: We present a proof-of-concept application of NNs to blazar science as the first step in a list of future applications involving hadronic processes and even larger parameter spaces. We offer this tool to the community through a public repository. The results of our work are available in GitHub; https://github.com/tzavellas/blazar_ml. This includes: (a) the NN and accompanied code produced to train them, (b) code for visualization of results in python and jupyter notebooks with instructions, and (c) part of the ATHEvA datasets that can be used for evaluation and plotting examples.
Fiorillo DFG, Petropoulou M, Comisso L, Peretti E, Sironi L. TeV Neutrinos and Hard X-Rays from Relativistic Reconnection in the Corona of NGC 1068. [Internet]. 2024;961:L14. WebsiteAbstract
The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of nonthermal protons within a compact "coronal" region close to the central black hole. The acceleration mechanism of these nonthermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and nonthermal protons is established in the reconnection region. Motivated by recent 3D particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is ${{dn}}_{p}/{{dE}}_{p}\propto {E}_{p}^{-1}$ below the break and ${{dn}}_{p}/{{dE}}_{p}\propto {E}_{p}^{-s}$ above the break, with IceCube neutrino observations suggesting s ≃ 3. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron-positron pairs that makes the corona moderately Compton thick.