Source author record

Apostolos Mastichiadis

Apostolos Mastichiadis 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

23works
2topics
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

23 published item(s)

preprint2025arXiv

A one-parameter two-zone leptonic model for the blazar sequence

Blazars, a subclass of radio-loud active galactic nuclei with relativistic jets aligned close to our line of sight, emit highly variable non-thermal radiation across the electromagnetic spectrum. The physical origin of their emission and the blazar sequence remain open questions. We present a self-consistent two-zone leptonic model in which relativistic electrons accelerate in a compact region, losing energy via synchrotron and inverse Compton processes, and escape into a larger zone permeated by an external photon field associated with magnetohydrodynamic winds from the accretion disk. By varying only the mass accretion rate onto the central black hole, the model naturally reproduces the blazar sequence, including Compton Dominance, $γ$-ray spectral indices, and the positions of synchrotron and inverse Compton peaks, while variations in secondary parameters account for the observed spread in the data. Flat Spectrum Radio Quasars exhibit strong external Compton emission from the extended zone, whereas BL Lac objects are dominated by synchrotron and synchrotron self-Compton emission from the compact acceleration region. This framework highlights the key role of accretion rate and spatially structured emission zones in shaping blazar spectra and provides a unified interpretation of their diverse phenomenology.

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

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

The Onset of Hadronic Supercriticality in Expanding Sources

An overlooked property of hadronic models is that they can become supercritical by abruptly transforming the energy stored in the relativistic protons into radiation. Supercriticality manifests itself when the proton density exceeds a critical value. We seek to map the complete parameter space of this behaviour in cases where a source, consisting of relativistic protons, is expanding so adiabatic losses become important. We search those critical values that lead the system to the supercritical regime and are closely related to Gamma Ray Bursts. For this reason we adopt the Relativistic Blast Wave (RBW) model which is thought to describe the production of a GRB after the initial explosion.

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

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.

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

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 statistics of BAT-to-XRT flux ratio in GRB: Evidence for a characteristic value and its implications

We present the statistics of the ratio, ${\cal R}$, between the prompt and afterglow "plateau" fluxes of GRB. This we define as the ratio between the mean prompt energy flux in {\em Swift} BAT and the {\em Swift} XRT one, immediately following the steep transition between these two states and the beginning of the afterglow stage referred to as the "plateau". Like the distribution of many other GRB observables, the histogram of ${\cal R}$~ is log-normal with maximum at a value ${\cal R}_m \simeq 2,000$, FWHM of about 2 decades and with the entire distribution spanning about 5 decades in the value of ${\cal R}$. We note that the peak of the distribution is close to the proton-to-electron mass ratio $({\cal R}_m \simeq m_p/m_e = 1836)$, as proposed to be the case in an earlier publication, on the basis of a specific model of the GRB dissipation process. It therefore appears that, in addition to the values of the energy of peak luminosity ${E_{\rm pk}\sim m_{e} c^2}$, GRB present us with one more quantity with an apparent characteristic value. The fact that the values of both these quantities (\Ep~and \R) are consistent with the same specific model invoked to account for the efficient conversion of their relativistic proton energies to electrons, argues favorably for its underlying assumptions.

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.

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.

preprint2013arXiv

The supercritical pile gamma-ray burst model: The GRB afterglow steep-decline-and-plateau phase

We present a process that accounts for the steep-decline-and-plateau phase of the Swift-XRT light curves, vexing features of GRB phenomenology. This process is an integral part of the "supercritical pile" GRB model, proposed a few years ago to provide an account for the conversion of the GRB kinetic energy into radiation with a spectral peak at $E_{\rm pk} \sim m_ec^2$. We compute the evolution of the relativistic blast wave (RBW) Lorentz factor $Γ$ to show that the radiation--reaction force due to the GRB emission can produce an abrupt, small ($\sim 25%$) decrease in $Γ$ at a radius which is smaller (depending on conditions) than the deceleration radius $R_D$. Because of this reduction, the kinematic criticality criterion of the "supercritical pile" is no longer fulfilled. Transfer of the proton energy into electrons ceases, and the GRB enters abruptly the afterglow phase at a luminosity smaller by $\sim m_p/m_e$ than that of the prompt emission. If the radius at which this slow-down occurs is significantly smaller than $R_D$, the RBW internal energy continues to drive the RBW expansion at a constant (new) $Γ$, and its X-ray luminosity remains constant until $R_D$ is reached, at which point it resumes its more conventional decay, thereby completing the "unexpected" XRT light curve phase. If this transition occurs at $R \simeq R_D$, the steep decline is followed by a flux decrease instead of a "plateau", consistent with the conventional afterglow declines. Besides providing an account of these peculiarities, the model suggests that the afterglow phase may in fact begin before the RBW reaches $R \simeq R_D$, thus introducing novel insights into the GRB phenomenology.

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.

preprint1999arXiv

Evidence for Intergalactic Absorption in the TeV Gamma-Ray Spectrum of Mkn 501

The recent HEGRA observations of the blazar Mkn 501 show strong curvature in the very high energy gamma-ray spectrum. Applying the gamma-ray opacity derived from an empirically based model of the intergalactic infrared background radiation field (IIRF), to these observations, we find that the intrinsic spectrum of this source is consistent with a power-law: dN/dE~ E^-alpha with alpha=2.00 +/- 0.03 over the range 500 GeV - 20 TeV. Within current synchrotron self-Compton scenarios, the fact that the TeV spectral energy distribution of Mkn 501 does not vary with luminosity, combined with the correlated, spectrally variable emission in X-rays, as observed by the BeppoSAX and RXTE instruments, also independently implies that the intrinsic spectrum must be close to alpha=2. Thus, the observed curvature in the spectrum is most easily understood as resulting from intergalactic absorption.