View unanswered posts | View active topics
|
Page 1 of 1
|
[ 8 posts ] |
|
Author |
Message |
strong
|
Posted: Sat Jan 13, 2007 11:34 am |
|
Joined: Wed May 17, 2006 7:02 am Posts: 285 Location: MPE Garching
|
I recently noticed that electron inverse Compton energy losses are treated in the Thompson limit using the total ISRF energy density.
In my original development of the code around 2000 this was done as a provisional solution with the idea to upgrade this but somehow it got forgotten.
It means the energy losses at high energies (TeV) on optical radiation fields are incorrect as it enters the Klein-Nishina regime.
The correct way to do it is to integrate over the ISRF using the KN formula,
just as is done for gamma-ray production.
Also note that only continuous losses are handled by galprop,
while at high energies the losses become catastropic. This will require
a more fundamental development of the propagation code.
_________________ Andy Strong, MPE
|
|
Top |
|
|
tporter
|
Posted: Tue Nov 20, 2007 11:46 am |
|
Joined: Wed May 17, 2006 2:35 pm Posts: 23 Location: Santa Cruz, California
|
I went and checked this. The KN losses are treated correctly (in create_transport_arrays.cc when electrons/positrons are created, e_KN_loss is called to generate the IC losses for all points in the spatial grid).
_________________ Dr. Troy A. Porter Hansen Experimental Physics Laboratory 452 Lomita Mall Stanford University Stanford, CA 94305
|
|
Top |
|
|
strong
|
Posted: Tue Nov 20, 2007 12:18 pm |
|
Joined: Wed May 17, 2006 7:02 am Posts: 285 Location: MPE Garching
|
good to know.
In fact it was done by Igor in October 2006
(the other method is still there but only to test analytical cases)
But catastrophic losses are anyway not handled correcly,
a problem at TeV.
_________________ Andy Strong, MPE
|
|
Top |
|
|
tporter
|
Posted: Fri Jan 18, 2008 8:59 am |
|
Joined: Wed May 17, 2006 2:35 pm Posts: 23 Location: Santa Cruz, California
|
Right. To do it properly we would probably need to use a Monte Carlo solution for the transport equations. We could put this in but it would require some redesign (and not just a small amount of it, either).
_________________ Dr. Troy A. Porter Hansen Experimental Physics Laboratory 452 Lomita Mall Stanford University Stanford, CA 94305
|
|
Top |
|
|
imos
|
Posted: Sat Jan 19, 2008 12:46 am |
|
Joined: Tue Nov 29, 2005 6:07 pm Posts: 34 Location: Stanford
|
The effect of treating this properly is negligible, cf. secondary protons vs primary protons. I did the proper treatment of the catastrophic losses, here is why "secondary protons."
_________________ Igor Moskalenko
Stanford University
|
|
Top |
|
|
Beischer
|
Posted: Fri Nov 14, 2008 3:07 am |
|
Joined: Thu Apr 17, 2008 8:26 am Posts: 4 Location: Aachen
|
Examining the code, I found that right now the energy losses for electrons are treated in the following way:
First of all the method electron_loss(...) gets called, in which the inverse compton losses in the Thompson limit are included. Later in the code the e_KN_losses(...) function calculates the Klein-Nishina losses. Is it correct to have both methods in the code? Wouldn't energy losses due to inverse compton interactions be considered twice now?
Also, since I am interested in doing calculations with TeV electrons to study the effect of local sources on the spectrum, could you please elaborate on the problem with catastrophic losses?
Cheers Bastian Beischer
|
|
Top |
|
|
imos
|
Posted: Fri Nov 14, 2008 5:32 pm |
|
Joined: Tue Nov 29, 2005 6:07 pm Posts: 34 Location: Stanford
|
The IC losses are NOT taken into account twice. Before calling the routine electron_loss, the energy density of photon field is set to zero: double uevcm3=0.
The effect of catastrophic losses is negligible due to the very steep spectrum of VHE electrons.
_________________ Igor Moskalenko
Stanford University
|
|
Top |
|
|
strong
|
Posted: Wed Nov 19, 2008 7:54 am |
|
Joined: Wed May 17, 2006 7:02 am Posts: 285 Location: MPE Garching
|
But the steep spectrum is itself the result of the losses. Maybe it's less steep if they are catastrophic. Anyway eventually it should be done correctly if we want to do TeV and above. Meanwhile a simple analytical comparison of the approaches would shed light on the matter.
|
|
Top |
|
|
|
Page 1 of 1
|
[ 8 posts ] |
|
Who is online |
Users browsing this forum: No registered users and 2 guests |
|
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
|