search

agnwinds / python

This is the repository for Python, the radiative transfer code used to winds in AGN and other syatems
GNU General Public License v3.0
24 stars 25 forks source link

Improve dielectronic recombination #186

Open jhmatthews opened 8 years ago

jhmatthews commented 8 years ago

This is part of trying to keep track of outstanding problems before they come back to bite us, as discussed in the last telecon.

Currently dielectronic recombination is turned off in dev. This is because the current method adds kT to the cooling estimator everytime you get a dielectronic recomb, when in reality the cross section for DR is energy dependent- i.e. cooler electrons would preferentially recombine.

screenshot 2015-11-06 17 06 14

Cox and Tucker's (1969) calculation of radiative cooling of an optically thin plasma suggests that dielectronic may be important on the ~20% level -- this implies we either need it in or need to understand in what regimes it is important.

kslong commented 4 years ago

This still needs to be resolved.

Higginbottom commented 2 years ago

As part of the effort to properly incorporate Auger ionization, we now import inner shell PI rates. These are now used, via the Milne relation, so calculate dielectric recombination cooling rates. However, there is still no mechanism top account for the photons which will be produced. The sn_1d.pf example file exhibits very high DR cooling, and ends up very cold. It could be that this issue is to blame