arXiv:1402.0383v1 [astro-ph.HE] 3 Feb 2014 i hog ycrto msin u loin also in but notably emission, emissions, synchrotron popula extended large-scale through various such by dio of attested existence is The tions Way). Milky the for hycnha n oietegs seilya h eteof centre the at especially c gas, formation star the flowing. impacting thereby are ionise sig- clouds, molecular they and can dense which and heat in ISM can environment the They more the with are interact modify rays nificantly they e cosmic because Galactic side tion. mostly simple is a than This astrophysics. 200 energy al. et Strong (see CRs of cycle wi review). life us a provide the for can of data picture astronomical better other a of inter frame their the and in transport, tion and com acceleration the CR of of signatures processes are emissions These inelast collisions. clear and Bremsstrahlung, scattering, inverse-Compton ag ubro ceeaineioe o h re of order the (of episodes acceleration ov of f volume number galactic responsible large the are a in accumulation and and medium, of confinement intergalactic their transport to the m sources determine and from interactions radiation, These turbulence, fields. ray the gas, netic of cosmic (ISM): components galactic various medium the terstellar These with nebulae. interact CRs) wind (hereafter pulsar collidin explosions, or supernova binaries, as such evolution, associat phenomena stellar violent in of accelerated populations particles getic large-scale harbour galaxies Star-forming Introduction 1. uut1,2018 13, August srnm Astrophysics & Astronomy nesadn R srlvn omn ed eodhigh- beyond fields many to relevant is CRs Understanding ntttd ehrh nAtohsqee Plan´etologie, U et Astrophysique en Recherche de Institut odtos hr s oee,n nvra eainbetwe relation universal no however, is, There conditions. di ik a hnetaoae oohrsses rbbybeca probably systems, other to extrapolated when Way Milky ytesatritoue ydi by introduced scatter the by eevd 3Dc2013 Dec 23 Received: h di the Conclusions. the from inferred trend and normalisation the distributi gas several kpc, 40 to Results. 4 ta from sizes were with parameters simulated CR-related Most galaxies. star-forming Methods. Aims. Context. lblglci rpris ihtega fetbihn w establishing of goal the with properties, galactic global xetdfo obigthe doubling from expected as(R)i aiu lblcniin a il informati yield can conditions global various between in correlation (CRs) a rays for found was Evidence e words. Key ff so cee igedi single a scheme, usion ff hswr drse h usino h cln asta can that laws scaling the of question the addresses work This rn msinpoessadpril ouain,ado t on and populations, particle and processes emission erent h vlto fthe of evolution The Fermi eeoe Dmdlfrtenntemleisosfo ste from emissions non-thermal the for model 2D a developed I neselrgmaryeiso rmcsi rays cosmic from emission gamma-ray Interstellar ceeaino atce omcry am as galaxi rays: Gamma – rays Cosmic – particles of Acceleration h current The / A bevtoso trfrigglxe nthe in galaxies star-forming of observations LAT ff c fteedsae fselrevolu- stellar of stages end the of ect / cetd a 2014 Jan 8 Accepted: aucitn.aa23329 no. manuscript Fermi Fermi ff ff γ so coe usion rn aatcgoa rprisadtedwtr in downturn the and properties global galactic erent rylmnst vrte10kV10TVrnei presented, is range TeV keV-100 100 the over luminosity -ray / A ouainsuyde o alfrmjrmdfiain of modifications major for call not does study population LAT / A xouetm n ae rmosriga e nriswi energies TeV at observing from later and time exposure LAT nsa-omn galaxies star-forming in ffi in,adteasmto htCseprec ag-cl v large-scale experience CRs that assumption the and cient, Fermi γ γ rylmnst n rcr ftesa omto activity. formation star the of tracers and luminosity -ray ry through -rays / A ouainsuyoe oto h F ag.Ti sobtai is This range. SFR the of most over study population LAT dwith ed g-wind ehrtecreteprmna aai h e ag a ec be can range GeV the in data experimental current the hether ∼ cnu- ic PS preta- nhigh-energy en bined naotteroii n rnpr nteitrtla medi interstellar the in transport and origin their about on ondi- ener- n,adsa omto ae SR oeigsxodr fma of orders six covering (SFR) rates formation star and ons, e rmMlyWysuis n ag ubro aaiswer galaxies of number large a and studies, Way Milky from ken ABSTRACT CRs .Martin P. 10 / ag- NS M57,308Tuos ee ,France 4, cedex Toulouse 31028 UMR5277, CNRS, ra- in- or th er 7, s h netite r tl o ag.Adtoa const Additional large. too still are uncertainties the use 5 s - etasto ewe rnpr eie.Temdlcnrepr can model The regimes. transport between transition he ∼ rletra aaista hn tpoo energies photon at shine that galaxies external eral nry(VHE) energy crucial. is galaxies star-forming of emission non-therma and the a tent understanding et reasons, (Delahaye these matter all ex dark For more as 2010). of such Universe signatures the the class of hide from components CRs can last, objects and throu astrophysical 2004); fields cal play al. magnetic et may galactic (Hanasz CRs process large-scale 2013); dynamo shaping al. in et role (Hanasz a times cosmic contri over e can feedback mation they important 2006); an al. outflows, large-scale et to (Ptuskin consequ transport with own medium, their interstellar on the of turbulence di the ter their through 2010); (Papadopoulos tions as facnrlbakhl ciiy si h aefrthe for case be- the not is (and as ISM activity, their hole detected with black of interacting majority central vast a CRs of of cause result a as 02 cir ta.20) trfrigglxe a o be now of can class al. emerging galaxies et Abramowski an Star-forming as 2009; 2009). considered al. al. et et Acciari (Acero 2012; respectively NGC25 of M82, energies tele- TeV at Cherenkov and same detection ground-based the the VERITAS permitted see In have and activity; scopes 2012a). HESS hole al. the et black Ackermann time, central 2010c,b,d,a; by al. et be- contaminated Abdo systems be emis- other whose may NGC1068, two sion (SMC, and with (NGC4945 energies candidates NGC253), GeV possible and ing at M82, galaxies M31, external LMC, five of detection the Fermi .-0 e ag aemd osbeafis ouainstudy population first a possible made have range GeV 0.1-100 eepce o h interstellar the for expected be h urn eeaino iheeg H)advery-high- and (HE) high-energy of generation current The γ / es A pc-on arceto eecp a permitted has telescope pair-creation space-borne LAT rylmnst n trfrainatvt,a illustrate as activity, formation star and luminosity -ray d-tt Rppltositrcigwt h S in ISM the with interacting populations CR ady-state γ ryeiso ttelwend. low the at emission -ray γ rytlsoe a nbe st eetsev- detect to us enabled has telescopes -ray ihepai ntecnrbto of contribution the on emphasis with γ h rnpr ceefrCsi the in CRs for scheme transport the ryeitn xenlglxe) The galaxies). external -ray-emitting hthe th γ ryeiso safnto of function a as emission -ray hrno eecp Array Telescope Cherenkov lm-vrgdinterstellar olume-averaged tdigglci cosmic galactic Studying γ rysucs ihano- a with sources, -ray m(ISM). um e ihaplain a with ned ff so,Cscnal- can CRs usion, ansmybe may raints onstraining. gnitude. ff c nsa for- star on ect
c > S 2018 ESO then e oduce 0 MeV 100 d . . con- l ences ha gh bute otic i- l. 3 1 Martin P.: Interstellar gamma-ray emission from star-forming galaxies table increase of the number of detected objects from the previ- the case of the high-energy band. Last, the impact of the model ous generation of instruments (only the LMC was detected by parameters is discussed, and the connection to the far-infrared - CGRO/EGRET; see Sreekumar et al. 1992). radio correlation is addressed. Even with such a limited sample, it is tempting to perform a population study to find clues about the main parameters regu- 2. Model lating CR populations. From the Fermi/LAT sample of five ob- jects plus an estimate for the global emission of the Milky Way The model is based on a modified version of the GALPROP pub- and an upper limit for M33, the work of Abdoetal. (2010a) lic code1, originally designed for the Milky Way and adapted showed first hints for a slightly non-linear correlation between here to model any star-forming galaxy. For the most part, this γ-ray luminosity and star formation rate (SFR). A more system- was achieved by implementing general models for the interstel- atic study by Ackermann et al. (2012a), including upper limits lar medium components, as described below. This generic model for about sixty luminous and ultra-luminous infrared galaxies, uses as input and for calibration some characteristics of a typical found further evidence for a quasi-linear correlation between γ- GALPROP modelling of the Milky Way in the plain diffusion ray luminosity and tracers of the star formation activity. Such a case: Model z04LMPDS in Table 1 of Strong et al. (2010), re- relationship can be expected to first order, because the produc- ferred to in the following as the reference GALPROP run. tion of CRs is expected to scale with the rate of massive star All the synthetic galaxies discussed in this work are two- formation; but it is surprising that the correlation is so close to dimensional in space with cylindrical symmetry, and are defined linear, over such a wide range of galactic properties, from dwarfs by a disk radius Rmax and a halo height zmax. The model is there- like the SMC to interacting systems like Arp 220. fore better suited to simulate grand-design spiral galaxies than The interpretation of this correlation therefore raised the strongly irregular or interacting galaxies, but it can provide in- question of the scaling laws that can be expected for the dif- sights into the effects of galaxy size and star formation condi- fuse γ-ray emission as a function of global galactic properties. tions on non-thermal emissions. Clarifying this at the theoretical level is required to establish whether the current experimental data can be constraining in any 2.1. Gas way. This is all the more needed because the sample of detected objects is small, the coverageover the GeV-TeV range is uneven, The interstellar gas distribution was modelled through simple and the spectral characterization is still rather limited. In addi- analytical functions of the galactocentric radius R and height tion, the determination of global properties of external systems, above the galactic plane z. I considered only atomic and molec- such as total mass or star formation rate, is in itself a challenge, ular gas, and neglected ionised gas. A ratio of He/H=0.11 was and there is most likely some intrinsic scatter in the actual phys- assumed for the interstellar gas. ical correlations. Three types of gas distributions are discussed in the follow- That star-forming galaxies might become a class of γ- ing: a disk of atomic hydrogen with exponential density profile ray sources was anticipated early, well before the launch in radius and height (three parameters: scale radius, scale height, of the Fermi satellite (see for instance V¨olk et al. 1996; and peak atomic gas density); a central core of molecular gas Pavlidou & Fields 2001; Torreset al. 2004). The interest in with uniform density (three parameters: core radius, core height, the subject extended beyond the study of nearby objects, be- and molecular gas density); a ring of molecular gas with uniform cause this emission cumulated over cosmic time might ac- density (four parameters: ring inner and outer radii, ring height, count for a significant fraction of the isotropic diffuse γ-ray and molecular gas density). These profiles are intended to re- background (Pavlidou & Fields 2002; Thompson et al. 2007; produce a typical star-forming galaxy, where star formation can Ackermannet al. 2012a). Before and after the launch of the operate at a low surface density in a disk, and/or be concentrated Fermi satellite, the γ-ray emission from individual starbursts in more compact regions like central cores (such as in M82 and was studied by several authors, with the focus being mainly NGC253),or in spiral arms or rings(such as in the Milky Way or on M82 and NGC253 (Persic et al. 2008; de Cea del Pozo et al. M31).The different models and their parameters are summarised 2009; Lacki et al. 2011; Paglione & Abrahams 2012). below. The purpose of this paper is to provide a more global pic- An important aspect of the model, as becomes obviousin the ture of the evolution of the γ-ray interstellar emission from star- following, is that the gas content and distribution drives many forming galaxies, in the context of the recent population study other components, such as the magnetic field model, the infrared conducted on the basis of Fermi/LAT observations. Instead of energy density of the interstellar radiation field, or the CR source trying to fit a model to a particular object, I defined a global distribution. framework by fixing as many cosmic-ray-related parameters as possible from Milky Way studies, leaving only a few global quantities as independent variables. Here, I present a generic 2.2. Magnetic field model for the transport of CRs and the associated non-thermal The magnetic field strength is assumed to depend at each point interstellar emissions in star-forming galaxies. Simulating sys- on the vertical gas column density, according to tems with various sizes and gas distributions thus allows inves- a Σgas(R) −|z| tigating the impact of global galactic parameters on the γ-ray z B(R, z) = B0 e B , (1) output. In a first part, I introduce the different components and Σ0 ! related assumptions of the model, together with the set of syn- 1 thetic galaxies that was considered. Then, the γ-ray emission GALPROP is a software package for modelling the propaga- tion of CRs in the Galaxy and their interactions in the ISM. It al- over 100keV-100TeV for this series of synthetic galaxies is pre- lows simultaneous and coherent predictions of the CR population sented, and its evolution with global properties is analysed. In in the interplanetary and interstellar medium, gamma-ray emission, a subsequent part, the scaling of γ-ray luminosity with global and recently also synchrotron radiation intensity and polarization. See galactic properties is derived for soft, high-energy, and very- http://galprop.stanford.edu/ and http://sourceforge.net/projects/galprop/ high-energy γ-ray ranges, and is compared with observations in for the latest developments.
2 Martin P.: Interstellar gamma-ray emission from star-forming galaxies where Σgas is the gas surface density. Such a scaling of the mid- zIR=2kpc above and below the galactic plane. The latter value plane magnetic field value can be justified by different theoreti- was set by comparison with the reference GALPROP run, as cal expectations, from hydrostatic equilibrium to turbulent mag- explained in Sect. 2.6. The approximation adopted implements netic field amplification (see Schleicher & Beck 2013). The ver- an IR energy density in the galactic plane that has a minimum tical profile is assumed to be exponential, and I adopted a mag- value for the threshold of star formation, and then scales with netic field scale height equal to the halo size (following Sun et al. the star formation rate according to the term in parentheses2.I 2008, and quite similar to the reference GALPROP run; the influ- neglected the cirrus emission arising from evolved star heating, ence of this assumption is examined in Sect. 3.4). The exponent which can increase the IR energy density by up to a factor of 2 a was set to a value of 0.6 to match the observed far-infrared - (Kennicutt et al. 2009), but I tested the impact of varying the IR radio continuum correlation, as discussed in Sect. 3.6. The nor- field intensity, as discussed in Sect. 3.4. −2 malisation doublet (B0, Σ0) was set to (5 µG, 4M⊙ pc ), from a The mid-plane value for the starlight energy density distribu- comparison with the reference GALPROP run, as explained in tion was assumed to be the IM83 field scaled by the dust-weighted Sect. 2.6. mean starlight intensity hUi, I did not consider a global topology model for the magnetic field, or a separation into regular and turbulent components, be- 4π −|z| u (R, z) = hUi E e zUVO , (4) cause the diffusion of CRs is treated in a homogeneous and UVO c M83 isotropic way. A recent GALPROP modelling using a more com- where EM83 = IM83(λ)dλ (5) plete magnetic field model can be found in Orlando & Strong Z (2013). Theoretically, anisotropic diffusion along field lines can be expected, depending on the ratio of turbulent field to regu- With the parameters adopted (see the appendix), this mean lar field strength, and lead for instance to channelling and clus- starlight intensity is hUi=3.4, so it is about 3 times stronger ev- tering of particles preferentially along spiral arms (Dosch et al. erywhere in the galactic disk than in the local Galactic environ- 2011). But, the reality and impact of this compared with the ho- ment. Getting the same field everywhere in the galactic disk is mogeneous and isotropic approximation is poorly constrained probably not realistic. The stellar photon field experienced by at the experimental level (see the references and discussion in CRs is expected to be higher in regions of intense star forma- Jaffe et al. 2013). tion, especially because this is where they are released into the ISM, but these effects occur on scales that are similar to or below the typical spatial resolution of the model. The energy density of 2.3. Radiation field the starlight component is taken to decrease over a scale height The interstellar radiation field model has three main compo- zUVO=1kpc above and below the galactic plane, according to a nents: the cosmic microwave background (CMB), the infrared comparison with the reference GALPROP run (see Sect. 2.6). dust emission (IR), and the stellar emission (UVO). The first contribution is implemented straightforwardly as a blackbody 2.4. Cosmic-ray source with a temperature of TCMB=2.7K. The other two components are dependent on the particular galactic setting and ideally result The main sources of CRs are thought to be supernova remnants, from a complex radiation transfer problem involving a variety of possibly with some contribution of other objects such as pulsars, stellar sources and a mixture of dust grains, each with specific pulsar wind nebula, or colliding-wind binaries. Since most of spatial distributions. these are end products of massive star evolution with lifetimes The interstellar radiation field model (ISRF) is based on the shorter than that of their progenitors, the CR source distribution work by Draine & Li (2007). Interstellar dust is irradiated bya was assumed by default to follow the star formation distribution. distribution of starlight intensities with the IM83(λ) spectrum de- The latter is determined from the gas distribution using the em- termined by Mathis et al. (1983) for the local Galactic environ- pirical Schmidt-Kennicutt relation ment and scaled by a dimensionless factor U. The resulting dust Σ = ΣN Σ Σ emission Idust(λ) in the 1-10000 µm wavelength band and nor- SFR A gas for gas > thres, (6) malized per hydrogen atom is available for various combinations of the model parameters (see the appendix). where ΣSFR and Σgas are the star formation rate and −1 −2 The infrared intensity distribution over the galactic volume atomic+molecular gas surface density (in M⊙ yr kpc and −2 −4 ideally results from integrating at each point the contribution M⊙ pc units), and with A = 2.510 and N = 1.4 the canonical −2 of all sources, possibly affected by absorption, reemission, and values (Kennicutt 1998). I used a threshold at Σthres =4M⊙ pc , scattering (see for instance Ackermann et al. 2012b). As a sim- which is a lower limit from spatially resolved observations of plification, the energy density of the infrared dust component in galaxies in Kennicutt (1998). The relation is observed to hold the ISRF model was computed as from galactic averages down to kpc averages. At each galacto- centric radius, the vertical distribution of the star formation rate N 4π Σgas(R) −|z| was taken to be proportional to the gas density. z uIR(R, z) = Edust Σthres e IR , (2) The total CR source luminosity was assumed by default to be c Σthres ! proportional to the total star formation rate, using as normalisa- where Edust = Idust(λ)dλ, (3) tion the reference GALPROP run. The following CR input power Z 2 where Σ is a threshold gas surface density for star formation With the adopted parameterization for the infrared dust emission, thres the conversion factor between star formation rate and infrared lumi- and N the index of the Schmidt-Kennicutt relation (see the defi- −10 −1 nosity in the 8-1000 µm band is 1.3 10 M⊙ yr / L⊙. This can be nitions below). The quantity c is the speed of light. The intensity −10 −1 compared with 1.7 10 M⊙ yr / L⊙ in the original formulation by at each point (R, z) in the volume is thus assumed to be domi- Kennicutt (1998) using a Salpeter initial mass function (Salpeter 1955), −10 −1 nated by dust emission from the same galactocentric radius R, or with 1.3 10 M⊙ yr / L⊙ when using a Chabrier initial mass func- and the energy density is taken to decrease over a scale height tion (Chabrier 2003).
3 Martin P.: Interstellar gamma-ray emission from star-forming galaxies integrated over 100MeV-100GeV was used: Table 1. List of the galaxy configurations in the initial sample SFR = , QCRp,e QMW −1 (7) 1.9M⊙ yr Rmax zmax HI H2 40 −1 where QMW = 7.1010 erg s for CRp 2.0 1.0 disk with RH=1.0 core with Rcore=0.1 39 −1 - - - core with Rcore=0.2 QMW = 1.0510 erg s for CRe , 5.0 1.0 disk with RH=2.5 core with Rcore=0.2 where CRp stands for cosmic-ray hydrogen and helium nuclei, - - - core with Rcore=0.5 and CRe stands for cosmic-ray electrons. The SFR normalisation - - - core with Rcore=1.0 was obtained by integrating Eq. 6 over the Galaxy given the gas - - - ring with (Ri, Ro) = (1.5, 2.0) = distribution of the original GALPROP model; the result lies in - - - ring with (Ri, Ro) (2.0, 3.0) 10.0 2.0 disk with R =5.0 core with R =0.5 the range of current estimates (Diehl et al. 2006), especially if H core - - - core with R =1.0 one considers that the scatter in the Schmidt-Kennicutt relation core - - - ring with (Ri, Ro) = (2.0, 3.0) can reach an order of magnitude (Kennicutt 1998). The relation - - - ring with (R , R ) = (4.0, 6.0) between total star formation rate and CR luminosity is probably i o 20.0 4.0 disk with RH=10.0 core with Rcore=1.0 far from straightforward because of a possible dependence of the - - - ring with (Ri, Ro) = (4.0, 6.0) stellar initial mass function or particle acceleration process on - - - ring with (Ri, Ro) = (8.0, 12.0) the interstellar conditions, but taking this into account is beyond the scope of the present work. Note to the Table: The first two columns are the maximum radius and The CR injection spectra used in this work were consistently half-height of the galactic volume, respectively. The third column gives taken from the reference GALPROP run. Particles at source have the scale radius of the exponential distribution of atomic gas (the verti- a power-law distribution in rigidity with a break, following cal distribution is exponential with a scale height of 100 pc). The fourth column indicates whether the molecular gas distribution is a uniform −q dQCRp,e R core or uniform ring, and gives the corresponding core radius or in- = Q0 , (8) ner/outer radii (the vertical distribution is uniform with a half-height of dR R0 ! 100 pc). All lengths are in kpc units. For each model, 7 molecular gas −3 with q = 1.80 for R ≤ Rb densities were tested: nH2 =5, 10, 20, 50, 100, 200, 500 H2 cm , for the same peak atomic gas density of n =2 H cm−3. q = 2.25 for R > Rb H
Rb = 9GV for CRp Rb = 4GV for CRe , can be resolved by adopting inhomogeneousand anisotropic dif- fusion schemes (Evoli et al. 2012). This means that the same electron-to-proton ratio at injection For simplicity and as a starting point, the spatial transportof was used for all synthetic galaxies. This may introduceerrors be- CRs in the present model of star-forming galaxy was assumed cause acceleration can be expected to proceed differently in en- to occur by energy-dependent diffusion in a homogeneous and vironments where magnetic and radiation energy densities are a isotropic way. For the diffusion coefficient, I used the value from factor of 100-1000 higher (as is the case in the cores of starburst the reference GALPROP run galaxies). But these differences are currently poorly constrained 0.5 at the experimental level, and it is interesting enough to assess 28 R 2 −1 Dxx = 3.410 β cm s for R ≥ 4GV the effects of global galactic properties for a given CR injection 4GV spectrum. Moreover, the assumption that particle acceleration in 28 2 −1 Dxx = 3.410 cm s for R < 4GV, (9) other galaxies works similarly to what we think takes place in the Milky Way can be tested by comparison with the emerging where R is the particle rigidity and β = v/c, the ratio of particle class of γ-ray star-forming galaxies. velocity to speed of light. The spatial transport of CRs in other galaxies than the Milky Way may take place in a different way. In systems harbouring a 2.5. Transport large starburst region, the interstellar turbulence may be strongly At galactic scales, the spatial transport of CRs in the ISM is cur- enhanced and result in different diffusion conditions, or strong rently understood as being ruled by two main processes: advec- outflows out of the disk may be generated and advect particles tion in large-scale outflows such as galactic winds or superbub- away. I discuss these two possibilities in Sect. 3.4. bles breaking out into the halo, and diffusion on interstellar mag- netohydrodynamic turbulence (Strong et al. 2007). While the 2.6. Test and calibration former implies adiabatic energy losses, the latter may be asso- ciated with reacceleration of the CRs. Recent studies seem to The star-forming galaxy model described above has a total of dismiss the need for reacceleration, and show that a large set five undetermined parameters. One of these is the index of the of cosmic-ray-related observations can be reasonably well ac- magnetic field model, and it was set using the constraint of the counted for from a pure diffusion scheme with isotropic and ho- far-infrared - radio correlation (see Sect. 3.6). The other four mogeneous properties (Strong et al. 2011, note, however, that parameters are the normalisation doublet (B0, Σ0) for the mag- the good fit to the data is obtained by simultaneously tuning netic field model, and the scale heights zIR and zUVO for the ISRF the source and transport parameters). Several discrepancies be- model. These parameters were set by a comparison with the ref- tween such predictions and observations remain, such as the so- erence GALPROP modelling of the Milky Way, which is basi- called gradient problem, connected with the longitude distribu- cally what the code was optimised for. From the built-in axisym- tion of the diffuse interstellar γ-ray emission, and the large-scale metric atomic, molecular, and ionised gas distributions, models cosmic-ray anisotropy, which is observed to be smaller than ex- for the ISRF, magnetic field, and source distribution were de- pected at high energies. Several authors have shown that these rived using the assumptions detailed above. A typical run was
4 Martin P.: Interstellar gamma-ray emission from star-forming galaxies then performed with these modified conditions, and its results butions extending 100pc on either side of the galactic plane (it were compared with those of the original GALPROP calcula- was verified that using a smaller molecular gas scale height of tion. 50pc does not alter the trends in the evolution of the γ-ray lumi- The comparison was first made on the ISM conditions, using nosity with galactic properties). Not all combinations of molec- source-weighted and volume averages for the energy densities of ular gas distributions and galaxy sizes were tested, see the list the radiation and magnetic ISM components (these averages are of models summarised in Table 1. Then, for each molecular gas relevant to the conditions experienced by CRs at injection and profile, seven gas densities were tested: nH2 =5, 10, 20, 50, 100, −3 over their long-term transport, respectively). From this, the four 200, and 500H2 cm . The initial sample of models includes undefined parameterslisted above were set. Applied to the Milky 14 × 7 = 88 synthetic star-forming galaxies, not counting varia- Way, the generic model for star-forming galaxies results in radi- tions of some parameters to assess their importance. While some ation (respectively magnetic) energy densities with volume and configurations may be considered as good models for existing source-weighted averages that are different by about ∼20-30% systems such as the Milky Way or M82, others may be less real- (respectively ∼5-10%) from those of the original GALPROP istic, such as the (Rmax, zmax)=(20kpc, 4kpc) model with a large −3 4 −1 setup. Then, the total Galactic non-thermal outputs were com- ring of nH2 =500H2 cm gas, yielding an SFR of ∼ 10 M⊙ yr . pared and were found to agree very well, without any need for These were included to illustrate the effects of changing galactic further tuning of the model. The radio emission over 1MHz- parameters. 1THz differs at most by 6% at low frequencies3, while the γ-ray emission over 100keV-100TeV shows a maximum deviation of 2.8. Time scales 5%. When integrating the emission over >100MeV, the differ- ences in total luminosity are even smaller. The similarity of the To facilitate the interpretation of the results presented in the fol- results to that of the original GALPROP calculation applies not lowing,I providethe typical time scales for the variousprocesses only to the total emission but also to its components: inverse- involved in the transport of CRs. These hold for the assumptions Compton, Bremsstrahlung, and pion decay. used in the star-forming galaxy model (such as the spatial diffu- sion coefficient). −0.5 2 2.7. Setup 14 E z τdiff = 2.810 s (10) All the models discussed in this work are two-dimensional in 4GeV 1kpc! space with cylindrical symmetry, and are defined by a disk z V −1 τ = 3.11014 s (11) radius R and a halo half-height z , with a cell size of adv −1 max max 1kpc! 100kms 50pc in both dimensions. The particle energy grid runs from n −1 100MeV to 1PeV with ten bins per decade. The transport of = 15 H τpp 2.210 − s (12) CRs arises from diffusion for the spatial part, and includes ion- 1cm 3 n −1 isation, synchrotron, Bremsstrahlung, inverse-Compton scatter- = 14 H τBr 8.910 − s (13) ing, and hadronic interaction losses for the momentum part. For 1cm 3 the latter, the production of secondary electrons and positrons is −1 −1 15 E Urad taken into account and included in the complete transport and τIC = 9.910 s (14) 1GeV 1eVcm−3 radiation calculation. All solutions correspond to a steady state. E −1 U −1 For more details about the implementation of the physical pro- = 15 mag τSynch 9.910 − s (15) cesses and the inner workings of the code, see the GALPROP 1GeV 1eVcm 3 ! Explanatory Supplement available at the GALPROP website. E n −1 = 15 H Four different sizes were considered for the synthetic τion 4.810 − s (16) 1GeV 1cm 3 star-forming disk galaxies: (Rmax, zmax)=(2kpc, 1kpc), (5kpc, 1kpc), (10kpc, 2kpc), and (20kpc, 4kpc). The larger size cor- The time scales are those of the diffusion, advection, hadronic responds to the dimensions of the Milky Way in the reference interaction, Bremsstrahlung, inverse-Compton, synchrotron, and GALPROP run, and I kept the same aspect ratio for the smaller ionisation processes. For the spatial transport processes, the typ- sizes, except for the smaller model where it would have resulted ical length was taken at 1kpc, which is the order of magnitude in a too small halo. The size of the halo for the Milky Way and for the galactic halo size. The dependence on energy for the dif- how it evolves with galaxy size are unsettled questions, there- fusion time scale applies only for E > 4GeV. The hadronic in- fore this assumption of a nearly constant aspect ratio in galaxy teraction time scale was computed for an energy of 1GeV, with dimensions is challenged below (see Sects. 3.4, 3.5, and 3.6). a cross-section of 30mbarn and an inelasticity parameter of 0.5 For each size, the galaxy consisted of a disk of atomic gas (note that the cross-section increases slightly to reach 40 mbarn with an exponential distribution in R and z, with scale lengths at 1TeV, see Kelner et al. 2006). The Bremsstrahlung-loss time scale corresponds to an electron energy of 1GeV, and typically of RH = Rmax/2 and zH =100pc respectively, and a peak gas density of n =2Hcm−3. On top of that, different molecular has a logarithmic dependence on energy that was neglected in H the above formula (Schlickeiser 2002). The inverse-Compton- gas distributions were implemented: cores of radius Rcore =100, 200, 500pc, and 1kpc, and two rings with inner and outer radius loss time scale was computed in the Thomson approximation, with Urad corresponding to the radiation energy density. The (Ri, Ro) = (0.2 Rmax, 0.3 Rmax) and (0.4 Rmax, 0.6 Rmax), all distri- synchrotron-loss time scale is the same formula, with Umag cor- responding to the magnetic energy density. The ionisation-loss 3 Note that we compare the output of two different models for the unabsorbed synchrotron emission over a very large frequency range. At time scale was computed for protons in the ISM, using the low- the observational level, free-free absorption causes a turnover at low est value of the loss rate (which is obtained for a Lorentz fac- frequencies, and free-free and thermal dust emission dominate at high tor of ∼ 3) and neglecting its energy dependence (which leads frequencies. Eventually, unabsorbed non-thermal emission dominates to variation by a factor of at most 3 over the 100MeV-10GeV in the range ∼100 MHz-10 GHz. range).
5 Martin P.: Interstellar gamma-ray emission from star-forming galaxies
Because the ISRF and magnetic field components are de- 500pc core configuration exhibits the highest contrasts between −3 duced from the gas distribution in my model, it is useful to re- the low-density case (nH2 = 5H2 cm ) and the high-density −3 cast some of the above formulae as a function of molecular gas case (nH2 = 500H2 cm ): the luminosities increase because of density. For a given gas distribution, such as a molecular core or a higher conversion efficiency of particle energy to radiation, a ring of a given extent on top of the atomic disk, the molecu- which arises from a higher gas density and from a larger pro- lar gas density is the independent variable that drives the non- portion of CRs being injected in a higher-density medium. thermal properties. Accordingly, from Eqs. 1 and 2, and usinga The evolution of the luminosity in all three γ-ray bands, for thickness of 200pc for the molecular gas, one obtains an increase of a factor 100 of the molecular gas density, shows the following trends (I recall that this discussion holds for the n −1 = 15 H2 same input CR luminosity): τpp 1.110 − s (17) 1cm 3 Band I (100keV-10MeV): The emission is always domi- n −1 nated by IC, with Bremsstrahlung contributing at the 15-25% τ = 4.51014 H2 s (18) Br −3 level depending on gas density. The IC from primaries is al- 1cm E −1 n −1.4 most constant, which is consistent with primaries being IC loss- τ = . 16 H2 IC 1 810 −3 s (19) dominated (the increase in the radiation field density is balanced 1GeV 1cm by a decrease in the steady-state primaries population). The IC E −1 n −1.2 = 15 H2 from secondaries increases by a factor of about 10, which is due τSynch 5.010 − s (20) 1GeV 1cm 3 to a rise of the production rate of secondaries (the production E n −1 rate increases not as much as the gas density because low-energy τ = 2.41015 H2 s, (21) ion 1GeV 1cm−3 CRp progressively move from a diffusion-dominated regime to a loss-dominated regime). The net result is that emission in band where only the infrared component of the ISRF was considered I gradually becomes dominated by secondaries, at a level of 70- in the inverse-Compton time scale (it is by far the dominant com- 80% for the highest gas density (for observations and interpreta- = −3 ponent above nH2 10-20H2 cm ). tion of the diffuse soft γ-ray emission from the Milky Way, see Porter et al. 2008; Bouchet et al. 2011). Overall, this contribu- 3. Results tion of secondaries increases from the very-high-energy band to the soft-gamma-ray band (band III to I), which reflects mainly In this section, I present the results for the synthetic star-forming that secondaries are the particles with the steepest spectrum. galaxies introduced above. To allow the connection with obser- Band II (100MeV-100GeV): The emission is always dom- vations, the γ-ray luminosities are often given in three spectral inated by pion decay, while leptonic processes contribute at the bands: 100keV-10MeV (band I, soft gamma-rays), 100MeV- 20-30% level, with an increasing fraction of it coming from sec- 100GeV (band II, high-energy gamma-rays), 100GeV-100TeV ondaries (eventually reaching 75% for the highest gas densities). (band III, very-high-energy gamma-rays). Examples of cur- The evolution of the pion decay and Bremsstrahlung emission rent instruments covering these bands are INTEGRAL/SPI, with gas density flattens, and this is due to a transition in the Fermi/LAT, and HESS, respectively. transport regime. At the lowest gas densities, diffusion is an im- I first focuson a subset ofthe results to illustrate how gasdis- portant process at intermediate particle energies ∼1-10GeV, so tribution and density affect the γ-ray emission, and how the latter the steady-state particle population is not strongly affected by scales with global parameters. This is followed by a discussion gas density and the luminosity of gas-related components in- of the impact of changing some parameters of the model. I then creases with gas density, although slower than linearly. At the discuss the complete set of results in the context of the popula- highest gas densities, Bremsstrahlung losses (for CR leptons) tion study performed from Fermi/LAT observations, and of the and hadronic losses (for CR nuclei) seriously compete with dif- far-infrared - radio correlation. fusion, therefore the luminosity of gas-related components rises more weakly with gas density as the transport is moving towards a loss-dominated regime. Note that the shape of this flattening 3.1. Gamma-ray luminosities and spectra probably depends on the diffusion conditions (a higher diffusion To illustrate the effects of the gas distribution and density on γ- coefficient at low energies would imply a flattening occurring ray emission in a star-forming galaxy, I restrict the discussion to later, at a higher density). the case of a disk galaxy with Rmax =10kpc for the four molec- Band III (100GeV-100TeV): The emission is dominated by ular gas profiles and seven molecular gas densities introduced pion decay, except for the lowest average gas densities where it above. To allow a clear comparison, the runs were performed for is either dominated by IC or where IC is within a factor of 2- the same input CR luminosity (which is that of the Milky Way, 3. Inverse-Compton emission is here always dominated by pri- see Sect. 2.4). The plots of Fig. 1 show the evolution of the γ- maries, and this contribution is constant with gas density because ray luminosity in three different bands for this particular setup, the steady-state primary CRe population is loss-dominated. In while the plots of Fig. 2 show the distribution of the emission contrast, the contribution from secondaries will rise with gas in terms of physical processes and emitting particles (primary density, following the increase of the production rate of secon- electrons versus secondary electrons and positrons, hereafter pri- daries from high-energy CRp. But because secondaries have a maries and secondaries for short). softer injection spectrum, and because the parent CRp popula- These plotsdonotonlyshowthe effect of increasing gas den- tion gradually suffers from losses, secondaries never dominate sity, and correspondingly of the enhanced energy losses, but also primaries at these high energies, and the increase of the total IC include an effect of the distribution of CR sources. Indeed, from emission remains moderate, only by a factor of 2-3. In contrast, the smaller core to the larger ring, and with increasing density, pion decay emission increases faster with gas density. the ratio of star formation occurring in dense molecular regions is higher than that distributed in the lower-density atomic disk. This is why for a same increase in molecular gas density, the
6 Martin P.: Interstellar gamma-ray emission from star-forming galaxies
39 10 Pion decay Inverse Compton from primaries 1.2 Inverse Compton from secondaries Bremsstrahlung from primaries Bremsstrahlung from secondaries
1.0 ) (erg/s)
0.8
1038
0.6
0.4 Luminosity in 100.0keV - 10.0MeV 10.0MeV - 100.0keV in Luminosity (erg/s 10.0MeV - 100.0keV in Luminosity
0.2 Total 1037 Bremsstrahlung Inverse Compton Pion decay 0.0 Core 500pc Core 1000pc Ring 2-3kpc Ring 4-6kpc Core 500pc Core 1000pc Ring 2-3kpc Ring 4-6kpc
Pion decay 1040 Inverse Compton from primaries 1.2 Inverse Compton from secondaries Bremsstrahlung from primaries Bremsstrahlung from secondaries
1.0 (erg/s) (erg/s)
0.8
1039
0.6
0.4 Luminosity in 100.0MeV - 100.0GeV 100.0GeV - 100.0MeV in Luminosity 100.0GeV - 100.0MeV in Luminosity
0.2 1038 Total Bremsstrahlung Inverse Compton Pion decay 0.0 Core 500pc Core 1000pc Ring 2-3kpc Ring 4-6kpc Core 500pc Core 1000pc Ring 2-3kpc Ring 4-6kpc
Pion decay Inverse Compton from primaries 39 1.2 10 Inverse Compton from secondaries Bremsstrahlung from primaries Bremsstrahlung from secondaries
1.0 (erg/s) (erg/s)
0.8
1038
0.6
0.4
Luminosity in 100.0GeV - 100.0TeV 100.0TeV - 100.0GeV in Luminosity 1037 100.0TeV - 100.0GeV in Luminosity
0.2 Total Bremsstrahlung Inverse Compton Pion decay 0.0 Core 500pc Core 1000pc Ring 2-3kpc Ring 4-6kpc Core 500pc Core 1000pc Ring 2-3kpc Ring 4-6kpc
Fig. 1. Luminosity in 3 γ-ray bands for a star-forming disk Fig. 2. Distribution of the luminosity in 3 γ-ray bands, in terms galaxy with (Rmax, zmax)=(10kpc, 2kpc), 4 profiles of molec- of physical processes and emitting particles, for a star-forming ular gas (given in abscissa), and 7 molecular gas densities for disk galaxy with (Rmax, zmax)=(10kpc, 2kpc), 4 profiles of −3 each profile (nH2 =5, 10, 20, 50, 100, 200, 500H2 cm ). To il- molecular gas (given in abscissa), and 7 molecular gas densities −3 lustrate the effects of the gas distribution, the γ-ray luminosities for each profile (nH2 =5, 10, 20, 50, 100, 200, 500H2 cm ). are given for the same input cosmic-ray luminosity.
7 Martin P.: Interstellar gamma-ray emission from star-forming galaxies
Total Pion decay