5), with a constraint of $X_{\mathrm{CR}}<0.1$ for $\alpha=2.3$ (close to the spectral index predicted by the simulations of @article:PinzkePfrommer:2010 around 220 GeV). Constraints on $X_{\mathrm{CR}}$ with -LAT limits (99% CL) depend only weakly on $\alpha$ because GeV-band gamma rays are produced by CRs with energies near the relativistic transition, that dominantly contribute to the CR pressure. $X_{\mathrm{CR}}$-constraints with -LAT limits are most constraining for an aperture of 0.4$^{\circ}$; despite the slightly weaker flux upper limits in comparison to the smaller radii of integration, we expect a considerably larger gamma-ray luminosity due to the increasing volume in this model. The best limit of $X_{\mathrm{CR}}< 0.012$ is achieved for $\alpha=2.3$, while the limit for $\alpha=2.1$ is only slightly worse $(X_{\mathrm{CR}}<0.017)$.

Simulation-based approach {#sec:simulation}
-------------------------

We complement the simplified analytical analysis with a more realistic and predictive approach derived from cosmological hydrodynamical simulations. We adopt the universal spectral and spatial gamma-ray model developed by @article:PinzkePfrommer:2010 to estimate the emission from decaying neutral pions which in clusters dominates over the inverse-Compton (IC) emission above 100 MeV. Given a density profile as, e.g., inferred by cosmological simulations or X-ray observations, the analytic approach models the CR distribution and the associated radiative emission processes from radio to the gamma-ray band. This formalism was derived from high-resolution simulations of clusters of galaxies that included radiative hydrodynamics, star formation and supernova feedback, and it followed the CR physics by tracing the most important injection and loss processes self-consistently while accounting for the CR pressure in the equation of motion [@article:PfrommerSpringelEnsslinJubelgas; @article:EnsslinPfrommerSpringelJubelgas:2007; @article:JubelgasSpringelEnsslinPfrommer:2008]. The results are in line with earlier numerical results on some of the overall characteristics of the CR distribution and the associated radiative emission processes [@article:DolagEnsslin:2000; @article:MiniatiRyuKangJones:2001
5), with a constraint of $ X_{\mathrm{CR}}<0.1 $ for $ \alpha=2.3 $ (close to the spectral index predict by the model of @article: PinzkePfrommer:2010 around 220 GeV). Constraints on $ X_{\mathrm{CR}}$ with -LAT limits (99%   CL) depend entirely weakly on $ \alpha$ because GeV - band gamma beam are produce by CRs with energy near the relativistic transition, that dominantly contribute to the CR atmospheric pressure. $ X_{\mathrm{CR}}$-constraints with -LAT limits are most restrain for an aperture of 0.4$^{\circ}$; despite the slightly weaker flux upper limits in comparison to the smaller radius of integration, we expect a considerably bombastic gamma - ray luminosity due to the increasing volume in this exemplar. The best limit of $ X_{\mathrm{CR } } < 0.012 $ is achieved for $ \alpha=2.3 $, while the limit for $ \alpha=2.1 $ is only slenderly worse $ (X_{\mathrm{CR}}<0.017)$. 

 Simulation - based approach { # sec: simulation } 
 ------------------------- 

 We complement the simplified analytical analysis with a more naturalistic and predictive access derived from cosmological hydrodynamical simulations. We adopt the universal spectral and spatial gamma - ray model developed by @article: PinzkePfrommer:2010 to calculate the discharge from decaying achromatic pions which in clusters dominates over the inverse - Compton (IC) emission above 100 MeV. give a density profile as, e.g., inferred by cosmologic simulations or X - ray observations, the analytic approach models the CR distribution and the consort radiative emission processes from radio to the gamma - re band. This formalism was derive from high - resolution simulations of clusters of galaxies that include radiative hydrodynamics, star formation and supernova feedback, and it followed the CR purgative by tracing the about important injection and loss processes self - consistently while accounting for the CR imperativeness in the equation of motion [ @article: PfrommerSpringelEnsslinJubelgas; @article: EnsslinPfrommerSpringelJubelgas:2007; @article: JubelgasSpringelEnsslinPfrommer:2008 ]. The consequence are in line with early numerical consequence on some of the overall characteristics of the CR distribution and the consociate radiative emission processes [ @article: DolagEnsslin:2000; @article: MiniatiRyuKangJones:2001
5), wihh a constraint of $X_{\mathvm{CR}}<0.1$ for $\alpha=2.3$ (close to the slectral kndex predicted by the simuletiobs of @article:PinzkePfrommer:2010 around 220 HeV). Consrraiits on $X_{\mathrm{CR}}$ with -LAB limjbs (99% CL) vepend only weally on $\alpva$ because GeV-tavd gamma rays are produced by CRs witr energoed near the relwtivpseic fgaksition, that dominantly contribufe to tie CR pressure. $C_{\mathrm{CR}}$-constraints with -PAT pimits are most cojstraining dor wb aperture ow 0.4$^{\circ}$; deskice the sliggtly weaker flux upper limits iv com'arison to jkw sldller radii of igtegration, wc expecd a conxiderably largcr gakma-eay luminosity due to the increasing vojume in tvia model. The best oinit ox $X_{\mdthro{XR}}< 0.012$ is arhisved flr $\elpha=2.3$, while the limit dor $\alpha=2.1$ is only skidytly worse $(X_{\mzthrm{CW}}<0.017)$.

Symulation-based approach {#sec:simulation}
-------------------------

We cojplement the simplified analytical analysis aith a mowe realistic and predictive approach derived from cosmkuogngwu hjdrodynamical simulations. We adopt the univerfzl siectral and spaticl gamma-ray modek fefgloped by @artizle:PinvieLfrommer:2010 to estimahe the gmissiin from dtcayimg neutral pions which in cousters domiuatws over the inversz-Compton (IC) zmissipn abpve 100 MeV. Given a densitv profjle as, e.g., ijferred bg cosmological sioulsthons or X-ray observations, ehe analyvic a'proach oodeks the CR distrihutiok and the associatef radnativa emission processes from radio to the gamma-ray band. Thos fogmalism wcs derlved from high-rqsolution simukations of cuusters of galaxixs that incltded radiativa hydrodynamirs, star fjrmarion and suodrnova feedbacl, and it yjlloqed the CR physics by fracing the mosc umportant injecyiov agd lpsf processes salf-cunskxtentuy whilt agcojntimg for the CR pressuse ih the equation of kobion [@artixle:PfromierSpringelEnxslinJubelgas; @artifle:EnvslmnPfrokmetSpringelJubelgas:2007; @article:JubelgzsSpringepEnfslinPfrommew:2008]. Thc refults are nn line with earlier numerical results oi some of the overall cyaracteristics of tke CR distribntion wnd the avsociated radiative wmission processex [@article:DolagEnsslin:2000; @articne:MinlatiRyuKangJones:2001
5), with a constraint of $X_{\mathrm{CR}}<0.1$ for to spectral index by the simulations Constraints $X_{\mathrm{CR}}$ with -LAT (99% CL) depend weakly on $\alpha$ because GeV-band gamma are produced by CRs with energies near the relativistic transition, that dominantly contribute the CR pressure. $X_{\mathrm{CR}}$-constraints with -LAT limits are most constraining for an aperture 0.4$^{\circ}$; the weaker upper limits in comparison to the smaller radii of integration, we expect a considerably larger gamma-ray due to the increasing volume in this model. best limit of $X_{\mathrm{CR}}< is achieved for $\alpha=2.3$, while limit $\alpha=2.1$ is slightly $(X_{\mathrm{CR}}<0.017)$. approach {#sec:simulation} ------------------------- complement the simplified analytical analysis with a more realistic and predictive approach derived from cosmological hydrodynamical simulations. adopt the and spatial model by to estimate the decaying neutral pions which in clusters inverse-Compton (IC) emission above 100 MeV. Given a profile as, inferred by cosmological simulations or X-ray the analytic approach models the CR distribution and associated radiative emission processes from radio to the gamma-ray band. This formalism was derived from of clusters of galaxies included radiative hydrodynamics, formation supernova and followed the physics by tracing the most important injection and loss processes self-consistently accounting for the CR pressure in the equation of motion @article:JubelgasSpringelEnsslinPfrommer:2008]. results are in with earlier numerical results some the overall characteristics of distribution the processes @article:MiniatiRyuKangJones:2001
5), with a constraint of $X_{\mathrm{Cr}}<0.1$ for $\alpha=2.3$ (cLose tO thE spEcTral IndeX predicted by thE SimuLations of @article:PinzkepfromMeR:2010 ArouND 220 GEV). ConStraintS On $x_{\MAthRm{cR}}$ WitH -Lat lImits (99% cL) dEpend onLy weakly on $\AlpHa$ Because GeV-baND gAmma rays arE prOduced by CRs wIth EnergiEs NeaR The reLatIvistIc tranSItion, tHat dominaNtLY contrIBute to tHE cR PresSure. $X_{\mathrm{CR}}$-consTRaINts with -LAT limiTs are mOsT CoNSTraIniNg for an apeRtUre of 0.4$^{\CIrc}$; despITe THE SliGHtly weaker fluX upper limitS In cOmpariSoN to THe smalLer raDiI Of iNtegration, wE expEct a consiDerablY Larger gAMma-ray lUminosIty Due To thE InCrEasInG VolUMe In tHIs mOdel. The bEsT lImit oF $X_{\maTHRM{cR}}< 0.012$ is AchIeveD for $\aLpha=2.3$, while the lImiT for $\ALphA=2.1$ is onLy sliGhtlY wOrse $(X_{\Mathrm{cR}}<0.017)$.

SimUlAtion-based approAch {#sEc:simulatIon}
-------------------------

we ComPlEment THe simpLifIed AnalytiCal analYSis WiTH A MoRe realistic and predIcTIVe Approach DeriveD FrOm COsmologiCaL hyDrodYNAmicaL simULaTions. We aDopt thE UnIvErsal spEcTral anD sPatIal Gamma-RAy moDel devEloped by @ArticLE:PinzkePfrommeR:2010 To estimate the EMiSSIoN From DecAying neutraL pioNS whiCh in CLuSteRS domiNates OvER tHE inverse-Compton (IC) emIsSion abOve 100 Mev. Given a densitY profile as, E.G., INferred bY cosMOlOGical simulatioNs or X-Ray observaTIons, the aNalytIc approaCh models tHE cR distriButIon And The ASSoCiated radiatiVE EmisSiOn proceSseS from raDio To tHe gAmmA-rAy band. ThiS formaliSm WaS dErIveD from HIgh-resolUtIon SiMulAtionS Of clusTers oF galAxIeS ThaT includED rADIatiVe HyDrodYnaMiCs, staR forMAtiOn and suPernova feEdbACk, anD iT fOllowed The CR physics bY tRacing the mOsT imPortanT INjection And loss processes self-conSIstentlY whIle acCounTing for thE CR PressuRe iN The equAtion oF motiOn [@ArtICLe:PfrOMMeRSpRiNgelEnssliNjUbeLgas; @aRtIcle:ensslinpfrommerSpringelJuBElgAs:2007; @article:JubeLgaSSprINGeLEnSSlINPfRoMMer:2008]. tHE results are in liNe with earlIeR NuMerical resULts On Some of tHe overaLl chaRActerisTics of the cR distribUtIon aND The Associated RadiativE emission PRocesSEs [@ArticLe:DOlagEnSsLin:2000; @ArticLe:MiniATiRYuKanGJones:2001
5), with a constraint of $ X_{\mathrm {CR}} <0. 1$fo r $\ alph a=2.3$ (closet o th e spectral index predi ctedby thes im ulati ons of@ ar t i cle :P in zke Pf r om mer:2 010  around  220 GeV).  Co ns traints on $ X _{ \mathrm{CR }}$  with -LAT l imi ts (99 %CL) depen d o nly w eaklyo n $\al pha$ beca us e  GeV-b a nd gamm a ra ys a re produced by CR s  w i th energies ne ar the  r e la t i vis tic  transitio n,  that dominan t ly c o ntr i bute to the C R pressure. $X_ {\math rm {CR } }$-con strai nt s  wi th -LAT lim itsare mostconstr a ining f o r an ap erture  of  0. 4$^{ \ ci rc }$;  d e spi t ethe sli ghtly we ak er  flux  upp e r l imit s i n co mpari son to the sm all er r a dii  of i ntegr atio n,  we e xpecta con si derably largergamm a-ray lum ino si tydu e tot he inc rea sin g volum e in th i s m od e l .  T he best limit of $ X_ { \ ma thrm{CR} }< 0.0 1 2$  i s  achieve dfor  $\a l p ha=2. 3$,w hi le the l imit f o r$\ alpha=2 .1 $ is o nl y s lig htlyw orse  $(X_{ \mathrm{ CR}}< 0 .017)$.

Simul a tion-based ap p ro a c h{ #sec :si mulation}
- ---- - ---- ---- - -- --- - ----
We c om p le m ent the simplifiedan alytic al an alysis with a  more real i s t ic and p redi c ti v e approach der ivedfrom cosmo l ogical h ydrod ynamical  simulati o n s. We ad opt  th e u niv e r sa l spectral an d spat ia l gamma -ra y model  de vel ope d b y@article: PinzkePf ro mm er :2 010  to e s timate t he  em is sio n fro m  decay ing n eutr al  p i ons  whichi nc l uste rs  d omin ate soverthei nve rse-Com pton (IC)  em i ssio nab ove 100  MeV. Given a  d ensity pro fi leas, e. g . , inferr ed by cosmological simu l ationsorX-ray  obs ervations , t he ana lyt i c appr oach m odels  t heC R  dist r i bu tio nand the as s o cia ted r ad iati ve emis sion processes fro m  ra dio to the ga mma -ray b an d.T hi s  fo rm a lis m was derived fro m high-res ol u ti on simulat i ons  o f clust ers ofgalax i es that  included  radiativ ehydr o d yna mics, star  formati on and su p ernov a  f eedba ck,  and i tfol lowed  the C R  ph ysics  by tr ac ing th e mos timportan t injection and loss pr ocesse s sel f-c onsistent lyw hil e account ingfor the CR  pr ess ure i n t h e equ atio n  o f m o tion[@ar t icle:Pfro m me rSp r i ng elEnsslinJu b e l gas ; @ar tic l e:Enss linP frommerSpringelJu b elgas:2007; @a rtic l e :Ju bel g asSp ri ngelEnsslinPfr omm er : 2 008]. Th eresults are  in line  w i th ea rliernumeri cal res u l ts on som e of  th e overall  ch ar a cterist ic so f theCR d is tribut ion an d  the a ssociated radiat ive e m i ssion pro cesse s[@artic l e:Do lagEnsslin :2000; @art icle:M inia tiRyu KangJon es :2001
5), with_a constraint_of $X_{\mathrm{CR}}<0.1$ for $\alpha=2.3$_(close to_the_spectral index_predicted_by the simulations_of @article:PinzkePfrommer:2010 around_220 GeV). Constraints on_$X_{\mathrm{CR}}$ with -LAT_limits_(99% CL) depend only weakly on $\alpha$ because GeV-band gamma rays are produced by CRs_with_energies near_the_relativistic_transition, that dominantly contribute to_the CR pressure. $X_{\mathrm{CR}}$-constraints with_-LAT limits_are most constraining for an aperture of 0.4$^{\circ}$;_despite_the slightly weaker_flux upper limits in comparison to the smaller radii_of integration, we expect a considerably_larger gamma-ray luminosity_due_to_the increasing volume in_this model. The best limit of_$X_{\mathrm{CR}}< 0.012$ is achieved for $\alpha=2.3$,_while the limit for $\alpha=2.1$ is only_slightly worse $(X_{\mathrm{CR}}<0.017)$.

Simulation-based approach {#sec:simulation}
-------------------------

We complement_the simplified analytical analysis with_a more_realistic and predictive approach derived_from cosmological hydrodynamical_simulations. We_adopt the universal_spectral and spatial gamma-ray model developed_by @article:PinzkePfrommer:2010 to_estimate the emission from decaying neutral_pions_which in clusters_dominates_over_the inverse-Compton_(IC) emission above_100_MeV. Given_a_density profile as, e.g., inferred by_cosmological_simulations or X-ray observations, the analytic approach_models the CR distribution_and_the associated radiative emission_processes from radio to the_gamma-ray band. This formalism was derived_from high-resolution_simulations of_clusters of galaxies that included radiative hydrodynamics, star formation and supernova_feedback, and it followed the CR_physics by tracing the_most important_injection_and loss processes_self-consistently_while accounting_for the CR pressure in the equation_of motion_[@article:PfrommerSpringelEnsslinJubelgas; @article:EnsslinPfrommerSpringelJubelgas:2007; @article:JubelgasSpringelEnsslinPfrommer:2008]. The results are_in line with earlier_numerical_results on some of the overall_characteristics of the CR distribution and_the associated radiative emission processes_[@article:DolagEnsslin:2000;_@article:MiniatiRyuKangJones:2001