Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCosmic ray feedback in the FIRE simulations: constraining cosmic ray propagation with GeV gamma ray emission
202
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present the implementation and the first results of cosmic ray (CR) feedback in the Feedback In Realistic Environments (FIRE) simulations. We investigate CR feedback in non-cosmological simulations of dwarf, sub-$Lstar$ starburst, and $Lstar$ galaxies with different propagation models, including advection, isotropic and anisotropic diffusion, and streaming along field lines with different transport coefficients. We simulate CR diffusion and streaming simultaneously in galaxies with high resolution, using a two moment method. We forward-model and compare to observations of $gamma$-ray emission from nearby and starburst galaxies. We reproduce the $gamma$-ray observations of dwarf and $Lstar$ galaxies with constant isotropic diffusion coefficient $kappa sim 3times 10^{29},{rm cm^{2},s^{-1}}$. Advection-only and streaming-only models produce order-of-magnitude too large $gamma$-ray luminosities in dwarf and $Lstar$ galaxies. We show that in models that match the $gamma$-ray observations, most CRs escape low-gas-density galaxies (e.g. dwarfs) before significant collisional losses, while starburst galaxies are CR proton calorimeters. While adiabatic losses can be significant, they occur only after CRs escape galaxies, so they are only of secondary importance for $gamma$-ray emissivities. Models where CRs are ``trapped in the star-forming disk have lower star formation efficiency, but these models are ruled out by $gamma$-ray observations. For models with constant $kappa$ that match the $gamma$-ray observations, CRs form extended halos with scale heights of several kpc to several tens of kpc.
rate research
Read More
We review numerical methods for simulations of cosmic ray (CR) propagation on galactic and larger scales. We present the development of algorithms designed for phenomenological and self-consistent models of CR propagation in kinetic description based on numerical solutions of the Fokker-Planck equation. The phenomenological models assume a stationary structure of the galactic interstellar medium and incorporate diffusion of particles in physical and momentum space together with advection, spallation, production of secondaries and various radiation mechanisms. The self-consistent propagation models of CRs include the dynamical coupling of the CR population to the thermal plasma. The CR transport equation is discretized and solved numerically together with the set of magneto-hydrodynamic (MHD) equations in various approaches treating the CR population as a separate relativistic fluid within the two-fluid approach or as a spectrally resolved population of particles evolving in physical and momentum space. The relevant processes incorporated in self-consistent models include advection, diffusion and streaming well as adiabatic compression and several radiative loss mechanisms. We discuss applications of the numerical models for the interpretation of CR data collected by various instruments. We present example models of astrophysical processes influencing galactic evolution such as galactic winds, the amplification of large-scale magnetic fields and instabilities of the interstellar medium.
The secondary-to-primary B/C ratio is widely used to study Galactic cosmic-ray propagation processes. The 2H/4He and 3He/4He ratios probe a different Z/A regime, therefore testing the `universality of propagation. We revisit the constraints on diffusion-model parameters set by the quartet (1H, 2H, 3He, 4He), using the most recent data as well as updated formulae for the inelastic and production cross-sections. The analysis relies on the USINE propagation package and a Markov Chain Monte Carlo technique to estimate the probability density functions of the parameters. Simulated data are also used to validate analysis strategies. The fragmentation of CNO cosmic rays (resp. NeMgSiFe) on the ISM during their propagation contributes to 20% (resp. 20%) of the 2H and 15% (resp. 10%) of the 3He flux at high energy. The C to Fe elements are also responsible for up to 10% of the 4He flux measured at 1 GeV/n. The analysis of 3He/4He (and to a less extent 2H/4He) data shows that the transport parameters are consistent with those from the B/C analysis: the diffusion model with delta~0.7 (diffusion slope), Vc~20 km/s (galactic wind), Va~40 km/s (reacceleration) is favoured, but the combination delta~0.2, Vc~0, and Va~80 km/s is a close second. The confidence intervals on the parameters show that the constraints set by the quartet data are competitive with those brought by the B/C data. These constraints are tighter when adding the 3He (or 2H) flux measurements, and the tightest when further adding the He flux. For the latter, the analysis of simulated and real data show an increased sensitivity to biases. Using secondary-to-primary ratio along with a loose prior on the source parameters is recommended to get the most robust constraints on the transport parameters.
The Tibet ASgamma experiment just reported their measurement of sub-PeV diffuse gamma ray emission from the Galactic disk, with the highest energy up to 957 TeV. These gamma-rays are most likely the hadronic origin by cosmic ray interaction with interstellar gas in the Galaxy. This measurement provides direct evidence to the hypothesis that the Galactic cosmic rays can be accelerated beyond PeV energies. In this work, we try to explain the sub-PeV diffuse gamma-ray spectrum within cosmic rays diffusive propagation model. We find there is a tension between the sub-PeV diffuse gamma rays and the local cosmic ray spectrum. To describe the sub-PeV diffuse gamma-ray flux, it generally requires larger local cosmic-ray flux than measurement in the knee region. We further calculate the PeV neutrino flux from the cosmic ray propagation model. Even all of these sub-PeV diffuse gamma rays originate from the propagation, the Galactic neutrinos only account for less than ~15% of observed flux, most of which are still from extragalactic sources.
201 -
Marco Padovani Laboratoire den Radioastronomie Millimetrique
, Ecole Normale Superieure et Observatoiren de Paris
2013
Cosmic-rays constitute the main ionising and heating agent in dense, starless, molecular cloud cores. We reexamine the physical quantities necessary to determine the cosmic-ray ionisation rate (especially the cosmic ray spectrum at E < 1 GeV and the ionisation cross sections), and calculate the ionisation rate as a function of the column density of molecular hydrogen. Available data support the existence of a low-energy component (below about 100 MeV) of cosmic-ray electrons or protons responsible for the ionisation of diffuse and dense clouds. We also compute the attenuation of the cosmic-ray flux rate in a cloud core taking into account magnetic focusing and magnetic mirroring, following the propagation of cosmic rays along flux tubes enclosing different amount of mass and mass-to-flux ratios. We find that mirroring always dominates over focusing, implying a reduction of the cosmic-ray ionisation rate by a factor of 3-4 depending on the position inside the core and the magnetisation of the core.
Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $gamma$-ray flux from individual clusters in our sample.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
