Lepto-hadronic Model Research Articles

Model for Population III microquasars

Context.Current simulations indicate that the first stars were formed predominantly in binary systems. Studying the contribution of the first accreting binary systems to the reionization and heating of the intergalactic medium requires formulating a concrete model for Population III microquasars.Aims.We aim at constructing a complete model for microquasars with a Population III donor star.Methods.We considered that the mass loss of the star is exclusively caused by the spill of matter through the Roche lobe toward the black hole. We calculated the spectral energy distribution of the radiation produced by the accretion disk, the radiation-pressure-driven wind, and the relativistic particles in the jets, within the framework of a lepto-hadronic model. In addition, we estimated the effect on the reionization by the first microquasars.Results.We determine that Population III microquasars are powerful sources of ultraviolet radiation produced by the winds of their super-critical disks, and that they generate a broadband nonthermal emission in jets.Conclusions.Our results indicate that microquasars in the early Universe could have been important for the reionization and heating of the intergalactic medium.

Leptohadronic Blazar Models Applied to the 2014–2015 Flare of TXS 0506+056

We investigate whether the emission of neutrinos observed in 2014–2015 from the direction of the blazar TXS 0506+056 can be accommodated with leptohadronic multiwavelength models of the source commonly adopted for the 2017 flare. While multiwavelength data during the neutrino flare are sparse, the large number of neutrino events (13 ± 5) challenges the missing activity in gamma-rays. We illustrate that two to five neutrino events during the flare can be explained with leptohadronic models of different categories: a one-zone model, a compact-core model, and an external radiation field model. If, however, significantly more events were to be accommodated, the predicted multiwavelength emission levels would be in conflict with observational X-ray constraints, or with the high-energy gamma-ray fluxes observed by the Fermi Large Area Telescope, depending on the model. For example, while the external radiation field model can predict up to five neutrino events without violating X-ray constraints, the absorption of high-energy gamma-rays is in minor tension with the data. We therefore do not find any model that can simultaneously explain the high event number quoted by IceCube and the (sparse) electromagnetic data during the neutrino flare.

Detection of a flaring blazar coincident with an IceCube high-energy neutrino

In September 22, 2017, IceCube released a public alert announcing the detection of a 290 TeV neutrino track event with an angular uncertainty of one square degree (90% containment). A multi-messenger follow-up campaign was initiated resulting in the detection of a GeV gamma-ray flare by the Fermi Large Area Telescope positionally consistent with the location of the known Bl Lac object, TXS 0506+056 , located only 0.1 degrees from the best-fit neutrino position. The probability of finding a GeV gamma-ray flare in coincidence with a high-energy neutrino event assuming a correlation of the neutrino flux with the gamma-ray energy flux in the energy band between 1 and 100 GeV was calculated to be 3σ (after trials correction). Following the detection of the flaring blazar the imaging air Cherenkov telescope MAGIC detected the source for the first time in the > 100 GeV gamma-ray band. The activity of the source was confirmed in X-ray, optical and radio wavelength. Several groups have developed lepto-hadronic models which succeed to explain the multi-messenger spectral energy distribution.

The Blazar TXS 0506+056 Associated with a High-energy Neutrino: Insights into Extragalactic Jets and Cosmic-Ray Acceleration

Abstract A neutrino with energy ∼290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS 0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at ∼3σ level. We monitored the object in the very-high-energy (VHE) band with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes for ∼41 hr from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of . The steep spectrum observed by MAGIC is concordant with internal γγ absorption above ∼100 GeV entailed by photohadronic production of a ∼290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton upscattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet comoving frame can be in the range ∼1014 – 1018 eV.

Population III microquasars

We present the first results obtained in the elaboration of a complete model of a microquasar where the donor star is from Population III. These stars do not produce stellar winds so we consider that the mass loss is due exclusively to matter overflowing the Roche lobe towards the compact object, a maximally rotating black hole. The rate of accretion is extremely super-Eddington, with an intense mass loss from the system in the form of winds and jets. We calculate the relativistic particle content of the jet and the corresponding spectral energy distribution (SED) considering a lepto-hadronic model. Prospects for the cosmological implications of these objects are briefly discussed.

Is the GeV-TeV emission of PKS 0447-439 from the proton synchrotron radiation?

We study the multi-wavelength emission features of PKS 0447-439 in the frame of the one-zone homogeneous lepto-hadronic model. In this model, we assumed that the steady power-laws with exponential cut-offs distributions of protons and electrons are injected into the source. The non-linear time-dependent kinematic equations, describing the evolution of protons, electrons and photons, are defined; these equations self-consistently involve synchrotron radiation of protons, photon-photon interaction, synchrotron radiation of electron/positron pairs, inverse Compton scattering and synchrotron self-absorption. The model is applied to reproduce the multi-wavelength spectrum of PKS 0447-439. Our results indicate that the spectral energy distribution (SED) of PKS 0447-439 can be reproduced well by the model. In particular, the GeV-TeV emission is produced by the synchrotron radiation of relativistic protons. The physically plausible solutions require the magnetic strength $10~\text{G}\lesssim B \lesssim 100~\text{G}$ . We found that the observed spectrum of PKS 0447-439 can be reproduced well by the model whether $z = 0.16$ or $z = 0.2$ , and the acceptable upper limit of redshift is $z=0.343$ .

Lepto-hadronic model of gamma rays from Eta Carinae and prospects for neutrino telescopes

The stellar binary $\eta$ Carinae has been observed during its full orbital period in gamma rays by the Fermi-Large Area Telescope (LAT). The shock-accelerated electrons in the colliding winds of the two stars radiate synchrotron photons in the magnetic field of the shocked region and inverse Compton photons, where the target photons are from the thermal emissions by the more massive and luminous of the two stars. The inverse Compton emission dominates the gamma-ray flux data from the $\eta$ Carinae, however the spectral energy distribution shows signature of a hadronic component in the $\sim 10$-300 GeV range during the periastron passage. Current and future air Cherenkov telescopes will be able to constrain this component at TeV energies. Acceleration of cosmic-ray protons to $\gg 1$ TeV energies in the colliding winds, required to explain the hadronic emission component through photopion interactions, can lead to detectable signal of $\gtrsim 10$ TeV neutrino events in large kilometer scale neutrino telescopes.

On the Direct Correlation between Gamma-Rays and PeV Neutrinos from Blazars

Abstract We study the frequently used assumption in multi-messenger astrophysics that the gamma-ray and neutrino fluxes are directly connected because they are assumed to be produced by the same photohadronic production chain. An interesting candidate source for this test is the flat-spectrum radio quasar PKS B1424-418, which recently called attention to a potential correlation between an IceCube PeV neutrino event and its burst phase. We simulate both the multi-waveband photon and the neutrino emission from this source using a self-consistent radiation model. We demonstrate that a simple hadronic model cannot adequately describe the spectral energy distribution for this source, but a lepto-hadronic model with a subdominant hadronic component can reproduce the multi-waveband photon spectrum observed during various activity phases of the blazar. As a conclusion, up to about 0.3 neutrino events may coincide with the burst, which implies that the leptonic contribution dominates in the relevant energy band. We also demonstrate that the time-wise correlation between the neutrino event and burst phase is weak.

3C 279 IN OUTBURST IN 2015 JUNE: A BROADBAND SED STUDY BASED ON THE INTEGRAL DETECTION

ABSTRACT Blazars radiate from radio through gamma-ray frequencies and thereby make ideal targets for multifrequency studies. Such studies allow the properties of the emitting jet to be constrained. 3C 279 is among the most notable blazars and therefore subject to extensive multifrequency campaigns. We report the results of a campaign ranging from near-IR to gamma-ray energies that targeted an outburst of 3C 279 in 2015 June. The campaign pivots around the detection in only 50 ks by INTEGRAL, whose IBIS/ISGRI data pin down the high-energy component of the spectral energy distribution (SED) between Swift-XRT data and Fermi-LAT data. The overall SED from near-IR to gamma rays can be well represented by either a leptonic or a lepto-hadronic radiation transfer model. Even though the data are equally well represented by the two models, their inferred parameters challenge the physical conditions in the jet. In fact, the leptonic model requires parameters with a magnetic field far below equipartition with the relativistic particle energy density. In contrast, equipartition may be achieved with the lepto-hadronic model, although this implies an extreme total jet power close to the Eddington luminosity.

NEUTRINO, γ-RAY, AND COSMIC-RAY FLUXES FROM THE CORE OF THE CLOSEST RADIO GALAXIES

ABSTRACT The closest radio galaxies; Centaurus A (Cen A), M87, and NGC 1275, have been detected from radio wavelengths to TeV γ-rays, and also studied as high-energy neutrino and ultra-high-energy cosmic-ray (UHECR) potential emitters. Their spectral energy distributions (SEDs) show a double-peak feature, which is explained by a synchrotron self-Compton (SSC) model. However, TeV γ-ray measured spectra could suggest that very-high-energy γ-rays might have a hadronic origin. We introduce a lepto-hadronic model to describe the broadband SED; from radio to sub-GeV photons as synchrotron SSC emission and TeV γ-ray photons as neutral pion decay resulting from pγ interactions occurring close to the core. These photo-hadronic interactions take place when Fermi-accelerated protons interact with the seed photons around synchrotron SSC peaks. Obtaining a good description of the TeV γ-ray fluxes, first, we compute neutrino fluxes and events expected in the IceCube detector and, second, we estimate UHECR fluxes and the event rate expected in Telescope Array, Pierre Auger, and HiRes observatories. Within this scenario, we show that the expected high-energy neutrinos cannot explain the astrophysical flux observed by IceCube, and the connection with UHECRs observed by Auger experiment around Cen A might be possible only considering a heavy nuclei composition in the observed events.

Spectral and Polarization Signatures of Relativistic Shocks in Blazars

Relativistic shocks are one of the most plausible sites of the emission of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars, via the diffusive shock acceleration (DSA) of relativistic particles. This paper summarizes recent results on a self-consistent coupling of diffusive shock acceleration and radiation transfer in blazar jets. We demonstrate that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of hydromagnetic turbulence responsible for pitch-angle scattering by requiring a strongly energy-dependent pitch-angle mean free path. The prominent soft X-ray excess (“Big Blue Bump”) in the SED of the BL Lac object AO 0235+164 can be modelled as the signature of bulk Compton scattering of external radiation fields by the thermal electron population, which places additional constraints on the level of hydromagnetic turbulence. It has further been demonstrated that internal shocks propagating in a jet pervaded by a helical magnetic field naturally produce polarization-angle swings by 180 ∘ , in tandem with multi-wavelength flaring activity, without requiring any helical motion paths or other asymmetric jet structures. The specific application of this model to 3C279 presents the first consistent simultaneous modeling of snap-shot SEDs, multi-wavelength light curves, and time-dependent polarization signatures of a blazar during a polarization-angle (PA) rotation. This model has recently been generalized to a lepto-hadronic model, in which the high-energy emission is dominated by proton synchrotron radiation. It is shown that in this case, the high-energy (X-ray and γ-ray) polarization signatures are expected to be significantly more stable (not showing PA rotations) than the low-energy (electron-synchrotron) signatures.

LEPTONIC AND LEPTO-HADRONIC MODELING OF THE 2010 NOVEMBER FLARE FROM 3C 454.3

In this study, we use a one-zone leptonic and a lepto-hadronic model to investigate the multi-wavelength emission and the prominent flare of the flat spectrum radio quasar 3C 454.3 in Nov 2010. We perform a parameter study with both models to obtain broadband fits to the spectral energy distribution of 3C 454.3. Starting with the baseline parameters obtained from the fits, we then investigate different flaring scenarios for both models to explain an extreme outburst and spectral hardening of 3C 454.3 that occurred in Nov 2010. We find that the one zone lepto-hadronic model can successfully explain both the broadband multi-wavelength spectral energy distribution and light curves in the optical R, Swift XRT and Fermi gamma-ray bandpasses for 3C 454.3 during quiescence and the peak of the Nov. 2010 flare. We also find that the one-zone leptonic model produces poor fits to the broadband spectra in the X-ray and high energy gamma-ray band passes for the Nov. 2010 flare.

On the origin of the soft photons of the high-synchrotron-peaked blazar PKS 1424+240

PKS 1424+240 is a distant very high energy gamma-ray BL Lac object with redshift $z=0.601$. It was found that pure synchrotron self-Compton (SSC) process normally need extreme input parameters (e.g., very low magnetic field intensity and extraordinarily large Doppler factor) to explain its multi-wavelength spectral energy distributions (SEDs). To avoid the extreme model parameters, different models have been proposed (e.g., two-zone SSC model or lepto-hadronic model). In this work, we employ the traditional one-zone leptonic model after including a weak external Compton component to re-explore the simultaneous multi-wavelength SEDs of PKS 1424+240 in both high (2009) and low (2013) states. We find that the input parameters of magnetic field and Doppler factor are roughly consistent with those of other BL Lacs if a weak external photon field from either broad line region (BLR) or the dust torus. However, the required energy density of seed photons from BLR or torus is about 3 orders of magnitude less than that constrained in luminous quasars (e.g., flat-spectrum radio quasars, FSRQs). This result suggests that the BLR/torus in BL Lacs is much weaker than that of luminous FSRQs (but not fully disappear), and the inverse-Compton of external photons from BLR/torus may still play a role even in high synchrotron peaked blazars.

Lepto-hadronic model for the broadband emission of Cygnus X-1

Cygnus X-1 is a well observed microquasar. Broadband observations at all wavelengths have been collected over the years. The origin of the MeV tail observed with COMPTEL and INTEGRAL is still under debate and it has mostly been attributed to the corona, although its high degree of polarization suggests it is synchrotron radiation from a jet. The origin of the transient emission above $\sim 100$ GeV is also unclear. We aim to disentangle the origin of the broadband spectral energy distribution (SED) of Cygnus X-1, focusing particularly on the gamma-ray emission, and to gain information on the physical conditions inside the jets. We develop and apply a lepto-hadronic, inhomogeneous jet model to the non-thermal SED of Cygnus X-1. We calculate the contributions to the SED of both protons and electrons accelerated in an extended region of the jet. We also estimate the radiation of charged secondaries produced in hadronic interactions, through several radiative processes. Absorption effects are considered. We produce synthetic maps of the jets at radio wavelengths. We find two sets of model parameters that lead to good fits of the SED. One of the models fits all the observations, including the MeV tail. This model also predicts hadronic gamma-ray emission slightly below the current upper limits. The flux predicted at 8.4 GHz is in agreement with the observations available in the literature, although the synthetic source is more compact than the imaged radio jet. Our results show that the MeV emission in Cygnus X-1 may be jet synchrotron radiation. This depends mainly on the strength of the jet magnetic field and the location of the injection region of the relativistic particles. Our calculations show that there must be energetic electrons in the jets quite far from the black hole.

TIME DEPENDENT HADRONIC MODELING OF FLAT SPECTRUM RADIO QUASARS

We introduce a new time-dependent lepto-hadronic model for blazar emission that takes into account the radiation emitted by secondary particles, such as pions and muons, from photo hadronic interactions. Starting from a baseline parameter set guided by a fit to the spectral energy distribution of the blazar 3C 279, we perform a parameter study to investigate the effects of perturbations of the input parameters to mimic different flaring events to study the resulting lightcurves in the optical, X-ray, high energy (HE: E > 100 MeV) and very-high-energy (VHE: E > 100 GeV) gamma-rays as well as the neutrino emission associated with charged-pion and muon decay. We find that flaring events from an increase in the efficiency of Fermi II acceleration will produce a positive correlation between all bandpasses and a marked plateau in the HE gamma-ray lightcurve. We also predict a distinctive dip in the HE lightcurve for perturbations caused by a change in the proton injection spectral index. These plateaus / dips could be a tell tale signature of hadronic models for perturbations that lead to more efficient acceleration of high energy protons in parameter regimes where pion and muon synchrotron emission is non-negligible.

Photohadronic origin of $\boldsymbol {\gamma }$-ray BL Lac emission: implications for IceCube neutrinos

The recent IceCube discovery of 0.1-1 PeV neutrinos of astrophysical origin opens up a new era for high-energy astrophysics. Although there are various astrophysical candidate sources, a firm association of the detected neutrinos with one (or more) of them is still lacking. A recent analysis of plausible astrophysical counterparts within the error circles of IceCube events showed that likely counterparts for nine of the IceCube neutrinos include mostly BL Lacs, among which Mrk 421. Motivated by this result and a previous independent analysis on the neutrino emission from Mrk 421, we test the BL Lac-neutrino connection in the context of a specific theoretical model for BL Lac emission. We model the spectral energy distribution (SED) of the BL Lacs selected as counterparts of the IceCube neutrinos using a one-zone leptohadronic model and mostly nearly simultaneous data. The neutrino flux for each BL Lac is self-consistently calculated, using photon and proton distributions specifically derived for every individual source. We find that the SEDs of the sample, although different in shape and flux, are all well fitted by the model using reasonable parameter values. Moreover, the model-predicted neutrino flux and energy for these sources are of the same order of magnitude as those of the IceCube neutrinos. In two cases, namely Mrk 421 and H 1914-194, we find a suggestively good agreement between the model prediction and the detected neutrino flux. Our predictions for all the BL Lacs of the sample are in the range to be confirmed or disputed by IceCube in the next few years of data sampling.

Bethe–Heitler emission in BL Lacs: filling the gap between X-rays and γ-rays

We present the spectral signatures of the Bethe-Heitler pair production ($pe$) process on the spectral energy distribution (SED) of blazars, in scenarios where the hard $\gamma$-ray emission is of photohadronic origin. If relativistic protons interact with the synchrotron blazar photons producing $\gamma$ rays through photopion processes, we show that, besides the $2-20$ PeV neutrino emission, the typical blazar SED should have an emission feature due to the synchrotron emission of $pe$ secondaries that bridges the gap betweeen the low-and high-energy humps of the SED, namely in the energy range 40 keV-40 MeV. We first present analytical expressions for the photopion and $pe$ loss rates in terms of observable quantities of blazar emission. For the $pe$ loss rate in particular, we derive a new approximate analytical expression for the case of a power-law photon distribution, which has an excellent accuracy with the numerically calculated exact one, especially at energies above the threshold for pair production. We show that for typical blazar parameters, the photopair synchrotron emission emerges in the hard X-ray/soft $\gamma$-ray energy range with a characteristic spectral shape and non negligible flux, which may be even comparable to the hard $\gamma$-ray flux produced through photopion processes. We argue that the expected "$pe$ bumps" are a natural consequence of leptohadronic models, and as such, they may indicate that blazars with a three-hump SED are possible emitters of high-energy neutrinos.

Lepton-hadronic processes and high-energy neutrinos in NGC 1275

AbstractThe nearby active galaxy NGC 1275, has been widely detected from radio to gamma rays. Its spectral energy distribution (SED) shows a double-peak feature, which is well explained by synchrotron self-Compton (SSC) model. However, recent TeV detections might suggest that very-high-energy γ-rays (E⩾100 GeV) may not have a leptonic origin. We test a lepto-hadronic model to describe the whole SED through SSC emission and neutral pion decay resulting from pγ interactions. Also, we estimate the neutrino events expected in a km3 Cherenkov telescope.

Origin of gamma-ray emission in the shell of Cassiopeia A

Non-thermal X-ray emission from the shell of Cassiopeia A (Cas A) has been an interesting subject of study, as it provides information about relativistic electrons and their acceleration mechanisms in the shocks. Chandra X-ray observatory revealed the detailed spectral and spatial structure of this SNR in X-rays. The spectral analysis of Chandra X-ray data of Cas A shows unequal flux levels for different regions of the shell, which can be attributed to different magnetic fields in those regions. Additionally, the GeV gamma-ray emission observed by Large Area Telescope on board Fermi Gamma Ray Space Telescope showed that the hadronic processes are dominating in Cas A, a clear signature of acceleration of protons. In this paper we aim to explain the GeV-TeV gamma-ray data in the context of both leptonic and hadronic scenario. We modeled the multi-wavelength spectrum of Cas A. We use synchrotron emission process to explain the observed non-thermal X-ray fluxes from different regions of the shell. These result in estimation of the model parameters, which are then used to explain TeV gamma-ray emission spectrum. We also use hadronic scenario to explain both GeV and TeV fluxes simultaneously. We show that a leptonic model alone cannot explain the GeV-TeV data. Therefore, we need to invoke a hadronic model to explain the observed GeV-TeV fluxes. We found that although pure hadronic model is able to explain the GeV-TeV data, a lepto-hadronic model provides the best fit to the data.

NEUTRINO AND UHECR SPECTRA FROM MRK 421

We present the neutrino and UHECR spectra obtained from a detailed fitting of the spectral energy distribution (SED) of Mrk 421 (March 2001) using two variations of the leptohadronic model. In particular, while the low-energy component (optical to X-rays) of the SED is fitted by synchrotron emission of primary electrons in both models, the high-energy one (GeV-TeV gamma-rays) is synchrotron emission attributed either to ultra-high energy protons (LHs model) or to secondary electrons produced by the decay of charged pions (LHπ model). In the LHπ case we find that the produced neutrino spectra are sharply peaked at Eν ~ 30 PeV with a peak flux slightly below the IC-40 sensitivity limit for Mrk 421. In the LHs model, on the other hand, the neutrino spectra fall well outside the PeV energy range, but the calculated E ~ 30 EeV — UHECR flux at earth is close to that observed by HiresI, Telescope Array and Pierre Augere experiments.