B. T. Draine

Discrete-Dipole Approximation For Scattering Calculations

The discrete-dipole approximation (DDA) for scattering calculations, including the relationship between the DDA and other methods, is reviewed. Computational considerations, i.e., the use of complex-conjugate gradient algorithms and fast-Fourier-transform methods, are discussed. We test the accuracy of the DDA by using the DDA to compute scattering and absorption by isolated, homogeneous spheres as well as by targets consisting of two contiguous spheres. It is shown that, for dielectric materials (|m| ≲ 2), the DDA permits calculations of scattering and absorption that are accurate to within a few percent.

Dust Grain-Size Distributions and Extinction in the Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud

We construct size distributions for carbonaceous and silicate grain populations in different regions of the Milky Way, LMC, and SMC. The size distributions include sufficient very small carbonaceous grains (including polycyclic aromatic hydrocarbon molecules) to account for the observed infrared and microwave emission from the diffuse interstellar medium. Our distributions reproduce the observed extinction of starlight, which varies depending on the interstellar environment through which the light travels. As shown by Cardelli, Clayton, and Mathis in 1989, these variations can be roughly parameterized by the ratio of visual extinction to reddening, RV. We adopt a fairly simple functional form for the size distribution, characterized by several parameters. We tabulate these parameters for various combinations of values for RV and bC, the C abundance in very small grains. We also find size distributions for the line of sight to HD 210121 and for sight lines in the LMC and SMC. For several size distributions, we evaluate the albedo and scattering asymmetry parameter and present model extinction curves extending beyond the Lyman limit.We construct size distributions for carbonaceous and silicate grain populations in different regions of the Milky Way, LMC, and SMC. The size distributions include sufficient very small carbonaceous grains (including polycyclic aromatic hydrocarbon molecules) to account for the observed infrared and microwave emission from the diffuse interstellar medium. Our distributions reproduce the observed extinction of starlight, which varies depending upon the interstellar environment through which the light travels. As shown by Cardelli, Clayton & Mathis in 1989, these variations can be roughly parameterized by the ratio of visual extinction to reddening, R_V. We adopt a fairly simple functional form for the size distribution, characterized by several parameters. We tabulate these parameters for various combinations of values for R_V and b_C, the C abundance in very small grains. We also find size distributions for the line of sight to HD 210121, and for sightlines in the LMC and SMC. For several size distributions, we evaluate the albedo and scattering asymmetry parameter, and present model extinction curves extending beyond the Lyman limit.

Interstellar Dust Grains

▪ Abstract This review surveys the observed properties of interstellar dust grains: the wavelength-dependent extinction of starlight, including absorption features, from UV to infrared; optical luminescence; infrared emission; microwave emission; optical, UV, and X-ray scattering by dust; and polarization of starlight and of infrared emission. The relationship between presolar grains in meteorites and the interstellar grain population is discussed. Candidate grain materials and abundance constraints are considered. A dust model consisting of amorphous silicate grains, graphite grains, and polycyclic aromatic hydrocarbons is compared with observed emission and scattering. Some issues concerning evolution of interstellar dust are discussed.

The discrete-dipole approximation and its application to interstellar graphite grains

The discrete dipole approximation (DDA), a flexible method for computing scattering of radiation by particles of arbitrary shape, is extended to incorporate the effects of radiative reaction and to allow for possible anisotropy of the dielectric tensor of the material. Formulas are given for the evaluation of extinction, absorption, scattering, and polarization cross sections. A simple numerical algorithm based on the method of conjugate gradients is found to provide an efficient and robust method for obtaining accurate solutions to the scattering problem. The method works well for absorptive, as well as dielectric, grain materials. Two validity criteria for the DDA are presented. The DDA is then used to compute extinction cross sections for spherical graphite grains and to calculate extinction cross sections for nonspherical graphite grains with three different geometries. It is concluded that the interstellar 2175 A extinction feature could be produced by small graphite grains which should have aspect ratios not far from unity. 35 references.

Infrared Emission from Interstellar Dust. IV. The Silicate-Graphite-PAH Model in the Post-Spitzer Era

IR emission spectra are calculated for dust heated by starlight, for mixtures of amorphous silicate and graphitic grains, including varying amounts of PAH particles. The models are constrained to reproduce the average Milky Way extinction curve. The calculations include the effects of single-photon heating. Updated IR absorption properties for the PAHs are presented that are consistent with observed emission spectra, including those newly obtained by Spitzer. We find a size distribution for the PAHs giving emission band ratios consistent with the observed spectra of the Milky Way and other galaxies. Emission spectra are presented for a wide range of starlight intensities. We calculate how the efficiency of emission into different IR bands depends on PAH size; the strong 7.7 μm emission feature is produced mainly by PAH particles containing Umin. We present graphical procedures using Spitzer IRAC and MIPS photometry to estimate the parameters qPAH, Umin, and γ, the fraction fPDR of the dust luminosity coming from photodissociation regions with U > 100, and the total dust mass Mdust.

Infrared Emission from Interstellar Dust. II. The Diffuse Interstellar Medium

We present a quantitative model for the infrared emission from dust in the diffuse interstellar medium. The model consists of a mixture of amorphous silicate grains and carbonaceous grains, each with a wide size distribution ranging from molecules containing tens of atoms to large grains 1 μm in diameter. We assume that the carbonaceous grains have properties like polycyclic aromatic hydrocarbons (PAHs) at very small sizes and graphitic properties for radii a 50 A. On the basis of recent laboratory studies and guided by astronomical observations, we propose astronomical absorption cross sections for use in modeling neutral and ionized PAHs from the far-ultraviolet to the far-infrared. We also propose modifications to the far-infrared emissivity of astronomical silicate. We calculate energy distribution functions for small grains undergoing temperature spikes caused by stochastic absorption of starlight photons using realistic heat capacities and optical properties. Using a grain-size distribution consistent with the observed interstellar extinction, we are able to reproduce the near-IR to submillimeter emission spectrum of the diffuse interstellar medium, including the PAH emission features at 3.3, 6.2, 7.7, 8.6, and 11.3 μm. The model is compared with the observed emission at high Galactic latitudes as well as in the Galactic plane, as measured by the COBE/DIRBE, COBE/FIRAS, IRTS/MIRS, and IRTS/NIRS instruments. The model has 60 × 10-6 of C (relative to H) locked up in PAHs, with 45 × 10-6 of C in a component peaking at ~6 A (NC ≈ 100 carbon atoms) to account for the PAH emission features and with 15 × 10-6 of C in a component peaking at ~50 A to account for the 60 μm flux. The total infrared emission is in excellent agreement with COBE/DIRBE observations at high Galactic latitudes, just as the albedo for our grain model is in accord with observations of the diffuse Galactic light. The aromatic absorption features at 3.3 and 6.2 μm predicted by our dust model are consistent with observations. We calculate infrared emission spectra for our dust model heated by a range of starlight intensities, from 0.3 to 104 times the local interstellar radiation field, and we tabulate the intensities integrated over the SIRTF/IRAC and MIPS bands. We also provide dust opacities tabulated from the extreme-ultraviolet to submillimeter wavelengths.

SINGS: The SIRTF Nearby Galaxies Survey

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multiwavelength library of ancillary and complementary observations, including radio continuum, H i, CO, submillimeter, BVRIJHK ,H a ,P aa, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.

The calibration of mid-infrared star formation rate indicators

With the goal of investigating the degree to which the MIR emission traces the SFR, we analyze Spitzer 8 and 24 μm data of star-forming regions in a sample of 33 nearby galaxies with available HST NICMOS images in the Paα (1.8756 μm) emission line. The galaxies are drawn from the SINGS sample and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and LIRGs are also included in the analysis. Both the stellar continuum-subtracted 8 μm emission and the 24 μm emission correlate with the extinction-corrected Paα line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed nonlinear trend of the 24 μm emission versus number of ionizing photons, including the modest deficiency of 24 μm emission in the low-metallicity regions, which results from a combination of decreasing dust opacity and dust temperature at low luminosities. Conversely, the trend of the 8 μm emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 μm emission is contributed, in larger measure than the 24 μm emission, by dust heated by nonionizing stellar populations, in addition to the ionizing ones, in agreement with previous findings. Two SFR calibrations, one using the 24 μm emission and the other using a combination of the 24 μm and Hα luminosities (Kennicutt and coworkers), are presented. No calibration is presented for the 8 μm emission because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by ongoing star formation.

THE MID-INFRARED SPECTRUM OF STAR-FORMING GALAXIES: GLOBAL PROPERTIES OF POLYCYCLIC AROMATIC HYDROCARBON EMISSION

We present a sample of low-resolution 5-38 μm Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths that indicate that most normal galaxies undergo A_V ≲ 3 mag averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20% and is suppressed at metallicities Z ≲ Z_☉/4, as well as in low-luminosity AGN environments. Strong interband PAH feature strength variations (2-5 times) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8 μm features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGNs. The AGNs may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGNs may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled and points to more than a factor of 2 uncertainty in results that assume a fixed PAH emission spectrum, for redshifts z = 0-2.5.

Structure of stationary photodissociation fronts

The structure of stationary photodissociation fronts is revisited. H_2 self- shielding is discussed, including the effects of line overlap. We find that line overlap is important for N(H_2) > 10^{20} cm^{-2}. We compute multiline UV pumping models, and compare these with simple analytic approximations for the effects of self-shielding. The overall fluorescent efficiency of the photodissociation front is obtained for different ratios of chi/n_H (where chi characterizes the intensity of the incident UV) and different dust extinction laws. The dust optical depth tau_{pdr} to the point where 50% of the H is molecular is found to be a simple function of a dimensionless quantity phi_0 depending on chi/n_H, the rate coefficient for H_2 formation on grains, and the UV dust opacity. The fluorescent efficiency of the PDR also depends primarily on phi_0 for chi 10^4K, but shows some sensitivity to the v-J distribution of newly-formed H_2. The 1-0S(1)/2-1S(1) and 2-1S(1)/6-4Q(1) intensity ratios, the ortho/para ratio, and the rotational temperature in the