Source author record

Maria Petropoulou

Maria Petropoulou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

preprint2022arXiv

Model selection for component network meta-analysis in connected and disconnected networks: a simulation study

Network meta-analysis (NMA) is widely used in evidence synthesis to estimate the effects of several competing interventions for a given clinical condition. One of the challenges is that it is not possible in disconnected networks. Component network meta-analysis (CNMA) allows technically 'reconnecting' a disconnected network with multicomponent interventions. The additive CNMA model assumes that the effect of any multicomponent intervention is the additive sum of its components. This assumption can be relaxed by adding interaction component terms, which improves the goodness of fit but decreases the network connectivity. Model selection aims at finding the model with a reasonable balance between the goodness of fit and connectivity (selected CNMA model). We aim to introduce a forward model selection strategy for CNMA models and to investigate the performance of CNMA models for connected and disconnected networks. We applied the methods to a real Cochrane review dataset and simulated data with additive, mildly, or strongly violated intervention effects. We started with connected networks, and we artificially constructed disconnected networks. We compared the results of the additive and the selected CNMAs from each connected and disconnected network with the NMA using the mean squared error and coverage probability. CNMA models provide good performance for connected networks and can be an alternative to standard NMA if additivity holds. On the contrary, model selection does not perform well for disconnected networks, and we recommend conducting separate analyses of subnetworks.

preprint2021arXiv

A unified model for orphan and multi-wavelength blazar flares

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.

preprint2021arXiv

Circumnuclear Dust in AP Librae and the source of its VHE emission

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.

preprint2021arXiv

Multi-messenger emission from the parsec-scale jet of the flat-spectrum radio quasar PKS 1502+106 coincident with high-energy neutrino IceCube-190730A

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.

preprint2020arXiv

A roadmap to hadronic supercriticalities: a comprehensive study of the parameter space for high-energy astrophysical sources

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.

preprint2020arXiv

Comprehensive Multimessenger Modeling of the Extreme Blazar 3HSP J095507.9+355101 and Predictions for IceCube

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.

preprint2020arXiv

Hadronic X-ray Flares from Blazars

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.

preprint2020arXiv

Inverse Compton Signatures of Gamma-Ray Burst Afterglows

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.

preprint2020arXiv

Multi-Epoch Modeling of TXS 0506+056 and Implications for Long-Term High-Energy Neutrino Emission

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.

preprint2020arXiv

Ready, set, launch: time interval between BNS merger and short GRB jet formation

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.

preprint2019arXiv

A Neutral Beam Model for High-Energy Neutrino Emission from the Blazar TXS 0506+056

The IceCube collaboration reported a $\sim 3.5σ$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 0506+56 during a $\sim$6 month period in 2014-2015, as well as the ($\sim3σ$) detection of a high-energy muon neutrino during an electromagnetic flare in 2017. We explore the possibility that the 2014-2015 neutrino excess and the 2017 multi-messenger flare are both explained in a common physical framework that relies on the emergence of a relativistic neutral beam in the blazar jet due to interactions of accelerated cosmic rays (CRs) with photons. We demonstrate that the neutral beam model provides an explanation for the 2014-2015 neutrino excess without violating X-ray and $γ$-ray constraints, and also yields results consistent with the detection of one high-energy neutrino during the 2017 flare. If both neutrino associations with TXS 05065+056 are real, our model requires that (i) the composition of accelerated CRs is light, with a ratio of helium nuclei to protons $\gtrsim5$, (ii) a luminous external photon field ($\sim 10^{46}$ erg s$^{-1}$) variable (on year-long timescales) is present, and (iii) the CR injection luminosity as well as the properties of the dissipation region (i.e., Lorentz factor, magnetic field, and size) vary on year-long timescales.

preprint2019arXiv

A two-zone model for blazar emission: implications for TXS 0506+056 and the neutrino event IceCube-170922A

A high-energy muon neutrino event, IceCube-170922A, was recently discovered in both spatial and temporal coincidence with a gamma-ray flare of the blazar TXS 0506+056. It has been shown, with standard one-zone models, that neutrinos can be produced in the blazar jet via hadronic interactions, but with a flux which is mostly limited by the X-ray data. In this work, we explore the neutrino production from TXS 0506+056 by invoking two physically distinct emission zones in the jet, separated by the broad line region (BLR). Using the Doppler-boosted radiation of the BLR as the target photon field, the inner zone accounts for the neutrino and gamma-ray emission via $pγ$ interactions and inverse Compton scattering respectively, while the outer zone produces the optical and X-ray emission via synchrotron and synchrotron self-Compton processes. The different conditions of the two zones allow us to suppress the X-ray emission from the electromagnetic cascade, and set a much higher upper limit on the muon neutrino flux (i.e., $\sim 10^{-11}\rm erg~cm^{-2}s^{-1}$) than in one-zone models. We compare, in detail, our scenario with one-zone models discussed in the literature, and argue that differentiating between such scenarios will become possible with next generation neutrino telescopes, such as IceCube-Gen2.

preprint2016arXiv

Blazar flares powered by plasmoids in relativistic reconnection

Powerful flares from blazars with short ($\sim$ min) variability timescales are challenging for current models of blazar emission. Here, we present a physically motivated ab initio model for blazar flares based on the results of recent particle-in-cell (PIC) simulations of relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids filled with high-energy particles and magnetic fields are a self-consistent by-product of the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration profile, particle and magnetic content) with a kinetic equation for the evolution of the electron distribution function we demonstrate that relativistic reconnection in blazar jets can produce powerful flares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron self-Compton flares of duration of several hours--days powered by the largest and slowest moving plasmoids that form in the reconnection layer. Smaller and faster plasmoids produce flares of sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed fluence in both types of flares is similar. Multiple flares with a range of flux-doubling timescales (minutes to several hours) observed over a longer period of flaring activity (days or longer) may be used as a probe of the reconnection layer's orientation and the jet's magnetization. Our model shows that blazar flares are naturally expected as a result of magnetic reconnection in a magnetically-dominated jet.

preprint2016arXiv

Early-time signatures of γ-ray emission from supernovae in dense circumstellar media

We present our results on the γ-ray emission from interaction-powered supernovae (SNe), a recently discovered SN type that is suggested to be surrounded by a circumstellar medium (CSM) with densities 10^7-10^12~ cm^-3. Such high densities favor inelastic collisions between relativistic protons accelerated in the SN blast wave and CSM protons and the production of γ-ray photons through neutral pion decays. Using a numerical code that includes synchrotron radiation, adiabatic losses due to the expansion of the source, photon-photon interactions, proton-proton collisions and proton-photon interactions, we calculate the multi-wavelength non-thermal photon emission soon after the shock breakout and follow its temporal evolution until 100-1000 days. Focusing on the γ-ray emission at >100 MeV, we show that this could be detectable by the Fermi-LAT telescope for nearby (<10 Mpc) SNe with dense CSM (>10^11 cm^-3).

preprint2016arXiv

Radio synchrotron emission from secondary electrons in interaction-powered supernovae

Several supernovae (SNe) with an unusually dense circumstellar medium (CSM) have been recently observed at radio frequencies. Their radio emission is powered by relativistic electrons that can be either accelerated at the SN shock (primaries) or injected as a by-product (secondaries) of inelastic proton-proton collisions. We investigate the radio signatures from secondary electrons, by detailing a semi-analytical model to calculate the temporal evolution of the distributions of protons, primary and secondary electrons. With our formalism, we track the cooling history of all the particles that have been injected into the emission region up to a given time, and calculate the resulting radio spectra and light curves. For a SN shock propagating through the progenitor wind, we find that secondary electrons control the early radio signatures, but their contribution decays faster than that of primary electrons. This results in a flattening of the light curve at a given radio frequency that depends only upon the radial profiles of the CSM density and of the shock velocity, $v_0$. The relevant transition time at the peak frequency is $\sim 190 \, K_{\rm ep,-3}^{-1} A_{\rm w, 16} β_{0, -1.5}^{-2}\,{\rm d}$, where $A_{\rm w}$ is the wind mass-loading parameter, $β_0=v_0/c$ and $K_{\rm ep}$ is the electron-to-proton ratio of accelerated particles. We explicitly show that late peak times at 5 GHz (i.e., $t_{\rm pk} \gtrsim 300-1000$ d) suggest a shock wave propagating in a dense wind ($A_{\rm w} \gtrsim 10^{16}-10^{17}$ gr cm$^{-1}$), where secondary electrons are likely to power the observed peak emission.

preprint2016arXiv

The many faces of blazar emission in the context of hadronic models

We present two ways of modeling the spectral energy distribution of blazars in the hadronic context and discuss the predictions of each "hadronic variant" on the spectral shape, the multi-wavelength variability, the cosmic-ray flux, and the high-energy neutrino emission. Focusing on the latter, we then present an application of the hadronic model to individual BL Lacs that were recently suggested to be the counterparts of some of the IceCube neutrinos, and conclude by discussing the contribution of the whole BL Lac class to the observed neutrino background.

preprint2016arXiv

The TeV emission of Ap Librae: a hadronic interpretation and prospects for CTA

Ap Librae is one out of a handful of low-frequency peaked blazars to be detected at TeV $γ$-rays and the only one with an identified X-ray jet. Combined observations of Fermi-LAT at high energies (HE) and of H.E.S.S at very high energies (VHE) revealed a striking spectral property of Ap Librae; the presence of a broad high-energy component that extends more than nine orders of magnitude in energy and is, therefore, hard to be explained by the usual single-zone synchrotron self-Compton model. We show that the superposition of different emission components related to photohadronic interactions can explain the $γ$-ray emission of Ap Librae without invoking external radiation fields. We present two indicative model fits to the spectral energy distribution of Ap Librae where the VHE emission is assumed to originate from a compact, sub-pc scale region of the jet. A robust prediction of our model is VHE flux variability on timescales similar to those observed at X-rays and HE $γ$-rays, which can be further used to distinguish between a sub-pc or kpc scale origin of the TeV emission. We thus calculate the expected variability signatures at X-rays, HE and VHE $γ$-rays and show that quasi-simultaneous flares are expected, with larger amplitude flares appearing at $γ$-rays. We assess the detectability of VHE variability from Ap Librae with CTA, next generation of IACTs. We show that $\sim$hr timescale variability at $E_γ>0.1$ TeV could be detectable at high significance with shorter exposure times than current Cherenkov telescopes.

preprint2016arXiv

Time-dependent neutrino emission from Mrk 421 during flares and predictions for IceCube

Blazars are prime candidate sources for the high energy neutrinos recently detected by IceCube. Being one of the brightest sources in the extragalactic X-ray and $γ$-ray sky as well as one of the nearest blazars to Earth, Mrk 421 is an excellent source for testing the scenario of the blazar-neutrino connection. Here, we model the spectral energy distribution of Mrk 421 during a 13-day flare in 2010 with unprecedented multi-wavelength coverage, and calculate the respective neutrino flux. We find a correlation between the $>1$ PeV neutrino and photon fluxes, in all energy bands. Using typical IceCube through-going muon event samples with good angular resolution and high statistics, we derive the mean event rate above 100 TeV ($\sim0.57$ evt/yr) and show that it is comparable to that expected from a four-month quiescent period in 2009. Due to the short duration of the flare, an accumulation of similar flares over several years would be necessary to produce a meaningful signal for IceCube. To better assess this, we apply the correlation between the neutrino and $γ$-ray fluxes to the 6.9 yr Fermi-LAT light curve of Mrk 421. We find that the mean event count above 1 PeV for the full IceCube detector livetime is $3.59\pm0.60$ ($2.73\pm0.38$) $ν_μ+\barν_μ$ with (without) major flares included in our analysis. This estimate exceeds, within the uncertainties, the $95\%$ ($90\%$) threshold value for the detection of one or more muon (anti-)neutrinos. Meanwhile, the most conservative scenario, where no correlation of $γ$-rays and neutrinos is assumed, predicts $1.60\pm0.16$ $ν_μ+\barν_μ$ events. We conclude that a non-detection of high-energy neutrinos by IceCube would probe the neutrino/$γ$-ray flux correlation during major flares or/and the hadronic contribution to the blazar emission.

preprint2015arXiv

Constraints of Flat Spectrum Radio Quasars in the hadronic model: the case of 3C 273

We present a method of constraining the properties of the $γ$-ray emitting region in flat spectrum radio quasars (FSRQs) in the one-zone proton synchrotron model, where the $γ$-rays are produced by synchrotron radiation of relativistic protons. We show that for low enough values of the Doppler factor $δ$, the emission from the electromagnetic (EM) cascade which is initiated by the internal absorption of high-energy photons from photohadronic interactions may exceed the observed $\sim$GeV flux. We use that effect to derive an absolute lower limit of $δ$; first, an analytical one, in the asymptotic limit where the external radiation from the broad line region (BLR) is negligible, and then a numerical one in the more general case that includes BLR radiation. As its energy density in the emission region depends on $δ$ and the region's distance from the galactic center, we use the EM cascade to determine a minimum distance for each value of $δ$. We complement the EM cascade constraint with one derived from variability arguments and apply our method to the FSRQ 3C 273. We find that $δ\gtrsim 18-20$ for $B \lesssim 30$ G and $\sim$day timescale variability; the emission region is located outside the BLR, namely at $r \gtrsim 10 R_{\rm BLR} \sim 3$ pc; the model requires at pc-scale distances stronger magnetic fields than those inferred from core shift observations; while the jet power exceeds by at least one order of magnitude the accretion power. In short, our results disfavour the proton synchrotron model for the FSRQ 3C 273.

preprint2015arXiv

Photohadronic origin of $γ$-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.

preprint2015arXiv

Relativistic Jets Shine through Shocks or Magnetic Reconnection?

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative efficiency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are efficient at dissipating energy, they typically do not accelerate particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger than the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating particles, which naturally accounts for a commonly observed property of blazar jets.

preprint2015arXiv

Spectral signatures of compact sources in the inverse Compton catastrophe limit

The inverse Compton catastrophe is defined as a dramatic rise in the luminosity of inverse Compton scattered photons. It is described by a non-linear loop of radiative processes that sets in for high values of the electron compactness and is responsible for the efficient transfer of energy from electrons to photons, predominantly through inverse Compton scatterings. We search for the conditions that drive a magnetized non-thermal source to the inverse Compton catastrophe regime and study its multi-wavelength (MW) photon spectrum. We develop a generic analytical framework and use numerical calculations as a backup to the analytical predictions. We find that the escaping radiation from a source in the Compton catastrophe regime bears some unique features. The MW photon spectrum is a broken power law with a break at $\sim m_e c^2$ due to the onset of the Klein-Nishina suppression. The spectral index below the break energy depends on the electron and magnetic compactnesses logarithmically, while it is independent of the electron power-law index ($s$). The maximum radiating power emerges typically in the $γ$-ray regime, at energies $\sim m_e c^2$ ($\sim γ_{\max} m_e c^2$ ) for $s>2$ ($s\lesssim 2$), where $γ_{\max}$ is the maximum Lorentz factor of the injected electron distribution. We apply the principles of the inverse Compton catastrophe to blazars and $γ$-ray bursts using the analytical framework we developed, and show how these can be used to impose robust constraints on the source parameters.

preprint2015arXiv

The X-ray dust scattered rings of the black hole low mass binary V404 Cyg

We report on the first detection of X-ray dust scattered rings from the Galactic low mass X-ray binary V404 Cyg. The observation of the system with Swift/XRT on June 30 2015 revealed the presence of five concentric ring-like structures centred at the position of V404 Cyg. Follow-up Swift/XRT observations allowed a time-dependent study of the X-ray rings. Assuming that these are the result of small-angle, single X-ray scattering by dust grains along the line of sight, we find that their angular size scales as $θ\propto\sqrt{t}$ in agreement with theoretical predictions. The dust grains are concentrated in five dust layers located at about 2.12, 2.05, 1.63, 1.50 and 1.18 kpc from the observer. These coincide roughly with locations of enhanced extinction as determined by infrared photometry. Assuming that the grain size distribution is described by a generalized Mathis-Rumpl-Nordsieck model, we find that the power-law index of the most distant cloud is $q\sim 4.4$, while $q \sim 3.5-3.7$ in all other clouds. We constrain at a $3σ$ level the maximum grain size of the intermediate dust layers in the range $0.16-0.20\,μ$m and set a lower limit of $\sim 0.2\,μ$m in the other clouds. Hints of an exponential cutoff at the angular intensity profile of the outermost X-ray ring suggest that the smallest grains have sizes $0.01 μ{\rm m}\le α_{\min} \lesssim 0.03\,μ$m. Based on the relative ratios of dust column densities we find the highest dust concentration at $\sim 1.6$ kpc. Our results indicate a gradient in the dust properties within 1 kpc from V404 Cyg.

preprint2014arXiv

Hadronic supercriticality as a trigger for GRB emission

We explore a one-zone hadronic model that may be able to reproduce $γ$-ray burst (GRB) prompt emission with a minimum of free parameters. Assuming only that GRBs are efficient high-energy proton accelerators and without the presence of an ab initio photon field, we investigate the conditions under which the system becomes supercritical, i.e. there is a fast, non-linear transfer of energy from protons to secondary particles initiated by the spontaneous quenching of proton-produced $γ$-rays. We first show analytically that the transition to supercriticality occurs whenever the proton injection compactness exceeds a critical value, which favours high proton injection luminosities and a wide range of bulk Lorentz factors. The properties of supercriticality are then studied with a time-dependent numerical code that solves concurrently the coupled equations of proton, photon, electron, neutron and neutrino distributions. For conditions that drive the system deep into the supercriticality we find that the photon spectra obtain a Band-like shape due to Comptonization by cooled pairs and that the energy transfer efficiency from protons to $γ$-rays and neutrinos is high reaching $\sim 0.3$. Although some questions concerning its full adaptability to the GRB prompt emission remain open, supercriticality is found to be a promising process in that regard.

preprint2014arXiv

Implications of a PeV neutrino spectral cutoff in GRB models

The recent discovery of extragalactic PeV neutrinos opens a new window to the exploration of cosmic-ray accelerators. The observed PeV neutrino flux is close to the Waxman-Bahcall upper bound implying that gamma-ray bursts (GRBs) may be the source of ultra-high energy cosmic rays (UHECRs). Starting with the assumption of the GRB-UHECR connection, we show using both analytical estimates and numerical simulations that the observed neutrinos can originate at the jet as a result of photopion interactions with the following implications: the neutrino spectra are predicted to have a cutoff at energy $\le 10$ PeV; the dissipation responsible for the GRB emission and cosmic-ray acceleration takes place at distances $r_{\rm diss} \approx 3 \times 10^{11}-3 \times 10^{13} \ {\rm cm}$ from the central engine; the Thomson optical depth at the dissipation region is $τ_{\rm T} \sim 1$; the jet carries a substantial fraction of its energy in the form of Poynting flux at the dissipation region, and has a Lorentz factor $Γ\simeq 100-500$. The non-detection of $\sim$PeV neutrinos coincident with GRBs will indicate that GRBs are either poor cosmic accelerators or the dissipation takes place at small optical depths in the jet.

preprint2014arXiv

The role of hadronic cascades in GRB models of efficient neutrino production

We investigate the effects of hadronic cascades on the gamma-ray burst (GRB) prompt emission spectra in scenarios of efficient neutrino production. By assuming a fiducial GRB spectrum and a power-law proton distribution extending to ultra-high energies, we calculate the proton cooling rate and the neutrino emission produced through photopion processes. For this, we employ a numerical code that follows the formation of the hadronic cascade by taking into account non-linear feedback effects, such as the evolution of the target photon field itself due to the contribution of secondary particles. We show that in cases of efficient proton cooling and subsequently efficient high-energy neutrino production, the emission from the hadronic cascade distorts and may even dominate the GRB spectrum. Taking this into account, we constrain the allowable values of the ratio $η_p=L_p/L_γ$, where $L_{p}$ and $L_γ$ are the isotropic equivalent proton and prompt gamma-ray luminosities. For the highest value of $η_{p}$ that does not lead to the dominance of the cascading emission, we then calculate the maximum neutrino luminosity from a single burst and show that it ranges between $(0.01- 0.6)L_{p}$ and $(0.5-1.4)L_γ$ for various parameter sets. We discuss possible implications of other parameters, such as the magnetic field strength and the shape of the initial gamma-ray spectrum, on our results. Finally, we compare the upper limit on $η_p$ derived here with various studies in the field, and we point out the necessity of a self-consistent treatment of the hadronic emission in order to avoid erroneously high neutrino fluxes from GRB models.

preprint2014arXiv

Time-dependent modelling of PKS 2155-304 in a low state

We apply both leptonic and leptohadronic emission scenarios for modelling the multiwavelength photon spectra and the observed variability in the optical, X-ray, and TeV gamma-ray energy bands of blazar PKS 2155-304 while being in a low state between 25 August and 6 September 2008. We consider three emission models, i.e. a one-component synchrotron self-Compton model (1-SSC), a one-zone proton synchrotron model (LHs), and a two-component SSC model (2-SSC). Only in the first scenario can the emission from the optical up to TeV gamma-rays be attributed to a single particle population from one emission region. Using a time-dependent numerical code that solves the kinetic equations for each particle species, we derived, in all cases, acceptable fits to the time-averaged SED. By imposing variations to one (or more) model parameters according to observed variability pattern in one (or more) frequencies we calculated the respective light curves and compared them with the observations. We show that the 1-SSC model cannot account for the anticorrelation observed between the X-rays and VHE gamma-rays, although it can explain the time-averaged SED. The anticorrelation can be more naturally explained by the two-component emission models. Both of them reproduce satisfactorily the optical, X-ray, and TeV variability but at the cost of additional free parameters, which from four in the 2-SSC model increase to six in the LHs model. Although the results of our time-resolved analysis do not favour one of the aforementioned models, they suggest that a two-component scenario is more adequate for the emission of PKS 2155-304 in the low state of 2008, which agrees with a recent independent analysis. This suggests that the quiescent blazar radiation might result from a superposition of the radiation from different components, while a flare might still be the result of a single component.

preprint2013arXiv

Mrk 421 as a case study for TeV and X-ray variability in leptohadronic models

We investigate the origin of high-energy emission in blazars within the context of the leptohadronic one-zone model. We find that $γ$-ray emission can be attributed to synchrotron radiation either from protons or from secondary leptons produced via photohadronic processes. These possibilities imply differences not only in the spectral energy distribution (SED) but also in the variability signatures, especially in the X- and $γ$-ray regime. Thus, the temporal behavior of each leptohadronic scenario can be used to probe the particle population responsible for the high-energy emission as it can give extra information not available by spectral fits. In the present work we apply these ideas to the non-thermal emission of Mrk 421. We focus on the observations of March 2001, since during that period Mrk 421 showed multiple flares that have been observed in detail both in X-rays and $γ$-rays. First, we obtain pre-flaring fits to the SED using the different types of leptohadronic scenarios. Then, we introduce random-walk type, small-amplitude variations on the injection compactness or on the maximum energy of radiating particles and follow the subsequent response of the radiated photon spectrum. For each leptohadronic scenario, we calculate the X-ray and $γ$-ray fluxes and investigate their possible correlation. Whenever the `input' variations lead, apart from flux variability, also to spectral variability, we present the resulting relations between the spectral index and the flux, both in X-rays and $γ$-rays. We find that proton synchrotron models are favoured energetically but require fine tuning between electron and proton parameters to reproduce the observed quadratic behaviour between X-rays and TeV $γ-$rays. On the other hand, models based on pion-decay can reproduce this behaviour in a much more natural way.

preprint2013arXiv

One-zone SSC model for the core emission of Centaurus A revisited

Aims: We investigate the role of the second synchrotron self-Compton (SSC) photon generation to the multiwavelength emission from the compact regions of sources that are characterized as misaligned blazars. For this, we focus on the nearest high-energy emitting radio galaxy Centaurus A and we revisit the one-zone SSC model for its core emission. Methods: We have calculated analytically the peak luminosities of the first and second SSC components by, first, deriving the steady-state electron distribution in the presence of synchrotron and SSC cooling and, then, by using appropriate expressions for the positions of the spectral peaks. We have also tested our analytical results against those derived from a numerical code where the full emissivities and cross-sections were used. Results: We show that the one-zone SSC model cannot account for the core emission of Centaurus A above a few GeV, where the peak of the second SSC component appears. We, thus, propose an alternative explanation for the origin of the high energy ($\gtrsim 0.4$ GeV) and TeV emission, where these are attributed to the radiation emitted by a relativistic proton component through photohadronic interactions with the photons produced by the primary leptonic component. We show that the required proton luminosities are not extremely high, e.g. $\sim 10^{43}$ erg/s, provided that the injection spectra are modelled by a power-law with a high value of the lower energy cutoff. Finally, we find that the contribution of the core emitting region of Cen A to the observed neutrino and ultra-high energy cosmic-ray fluxes is negligible.

preprint2013arXiv

Self-consistent neutrino and UHE cosmic ray spectra from Mrk 421

We examine the neutrino and cosmic ray spectra resulting from two models of fitting the spectral energy distribution (SED) of the blazar Mrk~421 using a self-consistent leptohadronic code. The $γ$-ray emission is attributed to either synchrotron radiation of ultra-high energy protons (LHs model) or to synchrotron radiation from electrons that result from photopion interactions of lower energy protons (LH$π$ model). Although both models succeed in fitting satisfactorily the SED, the parameter values that they use result in significantly different neutrino and cosmic-ray spectra. For the LH$π$ model, which requires high proton energy density, we find that the neutrino spectrum peaks at an energy $E_{\rm ν,peak}= 3.3$ PeV which falls well within the energy range of recent neutrino observations. While at the same time its peak flux is just under the sensitivity limit of IC-40 observations, it cannot produce ultra-high energy cosmic rays. In the LHs model, on the other hand, neutrinos are far from being detectable because of their low flux and peak energy at $E_{\rm ν,peak} \simeq 100$ PeV. However, the propagation of protons produced by the decay of escaping neutrons results in an ultra-high energy cosmic ray flux close to that observed by Pierre Augere, HiRes and Telescope Array at energies $E_p \simeq 30$ EeV.

preprint2013arXiv

Spontaneously quenched gamma-ray spectra from compact sources

We study a mechanism for producing intrinsic broken power-law gamma-ray spectra in compact sources. This is based on the principles of automatic photon quenching, according to which, gamma-rays are being absorbed on spontaneously produced soft photons, whenever the injected luminosity in gamma-rays lies above a certain critical value. We derive an analytical expression for the critical gamma-ray compactness in the case of power-law injection. For the case where automatic photon quenching is relevant, we calculate analytically the emergent steady-state gamma-ray spectra. We perform also numerical calculations in order to back up our analytical results. We show that a spontaneously quenched power-law gamma-ray spectrum obtains a photon index 3Γ/2, where Γ is the photon index of the power-law at injection. Thus, large spectral breaks of the gamma-ray photon spectrum, e.g. $ΔΓ\gtrsim 1$, can be obtained by this mechanism. We also discuss additional features of this mechanism that can be tested observationally. Finally, we fit the multiwavelength spectrum of a newly discovered blazar (PKS 0447-439) by using such parameters, as to explain the break in the gamma-ray spectrum by means of spontaneous photon quenching, under the assumption that its redshift lies in the range 0.1<z<0.24.

preprint2012arXiv

On proton synchrotron blazar models: the case of quasar 3C 279

In the present work we propose an innovative estimation method for the minimum Doppler factor and energy content of the gamma-ray emitting region of quasar 3C 279, using a standard proton synchrotron blazar model and the principles of automatic photon quenching. The latter becomes relevant for high enough magnetic fields and results in spontaneous annihilation of gamma-rays. The absorbed energy is then redistributed into electron-positron pairs and soft radiation. We show that as quenching sets an upper value for the source rest-frame gamma-ray luminosity, one has, by neccessity, to resort to Doppler factors that lie above a certain value in order to explain the TeV observations. The existence of this lower limit for the Doppler factor has also implications on the energetics of the emitting region. In this aspect, the proposed method can be regarded as an extension of the widely used one for estimating the equipartition magnetic field using radio observations. In our case, the leptonic synchrotron component is replaced by the proton synchrotron emission and the radio by the VHE gamma-ray observations. We show specifically that one can model the TeV observations by using parameter values that minimize both the energy density and the jet power at the cost of high-values of the Doppler factor. On the other hand, the modelling can also be done by using the minimum possible Doppler factor; this, however, leads to a particle dominated region and high jet power for a wide range of magnetic field values. Despite the fact that we have focused on the case of 3C 279, our analysis can be of relevance to all TeV blazars favoring hadronic modelling that have, moreover, simultaneous X-ray observations.

preprint2012arXiv

Temporal signatures of leptohadronic feedback mechanisms in compact sources

The hadronic model of Active Galactic Nuclei and other compact high energy astrophysical sources assumes that ultra-relativistic protons, electron-positron pairs and photons interact via various hadronic and electromagnetic processes inside a magnetized volume, producing the multiwavelength spectra observed from these sources. A less studied property of such systems is that they can exhibit a variety of temporal behaviours due to the operation of different feedback mechanisms. We investigate the effects of one possible feedback loop, where γ-rays produced by photopion processes are being quenched whenever their compactness increases above a critical level. This causes a spontaneous creation of soft photons in the system that result in further proton cooling and more production of γ-rays, thus making the loop operate. We perform an analytical study of a simplified set of equations describing the system, in order to investigate the connection of its temporal behaviour with key physical parameters. We also perform numerical integration of the full set of kinetic equations verifying not only our analytical results but also those of previous numerical studies. We find that once the system becomes `supercritical', it can exhibit either a periodic behaviour or a damped oscillatory one leading to a steady state. We briefly point out possible implications of such a supercriticality on the parameter values used in Active Galactic Nuclei spectral modelling, through an indicative fitting of the VHE emission of blazar 3C 279.

preprint2011arXiv

Effects of a low electron distribution cutoff on multiwavelength spectra and light curves of GRB afterglows

Aims: We investigate the behavior of the frequency-centered light curves expected within the standard model of Gamma Ray Bursts allowing the maximum electron energy to be a free parameter permitted to take low values. Methods: We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and photons, obtaining the multi-wavelength spectra as a function of time. From these we construct the frequency-centered light curves giving emphasis in the X-ray and optical bands. Results: We show that in cases where the maximum electron energy takes low values, the produced X-ray light curves show a plateau as the synchrotron component gives its place to the Synhro Self-Compton one in the X-ray band.

preprint2011arXiv

Implications of automatic photon quenching on compact gamma-ray sources

Aims: We investigate photon quenching in compact non-thermal sources. This involves photon-photon annihilation and lepton synchrotron radiation in a network that can become non-linear. As a result the gamma-ray luminosity of a source cannot exceed a critical limit that depends only on the radius of the source and on the magnetic field. Methods: We perform analytic and numerical calculations that verify previous results and extend them so that the basic properties of photon quenching are investigated. Results: We apply the above to the 2006 TeV observations of quasar 3C279 and obtain the parameter space of allowed values for the radius of the emitting source, its magnetic field strength and the Doppler factor of the flow. We argue that the TeV observations favour either a modest Doppler factor and a low magnetic field or a high Doppler factor and a high magnetic field.

preprint2011arXiv

The time-dependent one-zone hadronic model: first principles

We present some results on the radiative signatures of the one zone hadronic model. For this we have solved five spatially averaged, time-dependent coupled kinetic equations which describe the evolution of relativistic protons, electrons, photons, neutrons and neutrinos in a spherical volume containing a magnetic field. Protons are injected and lose energy by synchrotron, photopair and photopion production. We model photopair and photopion using the results of relevant MC codes, like the SOPHIA code in the case of photopion, which give accurate description for the injection of secondaries which then become source functions in their respective equations. This approach allows us to calculate the expected photon and neutrino spectra simultaneously in addition to examining questions like the efficiency and the temporal behaviour of the hadronic models.