Surface brightness profiles: sech^2 fit in mag/arcsec^2 units?

AI Thread Summary
The discussion revolves around fitting a sech^2 profile to surface brightness values in mag/arcsec^2 for an edge-on galactic disk. The original poster faced challenges due to the inverted nature of magnitude units, which caused the profile to appear upside-down. They sought a method to convert the density profile into magnitude units while maintaining accuracy. Ultimately, they found success by converting mag/arcsec^2 values into luminosity density (L_sun/pc^2), fitting the data, and then translating the fitted curve back into magnitude units. The solution, while awkward, effectively preserved the necessary information for analysis.
mdemone
Messages
5
Reaction score
0
Hi everybody, I have a very quick question about surface brightness profiles for an edge-on galactic disk.

I have surface brightness values (V-band, not that it matters) in units of mag/arcsec^2 as usual, as a function of z, which is the distance from the galactic disk's plane. I'd like to fit a typical sech^2 profile, usually quoted as:

I(z) = I(0)*(sech(z/z0))^2

where I(0) is the peak luminosity density and z0 is some scale height. Now this would be fine for me, except I'm working in magnitude units, i.e. the peak is something like 21 mag/arcsec^2 and it drops off quickly (for a thin disk) to ~30 mag/arcsec^2. Thus when I plot the above profile on a y-axis inverted for magnitude units, I get an upside-down profile that disappears off in the positive y-direction.

So my question, as dumb as it may be, is how do I convert the above density profile into magnitude units, i.e. I need mu(z) = mu(0) * some sech^2 function.

I know I've seen this form floating around somewhere before, but I can't find it and my first few attempts have been miserable. I really would appreciate just a quick answer for this simple problem.

Thanks in advance! Please mention a source if at all possible, either textbook or article or whatever, just for my own future benefit.
 
Astronomy news on Phys.org
I solved it -- of course the easiest thing to do turned out to be converting the mag/arcsec^2 values into real luminosity density = L_sun/pc^2, then fitting the result and translating the fitted curve back into mag/arcsec^2 units.

Sort of awkward, but it seems to have worked without losing any information. I still could swear I've seen a magnitude version of the vertical disk sech^2 profile somewhere; ah well...
 
This thread is dedicated to the beauty and awesomeness of our Universe. If you feel like it, please share video clips and photos (or nice animations) of space and objects in space in this thread. Your posts, clips and photos may by all means include scientific information; that does not make it less beautiful to me (n.b. the posts must of course comply with the PF guidelines, i.e. regarding science, only mainstream science is allowed, fringe/pseudoscience is not allowed). n.b. I start this...
Asteroid, Data - 1.2% risk of an impact on December 22, 2032. The estimated diameter is 55 m and an impact would likely release an energy of 8 megatons of TNT equivalent, although these numbers have a large uncertainty - it could also be 1 or 100 megatons. Currently the object has level 3 on the Torino scale, the second-highest ever (after Apophis) and only the third object to exceed level 1. Most likely it will miss, and if it hits then most likely it'll hit an ocean and be harmless, but...
Back
Top