https://vmcoolwiki.ipac.caltech.edu/index.php?title=Making_SEDs_-_Overview&feed=atom&action=history Making SEDs - Overview - Revision history 2024-03-29T08:02:50Z Revision history for this page on the wiki MediaWiki 1.34.2 https://vmcoolwiki.ipac.caltech.edu/index.php?title=Making_SEDs_-_Overview&diff=13371&oldid=prev Rebull at 17:14, 31 July 2020 2020-07-31T17:14:17Z <p></p> <table class="diff diff-contentalign-left" data-mw="interface"> <col class="diff-marker" /> <col class="diff-content" /> <col class="diff-marker" /> <col class="diff-content" /> <tr class="diff-title" lang="en"> <td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td> <td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 17:14, 31 July 2020</td> </tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l38" >Line 38:</td> <td colspan="2" class="diff-lineno">Line 38:</td></tr> <tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>=The next step =</div></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>=The next step =</div></td></tr> <tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td></tr> <tr><td class='diff-marker'>−</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>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 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)</div></td><td class='diff-marker'>+</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>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 <ins class="diffchange diffchange-inline">or Herschel or WISE </ins>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)</div></td></tr> <tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>   &lt;math&gt;B_{\lambda} = \left(\frac{2hc^2/\lambda^5}{\exp(hc/\lambda kT)-1)}\right)&lt;/math&gt;      (eqn 3)</div></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>   &lt;math&gt;B_{\lambda} = \left(\frac{2hc^2/\lambda^5}{\exp(hc/\lambda kT)-1)}\right)&lt;/math&gt;      (eqn 3)</div></td></tr> <tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>but of course we want to plot &lt;math&gt;\lambda \times B_{\lambda}&lt;/math&gt;:</div></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>but of course we want to plot &lt;math&gt;\lambda \times B_{\lambda}&lt;/math&gt;:</div></td></tr> </table> Rebull https://vmcoolwiki.ipac.caltech.edu/index.php?title=Making_SEDs_-_Overview&diff=13370&oldid=prev Rebull: Created page with " =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 se..." 2020-07-31T17:13:51Z <p>Created page with &quot; =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 se...&quot;</p> <p><b>New page</b></p><div><br /> =Making Spectral Energy Distributions (SEDs)=<br /> <br /> 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.<br /> Now it's really going to get a little hairy! Steel your nerves and plunge onwards... it really all comes down to unit conversion.<br /> <br /> The units that are used in Spitzer and Herschel data are Janskys. 1 Jy = &lt;math&gt;10^{-23}&lt;/math&gt; 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 &quot;flux density,&quot; represented by &lt;math&gt;F_{\nu}&lt;/math&gt;. We want to plot &quot;energy density&quot;, so one way to do this (''DON'T DO THIS TO YOUR DATA YET'') is to <br /> get rid of the &quot;per Hz&quot;, e.g., multiply the Jy by the frequency (&lt;math&gt;\nu&lt;/math&gt;) of the bandpass center, or &lt;math&gt;\nu F_{\nu}&lt;/math&gt;. You could just stop here, and plot &lt;math&gt;\nu F_{\nu}&lt;/math&gt; vs. &lt;math&gt;\nu&lt;/math&gt; 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 &lt;math&gt;\nu F_{\nu}&lt;/math&gt; against &lt;math&gt;\nu&lt;/math&gt; (nu). The units of &lt;math&gt;F_{\nu}&lt;/math&gt; really are Janskys. BUT, astronomers coming from the shorter wavelengths (optical, etc.), because they think in units of wavelengths, will tend to plot instead &lt;math&gt;\lambda F_\lambda&lt;/math&gt; , where &lt;math&gt;\lambda&lt;/math&gt; (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 &lt;math&gt;F_{\lambda}&lt;/math&gt; are NOT Janskys.<br /> <br /> &lt;math&gt;\lambda \times \nu = c&lt;/math&gt;, the speed of light. In order to convert the Janskys into units of &lt;math&gt;F_{\lambda}&lt;/math&gt;, you need to take into account the differentials (ah-HA, calculus being used here!), e.g., the fact that <br /> &lt;math&gt;\frac{dF}{d\lambda} = \frac{dF}{d\nu} \frac{d\nu}{d\lambda}&lt;/math&gt; and &lt;math&gt;d\nu = \frac{c}{\lambda^2}d\lambda&lt;/math&gt;<br /> So you need to multiply the &lt;math&gt;F_{\nu}&lt;/math&gt; by &lt;math&gt;c/\lambda^2&lt;/math&gt; to convert it into &lt;math&gt;F_{\lambda}&lt;/math&gt;. But we are not done yet! Recalling from above, the units of &lt;math&gt;F_{\lambda}&lt;/math&gt; are not an energy density. You need to get another factor of &lt;math&gt;\lambda&lt;/math&gt; in there to make the units work out to be energy density: calculate &lt;math&gt;\lambda F_{\lambda}&lt;/math&gt; to get units of ergs/s/cm^2. <br /> <br /> '''SO, IN SUMMARY''': Take your &lt;math&gt;F_{\nu}&lt;/math&gt; 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 &lt;math&gt;10^{-23}&lt;/math&gt; to get them into cgs units. Multiply these &lt;math&gt;F_{\nu}&lt;/math&gt; values by &lt;math&gt;c/\lambda^2&lt;/math&gt; to get them into &lt;math&gt;F_{\lambda}&lt;/math&gt;. Multiply them by &lt;math&gt;\lambda&lt;/math&gt; to get them into &lt;math&gt;\lambda F_\lambda&lt;/math&gt;. WATCH YOUR UNITS.<br /> NB: c = 2.997924d10 cm/sec<br /> <br /> In other words, in order to convert our data from photometry to SEDs, we need to do the following:<br /> <br /> #Read in the catalogs you have.<br /> #Convert any magnitudes (and errors) into flux densities (if necessary).<br /> #Make SED plots for individual objects, but converting numbers first into the right units:<br /> ##Create an array of the wavelengths of each measurement, keeps a copy of the version in microns, and convert to cm.<br /> ##For any real measurements, convert the flux densities (probably in microJanskys) into cgs units.<br /> ##For any real measurements, convert &lt;math&gt;F_{\nu}&lt;/math&gt; into &lt;math&gt;F_{\lambda}&lt;/math&gt; by multiplying the &lt;math&gt;F_{\nu}&lt;/math&gt; values by the &lt;math&gt;d\nu/d\lambda&lt;/math&gt; corresponding to the wavelength of each bandpass.<br /> ##For any real measurements, multiply &lt;math&gt;F_{\lambda}&lt;/math&gt; by the lambda corresponding to the wavelength of each bandpass to get &lt;math&gt;\lambda F_{\lambda}&lt;/math&gt;.<br /> #For any real measurements, plot the log of the &lt;math&gt;\lambda F_{\lambda}&lt;/math&gt; 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).<br /> <br /> <br /> <br /> =Notes on plotting=<br /> <br /> Why are we plotting &lt;math&gt;\lambda F_{\lambda}&lt;/math&gt; vs. &lt;math&gt;\lambda&lt;/math&gt; instead of &lt;math&gt;\nu&lt;/math&gt;? 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. <br /> <br /> Why are we plotting &lt;math&gt;\lambda F_{\lambda}&lt;/math&gt; instead of &lt;math&gt;\nu F_{\nu}&lt;/math&gt;? Well, only for internal consistency. Since one axis is in wavelength units, it makes sense to have the other axis also in wavelength units.<br /> <br /> 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(&lt;math&gt;\lambda F_{\lambda}&lt;/math&gt;) vs. log(&lt;math&gt;\lambda&lt;/math&gt;).<br /> <br /> =The next step =<br /> <br /> 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 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)<br /> &lt;math&gt;B_{\lambda} = \left(\frac{2hc^2/\lambda^5}{\exp(hc/\lambda kT)-1)}\right)&lt;/math&gt; (eqn 3)<br /> but of course we want to plot &lt;math&gt;\lambda \times B_{\lambda}&lt;/math&gt;:<br /> &lt;math&gt;\lambda B_{\lambda} = \left(\frac{2hc^2/\lambda^4}{\exp(hc/\lambda kT)-1)}\right)&lt;/math&gt; (eqn 4)<br /> Values of these constants all in cgs units:<br /> *h = 6.6260755d-27 erg*sec<br /> *c = 2.997924d10 cm/sec<br /> *k = 1.380658d-16 erg/deg<br /> <br /> So, in summary, for the list of things to do above, add this one:<br /> #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.</div> Rebull