Infrared Astronomy: Studying Dust Around Stars

Infrared astronomy lets you see the heat glow of dust grains orbiting distant stars, revealing the raw material for planet formation. By measuring that glow, you can estimate how much dust is present and where it sits. This guide walks through the process, from choosing a telescope to interpreting the brightness you measure.

What You’ll Need

• A telescope with infrared capability – either ground-based with an infrared camera (e.g., on a 3–10 m class observatory) or space-based like JWST or Spitzer (archival data works too).
• Access to data archives – for example, the NASA/IPAC Infrared Science Archive (IRSA) or the European Space Agency’s archives, where you can download calibrated images and spectra.
• Data analysis software – Python with astropy, photutils, and matplotlib, or IRAF if you prefer a classic tool.
• Knowledge of the target star – its distance, spectral type, and approximate age (from catalogs like SIMBAD or Gaia).
• Calibration files – flat fields, dark frames, and standard star observations if you’re taking your own data.

Step-by-Step

1. Select a target star. Pick a star that is not too bright in the infrared (to avoid saturation) and not too faint (to get a good signal). Young stars with known debris disks, like Beta Pictoris or Fomalhaut, are good starting points. Check the star’s coordinates and distance from the SIMBAD database.
2. Choose the right infrared wavelength. Dust grains emit most strongly at wavelengths comparable to their size. For typical dust (micrometer-sized), observe at 10–30 microns (mid-infrared) or 70–100 microns (far-infrared) to catch the thermal emission. Near-infrared (1–5 microns) is less sensitive to cold dust but can trace hot dust close to the star.
3. Obtain the data. If using archival data, download the calibrated images from IRSA or the JWST archive. If observing yourself, take images with your infrared camera, using a standard dither pattern to remove bad pixels. For spectroscopy, take a spectrum of the star and the surrounding region.
4. Subtract the background. In your analysis software, measure the sky background in an annulus around the star (e.g., 5–10 arcseconds away) and subtract it from the image. This removes the zodiacal light and telescope thermal emission, which are strong in the infrared.
5. Perform aperture photometry. Place an aperture (typically 2–3 times the full width at half maximum of the star’s point spread function) over the star and sum the counts inside. Use a sky annulus to estimate the background level and subtract it. Record the total flux in counts per second.
6. Calibrate the flux. Convert counts to physical units (Jy or mJy) using the calibration factor from the instrument’s manual or from standard star observations. For archival data, this is often already done – check the header keywords like FLUXCON or BUNIT.
7. Measure the excess emission. Compare the observed infrared flux to the expected photospheric emission from the star. You can model the star’s spectrum using a stellar atmosphere model (e.g., from the PHOENIX or Kurucz grids) based on its effective temperature and luminosity. The excess above the photosphere is the dust emission.
8. Model the dust. Fit a blackbody or a more realistic dust model (e.g., using the dusty code or a simple modified blackbody) to the excess spectrum. The temperature of the dust tells you its distance from the star (using the star’s luminosity and the equilibrium temperature equation). The total infrared luminosity of the dust gives you the dust mass, if you assume a typical grain size and composition.
9. Check for spatial structure. If your image has high resolution (like JWST), look for asymmetries, rings, or clumps in the dust emission. These can indicate planets or recent collisions.
10. Interpret your results. Compare the dust temperature and mass to models of debris disks. A warm disk (100–200 K) likely lies at 1–10 AU, while a cold disk (30–50 K) is farther out. The dust mass, combined with the star’s age, tells you about the collisional cascade rate.

Common Mistakes

• Ignoring the stellar photosphere. The star itself emits infrared light, especially in the near-infrared. If you don’t subtract it properly, you’ll overestimate the dust emission.
• Using the wrong calibration. Infrared detectors have wavelength-dependent responses. Using a calibration from a different filter or instrument will give wrong fluxes.
• Underestimating background. The infrared sky is bright and variable. If you don’t subtract it accurately, your photometry will be off.
• Assuming a single temperature. Real dust disks have a range of temperatures. Fitting a single blackbody can give misleading results.
• Forgetting the distance. The flux you measure depends on the distance to the star. Always use the distance from Gaia or a reliable catalog.

Tips

• Use archival data from Spitzer or JWST for your first attempt – the calibration is already done, and you can focus on the analysis.
• Work in Janskys (Jy) for flux, as most infrared catalogs use this unit.
• For a quick estimate, assume the dust is a blackbody and use the formula: dust temperature T_d = 278 * (L_star / (r^2))^(1/4) K, where L_star is in solar luminosities and r is in AU.
• Plot the spectral energy distribution (SED) – flux vs. wavelength – to see the excess clearly.
• If you’re using Python, the photutils package has a aperture_photometry function that makes step 5 easy.

Troubleshooting

• No excess detected. Your star may have no dust, or the dust is too faint. Try a longer wavelength or a star known to have a disk.
• Negative flux after background subtraction. This happens if you over-subtracted the background. Reduce the sky annulus size or use a median instead of mean.
• Calibration factor not in header. Check the instrument documentation or the data archive’s pipeline description. For Spitzer, the calibration is in the BUNIT header keyword.
• Image saturated. If the star is too bright, use a shorter exposure or a smaller aperture. For archival data, look for a shallower observation.
• Confused with background galaxies. In deep images, distant galaxies can mimic dust emission. Use a color cut (e.g., [3.6]-[4.5] color) to distinguish them.

Frequently Asked Questions

Why use infrared instead of visible light?

Dust grains are only a few microns in size, so they scatter visible light poorly but emit thermal radiation in the infrared. The dust is also cooler than the star, so its peak emission falls in the infrared.

How do you know the dust is around the star and not in the foreground?

Foreground dust in our solar system (zodiacal light) is diffuse and smooth. Dust around another star appears as a point source or a small resolved disk that moves with the star. You can also check the star’s proper motion – if the dust moves with it, it’s associated.

Can you see planets with this method?

Not directly – planets are too faint. But the dust disk’s shape and gaps can indicate the presence of planets that sculpt the disk gravitationally.

What does the dust mass tell us?

Dust mass is a proxy for the amount of solid material available for planet formation. It also tells you about the collisional history – more dust means more recent collisions.

Do all stars have dust?

No. Young stars often have disks, but old stars may have cleared their dust. About 20% of main-sequence stars show infrared excess from debris disks.

How accurate is the dust mass estimate?

It’s uncertain by a factor of a few, because you assume a grain size distribution and composition. But it’s still useful for comparing stars.

Related

• COUNTIFS with Variable Table Column?
• Forza Horizon 6 Tokyo Tower Photo Challenge Location
• Tech Boom Vegas Bust: Nasdaq Gains and Las Vegas Housing Dips