search

rwesson / NEAT

Empirical analysis of ionised nebulae, with uncertainty propagation
GNU General Public License v3.0
9 stars 5 forks source link

Recombination line temperature diagnostics #18

Closed rwesson closed 9 years ago

rwesson commented 12 years ago

The code should calculate the Balmer jump temperature. To do this we would have to allow the user to specify continuum fluxes either side of the jump. Then we can just code in equation 3 from Liu et al. 2001:

http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?bibcode=2001MNRAS.327..141L&db_key=AST&page_ind=14&data_type=GIF&type=SCREEN_VIEW&classic=YES

It would also be very simple to make it calculate temperatures from Helium line ratios e.g. 6678/4471, 5876/4471.

pscicluna commented 12 years ago

Thought I should make it look like we're doing something, so I added a short routine that calculates the Balmer jump temperature. The wavelengths for continuum fluxes have been added to the complete line list, but this seems like a slightly opaque way of dealing with it. Will try to think of a better way.

rwesson commented 12 years ago

Nice work. This was exactly how I'd have implemented it but if you think of any better way that would be great. I fixed the dereddening routine to deredden the whole Balmer series array instead of just the first 4 lines, so that the correct dereddened H11 flux is used for T(BJ).

pscicluna commented 12 years ago

Great, thanks for spotting the bugs.

rwesson commented 9 years ago

Just added routines to get temperatures and densities from OII line ratios, now working on getting helium temperatures. 8754d3b4beb3e392e0cd7bd512771998b0e8a542

rwesson commented 9 years ago

Helium temperatures included. Both OII and HeI calculations can be significantly improved but what I've added basically works so I'll close the issue. 16fa1a1be075d4bfa398946229f80ce956ff60cb