We examine a model for the observed temporal variability of powerful blazars in the $\gamma$-ray band in which the dynamics is described in terms of a stochastic differential equation, including the contribution of a deterministic drift and a stochastic term. The form of the equation is motivated by the current astrophysical framework, accepting that jets are powered through the extraction of the rotational energy of the central supermassive black hole mediated by magnetic fields supported by a so-called \emph{magnetically arrested} accretion disk. We apply the model to the $\gamma$-ray light curves of several bright blazars and we infer the parameters suitable to describe them. In particular, we examine the differential distribution of fluxes ($dN/dF_{\gamma}$) and we show that the predicted probability density function for the assumed stochastic equation naturally reproduces the observed power law shape at large fluxes $dN/dF_{\gamma} \propto F_{\gamma}^{-\alpha}$ with $\alpha>2$.

We consider the largest observed sample including all intermediate-frequency peaked (IBL) and high-frequency peaked (HBL) flaring blazars above 100 GeV up to redshift $z = 0.6$. We show that the best-fit regression line of the emitted spectral indices $\Gamma_{\rm em} (z)$ is a concave parabola decreasing as $z$ increases, thereby implying a statistical correlation between the $\{\Gamma_{\rm em} (z) \}$ distribution and $z$. This result contradicts our expectation that such a distribution should be $z$-independent. We argue that the above correlation does not arise from any selection bias. We show that our expectation naturally emerges provided that axion-like particles (ALPs) are put into the game. Moreover, ALPs can also explain why flat spectrum radio quasars emit up to 400 GeV, in sharp contradiction with conventional physics. So, the combination of the two very different but consistent results -- taken at face value -- leads to a hint at an ALP with mass $m = {\cal O} (10^{-10} \, {\rm eV})$ and two-photon coupling in the range $2.94 \times 10^{- 12} \, {\rm GeV}^{- 1} < g_{a \gamma \gamma} < 0.66 \times 10^{- 10} \, {\rm GeV}^{- 1}$. As a bonus, the Universe would become considerably more transparent above energies $E \gtrsim 1 \, {\rm TeV}$ than dictated by conventional physics. Our prediction can be checked not only by the new generation of observatories like CTA, HAWC, GAMMA-400, LHAASO, TAIGA-HiSCORE and HERD, but also thanks to the planned laboratory experiments ALPS II (upgraded), STAX, IAXO and with other techniques now being developed by Avignone and collaborators.

Within the classical emission model, where the emission region is placed within the broad line region (BLR), flat spectrum radio quasars (FSRQs) were believed not to emit photons with energies above few tens of GeV because of the absorption with the optical-UV photons from the BLR. However, photons with observed energies up to about $300 \, \rm GeV$ have been detected for few FSRQs, whose most iconic example is PKS 1441+25 at redshift $z = 0.94$. The most conservative explanation for these observations is that the emission occurs at distances comparable to the size of the dusty torus. In this case, absorption of high-energy gamma-ray photons for energies above $200-300 \, {\rm GeV}$ is dominated by the interaction with infrared radiation emitted by the torus. We investigate if current observational data about FSRQs in flaring state can give us information about: (i) the importance of the torus absorption and (ii) the properties of the torus i.e. its temperature and its geometry. We find that present data do not arrive at energies where the torus influence is prominent and as a result it is currently hardly possible to infer torus properties from observations. However, with dedicated simulations, we demonstrate that observations with the forthcoming Cherenkov Telescope Array (CTA) will be able to constrain the torus parameters (temperature and geometry).

Surprisingly, the contribution to the cosmic opacity to UHE ($E>10^{18} \, \rm eV$) $\gamma$-rays has been systematically computed so far for the $\gamma \gamma \to e^+e^-$ and $\gamma \gamma \to e^+ e^- e^+ e^-$ processes alone. We go a step further by systematically evaluating the additional opacity brought about by other leptons and hadrons. We find that the dominant channels are those leading to the production of $\mu^{\pm}$ and hadrons (mainly $\pi^{\pm}$, $\pi^0$, $K^{\pm}$, $K^0$, $\eta$). For nearly the GZK radius, the photon survival probability becomes smaller by about a factor of two with respect to current estimates.

Axion-like particles (ALPs) are light, neutral, pseudo-scalar bosons predicted by several extensions of the Standard Model of particle physics -- such as the String Theory -- and are supposed to interact primarily only with two photons. In the presence of an external magnetic field, photon-ALP oscillations occur and can produce sizable astrophysical effects in the very-high energy (VHE) band ($100 \, {\rm GeV} - 100 \, {\rm TeV}$). Photon-ALP oscillations increase the transparency of the Universe to VHE photons partially preventing the gamma-gamma absorption caused by the Extragalactic Background Light (EBL). Moreover, they have important implications for active galactic nuclei (AGN) by modifying their observed spectra both for flat spectrum radio quasars (FSRQs) and BL Lacs. Many attempts have been made in order to constrain the ALP parameter space by studying irregularities in spectra due to photon-ALP conversion in galaxy clusters and the consequences of ALP emission by stars. Upcoming new VHE photon detectors like CTA, HAWC, GAMMA-400, LHAASO, TAIGA-HiSCORE, HERD and ASTRI will settle the issue.

Very-high-energy (VHE) BL Lac spectra extending above $10 \, \rm TeV$ provide a unique opportunity for testing physics beyond the standard model of elementary particle and alternative blazar emission models. We consider the hadron beam, the photon to axion-like particle (ALP) conversion, and the Lorentz invariance violation (LIV) scenarios by analyzing their consequences and induced modifications to BL Lac spectra. In particular, we consider how different processes can provide similar spectral features (e.g. hard tails) and we discuss the ways they can be disentangled. We use HEGRA data of a high state of Markarian 501 and the HESS spectrum of the extreme BL Lac (EHBL) 1ES 0229+200. In addition, we consider two hypothetical EHBLs similar to 1ES 0229+200 located at redshifts $z=0.3$ and $z=0.5$. We observe that both the hadron beam and the photon-ALP oscillations predict a hard tail extending to energies larger than those possible in the standard scenario. Photon-ALP interaction predicts a peak in the spectra of distant BL Lacs at about $20-30 \, \rm TeV$, while LIV produces a strong peak in all BL Lac spectra around $\sim 100 \, \rm TeV$. The peculiar feature of the photon-ALP conversion model is the production of oscillations in the spectral energy distribution, so that its detection/absence can be exploited to distinguish among the considered models. The above mentioned features coming from the three models may be detected by the upcoming Cherenkov Telescope Array (CTA). Thus, future observations of BL Lac spectra could eventually shed light about new physics and alternative blazar emission models, driving fundamental research towards a specific direction.

Axion-like particles (ALPs) are attracting increasing interest since, among other things, they are a prediction of many extensions of the standard model of elementary particles physics and in particular of superstrings and superbranes. Remarkably, depending on the set of their parameter space, they strongly increase the photon transparency in the very-high energy band. The recent discovery of photon dispersion on the CMB requires a substantial modification of the previous picture: this is indeed the goal of the present paper. We compute the photon survival probability from a blazar to us exactly, and we plot it versus the observed energy for 7 simulated blazars at different $z$ and 4 values of a model parameter. Our predictions can be tested by the new generation of $\gamma$-ray observatories like CTA, HAWC, GAMMA-400, LHAASO, TAIGA-HiSCORE and HERD. Finally, for our guessed values of $m_a$ and $g_{\gamma \gamma a}$ our ALP can be detected in the upgrade of ALPS II at DESY, the planned experiments IAXO, STAX and ABRACADABRA as well as with other techniques.

Basically, in certain circumstances axion-like particles (ALPs) substantially enhance the photon survival probability $P_{\gamma \to \gamma} ({\cal E})$ of a beam emitted by a far-away source through the mechanism of photon-ALP oscillation (${\cal E}$ denotes the energy). But in order for this to work, an external magnetic field ${\bf B}$ must be present. In several cases ${\bf B}$ is modeled as a domain-like network with `sharp edges': all domains have the same size $L_{\rm dom}$ (set by the ${\bf B}$ coherence length) and the same strength B, but the direction of ${\bf B}$ changes randomly and abruptly from one domain to the next. It is obviously a highly mathematical idealization wherein the components of ${\bf B}$ are discontinuous across the edges (whence the name sharp edges). It is therefore highly desirable to go a step further, and to find out what happens when the edges are smoothed out, namely when the abrupt change of ${\bf B}$ is replaced by a smooth one. Moreover, this step becomes compelling when the photon-ALP oscillation length $l_{\rm osc}$ turns out to be comparable to -- or smaller than -- $L_{\rm dom}$, because in this case the photon survival probability $P_{\gamma \to \gamma} ({\cal E})$ critically depends on the domain shape. In the present paper we propose a smoothed out version of the previous domain-like structure of ${\bf B}$ which incorporates the above changes, and we work out its implications. Even in the present case we are able to solve analytically and exactly the photon/ALP beam propagation equation inside a single smoothed-out domain, thereby evaluating the corresponding $P_{\gamma \to \gamma} ({\cal E})$ exactly. Our results is of particular importance in view of the new generation of gamma-ray detectors, since in such a situation $l_{\rm osc} \lesssim L_{\rm dom}$ occurs just above the TeV scale.

Prompted by the increasing interest of axion-like particles (ALPs) for very-high-energy (VHE) astrophysics, we have considered a full scenario for the propagation of a VHE photon/ALP beam emitted by a BL Lac and reaching us in the light of the most up-to-date astrophysical information and for energies up to above $100 \, \rm TeV$. During its trip, the beam -- generated in a small region of a BL Lac jet -- crosses a variety of magnetic structures in very different astronomical environments: the BL Lac jet, the host elliptical galaxy, the extragalactic space and the Milky Way. We have taken an effort to model all these magnetic fields in the most realistic fashion and using a new model developed by us concerning the extragalactic magnetic field. Assuming an intrinsic spectrum with a power law exponentially truncated at a fixed cut-off energy, we have evaluated the resulting observed spectra of Markarian 501, the extreme BL Lac 1ES 0229+200 and a similar source located at $z = 0.6$ up to above $100 \, \rm TeV$. We obtain interesting results: the model with photon-ALP oscillations possesses features (spectral energy oscillatory behaviour and photon excess above $20 \, \rm TeV$) which can be tested by $\gamma$-ray observatories like CTA, HAWC, GAMMA 400, LHAASO, TAIGA-HiSCORE and HERD. In addition, our ALP can be detected in dedicated laboratory experiments like the upgrade of ALPS II at DESY, the planned IAXO and STAX experiments, as well as with other techniques developed by Avignone and collaborators.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

A strong bound on the properties of axion-like particles (ALPs) has been set by assuming that ALPs are emitted by the protoneutron star just before the core-bounce in Galactic core-collapse supernovae, and that these ALPs subsequently convert to $\gamma$-ray photons which ought to be detected by a $\gamma$-ray mission. This argument has been applied to supernova 1987A to derive the bound on the ALP-photon coupling $g_{a \gamma \gamma} \lesssim 5.3 \cdot 10^{- 12} \, {\rm GeV}^{- 1}$ for an ALP mass $m_a \lesssim 4.4 \cdot 10^{- 10} \, {\rm eV}$, and can be applied to the next Galactic supernova to derive the even stronger bound $g_{a \gamma \gamma} \lesssim 2 \cdot 10^{- 13} \, {\rm GeV}^{- 1}$ for an ALP mass $m_a \lesssim 10^{- 9} \, {\rm eV}$. We carefully analyze the considered ALP production mechanism and find that it is oversimplified to an unacceptable extent. By taking into account the minimal ingredients required by a realistic analysis, we conclude that the previous results are doomed to failure. As a consequence, all papers quoting the above bound should be properly revised. Yet, since we are unable to rule out the possibility that protoneutron stars emit ALPs, in case a core-collapse supernova explodes in the Galaxy the $\gamma$-ray satellite missions active at that time should look for photons possibly coming from the supernova.

e-ASTROGAM is a space mission with unprecedented sensitivity to photons in the MeV range, proposed within the ESA M5 call. In this note we describe some measurements sensitive to axion-like particles, for which performance in the MeV/GeV range is of primary importance, and e-ASTROGAM could be the key for discovery. Keywords: Axions. ALPs. Dark-matter candidates. Gamma-ray astrophysics.

The idea that photons can convert into axion-like particles (ALPs) $\gamma \to a$ in or around an AGN and reconvert back into photons $a \to \gamma$ in the Milky Way magnetic field has been put forward in 2008 and has recently attracted growing interest. Yet, so far nobody has estimated the conversion probability $\gamma \to a$ as carefully as allowed by present-day knowledge. Our aim is to fill this gap. We first remark that AGN that can be detected above 100 GeV are blazars, namely AGN with jets, with one of them pointing towards us. Moreover, blazars fall into two well defined classes: BL Lac objects (BL Lacs) and Flat Spectrum Radio Quasars (FSRQs), with drastically different properties. In this Letter we report a preliminary evaluation of the $\gamma \to a$ conversion probability inside these two classes of blazars. Our findings are surprising. Indeed, while in the case of BL Lacs the conversion probability turns out to be totally unpredictable due to the strong dependence on the values of the somewhat uncertain position of the emission region along the jet and strength of the magnetic field therein, for FSRQs we are able to make a clear-cut prediction. Our results are of paramount importance in view of the planned very-high-energy photon detectors like the CTA, HAWK, GAMMA-400 and HISCORE.

Flat spectrum radio quasars (FSRQs) are a particular class of blazars rich of optical/ultraviolet photons inside the broad line region (BLR), necessarily implying a huge optical depth for $\gamma$ rays above 20 GeV. As a consequence, photons with energy above such a threshold should not be emitted. However, photons in the energy range 70 - 400 GeV have been observed by MAGIC from the FSRQ PKS 1222+216. Several astrophysical explanations exist in the literature, but they are all ad hoc, namely devised only for that specific purpose. We show that such a surprising discovery can be explained within standard blazar models by adding the new possibility that photons oscillate into axion-like particles (ALPs) and vice-versa inside the source. Through the photon-ALP oscillation mechanism a sizable fraction of very-high energy (VHE) photons can escape absorption from the BLR in a similar fashion as they largely avoid absorption from the extragalactic background light (EBL) in the intergalactic space. Actually, we show that not only are VHE photon indeed emitted, but also that their spectral energy distribution (SED) is such that they lie on the same Compton peak to which also lower energy photons simultaneously detected by Fermi/LAT belong.

An intriguing possibility to partially circumvent extragalactic background light (EBL) absorption in very-high-energy (VHE) observations of blazars is that photons convert into axion-like particles (ALPs) $\gamma \to a$ inside or close to a blazar and reconvert into photons $a \to \gamma$ in the Milky Way magnetic field. This idea has been put forward in 2008 and has attracted a considerable interest. However, while the probability for the back-conversion $a \to \gamma$ has been computed in detail (using the maps of the Galatic magnetic field), regretfully no realistic estimate of the probability for the conversion $\gamma \to a$ inside a blazar has been performed, in spite of the fact that the present-day knowledge allows this task to be accomplished in a reliable fashion. We present a detailed calculation that fills this gap, considering both types of blazars, namely BL Lac objects (BL Lacs) and flat spectrum radio quasars (FSRQ) with their specific structural and environmental properties. We also include the host elliptical galaxy into account. Our somewhat surprising results show that the conversion probability in BL Lacs is strongly dependent on the source parameters -- like the position of the emission region along the jet and the strength of the magnetic field therein -- making it effectivelly unpredictable. On the other hand, the lobes at the termination of FSRQ jets lead to an effective "equipartition" between photons and ALPs due to its chaotic nature, thereby allowing us to make a clear-cut prediction. These results are quite important in view of the planned VHE detectors like the CTA, HAWK and HiSCORE.

In this paper we explore a possible route of black hole seed formation that appeal to a model by Davies, Miller & Bellovary who considered the case of the dynamical collapse of a dense cluster of stellar black holes subjected to an inflow of gas. Here, we explore this case in a broad cosmological context. The working hypotheses are that (i) nuclear star clusters form at high redshifts in pre-galactic discs hosted in dark matter halos, providing a suitable environment for the formation of stellar black holes in their cores, (ii) major central inflows of gas occur onto these clusters due to instabilities seeded in the growing discs and/or to mergers with other gas-rich halos, and that (iii) following the inflow, stellar black holes in the core avoid ejection due to the steepening to the potential well, leading to core collapse and the formation of a massive seed of $<~ 1000\, \rm M_\odot$. We simulate a cosmological box tracing the build up of the dark matter halos and there embedded baryons, and explore cluster evolution with a semi-analytical model. We show that this route is feasible, peaks at redshifts $z <~ 10$ and occurs in concomitance with the formation of seeds from other channels. The channel is competitive relative to others, and is independent of the metal content of the parent cluster. This mechanism of gas driven core collapse requires inflows with masses at least ten times larger than the mass of the parent star cluster, occurring on timescales shorter than the evaporation/ejection time of the stellar black holes from the core. In this respect, the results provide upper limit to the frequency of this process.

One of the mysteries of very-high-energy (VHE) astrophysics is the observation of flat spectrum radio quasars (FSRQs) above about 30 GeV, because at those energies their broad line region should prevent photons produced by the central engine to escape. Although a few astrophysical explanations have been put forward, they are totally ad hoc. We show that a natural explanation emerges within the conventional models of FSRQs provided that photon-ALP oscillations take place inside the source for the model parameters within an allowed range.

D. Wouters and P. Brun in Phys. Rev. D 86, 043005 (2012) claim that an observable effect in the spectra of distant very-high-energy blazars arises as a consequence of oscillations of photons into axion-like particles (ALPs) in the presence of turbulent extra-galactic magnetic fields. The main objective of this comment is to demonstrate that such a result is physically incorrect. We also show that a physically correct treatment of the same issue leads to a much less relevant conclusion, which makes the effect pointed out by WB likely unobservable with the present capabilities.

Using the most recent observational data concerning the Extragalactic Background Light and the Radio Background, for a source at a redshift z_s < 3 we compute the energy E_0 of an observed gamma-ray photon in the range 10 GeV < E_0 < 10^13 GeV such that the resulting optical depth tau_gamma(E_0,z_s) takes the values 1, 2, 3 and 4.6, corresponding to an observed flux dimming of e^-1 = 0.37, e^-2 = 0.14, e^-3 = 0.05 and e^-4.6 = 0.01, respectively. Below a source distance D = 8 kpc we find that tau_gamma(E_0,DH_0/c) < 1 for any value of E_0. In the limiting case of a local Universe (z_s = 0) we compare our result with the one derived in 1997 by Coppi and Aharonian. The present achievement is of paramount relevance for the planned ground-based detectors like CTA, HAWC and HiSCORE.

Many extensions of the Standard Model predict the existence of ALPs, which are very light spin-zero bosons with a two-photon coupling. Photon-ALP oscillations occur in the presence of an external magnetic field, and ALPs can lead to observable effects on the measured photon spectrum of astrophysical sources. An intriguing situation arises when blazars are observed with the Cherenkov Telescopes H.E.S.S., MAGIC, CANGAROO III and VERITAS. The extragalactic background light (EBL) produced by galaxies during cosmic evolution gives rise to a source dimming which becomes important in the VHE band. This dimming can be considerably reduced by photon-ALP oscillations in the large-scale magnetic fields, and the resulting blazar spectra become harder than expected. We find that for ALPs lighter than 5 x 10^{-10} eV the photon survival probability is larger than predicted by conventional physics above a few hundred GeV. This is a clear-cut prediction which can be tested with the planned Cherenkov Telescope Array and HAWC. Moreover, we offer a new interpretation of the VHE blazars detected so far, according to which the large spread in the values of the observed spectral index is mainly due to the wide spread in the source distances rather than to large variations of their internal physical properties.

The surprising discovery by MAGIC of an intense, rapidly varying emission in the energy range 70 - 400 GeV from the flat spectrum radio quasar PKS 1222+216 represents a challenge for all interpretative scenarios. Indeed, in order to avoid absorption of \gamma rays in the dense ultraviolet radiation field of the broad line region (BLR), one is forced to invoke some unconventional astrophysical picture, like for instance the existence of a very compact (r\sim 10^{14} cm) emitting blob at a large distance (R \sim10^{18} cm) from the jet base. We offer the investigation of a scenario based on the standard blazar model for PKS 1222+216 where \gamma rays are produced close to the central engine, but we add the new assumption that inside the source photons can oscillate into axion-like particles (ALPs), which are a generic prediction of several extensions of the Standard Model of elementary particle interactions. As a result, a considerable fraction of very-high-energy photons can escape absorption from the BLR through the mechanism of photon-ALP oscillations much in the same way as they largely avoid absorption from extragalactic background light when propagating over cosmic distances in the presence of large-scale magnetic fields in the nG range. In addition we show that the above MAGIC observations and the simultaneous Fermi/LAT observations in the energy range 0.3 - 3 GeV can both be explained by a standard spectral energy distribution for experimentally allowed values of the model parameters. In particular, we need a very light ALP just like in the case of photon-ALP oscillations in cosmic space. Moreover, we find it quite tantalizing that the most favorable value of the photon-ALP coupling happens to be the same in both situations. Although our ALPs cannot contribute to the cold dark matter, they are a viable candidate for the quintessential dark energy. [abridged]

Several extensions of the Standard Model predict the existence of Axion-Like Particles (ALPs), very light spin-zero bosons with a two-photon coupling. ALPs can give rise to observable effects in very-high-energy astrophysics. Above roughly 100 GeV the horizon of the observable Universe progressively shrinks as the energy increases, due to scattering of beam photons off background photons in the optical and infrared bands, which produces e+e- pairs. In the presence of large-scale magnetic fields photons emitted by a blazar can oscillate into ALPs on the way to us and back into photons before reaching the Earth. Since ALPs do not interact with background photons, the effective mean free path of beam photons increases, enhancing the photon survival probability. While the absorption probability increases with energy, photon-ALP oscillations are energy-independent, and so the survival probability increases with energy compared to standard expectations. We have performed a systematic analysis of this effect, interpreting the present data on very-high-energy photons from blazars. Our predictions can be tested with presently operating Cherenkov Telescopes like H.E.S.S., MAGIC, VERITAS and CANGAROO III as well as with detectors like ARGO-YBJ and MILAGRO and with the planned Cherenkov Telescope Array and the HAWC-ray observatory. ALPs with the right properties to produce the above effects can possibly be discovered by the GammeV experiment at FERMILAB and surely by the planned photon regeneration experiment ALPS at DESY.

Observations with the Imaging Atmospheric Cherenkov Telescopes H.E.S.S., MAGIC, CANGAROO III and VERITAS have shown that the Universe is more transparent than expected to gamma rays above 100GeV. As a natural explanation, the DARMA scenario has previously been proposed, wherein photons can oscillate into a new very light axionlike particle and vice-versa in the presence of cosmic magnetic fields. Here we demonstrate that the most recent observations further support the DARMA scenario, thereby making the existence of a very light axion-like particle more likely.