-
A magnetized strongly turbulent corona as the source of neutrinos from NGC 1068
Authors:
Damiano F. G. Fiorillo,
Luca Comisso,
Enrico Peretti,
Maria Petropoulou,
Lorenzo Sironi
Abstract:
The cores of active galactic nuclei (AGN) 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 which is the site and mechanism of cosmic ray acceleration. One candidate is the magnetized turbulence surrounding…
▽ More
The cores of active galactic nuclei (AGN) 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 which 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 non-resonant, 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.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
The Spectra of IceCube Neutrino (SIN) candidate sources -- V. Modeling and interpretation of multiwavelength and neutrino data
Authors:
Xavier Rodrigues,
Martina Karl,
Paolo Padovani,
Paolo Giommi,
Simona Paiano,
Renato Falomo,
Maria Petropoulou,
Foteini Oikonomou
Abstract:
A correlation has been reported between the arrival directions of high-energy IceCube events and gamma-ray blazars classified as intermediate- and high-synchrotron-peaked BL Lacs. Subsequent studies have investigated the optical properties of these sources, analyzed public multiwavelength data, and constrained their individual neutrino emission based on public IceCube data. In this work, we provid…
▽ More
A correlation has been reported between the arrival directions of high-energy IceCube events and gamma-ray blazars classified as intermediate- and high-synchrotron-peaked BL Lacs. Subsequent studies have investigated the optical properties of these sources, analyzed public multiwavelength data, and constrained their individual neutrino emission based on public IceCube data. In this work, 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 have 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, the multiwavelength fluxes, and the IceCube point source constraints. We show that a single-zone leptohadronic model can describe the multiwavelength fluxes from the 32 IceCube candidates. In some cases, the model suggests that hadronic emission may contribute a considerable fraction of the gamma-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. The predicted neutrino flux peaks between a few PeV and 100 PeV and scales positively with the flux in the GeV, MeV, 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.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
The role of magnetospheric current sheets in pair enrichment and ultra-high energy proton acceleration in M87*
Authors:
S. I. Stathopoulos,
M. Petropoulou,
L. Sironi,
D. Giannios
Abstract:
Recent advances in numerical simulations of magnetically arrested accretion onto supermassive black holes have shed light on the formation and dynamics of magnetospheric current sheets near the black hole horizon. By considering the pair magnetization $σ_{\rm e}$ in the upstream region and the mass accretion rate $\dot{m}$ (in units of the Eddington mass accretion rate) as free parameters we estim…
▽ More
Recent advances in numerical simulations of magnetically arrested accretion onto supermassive black holes have shed light on the formation and dynamics of magnetospheric current sheets near the black hole horizon. By considering the pair magnetization $σ_{\rm e}$ in the upstream region and the mass accretion rate $\dot{m}$ (in units of the Eddington mass accretion rate) as free parameters we estimate the strength of the magnetic field and develop analytical models, motivated by recent three-dimensional particle-in-cell simulations, to describe the populations of relativistic electron and positrons (pairs) in the reconnection region. Applying our model to M87*, we numerically compute the non-thermal photon spectra for various values of $σ_e$. We show that pairs that are accelerated up to the synchrotron radiation-limited energy while meandering across both sides of the current sheet, can produce MeV flares with luminosity of $\sim 10^{41}$~erg s$^{-1}$ -- independent of $σ_e$ -- for a black hole accreting at $\dot{m}=10^{-5}$. Pairs that are trapped in the transient current sheet can produce X-ray counterparts to the MeV flares, lasting about a day for current sheets with length of a few gravitational radii. We also show that the upstream plasma can be enriched due to photon-photon pair creation, and derive a new equilibrium magnetization of $σ_e \sim 10^3-10^4$ for $\dot{m}=10^{-6} - 10^{-5}$. Additionally, we explore the potential of magnetospheric current sheets to accelerate protons to ultra-high energies, finding that while acceleration to such energies is limited by various loss mechanisms, such as synchrotron and photopion losses from the non-thermal emission from pairs, maximal proton energies in the range of a few EeV are attainable in magnetospheric sheets forming around supermassive sub-Eddington accreting black holes.
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
A closer look at the electromagnetic signatures of Bethe-Heitler pair production process in blazars
Authors:
Despina Karavola,
Maria Petropoulou
Abstract:
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…
▽ More
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. 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 low- and intermediate-peaked blazars may arise from Bethe-Heitler pairs in regions of the jet with typical transverse size $\sim 10^{15}$ cm and co-moving magnetic field $5-500$ G.
△ Less
Submitted 30 May, 2024; v1 submitted 10 January, 2024;
originally announced January 2024.
-
Application of neural networks to synchro-Compton blazar emission models
Authors:
A. Tzavellas,
G. Vasilopoulos,
M. Petropoulou,
A. Mastichiadis,
S. I. Stathopoulos
Abstract:
Jets from supermassive black holes in the centers of active galaxies are the most powerful 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 of their spectral energy distribution (SED). The calculation of radiative models for the jet non-thermal emission usually re…
▽ More
Jets from supermassive black holes in the centers of active galaxies are the most powerful 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 of their spectral energy distribution (SED). The calculation of radiative models for the jet non-thermal emission usually relies on numerical solvers of coupled partial differential equations. In this work machine learning is used to tackle the problem of high computational complexity in order to significantly reduce the SED model evaluation time, which is needed for SED fitting with Bayesian inference methods. We compute SEDs based on the synchrotron self-Compton model for blazar emission using the radiation code ATHE$ν$A, and use them to train Neural Networks exploring whether these can replace the original computational expensive code. We find that a Neural Network with Gated Recurrent Unit neurons can effectively replace the ATHE$ν$A leptonic code for this application, while it can be efficiently coupled with MCMC and nested sampling algorithms for fitting purposes. We demonstrate this through an application to simulated data sets and with an application to observational data. We offer this tool in the community through a public repository. We present a proof-of-concept application of neural networks to blazar science. This is the first step in a list of future applications involving hadronic processes and even larger parameter spaces.
△ Less
Submitted 10 November, 2023;
originally announced November 2023.
-
The High Energy X-ray Probe (HEX-P): the most powerful jets through the lens of a superb X-ray eye
Authors:
Lea Marcotulli,
Marco Ajello,
Markus Böttcher,
Paolo Coppi,
Luigi Costamante,
Laura Di Gesu,
Manel Errando,
Javier A. García,
Andrea Gokus,
Ioannis Liodakis,
Greg Madejski,
Kristin Madsen,
Alberto Moretti,
Riccardo Middei,
Felicia McBride,
Maria Petropoulou,
Bindu Rani,
Tullia Sbarrato,
Daniel Stern,
Georgios Vasilopoulos,
Michael Zacharias,
Haocheng Zhang,
the HEX-P Collaboration
Abstract:
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 several decades at all accessible wavelengths. However, many open questions remain unanswered about the pro…
▽ More
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 several decades at all accessible wavelengths. However, many 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\sim0.1\,\rm keV$ up to $E>100\,\rm keV$. 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 ($F_{8-24\,\rm keV}<10^{-15}\,\rm erg~cm^{-2}~s^{-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.
△ Less
Submitted 8 November, 2023;
originally announced November 2023.
-
TeV neutrinos and hard X-rays from relativistic reconnection in the corona of NGC 1068
Authors:
Damiano F. G. Fiorillo,
Maria Petropoulou,
Luca Comisso,
Enrico Peretti,
Lorenzo Sironi
Abstract:
The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of non-thermal protons within a compact "coronal" region close to the central black hole. The acceleration mechanism of these non-thermal 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 su…
▽ More
The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of non-thermal protons within a compact "coronal" region close to the central black hole. The acceleration mechanism of these non-thermal 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 non-thermal protons is established in the reconnection region. Motivated by recent three-dimensional 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 \simeq 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.
△ Less
Submitted 19 January, 2024; v1 submitted 27 October, 2023;
originally announced October 2023.
-
An X-rays-to-radio investigation of the nuclear polarization from the radio-galaxy Centaurus A
Authors:
Frédéric Marin,
Thibault Barnouin,
Steven R. Ehlert,
Abel Lawrence Peirson,
Enrique Lopez-Rodriguez,
Maria Petropoulou,
Kinwah Wu,
Iván Martí-Vidal
Abstract:
Centaurus A is one of the closest radio-galaxies to Earth. Its proximity allowed us to extensively study its active galactic nucleus but the core emission mechanism remains elusive because of local strong dust and gas obscuration. The capability of polarimetry to shave-off contaminating emission has been exploited without success in the near-infrared by previous studies but the very recent measure…
▽ More
Centaurus A is one of the closest radio-galaxies to Earth. Its proximity allowed us to extensively study its active galactic nucleus but the core emission mechanism remains elusive because of local strong dust and gas obscuration. The capability of polarimetry to shave-off contaminating emission has been exploited without success in the near-infrared by previous studies but the very recent measurement of the 2 - 8 keV polarization by the Imaging X-ray Polarimetry Explorer (IXPE) brought the question back to the fore. To determine what is the prevalent photon generation mechanism to the multi-wavelength emission from the core of Centaurus A, we retrieved from the archives the panchromatic polarization measurements of the central compact component. We built the total and polarized flux spectral energy distributions of the core and demonstrated that synchrotron self-Compton models nicely fit the polarized flux from the radio to the X-ray band. The linear polarization of the synchrotron continuum is perpendicular to the jet radio axis from the optical to the radio band, and parallel to it at higher energies. The observed smooth rotation of the polarization angle in the ultraviolet band is attributed to synchrotron emission from regions that are getting closer to the particle acceleration site, where the orientation of the jet's magnetic fields become perpendicular to the jet axis. This phenomenon support the shock acceleration mechanism for particle acceleration in Centaurus A, in line with IXPE observations of several high-synchrotron peak blazars.
△ Less
Submitted 6 October, 2023;
originally announced October 2023.
-
Bright blazar flares with CTA
Authors:
M. Cerruti,
J. Finke,
G. Grolleron,
J. P. Lenain,
T. Hovatta,
M. Joshi,
E. Lindfors,
P. Morris,
M. Petropoulou,
P. Romano,
S. Vercellone,
M. Zacharias
Abstract:
The TeV extragalactic sky is dominated by blazars, radio-loud active galactic nuclei with a relativistic jet pointing towards the Earth. Blazars show variability that can be quite exceptional both in terms of flux (orders of magnitude of brightening) and time (down to the minute timescale). This bright flaring activity contains key information on the physics of particle acceleration and photon pro…
▽ More
The TeV extragalactic sky is dominated by blazars, radio-loud active galactic nuclei with a relativistic jet pointing towards the Earth. Blazars show variability that can be quite exceptional both in terms of flux (orders of magnitude of brightening) and time (down to the minute timescale). This bright flaring activity contains key information on the physics of particle acceleration and photon production in the emitting region, as well as the structure and physical properties of the jet itself. The TeV band is accessed from the ground by Cherenkov telescopes that image the pair cascade triggered by the interaction of the gamma ray with the Earth's atmosphere. The Cherenkov Telescope Array (CTA) represents the upcoming generation of imaging atmospheric Cherenkov telescopes, with a significantly higher sensitivity and larger energy coverage with respect to current instruments. It will thus provide us with unprecedented statistics on blazar light-curves and spectra. In this contribution we present the results from realistic simulations of CTA observations of bright blazar flares, taking as input state-of-the-art numerical simulations of blazar emission models and including all relevant observational constraints.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
TeV pion bumps in the gamma-ray spectra of flaring blazars
Authors:
M. Petropoulou,
A. Mastichiadis,
G. Vasilopoulos,
D. Paneque,
J. Becerra González,
F. Zanias
Abstract:
(Abridged) 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. Motivated by this rare observation, we prop…
▽ More
(Abridged) 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. 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. We explored the conditions needed for the emergence of narrow $π^0$ bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies $\sim 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. We find that synchrotron-dominated flares with peak energies $\gtrsim100$ 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.
△ Less
Submitted 5 March, 2024; v1 submitted 27 August, 2023;
originally announced August 2023.
-
LeHaMoC: a versatile time-dependent lepto-hadronic modeling code for high-energy astrophysical sources
Authors:
S. I. Stathopoulos,
M. Petropoulou,
G. Vasilopoulos,
A. Mastichiadis
Abstract:
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 soon fast numerical models that may be applied to large source samples. We develop a time-dependent leptohadronic code, LeHaMoC, that offers several no…
▽ More
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 soon fast numerical models that may be applied to large source samples. We develop a time-dependent leptohadronic code, LeHaMoC, that offers several notable benefits compared to other existing codes, such as versatility and speed. 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. We first present test cases where we compare the numerical results obtained with LeHaMoC against exact analytical solutions and numerical results computed with ATHE$ν$A, a well-tested code of similar philosophy but different numerical implementation. We find a good agreement (within 10-30%) with the numerical results obtained with ATHE$ν$A 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. 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.
△ Less
Submitted 5 July, 2024; v1 submitted 11 August, 2023;
originally announced August 2023.
-
Quantitative comparisons of VHE gamma-ray blazar flares with relativistic reconnection models
Authors:
J. Jormanainen,
T. Hovatta,
I. M. Christie,
E. Lindfors,
M. Petropoulou,
I. Liodakis
Abstract:
The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute time scales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection…
▽ More
The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute time scales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection layer, is a candidate solution to this problem. In this work, we employ state-of-the-art particle-in-cell simulations in a statistical comparison with observations of a flaring episode of a well-known blazar, Mrk 421, at very high energy (VHE, E > 100 GeV). We tested the predictions of our model by generating simulated VHE light curves that we compared quantitatively with methods that we have developed for a precise evaluation of theoretical and observed data. With our analysis, we can constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma, viewing angle and the reconnection layer orientation in the blazar jet. Our analysis favours parameter spaces with magnetic field strength 0.1 G, rather large viewing angles (6-8 degrees), and misaligned layer angles, offering a strong candidate explanation for the Doppler crisis often observed in the jets of high synchrotron peaking blazars.
△ Less
Submitted 11 August, 2023;
originally announced August 2023.
-
Possible contribution of X-ray binary jets to the Galactic cosmic ray and neutrino flux
Authors:
Dimitrios Kantzas,
Sera Markoff,
Alex J. Cooper,
Daniele Gaggero,
Maria Petropoulou,
Pedro De La Torre Luque
Abstract:
For over a century, the identification of high-energy cosmic ray (CR) sources remains an open question. For Galactic CRs with energy up to $10^{15}$ eV, supernova remnants (SNRs) have traditionally been thought the main candidate source. However, recent TeV gamma-ray observations have questioned the SNR paradigm. Propagating CRs are deflected by the Galactic magnetic field, hence, gamma-rays and n…
▽ More
For over a century, the identification of high-energy cosmic ray (CR) sources remains an open question. For Galactic CRs with energy up to $10^{15}$ eV, supernova remnants (SNRs) have traditionally been thought the main candidate source. However, recent TeV gamma-ray observations have questioned the SNR paradigm. Propagating CRs are deflected by the Galactic magnetic field, hence, gamma-rays and neutrinos produced via inelastic hadronic interactions are the only means for unveiling the CR sources. In this work, we study the gamma-ray and neutrino emission produced by CRs accelerated inside Galactic jets of stellar-mass black holes in X-ray binaries (BHXBs). We calculate the intrinsic neutrino emission of two prototypical BHXBs, Cygnus X-1 and GX 339-4, for which we have high-quality, quasi-simultaneous multiwavelength spectra. Based on these prototypical sources, we discuss the likelihood of the 35 known Galactic BHXBs to be efficient CR accelerators. Moreover, we estimate the potential contribution to the CR spectrum of a viable population of BHXBs that reside in the Galactic plane. When these BHXBs go into outburst, they may accelerate particles up to 100s of TeV that contribute to the diffuse gamma-ray and neutrino spectra while propagating in the Galactic medium. Using HERMES, an open-source code that calculates the hadronic processes along the line of sight, we discuss the contribution of BHXBs to the diffuse gamma-ray and neutrino fluxes, and compare these to their intrinsic gamma-ray and neutrino emissions. Finally, we discuss the contribution of BHXBs to the observed spectrum of Galactic CRs.
△ Less
Submitted 22 June, 2023;
originally announced June 2023.
-
The origin of power-law spectra in relativistic magnetic reconnection
Authors:
Hao Zhang,
Lorenzo Sironi,
Dimitrios Giannios,
Maria Petropoulou
Abstract:
Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model - supported and benchmarked with large-scale three-dimensional particle-in-cell simulations - that elucidates the physics governing the generation of power-law energy spect…
▽ More
Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model - supported and benchmarked with large-scale three-dimensional particle-in-cell simulations - that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factor $γ\gtrsim 3σ$ (here, $σ$ is the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is $t_{\rm acc}\sim γ\,η_{\rm rec}^{-1}ω_{\rm c}^{-1}\simeq 20\,γ\,ω_{\rm c}^{-1}$, where $η_{\rm rec}\simeq0.06$ is the inflow speed in units of the speed of light and $ω_{\rm c}=eB_0/mc$ is the gyrofrequency in the upstream magnetic field. They leave the region of active energization after $t_{\rm esc}$, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measure $t_{\rm esc}$ in our simulations and find that $t_{\rm esc}\sim t_{\rm acc}$ for $σ\gtrsim {\rm few}$. This leads to a universal (i.e., $σ$-independent) power-law spectrum $dN_{\rm free}/dγ\propto γ^{-1}$ for the particles undergoing active acceleration, and $dN/dγ\propto γ^{-2}$ for the overall particle population. Our results help shedding light on the ubiquitous presence of power-law particle and photon spectra in astrophysical non-thermal sources.
△ Less
Submitted 23 February, 2023;
originally announced February 2023.
-
Dust-scattering rings of GRB 221009A as seen by the Neil Gehrels Swift satellite: can we count them all?
Authors:
Georgios Vasilopoulos,
Despina Karavola,
Stamatios I. Stathopoulos,
Maria Petropoulou
Abstract:
We present the first results for the dust-scattering rings of GRB 221009A, coined as the GRB of the century, as observed by the Neil Gehrels Swift satellite. We perform analysis of both time resolved observations and stacked data. The former approach enable us to study the expansion of the most prominent rings, associate their origin with the prompt X-ray emission of the GRB and determine the loca…
▽ More
We present the first results for the dust-scattering rings of GRB 221009A, coined as the GRB of the century, as observed by the Neil Gehrels Swift satellite. We perform analysis of both time resolved observations and stacked data. The former approach enable us to study the expansion of the most prominent rings, associate their origin with the prompt X-ray emission of the GRB and determine the location of the dust layers. The stacked radial profiles increase the signal-to-noise ratio of the data and allows detection of fainter and overlapping peaks in the angular profile. We find a total of 16 dust concentrations (with hints of even more) that span about 15 kpc in depth and could be responsible for the highly structured X-ray angular profiles. By comparing the relative scattered fluxes of the five most prominent rings we show that the layer with the largest amount of dust is located at about 0.44 kpc away from us. We finally compare the location of the dust layers with results from experiments that study the 3D structure of our Galaxy via extinction or CO radio observations, and highlight the complementarity of dust X-ray tomography to these approaches.
△ Less
Submitted 8 February, 2023; v1 submitted 5 February, 2023;
originally announced February 2023.
-
Multi-messenger model for the prompt emission from GRB 221009A
Authors:
Annika Rudolph,
Maria Petropoulou,
Walter Winter,
Željka Bošnjak
Abstract:
We present a multi-messenger model for the prompt emission from GRB 221009A within the internal shock scenario. We consider the time-dependent evolution of the outflow with its impact on the observed light curve from multiple collisions, and the self-consistent generation of the electromagnetic spectrum in synchrotron and inverse Compton-dominated scenarios. Our leptohadronic model includes UHE pr…
▽ More
We present a multi-messenger model for the prompt emission from GRB 221009A within the internal shock scenario. We consider the time-dependent evolution of the outflow with its impact on the observed light curve from multiple collisions, and the self-consistent generation of the electromagnetic spectrum in synchrotron and inverse Compton-dominated scenarios. Our leptohadronic model includes UHE protons potentially accelerated in the outflow, and their feedback on spectral energy distribution and on the neutrino emission. We find that we can roughly reproduce the observed light curves with an engine with varying ejection velocity of ultra-relativistic material, which has an intermediate quiescent period of about 200 seconds and a variability timescale of $\sim1$~s. We consider baryonic loadings of 3 and 30 that are compatible with the hypothesis that the highest-energetic LHAASO photons might come from UHECR interactions with the extragalactic background light, and the paradigm that energetic GRBs may power the UHECR flux. For these values and the high dissipation radii considered we find consistency with the non-observation of neutrinos and no significant signatures on the electromagnetic spectrum. Inverse Compton-dominated scenarios from the prompt emission are demonstrated to lead to about an order of magnitude higher fluxes in the HE-range; this enhancement is testable by its spectral impact in the Fermi-GBM and LAT ranges.
△ Less
Submitted 1 April, 2023; v1 submitted 1 December, 2022;
originally announced December 2022.
-
Multi-collision internal shock lepto-hadronic models for energetic GRBs
Authors:
Annika Rudolph,
Maria Petropoulou,
Željka Bošnjak,
Walter Winter
Abstract:
For a sub-population of energetic Gamma-Ray Bursts (GRBs), a moderate baryonic loading may suffice to power Ultra-High-Energy Cosmic Rays (UHECRs). Motivated by this, we study the radiative signatures of cosmic-ray protons in the prompt phase of energetic GRBs. Our framework is the internal shock model with multi-collision descriptions of the relativistic ejecta (with different emission regions al…
▽ More
For a sub-population of energetic Gamma-Ray Bursts (GRBs), a moderate baryonic loading may suffice to power Ultra-High-Energy Cosmic Rays (UHECRs). Motivated by this, we study the radiative signatures of cosmic-ray protons in the prompt phase of energetic GRBs. Our framework is the internal shock model with multi-collision descriptions of the relativistic ejecta (with different emission regions along the jet), plus time-dependent calculations of photon and neutrino spectra. Our GRB prototypes are motivated by {\em Fermi}-LAT detected GRBs (including GRB~221009A) for which further, owing to the large energy flux, neutrino non-observation of single events may pose a strong limit on the baryonic loading. We study the feedback of protons on electromagnetic spectra in synchrotron- and inverse Compton-dominated scenarios to identify the multi-wavelength signatures, to constrain the maximally allowed baryonic loading, and to point out the differences between hadronic and inverse Compton signatures. We find that hadronic signatures appear as correlated flux increases in the optical-UV to soft X-ray and GeV to TeV gamma-ray ranges in the synchrotron scenarios, whereas they are difficult to identify in inverse Compton-dominated scenarios. We demonstrate that baryonic loadings around 10, which satisfy the UHECR energetic requirements, do not distort the predicted photon spectra in the {\em Fermi}-GBM range and are consistent with constraints from neutrino data if the collision radii are large enough (i.e., the time variability is not too short). It therefore seems plausible that under the condition of large dissipation radii a population of energetic GRBs can be the origin of the UHECRs.
△ Less
Submitted 1 April, 2023; v1 submitted 1 December, 2022;
originally announced December 2022.
-
A Bayesian approach for torque modelling of BeXRB pulsars with application to super-Eddington accretors
Authors:
A. S. Karaferias,
G. Vasilopoulos,
M. Petropoulou,
P. A. Jenke,
C. A. Wilson-Hodge,
C. Malacaria
Abstract:
In this study we present a method to estimate posterior distributions for standard accretion torque model parameters and binary orbital parameters for X-ray binaries using a nested sampling algorithm for Bayesian Parameter Estimation. We study the spin evolution of two Be X-ray binary systems in the Magellanic Clouds, RX J0520.5-6932 and RX J0209-7427, during major outbursts, in which they surpass…
▽ More
In this study we present a method to estimate posterior distributions for standard accretion torque model parameters and binary orbital parameters for X-ray binaries using a nested sampling algorithm for Bayesian Parameter Estimation. We study the spin evolution of two Be X-ray binary systems in the Magellanic Clouds, RX J0520.5-6932 and RX J0209-7427, during major outbursts, in which they surpassed the Eddington-limit. Moreover, we apply our method to the recently discovered Swift J0243.6+6124; the only known Galactic pulsating ultra-luminous X-ray source. This is an excellent candidate for studying the disc evolution at super-Eddington accretion rates, for its luminosity span several orders of magnitude during its outburst, with a maximum $L_{\rm X}$ that exceeded the Eddington limit by a factor of $\sim 10$. Our method, when applied to RX J0520.5-6932 and RX J0209-7427, is able to identify the more favourable torque model for each system, while yielding meaningful ranges for the NS and orbital parameters. Our analysis for Swift J0243.6+6124 illustrates that, contrary to the standard torque model predictions, the magnetospheric radius and the Alfvén radius are not proportional to each other when surpassing the Eddington limit. Reported distance estimates of this source range between 5 and 7 kpc. Smaller distances require non-typical neutron star properties (i.e. mass and radius) and possibly lower radiative efficiency of the accretion column.
△ Less
Submitted 29 November, 2022;
originally announced November 2022.
-
Hadronic signatures from magnetically dominated baryon-loaded AGN jets
Authors:
Maria Petropoulou,
Filippos Psarras,
Dimitrios Giannios
Abstract:
Blazars are a rare class of active galactic nuclei (AGN) with relativistic jets pointing towards the observer. Jets are thought to be launched as Poynting-flux dominated outflows that accelerate to relativistic speeds at the expense of the available magnetic energy. In this work, we consider electron-proton jets and assume that particles are energized via magnetic reconnection in parts of the jet…
▽ More
Blazars are a rare class of active galactic nuclei (AGN) with relativistic jets pointing towards the observer. Jets are thought to be launched as Poynting-flux dominated outflows that accelerate to relativistic speeds at the expense of the available magnetic energy. In this work, we consider electron-proton jets and assume that particles are energized via magnetic reconnection in parts of the jet where the magnetization is still high ($σ\ge 1$). The magnetization and bulk Lorentz factor $Γ$ are related to the available jet energy per baryon as $μ=Γ(1+σ)$. We adopt an observationally motivated relation between $Γ$ and the mass accretion rate into the black hole $\dot{m}$, which also controls the luminosity of external radiation fields. We numerically compute the photon and neutrino jet emission as a function of $μ$ and $σ$. We find that the blazar SED is produced by synchrotron and inverse Compton radiation of accelerated electrons, while the emission of hadronic-related processes is subdominant except for the highest magnetization considered. We show that low-luminosity blazars ($L_γ \lesssim 10^{45}$ erg s$^{-1}$) are associated with less powerful, slower jets with higher magnetizations in the jet dissipation region. Their broadband photon spectra resemble those of BL Lac objects, and the expected neutrino luminosity is $L_{ν+\barν}\sim (0.3-1)\, L_γ$. High-luminosity blazars ($L_γ \gg 10^{45}$ erg s$^{-1}$) are associated with more powerful, faster jets with lower magnetizations. Their broadband photon spectra resemble those of flat spectrum radio quasars, and they are expected to be dim neutrino sources with $L_{ν+\barν}\ll L_γ$.
△ Less
Submitted 3 November, 2022;
originally announced November 2022.
-
Flare Duty Cycle of Gamma-Ray Blazars and Implications for High-Energy Neutrino Emission
Authors:
Kenji Yoshida,
Maria Petropoulou,
Kohta Murase,
Foteini Oikonomou
Abstract:
Gamma-ray flares of blazars may be accompanied by high-energy neutrinos due to interactions of high-energy cosmic rays in the jet with photons, as suggested by the detection of the high-energy neutrino IceCube-170922A during a major gamma-ray flare from blazar TXS 0506+056 at the $\sim3σ$ significance level. In this work, we present a statistical study of gamma-ray emission from blazars to constra…
▽ More
Gamma-ray flares of blazars may be accompanied by high-energy neutrinos due to interactions of high-energy cosmic rays in the jet with photons, as suggested by the detection of the high-energy neutrino IceCube-170922A during a major gamma-ray flare from blazar TXS 0506+056 at the $\sim3σ$ significance level. In this work, we present a statistical study of gamma-ray emission from blazars to constrain the contribution of gamma-ray flares to their neutrino output. We construct weekly binned light curves for 145 gamma-ray bright blazars in the {\it Fermi} Large Area Telescope (LAT) Monitored Source List adding TXS 0506+056. We derive the fraction of time spent in the flaring state (flare duty cycle) and the fraction of energy released during each flare from the light curves with a Bayesian blocks algorithm. We find that blazars with lower flare duty cycles and energy fractions are more numerous among our sample. We identify a significant difference in flare duty cycles between blazar sub-classes at a significance level of 5~\%. Then using a general scaling relation for the neutrino and gamma-ray luminosities, $L_ν \propto (L_γ)^γ$ with a weighting exponent of $γ = 1.0 - 2.0$, normalized to the quiescent gamma-ray or X-ray flux of each blazar, we evaluate the neutrino energy flux of each gamma-ray flare. The gamma-ray flare distribution indicates that blazar neutrino emission may be dominated by flares for $γ\gtrsim1.5$. The neutrino energy fluxes for one-week and 10-year bins are compared with the declination-dependent IceCube sensitivity to constrain the standard neutrino emission models for gamma-ray flares. Finally, we present the upper-limit contribution of blazar gamma-ray flares to the isotropic diffuse neutrino flux.
△ Less
Submitted 8 July, 2023; v1 submitted 18 October, 2022;
originally announced October 2022.
-
The hunt for extraterrestrial high-energy neutrino counterparts
Authors:
I. Liodakis,
T. Hovatta,
V. Pavlidou,
A. C. S. Readhead,
R. D. Blandford,
S. Kiehlmann,
E. Lindfors,
W. Max-Moerbeck,
T. J. Pearson,
M. Petropoulou
Abstract:
The origin of Petaelectronvolt (PeV) astrophysical neutrinos is fundamental to our understanding of the high-energy Universe. Apart from the technical challenges of operating detectors deep below ice, oceans, and lakes, the phenomenological challenges are even greater than those of gravitational waves; the sources are unknown, hard to predict, and we lack clear signatures. Neutrino astronomy there…
▽ More
The origin of Petaelectronvolt (PeV) astrophysical neutrinos is fundamental to our understanding of the high-energy Universe. Apart from the technical challenges of operating detectors deep below ice, oceans, and lakes, the phenomenological challenges are even greater than those of gravitational waves; the sources are unknown, hard to predict, and we lack clear signatures. Neutrino astronomy therefore represents the greatest challenge faced by the astronomy and physics communities thus far. The possible neutrino sources range from accretion disks and tidal disruption events, to relativistic jets and galaxy clusters with blazar TXS~0506+056 the most compelling association thus far. Since that association, immense effort has been put into proving or disproving that jets are indeed neutrino emitters, but to no avail. By generating simulated neutrino counterpart samples, we explore the potential of detecting a significant correlation of neutrinos with jets from active galactic nuclei. We find that, given the existing challenges, even our best experiments could not have produced a $>3σ$ result. Larger programs over the next few years will be able to detect a significant correlation only if the brightest radio sources, rather than all jetted active galactic nuclei, are neutrino emitters. We discuss the necessary strategies required to steer future efforts into successful experiments.
△ Less
Submitted 15 August, 2022;
originally announced August 2022.
-
A study of natural frequencies in a dynamic corona-disk system
Authors:
A. Mastichiadis,
M. Petropoulou,
N. D. Kylafis
Abstract:
Black-hole X-ray binaries (BHXRBs) in the hard and hard-intermediate spectral (and temporal) states exhibit in their power spectra characteristic frequencies called type-C quasi-periodic oscillations (QPOs). Various models that can explain them with various degrees of success have been proposed, but a definitive answer is still missing. The hot Comptonizing corona interacting with the cold accreti…
▽ More
Black-hole X-ray binaries (BHXRBs) in the hard and hard-intermediate spectral (and temporal) states exhibit in their power spectra characteristic frequencies called type-C quasi-periodic oscillations (QPOs). Various models that can explain them with various degrees of success have been proposed, but a definitive answer is still missing. The hot Comptonizing corona interacting with the cold accretion disk, both of which are central in understanding BHXRBs, is essentially a dynamical system. Our aim is to investigate if the radiative coupling between the two components can produce QPOs. We write and solve the time-dependent equations that describe energy conservation in the system corona - accretion disk. We examine both constant and variable mass accretion rates. By necessity, in this first investigation we use a simple model, but it contains all the essential ingredients. For a constant mass accretion rate and certain justifiable conditions, the dynamic corona - disk system exhibits oscillations, which die out after a few cycles. The characteristic frequencies of these oscillations are similar to the ones observed in the power spectra of BHXRBs. For most parameters, the natural frequencies persist even in the case of variable accretion rates. We argue that type-C QPOs in BHXRBs could, in principle, arise from the interaction of the hot Comptonizing corona with the much colder accretion disk. If this picture is correct, it has immediate implications for other systems that contain the above constituents, such as active galactic nuclei.
△ Less
Submitted 7 May, 2022; v1 submitted 29 April, 2022;
originally announced April 2022.
-
Hadronic supercriticality in spherically expanding sources: application to GRB prompt emission
Authors:
Ioulia Florou,
Apostolos Mastichiadis,
Maria Petropoulou
Abstract:
Relativistic hadronic plasmas can become under certain conditions supercritical, abruptly and efficiently releasing the energy stored in protons through photon outbursts. Past studies have tried to relate the features of such hadronic supercriticalities (HSC) to the phenomenology of Gamma-Ray Burst (GRB) prompt emission. In this work we investigate, for the first time, HSC in adiabatically expandi…
▽ More
Relativistic hadronic plasmas can become under certain conditions supercritical, abruptly and efficiently releasing the energy stored in protons through photon outbursts. Past studies have tried to relate the features of such hadronic supercriticalities (HSC) to the phenomenology of Gamma-Ray Burst (GRB) prompt emission. In this work we investigate, for the first time, HSC in adiabatically expanding sources. We examine the conditions required to trigger HSC, study the role of expansion velocity, and discuss our results in relation to GRB prompt emission. We find multi-pulse light curves from slowly expanding regions ($u_{\rm exp}\lesssim 0.01 c)$ that are a manifestation of the natural HSC quasi-periodicity, while single-pulse light curves with a fast rise and slow decay are found for higher velocities. The formation of the photon spectrum is governed by an in-source electromagnetic cascade. The peak photon energy is $\sim 1$ MeV ($\sim 1$ GeV) for maximum proton energies $\sim 1-10$ PeV ($1-10$ EeV) assuming a jet Lorentz factor 100. Peak $γ$-ray luminosities are in the range $10^{49}-10^{52}$ erg s$^{-1}$, with the MeV-peaked spectra being $\sim 100-300$ times more luminous than their GeV-peaked analogues. HSC bursts peaking in the MeV are also copious $\sim 10$ TeV neutrino emitters, with an all-flavour fluence $\sim 10\%$ of the $γ$-ray one. The hypothesis that typical long-duration GRBs are powered by HSC could be tested in the near future with more sensitive neutrino telescopes like IceCube-Gen2.
△ Less
Submitted 17 April, 2022;
originally announced April 2022.
-
A multi-messenger study of the blazar PKS 0735+178: a new major neutrino source candidate
Authors:
N. Sahakyan,
P. Giommi,
P. Padovani,
M. Petropoulou,
D. Bégué,
B. Boccardi,
S. Gasparyan
Abstract:
The blazar PKS 0735+178 is possibly associated with multiple neutrino events observed by the IceCube, Baikal, Baksan, and KM3NeT neutrino telescopes while it was flaring in the $γ$-ray, X-ray, ultraviolet and optical bands. We present a detailed study of this peculiar blazar to investigate the temporal and spectral changes in the multi-wavelength emission when the neutrino events were observed. Th…
▽ More
The blazar PKS 0735+178 is possibly associated with multiple neutrino events observed by the IceCube, Baikal, Baksan, and KM3NeT neutrino telescopes while it was flaring in the $γ$-ray, X-ray, ultraviolet and optical bands. We present a detailed study of this peculiar blazar to investigate the temporal and spectral changes in the multi-wavelength emission when the neutrino events were observed. The analysis of Swift-XRT snapshots reveal a flux variability of more than a factor 2 in about $5\times10^3$ seconds during the observation on December 17, 2021. In the $γ$-ray band, the source was in its historical highest flux level at the time of the arrival of the neutrinos. The observational comparison between PKS 0735+178 and other neutrino source candidates, such as TXS 0506+056, PKS 1424+240, and GB6 J1542+6129, shows that all these sources share similar spectral energy distributions, very high radio and $γ$-ray powers, and parsec scale jet properties. Moreover, we present strong supporting evidence for PKS 0735+178 to be, like all the others, a masquerading BL Lac. We perform comprehensive modelling of the multiwavelength emission from PKS 0735+178 within one-zone lepto-hadronic models considering both internal and external photon fields and estimate the expected accompanying neutrino flux. The most optimistic scenario invokes a jet with luminosity close to the Eddington value and the interactions of $\sim$ PeV protons with an external UV photon field. This scenario predicts $\sim 0.067$ muon and antimuon neutrinos over the observed 3-week flare. Our results are consistent with the detection of one very-high-energy neutrino like IceCube-211208A.
△ Less
Submitted 6 December, 2022; v1 submitted 11 April, 2022;
originally announced April 2022.
-
The Spectra of IceCube Neutrino (SIN) candidate sources -- II. Source Characterisation
Authors:
P. Padovani,
P. Giommi,
R. Falomo,
F. Oikonomou,
M. Petropoulou,
T. Glauch,
E. Resconi,
A. Treves,
S. Paiano
Abstract:
Eight years after the first detection of high-energy astrophysical neutrinos by IceCube we are still almost clueless as regards to their origin, although the case for blazars being neutrino sources is getting stronger. After the first significant association at the $3 - 3.5\,σ$ level in time and space with IceCube neutrinos, i.e. the blazar TXS 0506+056 at $z=0.3365$, some of us have in fact selec…
▽ More
Eight years after the first detection of high-energy astrophysical neutrinos by IceCube we are still almost clueless as regards to their origin, although the case for blazars being neutrino sources is getting stronger. After the first significant association at the $3 - 3.5\,σ$ level in time and space with IceCube neutrinos, i.e. the blazar TXS 0506+056 at $z=0.3365$, some of us have in fact selected a unique sample of 47 blazars, out of which $\sim 16$ could be associated with individual neutrino track events detected by IceCube. Building upon our recent spectroscopy work on these objects, here we characterise them to determine their real nature and check if they are different from the rest of the blazar population. For the first time we also present a systematic study of the frequency of masquerading BL Lacs, i.e. flat-spectrum radio quasars with their broad lines swamped by non-thermal jet emission, in a $γ$-ray- and IceCube-selected sample, finding a fraction $>$ 24 per cent and possibly as high as 80 per cent. In terms of their broad-band properties, our sources appear to be indistinguishable from the rest of the blazar population. We also discuss two theoretical scenarios for neutrino emission, one in which neutrinos are produced in interactions of protons with jet photons and one in which the target photons are from the broad line region. Both scenarios can equally account for the neutrino-blazar correlation observed by some of us. Future observations with neutrino telescopes and X-ray satellites will test them out.
△ Less
Submitted 10 December, 2021;
originally announced December 2021.
-
A unified model for orphan and multi-wavelength blazar flares
Authors:
Ze-Rui Wang,
Ruo-Yu Liu,
Maria Petropoulou,
Foteini Oikonomou,
Rui Xue,
Xiang-Yu Wang
Abstract:
Blazars are a class of active galactic nuclei which host relativistic jets oriented close to the observer's line of sight. Blazars have very complex variability properties. Flares, namely flux variations around the mean value with a well-defined shape and duration, are one of the identifying properties of the blazar phenomenon. Blazars are known to exhibit multi-wavelength flares, but also "orphan…
▽ More
Blazars are a class of active galactic nuclei which host relativistic jets oriented close to the observer's line of sight. Blazars have very complex variability properties. Flares, namely flux variations around the mean value with a well-defined shape and duration, are one of the identifying properties of the blazar phenomenon. Blazars are known to exhibit multi-wavelength flares, but also "orphan" flares, namely flux changes that appear only in a specific energy range. Various models, sometimes at odds with each other, have been proposed to explain specific flares even for a single source, and cannot be synthesized into a coherent picture. In this paper, we propose a unified model for explaining orphan and multi-wavelength flares from blazars in a common framework. We assume that the blazar emission during a flare consists of two components: (i) a quasi-stable component that arises from the superposition of numerous but comparatively weak dissipation zones along the jet, forming the background (low-state) emission of the blazar, and (ii) a transient component, which is responsible for the sudden enhancement of the blazar flux, forming at a random distance along the jet by a strong energy dissipation event. Whether a multi-wavelength or orphan flare is emitted depends on the distance from the base of the jet where the dissipation occurs. Generally speaking, if the dissipation occurs at a small/large distance from the supermassive black hole, the inverse Compton/synchrotron radiation dominates and an orphan gamma-ray/optical flare tends to appear. On the other hand, we may expect a multi-wavelength flare if the dissipation occurs at an intermediate distance. We show that the model can successfully describe the spectral energy distribution of different flares from the flat spectrum radio quasar 3C 279 and the BL Lac object PKS 2155-304.
△ Less
Submitted 3 December, 2021;
originally announced December 2021.
-
High-energy neutrinos from X-rays flares of blazars frequently observed by the Swift X-Ray Telescope
Authors:
S. I. Stathopoulos,
M. Petropoulou,
P. Giommi,
G. Vasilopoulos,
P. Padovani,
A. Mastichiadis
Abstract:
Blazar flares have been suggested as ideal candidates for enhanced neutrino production. While the neutrino signal of $γ$-ray flares has been widely discussed, the neutrino yield of X-ray flares has received less attention. Here, we compute the predicted neutrino signal from X-ray flares detected in 66 blazars observed more than 50 times with the X-ray Telescope (XRT) on board the Neil Gehrels Swif…
▽ More
Blazar flares have been suggested as ideal candidates for enhanced neutrino production. While the neutrino signal of $γ$-ray flares has been widely discussed, the neutrino yield of X-ray flares has received less attention. Here, we compute the predicted neutrino signal from X-ray flares detected in 66 blazars observed more than 50 times with the X-ray Telescope (XRT) on board the Neil Gehrels Swift Observatory. We consider a scenario where X-ray flares are powered by synchrotron radiation of relativistic protons, and neutrinos are produced through photomeson interactions between protons with their own synchrotron X-ray photons. Using the 1 keV X-ray light curves for flare identification, the 0.5-10 keV fluence of each flare as a proxy for the all-flavour neutrino fluence, and the IceCube point-source effective area for different detector configurations, we calculate the number of muon and antimuon neutrinos above 100 TeV expected for IceCube from each flaring source. The bulk of the neutrino events from the sample originates from flares with durations $\sim 1-10$ d. Accounting for the X-ray flare duty cycle of the sources in the sample, which ranges between $\sim2$ and 24 per cent, we compute an average yearly neutrino rate for each source. The median of the distribution (in logarithm) is $\sim0.03$ yr$^{-1}$, with Mkn 421 having the highest predicted rate $1.2\pm 0.3$ yr$^{-1}$, followed by 3C 273 $(0.33\pm0.03$ yr$^{-1})$ and PG 1553+113 ($0.25\pm0.02$ yr$^{-1}$). Next-generation neutrino detectors together with regular X-ray monitoring of blazars could constrain the duty cycle of hadronic X-ray flares.
△ Less
Submitted 17 November, 2021;
originally announced November 2021.
-
Circumnuclear Dust in AP Librae and the source of its VHE emission
Authors:
Agniva Roychowdhury,
Eileen T. Meyer,
Markos Georganopoulos,
Peter Breiding,
Maria Petropoulou
Abstract:
The broad high-energy spectral component in blazars is usually attributed to various inverse Compton scattering processes in the relativistic jet, but has not been clearly identified in most cases due to degeneracies in physical models. AP Librae, a low-synchrotron-peaking BL Lac object (LBL) detected in 2015 by H.E.S.S. at very high energies (VHE; $>$ 0.5 TeV), has an extremely broad high-energy…
▽ More
The broad high-energy spectral component in blazars is usually attributed to various inverse Compton scattering processes in the relativistic jet, but has not been clearly identified in most cases due to degeneracies in physical models. AP Librae, a low-synchrotron-peaking BL Lac object (LBL) detected in 2015 by H.E.S.S. at very high energies (VHE; $>$ 0.5 TeV), has an extremely broad high-energy spectrum, covering $\sim$ 9 decades in energy. Standard synchrotron self-Compton models generally fail to reproduce the VHE emission, which has led to the suggestion that it might arise not from the blazar core, but on kiloparsec scales from inverse Compton scattering of cosmic microwave background (CMB) photons by a still-relativistic jet (IC/CMB). IC/CMB models for the TeV emission of AP Librae in prior works have implied a high level of infrared emission from the kpc-scale jet. With newly obtained Hubble Space Telescope imaging, we obtain a deep upper limit on the kpc-scale jet emission at 1.6 $μ$m, well below the expected level. High-resolution ALMA imaging in bands 3-9 reveals a residual dust disk signature after core subtraction, with a clearly thermal spectrum, and an extent ($\sim$500 pc) which matches with a non-jet residual emission seen after PSF subtraction in our 1.6 $μ$m HST imaging. We find that the unusually broad GeV and VHE emission in AP Librae can be reproduced through the combined IC scattering of photons from the CMB and the dust disk, respectively, by electrons in both the blazar core and sub-kpc jet.
△ Less
Submitted 28 October, 2021; v1 submitted 22 October, 2021;
originally announced October 2021.
-
Radio emission from colliding outflows in high-mass X-ray binaries with strongly magnetized neutron stars
Authors:
Margaritis Chatzis,
Maria Petropoulou,
Georgios Vasilopoulos
Abstract:
We present a toy model for radio emission in HMXBs with strongly magnetized neutron stars (NS) where a wind-collision region is formed by the NS outflow and the stellar wind of the massive companion. Radio emission is expected from the synchrotron radiation of shock-accelerated electrons and the free-free emission of the stellar wind. We found that the predicted relation between the GHz luminosity…
▽ More
We present a toy model for radio emission in HMXBs with strongly magnetized neutron stars (NS) where a wind-collision region is formed by the NS outflow and the stellar wind of the massive companion. Radio emission is expected from the synchrotron radiation of shock-accelerated electrons and the free-free emission of the stellar wind. We found that the predicted relation between the GHz luminosity ($L_R$) and the accretion X-ray luminosity ($L_X$) can be written as $L_R \propto L_X^b$ for most parameters. No correlation with X-rays is expected ($b=0$) when the thermal emission of the stellar wind dominates in radio. We typically find a steep correlation ($b=12/7$) for sub-Eddington X-ray luminosities and a more shallow one ($b=2(p-1)/7$) for super-Eddington X-ray luminosities, where $p$ is the power-law index of accelerated electrons. The maximum predicted radio luminosity is independent of the NS properties, while it depends on the stellar wind momentum, binary separation distance, and the minimum electron Lorentz factor. Using a Bayesian approach we modelled the radio observations of \sj that cover a wide range of mass accretion rates. Our results support a shock origin for the radio detections at sub-Eddington X-ray luminosities. However, no physically meaningful parameters could be found for the super-Eddington phase of the outburst, suggesting a different origin. Future observations with more sensitive instruments might reveal a large number of HMXBs with strongly magnetized NSs in radio, allowing determination of the slope in the $L_R-L_X$ relation, and putting the wind-collision scenario into test.
△ Less
Submitted 22 October, 2021;
originally announced October 2021.
-
Confronting observations of VHE gamma-ray blazar flares with reconnection models
Authors:
J. Jormanainen,
T. Hovatta,
E. Lindfors,
I. Christie,
M. Petropoulou,
I. Liodakis
Abstract:
Several models have been suggested to explain the fast gamma-ray variability observed in blazars, but its origin is still debated. One scenario is magnetic reconnection, a process that can efficiently convert magnetic energy to energy of relativistic particles accelerated in the reconnection layer. In our study, we compare results from state-of-the-art particle-in-cell simulations with observation…
▽ More
Several models have been suggested to explain the fast gamma-ray variability observed in blazars, but its origin is still debated. One scenario is magnetic reconnection, a process that can efficiently convert magnetic energy to energy of relativistic particles accelerated in the reconnection layer. In our study, we compare results from state-of-the-art particle-in-cell simulations with observations of blazars at Very High Energy (VHE, E > 100 GeV) gamma-rays. Our goal is to test our model predictions on fast gamma-ray variability with data and to constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma and the reconnection layer orientation in the blazar jet. For this first comparison, we used the remarkably well-sampled VHE gamma-ray light curve of Mrk 421 observed with the MAGIC and VERITAS telescopes in 2013.The simulated VHE light curves were generated using the observable parameters of Mrk 421, such as the jet power, bulk Lorentz factor, and the jet viewing angle, and sampled as real data. Our results pave the way for future model-to-data comparison with next-generation Cherenkov telescopes, which will help further constrain the different variability models.
△ Less
Submitted 17 September, 2021;
originally announced September 2021.
-
Probing Neutrino Emission from X-ray Blazar Flares observed with Swift-XRT
Authors:
Stamatios I. Stathopoulos,
Maria Petropoulou,
Paolo Giommi,
Georgios Vasilopoulos,
Paolo Padovani,
Apostolos Mastichiadis
Abstract:
Blazars are the most extreme subclass of active galactic nuclei with relativistic jets emerging from a super-massive black hole and forming a small angle with respect to our line of sight. Blazars are also known to be related to flaring activity as they exhibit large flux variations over a wide range of frequency and on multiple timescales, ranging from a few minutes to several months. The detecti…
▽ More
Blazars are the most extreme subclass of active galactic nuclei with relativistic jets emerging from a super-massive black hole and forming a small angle with respect to our line of sight. Blazars are also known to be related to flaring activity as they exhibit large flux variations over a wide range of frequency and on multiple timescales, ranging from a few minutes to several months. The detection of a high-energy neutrino from the flaring blazar TXS 0506+056 and the subsequent discovery of a neutrino excess from the same direction have naturally strengthened the hypothesis that blazars are cosmic neutrino sources. While neutrino production during gamma-ray flares has been widely discussed, the neutrino yield of X-ray flares has received less attention. Motivated by a theoretical scenario where high energy neutrinos are produced by energetic protons interacting with their own X-ray synchrotron radiation, we make neutrino predictions over a sample of a sample of X-ray blazars. This sample consists of all blazars observed with the X-ray Telescope (XRT) on board Swift more than 50 times from November 2004 to November 2020. The statistical identification of a flaring state is done using the Bayesian Block algorithm to the 1 keV XRT light curves of frequently observed blazars. We categorize flaring states into classes based on their variation from the time-average value of the data points. During each flaring state, we compute the expected muon plus anti-muon neutrino events as well as the total signal for each source using the point-source effective area of Icecube for different operational seasons. We find that the median of the total neutrino number (in logarithm) from flares with duration $<30$ d is $\mathcal{N}^{(\rm tot)}_{ν_μ+\barν_μ} \sim 0.02$.
△ Less
Submitted 30 July, 2021;
originally announced July 2021.
-
Multi-messenger emission from the parsec-scale jet of the flat-spectrum radio quasar PKS 1502+106 coincident with high-energy neutrino IceCube-190730A
Authors:
Foteini Oikonomou,
Maria Petropoulou,
Kohta Murase,
Aaron Tohuvavohu,
Georgios Vasilopoulos,
Sara Buson,
Marcos Santander
Abstract:
On July 30th, 2019 IceCube detected a high-energy astrophysical muon neutrino candidate, IC-190730A, with a $67\%$ probability of astrophysical origin. The flat spectrum radio quasar (FSRQ) PKS 1502+106 is in the error circle of the neutrino. Motivated by this observation, we study PKS 1502+106 as a possible source of IC-190730A. PKS 1502+106 was in a quiet state in terms of UV/optical/X-ray/gamma…
▽ More
On July 30th, 2019 IceCube detected a high-energy astrophysical muon neutrino candidate, IC-190730A, with a $67\%$ probability of astrophysical origin. The flat spectrum radio quasar (FSRQ) PKS 1502+106 is in the error circle of the neutrino. Motivated by this observation, we study PKS 1502+106 as a possible source of IC-190730A. PKS 1502+106 was in a quiet state in terms of UV/optical/X-ray/gamma-ray flux at the time of the neutrino alert, we therefore model the expected neutrino emission from the source during its average long-term state, and investigate whether the emission of IC-190730A as a result of the quiet long-term emission of PKS 1502+106 is plausible. We analyse UV/optical and X-ray data and collect additional observations from the literature to construct the multi-wavelength spectral energy distribution of PKS 1502+106. We perform leptohadronic modelling of the multi-wavelength emission of the source and determine the most plausible emission scenarios and the maximum expected accompanying neutrino flux. A model in which the multi-wavelength emission of PKS 1502+106 originates beyond the broad-line region and inside the dust torus is most consistent with the observations. In this scenario, PKS 1502+106 can have produced up to of order one muon neutrino with energy exceeding 100 TeV in the lifetime of IceCube. An appealing feature of this model is that the required proton luminosity is consistent with the average required proton luminosity if blazars power the observed ultra-high-energy-cosmic-ray flux and well below the source's Eddington luminosity. If such a model is ubiquitous among FSRQs, additional neutrinos can be expected from other bright sources with energy $\gtrsim 10$ PeV.
△ Less
Submitted 3 November, 2021; v1 submitted 23 July, 2021;
originally announced July 2021.
-
The Blazar Hadronic Code Comparison Project
Authors:
Matteo Cerruti,
Michael Kreter,
Maria Petropoulou,
Annika Rudolph,
Foteini Oikonomou,
Markus Böttcher,
Stavros Dimitrakoudis,
Anton Dmytriiev,
Shan Gao,
Apostolos Mastichiadis,
Susumu Inoue,
Kohta Murase,
Anita Reimer,
Joshua Robinson,
Xavier Rodrigues,
Walter Winter,
Andreas Zech,
Natalia Żywucka
Abstract:
Blazar hadronic models have been developed in the past decades as an alternative to leptonic ones. In hadronic models the gamma-ray emission is associated with synchrotron emission by protons, and/or secondary leptons produced in proton-photon interactions. Together with photons, hadronic emission models predict the emission of neutrinos that are therefore the smoking gun for acceleration of relat…
▽ More
Blazar hadronic models have been developed in the past decades as an alternative to leptonic ones. In hadronic models the gamma-ray emission is associated with synchrotron emission by protons, and/or secondary leptons produced in proton-photon interactions. Together with photons, hadronic emission models predict the emission of neutrinos that are therefore the smoking gun for acceleration of relativistic hadrons in blazar jets. The simulation of proton-photon interactions and all associated radiative processes is a complex numerical task, and different approaches to the problem have been adopted in the literature. So far, no systematic comparison between the different codes has been performed, preventing a clear understanding of the underlying uncertainties in the numerical simulations. To fill this gap, we have undertaken the first comprehensive comparison of blazar hadronic codes, and the results from this effort will be presented in this contribution.
△ Less
Submitted 13 July, 2021;
originally announced July 2021.
-
A numerical study of long-term multi-wavelength blazar variability
Authors:
M. Polkas,
M. Petropoulou,
G. Vasilopoulos,
A. Mastichiadis,
M. C. Urry,
P. Coppi,
C. Bailyn
Abstract:
Decade-long monitoring of blazars at optical and infrared (OIR) wavelengths with the Small and Moderate Aperture Research Telescope System (SMARTS) in Chile and in $γ$-rays with the Fermi Large Area Telescope (LAT) has enabled the systematic study of their multi-wavelength long-term variability. In this work we investigate, from a theoretical perspective, the long-term variability properties of bl…
▽ More
Decade-long monitoring of blazars at optical and infrared (OIR) wavelengths with the Small and Moderate Aperture Research Telescope System (SMARTS) in Chile and in $γ$-rays with the Fermi Large Area Telescope (LAT) has enabled the systematic study of their multi-wavelength long-term variability. In this work we investigate, from a theoretical perspective, the long-term variability properties of blazar emission by introducing an observationally motivated time-dependence to four main parameters of the one-zone leptonic model: injection luminosity of relativistic electrons, strength of magnetic field, Doppler factor, and external photon field luminosity. For the first time, we use both the probability density function and the power spectral density of the 10 year-long Fermi-LAT light curves to create variation patterns for the model parameters. Using as test beds two bright blazars from the SMARTS sample (PKS 2155-304 and 3C 273), we compute 10 year-long OIR, X-ray, and $γ$-ray model light curves for different varying parameters. We compare the findings of our theoretical investigation with multi-wavelength observations using various measures of variability. While no single-varying parameter simulation can explain all multi-wavelength variability properties, changes in the electron luminosity and external radiation field in PKS 2155-304 and 3C 273, respectively, can account for most of them. Our results motivate future time-dependent studies with coupling between two or more physical parameters to describe the multi-wavelength long-term blazar variability.
△ Less
Submitted 8 June, 2021; v1 submitted 14 May, 2021;
originally announced May 2021.
-
A marginally fast-cooling proton-synchrotron model for prompt GRBs
Authors:
Ioulia Florou,
Maria Petropoulou,
Apostolos Mastichiadis
Abstract:
A small fraction of GRBs with available data down to soft X-rays ($\sim0.5$ keV) have been shown to feature a spectral break in the low-energy part ($\sim1-10$ keV) of their prompt emission spectrum. The overall spectral shape is consistent with optically thin synchrotron emission from a population of particles that have cooled on a timescale comparable to the dynamic time to energies that are sti…
▽ More
A small fraction of GRBs with available data down to soft X-rays ($\sim0.5$ keV) have been shown to feature a spectral break in the low-energy part ($\sim1-10$ keV) of their prompt emission spectrum. The overall spectral shape is consistent with optically thin synchrotron emission from a population of particles that have cooled on a timescale comparable to the dynamic time to energies that are still much higher than their rest mass energy (marginally fast cooling regime). We consider a hadronic scenario and investigate if the prompt emission of these GRBs can originate from relativistic protons that radiate synchrotron in the marginally fast cooling regime. Using semi-analytical methods, we derive the source parameters, such as magnetic field strength and proton luminosity, and calculate the high-energy neutrino emission expected in this scenario. We also investigate how the emission of secondary pairs produced by photopion interactions and $γγ$ pair production affect the broadband photon spectrum. We support our findings with detailed numerical calculations. Strong modification of the photon spectrum below the break energy due to the synchrotron emission of secondary pairs is found, unless the bulk Lorentz factor is very large ($Γ\gtrsim 10^3$). Moreover, this scenario predicts unreasonably high Poynting luminosities because of the strong magnetic fields ($10^6-10^7$ G) that are necessary for the incomplete proton cooling. Our results strongly disfavour marginally fast cooling protons as an explanation of the low-energy spectral break in the prompt GRB spectra.
△ Less
Submitted 30 April, 2021; v1 submitted 4 February, 2021;
originally announced February 2021.
-
Neutrino signal dependence on gamma-ray burst emission mechanism
Authors:
Tetyana Pitik,
Irene Tamborra,
Maria Petropoulou
Abstract:
Long duration gamma-ray bursts (GRBs) are among the least understood astrophysical transients powering the high-energy universe. To date, various mechanisms have been proposed to explain the observed electromagnetic GRB emission. In this work, we show that, although different jet models may be equally successful in fitting the observed electromagnetic spectral energy distributions, the neutrino pr…
▽ More
Long duration gamma-ray bursts (GRBs) are among the least understood astrophysical transients powering the high-energy universe. To date, various mechanisms have been proposed to explain the observed electromagnetic GRB emission. In this work, we show that, although different jet models may be equally successful in fitting the observed electromagnetic spectral energy distributions, the neutrino production strongly depends on the adopted emission and dissipation model. To this purpose, we compute the neutrino production for a benchmark high-luminosity GRB in the internal shock model, including a dissipative photosphere as well as three emission components, in the jet model invoking internal-collision-induced magnetic reconnection and turbulence (ICMART), in the case of a magnetic jet with gradual dissipation, and in a jet with dominant proton synchrotron radiation. We find that the expected neutrino fluence can vary up to three orders of magnitude in amplitude and peak at energies ranging from $10^4$ to $10^8$ GeV. For our benchmark input parameters, none of the explored GRB models is excluded by the targeted searches carried out by the IceCube and ANTARES Collaborations. However, our work highlights the potential of high-energy neutrinos of pinpointing the underlying GRB emission mechanism and the importance of relying on different jet models for unbiased stacking searches.
△ Less
Submitted 22 April, 2021; v1 submitted 3 February, 2021;
originally announced February 2021.
-
The observability of plasmoid-powered $γ$-ray flares with the Fermi Large Area Telescope
Authors:
Manuel Meyer,
Maria Petropoulou,
Ian Christie
Abstract:
The exact mechanism for the production of fast $γ$-ray variability in blazars remains debated. Magnetic reconnection, in which plasmoids filled with relativistic particles and magnetic fields are formed, is a viable candidate to explain the broadband electromagnetic spectrum and variability of these objects. Using state-of-the-art magnetic reconnection simulations, we generate realistic $γ$-ray li…
▽ More
The exact mechanism for the production of fast $γ$-ray variability in blazars remains debated. Magnetic reconnection, in which plasmoids filled with relativistic particles and magnetic fields are formed, is a viable candidate to explain the broadband electromagnetic spectrum and variability of these objects. Using state-of-the-art magnetic reconnection simulations, we generate realistic $γ$-ray light curves that would be observed with the Fermi Large Area Telescope. A comparison with observed $γ$-ray flares from flat spectrum radio quasars (FSRQs) reveals that magnetic reconnection events lead to comparable flux levels and variability patterns, in particular when the reconnection layer is slightly misaligned with the line of sight. Emission from fast plasmoids moving close to the line of sight could explain fast variability on the time scales of minutes for which evidence has been found in observations of FSRQs. Our results motivate improvements in existing radiative transfer simulations as well as dedicated searches for fast variability as evidence for magnetic reconnection events.
△ Less
Submitted 5 May, 2021; v1 submitted 17 December, 2020;
originally announced December 2020.
-
Hadronic X-ray Flares from Blazars
Authors:
Apostolos Mastichiadis,
Maria Petropoulou
Abstract:
The detection of a high-energy neutrino from the flaring blazar TXS 0506+056 and the subsequent discovery of a neutrino excess from the same direction have strengthened the hypothesis that blazars are cosmic neutrino sources. The lack, however, of $γ$-ray flaring activity during the latter period challenges the standard scenario of correlated $γ$-ray and high-energy neutrino emission in blazars. W…
▽ More
The detection of a high-energy neutrino from the flaring blazar TXS 0506+056 and the subsequent discovery of a neutrino excess from the same direction have strengthened the hypothesis that blazars are cosmic neutrino sources. The lack, however, of $γ$-ray flaring activity during the latter period challenges the standard scenario of correlated $γ$-ray and high-energy neutrino emission in blazars. We propose instead that TeV-PeV neutrinos are produced in coincidence with X-ray flares that are powered by proton synchrotron radiation. In this case, neutrinos are produced by photomeson interactions of protons with their own synchrotron radiation, while MeV to GeV $γ$-rays are the result of synchrotron-dominated electromagnetic cascades developed in the source. Using a time-dependent approach, we find that this "pure hadronic flaring" hypothesis has several interesting consequences. The X-ray flux is a good proxy for the all-flavor neutrino flux, while certain neutrino-rich X-ray flares may be dark in GeV-TeV $γ$-rays. Lastly, hadronic X-ray flares are accompanied by an equally bright MeV component that is detectable by proposed missions like e-ASTROGAM and AMEGO. We then applied this scenario to the extreme blazar 3HSP J095507.9+355101 that has been associated with IceCube-200107A while undergoing an X-ray flare. We showed that the number of muon and antimuon neutrinos above 100 TeV during hadronic flares can be up to $\sim3-10$ times higher than the expected number in standard leptohadronic models. Still, frequent hadronic flaring activity is necessary for explaining the detected neutrino event IceCube-200107A.
△ Less
Submitted 14 November, 2020; v1 submitted 25 September, 2020;
originally announced September 2020.
-
Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas
Authors:
H. Ji,
J. Karpen,
A. Alt,
S. Antiochos,
S. Baalrud,
S. Bale,
P. M. Bellan,
M. Begelman,
A. Beresnyak,
A. Bhattacharjee,
E. G. Blackman,
D. Brennan,
M. Brown,
J. Buechner,
J. Burch,
P. Cassak,
B. Chen,
L. -J. Chen,
Y. Chen,
A. Chien,
L. Comisso,
D. Craig,
J. Dahlin,
W. Daughton,
E. DeLuca
, et al. (83 additional authors not shown)
Abstract:
Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen…
▽ More
Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.
△ Less
Submitted 16 September, 2020;
originally announced September 2020.
-
Secondary Energization in Compressing Plasmoids during Magnetic Reconnection
Authors:
Hayk Hakobyan,
Maria Petropoulou,
Anatoly Spitkovsky,
Lorenzo Sironi
Abstract:
Plasmoids -- magnetized quasi-circular structures formed self-consistently in reconnecting current sheets -- were previously considered to be the graveyards of energetic particles. In this paper, we demonstrate the important role of plasmoids in shaping the particle energy spectrum in relativistic reconnection (i.e., with upstream magnetization $σ_{\rm up} \gg 1$). Using two dimensional particle-i…
▽ More
Plasmoids -- magnetized quasi-circular structures formed self-consistently in reconnecting current sheets -- were previously considered to be the graveyards of energetic particles. In this paper, we demonstrate the important role of plasmoids in shaping the particle energy spectrum in relativistic reconnection (i.e., with upstream magnetization $σ_{\rm up} \gg 1$). Using two dimensional particle-in-cell simulations in pair plasmas with $σ_{\rm up}=10$ and $100$, we study a secondary particle energization process that takes place inside compressing plasmoids. We demonstrate that plasmoids grow in time, while their interiors compress, amplifying the internal magnetic field. The magnetic field felt by particles injected in an isolated plasmoid increases linearly with time, which leads to particle energization as a result of magnetic moment conservation. For particles injected with a power-law distribution function, this energization process acts in such a way that the shape of the injected power law is conserved, while producing an additional non-thermal tail $f(E)\propto E^{-3}$ at higher energies followed by an exponential cutoff. The cutoff energy, which increases with time as $E_{\rm cut}\propto\sqrt{t}$, can greatly exceed $σ_{\rm up} m_e c^2$. We analytically predict the secondary acceleration timescale and the shape of the emerging particle energy spectrum, which can be of major importance in certain astrophysical systems, such as blazar jets.
△ Less
Submitted 4 May, 2021; v1 submitted 22 June, 2020;
originally announced June 2020.
-
Deciphering the properties of the central engine in GRB collapsars
Authors:
M. Petropoulou,
P. Beniamini,
G. Vasilopoulos,
D. Giannios,
R. Barniol Duran
Abstract:
The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic pr…
▽ More
The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity timescale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout timescale on the engine luminosity and the effects of the detector's flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black-hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model.
△ Less
Submitted 12 June, 2020;
originally announced June 2020.
-
High-Energy Neutrino and Gamma-Ray Emission from Tidal Disruption Events
Authors:
Kohta Murase,
Shigeo S. Kimura,
B. Theodore Zhang,
Foteini Oikonomou,
Maria Petropoulou
Abstract:
Tidal disruption events (TDE) have been considered as cosmic-ray and neutrino sources for a decade. We suggest two classes of new scenarios for high-energy multi-messenger emission from TDEs that do not have to harbor powerful jets. First, we investigate high-energy neutrino and gamma-ray production in the core region of a supermassive black hole. In particular, we show that about 1-100 TeV neutri…
▽ More
Tidal disruption events (TDE) have been considered as cosmic-ray and neutrino sources for a decade. We suggest two classes of new scenarios for high-energy multi-messenger emission from TDEs that do not have to harbor powerful jets. First, we investigate high-energy neutrino and gamma-ray production in the core region of a supermassive black hole. In particular, we show that about 1-100 TeV neutrinos and MeV gamma-rays can efficiently be produced in hot coronae around an accretion disk. We also study the consequences of particle acceleration in radiatively inefficient accretion flows (RIAFs). Second, we consider possible cosmic-ray acceleration by sub-relativistic disk-driven winds or interactions between tidal streams, and show that subsequent hadronuclear and photohadronic interactions inside the TDE debris lead to GeV-PeV neutrinos and sub-GeV cascade gamma-rays. We demonstrate that these models should be accompanied by soft gamma-rays or hard X-rays as well as optical/UV emission, which can be used for future observational tests. Although this work aims to present models of non-jetted high-energy emission, we discuss the implications of the TDE AT2019dsg that might coincide with the high-energy neutrino IceCube-191001A, by considering the corona, RIAF, hidden sub-relativistic wind, and hidden jet models. It is not yet possible to be conclusive about their physical association and the expected number of neutrinos is typically much less than unity. We find that the most optimistic cases of the corona and hidden wind models could be consistent with the observation of IceCube-191001A, whereas jet models are unlikely to explain the multi-messenger observations.
△ Less
Submitted 28 August, 2020; v1 submitted 18 May, 2020;
originally announced May 2020.
-
Comprehensive Multimessenger Modeling of the Extreme Blazar 3HSP J095507.9+355101 and Predictions for IceCube
Authors:
Maria Petropoulou,
Foteini Oikonomou,
Apostolos Mastichiadis,
Kohta Murase,
Paolo Padovani,
Georgios Vasilopoulos,
Paolo Giommi
Abstract:
3HSP J095507.9+355101 is an extreme blazar which has been possibly associated with a high-energy neutrino (IceCube-200107A) detected one day before the blazar was found to undergo a hard X-ray flare. We perform a comprehensive study of the predicted multimessenger emission from 3HSP J095507.9+355101 during its recent X-ray flare, but also in the long term. We focus on one-zone leptohadronic models…
▽ More
3HSP J095507.9+355101 is an extreme blazar which has been possibly associated with a high-energy neutrino (IceCube-200107A) detected one day before the blazar was found to undergo a hard X-ray flare. We perform a comprehensive study of the predicted multimessenger emission from 3HSP J095507.9+355101 during its recent X-ray flare, but also in the long term. We focus on one-zone leptohadronic models, but we also explore alternative scenarios: (i) a blazar-core model, which considers neutrino production in the inner jet, close to the supermassive black hole; (ii) a hidden external-photon model, which considers neutrino production in the jet through interactions with photons from a weak broad line region; (iii) a proton synchrotron model, where high-energy protons in the jet produce $γ$-rays via synchrotron; and (iv) an intergalactic cascade scenario, where neutrinos are produced in the intergalactic medium by interactions of a high-energy cosmic-ray beam escaping the jet. The Poisson probability to detect one muon neutrino in ten years from 3HSP J095507.9+355101 with the real-time IceCube alert analysis is $\sim 1\%$ ($3\%$) for the most optimistic one-zone leptohadronic model (the multi-zone blazar-core model). Meanwhile, detection of one neutrino during the 44-day-long high X-ray flux-state period following the neutrino detection is $0.06\%$, according to our most optimistic leptohadronic model. The most promising scenarios for neutrino production also predict strong intra-source $γ$-ray attenuation above $\sim100$ GeV. If the association is real, then IceCube-Gen2 and other future detectors should be able to provide additional evidence for neutrino production in 3HSP J095507.9+355101 and other extreme blazars.
△ Less
Submitted 21 July, 2020; v1 submitted 14 May, 2020;
originally announced May 2020.
-
Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena throughout the Universe
Authors:
H. Ji,
A. Alt,
S. Antiochos,
S. Baalrud,
S. Bale,
P. M. Bellan,
M. Begelman,
A. Beresnyak,
E. G. Blackman,
D. Brennan,
M. Brown,
J. Buechner,
J. Burch,
P. Cassak,
L. -J. Chen,
Y. Chen,
A. Chien,
D. Craig,
J. Dahlin,
W. Daughton,
E. DeLuca,
C. F. Dong,
S. Dorfman,
J. Drake,
F. Ebrahimi
, et al. (75 additional authors not shown)
Abstract:
This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.
This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.
△ Less
Submitted 31 March, 2020;
originally announced April 2020.
-
Proton Synchrotron $γ$-rays and the Energy Crisis in Blazars
Authors:
I. Liodakis,
M. Petropoulou
Abstract:
The origin of high-energy emission in blazars jets (i.e., leptonic versus hadronic) has been a long-standing matter of debate. Here, we focus on one variant of hadronic models where proton synchrotron radiation accounts for the observed steady $γ$-ray blazar emission. Using analytical methods, we derive the minimum jet power ($P_{j,\min}$) for the largest blazar sample analyzed to date (145 source…
▽ More
The origin of high-energy emission in blazars jets (i.e., leptonic versus hadronic) has been a long-standing matter of debate. Here, we focus on one variant of hadronic models where proton synchrotron radiation accounts for the observed steady $γ$-ray blazar emission. Using analytical methods, we derive the minimum jet power ($P_{j,\min}$) for the largest blazar sample analyzed to date (145 sources), taking into account uncertainties of observables and jet's physical parameters. We compare $P_{j,\min}$ against three characteristic energy estimators for accreting systems, i.e., the Eddington luminosity, the accretion disk luminosity, and the power of the Blandford-Znajek process, and find that $P_{j,\min}$ is about 2 orders of magnitude higher than all energetic estimators for the majority of our sample. The derived magnetic field strengths in the emission region require either large amplification of the jet's magnetic field (factor of 30) or place the $γ$-ray production site at sub-pc scales. The expected neutrino emission peaks at $\sim 0.1-10$ EeV, with typical peak neutrino fluxes $\sim 10^{-4}$ times lower than the peak $γ$-ray fluxes. We conclude that if relativistic hadrons are present in blazar jets, they can only produce a radiatively subdominant component of the overall spectral energy distribution of the blazar's steady emission.
△ Less
Submitted 26 April, 2020; v1 submitted 23 March, 2020;
originally announced March 2020.
-
A roadmap to hadronic supercriticalities: a comprehensive study of the parameter space for high-energy astrophysical sources
Authors:
Apostolos Mastichiadis,
Ioulia Florou,
Elina Kefala,
Stella S. Boula,
Maria Petropoulou
Abstract:
Hadronic supercriticalities are radiative instabilities that appear when large amounts of energy are stored in relativistic protons. When the proton energy density exceeds some critical value, a runaway process is initiated resulting in the explosive transfer of the proton energy into electron-positron pairs and radiation. The runaway also leads to an increase of the radiative efficiency, namely t…
▽ More
Hadronic supercriticalities are radiative instabilities that appear when large amounts of energy are stored in relativistic protons. When the proton energy density exceeds some critical value, a runaway process is initiated resulting in the explosive transfer of the proton energy into electron-positron pairs and radiation. The runaway also leads to an increase of the radiative efficiency, namely the ratio of the photon luminosity to the injected proton luminosity. We perform a comprehensive study of the parameter space by investigating the onset of hadronic supercriticalities for a wide range of source parameters (i.e., magnetic field strengths of 1 G$- 100$ kG and radii of $10^{11}$ cm$-10^{16}$ cm) and maximum proton Lorentz factors ($10^3-10^9$). We show that supercriticalities are possible for the whole range of source parameters related to compact astrophysical sources, like gamma-rays bursts, cores and jets of active galactic nuclei. We also provide an in-depth look at the physical mechanisms of hadronic supercriticalities and show that magnetized relativistic plasmas are excellent examples of non-linear dynamical systems in high-energy astrophysics.
△ Less
Submitted 12 June, 2020; v1 submitted 15 March, 2020;
originally announced March 2020.
-
Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013
Authors:
MAGIC collaboration,
V. A. Acciari,
S. Ansoldi,
L. A. Antonelli,
A. Arbet Engels,
D. Baack,
A. Babic,
B. Banerjee,
U. Barres de Almeida,
J. A. Barrio,
J. Becerra Gonzalez,
W. Bednarek,
L. Bellizzi,
E. Bernardini,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
A. Biland,
O. Blanch,
G. Bonnoli,
Z. Bosnjak,
G. Busetto,
R. Carosi,
G. Ceribella
, et al. (215 additional authors not shown)
Abstract:
We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on…
▽ More
We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.
△ Less
Submitted 23 January, 2020;
originally announced January 2020.
-
Ready, set, launch: time interval between BNS merger and short GRB jet formation
Authors:
Paz Beniamini,
Rodolfo Barniol Duran,
Maria Petropoulou,
Dimitrios Giannios
Abstract:
The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents…
▽ More
The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents. Combining this knowledge with the observed luminosities and durations of cosmological sGRBs detected by the Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift Observatory, we revisit the breakout conditions of sGRB jets. Assuming self-collimation of sGRB jets does not play a critical role, we find that the time interval between the binary merger and the launching of a typical sGRB jet is $\lesssim0.1$~s. We also show that for a fraction of at least $\sim 30\%$ of sGRBs, the usually adopted assumption of static ejecta is inconsistent with observations, even if the polar ejecta mass is an order of magnitude smaller than the one in GRB 170817. Our results disfavour magnetar central engines for powering cosmological sGRBs, limit the amount of energy deposited in the cocoon prior to breakout, and suggest that the observed delay of $\sim 1.$7~s in GW 170817 /GRB 170817 between the gravitational wave and $γ$-ray signals is likely dominated by the propagation time of the jet to the $γ$-ray production site.
△ Less
Submitted 18 May, 2020; v1 submitted 3 January, 2020;
originally announced January 2020.
-
Multi-Epoch Modeling of TXS 0506+056 and Implications for Long-Term High-Energy Neutrino Emission
Authors:
Maria Petropoulou,
Kohta Murase,
Marcos Santander,
Sara Buson,
Aaron Tohuvavohu,
Taiki Kawamuro,
Georgios Vasilopoulos,
Hiroshi Negoro,
Yoshihiro Ueda,
Michael H. Siegel,
Azadeh Keivani,
Nobuyuki Kawai,
Apostolos Mastichiadis,
Stavros Dimitrakoudis
Abstract:
The IceCube report of a $\sim 3.5σ$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 05056+056 in 2014-2015 and the 2017 detection of a high-energy neutrino, IceCube-170922A, during a gamma-ray flare from the same blazar, have revived the interest in scenarios for neutrino production in blazars. We perform comprehensive analyses on the long-term electromagnetic emission of TXS…
▽ More
The IceCube report of a $\sim 3.5σ$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 05056+056 in 2014-2015 and the 2017 detection of a high-energy neutrino, IceCube-170922A, during a gamma-ray flare from the same blazar, have revived the interest in scenarios for neutrino production in blazars. We perform comprehensive analyses on the long-term electromagnetic emission of TXS 05056+056 using optical, X-ray, and gamma-ray data from the All-Sky Automated Survey for Supernovae (ASAS-SN), the Neil Gehrels Swift Observatory (Swift), the Monitor of All-sky X-ray Image (MAXI), and the Fermi Large Area Telescope (Fermi-LAT). We also perform numerical modeling of the spectral energy distributions (SEDs) in four epochs prior to 2017 with contemporaneous gamma-ray and lower energy (optical and/or X-ray) data. We find that the multi-epoch SEDs are consistent with a hybrid leptonic scenario, where the gamma-rays are produced in the blazar zone via external inverse Compton scattering of accelerated electrons, and high-energy neutrinos are produced via the photomeson production process of co-accelerated protons. The multi-epoch SEDs can be satisfactorily explained with the same jet parameters and variable external photon density and electron luminosity. Using the maximal neutrino flux derived for each epoch, we put an upper limit of $\sim0.4-2$ on the muon neutrino number in ten years of IceCube observations. Our results are consistent with the IceCube-170922A detection, which can be explained as an upper fluctuation from the average neutrino rate expected from the source, but in strong tension with the 2014-2015 neutrino flare.
△ Less
Submitted 16 February, 2020; v1 submitted 10 November, 2019;
originally announced November 2019.
-
Inverse Compton Signatures of Gamma-Ray Burst Afterglows
Authors:
Hao Zhang,
Ian Christie,
Maria Petropoulou,
Jesus M. Rueda-Becerril,
Dimitrios Giannios
Abstract:
The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~G…
▽ More
The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~GeV) remains uncertain. The recent detection of sub-TeV emission from GRB~190114C by MAGIC raises further debate on what powers the very high-energy (VHE; $\gtrsim 300$GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multi-wavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB~190114C, we find that its afterglow emission in the \fermi-LAT band is synchrotron-dominated.The late-time \fermi-LAT measurement (i.e., $t\sim 10^4$~s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. $\lesssim 3\times 10^{-9}\,{\rm erg\,cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies.
△ Less
Submitted 17 June, 2020; v1 submitted 30 October, 2019;
originally announced October 2019.