Free Access
Issue
A&A
Volume 607, November 2017
Article Number A18
Number of page(s) 12
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201731452
Published online 30 October 2017

© ESO, 2017

1. Introduction

The extragalactic gamma-ray background (EGB) represents a superposition of all γ-ray sources, both individual and diffuse, from the Milky Way to the edge of the observable universe, and provides a view of the high-energy processes in the universe. Here we consider the total γ-ray photon flux produced outside of the Milky Way, including both resolved and unresolved sources. Indeed, the diffuse Galactic emission produced by the interaction of Galactic cosmic rays (CR), mainly protons and electrons, with the Galactic interstellar medium (ISM) and interstellar radiation field, is comparable to the EGB intensity and represents a strong foreground to the EGB measurement. The latter have recently been measured by the Large Area Telescope (LAT, Atwood et al. 2009 ) on board the Fermi Gamma-ray Space Telescope (Fermi), in the range between 100 Mev and 820 GeV (Ackermann et al. 2015). The EGB spectrum is well described by a power-law with exponential cut-off having a spectral index of ~2.3 and cut-off energy greater than 300 GeV.

How much different source classes contribute to the EGB remains one of the main unanswered questions of γ-ray astrophysics. Well-established astrophysical populations, whose brightest members have been robustly detected, represent guaranteed components to the EGB. Among these, the extragalactic components are blazars, radio galaxies, and star-forming galaxies (see Fornasa & Sánchez-Conde 2015, for a review).

Blazars are among the brightest γ-ray emitters in the sky. They account for ~50% of the EGB in the energy interval Eγ ≲ 10 GeV, and for ~85% at Eγ ≳ 10 GeV (Ajello et al. 2015). They are interpreted as active galactic nuclei (AGN) with the relativistic jet directed towards the observer. The γ-ray emission in blazars is produced by inverse Compton (IC) scattering of the electrons accelerated in the jet and either the synchrotron photons emitted by the same leptonic population (synchrotron-self Compton), or from accretion disk photons (external Compton).

According to the AGN unification model (Antonucci & Miller 1985) the viewing angle discriminates among blazars and radio galaxies. With no doppler boost, radio galaxies are expected to be less bright but more abundant than blazars (blazars represent ~10% of the AGN population), making them potentially important contributors to the EGB. However their contribution to EGB is not well constrained ranging from ~7% to ~30% of the EGB intensity measured by Fermi-LAT at Eγ ≲ 10 GeV (Inoue 2011; Di Mauro et al. 2014; Wang & Loeb 2016a).

Recently, the AGN population that do not exhibit relativistic jets have also been considered as possible astrophysical contributors to the EGB (Tamborra et al. 2014; Wang & Loeb 2016a). In fact, several pieces of observational evidence indicate that AGN produce wide-angle winds with velocities of v ~ 0.1–0.3 c (e.g. Chartas et al. 2002; Pounds et al. 2003; Reeves et al. 2003; Tombesi et al. 2010, 2015). The shocks produced by the interaction of AGN winds with the ambient medium are expected to accelerate particles to relativistic energies, and the interactions of shock-accelerated particles with surrounding ISM and interstellar radiation field can produce non-thermal emission in the γ-ray band (Nims et al. 2015; Wang & Loeb 2016a; Lamastra et al. 2016).

The same emission mechanisms are expected to produce γ-rays in star-forming galaxies. In this case, the shocks are produced by supernovae (SN) explosions following star formation. Two modes of star formation have been observationally identified: a quiescent mode where the star formation is extended over the whole galactic disk and occurs on time scales of (1–2) Gyr; and a starburst mode where the star formation is concentrated in the dense, nuclear region of galaxies, and it is sustained at an enhanced rate in comparison to quiescently star-forming galaxies. There are several studies that derive the contribution to the EGB from star-forming galaxies (e.g. Fields et al. 2010; Stecker & Venters 2011; Makiya et al. 2011; Ackermann et al. 2012; Chakraborty & Fields 2013; Lacki et al. 2014; Tamborra et al. 2014). The studies that analyse quiescent and starburst galaxies separately find that the starburst contribution is always comparatively minor and the total (quiescent plus starburst) γ-ray emission is between 10% and 50% of the EGB intensity.

However, large uncertainties remain for the contribution to the EGB of the above source classes. In this paper we present improved modelling of the integrated γ-ray emission from AGN winds and star-forming galaxies by incorporating a physical model for the γ-ray emission produced by particles accelerated in AGN-driven and SN-driven shocks (Lamastra et al. 2016) into a state-of-the-art semi-analytic model (SAM) of galaxy formation (Menci et al. 2014). Our SAM includes a physical description of starburst and AGN activities triggered by galaxy interactions during their merging histories, and is ideally suited for this goal as it has been tested against several observational properties of the AGN and galaxy population both in the local and high-redshift universe, and in different electromagnetic bands (e.g. Menci et al. 2005, 2006; Lamastra et al. 2010, 2013a,b; Menci et al. 2014; Gatti et al. 2015). Moreover, galaxy and AGN number densities, and galaxy properties that determine the γ-ray emission, like gas mass, star formation rate (SFR), and AGN bolometric luminosity, can be calculated self-consistently by our SAM. This represents an advantage of the semi-analytic approach with respect to previous studies based on parametric expressions for the evolution of the AGN and galaxy populations (derived from observations in a particular electromagnetic band), and on simple scaling laws to relate the γ-ray luminosity with the properties of the host galaxies.

The paper is organised as follows. Section 2 describes the physical processes producing γ-ray emission in AGN-driven and SN-driven shocks. A description of the SAM is given in Sect. 3. Section 4 describes how we model the γ-ray emission from AGN winds and star-forming galaxies. In Sect. 5 we derive the contribution to the EGB from AGN winds and star-forming galaxies, and the Discussion and Conclusions follow in Sects. 6 and 7.

2. Gamma-ray emission from astrophysical shocks

The standard paradigm for the origin of the γ-ray emission in star-forming galaxies is non-thermal emission from relativistic particles accelerated in the shocks produced by SN explosions. Similarly to the shocks surrounding SN remnants, the shocks produced by the interaction of AGN winds with the surrounding ISM are expected to accelerate particles to relativistic energies (Nims et al. 2015; Wang & Loeb 2016a; Lamastra et al. 2016). In fact, outflows of ionised, neutral, and molecular gas, extended from a few milli-pc to kpc scales from the central supermassive black hole (SMBH), are now commonly observed in local and high-redshift AGN (see Fiore et al. 2017, and references therein). The most powerful of these AGN winds are made by fast (v ~ 0.1 − 0.3 c) highly ionised gas particles that are likely accretion disc particles accelerated by the AGN radiation field. The shock pattern resulting from the impact of an AGN wind on the ISM gas is similar to that of the stellar wind hitting the ISM around it (e.g. Weaver et al. 1977; King & Pounds 2015; King 2003; King & Pounds 2003; King et al. 2011; Lapi et al. 2005; Faucher-Giguère & Quataert 2012; Zubovas & King 2012; Zubovas & Nayakshin 2014). The wind-ISM interaction is expected to drive an outer forward shock into the ISM accelerating the swept-up material, and an inner reverse shock into the wind decelerating itself, separated by a contact discontinuity. The cooling properties of the shocked wind gas determines whether the outflow is energy- or momentum-driven. In the limit of efficient cooling of the shocked wind gas, most of the pre-shock kinetic energy is radiated away, and only its momentum flux is transferred to the ISM (momentum-driven). In contrast, if the shocked wind gas does not cool, all the energy initially provided by the shock is retained within the systems, and the shocked wind gas expands adiabatically, pushing the ISM gas away (energy-driven).

Inelastic collisions between CR protons accelerated by AGN-driven and SN-driven shocks with ambient protons may produce a significant γ-ray emission. In fact, inelastic proton-proton collisions produce neutral and charged pions. Neutral pions decay into two γ-rays: π0γ + γ; while charged pions decay into secondary electrons and positrons and neutrinos: π+μ+ + νμ and ; and . CR electrons can also produce γ-ray emission either through interaction with ISM gas (bremsstrahlung) or interstellar radiation field (IC scattering).

In our previous paper (Lamastra et al. 2016) we developed a physical model for the γ-ray emission from relativistic protons and electrons accelerated by astrophysical shocks. This model was used to predict the γ-ray spectrum produced by CR particles accelerated by the shocks observed in the molecular disk of the Seyfert galaxy NGC 1068. In this paper we derive the γ-ray emission from AGN winds and star-forming galaxies in a cosmological context by including the physical model for the γ-ray emission into a semi-analytic model of hierarchical galaxy formation. Our aim is to compare the model predictions with the measurement of the EGB intensity performed by the Fermi-LAT in the range between 100 MeV and 820 GeV (Ackermann et al. 2015). In this energy range leptonic γ-ray emission is expected to be lesser than the hadronic one; thus we limit the calculation of the γ-ray spectrum to the hadronic component.

Here we briefly recall the basic points of our model. We assume that protons are accelerated by diffusive shock acceleration (DSA) to relativistic energies in the forward outflow shock. The resulting proton number density per unit volume can be expressed as a power-law with spectral index p ≃ 2 and an exponential high-energy cut-off (Bell 1978a,b; Blandford & Ostriker 1978; Drury 1983): (1)The normalisation constant Ap is determined by the total energy supplied to relativistic protons at the shock, and Emax is the maximum energy of accelerated protons. The latter can be obtained by equating the proton acceleration time , where e is the electron charge, vs is the shock velocity, and B is the magnetic field strength, to either the time scale of proton-proton collisions τpp ≈ 5 × 107 yr / (nH/ cm-3), where nH in the ISM number density, or the outflow time scale τs = Rs/vs. Thus: (2)where vs,8 is the shock velocity in units of 108 cm/s, τage,3 is the age of the accelerator in units of 103 yr, and BμG is the magnetic field strength in units of μG (Reynolds 2008).

We constrain the normalisation constant Ap in Eq. (1) as: (3)where Emin = mpc2 is the minimum energy of an accelerated proton which is set to be proton rest mass, Ekin is the kinetic energy of the shocked particles, and ηp is the fraction of the kinetic energy transferred to protons. For the latter, we adopt ηp ≃ 0.1, that is, the value assumed in standard SN-driven shocks (Keshet et al. 2003; Thompson et al. 2006; Tatischeff 2008; Lacki et al. 2010).

We compute the γ-ray spectrum produced by neutral pion decay using the δ-functional approximation (Stecker 1970; Aharonian & Atoyan 2000): (4)where V is the volume of the outflow, , and Eπ and mπ are the energy and mass of the neutral pion. The emissivity of π0 is given by: (5)where x = mpc2 + Eπ/kpp, kpp = 0.17 is the fraction of the accelerated proton energy that goes to neutral pions in each interactions, σpp is the inelastic cross-section of a proton-proton collision, and N(x) is the accelerated proton energy distribution.

thumbnail Fig. 1

Schematic representation of the interaction of CR protons accelerated in the AGN blast wave with ISM protons. The shock radius Rs(t) expands outwards compressing the swept gas into a thin shell, and leaving a cavity inside. The γ-ray emission from neutral pion decays occurs in the regions outside the yellow cones where the lines of sight intercept the galactic gas that has not been swept out by the blast wave.

3. The semi-analytic model

In order to connect the above modelling of source emission to the statistical description of galaxy and AGN populations in a cosmological framework, we use the SAM described in detail in Menci et al. (2014; see also Gatti et al. 2015).

3.1. Hierarchical galaxy formation and SMBH growth

The SAM connects the cosmological evolution of dark matter halos with the processes involving their baryonic content. An accurate Monte Carlo procedure is used to generate the merging trees of dark matter halos following the Press & Schecter formalism (Bond et al. 1991; Lacey & Cole 1993), and to describe the gradual inclusion of sub-halos and their dynamical friction processes and binary interactions (major and minor merging and fly-by events).

We assume a cold dark matter power spectrum of perturbations and we adopted a Hubble constant h = 0.7 in units of 100 km s-1 Mpc-1, a dark energy density parameter ΩΛ = 0.7, a matter density parameter ΩM = 0.3, and a baryon density parameter Ωb = 0.035.

The baryonic processes taking place in each dark matter halo are computed following the standard recipes commonly adopted in SAMs. Starting from an initial amount, Mgas = MhaloΩb/ ΩM, of gas at the virial temperature in each dark matter halos, we compute the mass of cold baryons which are able to radiatively cool. The cooled gas settles into a rotationally supported disk with mass Mc, disk circular velocity vd, and disk radius rd (typically ranging from 1 to 5 kpc) computed as in Mo et al. (1998).

The cooled gas mass Mc is converted into stars through two different channels: i) quiescent star formation, gradually converting the gas into stars with a rate SFR = Mc/τ given by the Schmidt-Kennicutt law with τ = 1 Gyr; or ii) starbursts following galaxy interactions (merging and fly-bys) occurring on time scales of ~107 − 108 yr given by the duration of the interaction.

We assume that all stars with masses in the SN regime explode together, giving rise to a single bubble, and that a fraction of the total energy released by SN explosions is fed back onto the galactic gas. Thus, the effect of SN feedback is to return part of the cooled gas to the hot phase. The mass Δmh returned from the cold gas content of the disk to the hot gas phase is estimated, at each time-step from canonical energy balance arguments (Kauffmann 1996; Kauffmann & Charlot 1998) as: (6)where ESN = 1051 erg is the energy of ejecta of each SN, ϵSN = 0.01–0.5 is the efficiency for the coupling of the emitted energy with the cold ISM, is the number of SN per unit solar mass, depending on the assumed initial mass function (IMF), and vc is the circular velocity of the galactic halo. The model free parameter ϵSN = 0.1 is chosen as to match the local B-band luminosity function and the Tully-Fisher relation adopting a Salpeter IMF. Although our simple modelling of SN feedback does not include a detailed treatment of the gas kinematics, including the dynamics of superbubbles (Ferrara et al. 2000), it provides a good match to the observed correlations between the outflow velocity with the galactic circular velocity, and the SFR (Calura & Menci 2009).

The luminosity produced by the stellar population of the galaxies is computed by convolving the star-formation histories of the galaxy progenitors with a synthetic spectral energy distribution (SED, Bruzual & Charlot 2003). The dust extinction affecting the above luminosities is computed assuming the dust optical depth to be proportional to the metallicity Zcold of the cold phase (computed assuming a constant effective yield) and to the disk surface density, so that for the V band . The proportionality constant is taken as to match the bright end of the local luminosity function. To compute the extinction at other wavelengths we applied a proper extinction curve (see Menci et al. 2002, 2005).

The SAM includes the growth of SMBHs from primordial seeds. The latter are assumed to be the end-product of PopIII stars with a mass Mseed = 100M (Madau & Rees 2001), and to be initially present in all galaxy progenitors. SMBHs grow by merging with other black holes following the coalescence of the host galaxies and by accretion of cold galactic gas. The latter gives rise to the AGN activity. The gas accretion is triggered by galaxy interactions. In particular, we assume the analytical description of the gas inflows induced by galaxy interactions derived by Cavaliere & Vittorini (2000; see also Menci et al. 2006, 2008; Lamastra et al. 2013b), and that in each galaxy interaction one quarter of the destabilised gas feeds the SMBH, while the remaining fraction feeds the circumnuclear starburst (Sanders & Mirabel 1996). These gas fractions are calibrated as to yield final SMBH masses matching the observed local correlations with the properties of the host galaxies. We converted the mass accretion rate BH into AGN bolometric luminosity as: (7)where η ≃ 0.1 is the efficiency for the conversion of gravitational energy into radiation (Yu & Tremaine 2002; Marconi et al. 2004). Our SAM also includes a model of AGN feedback described in detail in Sect. 3.2.

The SAM has been tested against the statistical properties of the galaxy and AGN populations at low and high redshift and in different electromagnetic bands. In particular, the model provides galaxy luminosity functions (LF) in the K-band (a proxy for the stellar mass content) and in the UV band (a proxy for the instantaneous SFR) that are in good agreement with the observed evolution of the galaxy LF in the K-band up to z ~ 3 and with the galaxy LF in the UV band up to z ~ 6 (see Menci et al. 2014). The model is also able to reproduce the well known bimodal distribution of galaxies in the colour-magnitude diagram (see Menci et al. 2014).

The model predicts AGN LFs in the UV band that are in good agreement with the observational estimates at intermediate and high luminosities up to z ~ 6. At all redshifts, the model tends to slightly overestimate the data at faint luminosities (see Menci et al. 2014). The observational scaling relations between the galaxy and AGN physical properties (such as stellar mass, SFR, SMBH mass, and BH) are also well described by the model (Menci et al. 2005, 2006; Lamastra et al. 2010, 2013a,b; Menci et al. 2014; Gatti et al. 2015).

3.2. The blast wave model for AGN feedback

Our SAM includes a physical model for AGN feedback which is related to the impulsive luminous AGN phase. As discussed in Sect. 2, mildly relativistic winds (v ~ 0.1–0.3 c ) are injected by AGN into the surrounding ISM (Chartas et al. 2002; Pounds et al. 2003; Reeves et al. 2003; Tombesi et al. 2010, 2015). As these winds propagate into the ISM, they compress the gas into a blast wave terminated by a leading shock front, which moves outward with a lower but still supersonic velocity and sweeps out the surrounding medium. The expansion of the blast wave into the ISM is described by hydrodynamical equations. Taking into account the effect of dark-matter gravity, upstream pressure, and initial density gradient, and assuming the Rankine-Hugoniot boundary condition at the shock, Lapi et al. (2005) derived an analytic expression for the radius Rs of the blast wave in the case of shock expansion in a gas with a power-law density profile ρrω, where the exponent ω is in the range 2 ≤ ω< 2.5 (see also Chevalier 1976, 1982; Weaver et al. 1977; Ostriker & McKee 1988; Franco et al. 1991).

In Menci et al. (2008) the expression for the shock radius is given in terms of the galactic disk radius, disk velocity, and Mach number ℳ = vs/cs(Rs(t)): (8)The Mach number is related to ratio between the energy ΔE injected by the AGN into the surrounding medium and the total thermal energy EMc of the ISM: (9)Thus the production of weak (ℳ ≃ 1) or strong shocks (ℳ ≫ 1) depends on the value of ΔE which is computed as: (10)here ϵAGN is the fraction of the AGN bolometric luminosity transferred to the gas in the form of kinetic energy, and τAGN is the duration of the AGN phase.

The blast wave model for AGN feedback was used in our previous papers to explain the distribution of hydrogen column densities in AGN as a function of luminosity and redshift, and to predict hydrogen phoionization rate as a function of redshift (Menci et al. 2008; Giallongo et al. 2012).

4. Model set up

In this section we describe the model parameters that we will use in the computation of the γ-ray emission from AGN winds and star-forming galaxies.

In particular, we define the parameters that describe the γ-ray spectra of individual AGN winds and star-forming galaxies, and the environment into which the shocks expand. We limit the shock expansion into galactic disks, for which we assume a constant scale height hd = 100 pc (Narayan & Jog 2002; van der Kruit & Freeman 2011), and an isothermal gas density profile nH = nH,0/r-2 (see Fig. 1). The constant nH,0 in the density profile can be constrained by the total gas content in the disk Mc.

4.1. Modelling γ-ray emission from AGN winds

To derive the γ-ray spectrum of an individual AGN wind, we need to determine the energy distribution of the particles accelerated in the shocks. The latter is determined by the CR particle spectral index p, the particle maximum energy Emax, and by the total energy supplied to relativistic particles at the shock. As discussed in Sect. 2, we assume that DSA is the mechanism which produces CR protons in AGN-driven shocks (but see Vazza et al. 2015, 2016; van Weeren et al. 2016, for results showing that DSA has difficulty in explaining the observed emissions of particles accelerated in some astrophysical shocks). DSA could result in the production of a power-law accelerated proton population with a power-law index p ≃ 2 (Bell 1978a,b; Blandford & Ostriker 1978; Drury 1983). As the particles diffuse from the acceleration region, energy-dependent diffusion losses can soften the source spectrum leading to larger values of the spectral index p. The accelerated particle energy distribution extends to energies as high as is permitted by various loss processes. Protons are accelerated up to a maximum energy that depends on the shock velocity, age of the accelerator, and on the magnetic field strength in the shock region (Eq. (2)). The first two parameters are determined by the hydrodynamics of the shocks and by the density of the galactic disk. For an isothermal gas disk, the blast wave shock radius is given by Rs = 0.9 ℳ vdt (Eq. (8)), which implies a constant outflow velocity vs = dRs/ dt = 0.9ℳvd. For the magnetic field B we assume values bracketed by a minimum value that is given by the volume average ISM magnetic field strength BISM = 6 × (Σgas/ 0.0025 g cm-2)aμG where a ≃ 0.4–1 and Σgas = 2nHmHhd is the disk gas surface density (Robishaw et al. 2008; Lacki et al. 2010; McBride et al. 2014); and a maximum value that is derived by assuming that a fraction ξB ≃ 0.1, based on observations of SN remnants (Chevalier 1998), of the post shock thermal energy is carried by the magnetic field Bshock = (8πξBnskTs)0.5 where ns and Ts are the post-shock density and temperature of the gas, respectively. To derive the maximum value of the magnetic field Bshock we assume the temperature and density jumps given by the approximations valid for very strong shocks as given in Lapi et al. (2005): ns ≃ 4nH and , irrespective of .

A constant outflow velocity corresponds to an energy-driven outflow in which the kinetic luminosity does not vary in time. In this case, the kinetic energy of accelerated protons is simply given by the product of the outflow kinetic luminosity Lkin and the residence time of the particles in the acceleration region τres:(11)In the case of AGN-driven winds, we assume that the outflow kinetic luminosity is a fraction ϵAGN of the AGN bolometric luminosity: (12)The ratio between the outflow kinetic power and AGN bolometric luminosity has recently been determined in a sample of 94 AGN by Fiore et al. (2017) to be in the range ϵAGN = 0.001–0.1. The parameter ϵAGN = 0.01 is chosen within the observational range as to provide a good fit of the bright end of the AGN luminosity function (see Menci et al. 2014).

The collisions between CR protons and ambient protons in galactic disks produce hadronic γ-ray emission. The blast wave model for AGN feedback allows us to self-consistently compute the fraction of accelerated protons that may interact with ambient protons. In fact γ-ray emission from neutral pion decays occurs along the lines of sight where the galactic gas has not been swept out by the blast wave produced by the AGN (see Fig. 1). This fraction is the complementary part of the escape fraction of ionising photons that we derived in Giallongo et al. (2012). The average of this fraction over the duration τAGN of the AGN activity is given by: (13)where τs is the time at which the shock radius first encompasses the width of the galactic disk: τs = hd/ 0.9ℳrd. For shocks with high Mach number, that is, for large AGN injected energies ΔEE, τsτAGN yielding small fractions of interacting protons. Figure 2 shows the distribution of simulated AGN as a function of and AGN bolometric luminosity for our fiducial model.

We assume that all the energy of the protons that interact with the protons in the galactic disk is converted into pion production. This corresponds to an assumtion that AGN winds act as proton calorimeters. The calorimetric limit has the maximum efficiency to convert AGN blast wave energy into γ-rays, and this corresponds to the assumption that τres = τpp in Eq. (11). The proton-proton collisional time scale τpp is inversely proportional to the density nH of target material, thus the resulting hadronic γ-ray luminosity is independent of nH, and it scales linearly with the outflow kinetic luminosity (see Eq. (4)). Figure 3 shows the predicted γ-ray spectrum of an AGN with LAGN = 7 × 1044 erg/s hosted in a halo of mass 1012M at z = 0.1. In our derivation of the γ-ray spectrum we neglect the γ-ray emission from IC and bremsstrahlung processes of the primary and secondary leptonic populations, as it is expected to be lesser than the hadronic one at Eγ ≳ 100 MeV (Lacki et al. 2014).

thumbnail Fig. 2

Distribution of simulated AGN as a function of and AGN bolometric luminosity for our fiducial model. The contours correspond to equally spaced values of the density (per Mpc3) of objects in a given LAGN bin in logarithmic scale: from 10-5 for the lightest filled region to 10-2 for the darkest.

thumbnail Fig. 3

γ-ray spectrum of an AGN with LAGN = 7 × 1044 erg/s hosted in a halo of mass 1012M at z = 0.1. Energy spectrum parameters are set at p = 2 (black lines), p = 2.2 (red lines), B = BISM (solid lines), and B = Bshock (dashed lines).

4.2. Modelling γ-ray emission from star-forming galaxies

Besides the Milky Way, ten external star-forming galaxies have been firmly detected in γ-rays with the Fermi-LAT (Ackermann et al. 2012, 2017). Among these, seven are starburst and active galaxies that are more luminous at γ-ray energies compared to quiescently star-forming galaxies by a factor larger than 10. The γ-ray spectra of starbursts look similar and can be described by a single power-law with spectral index p ~ 2.2. Because of the high density of star-forming regions in starburst galaxies, it is likely that they act as proton calorimeters (Wang & Fields 2016). Thus starburst galaxies could have harder spectra than quiescent galaxies as in their acceleration region proton energy losses are dominated by nearly energy-independent proton-proton collisions. On the contrary, in quiescent galaxies proton losses are thought to be dominated by energy-dependent diffusion (as in the Milky Way), and their γ-ray luminosity is set by equilibrium between proton injection rate, diffusion processes, and energy losses. Both kinds of star-forming galaxy are included in our SAM. However, for the above reasons, modelling the detailed escape and energy-loss process of accelerated protons in quiescent galaxies is a very difficult task. Thus, in the following we consider only the calorimetric regime, which provides a good description of starburst galaxies, and we discuss the implications for the γ-ray emission from the more numerous but less powerful quiescient galaxies.

To derive the energy distribution of particles accelerated in starburst galaxies we assume a power-law index p = 2.2 characteristic of the starbursts detected in the γ-ray band (Ackermann et al. 2012). We also assume a constant shock velocity equal to vs = 300 km s-1 (Lacki et al. 2010), and a minimum and maximum value of the magnetic field in the shock region given by BISM and Bshock, respectively, as in the case of AGN-driven shocks. The kinetic luminosity available in the form of accelerated protons for a starburst galaxy is: (14)where νSN is the SN rate, and ϵSN = 0.1 is the fraction of the SN energy transferred to the gas in the form of kinetic energy (see Sect. 3). SN rates can be estimated from SFR and the IMF as νSN = φ × SFR where φ ≃ 0.003 for a Salpeter IMF. We note that for starburst galaxies in the calorimetric regime the γ-ray luminosity scales linearly with the SFR.

Starburst galaxies are selected from the SAM as model galaxies dominated by the SFR triggered by galaxy interactions. The latter at the same time trigger AGN activity. The starburst-AGN connection predicted by the SAM (see Lamastra et al. 2013b) implies that the blast wave produced by the AGN could sweep out the disk gas when the starbursts are in action. This hampers the derivation of the fraction of protons accelerated in SN-driven shocks that interact with protons in the ISM. In fact, the latter depends on the time delay between the trigger of the AGN and starburst activities, and on the starburst spatial distribution quantity that cannot be provided by the SAM. For this reason, in the following, to derive the contribution from starburst galaxies to the EGB, we adopt an empirical calorimetric fraction, which is derived from observations of star-forming galaxies in the GeV band (Ackermann et al. 2012): (15)

5. Results

thumbnail Fig. 4

The contribution to the EGB from AGN winds (blue) and starburst galaxies (red) predicted by our SAM. Fiducial model parameters are set at: hd = 100 pc, ηp = 0.1, ϵAGN = 0.01, and ϵSN = 0.1. Accelerated proton spectral index equal to p = 2.2 is assumed for both AGN winds and starburst galaxies. The blue and red shaded bands represent the uncertainty related to the EBL model adopted (Stecker et al. 2016). The data points are the Fermi-LAT measurement of the EGB, with black and grey vertical error bars indicating EGB statistical and systematic uncertainties, respectively (Ackermann et al. 2015).

thumbnail Fig. 5

Upper panel: the contribution to the EGB from AGN winds predicted by our SAM is compared to the contribution provided by the blazar population (Ajello et al. 2015). Model predictions are shown for our fiducial model and for different accelerated proton spectral index: p = 2 (dotted blue lines), p = 2.2 (solid blue lines), and p = 2.4 (dashed blue lines). The red shaded region represent the sum of the AGN wind and blazar contributions to the EGB. The data points are as in Fig. 4. Lower panel: fraction of EGB provided by AGN winds and blazars.

In this section we derive the cumulative γ-ray emission from AGN winds and starburst galaxies predicted by the SAM. The contribution from AGN winds (starburst galaxies) to the EGB spectrum can be estimated as: (16)where is the intrinsic photon energy, φ(Lγ,z) is the comoving number density of AGN (starburst galaxies) per unit of γ-ray luminosity as a function of redshift, DL(z) is the luminosity distance, I(Eγ,Lγ,z) is the γ-ray spectrum of an individual AGN (starburst galaxy) with integral γ-ray luminosity Lγ at redshift z, the factor d2V/ dz represents the comoving volume element per unit redshift and unit solid angle, and τγγ is the diffuse extragalactic background light (EBL) optical depth for photons with energy Eγ at a redshift of z. In fact, the emitted γ-rays, while travelling through the inter galactic medium, interact with the photons of EBL and become absorbed through e+/e pair production. The absorption probability increases with energy and distance of the γ-ray source. In our calculation we assume zmax = 5 and adopt the EBL model of Stecker et al. (2016).

Figure 4 compares the cumulative γ-ray emission from AGN winds and starburst galaxies predicted by the SAM with the EGB spectrum measured by Fermi-LAT (Ackermann et al. 2015). The latter is well described by a power-law with a spectral index of ~2.3 and an exponential cut-off. Both the statistical and systematic uncertainties of the EGB measurement are shown in Fig. 4. The systematic uncertainty ranges between a factor ~15% and ~30% (depending on the energy range considered) and depends on the modelling of the Galactic diffuse emission. The model predictions are shown for our fiducial model corresponding to hd = 100 pc, ηp = 0.1, ϵAGN = 0.01, and ϵSN = 0.1. The accelerated proton spectral index is assumed to be p = 2.2 for both AGN winds and starburst galaxies. We find that our results do not depend on the exact value of the magnetic field. In fact, the magnetic field determines the cut-off energy of the EGB spectrum that is mainly affected by EBL attenuation at the bright end. Figure 4 shows that, assuming a comparable efficiency for accelerating protons in AGN-driven and SN-driven shocks (ηp = 0.1), hierarchical scenarios predict a contribution to the EGB from AGN winds larger by about a factor of 100 than that provided by starburst galaxies. This implies that the kinetic energy available to accelerate protons in AGN winds significantly exceeds that provided by star formation driven by galaxy interactions. However, in hierarchical clustering scenarios starburst galaxies account for a small fraction of the cosmic star formation rate density (SFRD), that is, the mass converted into stars per unit time and comoving volume. The starburst contribution to the SFRD was estimated to be ~5% at z ≃ 0.1 and ~20% at z ≃ 5 (Lamastra et al. 2013a). Thus, the cosmic SFRD is dominated by quiescent galaxies at all redshifts. If we assume that, at all redshifts, the γ-ray emission scales as the cosmic SFRD, we find that the contribution to the EGB from quiescent galaxies should be a factor of ~20 greater than that provided by starburst galaxies. This implies that, in hierarchical scenarios, the γ-ray emission from AGN winds dominates over that powered by star-forming galaxies (quiescent plus starburst). We note that the above order-of-magnitude estimate for the quiescent contribution to the EGB should be thought of as an upper bound. In fact, as discussed in Sect. 4.2, diffusive and advective losses cannot be neglected in quiescent galaxies, thus the hypothesis of proton calorimeter adopted for starburst galaxies, which has the maximum efficiency to convert SN blast wave into γ-rays, should not be valid in quiescent galaxies.

We find that AGN winds account for ~40% of the observed EGB in the energy interval Eγ = 0.1–1 GeV, for ~90% at Eγ = 1–10 GeV, and for ~70% at Eγ ≳ 10 GeV. Other classes of sources are known to contribute to the EGB. Among these, the major contribution is from blazars. The predicted cumulative emission of blazars is shown in Fig. 5 as a green dashed band (Ajello et al. 2015). The latter encompasses systematic uncertainties on blazar luminosity function models and energy spectrum models. We also show the integrated emission from AGN winds predicted by our SAM assuming different values of the accelerated proton spectral index. As shown in the bottom panel of Fig. 5, at energies Eγ ≲ 10 GeV, and for spectral index p> 2, the γ-ray emission from AGN winds dominates over that from blazars. The AGN wind contribution to the EGB is peaked at energies Eγ ≃ 1–10 GeV, depending on the value of p, while the blazar contribution reachs its maximum at larger energies Eγ ≃ 60 GeV. This analysis shows that AGN winds and blazars can account for the amplitude and spectral shape of the EGB, leaving only little room for other contributions.

6. Discussion

6.1. Robustness and comparison with previous works

Here we discuss our assumptions and compare our results with previous estimates of the γ-ray emission from AGN winds and star-forming galaxies.

The results shown in the previous section were obtained assuming: i) energy-driven winds powered by the active nucleus in AGN host galaxies; ii) the calorimetric regime; and iii) shock expansion limited to galactic disks.

As for the first point, the observations of winds are very common in AGN host galaxies. X-ray and UV emission and absorption line studies revealed outflows of highly ionised gas on an ~0.001–1 pc scale, with velocities v = 0.1 − 0.3 c, both at low and high redshift (see e.g. Tombesi et al. 2013, 2015). On larger scales (100–1000 pc), outflows of ionised, neutral atomic, and molecular gas, with velocities v = 102 − 103 km s-1, have been observed through deep optical/near-infrared spectroscopy and interferometric observations in the (sub)millimetre band (e.g. Rupke & Veilleux 2011; Feruglio et al. 2010, 2015; Cicone et al. 2014; Harrison et al. 2016; Shen 2016; Zakamska et al. 2016). A collection of AGN winds detected at different scales and ionisation states is given in Fiore et al. (2017). These observations indicate that the majority of the large-scale outflows are driven by the nuclear activity, and may be identified with the energy-driven phase. For two sources in particular, namely MrK 231 (Feruglio et al. 2015) and IRAS F11119 (Tombesi et al. 2015), both X-ray winds and molecular winds have been detected. The comparison between the momentum rate of the X-ray and molecular winds indicates that these winds are energy-driven. These findings seem to support our assumptions about the frequency and nature of AGN winds, however, these observations are sparse, and mostly limited to AGN selected to have high chances of being in an outflowing phase. In order to gain more insight into these topics, the measurements of the frequency and parameters of AGN winds in unbiased AGN samples over a large range of redshifts are necessary.

As for the calorimetric regime, our model assumes that advective and diffusive escape of accelerated protons in galactic disks are negligible. This condition is satisfied if the proton-proton collisional time scale τpp is less than the advective and diffusion time scales. The wind advection time is τwindhd/vs. Assuming the velocities measured in galactic scale AGN winds, and for hd = 100 pc, we expect that hadronic losses dominate advection losses when nH ≳ 50–500 cm-3, which are values observed in circumnuclear disks of active and starburst galaxies (e.g. García-Burillo et al. 2014; Yoast-Hull et al. 2014). The particle diffusion processes in the environments of active and starburst galaxies are poorly constrained. For this reason we neglect diffusive losses, and this could constitute a source of uncertainty in our computation.

thumbnail Fig. 6

EGB spectrum. The blue and red shaded regions correspond to AGN winds and starburst galaxies predicted by our SAM, respectively. The purple shaded region corresponds to the contribution to the EGB from AGN winds estimated by Wang & Loeb (2016a). The cyan and gold lines indicate the estimate for quiescent and starburst galaxies, respectively, predicted by the Makiya et al. (2011) SAM. The green shaded region shows the contribution to the EGB from blazars as estimated by Ajello et al. (2015). The data points are as in Fig. 4.

As for the environment into which the shocks expand, we limit the expansion to galactic disks, however it is also possible that shocks can propagate in gaseous haloes of galaxies. The lower densities of galactic haloes with respect to galactic disks imply a major role of escape (advective and diffusive) with respect to hadronic interactions, and therefore a low γ-ray-production rate. As a check, we can compare our predictions with other recent predictions of the EGB from AGN winds that assume shock expansion in both the disk and halo components of galaxies. Wang & Loeb (2016a) derived the cumulative γ-ray emission from AGN winds using a hydrodynamical model for AGN wind interaction with the ambient medium (Wang & Loeb 2015), and the empirical AGN bolometric luminosity function of Hopkins et al. (2007). In their computation, the γ-ray emission is assumed to be produced by the same hadronic processes considered in this work and described in Sect. 2.

Figure 6 compares the estimates for the cumulative γ-ray emission from AGN winds derived by Wang & Loeb (2016a) with that predicted by our SAM assuming the same fraction of the AGN bolometric luminosity that powers the winds, and the same energy spectrum. Figure 6 shows that the normalisation of the EGB spectrum predicted by the SAM is a factor of ~3 larger than that derived by Wang & Loeb (2016a). The SAM also predicts a different shape of the EGB spectrum at energies Eγ ≳ 1 GeV. As for the normalisation, we have checked whether the mismatch seen in Fig. 6 stems from the different AGN luminosity functions adopted. Wang & Loeb (2016a) used the AGN bolometric luminosity function of Hopkins et al. (2007) which is derived by combining measurements of AGN number density in IR, optical, and X-ray bands, in the redshift interval z = 0–6. Deriving the bolometric luminosity function from the observed luminosity functions in different electromagnetic bands is not a trivial procedure as it requires assumptions on the bolometric corrections, and corrections for obscured sources. Hopkins et al. (2007) use a luminosity-dependent bolometric correction, and they correct their data for extinction and the fraction of Compton-thick AGN missed in IR, optical, and X-ray data. They also give analytical approximations of the empirical bolometric luminosity function. In their calculation, Wang & Loeb (2016a) use the pure luminosity evolution (PLE) model given by Hopkins et al. (2007) to describe the evolution of the AGN population. As discussed by Hopkins et al. (2007), although the PLE model provides a reasonable lowest order of magnitude approximation to the data, it underpredicts the abundance of low-luminosity AGN at z ≲ 0.5. Such an underprediction of faint objects explains the normalisation difference seen in Fig. 6.

As for the shape of the EGB spectrum, we find that while the two spectral shapes are consistent at energies Eγ ≲ 1 GeV, at higher energies the SAM predicts a decline in the EGB spectrum that begins at larger energies and is shallower than that predicted by Wang & Loeb (2016a). The high-energy part of the predicted EGB spectrum is affected by γ-ray absorption in the intergalactic medium due to the EBL. There are large uncertainties regarding the EBL estimates and thus also the γ-ray optical depth used in Eq. (16) to calculate the EGB spectrum. We used τγγ given by Stecker et al. (2016), while Wang & Loeb (2016a) used the γ-ray optical depth estimated by the same authors in a previous paper (Stecker et al. 2007). The latter results in a larger absorption of γ-rays at energies Eγ ≳ 10 GeV, and this could in part explain the discrepancy in the shape of the EGB spectra seen in Fig. 6. In fact, we have verified that assuming τγγ given by Stecker et al. (2007), the SAM predicts a steeper decline of the cumulative γ-ray emission from AGN winds, however the predicted EGB spectrum still remains shallower that that predicted by Wang & Loeb (2016a).

It is also worth noting that the solutions of the hydrodynamical equations describing the outflow motion in an isothermal gas density profile give a constant outflow velocity (i.e. an energy-driven outflow) in the Lapi et al. (2005) model, and an outflow velocity that decreases with increasing outflow radius (i.e. a momentum-driven outflow) in the Wang & Loeb (2016a) model. The shock velocity determines the maximum energy of accelerated protons, which in turn shapes the high-energy part of the EGB spectrum. The maximum energy that can attain accelerated protons can be expressed in terms of the shock velocity, size of the accelerator, and magnetic field strength as EmaxvsRsB (see Eq. (2)). Since for an isothermal gas density profile, the maximum energy scales as Emaxvs. Thus in an energy-driven outflow the proton maximum energy remains constant during the outflow expansion, while in a momentum-driven outflow Emax decreases with time.

As for the contribution to the EGB from star-forming galaxies, Makiya et al. (2011) derived the contribution from both starburst and quiescent galaxies using a SAM of hierarchical galaxy formation. As in our model, the star formation in starburst galaxies is triggered by galaxy interactions, while in quiescent galaxies the star formation is determined by the cold gas reservoir and the galaxy dynamical time-scale. In the Makiya et al. (2011) model the emission of quiescent and starburst galaxies are modelled based on templates which are tuned to reproduce the γ-ray spectra of the Milky Way and M82, respectively. For quiescent galaxies, they assumed the so-called escape regime. In this regime the energy losses of accelerated protons are dominated by their escape from the diffuse region of the galaxy, and the γ-ray luminosity depends on the SFR and mass of ISM gas, that is, LγSFR × Mc. For the starburst galaxies, they assumed the calorimetric regime. In this case the γ-ray luminosity is no longer dependent on the gas mass, that is, LγSFR. The authors determined the relation between Lγ and SFR × Mc, and Lγ and SFR, by fitting the results of four star-forming galaxies detected by Fermi-LAT. The contributions to the EGB from quiescent and starburst galaxies derived by Makiya et al. (2011) are shown in Fig. 6. They found that the contribution to the EGB from quiescent galaxies is a factor of ~10 larger than that provided by starburst galaxies. The latter agrees reasonably well with our estimate. This supports the scenario discussed in Sect. 5 where the γ-ray emission from star-forming galaxies (quiescent plus starburst) predicted by hierarchical galaxy formation models is lesser than that provided by AGN winds.

6.2. Multi-messenger implications

In Sect. 5 we show that a sizeable fraction of the EGB can be accounted by hadronic γ-ray emission from AGN winds. The hadronic origin of the γ-ray flux will make the production of neutrinos unavoidable, creating a diffuse neutrino background. Indeed, in proton-proton interactions approximately two thirds of the pions produced are charged pions that decay into muons and neutrinos followed by electrons and positrons and more neutrinos: π+μ+ + νμ and ; and . Neutrinos can also be created in interactions of CR protons with the ambient radiation field. Here we focus on proton-proton interactions as they are the dominant process for AGN winds (Tamborra et al. 2014; Wang & Loeb 2016b; Lamastra et al. 2016).

thumbnail Fig. 7

Neutrino spectrum of an AGN with LAGN = 7 × 1044 erg/s hosted in a halo of mass 1012M at z = 0.1. Model parameter as in Fig. 3.

To calculate the energy spectra of neutrinos produced in AGN winds we use the parametrisations derived by Kelner et al. (2006; see Fig. 7). Then, the cumulative neutrino flux from AGN winds can be obtained by summing the neutrino emission over the entire AGN population, as described by Eq. (16). The absorption due to the EBL is ignored in the calculation of the neutrino background, and this introduces a dependence of the predicted spectrum on the values of the magnetic field in the shock region.

Figure 8 compares the cumulative neutrino background from AGN winds with the most recent IceCube data. The latter are fitted by two different models: a single power-law model, and a differential model with nine free parameters (Aartsen et al. 2015). We find that the estimated neutrino intensity is comparable, within the astrophysical uncertainties, to the IceCube measurement for spectral index p ~ 2.2–2.3. As shown in Fig. 5, these spectral indexes imply the largest contribution to the EGB from AGN winds at energies Eγ ≃ 1–10 GeV, which, when added to the contribution from the balazar population, lead to a slight overestimate of Fermi-LAT data. Of course, any plausible model for the IceCube neutrino background should not overpopulate the gamma bounds. This tension points to a class of sources that are opaque in the γ-rays (Chakraborty & Izaguirre 2016; Murase et al. 2016). As discussed in Chakraborty & Izaguirre (2016), the large photon number density present in the environment of starburst galaxies implies that γ-rays could interact inside the galaxy before escape. The AGN-starburst connection predicted by our SAM implies that the internal absorption of γ-rays could also be present in AGN host galaxies. Thus, the Fermi-LAT and IceCube data could be reproduced simultaneously by our SAM (see also Murase et al. 2013; Wang & Loeb 2016b).

The analysis of positional coincidence of IceCube neutrino events with known astrophysical sources is a difficult task owing to the poor angular resolution of the detector. For this reason, there are as yet no confirmed identifications for astrophysical sources of IceCube neutrino events (e.g. Aartsen et al. 2014; Adrián-Martínez et al. 2016b; The IceCube Collaboration et al. 2015). Recently, Padovani et al. (2016 see also Resconi et al. 2017) argued for a statistical significant correlation between IceCube neutrino events and high-energy peaked BL Lacertae (HBL) objects in the second catalogue of hard Fermi-LAT sources (Ackermann et al. 2016). Although HBL are promising neutrino-emitter candidates (Lucarelli et al. 2017), they can account for ~10–20% of the IceCube signal.

In the near future, the better angular resolution of KM3NeT (~0.2° for neutrinos with energy E ≳ 10 TeV, track-like events, Adrián-Martínez et al. 2016a) will allow us to effectively constrain the position of the possible counterparts of neutrino events, thus providing a possible direct test of neutrino background models.

thumbnail Fig. 8

Cumulative γ-ray (left) and neutrino (right) background from AGN winds predicted by our SAM. Model predictions are shown for our fiducial model and for different AGN wind spectral energy parameters: p = 2.2 (purple lines), p = 2.3 (blue lines), B = BISM (solid lines), and B = Bshock dashed lines. The circles represent Fermi-LAT data (Ackermann et al. 2015). The squares show the differential model of IceCube neutrino data, while purple dotted lines represent the power-law models (Aartsen et al. 2015).

7. Conclusions

We have incorporated the description of the hadronic γ-ray emission from relativistic protons accelerated in AGN-driven and SN-driven shocks into a state-of-the-art SAM of hierarchical galaxy formation. Our SAM has already proven to match the statistical properties of the galaxy and AGN populations at low and high redshift and in different electromagnetic bands. We have compared the predictions for the cumulative γ-ray emission from AGN winds and star-forming galaxies with the latest measurement of the EGB performed by the Fermi-LAT in the range between 100 Mev and 820 GeV (Ackermann et al. 2015). The main results of this paper are as follows:

  • In hierarchical clustering scenarios, connecting the physics ofAGN and starburst galaxies to the merging histories of the hostgalaxies, assuming a comparable efficiency for acceleratingprotons in AGN-driven and SN-driven shocks(ηp = 0.1), the contribution to the EGB from AGN winds dominates over that from starburst galaxies. If we consider also the contribution of the less powerful but more numerous quiescent galaxies, the contribution to the EGB from all star-forming galaxies is a factor ~3–5 lower than that provided by AGN winds.

  • The cumulative γ-ray emission from AGN winds and blazars can account for the amplitude and spectral shape of the EGB, assuming the standard acceleration theory, and AGN wind parameters that agree with observations. At energies lower and greater than Eγ ≃ 10 GeV the EGB is dominated by AGN winds and blazars, respectively. The transition between these two regimes could, in principle, give rise to breaks and features in the EGB energy spectrum.

  • The neutrino background resulting from charged pion decays following hadronic interactions can reproduce the IceCube data assuming accelerated proton spectral index p ~ 2.2–2.3. The Fermi-LAT data could be reproduced simultaneously, taking into account internal absorption of γ-rays.

References

  1. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2014, ApJ, 796, 109 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aartsen, M. G., Abraham, K., Ackermann, M., et al. 2015, ApJ, 809, 98 [NASA ADS] [CrossRef] [Google Scholar]
  3. Ackermann, M., Ajello, M., Allafort, A., et al. 2012, ApJ, 755, 164 [NASA ADS] [CrossRef] [Google Scholar]
  4. Ackermann, M., Ajello, M., Albert, A., et al. 2015, ApJ, 799, 86 [NASA ADS] [CrossRef] [Google Scholar]
  5. Ackermann, M., Ajello, M., Atwood, W. B., et al. 2016, ApJS, 222, 5 [NASA ADS] [CrossRef] [Google Scholar]
  6. Ackermann, M., Ajello, M., Albert, A., et al. 2017, ApJ, 836, 208 [NASA ADS] [CrossRef] [Google Scholar]
  7. Adrián-Martínez, S., Ageron, M., Aharonian, F., et al. 2016a, J. Phys. G Nucl. Phys., 43, 084001 [Google Scholar]
  8. Adrián-Martínez, S., Albert, A., André, M., et al. 2016b, ApJ, 823, 65 [Google Scholar]
  9. Aharonian, F. A., & Atoyan, A. M. 2000, A&A, 362, 937 [NASA ADS] [Google Scholar]
  10. Ajello, M., Gasparrini, D., Sánchez-Conde, M., et al. 2015, ApJ, 800, L27 [NASA ADS] [CrossRef] [Google Scholar]
  11. Antonucci, R. R. J., & Miller, J. S. 1985, ApJ, 297, 621 [NASA ADS] [CrossRef] [Google Scholar]
  12. Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bell, A. R. 1978a, MNRAS, 182, 147 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bell, A. R. 1978b, MNRAS, 182, 443 [NASA ADS] [CrossRef] [Google Scholar]
  15. Blandford, R. D., & Ostriker, J. P. 1978, ApJ, 221, L29 [Google Scholar]
  16. Bond, J. R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 440 [NASA ADS] [CrossRef] [Google Scholar]
  17. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  18. Calura, F., & Menci, N. 2009, MNRAS, 400, 1347 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cavaliere, A., & Vittorini, V. 2000, ApJ, 543, 599 [NASA ADS] [CrossRef] [Google Scholar]
  20. Chakraborty, N., & Fields, B. D. 2013, ApJ, 773, 104 [NASA ADS] [CrossRef] [Google Scholar]
  21. Chakraborty, S., & Izaguirre, I. 2016, ArXiv e-prints [arXiv:1607.03361] [Google Scholar]
  22. Chartas, G., Brandt, W. N., Gallagher, S. C., & Garmire, G. P. 2002, ApJ, 579, 169 [NASA ADS] [CrossRef] [Google Scholar]
  23. Chevalier, R. A. 1976, ApJ, 207, 872 [NASA ADS] [CrossRef] [Google Scholar]
  24. Chevalier, R. A. 1982, ApJ, 258, 790 [NASA ADS] [CrossRef] [Google Scholar]
  25. Chevalier, R. A. 1998, ApJ, 499, 810 [NASA ADS] [CrossRef] [Google Scholar]
  26. Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&A, 562, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Di Mauro, M., Calore, F., Donato, F., Ajello, M., & Latronico, L. 2014, ApJ, 780, 161 [NASA ADS] [CrossRef] [Google Scholar]
  28. Drury, L. 1983, Space Sci. Rev., 36, 57 [NASA ADS] [CrossRef] [Google Scholar]
  29. Faucher-Giguère, C.-A., & Quataert, E. 2012, MNRAS, 425, 605 [NASA ADS] [CrossRef] [Google Scholar]
  30. Ferrara, A., Pettini, M., & Shchekinov, Y. 2000, MNRAS, 319, 539 [NASA ADS] [CrossRef] [Google Scholar]
  31. Feruglio, C., Maiolino, R., Piconcelli, E., et al. 2010, A&A, 518, L155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Feruglio, C., Fiore, F., Carniani, S., et al. 2015, A&A, 583, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Fields, B. D., Pavlidou, V., & Prodanović, T. 2010, ApJ, 722, L199 [NASA ADS] [CrossRef] [Google Scholar]
  34. Fiore, F., Feruglio, C., Shankar, F., et al. 2017, A&A 601, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Fornasa, M., & Sánchez-Conde, M. A. 2015, Phys. Rep., 598, 1 [NASA ADS] [CrossRef] [Google Scholar]
  36. Franco, J., Tenorio-Tagle, G., Bodenheimer, P., & Rozyczka, M. 1991, PASP, 103, 803 [NASA ADS] [CrossRef] [Google Scholar]
  37. García-Burillo, S., Combes, F., Usero, A., et al. 2014, A&A, 567, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Gatti, M., Lamastra, A., Menci, N., Bongiorno, A., & Fiore, F. 2015, A&A, 576, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Giallongo, E., Menci, N., Fiore, F., et al. 2012, ApJ, 755, 124 [NASA ADS] [CrossRef] [Google Scholar]
  40. Harrison, C. M., Alexander, D. M., Mullaney, J. R., et al. 2016, MNRAS, 456, 1195 [NASA ADS] [CrossRef] [Google Scholar]
  41. Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 731 [NASA ADS] [CrossRef] [Google Scholar]
  42. Inoue, Y. 2011, ApJ, 733, 66 [NASA ADS] [CrossRef] [Google Scholar]
  43. Kauffmann, G. 1996, MNRAS, 281, 475 [NASA ADS] [CrossRef] [Google Scholar]
  44. Kauffmann, G., & Charlot, S. 1998, MNRAS, 294, 705 [NASA ADS] [CrossRef] [Google Scholar]
  45. Kelner, S. R., Aharonian, F. A., & Bugayov, V. V. 2006, Phys. Rev. D, 74, 034018 [NASA ADS] [CrossRef] [Google Scholar]
  46. Keshet, U., Waxman, E., Loeb, A., Springel, V., & Hernquist, L. 2003, ApJ, 585, 128 [NASA ADS] [CrossRef] [Google Scholar]
  47. King, A. 2003, ApJ, 596, L27 [NASA ADS] [CrossRef] [Google Scholar]
  48. King, A. R., & Pounds, K. A. 2003, MNRAS, 345, 657 [NASA ADS] [CrossRef] [Google Scholar]
  49. King, A., & Pounds, K. 2015, ARA&A, 53, 115 [Google Scholar]
  50. King, A. R., Zubovas, K., & Power, C. 2011, MNRAS, 415, L6 [NASA ADS] [CrossRef] [Google Scholar]
  51. Lacey, C., & Cole, S. 1993, MNRAS, 262, 627 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lacki, B. C., Thompson, T. A., & Quataert, E. 2010, ApJ, 717, 1 [NASA ADS] [CrossRef] [Google Scholar]
  53. Lacki, B. C., Horiuchi, S., & Beacom, J. F. 2014, ApJ, 786, 40 [NASA ADS] [CrossRef] [Google Scholar]
  54. Lamastra, A., Menci, N., Maiolino, R., Fiore, F., & Merloni, A. 2010, MNRAS, 405, 29 [NASA ADS] [Google Scholar]
  55. Lamastra, A., Menci, N., Fiore, F., & Santini, P. 2013a, A&A, 552, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Lamastra, A., Menci, N., Fiore, F., et al. 2013b, A&A, 559, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Lamastra, A., Fiore, F., Guetta, D., et al. 2016, A&A, 596, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Lapi, A., Cavaliere, A., & Menci, N. 2005, ApJ, 619, 60 [NASA ADS] [CrossRef] [Google Scholar]
  59. Lucarelli, F., Pittori, C., Verrecchia, F., et al. 2017, ApJ, 846, 121 [NASA ADS] [CrossRef] [Google Scholar]
  60. Madau, P., & Rees, M. J. 2001, ApJ, 551, L27 [NASA ADS] [CrossRef] [Google Scholar]
  61. Makiya, R., Totani, T., & Kobayashi, M. A. R. 2011, ApJ, 728, 158 [NASA ADS] [CrossRef] [Google Scholar]
  62. Marconi, A., Risaliti, G., Gilli, R., et al. 2004, MNRAS, 351, 169 [NASA ADS] [CrossRef] [Google Scholar]
  63. McBride, J., Quataert, E., Heiles, C., & Bauermeister, A. 2014, ApJ, 780, 182 [NASA ADS] [CrossRef] [Google Scholar]
  64. Menci, N., Cavaliere, A., Fontana, A., Giallongo, E., & Poli, F. 2002, ApJ, 575, 18 [NASA ADS] [CrossRef] [Google Scholar]
  65. Menci, N., Fontana, A., Giallongo, E., & Salimbeni, S. 2005, ApJ, 632, 49 [NASA ADS] [CrossRef] [Google Scholar]
  66. Menci, N., Fontana, A., Giallongo, E., Grazian, A., & Salimbeni, S. 2006, ApJ, 647, 753 [NASA ADS] [CrossRef] [Google Scholar]
  67. Menci, N., Fiore, F., Puccetti, S., & Cavaliere, A. 2008, ApJ, 686, 219 [NASA ADS] [CrossRef] [Google Scholar]
  68. Menci, N., Gatti, M., Fiore, F., & Lamastra, A. 2014, A&A, 569, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 [Google Scholar]
  70. Murase, K., Ahlers, M., & Lacki, B. C. 2013, Phys. Rev. D, 88, 121301 [NASA ADS] [CrossRef] [Google Scholar]
  71. Murase, K., Guetta, D., & Ahlers, M. 2016, Phys. Rev. Lett., 116, 071101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  72. Narayan, C. A., & Jog, C. J. 2002, A&A, 394, 89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Nims, J., Quataert, E., & Faucher-Giguère, C.-A. 2015, MNRAS, 447, 3612 [NASA ADS] [CrossRef] [Google Scholar]
  74. Ostriker, J. P., & McKee, C. F. 1988, Rev. Mod. Phys., 60, 1 [NASA ADS] [CrossRef] [Google Scholar]
  75. Padovani, P., Resconi, E., Giommi, P., Arsioli, B., & Chang, Y. L. 2016, MNRAS, 457, 3582 [NASA ADS] [CrossRef] [Google Scholar]
  76. Pounds, K. A., Reeves, J. N., King, A. R., et al. 2003, MNRAS, 345, 705 [NASA ADS] [CrossRef] [Google Scholar]
  77. Reeves, J. N., O’Brien, P. T., & Ward, M. J. 2003, ApJ, 593, L65 [NASA ADS] [CrossRef] [Google Scholar]
  78. Resconi, E., Coenders, S., Padovani, P., Giommi, P., & Caccianiga, L. 2017, MNRAS, 468, 597 [NASA ADS] [CrossRef] [Google Scholar]
  79. Reynolds, S. P. 2008, ARA&A, 46, 89 [Google Scholar]
  80. Robishaw, T., Quataert, E., & Heiles, C. 2008, ApJ, 680, 981 [NASA ADS] [CrossRef] [Google Scholar]
  81. Rupke, D. S. N., & Veilleux, S. 2011, ApJ, 729, L27 [NASA ADS] [CrossRef] [Google Scholar]
  82. Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 [NASA ADS] [CrossRef] [Google Scholar]
  83. Shen, Y. 2016, ApJ, 817, 55 [NASA ADS] [CrossRef] [Google Scholar]
  84. Stecker, F. W. 1970, Ap&SS, 6, 377 [NASA ADS] [CrossRef] [Google Scholar]
  85. Stecker, F. W., & Venters, T. M. 2011, ApJ, 736, 40 [NASA ADS] [CrossRef] [Google Scholar]
  86. Stecker, F. W., Malkan, M. A., & Scully, S. T. 2007, ApJ, 658, 1392 [NASA ADS] [CrossRef] [Google Scholar]
  87. Stecker, F. W., Scully, S. T., & Malkan, M. A. 2016, ApJ, 827, 6 [NASA ADS] [CrossRef] [Google Scholar]
  88. Tamborra, I., Ando, S., & Murase, K. 2014, J. Cosmol. Astropart. Phys., 9, 043 [Google Scholar]
  89. Tatischeff, V. 2008, ArXiv eprint [arXiv:0804.1004] [Google Scholar]
  90. The IceCube Collaboration, Aartsen, M. G., Abraham, K., et al. 2015, ArXiv eprint [arXiv:1510.05222] [Google Scholar]
  91. Thompson, T. A., Quataert, E., Waxman, E., Murray, N., & Martin, C. L. 2006, ApJ, 645, 186 [NASA ADS] [CrossRef] [Google Scholar]
  92. Tombesi, F., Sambruna, R. M., Reeves, J. N., et al. 2010, ApJ, 719, 700 [NASA ADS] [CrossRef] [Google Scholar]
  93. Tombesi, F., Cappi, M., Reeves, J. N., et al. 2013, MNRAS, 430, 1102 [Google Scholar]
  94. Tombesi, F., Meléndez, M., Veilleux, S., et al. 2015, Nature, 519, 436 [Google Scholar]
  95. van der Kruit, P. C., & Freeman, K. C. 2011, ARA&A, 49, 301 [NASA ADS] [CrossRef] [Google Scholar]
  96. van Weeren, R. J., Brunetti, G., Brüggen, M., et al. 2016, ApJ, 818, 204 [NASA ADS] [CrossRef] [Google Scholar]
  97. Vazza, F., Eckert, D., Brüggen, M., & Huber, B. 2015, MNRAS, 451, 2198 [NASA ADS] [CrossRef] [Google Scholar]
  98. Vazza, F., Brüggen, M., Wittor, D., et al. 2016, MNRAS, 459, 70 [NASA ADS] [CrossRef] [Google Scholar]
  99. Wang, X., & Fields, B. D. 2016, ApJ, submitted [arXiv:1612.07290] [Google Scholar]
  100. Wang, X., & Loeb, A. 2015, MNRAS, 453, 837 [NASA ADS] [CrossRef] [Google Scholar]
  101. Wang, X., & Loeb, A. 2016a, Nat. Phys., 12, 1116 [CrossRef] [Google Scholar]
  102. Wang, X., & Loeb, A. 2016b, J. Cosmol. Astropart. Phys., 12, 012 [NASA ADS] [CrossRef] [Google Scholar]
  103. Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, R. 1977, ApJ, 218, 377 [NASA ADS] [CrossRef] [Google Scholar]
  104. Yoast-Hull, T. M., Gallagher, III, J. S., Zweibel, E. G., & Everett, J. E. 2014, ApJ, 780, 137 [NASA ADS] [CrossRef] [Google Scholar]
  105. Yu, Q., & Tremaine, S. 2002, MNRAS, 335, 965 [NASA ADS] [CrossRef] [Google Scholar]
  106. Zakamska, N. L., Hamann, F., Pâris, I., et al. 2016, MNRAS, 459, 3144 [NASA ADS] [CrossRef] [Google Scholar]
  107. Zubovas, K., & King, A. 2012, ApJ, 745, L34 [Google Scholar]
  108. Zubovas, K., & Nayakshin, S. 2014, MNRAS, 440, 2625 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

Schematic representation of the interaction of CR protons accelerated in the AGN blast wave with ISM protons. The shock radius Rs(t) expands outwards compressing the swept gas into a thin shell, and leaving a cavity inside. The γ-ray emission from neutral pion decays occurs in the regions outside the yellow cones where the lines of sight intercept the galactic gas that has not been swept out by the blast wave.

In the text
thumbnail Fig. 2

Distribution of simulated AGN as a function of and AGN bolometric luminosity for our fiducial model. The contours correspond to equally spaced values of the density (per Mpc3) of objects in a given LAGN bin in logarithmic scale: from 10-5 for the lightest filled region to 10-2 for the darkest.

In the text
thumbnail Fig. 3

γ-ray spectrum of an AGN with LAGN = 7 × 1044 erg/s hosted in a halo of mass 1012M at z = 0.1. Energy spectrum parameters are set at p = 2 (black lines), p = 2.2 (red lines), B = BISM (solid lines), and B = Bshock (dashed lines).

In the text
thumbnail Fig. 4

The contribution to the EGB from AGN winds (blue) and starburst galaxies (red) predicted by our SAM. Fiducial model parameters are set at: hd = 100 pc, ηp = 0.1, ϵAGN = 0.01, and ϵSN = 0.1. Accelerated proton spectral index equal to p = 2.2 is assumed for both AGN winds and starburst galaxies. The blue and red shaded bands represent the uncertainty related to the EBL model adopted (Stecker et al. 2016). The data points are the Fermi-LAT measurement of the EGB, with black and grey vertical error bars indicating EGB statistical and systematic uncertainties, respectively (Ackermann et al. 2015).

In the text
thumbnail Fig. 5

Upper panel: the contribution to the EGB from AGN winds predicted by our SAM is compared to the contribution provided by the blazar population (Ajello et al. 2015). Model predictions are shown for our fiducial model and for different accelerated proton spectral index: p = 2 (dotted blue lines), p = 2.2 (solid blue lines), and p = 2.4 (dashed blue lines). The red shaded region represent the sum of the AGN wind and blazar contributions to the EGB. The data points are as in Fig. 4. Lower panel: fraction of EGB provided by AGN winds and blazars.

In the text
thumbnail Fig. 6

EGB spectrum. The blue and red shaded regions correspond to AGN winds and starburst galaxies predicted by our SAM, respectively. The purple shaded region corresponds to the contribution to the EGB from AGN winds estimated by Wang & Loeb (2016a). The cyan and gold lines indicate the estimate for quiescent and starburst galaxies, respectively, predicted by the Makiya et al. (2011) SAM. The green shaded region shows the contribution to the EGB from blazars as estimated by Ajello et al. (2015). The data points are as in Fig. 4.

In the text
thumbnail Fig. 7

Neutrino spectrum of an AGN with LAGN = 7 × 1044 erg/s hosted in a halo of mass 1012M at z = 0.1. Model parameter as in Fig. 3.

In the text
thumbnail Fig. 8

Cumulative γ-ray (left) and neutrino (right) background from AGN winds predicted by our SAM. Model predictions are shown for our fiducial model and for different AGN wind spectral energy parameters: p = 2.2 (purple lines), p = 2.3 (blue lines), B = BISM (solid lines), and B = Bshock dashed lines. The circles represent Fermi-LAT data (Ackermann et al. 2015). The squares show the differential model of IceCube neutrino data, while purple dotted lines represent the power-law models (Aartsen et al. 2015).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.