It is currently Wed Dec 04, 2024 10:11 pm

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 3 posts ] 
Author Message
PostPosted: Sun Aug 29, 2010 6:09 am 
Offline

Joined: Tue Aug 24, 2010 5:46 pm
Posts: 11
Hi,

I am trying to model the gamma-ray emission up to 100 TeV.

Here is my first attempt:
http://galprop.stanford.edu/wrxiv/0353/ ... uo2.tar.gz

As the WebRun Quick Plot (attached) shows, the inverse Compton spectrum is non-smooth above 1 TeV.

Why is that?
What parameters do I have to change to get a smooth spectrum?

Are my other parameters reasonable, or should I change some to get a good estimate of the emission at very high energies?
What parameters would correspond to a model where the high-energy emission is as bright as possible without violating any strong constraints, such as Fermi measurements?
Besides the spectrum I am also interested in the scale height in GLAT of the emission, i.e. in the region abs(GLAT) < 3.

Any help and pointers are much appreciated!

Christoph


Attachments:
gamma_54_03530002_1282898632.png
gamma_54_03530002_1282898632.png [ 12.5 KiB | Viewed 12000 times ]
Top
 Profile  
 
PostPosted: Tue Aug 31, 2010 8:55 pm 
Offline
Site Admin

Joined: Wed Oct 21, 2009 6:30 pm
Posts: 47
Location: Stanford U.
Christoph,

I made a similar run, setting the electron and photon energy grids to your values (Ekin_min=1e1, Ekin_max=1e8; E_gamma_min=1e4, E_gamma_max=1e8; long_min=0; long_max=360; lat_min=-5; lat_max=+5). I did not use anisotropic IC, just to save time. Results:
1) for E_gamma_factor=1.5 (like in your run) I get the same jagged spectrum of IC.
2) for E_gamma_factor=1.3 (just to make it smaller), IC came out with a smooth spectrum.
3) for E_gamma_factor=1.19 (a test), the IC spectrum is jagged again.

I think it is related to the fact that Ekin_factor=1.3 in these runs, and there must be some interpolation problem in GALPROP. Anyway, setting E_gamma_factor=Ekin_factor seems to work for now, and the other cases reveal a problem.

Please try re-running your configuration with a different Ekin_factor (the Import button in the Download Results section will copy all the parameters into a new form), and if you find the same, would you kindly file a bug report via http://galprop.stanford.edu/bugs/ ? I am asking you to do this instead of doing it myself because it will set a good example if users outside the GALPROP team participate in bug reporting :)

Andrey

_________________
Andrey Vladimirov
http://galprop.stanford.edu/


Top
 Profile  
 
PostPosted: Wed Sep 01, 2010 1:32 am 
Offline

Joined: Tue Aug 24, 2010 5:46 pm
Posts: 11
I ran using the parameters you suggested, got the same results and filed a bug report:
http://galprop.stanford.edu/bugs/show_bug.cgi?id=6

My other questions stand: if someone has suggestions about good GALPROP parameters to model 1 -- 100 TeV emission, please don't hold back. Preferably I would use a published model, so if someone can refer me to a paper that would be great!


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 3 posts ] 

All times are UTC - 8 hours [ DST ]


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB® Forum Software © phpBB Group