How To Use Thermal Dust Continuum Emission To Measure The Physical Properties Of Dusty Astrophysical Objects

TL;DR

This work provides a comprehensive, derivation-based framework for interpreting thermal dust continuum emission across astrophysical environments. By starting from the one-dimensional radiative transfer equation and progressing through optically thin and thick limits, it links observed fluxes to dust mass, column density, and luminosity while incorporating redshift, lensing, and CMB corrections. It also surveys dust opacity models, two common SED fitting approaches, and the use of radiative transfer codes to compute self-consistent dust temperature distributions. The resulting toolkit enables robust estimation of physical properties of dusty objects—from molecular clouds to galaxies—while highlighting degeneracies and the need for multi-wavelength, radiative-transfer-enabled analyses. The work emphasizes practical strategies for translating flux measurements into meaningful astrophysical quantities with well-defined assumptions and limits.

Abstract

Dust grains in the interstellar medium interact with photons across the electromagnetic spectrum. They are generally photon energy converters, absorbing short wavelength radiation and emitting long wavelength radiation. Sixty years ago in 1965, thermal emission from dust grains in the interstellar medium was discovered. This tutorial is a summary of the physics of thermal dust continuum emission and how to use observations of the intensity and flux density of dusty objects to calculate physical properties such as mass, column density, luminosity, dust temperature, and dust opacity spectral index. Equations are derived, when feasible, from first principles with all limits and assumptions explicitly stated. Properties of dust opacities appropriate for different astrophysical environments (e.g. diffuse ISM, dense cores, protoplanetary disks) are discussed and tabulated for the wavelengths of past, current, and future bolometer cameras. Corrections for observations at high redshift as well as the effects of telescope measurement limitations are derived. We also update the calculation of the mean molecular weight in different ISM environments and find that it is 1.404 per H atom, 2.809 per H2 molecule, and 2.351 per gas particle assuming protosolar metallicity and the latest values of the ISM gas phase abundances of metals.

Figures (21)

• Figure 1: LEFT: Optical color image of dust absorption and scattering in the Taurus Molecular Cloud. The optical image is a 21 hour exposure taken by Adam Block (Steward Observatory/University of Arizona) using the Pomenis Astrograph (Takahashi E-180 - Epsilon f2.8 180mm ED) with an APOGEE ALTA F9000 camera and equal integration time with R, G, B filters. RIGHT: Herschel Space Observatory image of thermal dust emission in the Taurus Molecular Cloud. The colors correspond to 160 $\mu$m (blue), 250 $\mu$m (green), 350 $\mu$m (split between green and red), and 500 $\mu$m (red). Regions with clear dust absorption in the optical image glow with thermal dust emission at far-infrared and submillimeter wavelengths. Credit: ESA/Herschel/NASA/JPL-Caltech, CC BY-SA 3.0 IGO; Acknowledgement: R. Hurt (JPL-Caltech). Note that some portions covered by the optical image were not imaged by the Herschel Space Observatory.
• Figure 2: LEFT: The Spectral Energy Distribution (SED) of a template Ultra Luminous InfraRed Galaxy (ULIRG) is dominated by thermal dust emission at far-infrared and submillimeter wavelengths (shaded yellow region from $\sim 20$$\mu$m to $3$ mm). The Arp 220 template (red curve) is from 2007ApJ...660..167D. RIGHT: Hubble Space Telescope image of Arp 220. The image is a combination of filters F435W (B) and F814W (I). Credit: NASA, ESA, the Hubble Heritage (STScI/AURA)-ESA/Hubble Collaboration, and A. Evans (University of Virginia, Charlottesville/NRAO/Stony Brook University)
• Figure 3: Geometry of an emitting dusty object with two lines-of-sight labeled A (observed line-of-sight) and C (line-of-sight through object center). The observer is located to the left and is looking toward the dusty object. A differential volume, $dV$, is located a position $\vec{r}$ from the center of the object. The $z$-direction is measured along the observed line-of-sight through the cloud that passes through $dV$ (labeled A). The coordinates $(x^{\prime},y^{\prime})$ are perpendicular coordinates in the plane of the sky of the observer ($x^{\prime}$ points out of the page). The impact parameter, $b$, is the shortest distance between the observed line-of-sight (A) and the cloud center (C). For an object located at a distance, $D$, that is very far away from the observer compared to the size of the object, the angle between lines-of-sight A and C, the angular impact parameter, is $\theta = b/D$ (and lines-of-sight A and C are very close to parallel). This Figure is drawn in the limit $D \gg \rm{max}\{| \vec{r}| \}$.
• Figure 4: Geometry of the optical depth coordinates in the radiative transfer problem. The arrow towards the observer along line-of-sight A indicates the direction of integration in the radiative transfer problem. $\tau_{\nu}$ is the total optical depth along A through the object. If $dV$ emits radiation into the line-of-sight A, then this radiation will see an optical depth of $\tau_{\nu} - \tau_{\nu}^{\prime}$ as it passes out of the object.
• Figure 5: Examples of thermal dust emission at $250$$\mu$m (left) and dust absorption against a $24$$\mu$m continuum (right) toward the cold ($T_d \sim 10$ K) starless core, L1689B, in the Ophiuchus Molecular Cloud. The $250$$\mu$m Herschel Space Observatory and the $24$$\mu$m Spitzer Space Telescope images were downloaded from the NASA/IPAC Infrared Science Archive.