Making SEDs - Overview
Making Spectral Energy Distributions (SEDs)
SEDs are energy plotted against some measure of the photon -- frequency or wavelength. The reason astronomers do this is to see how much energy is produced by the object as a function of frequency or wavelength. Now it's really going to get a little hairy! Steel your nerves and plunge onwards... it really all comes down to unit conversion.
The units that are used in Spitzer and Herschel data are Janskys. 1 Jy = erg/s/cm^2/Hz (in cgs units rather than mks units, sorry -- and just to be clear, I mean erg / (s*cm^2*Hz), it's just easier to read in the way I wrote it above). A Jansky is technically a unit of "flux density," represented by . We want to plot "energy density", so one way to do this (DON'T DO THIS TO YOUR DATA YET) is to get rid of the "per Hz", e.g., multiply the Jy by the frequency () of the bandpass center, or . You could just stop here, and plot vs. to get something that is technically a spectral energy distribution. BUT, now it starts to get a little hairy, because there are some cultural influences here. Do you know off the top of your head what the frequency of the IRAC-1 band is? I don't either, but I do know its wavelength -- 3.6 microns. Astronomers coming from the longer wavelengths (mm, radio, etc.), because they think in units of frequencies, will tend to plot up against (nu). The units of really are Janskys. BUT, astronomers coming from the shorter wavelengths (optical, etc.), because they think in units of wavelengths, will tend to plot instead , where (lambda) is the wavelength of the light. This is what we want to do here (because we have been thinking of the wavelengths of the Spitzer bands but not the frequencies). The catch here is that the units of are NOT Janskys.
, the speed of light. In order to convert the Janskys into units of , you need to take into account the differentials (ah-HA, calculus being used here!), e.g., the fact that
So you need to multiply the by to convert it into . But we are not done yet! Recalling from above, the units of are not an energy density. You need to get another factor of in there to make the units work out to be energy density: calculate to get units of ergs/s/cm^2.
SO, IN SUMMARY: Take your measurements that are in Jy. (Ensure they are in Jy! If they're in magnitudes, convert them to Jy first; see 'magnitude' discussion above.) Multiply by to get them into cgs units. Multiply these values by to get them into . Multiply them by to get them into . WATCH YOUR UNITS. NB: c = 2.997924d10 cm/sec
In other words, in order to convert our data from photometry to SEDs, we need to do the following:
- Read in the catalogs you have.
- Convert any magnitudes (and errors) into flux densities (if necessary).
- Make SED plots for individual objects, but converting numbers first into the right units:
- Create an array of the wavelengths of each measurement, keeps a copy of the version in microns, and convert to cm.
- For any real measurements, convert the flux densities (probably in microJanskys) into cgs units.
- For any real measurements, convert into by multiplying the values by the corresponding to the wavelength of each bandpass.
- For any real measurements, multiply by the lambda corresponding to the wavelength of each bandpass to get .
- For any real measurements, plot the log of the data points (in cgs units) against the log of the lambda data points (in microns, only because that makes it easier to read). Label the axes (with units)! Plot the error bars on top of the data points (also converted from uJy).
Notes on plotting
Why are we plotting vs. instead of ? Well, only because I think in wavelength, not frequency. I don't know off the top of my head the frequencies of the Spitzer bandpasses, but I do know their wavelengths.
Why are we plotting instead of ? Well, only for internal consistency. Since one axis is in wavelength units, it makes sense to have the other axis also in wavelength units.
If you have gotten this far using real data, you will find (if you have done the calculations correctly) that you have numbers that are very small, like 9.77237e-12, 1.99526e-11, etc. Any time you find yourself with these kinds of numbers, you should automatically plot in log space (or log/log space), NOT linear space. You want to actually plot log() vs. log().
The next step
IF you are highly motivated and ready to go on to the next step... It can be useful, in the course of analysis of the Spitzer or Herschel or WISE data, to pretend that the contribution from the star is a blackbody. It's not really, but it's awful close, especially in the infrared. A blackbody's flux density is given by (where T is temperature, and other constants are given below)
but of course we want to plot :
Values of these constants all in cgs units:
- h = 6.6260755d-27 erg*sec
- c = 2.997924d10 cm/sec
- k = 1.380658d-16 erg/deg
So, in summary, for the list of things to do above, add this one:
- For any real measurements, for any star with at least 2 fluxes, fit a model -- one very simple way to do that is to fit a blackbody to the energies derived from the three 2MASS and first 2 IRAC bands. There are two free parameters in this fit -- the temperature of the blackbody and an additive (in the log) offset related to the distance of the object. If we know the temperature of the star (via a spectral type) and the distance to the object, then we know the values for the temperature and the offset.