Difference between revisions of "Units"

From CoolWiki
Jump to navigationJump to search
 
 
(120 intermediate revisions by 5 users not shown)
Line 1: Line 1:
=General Units=
+
=General IR astronomy units=
 
Wavelengths in infrared astronomy are commonly expressed in microns = micrometers = µm (or um if you don't have a µ).
 
Wavelengths in infrared astronomy are commonly expressed in microns = micrometers = µm (or um if you don't have a µ).
  
Line 7: Line 7:
 
*30 µm to ~350 µm = Far-infrared (~salt grain)
 
*30 µm to ~350 µm = Far-infrared (~salt grain)
  
Brightnesses or fluxes are most likely to be given in Janskys (Jy) or mJy (milli Jy) or µJy (micro Jy). 1 Jansky = 10^(-26) Watts/m^(2)/Hz. Jy can be converted to magnitudes which are rarely used in the mid- or far-infrared.
+
Brightnesses or fluxes are most likely to be given in Janskys (Jy) or mJy (milli Jy) or µJy (micro Jy). 1 Jansky = <math>10^{-26}</math> Watts/m^2/Hz.  
  
=Units of Spitzer Images=
+
Jy can be converted to magnitudes which have historically been relatively rarely used in the mid- or far-infrared.
  
Optical data with which you are familiar may be in counts or photons, or possibly (like Hubble data) calibrated to be energies. These are all ''flux units''. Spitzer data comes in ''flux (density) per unit area'' instead, MegaJanskys per steradian (MJy/sr). 1 MJy = 10^(6) Jy, and a sr is a solid angle.
+
Because the unit is named for [http://en.wikipedia.org/wiki/Karl_Jansky Karl Jansky], the plural of the unit is really Janskys, not Janskies.
  
If you've done photometry before, and expect to do it exactly the same way again here, '''it won't work''', because '''this matters'''.
+
=Fluxes and flux densities=
  
1 square arcsec is 2.3504x10^(-11) sr. (1 degree = 60 arcmin = 3600 arcsec.)
+
Astronomically, it can be important to understand the difference between luminosity, flux, and flux density.
  
If you want to convert the image from MJy/sr to uJy/square arcsec, multiply the image by 23.5045. The units of this number are (uJy/arcsec)/(MJy/sr).
+
Colloquially, flux means the rate of something through something else, such as water through a pipe, or traffic on a highway.  In physics and astronomy, it means basically the same thing, but it has a very specific meaning.
  
If you want to take a Spitzer image and use your previous routines on it, the most efficient way to do this is probably to take the image in MJy/sr and multiply out the "per sr" part of it so that it is instead in MJy/px. The subtlety in this step is that each Spitzer array has slightly different pixel sizes, and the mosaics that we create have different sizes yet again from the original images. The information on the pixel sizes are in the FITS header of each image.
+
''Flux'' is a measurement of ''energy per unit area per unit time.''  Using a traffic analogy, this would be the number of cars per lane per second that pass under a bridge on a highway (or grams of water through the cross-sectional area of a pipe per second). In measuring energy from celestial objects, the units of flux are Joules per second per meter squared if you like mks (meters-kilograms-seconds) units, or ergs per second per centimeter squared if you like cgs (centimeters-grams-seconds) units.
  
The following paragraphs are a high-level summary of what to do; see [[image unit conversion cookbook|this page]] for a cookbook of the process for one mosaic.
+
''Luminosity'' is a measurement of ''energy per unit of time,'' such as Joules per second if you like mks units, or ergs per second if you like cgs units.  This would be, in our traffic analogy, the total number of cars on the highway passing under the bridge per second.  (The flux of cars is the luminosity per lane.)
  
Look in the FITS header of the mosaics for the keywords "CDELT1" and "CDELT2". These keywords are set to be the scale of the rows and columns in degrees per pixel. Using the values of these keywords, and the conversions above, you can figure out the number of square degrees per pixel, the number of square arcsec per pixel, and finally the number of steradians per pixel. Multiply the whole image in MJy/sr by the number of sr/px to get MJy/px.  
+
''Flux density'' is a measurement essentially of ''energy per unit area per unit time "per photon".'' In our traffic analogy, this would be the number of RED cars per lane per second that pass under the bridge on the highway.  In this analogy, the "per photon" is seen in the red cars.  In astronomy, the "per photon" manifests itself as a "per Hz" (unit of frequency) or "per cm" (unit of wavelength). A Jansky is proportional to Watts/m^2/Hz.  Recall that Watts are energy per second.  So this is energy per second per square meter per Hertz.
  
If you are instead working with the individual BCDs, you should look for keywords "PXSCAL1" and "PXSCAL2". NOTE that these pixels ARE NOT SQUARE, and this is more important for MIPS data. From here, you now have the same information as the "CDELT1" and "CDELT2" above, so you can follow the same procedure.
+
Now, just to further confuse things, the units of Spitzer and some Herschel mosaics are not just Janskys, but Janskys per area!  To make the numbers easier, they are in MJy/sr, but they could also be in uJy/square arcsecond.  Read on for more, including definitions and scale factors!
  
=Units of Spitzer Photometry=
+
For completeness, we note here that ''magnitudes'' are proportional to ''the log of the ratio of two fluxes''.  Most magnitudes with which you are most likely familiar are tied to the magnitude of Vega, so a magnitude of 0 means that the object has the same flux as Vega. 
  
==Introduction==
+
[[Magnitudes]] has a different wiki page.
  
The photometry software that people use at the SSC, APEX, produces fluxes in microJanskys. The final bandmerged catalog you can get has listed fluxes in microJanskys, as well as magnitudes.
+
=Units of Spitzer and Herschel Images=
  
Astronomers use magnitudes in color-color or color-magnitude plots. Astronomers use a variant on fluxes in Spectral Energy Distribution (SED) plots.
+
Optical data with which you are familiar may be in counts or photons, or possibly (like Hubble data) calibrated to be energies. That, combined with the exposure time of the image, gives you ''flux units''. Spitzer data (and some Herschel data) come in ''flux (density) per unit (pixel) area'' instead, MegaJanskys per steradian (MJy/sr). 1 MJy = <math>10^{6}</math> Jy, and a sr is a solid angle.
  
==Magnitudes==
+
If you've done photometry before, and expect to do it exactly the same way again here, '''it won't work''', because '''this matters'''.
  
A magnitude is really a flux ratio. It is defined as follows, where M's are magnitudes and F's are fluxes:
+
1 square arcsec is <math>2.3504 \times 10^{-11}</math> sr. (1 degree = 60 arcmin = 3600 arcsec, and a sphere subtends 4pi steradians.
M1 - M2 = 2.5 * log (F2/F1)      (eqn 1)
 
  
The magnitude system (in the optical) was defined to be referenced to Vega. In other words, Vega is defined to be zero magnitude, and you would then define magnitudes of anything else as follows:
+
If you want to convert the image from MJy/sr to uJy/square arcsec, multiply the image by 23.5045. The units of this number are (uJy/arcsec)/(MJy/sr).
M = 2.5 * log (Fvega/F)      (eqn 2)
 
 
 
When they looked at Vega with IRAS, they discovered that it did NOT look like they expected, and in fact it has a large infrared excess! Therefore, infrared magnitudes are defined with respect to what Vega would be, if it did not have an excess.
 
 
 
We have published the zero points (e.g., the "Vega flux") for most of our bandpasses. They are (copied from various places on the web):
 
*IRAC 1 : 280.9 Jy
 
*IRAC 2 : 179.7 Jy
 
*IRAC 3 : 115.0 Jy
 
*IRAC 4 : 64.13 Jy
 
*MIPS 1 : 7.14 Jy
 
*MIPS 2 : 0.775 Jy
 
*MIPS 3 : 0.159 Jy
 
 
 
Therefore, in order to convert the uJy that apex returns into magnitudes, use the equation 2 above, substituting these so-called "zero-point fluxes" in for "Fvega." Note that the zero-point fluxes are in Janskys and the fluxes returned by APEX are in microJanskys.
 
 
 
You can find the zeropoints for 2MASS magnitudes on the web as well:
 
*J : 1594 Jy
 
*H : 1024 Jy
 
*K : 666.7 Jy
 
 
 
Note that plain magnitudes get fainter (the number gets larger) as the distance of the object increases. BUT, colors (differences in magnitudes) are ratios of fluxes, and therefore independent of distance.
 
 
 
==Spectral Energy Distributions (SEDs)==
 
 
 
SEDs are energy plotted against some measure of the photon -- frequency or wavelength.
 
 
 
1Jy = 10^(-23) erg/s/cm2/Hz (in cgs units rather than mks units, sorry) A Jansky is technically a unit of "flux density." In order to get rid of the "per Hz", you need to multiply the Jy by the frequency of the bandpass center.
 
 
 
Astronomers coming from the longer wavelengths will tend to plot up nu * F(nu) (written as a nu followed by an F with a subscript nu) against nu, where "nu" is the frequency. The units of F(nu) are Janskys.
 
 
 
Astronomers coming from the shorter wavelengths (like me) will tend to plot up lambda * F(lambda) (written as a lambda followed by an F with a subscript lambda), where "lambda" is the wavelength of the light. The units of F(lambda) are NOT Janskys.
 
 
 
lambda* nu = c, the speed of light. In order to convert the Janskys into units of F(lambda), you need to take into account the differentials, e.g., the fact that dF/dlambda = (dF/dnu)*(dnu/dlambda) and dnu = (c/lambda^2)dlambda.
 
 
 
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
 
 
 
A blackbody is given by
 
B(lambda) = (2hc^2/lambda^5)/(exp(hc/lambda*k*T)-1)      (eqn 3)
 
but of course we want to plot lambda*B(lambda):
 
lambda*B(lambda) = ((2hc^2)/lambda^4)/(exp((hc)/(lambda*k*T))-1)      (eqn 4)
 
  
In words, in order to begin to analyze these data, we need to have something that does the following:
+
If you want to take a Spitzer or Herschel image and use your previous routines on it, the most efficient way to do this is probably to take the image in MJy/sr and multiply out the "per sr" part of it so that it is instead in MJy/px. The subtlety in this step is that each Spitzer (or Herschel) detector has different pixel sizes, and the mosaics that we create have different sizes yet again from the original images.  You can make mosaics with whatever size pixels you want, so if you get Spitzer(or Herschel) mosaics from more than one astronomer, or more than one Spitzer(or Herschel) wavelength, chances are excellent that the pixels will be slightly different sizes (and for Herschel, maybe slightly different units).  The information on the pixel sizes are in the [[FITS_format|FITS]] header of each image.
  
  1. Reads in the fluxes from the files.
+
To figure out the pixel size, look in the FITS header of the mosaics for keywords corresponding to pixel size -- they are often "CDELT1" and "CDELT2", but may be "PXSCAL1" and "PXSCAL2" or something else. These keywords are set to be the scale of the rows and columns in degrees per pixel. Using the values of these keywords, and the conversions above, you can figure out the number of square degrees per pixel, the number of square arcsec per pixel, and finally the number of steradians per pixel. Multiply the whole image in MJy/sr by the number of sr/px to get MJy/px.
  2. Converts the Spitzer fluxes (and errors) into magnitudes (if necessary).
 
  3. Converts the 2MASS magnitudes (and errors) into fluxes (if necessary).
 
  4. Makes color-color and color-magnitude plots for stars in our region using magnitudes.
 
  5. Makes SED plots for individual objects, but converting numbers first into the right units:
 
        1. Creates an array of the wavelengths of each measurement, keeps a copy of the version in microns, and converts to cm.
 
        2. For any real measurements, converts the microJanskys into cgs units.
 
        3. For any real measurements, converts F(nu) into F(lambda) by multiplying the F(nu) values by the dnu/dlambda corresponding to the wavelength of each bandpass.
 
        4. For any real measurements, multiplies F(lambda) by the lambda corresponding to the wavelength of each bandpass to get lambda*F(lambda).
 
        5. For any real measurements, plots the log of the lambda*F(lambda) data points (in cgs units) against the log of the lambda data points (in microns, only because that makes it easier to read). Labels the axes (with units)! Plots the error bars on top of the data points (also converted from uJy).
 
        6. For any real measurements, for any star with at least 2 fluxes, fits 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.
 

Latest revision as of 18:42, 11 August 2020

General IR astronomy units

Wavelengths in infrared astronomy are commonly expressed in microns = micrometers = µm (or um if you don't have a µ).

  • 5000 Å =500 nm =0.5 µm =Visible light
  • ~0.9 to 5 µm =Near-infrared (~smoke particles)
  • 5 µm to ~30 µm = Mid-infrared (~hair)
  • 30 µm to ~350 µm = Far-infrared (~salt grain)

Brightnesses or fluxes are most likely to be given in Janskys (Jy) or mJy (milli Jy) or µJy (micro Jy). 1 Jansky = Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 10^{-26}} Watts/m^2/Hz.

Jy can be converted to magnitudes which have historically been relatively rarely used in the mid- or far-infrared.

Because the unit is named for Karl Jansky, the plural of the unit is really Janskys, not Janskies.

Fluxes and flux densities

Astronomically, it can be important to understand the difference between luminosity, flux, and flux density.

Colloquially, flux means the rate of something through something else, such as water through a pipe, or traffic on a highway. In physics and astronomy, it means basically the same thing, but it has a very specific meaning.

Flux is a measurement of energy per unit area per unit time. Using a traffic analogy, this would be the number of cars per lane per second that pass under a bridge on a highway (or grams of water through the cross-sectional area of a pipe per second). In measuring energy from celestial objects, the units of flux are Joules per second per meter squared if you like mks (meters-kilograms-seconds) units, or ergs per second per centimeter squared if you like cgs (centimeters-grams-seconds) units.

Luminosity is a measurement of energy per unit of time, such as Joules per second if you like mks units, or ergs per second if you like cgs units. This would be, in our traffic analogy, the total number of cars on the highway passing under the bridge per second. (The flux of cars is the luminosity per lane.)

Flux density is a measurement essentially of energy per unit area per unit time "per photon". In our traffic analogy, this would be the number of RED cars per lane per second that pass under the bridge on the highway. In this analogy, the "per photon" is seen in the red cars. In astronomy, the "per photon" manifests itself as a "per Hz" (unit of frequency) or "per cm" (unit of wavelength). A Jansky is proportional to Watts/m^2/Hz. Recall that Watts are energy per second. So this is energy per second per square meter per Hertz.

Now, just to further confuse things, the units of Spitzer and some Herschel mosaics are not just Janskys, but Janskys per area! To make the numbers easier, they are in MJy/sr, but they could also be in uJy/square arcsecond. Read on for more, including definitions and scale factors!

For completeness, we note here that magnitudes are proportional to the log of the ratio of two fluxes. Most magnitudes with which you are most likely familiar are tied to the magnitude of Vega, so a magnitude of 0 means that the object has the same flux as Vega.

Magnitudes has a different wiki page.

Units of Spitzer and Herschel Images

Optical data with which you are familiar may be in counts or photons, or possibly (like Hubble data) calibrated to be energies. That, combined with the exposure time of the image, gives you flux units. Spitzer data (and some Herschel data) come in flux (density) per unit (pixel) area instead, MegaJanskys per steradian (MJy/sr). 1 MJy = Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 10^{6}} Jy, and a sr is a solid angle.

If you've done photometry before, and expect to do it exactly the same way again here, it won't work, because this matters.

1 square arcsec is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 2.3504 \times 10^{-11}} sr. (1 degree = 60 arcmin = 3600 arcsec, and a sphere subtends 4pi steradians.

If you want to convert the image from MJy/sr to uJy/square arcsec, multiply the image by 23.5045. The units of this number are (uJy/arcsec)/(MJy/sr).

If you want to take a Spitzer or Herschel image and use your previous routines on it, the most efficient way to do this is probably to take the image in MJy/sr and multiply out the "per sr" part of it so that it is instead in MJy/px. The subtlety in this step is that each Spitzer (or Herschel) detector has different pixel sizes, and the mosaics that we create have different sizes yet again from the original images. You can make mosaics with whatever size pixels you want, so if you get Spitzer(or Herschel) mosaics from more than one astronomer, or more than one Spitzer(or Herschel) wavelength, chances are excellent that the pixels will be slightly different sizes (and for Herschel, maybe slightly different units). The information on the pixel sizes are in the FITS header of each image.

To figure out the pixel size, look in the FITS header of the mosaics for keywords corresponding to pixel size -- they are often "CDELT1" and "CDELT2", but may be "PXSCAL1" and "PXSCAL2" or something else. These keywords are set to be the scale of the rows and columns in degrees per pixel. Using the values of these keywords, and the conversions above, you can figure out the number of square degrees per pixel, the number of square arcsec per pixel, and finally the number of steradians per pixel. Multiply the whole image in MJy/sr by the number of sr/px to get MJy/px.