GALPROP
https://galprop.stanford.edu/forum/

Units of DM antiproton routines
https://galprop.stanford.edu/forum/viewtopic.php?f=13&t=208
Page 1 of 1

Author:  erccarls [ Wed May 27, 2015 5:13 pm ]
Post subject:  Units of DM antiproton routines

Hi All,

I am trying to calculate propagation greens functions for antiprotons from dark matter using galprop. I have implemented a simple NFW profile and am using the gen_DM_Source.cc with the gen_DM_source routine containing the line:

particle.secondary_source_function.d2[ir][iz].s[ip]+=pow(DM_profile_av(galaxy.r[ir], galaxy.z[iz], galaxy.dr, galaxy.dz, dzz),2);

if we are in the energy bin of interest -- i.e. I inject antiprotons proportional to the NFW squared (in units GeV^2/cm^6) into a single energy bin and look at the resulting antiproton output, then repeat the galprop run for each energy bin.

I am trying to understand the resulting output units after propagation. The shape looks approximately correct if I compare to semi-analytic results for propagation functions (e.g. fig 9 of http://arxiv.org/pdf/1012.4515v4.pdf), but the normalization doesn't match. assuming a local rho_dm of 0.3 GeV/cm^3 I obtain a value at 1 GeV of around 10^21.

Any help would be much appreciated. Picture of resulting runs is also attached where one can see the Green functions for one of the models.

Attachments:
green_func.png
green_func.png [ 79.49 KiB | Viewed 29847 times ]

Page 1 of 1 All times are UTC - 8 hours [ DST ]
Powered by phpBB® Forum Software © phpBB Group
https://www.phpbb.com/