Understanding intermittent dissipation in the interstellar medium is integral to understanding the chemical abundances and line luminosities observed therein. In the cold neutral medium, continuous heating processes due to cosmic rays and photoelectric emission from dust grains determine a unique equilibrium temperature varying with density. Excursions from this are due to intermittent processes including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^+$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We carry out simulations of MHD turbulence in molecular clouds including heating and cooling, and post-process them to study H$_2$ line emission and hot-gas chemistry, particularly the formation of CH$^+$. We explore multiple magnetic field strengths and equations of state. We use a new H$_2$ cooling function for $n_{\rm H} \leq 10^5\,{\rm cm}^{-3}$, $T\leq 5000\,{\rm K}$, and variable H$_2$ fraction. Our models produce H$_2$ emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic r.m.s. magnetic field strengths ($\approx 10$ $\mu$G) and velocity dispersions, we reproduce observed CH$^+$ abundances. Comparison of predicted dust polarization with observations by ${\it Planck}$ suggests that the mean field $\gtrsim 5\,\mu$G, so that the turbulence is sub-Alfv\'enic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH$^+$ abundance.

The sightline toward the luminous blue hypergiant Cyg OB2-12 is widely used in studying interstellar dust on account of its large extinction ($A_V \simeq 10$ mag) and the fact that this extinction appears to be dominated by dust typical of the diffuse interstellar medium. We present a new analysis of archival ISO-SWS and Spitzer IRS observations of Cyg OB2-12 using a model of the emission from the star and its stellar wind to determine the total extinction $A_\lambda$ from 2.4--37 $\mu$m. In addition to the prominent 9.7 and 18 $\mu$m silicate features, we robustly detect absorption features associated with polycyclic aromatic hydrocarbons (PAHs), including the first identification of the 7.7 $\mu$m feature in absorption. The 3.3 $\mu$m aromatic feature is found to be much broader in absorption than is typically seen in emission. The 3.4 and 6.85 $\mu$m aliphatic hydrocarbon features are observed with relative strengths consistent with observation of these features on sightlines toward the Galactic Center. We identify and characterize more than sixty spectral lines in this wavelength range, which may be useful in constraining models of the star and its stellar wind. Based on this analysis, we present an extinction curve $A_\lambda/A_{2.2 \mu m}$ that extrapolates smoothly to determinations of the mean Galactic extinction curve at shorter wavelengths and to dust opacities inferred from emission at longer wavelengths, providing a new constraint on models of interstellar dust in the mid-infrared.

We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $\alpha_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $\alpha_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $\alpha_{\mathrm{[CI](1-0)}}$ to $\alpha_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $\alpha_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas.

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.

The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample.

The X-ray regime is a largely underused resource for constraining interstellar dust grain models and improving our understanding of the physical processes that dictate how grains evolve over their lifetimes. This is mostly due to current detectors' relatively low sensitivity and high background, limiting the targets to the brightest sources. The improved sensitivity of the next generation of X-ray detectors will allow studies of much fainter sources, at much higher angular resolution, expanding our sampled sightlines in both quality and quantity.

X-ray absorption fine structure (XAFS) in the 0.2-2 keV band is a crucial component in multi-wavelength studies of dust mineralogy, size, and shape -- parameters that are necessary for interpreting astronomical observations and building physical models across all fields, from cosmology to exoplanets. Despite its importance, many fundamental questions about dust remain open. What is the origin of the dust that suffuses the interstellar medium (ISM)? Where is the missing interstellar oxygen? How does iron, predominantly produced by Type Ia supernovae, become incorporated into dust? What is the main form of carbon in the ISM, and how does it differ from carbon in stellar winds? The next generation of X-ray observatories, employing microcalorimeter technology and $R \equiv \lambda/\Delta \lambda \geq 3000$ gratings, will provide pivotal insights for these questions by measuring XAFS in absorption and scattering. However, lab measurements of mineralogical candidates for astrophysical dust, with R > 1000, are needed to fully take advantage of the coming observations.

We present maps of the dust properties in the Small and Large Magellanic Clouds (SMC, LMC) from fitting Spitzer and Herschel observations with the \citet{DL07} dust model. We derive the abundance of the small carbonaceous grain (or polycyclic aromatic hydrocarbon; PAH) component. The global PAH fraction (q_pah, the fraction of the dust mass in the form of PAHs) is smaller in the SMC (1.0$^{+0.3}_{-0.3}$%) than in the LMC (3.3$^{+1.4}_{-1.3}$%). We measure the PAH fraction in different gas phases (H II regions, ionized gas outside of H II regions, molecular gas, and diffuse neutral gas). H II regions appear as distinctive holes in the spatial distribution of the PAH fraction. In both galaxies, the PAH fraction in the diffuse neutral medium is higher than in the ionized gas, but similar to the molecular gas. Even at equal radiation field intensity, the PAH fraction is lower in the ionized gas than in the diffuse neutral gas. We investigate the PAH life-cycle as a function of metallicity between the two galaxies. The PAH fraction in the diffuse neutral medium of the LMC is similar to that of the Milky Way ($\sim4.6$%), while it is significantly lower in the SMC. Plausible explanations for the higher PAH fraction in the diffuse neutral medium of the LMC compared to the SMC include: a more effective PAH production by fragmentation of large grains at higher metallicity, and/or the growth of PAHs in molecular gas.

We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR.

Some authors have proposed that electron energy distributions in H II regions and planetary nebulae may be significantly nonthermal, and kappa-distributions have been suggested as being appropriate. Here it is demonstrated that the electron energy distribution function is extremely close to a Maxwellian up to electron kinetic energies ~13 eV in HII regions, and up to ~16eV in planetary nebulae: kappa-distributions are inappropriate. The small departures from a Maxwellian have negligible effects on line ratios. When observed line ratios in H II regions deviate from models with a single electron temperature, it must arise from spatial variations in electron temperature, rather than local deviations from a Maxwellian.

The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ISM in many post-starbursts raise important questions. Are they truly quiescent and, if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infrared spectroscopic and photometric survey of 33 SDSS-selected E+A post-starbursts, aimed at resolving these questions. We find compact, warm dust reservoirs with high PAH abundances, and total gas and dust masses significantly higher than expected from stellar recycling alone. Both PAH/TIR and dust-to-burst stellar mass ratios are seen to decrease with post-burst age, indicative of the accumulating effects of dust destruction and an incipient transition to hot, early-type ISM properties. Their infrared spectral properties are unique, with dominant PAH emission, very weak nebular lines, unusually strong H$_{2}$ rotational emission, and deep ${\rm [C\, II]}$ deficits. There is substantial scatter among SFR indicators, and both PAH and TIR luminosities provide overestimates. Even as potential upper limits, all tracers show that the SFR has typically experienced a more than two order-of-magnitude decline since the starburst, and that the SFR is considerably lower than expected given both their stellar masses and molecular gas densities. These results paint a coherent picture of systems in which star formation was, indeed, rapidly truncated, but in which the ISM was $\textit{not}$ completely expelled, and is instead supported against collapse by latent or continued injection of turbulent or mechanical heating. The resulting aging burst populations provide a "high-soft" radiation field which seemingly dominates the E+As' unusual ISM energetics.

Spinning dust grains in front of the bright Galactic synchrotron background can produce a weak absorption signal that could affect measurements of high redshift 21 cm absorption. At frequencies near 80 MHz where the EDGES experiment has reported 21\,cm absorption at $z \approx 17$, absorption could be produced by interstellar nanoparticles with radii $a \approx 50\AA$ in the cold interstellar medium at temperature $T \approx 50$ K. Atmospheric aerosols could contribute additional absorption. The strength of the absorption depends on the abundance of such grains and on their dipole moments, which are uncertain. The breadth of the absorption spectrum of spinning dust limits its possible impact on measurement of a relatively narrow 21 cm absorption feature.

Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016.

We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from ~1.5 to ~1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value ~1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples.

Although interstellar grains are known to be aspherical, their actual shapes remain poorly constrained. We assess whether three distributions of ellipsoidal shapes from the literature, employed principally for their convenient mathematical properties, are suitable for describing the shapes of interstellar grains. Randomly-selected shapes from each distribution are shown as illustrations. The often-used BHCDE distribution includes a very large fraction of extreme shapes: fully 10% of random draws have axial ratio a_3/a_1 > 19.7, and 5% have a_3/a_1 > 33. Of the three distributions considered, the CDE2 appears to be most realistic. For each of the three CDEs considered, we derive shape-averaged cross sections for extinction and polarization. Finally, we describe a method for "synthesizing" a dielectric function for an assumed shape or shape distribution if the actual absorption cross sections per grain volume are known from observations. This synthetic dielectric function predicts the wavelength dependence of polarization, which can then be compared to observations to constrain the grain shape.

The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.

The ratio of the [CII] 158$\,\mu$m emission line over the total infrared emission (TIR) is often used as a proxy for the photoelectric (PE) heating efficiency ($\epsilon_{\rm PE}$) of the far-ultraviolet (FUV) photons absorbed by dust in the interstellar medium. In the nearby galaxy M31, we measure a strong radial variation of [CII]/TIR that we rule out as being due to an intrinsic variation in $\epsilon_{\rm PE}$. [CII]/TIR fails as a proxy for $\epsilon_{\rm PE}$, because the TIR measures all dust heating, not just the contribution from FUV photons capable of ejecting electrons from dust grains. Using extensive multiwavelength coverage from the FUV to far-infrared (FIR), we infer the attenuated FUV emission ($\rm UV_{att}$), and the total attenuated flux ($\rm TOT_{att}$). We find [CII]/TIR to be strongly correlated with $\rm UV_{att}$/$\rm TOT_{att}$, indicating that, in M31 at least, one of the dominant drivers for [CII]/TIR variation is the relative hardness of the absorbed stellar radiation field. We define $\rm{ \epsilon_{PE}^{UV}}$, [CII]/$\rm{ UV_{att}}$ which should be more closely related to the actual PE efficiency, which we find to be essentially constant ($1.85 \pm 0.8 \%$) in all explored fields in M31. This suggests that part of the observed variation of [CII]/TIR in other galaxies is likely due to a change in the relative hardness of the absorbed stellar radiation field, caused by a combination of variations in the stellar population, dust opacity and galaxy metallicity, although PE efficiency may also vary across a wider range of environments.

Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $\mu$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

In light of recent observational results indicating an apparent lack of correlation between the Anomalous Microwave Emission (AME) and mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs), we assess whether rotational emission from spinning silicate and/or iron nanoparticles could account for the observed AME without violating observational constraints on interstellar abundances, ultraviolet extinction, and infrared emission. By modifying the SpDust code to compute the rotational emission from these grains, we find that nanosilicate grains could account for the entirety of the observed AME, whereas iron grains could be responsible for only a fraction, even for extreme assumptions on the amount of interstellar iron concentrated in ultrasmall iron nanoparticles. Given the added complexity of contributions from multiple grain populations to the total spinning dust emission, as well as existing uncertainties due to the poorly-constrained grain size, charge, and dipole moment distributions, we discuss generic, carrier-independent predictions of spinning dust theory and observational tests that could help identify the AME carrier(s).

Interstellar iron in the form of metallic iron nanoparticles may constitute a component of the interstellar dust. We compute the stability of iron nanoparticles to sublimation in the interstellar radiation field, finding that iron clusters can persist down to a radius of $\simeq 4.5\,$\AA, and perhaps smaller. We employ laboratory data on small iron clusters to compute the photoelectric yields as a function of grain size and the resulting grain charge distribution in various interstellar environments, finding that iron nanoparticles can acquire negative charges particularly in regions with high gas temperatures and ionization fractions. If $\gtrsim 10\%$ of the interstellar iron is in the form of ultrasmall iron clusters, the photoelectric heating rate from dust may be increased by up to tens of percent relative to dust models with only carbonaceous and silicate grains.

We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas.

Laboratory measurements are used to constrain the dielectric tensor for graphite, from microwave to X-ray frequencies. The dielectric tensor is strongly anisotropic even at X-ray energies. The discrete dipole approximation is employed for accurate calculations of absorption and scattering by single-crystal graphite spheres and spheroids. For randomly-oriented single-crystal grains, the so-called 1/3 - 2/3 approximation for calculating absorption and scattering cross sections is exact in the limit a/lambda -> 0, provides better than ~10% accuracy in the optical and UV even when a/lambda is not small, but becomes increasingly inaccurate at infrared wavelengths, with errors as large as ~40% at lambda = 10 micron. For turbostratic graphite grains, the Bruggeman and Maxwell Garnett treatments yield similar cross sections in the optical and ultraviolet, but diverge in the infrared, with predicted cross sections differing by over an order of magnitude in the far-infrared. It is argued that the Maxwell Garnett estimate is likely to be more realistic, and is recommended. The out-of-plane lattice resonance of graphite near 11.5 micron may be observable in absorption with the MIRI spectrograph on JWST. Aligned graphite grains, if present in the ISM, could produce polarized X-ray absorption and polarized X-ray scattering near the carbon K edge.

The effect of line overlap in the Lyman and Werner bands, often ignored in galactic studies of the atomic-to-molecular transition, greatly enhances molecular hydrogen self-shielding in low metallicity environments, and dominates over dust shielding for metallicities below about 10% solar. We implement that effect in cosmological hydrodynamics simulations with an empirical model, calibrated against the observational data, and provide fitting formulae for the molecular hydrogen fraction as a function of gas density on various spatial scales and in environments with varied dust abundance and interstellar radiation field. We find that line overlap, while important for detailed radiative transfer in the Lyman and Werner bands, has only a minor effect on star formation on galactic scales, which, to a much larger degree, is regulated by stellar feedback.

The attenuation of starlight by dust in galactic environments is investigated through models of radiative transfer in a spherical, clumpy ISM. Extinction properties for MW, LMC, and SMC dust types are considered. It is illustrated that the attenuation curves are primarily determined by the wavelength dependence of absorption rather than by the underlying extinction (absorption+scattering) curve. Attenuation curves consistent with the "Calzetti attenuation curve" are found by assuming the silicate-carbonaceous dust model for the MW, but with the 2175A absorption bump suppressed or absent. The discrepancy between our results and previous work that claimed the SMC-type dust to be the most likely origin of the Calzetti curve is ascribed to the difference in adopted albedos; this study uses the theoretically calculated albedos whereas the previous ones adopted empirically derived albedos from observations of reflection nebulae. It is also found that the model attenuation curves calculated with the MW dust are well represented by a modified Calzetti curve with a varying slope and UV bump strength. The strong correlation between the slope and UV bump strength, with steeper curves having stronger bumps, as found in star-forming galaxies at 0.5<z<2.0, is well reproduced by our models if the abundance of the UV bump carriers or PAHs is assumed to be 30% or 40% of that of the MW-dust. The trend is explained by radiative transfer effects which lead to shallower attenuation curves with weaker UV bumps as the ISM is more clumpy and dustier. We also argue that at least some of the IUE local starburst galaxies may have a UV bump feature in their attenuation curves, albeit much weaker than that of the MW extinction curve.

The quantization of energy levels in very nanoparticles suppresses dissipative processes that convert grain rotational kinetic energy into heat. For grains small enough to have GHz rotation rates, the suppression of dissipation can be extreme. As a result, alignment of such grains is suppressed. This applies both to alignment of the grain body with its angular momentum J, and to alignment of J with the local magnetic field B_0. If the anomalous microwave emission is rotational emission from spinning grains, it will be negligibly polarized at GHz frequencies, with P < 10^{-6} at frequencies above 10 GHz.

The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.

To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quantifying what physical processes drive the variations between different hybrid estimators. Doing so, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation, taking into account the intrinsic physical properties of galaxies. We use the CIGALE code to model the spatially-resolved FUV to FIR SED of eight nearby star-forming galaxies drawn from the KINGFISH sample. This allows us to determine their local physical properties, and in particular their UV attenuation, average SFR, average specific SFR (sSFR), and their stellar mass. We then examine how hybrid estimators depend on said properties. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 and 100 micron) and on the sSFR (in particular at 24 micron and the TIR). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of hybrid estimators based on observed UV to near-IR colours and near-IR luminosity densities per unit area. We find that they can reliably be extended to entire galaxies. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators. Naturally taking into account the variable impact of dust heated by old stellar populations, they constitute a step towards universal estimators.

We investigate the claim that the ratio {\beta} of radiation pressure force to gravitational force on a dust grain in our solar system can substantially exceed unity for some grain sizes, provided that grain porosity is high enough. For model grains consisting of random aggregates of silicate spherules, we find that the maximum value of {\beta} is almost independent of grain porosity, but for small (<0.3 {\mu}m) grains, {\beta} actually decreases with increasing porosity. We also investigate the effect of metallic iron and amorphous carbon inclusions in the dust grains and find that while these inclusions do increase the radiation pressure cross-section, {\beta} remains below unity for grains with 3 pg of silicate material. These results affect the interpretation of the grain trajectories estimated from the {\it Stardust} mission, which were modeled assuming {\beta} values exceeding one. We find that radiation pressure effects are not large enough for particles Orion and Hylabrook captured by {\it Stardust} to be of interstellar origin given their reported impact velocities. We also consider the effects of solar radiation on transverse velocities and grain spin, and show that radiation pressure introduces both transverse velocities and equatorial spin velocities of several hundred meters per second for incoming interstellar grains at 2 AU. These transverse velocities are not important for modeling trajectories, but such spin rates may result in centrifugal disruption of aggregates.

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

Interstellar abundance determinations from fits to X-ray absorption edges often rely on the incorrect assumption that scattering is insignificant and can be ignored. We show instead that scattering contributes significantly to the attenuation of X-rays for realistic dust grain size distributions and substantially modifies the spectrum near absorption edges of elements present in grains. The dust attenuation modules used in major X-ray spectral fitting programs do not take this into account. We show that the consequences of neglecting scattering on the determination of interstellar elemental abundances are modest; however, scattering (along with uncertainties in the grain size distribution) must be taken into account when near-edge extinction fine structure is used to infer dust mineralogy. We advertise the benefits and accuracy of anomalous diffraction theory for both X-ray halo analysis and near edge absorption studies. An open source Fortran suite, General Geometry Anomalous Diffraction Theory (GGADT), is presented that calculates X-ray absorption, scattering, and differential scattering cross sections for grains of arbitrary geometry and composition.

We employ an all-sky map of the anomalous microwave emission (AME) produced by component separation of the microwave sky to study correlations between the AME and Galactic dust properties. We find that while the AME is highly correlated with all tracers of dust emission, the best predictor of the AME strength is the dust radiance. Fluctuations in the AME intensity per dust radiance are uncorrelated with fluctuations in the emission from polycyclic aromatic hydrocarbons (PAHs), casting doubt on the association between AME and PAHs. The PAH abundance is strongly correlated with the dust optical depth and dust radiance, consistent with PAH destruction in low density regions. We find that the AME intensity increases with increasing radiation field strength, at variance with predictions from the spinning dust hypothesis. Finally, the temperature-dependence of the AME per dust radiance disfavors the interpretation of the AME as thermal emission. A reconsideration of other AME carriers, such as ultrasmall silicates, and other emission mechanisms, such as magnetic dipole emission, is warranted.

We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and $E_{K}\sim6-20\times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAH's with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow.

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.

We investigate the composition of interstellar grains along the line of sight toward Zeta Ophiuchi, a well-studied environment near the diffuse-dense cloud transition. A spectral decomposition analysis of the solid-state absorbers is performed using archival spectroscopic observations from the Spitzer Space Telescope and Infrared Space Observatory. We find strong evidence for the presence of sub-micron-sized amorphous silicate grains, principally comprised of olivine-like composition, with no convincing evidence of H2O ice mantles. However, tentative evidence for thick H2O ice mantles on large (a ~ 2.8 microns) grains is presented. Solid-state abundances of elemental Mg, Si, Fe, and O are inferred from our analysis and compared to standard reference abundances. We find that nearly all of elemental Mg and Si along the line of sight are present in amorphous silicate grains, while a substantial fraction of elemental Fe resides in compounds other than silicates. Moreover, we find that the total abundance of elemental O is largely inconsistent with the adopted reference abundances, indicating that as much as ~156 ppm of interstellar O is missing along the line of sight. After taking into account additional limits on the abundance of elemental O in other O-bearing solids, we conclude that any missing reservoir of elemental O must reside on large grains that are nearly opaque to infrared radiation.

In this chapter, we will outline the scientific motivation for studying Anomalous Microwave Emission (AME) with the SKA. AME is thought to be due to electric dipole radiation from small spinning dust grains, although thermal fluctuations of magnetic dust grains may also contribute. Studies of this mysterious component would shed light on the emission mechanism, which then opens up a new window onto the interstellar medium (ISM). AME is emitted mostly in the frequency range $\sim 10$--100\,GHz, and thus the SKA has the potential of measuring the low frequency side of the AME spectrum, particularly in band 5. Science targets include dense molecular clouds in the Milky Way, as well as extragalactic sources. We also discuss the possibility of detecting rotational line emission from Poly-cyclic Aromatic Hydrocarbons (PAHs), which could be the main carriers of AME. Detecting PAH lines of a given spacing would allow for a definitive identification of specific PAH species.

The [CII] 157.74 $\mu$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $\Sigma_{\rm [CII]}-\Sigma_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $\Sigma_{\rm[CII]}-\Sigma_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $\Sigma_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields responsible for heating the dust. After fitting the dust and stellar radial profiles with exponential functions, we fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified black bodies using both variable (MBBV) and fixed (MBBF) emissivity indices beta, as well as with physically motivated dust models. Results show that while most SED parameters decrease with radius, the emissivity index beta also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average beta~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with beta<= 1. We find that flatter slopes (beta<= 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust-beta degeneracy. This trend is related to variations in Umin since beta and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat beta<=1. Both these results strongly suggest that the low apparent \beta values (flat slopes) in MBBV fits are caused by temperature mixing along the line-of-sight, rather than by intrinsic variations in grain properties. Abstract truncated for arXiv.

NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the central 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems.

Because the 8 {\mu}m polycyclic aromatic hydrocarbon (PAH) emission has been found to correlate with other well-known star formation tracers, it has widely been used as a star formation rate (SFR) tracer. There are, however, studies that challenge the accuracy and reliability of the 8 {\mu}m emission as a SFR tracer. Our study, part of the Herschel M33 Extended Survey (HERM33ES) open time key program, aims at addressing this issue by analyzing the infrared emission from the nearby spiral galaxy M33 at the high spatial scale of 75 pc. Combining data from the Herschel Space Observatory and the Spitzer Space Telescope we find that the 8 {\mu}m emission is better correlated with the 250 {\mu}m emission, which traces cold interstellar gas, than with the 24 {\mu}m emission. The L(8)/L(24) ratio is highly depressed in 24 {\mu}m luminous regions, which correlate with known HII regions. We also compare our results with the dust emission models by Draine & Li (2007). We confirm that the depression of 8 {\mu}m PAH emission near star-forming regions is higher than what is predicted by models; this is possibly an effect of increased stellar radiation from young stars destroying the dust grains responsible for the 8 {\mu}m emission as already suggested by other authors. We find that the majority of the 8 {\mu}m emission is fully consistent with heating by the diffuse interstellar medium, similar to what recently determined for the dust emission in M31 by Draine at al. (2013). We also find that the fraction of 8 {\mu}m emission associated with the diffuse interstellar radiation field ranges between 60% and 80% and is 40% larger than the diffuse fraction at 24 {\mu}m.

We model the infrared to submillimeter spectral energy distribution of 11 nearby galaxies of the KINGFISH sample using Spitzer and Herschel data and compare model extrapolations at 870um (using different fitting techniques) with LABOCA 870um observations. We investigate how the differences between predictions and observations vary with model assumptions or environment. At global scales, we find that modified blackbody models using realistic cold emissivity indices (beta_c=2 or 1.5) are able to reproduce the 870um observed emission within the uncertainties for most of the sample. Low values (beta_c<1.3) would be required in NGC0337, NGC1512 and NGC7793. At local scales, we observe a systematic 870um excess when using beta_=2.0. The beta_c=1.5 or the Draine and Li (2007) models can reconcile predictions with observations in part of the disks. Some of the remaining excesses occur towards the centres and can be partly or fully accounted for by non-dust contributions such as CO(3-2) or, to a lesser extent, free-free or synchrotron emission. In three non-barred galaxies, the remaining excesses rather occur in the disk outskirts. This could be a sign of a flattening of the submm slope (and decrease of the effective emissivity index) with radius in these objects.

Spitzer Space Telescope and Herschel Space Observatory imaging of M31 is used, with a physical dust model, to construct maps of dust surface density, dust-to-gas ratio, starlight heating intensity, and PAH abundance, out to R=25kpc. The global dust mass is M_d=5.4x10^7Msol, the global dust/H mass ratio is M_d/M_H=0.0081, and the global PAH abundance is <q_PAH>=0.039. The dust surface density has an inner ring at R=5.6kpc, a maximum at R=11.2kpc, and an outer ring at R=15.1kpc. The dust/gas ratio varies from M_d/M_H=0.026 at the center to ~0.0027 at R=25kpc. From the dust/gas ratio, we estimate the ISM metallicity to vary by a factor ~10, from Z/Zsol=3 at R=0 to ~0.3 at R=25kpc. The dust heating rate parameter <U> peaks at the center, with <U> approx 35, declining to <U> approx 0.25 at R=20kpc. Within the central kpc, the starlight heating intensity inferred from the dust modeling is close to what is estimated from the stars in the bulge. The PAH abundance reaches a peak q_PAH=0.045 at R=11.2kpc. When allowance is made for the different spectrum of the bulge stars, q_PAH for the dust in the central kpc is similar to the overall value of q_PAH in the disk. The silicate-graphite-PAH dust model used here is generally able to reproduce the observed dust spectral energy distribution across M31, but overpredicts 500um emission at R=2-6kpc, suggesting that at R=2-6kpc, the dust opacity varies more steeply with frequency (with beta approx 2.3 between 200 and 600um) than in the model

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

'Normal' galaxies observed at z>6, when the Universe was <1 billion years old, thus far show no evidence of the cold dust that accompanies star formation in the local Universe, where the dust-to-gas mass ratio is 1%. A prototypical example is 'Himiko' (z=6.6), which a mere 840 Myr after the Big Bang is forming stars at a rate of 30-100 Msun/yr, yielding a mass assembly time M^{star}/SFR 150x10^6 yr. Himiko is estimated to have a low fraction (2-3% of the Solar value) of elements heavier than helium (metallicity), and although its gas mass cannot be asserted at this time its dust-to-stellar mass ratio is constrained to be <0.05%. The local galaxy I Zw 18, with a metallicity 4% solar and forming stars less rapidly than Himiko but still vigorously for its mass (M^{star}/SFR 1.6x10^9 yr), is also very dust deficient and perhaps one of the best analogues of primitive galaxies accessible to detailed study. Here we report observations of dust emission from I Zw 18 from which we determine its dust mass to be 450-1800 Msun, yielding a dust-to-stellar mass ratio \approx 10^{-6}-10^{-5} and a dust-to-gas mass ratio 3.2-13x10^{-6}. If I Zw 18 is a reasonable analog of Himiko, then Himiko's dust mass is \approx 50,000 Msun, a factor of 100 below the current upper limit. These numbers are considerably uncertain, but if most high-z galaxies are more like Himiko than like the quasar host SDSS J114816.64+525150.3, then the prospects for detecting the gas and dust in them are much poorer than hitherto anticipated.

We present kiloparsec (kpc) spatial resolution maps of the CO-to-H2 conversion factor (alpha_co) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for alpha_co and DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we can combine maps of dust mass surface density, CO integrated intensity and HI column density to solve for both alpha_co and DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high resolution far-IR maps from the Herschel key program KINGFISH, 12CO J=(2-1) maps from the IRAM 30m large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our alpha_co results on the more typically used 12CO J=(1-0) scale and show using literature measurements that variations in the line ratio do not effect our results. In total, we derive 782 individual solutions for alpha_co and DGR. On average, alpha_co = 3.1 Msun pc^-2 (K km s^-1)^-1 for our sample with a standard deviation of 0.3 dex. Within galaxies we observe a generally flat profile of alpha_co as a function of galactocentric radius. However, most galaxies exhibit a lower alpha_co in the central kpc---a factor of ~2 below the galaxy mean, on average. In some cases, the central alpha_co value can be factors of 5 to 10 below the standard Milky Way (MW) value of alpha_co,MW =4.4 Msun pc^-2 (K km s^-1)^-1. While for alpha_co we find only weak correlations with metallicity, DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate alpha_co for studies of nearby galaxies.

DDSCAT 7.3 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures in Draine & Flatau (2008), and efficient near-field calculations in Flatau & Draine (2012). DDSCAT 7.3 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for user-specified scattering directions. DDSCAT 7.3 can efficiently calculate E and B throughout a user-specified volume containing the target. This User Guide explains how to use DDSCAT 7.3 to carry out electromagnetic scattering calculations, including use of DDPOSTPROCESS, a Fortran-90 code to perform calculations with E and B at user-selected locations near the target. A number of changes have been made since the last release, DDSCAT 7.2 .

Dust has long been identified as a barrier to measuring inherent galaxy properties. However, the link between dust and attenuation is not straightforward and depends on both the amount of dust and its distribution. Herschel imaging of nearby galaxies undertaken as part of the KINGFISH project allows us to map the dust as seen in emission with unprecedented sensitivity and ~1 kpc resolution. We present here new optical integral field unit spectroscopy for eight of these galaxies that provides complementary 100-200 pc scale maps of the dust attenuation through observation of the reddening in both the Balmer decrement and the stellar continuum. The stellar continuum reddening, which is systematically less than that observed in the Balmer decrement, shows no clear correlation with the dust, suggesting that the distribution of stellar reddening acts as a poor tracer of the overall dust content. The brightest HII regions are observed to be preferentially located in dusty regions, and we do find a correlation between the Balmer line reddening and the dust mass surface density for which we provide an empirical relation. Some of the high-inclination systems in our sample exhibit high extinction, but we also find evidence that unresolved variations in the dust distribution on scales smaller than 500 pc may contribute to the scatter in this relation. We caution against the use of integrated A_V measures to infer global dust properties.

We use the near--infrared Br\gamma hydrogen recombination line as a reference star formation rate (SFR) indicator to test the validity and establish the calibration of the {\it Herschel} PACS 70 \mu m emission as a SFR tracer for sub--galactic regions in external galaxies. Br\gamma offers the double advantage of directly tracing ionizing photons and of being relatively insensitive to the effects of dust attenuation. For our first experiment, we use archival CFHT Br\gamma and Ks images of two nearby galaxies: NGC\,5055 and NGC\,6946, which are also part of the {\it Herschel} program KINGFISH (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel). We use the extinction corrected Br\gamma emission to derive the SFR(70) calibration for H{\sc ii} regions in these two galaxies. A comparison of the SFR(70) calibrations at different spatial scales, from 200 pc to the size of the whole galaxy, reveals that about 50% of the total 70\mu m emission is due to dust heated by stellar populations that are unrelated to the current star formation. We use a simple model to qualitatively relate the increase of the SFR(70) calibration coefficient with decreasing region size to the star formation timescale. We provide a calibration for an unbiased SFR indicator that combines the observed H\alpha with the 70 \mu m emission, also for use in H{\sc ii} regions. We briefly analyze the PACS 100 and 160 \mu m maps and find that longer wavelengths are not as good SFR indicators as 70\mu m, in agreement with previous results. We find that the calibrations show about 50% difference between the two galaxies, possibly due to effects of inclination.

We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved SED modelling of these objects using the Draine and Li (2007) dust models and investigate the influence of the addition of SPIRE measurements in the estimation of LTIR. We find that using data up to 250 um leads to local LTIR values consistent with those obtained with a complete coverage (up to 500 um) within 10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the (70-160 um) band. The (24-70 um) fraction decreases with increasing metallicity. The (160-1100 um) submillimeter band can account for up to 25 per cent of the LTIR in metal-rich galaxies. We investigate the correlation between TIR surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three PACS bands can be used as reliable monochromatic estimators of the LTIR, the 100 um band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250um and LTIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated LTIR of nearby galaxies as well as LTIR of z=1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.

We derive the distribution of the synchrotron spectral index across NGC6946 and investigate the correlation between the radio continuum (synchrotron) and far-infrared (FIR) emission using the KINGFISH Herschel PACS and SPIRE data. The radio--FIR correlation is studied as a function of star formation rate, magnetic field strength, radiation field strength, and the total gas surface brightness. The synchrotron emission follows both star-forming regions and the so-called magnetic arms present in the inter-arm regions. The synchrotron spectral index is steepest along the magnetic arms ($\alpha_n \sim 1$), while it is flat in places of giant H{\sc ii} regions and in the center of the galaxy ($\alpha_n \sim 0.6-0.7$). The map of $\alpha_n$ provides an observational evidence for aging and energy loss of cosmic ray electrons propagating in the disk of the galaxy. Variations in the synchrotron--FIR correlation across the galaxy are shown to be a function of both star formation and magnetic fields. We find that the synchrotron emission correlates better with cold rather than with warm dust emission, when the interstellar radiation field is the main heating source of dust. The synchrotron--FIR correlation suggests a coupling between the magnetic field and the gas density. NGC6946 shows a power-law behavior between the total (turbulent) magnetic field strength B and the star formation rate surface density $\Sigma_{\rm SFR}$ with an index of 0.14\,(0.16)$\pm$0.01. This indicates an efficient production of the turbulent magnetic field with the increasing gas turbulence expected in actively star forming regions. The scale-by-scale analysis of the synchrotron--FIR correlation indicates that the ISM affects the propagation of old/diffused cosmic ray electrons, resulting in a diffusion coefficient of $D_0=4.6\times 10^{28}$\,cm$^2$\,s$^{-1}$ for 2.2\,GeV CREs.

Combining Ha and IRAC images of the nearby spiral galaxy NGC 628, we find that between 30-43% of its 8um dust emission is not related to recent star formation. Contributions from dust heated by young stars are separated by identifying HII regions in the Ha map and using these areas as a mask to determine the 8um dust emission that must be due to heating by older stars. Corrections are made for sub-detection-threshold HII regions, photons escaping from HII regions and for young stars not directly associated to HII regions (i.e. 10-100 Myr old stars). A simple model confirms this amount of 8um emission can be expected given dust and PAH absorption cross-sections, a realistic star-formation history, and the observed optical extinction values. A Fourier power spectrum analysis indicates that the 8um dust emission is more diffuse than the Ha emission (and similar to observed HI), supporting our analysis that much of the 8um-emitting dust is heated by older stars. The 8um dust-to-Ha emission ratio declines with galactocentric radius both within and outside of HII regions, probably due to a radial increase in disk transparency. In the course of this work, we have also found that intrinsic diffuse Ha fractions may be lower than previously thought in galaxies, if the differential extinction between HII regions and diffuse regions is taken into account.

The unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 micrometer, commonly attributed to polycyclic aromatic hydrocarbon (PAH) molecules, have been recently ascribed to coal- or kerogen-like organic nanoparticles with a mixed aromatic-aliphatic structure. However, we show in this Letter that this hypothesis is inconsistent with observations. We estimate the aliphatic fraction of the UIE carriers based on the observed intensities of the 3.4 and 6.85 micrometer emission features by attributing them exclusively to aliphatic C-H stretch and aliphatic C-H deformation vibrational modes, respectively. We derive the fraction of carbon atoms in aliphatic form to be <15%. We conclude that the UIE emitters are predominantly aromatic with aliphatic material at most a minor part of the UIE carriers. The PAH model is consistent with astronomical observations and PAHs dominate the strong UIE bands.

We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.

Taking advantage of the sensitivity and angular resolution of the Herschel Space Observatory at far-infrared and submm wavelengths, we aim to characterize the physical properties of cold dust within nearby galaxies and study the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH program, we perform 2 temperature fits of the Spitzer and Herschel photometric data (24 to 500um), with a warm and a cold component, globally and in each resolution element.At global scales, we observe ranges of values for beta_c(0.8 to 2.5) and Tc(19.1 to 25.1K).We compute maps of our parameters with beta fixed or free to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe temperature gradients as a function of radius.When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities.Gathering resolved elements in a Tc-beta_c diagram underlines an anti-correlation between the two parameters.It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar ends (coinciding with molecular gas density enhancements usually found in these locations).We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (beta_c < 2.0) lead to significantly lower dust masses compared to the beta_c < 2.0 case, with dust masses lower by up to 50% if beta_c=1.5 for instance.The working resolution affects our total dust mass estimates: masses increase from global fits to spatially-resolved fits.

We examine Herschel Space Observatory images of one nearby prototypical outer ring galaxy, NGC 1291, and show that the ring becomes more prominent at wavelengths longer than 160um. The mass of cool dust in the ring dominates the total dust mass of the galaxy, accounting for at least 70% of it. The temperature of the emitting dust in the ring (T=19.5+/-0.3K) is cooler than that of the inner galaxy (T=25.7+/-0.7K). We discuss several explanations for the difference in dust temperature, including age and density differences in the stellar populations of the ring versus the bulge.

The presence of ferromagnetic or ferrimagnetic nanoparticles in the interstellar medium would give rise to magnetic dipole radiation at microwave and submm frequencies. Such grains may account for the strong mm-wavelength emission observed from a number of low-metallicity galaxies, including the Small Magellanic Cloud. We show how to calculate the absorption and scattering cross sections for such grains, with particular attention to metallic Fe, magnetite Fe3O4, and maghemite gamma-Fe2O3, all potentially present in the interstellar medium. The rate of Davis-Greenstein alignment by magnetic dissipation is also estimated. We determine the temperature of free-flying magnetic grains heated by starlight and we calculate the polarization of the magnetic dipole emission from both free-fliers and inclusions. For inclusions, the magnetic dipole emission is expected to be polarized orthogonally relative to the normal electric dipole radiation. Finally, we present self-consistent dielectric functions for metallic Fe, magnetite Fe3O4, and maghemite gamma-Fe2O3, enabling calculation of absorption and scattering cross sections from microwave to X-ray wavelengths.

The Small Magellanic Cloud (SMC) has surprisingly strong submm and mm-wavelength emission that is inconsistent with standard dust models, including those with emission from spinning dust. Here we show that the emission from the SMC may be understood if the interstellar dust mixture includes magnetic nanoparticles, emitting magnetic dipole radiation resulting from thermal fluctuations in the magnetization. The magnetic grains can be metallic iron, magnetite Fe3O4, or maghemite gamma-Fe2O3. The required mass of iron is consistent with elemental abundance constraints. The magnetic dipole emission is predicted to be polarized orthogonally to the normal electric dipole radiation if the nanoparticles are inclusions in larger grains. We speculate that other low-metallicity galaxies may also have a large fraction of the interstellar Fe in magnetic materials.

DDSCAT 7.2 is a freely available open-source Fortran-90 software package applying the discrete dipole approximation (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", which can be used to study absorption, scattering, and electric fields around arrays of nanostructures. The DDA approximates the target by an array of polarizable points. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures in Draine & Flatau (2008). Efficient near-field calculations are enabled following Flatau & Draine (2012). DDSCAT 7.2 allows accurate calculations of electromagnetic scattering from targets with size parameters 2*pi*aeff/lambda < 25 provided the refractive index m is not large compared to unity (|m-1| < 2). DDSCAT 7.2 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT 7.2 supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import new target geometries into the code. DDSCAT 7.2 calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. DDSCAT 7.2 calculates E throughout a user-specified rectangular volume containing the target. A Fortran-90 code READNF to read E and P from files created by DDSCAT 7.2 is included in the distribution.

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158{\mu}m/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $\mu$m, [O \one] 63 $\mu$m, and [N \ii] 122 $\mu$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\lesssim$ 0.95. In the warmest regions, where $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $\mu$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.

We measure the optical spectrum of the Diffuse Galactic Light (DGL) -- the local Milky Way in reflection -- using 92,000 blank-sky spectra from the Sloan Digital Sky Survey (SDSS). We correlate the SDSS optical intensity in regions of blank sky against 100 micron intensity independently measured by the COsmic Background Explorer (COBE) and InfraRed Astronomy Satellite (IRAS) satellites, which provides a measure of the dust column density times the intensity of illuminating starlight. The spectrum of scattered light is very blue and shows a clear 4000 Angstrom break and broad Mg b absorption. This is consistent with scattered starlight, and the continuum of the DGL is well-reproduced by a simple radiative transfer model of the Galaxy. We also detect line emission in H\alpha, H\beta, [N II], and [S II], consistent with scattered light from the local interstellar medium (ISM). The strength of [N II] and [S II], combined with upper limits on [O III] and [He I], indicate a relatively soft ionizing spectrum. We find that our measurements of the DGL can constrain dust models, favoring a grain size distribution with relatively few large grains. We also estimate the fraction of high-latitude H\alpha which is scattered to be 19+/-4%.

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we have detected emission in the S(0), S(1), and S(2) pure-rotational (v=0-0) transitions of molecular hydrogen (H2) towards 6 positions in two translucent high Galactic latitude clouds, DCld 300.2-16.9 and LDN 1780. The detection of these lines raises important questions regarding the physical conditions inside low-extinction clouds that are far from ultraviolet radiation sources. The ratio between the S(2) flux and the flux from PAHs at 7.9 microns averages 0.007 for these 6 positions. This is a factor of about 4 higher than the same ratio measured towards the central regions of non-active Galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). Thus the environment of these translucent clouds is more efficient at producing rotationally excited H2 per PAH-exciting photon than the disks of entire galaxies. Excitation analysis finds that the S(1) and S(2) emitting regions are warm (T >300K), but comprise no more than 2% of the gas mass. We find that UV photons cannot be the sole source of excitation in these regions and suggest mechanical heating via shocks or turbulent dissipation as the dominant cause of the emission. The clouds are located on the outskirts of the Scorpius-Centaurus OB association and may be dissipating recent bursts of mechanical energy input from supernova explosions. We suggest that pockets of warm gas in diffuse or translucent clouds, integrated over the disks of galaxies, may represent a major source of all non-active galaxy H2 emission.

Multi-wavelength study of extended astronomical objects requires combining images from instruments with differing point spread functions (PSFs). We describe the construction of convolution kernels that allow one to generate (multi-wavelength) images with a common PSF, thus preserving the colors of the astronomical sources. We generate convolution kernels for the cameras of the Spitzer Space Telescope, Herschel Space Observatory, Galaxy Evolution Explorer (GALEX), Wide-field Infrared Survey Explorer (WISE), ground-based optical telescopes (Moffat functions and sum of Gaussians), and Gaussian PSFs. These kernels allow the study of the Spectral Energy Distribution (SED) of extended objects, preserving the characteristic SED in each pixel. The convolution kernels and the IDL packages used to construct and use them are made publicly available.

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.

Planck is expected to answer crucial questions on the early Universe, but it also provides further understanding on anomalous microwave emission. Electric dipole emission from spinning dust grains continues to be the favored interpretation of anomalous microwave emission. In this paper, we present a method to calculate the rotational emission from small grains of irregular shape with moments of inertia $I_{1}> I_{2}> I_{3}$. We show that a torque-free rotating irregular grain with a given angular momentum radiates at multiple frequency modes. The resulting spinning dust spectrum has peak frequency and emissivity increasing with the degree of grain shape irregularity, which is defined by $I_{1}:I_{2}:I_{3}$. We discuss how the orientation of dipole moment $\bmu$ in body coordinates affects the spinning dust spectrum for different regimes of internal thermal fluctuations. We show that the spinning dust emissivity for the case of strong thermal fluctuations is less sensitive to the orientation of $\bmu$ than in the case of weak thermal fluctuations. We calculate spinning dust spectra for a range of gas density and dipole moment. The effect of compressible turbulence on spinning dust emission intensity is investigated. We show that the emission intensity in a turbulent medium increases by a factor from 1.2-1.4 relative to that in a uniform medium, as sonic Mach number $M_{s}$ increases from 2-7. Finally, spinning dust parameters are constrained by fitting our improved model to five-year {\it Wilkinson Microwave Anisotropy Probe} cross-correlation foreground spectra, for both the H$\alpha$-correlated and 100 $\mu$m-correlated emission spectra.

Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of H\alpha+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, H\alpha+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

Radiation pressure acting on gas and dust causes HII regions to have central densities that are lower than the density near the ionized boundary. HII regions in static equilibrium comprise a family of similarity solutions, parametrized by 3 parameters: beta, gamma, and the product (Q_0 n_rms); beta characterizes the stellar spectrum, gamma characterizes the dust/gas ratio, Q_0 is the ionizing output from the star (photons/s), and n_rms is the rms density within the ionized region. Adopting standard values for beta and gamma, varying (Q_0 n_rms) generates a one-parameter family of density profiles, ranging from nearly uniform density (small Q_0 n_rms), to hollow-sphere HII regions (large Q_0 n_rms). When (Q_0 n_rms) exceeds 10^{52} cm^{-3} s^{-1}, dusty HII regions have conspicuous central cavities, even if no stellar wind is present. For given beta, gamma and (Q_0 n_rms), a fourth quantity, which can be Q_0, determines the overall size and density of the HII region. Examples of density and emissivity profiles are given. We show how quantities of interest -- such as the peak-to-center emissivity ratio, the rms-to-mean density ratio, the edge-to-rms density ratio, and the fraction of the ionizing photons absorbed by the gas -- depend on the 3 parameters beta, gamma, and (Q_0 n_rms). For dusty HII regions, compression of the gas and dust into an ionized shell results in a substantial increase in the fraction of the >13.6 eV photons that actually ionize H (relative to a uniform density HII region with the same dust/gas ratio and density n=n_rms). We discuss the extent to which radial drift of dust grains in HII regions can alter the dust-to-gas ratio. The applicability of these solutions to real HII regions is discussed.

We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

We report the discovery of a supernova (SN) with the highest apparent energy output to date and conclude that it represents an extreme example of the Type IIn subclass. The SN, which was discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, peaked at M_R = -21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this period, SN 2003ma had an integrated bolometric luminosity of 4 x 10^51 ergs, more than any other SN to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen lines, a signature of interaction between the SN and a dense circumstellar medium. The light curves show further evidence for circumstellar interaction, including a long plateau with a shape very similar to the classic SN IIn 1988Z -- however, SN 2003ma is ten times more luminous at all epochs. The fast velocity measured for the intermediate-width H_alpha component (~6000 km/s) points towards an extremely energetic explosion (> 10^52 ergs), which imparts a faster blast-wave speed to the post-shock material and a higher luminosity from the interaction than is observed in typical SNe IIn. Mid-infrared observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at a distance ~0.5 pc from the SN.

Recent observations have found the ratio of H-alpha to free-free radio continuum to be surprisingly high in the diffuse ionized ISM (the so-called WIM), corresponding to an electron temperature of only ~3000K. Such low temperatures were unexpected in gas that was presumed to be photoionized. We consider a 3-component model for the observed diffuse emission, consisting of a mix of (1) photoionized gas, (2) gas that is recombining and cooling, and (3) cool H I gas. This model can successfully reproduce the observed intensities of free-free continuum, H-alpha, and collisionally-excited lines such as NII 6583. To reproduce the low observed value of free-free to H-alpha, the PAH abundance in the photoionized regions must be lowered by a factor ~3, and ~20% of the diffuse H-alpha must be reflected from dust grains, as suggested by Wood & Reynolds (1999).

We study a 24\,$\mu$m selected sample of 330 galaxies observed with the Infrared Spectrograph for the 5\,mJy Unbiased Spitzer Extragalactic Survey. We estimate accurate total infrared luminosities by combining mid-IR spectroscopy and mid-to-far infrared photometry, and by utilizing new empirical spectral templates from {\em Spitzer} data. The infrared luminosities of this sample range mostly from 10$^9$L$_\odot$ to $10^{13.5}$L$_\odot$, with 83% in the range 10$^{10}$L$_\odot$$<$L$_{\rm IR}$$<10^{12}$L$_\odot$. The redshifts range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 6.2\,$\mu$m aromatic feature have a bimodal distribution. We use the 6.2\,$\mu$m PAH EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have steeper mid-IR to far-IR spectral slopes and lower L$_{\rm IR}$ and redshifts. The low EW objects (AGN-dominated) tend to have less steep spectral slopes and higher L$_{\rm IR}$ and redshifts. This dichotomy leads to a gross correlation between EW and slope, which does not hold within either group. AGN dominated sources tend to have lower log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios than star-forming galaxies, possibly due to preferential destruction of the smaller aromatics by the AGN. The log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios for star-forming galaxies are lower in our sample than the ratios measured from the nuclear spectra of nearby normal galaxies, most probably indicating a difference in the ionization state or grain size distribution between the nuclear regions and the entire galaxy. Finally, we provide a calibration relating the monochromatic 5.8, 8, 14 and 24um continuum or Aromatic Feature luminosity to L$_{\rm IR}$ for different types of objects.

We review the main advances brought by the Spitzer Space Telescope in the field of nearby galaxies studies, concentrating on a few subject areas, including: (1) the physics of the Polycyclic Aromatic Hydrocarbons that generate the mid-infrared features between ~3.5 micron and ~20 micron; (2) the use of the mid- and far-infrared emission from galaxies as star formation rate indicators; and (3) the improvement of mid-infrared diagnostics to discriminate between thermal (star-formation) and non-thermal (AGN) emission in galaxies and galaxy centers.

Observations continue to support the interpretation of the anomalous microwave foreground as electric dipole radiation from spinning dust grains as proposed by Draine and Lazarian (1998ab). In this paper we present a refinement of the original model by improving the treatment of a number of physical effects. First, we consider a disk-like grain rotating with angular velocity at an arbitrary angle with respect to the grain symmetry axis and derive the rotational damping and excitation coefficients arising from infrared emission, plasma-grain interactions and electric dipole emission. The angular velocity distribution and the electric dipole emission spectrum for grains is calculated using the Langevin equation, for cases both with and without fast internal relaxation. Our results show that, the peak emissivity of spinning dust, compared to earlier studies, increases by a factor of ~2 for the Warm Neutral Medium (WNM), the Warm Ionized Medium (WIM), the Cold Neutral Medium (CNM) and the Photodissociation Region (PDR), and by a factor ~4 for Reflection Nebulae (RN). The frequency at the emission peak also increases by factors ~1.4 to ~2 for these media. The increased emission and peak frequency result from the non-sphericity of grain shape and from the anisotropy in damping and excitation along directions parallel and perpendicular to the grain symmetry axis. Second, we provide a detailed numerical study including transient spin-up of grains by single-ion collisions. The impulses broaden the emission spectrum and increase the peak emissivity for the CNM, WNM and WIM. In addition, we present an improved treatment of rotational excitation and damping by infrared emission.

The physical sizes of supernova remnants (SNRs) in a number of nearby galaxies follow an approximately linear cumulative distribution, contrary to what is expected for decelerating shock fronts. This has been attributed to selection effects, or to a majority of SNRs propagating in "free expansion", at constant velocity, into a tenuous ambient medium. We compile a list of 77 known SNRs in the Magellanic Clouds (MCs), and argue that they are a fairly complete record of the SNe that have exploded over the last ~20kyr, with most now in the adiabatic, Sedov phase of their expansions. The roughly linear cumulative size distribution (uniform in a differential distribution) can result from the combination of a deceleration during this phase, a transition to a radiation-loss-dominated phase at a radius that depends on the local gas density, and a distribution of ambient densities varying roughly as rho^{-1}. This explanation is supported by the observed -1 power-law distributions of three independent tracers of density: HI column density, Halpha surface brightness, and star formation rate from resolved stellar populations. In this picture, the observed cutoff at r~30 pc in the SNR size distribution is due to a minimum in the mean ambient gas density in the regions where supernovae (SNe) explode. We show that M33 has a SNR size distribution similar to that of the MCs, suggesting these features, and their explanation, may be universal. In a companion paper (Maoz & Badenes 2010), we use our sample of SNRs as an effective "SN survey" to calculate the SN rate and delay time distribution in the MCs. The hypothesis that most SNRs are in free expansion, rather than in the Sedov phase of their evolution, would result in SN rates that are in strong conflict with independent measurements, and with basic stellar evolution theory.

Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.

Stellar density and bar strength should affect the temperatures of the cool (T ~ 20-30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 microns. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 +/- 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23(+/-3)% hotter for morphological types earlier than Sc, and only 9(+/-3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar.

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $\mu$m, [OIII] 88 $\mu$m, [NII] 122 $\mu$m, [CII] 158 $\mu$m and [NII] 205 $\mu$m. We observe that in the [OI] 63 $\mu$m, [OIII] 88 $\mu$m, and [NII] 122 $\mu$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $\mu$m/[CII] 158 $\mu$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $\mu$m and [OIII] 88 $\mu$m maps. Variations of the [NII] 122 $\mu$m/[NII] 205 $\mu$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

We present the results of modeling dust SEDs across the SMC with the aim of mapping the distribution of PAHs in a low-metallicity environment. Using Spitzer Survey of the SMC (S3MC) photometry from 3.6-160 um over the main star-forming regions of the Wing and Bar along with spectral mapping from 5-38 um from the Spitzer Spectroscopic Survey of the SMC (S4MC) in selected regions, we model the dust SED and emission spectrum to determine the fraction of dust in PAHs across the SMC. We use the regions of overlapping photometry and spectroscopy to test the reliability of the PAH fraction as determined from SED fits alone. The PAH fraction in the SMC is low compared to the Milky Way and variable--with relatively high fractions (q_PAH~1-2%) in molecular clouds and low fractions in the diffuse ISM (<q_PAH>=0.6%). We use the map of PAH fraction across the SMC to test a number of ideas regarding the production, destruction and processing of PAHs in the ISM. We find weak or no correlation between the PAH fraction and the distribution of carbon AGB stars, the location of supergiant H I shells and young SN remnants, or the turbulent Mach number. We find that the PAH fraction is correlated with CO intensity, peaks in the dust surface density and the molecular gas surface density as determined from 160 um emission. The PAH fraction is high in regions of active star-formation, as predicted by its correlation with molecular gas, but is suppressed in H II regions. Because the PAH fraction in the diffuse ISM is generally very low--in accordance with previous work on modeling the SED of the SMC--and the PAH fraction is relatively high in molecular regions, we suggest that PAHs are destroyed in the diffuse ISM of the SMC and/or PAHs are forming in molecular clouds. We discuss the implications of these observations for our understanding of the PAH life cycle, particularly in low-metallicity galaxies.

(Abridged) Spitzer data at 24, 70, and 160 micron and ground-based H-alpha images are analyzed for a sample of 189 nearby star-forming and starburst galaxies to investigate whether reliable star formation rate (SFR) indicators can be defined using the monochromatic infrared dust emission centered at 70 and 160 micron. We compare recently published recipes for SFR measures using combinations of the 24 micron and observed H-alpha luminosities with those using 24 micron luminosity alone. From these comparisons, we derive a reference SFR indicator for use in our analysis. Linear correlations between SFR and the 70 and 160 micron luminosity are found for L(70)>=1.4x10^{42} erg/s and L(160)>=2x10^{42} erg/s, corresponding to SFR>=0.1-0.3 M_sun/yr. Below those two luminosity limits, the relation between SFR and 70 micron (160 micron) luminosity is non-linear and SFR calibrations become problematic. The dispersion of the data around the mean trend increases for increasing wavelength, becoming about 25% (factor ~2) larger at 70 (160) micron than at 24 micron. The increasing dispersion is likely an effect of the increasing contribution to the infrared emission of dust heated by stellar populations not associated with the current star formation. The non-linear relation between SFR and the 70 and 160 micron emission at faint galaxy luminosities suggests that the increasing transparency of the interstellar medium, decreasing effective dust temperature, and decreasing filling factor of star forming regions across the galaxy become important factors for decreasing luminosity. The SFR calibrations are provided for galaxies with oxygen abundance 12+Log(O/H)>8.1. At lower metallicity the infrared luminosity no longer reliably traces the SFR because galaxies are less dusty and more transparent.

DDSCAT 7.1 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures (and near-field calculations) in Draine & Flatau (2008). DDSCAT 7.1 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix. This User Guide explains how to use DDSCAT 7.1 to carry out electromagnetic scattering calculations. DDfield, a Fortran-90 code to calculate E and B at user-selected locations near the target, is included in the distribution. A number of changes have been made since the last release, DDSCAT 7.0 .

We combine new ASTE/AzTEC 1.1 mm maps of the galaxy pair NGC 1512/1510 with archival Spitzer IRAC and MIPS images covering the wavelength range 3.6--160 um from the SINGS project to derive accurate dust masses in each galaxy, and in sub--galactic regions in NGC 1512. The two galaxies form a pair consisting of a large, high--metallicity spiral (NGC 1512) and a low metallicity, blue compact dwarf (NGC 1510). The derived total dust masses are (2.4+/-0.6) 10^7 Msun and (1.7+/-3.6) 10^5 Msun for NGC 1512 and NGC 1510, respectively. The derived ratio of dust mass to H I gas mass for the galaxy pair (0.0034) is much lower than expected, while regions within NGC 1512, specifically the central region and the arms, do not show such unusually low ratios; furthermore, the dust--to--gas ratio is within expectations for NGC 1510. These results suggest that a fraction of the H I included in the determination of the M_d/M_HI ratio of the NGC 1512/1510 pair is not associated with the star forming disks/regions of either galaxy. We also perform simple two--temperature (warm + cold) modified--blackbody fits to the far--infrared/mm data. The warm dust temperature of the low--metallicity NGC 1510 (~36 K) is substantially higher than that of the high--metallicity NGC 1512 (~24 K). In both galaxies, a substantial fraction (>93%) of the total dust mass is in a cool dust component, with temperatures 14--16 K for NGC 1512 and 15--24 K for NGC 1510. In contrast, the warm dust in the arms and central regions of NGC 1512 represents a much larger fraction of the total dust content, in agreement with the fact that all three regions have higher specific star formation rates than the average in the galaxy.

We present a detailed analysis of the radial distribution of dust properties in the SINGS sample, performed on a set of UV, IR and HI surface brightness profiles, combined with published molecular gas profiles and metallicity gradients. The internal extinction, derived from the TIR-to-FUV luminosity ratio, decreases with radius, and is larger in Sb-Sbc galaxies. The TIR-to-FUV ratio correlates with the UV spectral slope beta, following a sequence shifted to redder UV colors with respect to that of starbursts. The star formation history (SFH) is identified as the main driver of this departure. We have also derived radial profiles of the total dust mass surface density, the fraction of the dust mass contributed by PAHs, the fraction of the dust mass heated by very intense starlight and the intensity of the radiation field heating the grains. The dust profiles are exponential, their radial scale-length being constant from Sb to Sd galaxies (only ~10% larger than the stellar scale-length). Many S0/a-Sab galaxies have central depressions in their dust radial distributions. The PAH abundance increases with metallicity for 12+\log(O/H)<9, and at larger metallicities the trend flattens and even reverses, with the SFH being a plausible underlying driver for this behavior. The dust-to-gas ratio is also well correlated with metallicity and therefore decreases with galactocentric radius.

A photon of wavelength lambda ~1 micron interacting with a dust grain of radius a_p ~ 1 mm (a "pebble") undergoes scattering in the forward direction, largely within a small characteristic diffraction angle theta_s ~ lambda/a_p ~ 100". Though mm-size dust grains contribute negligibly to the interstellar medium's visual extinction, the signal they produce in scattered light may be detectable, especially for variable sources. Observations of light scattered at small angles allows for the direct measurement of the large grain population; variable sources can also yield tomographic information of the interstellar medium's mass distribution. The ability to detect brilliant pebble halos require a careful understanding of the instrument PSF.

GRB 080607 has provided the first strong observational signatures of molecular absorption bands toward any galaxy hosting a gamma-ray burst. Despite the identification of dozens of features as belonging to various atomic and molecular (H_2 and CO) carriers, many more absorption features remained unidentified. Here we report on a search among these features for absorption from vibrationally-excited H_2, a species that was predicted to be produced by the UV flash of a GRB impinging on a molecular cloud. Following a detailed comparison between our spectroscopy and static, as well as dynamic, models of H_2* absorption, we conclude that a column density of 10^{17.5+-0.2} cm^{-2} of H_2* was produced along the line of sight toward GRB 080607. Depending on the assumed amount of dust extinction between the molecular cloud and the GRB, the model distance between the two is found to be in the range 230--940 pc. Such a range is consistent with a conservative lower limit of 100 pc estimated from the presence of Mg I in the same data. These distances show that substantial molecular material is found within hundreds of pc from GRB 080607, part of the distribution of clouds within the GRB host galaxy.

We present a catalog of mid-infrared (MIR) spectra of 150 infrared (IR) luminous galaxies in the Spitzer extragalactic first look survey obtained with IRS on board Spitzer. The sample is selected to be brighter than ~0.9 mJy at 24 micron and it has a z distribution in the range [0.3,3.5] with a peak at z=1. It primarily comprises ultraluminous IR galaxies at z>1 and luminous IR galaxies at z<1, as estimated from their monochromatic 14 micron luminosities. The number of sources with spectra that are dominated by an active galactic nucleus (AGN) continuum is 49, while 39 sources have strong, star-formation related features. For this classification, we used the equivalent width (EW) of the 11.3 micron polycyclic aromatic hydrocarbon (PAH) feature. Several intermediate/high z starbursts have higher PAH EW than local ULIRGs. An increase in the AGN activity is observed with increasing z and luminosity, based on the decreasing EW of PAHs and the increasing [NeIII]/[NeII] ratio. Spectral stacking leads to the detection of the 3.3 micron PAH, the H2 0-0 S(1) and S(3) lines, and the [NeV] line. We observe differences in the flux ratios of PAHs in the stacked spectra of IR-luminous galaxies with z or luminosity, which are not due to extinction effects. When placing the observed galaxies on IR color-color diagrams, we find that the wedge defining AGN comprises most sources with continuum-dominated spectra, but also contains many starbursts. The comparison of the 11.3 micron PAH EW and the H-band effective radius, measured from HST data, indicates that sources with EW>2 micron, are typically more extended than ~3 kpc. However, there is no strong correlation between the MIR spectral type and the near-IR extent of the sources. [Abridged].

We present theoretical calculations of the X-ray scattering properties of porous grain aggregates with olivine monomers. The small and large angle scattering properties of these aggregates are governed by the global structure and substructure of the grain, respectively. We construct two diagnostics, R_X and T_X, based on the optical and X-ray properties of the aggregates, and apply them to a Chandra measurement of the dust halo around the Galactic binary GX13+1. Grain aggregates with porosities higher than about 0.55 are ruled out. Future high-precision observations of X-ray dust haloes together with detailed modeling of the X-ray scattering properties of porous grain mixtures will further constrain the presence of porous grain aggregates in a given dust population.

DDSCAT 7.0 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures (and near-field calculations) in Draine & Flatau (2008). DDSCAT 7.0 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix. This User Guide explains how to use DDSCAT 7.0 to carry out electromagnetic scattering calculations. DDfield, a Fortran-90 code DDfield to calculate E and B at user-selected locations near the target, is included in the distribution.

The wavelength dependences of interstellar extinction and polarization, supplemented by observed elemental abundances and the spectrum of infrared emission from dust heated by starlight, strongly constrain dust models. One dust model that appears to be consistent with observations is presented. To reproduce the observed extinction, the model consumes the bulk of interstellar Mg, Si, and Fe (in amorphous silicates), and a substantial fraction of C (in carbonaceous material), with size distributions and alignment adjusted to match observations. The composition, structure, and size distribution of interstellar grains is the result of injection of dust from stellar outflows into the interstellar medium (ISM), followed by destruction, growth, coagulation, and photoprocessing of interstellar grains. The balance among these poorly-understood processes is responsible for the mix of solid material present in the ISM. Most interstellar grain material present in the diffuse ISM must be grown in the ISM. The amorphous silicate and carbonaceous materials that form the bulk of interstellar dust must therefore be the result of grain growth in the presence of ultraviolet radiation. Dust in high-z systems such as J1148+5251 is also produced primarily in the ISM, with supernova-produced dust contributing only a small fraction of the total dust mass.

We study the polarized infrared emission by Polycyclic Aromatic Hydrocarbons (PAHs), when anisotropically illuminated by UV photons. PAH molecules are modeled as planar disks with in-plane and out-of-plane vibrational dipoles. As first pointed out by Leger (1988), infrared emission features resulting from in-plane and out-of-plane modes should have orthogonal polarization directions. We show analytically how the degree of polarization depends on the viewing geometry and the molecule's internal alignment between principal axis of inertia and angular momentum, which gets worse after photon absorption. Longer wavelength features, emitted after better internal alignment is recovered, should be more strongly polarized. The degree of polarization for uni-directional illumination (e.g., by a star) is larger than for diffuse illumination (e.g., by a disk galaxy), all else being equal. For PAHs in the Cold Neutral Medium, the predicted polarization is probably too small to distinguish from the contribution of linear dichroism by aligned foreground dust. The level of polarization predicted for PAH emission from the Orion Bar is only ~0.06% at 3.3 microns; Sellgren et al. (1988) report a much larger value, 0.86+-0.28%, which suggests that the smallest PAHs may have moderately suprathermal rotation rates. Future observations of (or upper limits on) the degree of polarization for the Orion Bar or for dust above edge-on galaxies (e.g., NGC 891 or M82) may constrain the internal alignment of emitting PAHs, thus providing clues to their rotational dynamics.

We study the light scattering properties of random ballistic aggregates constructed in Shen et al. (Paper I). Using the discrete-dipole-approximation, we compute the scattering phase function and linear polarization for random aggregates with various sizes and porosities, and with two different compositions: 100% silicate and 50% silicate-50% graphite. We investigate the dependence of light scattering properties on wavelength, cluster size and porosity using these aggregate models. We find that while the shape of the phase function depends mainly on the size parameter of the aggregates, the linear polarization depends on both the size parameter and the porosity of the aggregates, with increasing degree of polarization as the porosity increases. Contrary to previous studies, we argue that monomer size has negligible effects on the light scattering properties of ballistic aggregates, as long as the constituent monomer is smaller than the incident wavelength up to 2*pi*a_0/lambda\sim 1.6 where a_0 is the monomer radius. Previous claims for such monomer size effects are in fact the combined effects of size parameter and porosity. Finally, we present aggregate models that can reproduce the phase function and polarization of scattered light from the AU Mic debris disk and from cometary dust, including the negative polarization observed for comets at scattering angles 160<theta<180 deg. These aggregates have moderate porosities, P\sim 0.6, and are of sub-micron-size for the debris disk case, or micron-size for the comet case.

We present an analysis of far-infrared dust emission from diffuse cirrus clouds. This study is based on serendipitous observations at 160 microns at high galactic latitude with the Multiband Imaging Photometer (MIPS) onboard the Spitzer Space Telescope by the Spitzer Infrared Nearby Galaxies Survey (SINGS). These observations are complemented with IRIS data at 100 and 60 microns and constitute one of the most sensitive and unbiased samples of far infrared observations at small scale of diffuse interstellar clouds. Outside regions dominated by the cosmic infrared background fluctuations, we observe a substantial scatter in the 160/100 colors from cirrus emission. We compared the 160/100 color variations to 60/100 colors in the same fields and find a trend of decreasing 60/100 with increasing 160/100. This trend can not be accounted for by current dust models by changing solely the interstellar radiation field. It requires a significant change of dust properties such as grain size distribution or emissivity or a mixing of clouds in different physical conditions along the line of sight. These variations are important as a potential confusing foreground for extragalactic studies.

We report on our knowledge of Galactic foregrounds, as well as on how a CMB satellite mission aiming at detecting a primordial B-mode signal (CMBPol) will contribute to improving it. We review the observational and analysis techniques used to constrain the structure of the Galactic magnetic field, whose presence is responsible for the polarization of Galactic emissions. Although our current understanding of the magnetized interstellar medium is somewhat limited, dramatic improvements in our knowledge of its properties are expected by the time CMBPol flies. Thanks to high resolution and high sensitivity instruments observing the whole sky at frequencies between 30 GHz and 850 GHz, CMBPol will not only improve this picture by observing the synchrotron emission from our galaxy, but also help constrain dust models. Polarized emission from interstellar dust indeed dominates over any other signal in CMBPol's highest frequency channels. Observations at these wavelengths, combined with ground-based studies of starlight polarization, will therefore enable us to improve our understanding of dust properties and of the mechanism(s) responsible for the alignment of dust grains with the Galactic magnetic field. CMBPol will also shed new light on observations that are presently not well understood. Morphological studies of anomalous dust and synchrotron emissions will indeed constrain their natures and properties, while searching for fluctuations in the emission from heliospheric dust will test our understanding of the circumheliospheric interstellar medium. Finally, acquiring more information on the properties of extra-Galactic sources will be necessary in order to maximize the cosmological constraints extracted from CMBPol's observations of CMB lensing. (abridged)

We characterize spinning dust emission in the warm ionized medium by comparing templates of Galactic dust and Halpha with the 5-year maps from the Wilkinson Microwave Anisotropy Probe. The Halpha-correlated microwave emission deviates from the thermal bremsstrahlung (free-free) spectrum expected for ionized gas, exhibiting an additional broad bump peaked at ~40 GHz which provides ~20% of the peak intensity. We confirm that the bump is consistent with a modified Draine & Lazarian (1998) spinning dust model, though the peak frequency of the emission is somewhat lower than the 50 GHz previously claimed. This frequency shift results from systematic errors in the large-scale modes of the 3-year WMAP data which have been corrected in the 5-year data release. We show that the bump is not the result of errors in the Halpha template by analyzing regions of high free-free intensity, where the WMAP K-band map may be used as the free-free template. We rule out a pure free-free spectrum for the Halpha-correlated emission at high confidence: ~27sigma for the nearly full-sky fit, even after marginalizing over the CMB cross-correlation bias. We also extend the previous analysis by searching the parameter space of the Draine & Lazarian model but letting the amplitude float. The best fit for reasonable values of the characteristic electric dipole moment and density requires an amplitude factor of ~0.3. This suggests that small PAHs in the warm ionized medium are depleted by a factor of ~3.

Measurements by dust detectors on interplanetary spacecraft appear to indicate a substantial flux of interstellar particles with masses exceeding 10^{-12}gram. The reported abundance of these massive grains cannot be typical of interstellar gas: it is incompatible with both interstellar elemental abundances and the observed extinction properties of the interstellar dust population. We discuss the likelihood that the Solar System is by chance located near an unusual concentration of massive grains and conclude that this is unlikely, unless dynamical processes in the ISM are responsible for such concentrations. Radiation pressure might conceivably drive large grains into "magnetic valleys". If the influx direction of interstellar gas and dust is varying on a ~10 yr timescale, as suggested by some observations, this would have dramatic implications for the small-scale structure of the interstellar medium.

Polarized far-infrared and submillimeter emission is calculated for models of nonspherical dust grains that are constrained to reproduce the observed wavelength-dependent extinction and polarization of starlight. For emission from regions where the magnetic field is perpendicular to the line-of-sight, the far-infrared emission is expected to have substantial linear polarization at wavelengths longer than 100 microns, but the degree of linear polarization, and its variation with wavelength, is model-dependent. Models in which the starlight polarization is produced by both amorphous silicate and graphite grains have linear polarizations between 6% and 10% at wavelengths longer than 100 microns, but for some models in which only silicate grains are spheroidal, the linear polarization increases from about 3% at 100 microns to about 15% at 1 mm. We briefly discuss the implications of these results for removal of the polarized dust emission from maps of the Cosmic Microwave Background, as well as the possibility of discriminating among interstellar dust models based on observations of far-infrared and submillimeter linear polarization.

The discrete-dipole approximation (DDA) is a powerful method for calculating absorption and scattering by targets that have sizes smaller than or comparable to the wavelength of the incident radiation. The DDA can be extended to targets that are singly- or doubly-periodic. We generalize the scattering amplitude matrix and the 4 x 4 Mueller matrix to describe scattering by singly- and doubly-periodic targets, and show how these matrices can be calculated using the DDA. The accuracy of DDA calculations using the open-source code DDSCAT is demonstrated by comparison to exact results for infinite cylinders and infinite slabs. A method for using the DDA solution to obtain fields within and near the target is presented, with results shown for infinite slabs.

We investigate the scattering and absorption of light by random ballistic aggregates of spherical monomers. We present a general measure for the porosity of an irregular particle. Three different classes of ballistic aggregates are considered, with different degrees of porosity. Scattering and absorption cross sections are calculated, using the discrete dipole approximation (DDA), for grains of three compositions (50% silicate and 50% graphite; 50% silicate and 50% amorphous carbon; and 100% silicate), for wavelengths from 0.1 micron to 4 micron. For fixed particle mass, increased porosity increases the extinction at short wavelengths, but decreases the extinction at wavelengths long compared to the overall aggregate size. Scattering and absorption cross sections are insensitive to monomer size as long as the constituent monomers are small compared with the incident wavelength. We compare our accurate DDA results with two other approximations: the analytical multi-layer sphere (MLS) model and effective medium theory (EMT). For high porosity and/or absorptive materials, the MLS model does not provide a good approximation for scattering and absorption by ballistic aggregates. The EMT method provides a much better approximation than the MLS model for these aggregates, with a typical difference less than 20% in extinction and scattering cross sections compared with DDA results, for all types, compositions and wavelengths probed in this study.

We have modeled optical to far infrared images, photometry, and spectroscopy of the object known as Gomez's Hamburger. We reproduce the images and spectrum with an edge-on disk of mass 0.3M_sun and radius 1600AU, surrounding an A0 III star at a distance of 280pc. Our mass estimate is in excellent agreement with recent CO observations. However, our distance determination is more than an order of magnitude smaller than previous analyses which inaccurately interpreted the optical spectrum. To accurately model the infrared spectrum we have extended our Monte Carlo radiation transfer codes to include emission from polycyclic aromatic hydrocarbon (PAH) molecules and very small grains (VSG). We do this using pre-computed PAH/VSG emissivity files for a wide range of values of the mean intensity of the exciting radiation field. When Monte Carlo energy packets are absorbed by PAHs/VSGs we reprocess them to other wavelengths by sampling from the emissivity files, thus simulating the absorption and re-emission process without reproducing lengthy computations of statistical equilibrium, excitation and de-excitation in the complex many level molecules. Using emissivity lookup tables in our Monte Carlo codes gives the flexibility to use the latest grain physics calculations of PAH/VSG emissivity and opacity that are being continually updated in the light of higher resolution infrared spectra. We find our approach gives a good representation of the observed PAH spectrum from the disk of Gomez's Hamburger. Our models also indicate the PAHs/VSGs in the disk have a larger scaleheight than larger radiative equilibrium grains, providing evidence for dust coagulation and settling to the midplane.

We investigate the relations among the stellar continuum-subtracted 8 micron polycyclic aromatic hydrocarbon (PAH 8 micron) emission, 24 micron hot dust emission, and 160 micron cold dust emission in fifteen nearby face-on spiral galaxies in the Spitzer Infrared Nearby Galaxies Survey sample. The relation between PAH 8 and 24 micron emission measured in ~2 kpc regions is found to exhibit a significant amount of scatter, and strong spatial variations are observed in the (PAH 8 micron)/24 micron surface brightness ratio. In particular, the (PAH 8 micron)/24 micron surface brightness ratio is observed to be high in the diffuse interstellar medium and low in bright star-forming regions and other locations with high 24 micron surface brightness. PAH 8 micron emission is found to be well-correlated with 160 micron emission on spatial scales of ~2 kpc, and the (PAH 8 micron)/160 micron surface brightness ratio is generally observed to increase as the 160 micron surface brightness increases. These results suggest that the PAHs are associated with the diffuse, cold dust that produces most of the 160 micron emission in these galaxies, and the variations in the (PAH 8 micron)/160 micron ratio may generally be indicative of either the intensity or the spectrum of the interstellar radiation field that is heating both the PAHs and the diffuse interstellar dust.

We have studied the relationship between the star formation rate (SFR) surface density and gas surface density in the spiral galaxy M51a (NGC 5194), using multi-wavelength data obtained as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS). We introduce a new SFR index based on a linear combination of H-alpha emission-line and 24 micron continuum luminosities, that provides reliable extinction-corrected ionizing fluxes and SFR densities over a wide range of dust attenuations. The combination of these extinction-corrected SFR densities with aperture synthesis HI and CO maps has allowed us to probe the form of the spatially-resolved star formation law on scales of 0.5 to 2 kpc. We find that the resolved SFR vs gas surface density relation is well represented by a Schmidt power law, which is similar in form and dispersion to the disk-averaged Schmidt law. We observe a comparably strong correlation of the SFR surface density with the molecular gas surface density, but no significant correlation with the surface density of atomic gas. The best-fitting slope of the Schmidt law varies from N = 1.37 to 1.56, with zeropoint and slope that change systematically with the spatial sampling scale. We tentatively attribute these variations to the effects of areal sampling and averaging of a nonlinear intrinsic star formation law. Our data can also be fitted by an alternative parametrization of the SFR surface density in terms of the ratio of gas surface density to local dynamical time, but with a considerable dispersion.

The quantity of dust in a spiral disk can be estimated using the dust's typical emission or the extinction of a known source. In this paper, we compare two techniques, one based on emission and one on absorption, applied on sections of fourteen disk galaxies. The two measurements reflect, respectively the average and apparent optical depth of a disk section. Hence, they depend differently on the average number and optical depth of ISM structures in the disk. The small scale geometry of the cold ISM is critical for accurate models of the overall energy budget of spiral disks. ISM geometry, relative contributions of different stellar populations and dust emissivity are all free parameters in galaxy Spectral Energy Distribution (SED) models; they are also sometimes degenerate, depending on wavelength coverage. Our aim is to constrain typical ISM geometry. The apparent optical depth measurement comes from the number of distant galaxies seen in HST images through the foreground disk. We measure the IR flux in images from the {\it Spitzer} Infrared Nearby Galaxy Survey in the same section of the disk that was covered by HST. A physical model of the dust is fit to the SED to estimate the dust surface density, mean temperature, and brightness in these disk sections. The surface density is subsequently converted into the average optical depth estimate. The two measurements generally agree. The ratios between the measured average and apparent optical depths of the disk sections imply optically thin clouds in these disks. Optically thick disks, are likely to have more than a single cloud along the line-of-sight.

With the goal of investigating the degree to which the mid-infrared emission traces the star formation rate (SFR), we analyze Spitzer 8 um and 24 um data of star-forming regions in a sample of 33 nearby galaxies with available HST/NICMOS images in the Paschen-alpha (1.8756 um) emission line. The galaxies are drawn from the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample, and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and Luminous Infrared Galaxies are also included in the analysis. Both the stellar-continuum-subtracted 8 um emission and the 24 um emission correlate with the extinction-corrected Pa-alpha line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed non-linear trend of the 24 um emission versus number of ionizing photons, including the modest deficiency of 24 um 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 um emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 um emission is contributed, in larger measure than the 24 um emission, by dust heated by non-ionizing stellar populations, in agreement with previous findings. Two SFR calibrations, one using the 24 um emission and the other using a combination of the 24 um and H-alpha luminosities (Kennicutt et al. 2007), are presented. No calibration is presented for the 8 um emission, because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by on-going star formation.

(simplified) Results on the properties of warm H2 in 57 normal galaxies are derived from H2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of H2, the most abundant constituent of the molecular ISM, to more common systems (L_FIR = e7 to 6e10 L_sun) of all morphological and nuclear types. The S(1) transition is securely detected in the nuclear regions of 86% of SINGS galaxies with stellar masses above 10^9.5 M_sun. The derived column densities of warm H2 (T > ~100 K), even though averaged over kiloparsec-scale areas, are commensurate with those of resolved PDRs; the median of the sample is 3e20 cm-2. They amount to between 1% and >30% of the total H2. The power emitted in the sum of the S(0) to S(2) transitions is on average 30% of the [SiII] line power, and ~4e-4 of the total infrared power (TIR) within the same area for star-forming galaxies, which is consistent with excitation in PDRs. The fact that H2 emission scales tightly with PAH emission, even though the average radiation field intensity varies by a factor ten, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. A large fraction of the 25 LINER/Sy targets, however, strongly depart from the rest of the sample, in having warmer H2 in the excited states, and an excess of H2 emission with respect to PAHs, TIR and [SiII]. We propose a threshold in H2 to PAH power ratios, allowing the identification of low-luminosity AGNs by an excess H2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect, in nearly half the star-forming targets, non-equilibrium ortho to para ratios, consistent with FUV pumping combined with incomplete ortho-para thermalization by collisions, or possibly non-equilibrium PDR fronts advancing into cold gas.

Physical dust models are presented for 65 galaxies in the SINGS survey that are strongly detected in the four IRAC bands and three MIPS bands. For each galaxy we estimate (1) the total dust mass, (2) the fraction of the dust mass contributed by PAHs, and (3) the intensity of the starlight heating the dust grains. We find that spiral galaxies have dust properties resembling the dust in the local region of the Milky Way, with similar dust-to-gas ratio, and similar PAH abundance. The observed SEDs, including galaxies with SCUBA photometry, can be reproduced by dust models that do not require "cold" (T<10K) dust. The dust-to-gas ratio is observed to be dependent on metallicity. In the interstellar media of galaxies with A_O=12+log(O/H)>8.1, grains contain a substantial fraction of interstellar Mg, Si and Fe. Galaxies with A_O<8.1 and extended HI envelopes in some cases appear to have global dust-to-gas ratios that are low for their measured oxygen abundance, but in the regions where infrared emission is detected, the dust-to-gas ratio generally appears to be consistent with a substantial fraction of interstellar Mg, Si, and Fe being contained in dust. The PAH index q_PAH -- the fraction of the dust mass in the form of PAHs -- correlates with metallicity. The nine galaxies in our sample with A_O<8.1 have a median q_PAH=1.0%, whereas galaxies with A_O>8.1 have a median q_PAH=3.55%. The derived dust masses favor a value X_CO approx 4e20 cm^{-2}(K kms)^{-1} for the CO to H_2 conversion factor. Except for some starbursting systems (Mrk33, Tolo89, NGC3049), dust in the diffuse ISM dominates the IR power.

[abridged] We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the M81 group from the SINGS and THINGS surveys. The Spitzer observations provide a first glimpse of the nature of the non-atomic ISM in these metal-poor (Z~0.1 Z_sun), quiescent (SFR~0.001-0.1 M_sun/yr) dwarf galaxies. Dust emission is detected in five out of the seven targets. Most detected dust emission is restricted to HI column densities >1x10^21 cm^-2. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral galaxy by a factor of \~7 . The aromatic features in Holmberg~II (which has only a slightly lower gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises that the strength of the aromatic features is not a simple linear function of metallicity. We estimate dust masses of ~10^4-10^6 M_sun for the M81 dwarf galaxies, resulting in an average dust--to--gas ratio (M_dust/M_HI) of ~3x10^-4 (1.5x10^-3 if only the HI that is associated with dust emission is considered); this is an order of magnitude lower than the typical value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation rate at 70um as compared to the more massive galaxies in SINGS by a factor of ~2. However, the average 70um/160um ratio in the sample dwarf galaxies is higher than what is found in the other galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in conjunction with a higher dust temperature in the dwarfs.

Infrared (IR) emission spectra are calculated for dust composed of mixtures of amorphous silicate and graphitic grains, including varying amounts of polycyclic aromatic hydrocarbon (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 Space Telescope. We find a size distribution for the PAHs that results in emission band ratios consistent with observed spectra of the Milky Way and other galaxies. Emission spectra are presented for various intensities of the illuminating starlight. We calculate how the efficiency of emission into different IR bands depends on PAH size; the strong 7.7um emission feature is produced mainly by PAH particles containing <1000 C atoms. We show how the emission spectrum depends on U, the starlight intensity relative to the local interstellar radiation field. The submm and far-infrared emission is compared to the observed emission from the local interstellar medium. Using a simple distribution function, we calculate the emission spectrum for dust heated by a distribution of starlight intensities, such as occurs within galaxies. The models are parameterized by the PAH mass fraction qpah, the lower cutoff Umin, and the fraction gamma of the dust heated by starlight with U>Umin. We present graphical procedures using IRAC and MIPS photometry to estimate qpah, Umin, and gamma, the fraction f_PDR of the dust luminosity coming from photodissociation regions with U>100, and the total dust mass Mdust.

DDSCAT 6.1 is a software package which applies the discrete dipole approximation (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. DDSCAT 6.1 allows accurate calculations of electromagnetic scattering from targets with size parameters 2 pi a_eff/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 2). DDSCAT 6.1 includes support for MPI and FFTW. We also make available a "plain" distribution of DDSCAT 6.1 that does not include support for MPI, FFTW, or netCDF, but is much simpler to install than the full distribution. The DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions.

We present a sample of low-resolution 5-38um 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 which indicate most normal galaxies undergo A_V < ~3 magnitudes 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_sun/4, as well as in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5x) 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-8um features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGN. The AGN 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 AGN 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 two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z=0-2.5.

Using data obtained for twelve galaxies as part of the {\it Spitzer} Infrared Nearby Galaxies Survey (SINGS) and the Westerbork Synthesis Radio Telescope (WSRT)-SINGS radio continuum survey, we study how star formation activity affects the far-infrared (FIR)--radio correlation {\it within} galaxies by testing a phenomenological model, which describes the radio image as a smeared version of the FIR image. The physical basis of this description is that cosmic-ray (CR) electrons will diffuse measurably farther than the mean free path of dust-heating photons before decaying by synchrotron radiation. This description works well in general. Galaxies with higher infrared surface brightnesses have best-fit smoothing scale-lengths of a few hundred parsecs, substantially shorter than those for lower surface brightness galaxies. We interpret this result to suggest that galaxies with higher disk averaged star formation rates have had a recent episode of enhanced star formation and are characterized by a higher fraction of young CR electrons that have traveled only a few hundred parsecs from their acceleration sites in supernova remnants compared to galaxies with lower star formation activity.

Analyses of spectra obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE) satellite, together with spectra from the Copernicus and IMAPS instruments, reveal an unexplained very wide range in the observed deuterium/hydrogen (D/H) ratios for interstellar gas in the Galactic disk beyond the Local Bubble. We argue that spatial variations in the depletion of deuterium onto dust grains can explain these local variations in the observed gas-phase D/H ratios. We present a variable deuterium depletion model that naturally explains the constant measured values of D/H inside the Local Bubble, the wide range of gas-phase D/H ratios observed in the intermediate regime (log N(H I} = 19.2-20.7), and the low gas-phase D/H ratios observed at larger hydrogen column densities. We consider empirical tests of the deuterium depletion hypothesis: (i) correlations of gas-phase D/H ratios with depletions of the refractory metals iron and silicon, and (ii) correlation with the molecular hydrogen rotational temperature. Both of these tests are consistent with deuterium depletion from the gas phase in cold, not recently shocked, regions of the ISM, and high gas-phase D/H ratios in gas that has been shocked or otherwise heated recently. We argue that the most representative value for the total (gas plus dust) D/H ratio within 1 kpc of the Sun is >=23.1 +/- 2.4 (1 sigma) parts per million (ppm). This ratio constrains Galactic chemical evolution models to have a very small deuterium astration factor, the ratio of primordial to total (D/H) ratio in the local region of the Galactic disk, which we estimate to be f_d <= 1.19 +/-0.16 (1 sigma) or <= 1.12 +/- 0.14 (1 sigma) depending on the adopted light element nuclear reaction rates.

We present Spitzer imaging of the metal-deficient (Z ~30% Z_sun) Local Group dwarf galaxy NGC 6822. On spatial scales of ~130 pc, we study the nature of IR, H alpha, HI, and radio continuum emission. Nebular emission strength correlates with IR surface brightness; however, roughly half of the IR emission is associated with diffuse regions not luminous at H alpha (as found in previous studies). The global ratio of dust to HI gas in the ISM, while uncertain at the factor of ~2 level, is ~25 times lower than the global values derived for spiral galaxies using similar modeling techniques; localized ratios of dust to HI gas are about a factor of five higher than the global value in NGC 6822. There are strong variations (factors of ~10) in the relative ratios of H alpha and IR flux throughout the central disk; the low dust content of NGC 6822 is likely responsible for the different H alpha/IR ratios compared to those found in more metal-rich environments. The H alpha and IR emission is associated with high-column density (> ~1E21 cm^-2) neutral gas. Increases in IR surface brightness appear to be affected by both increased radiation field strength and increased local gas density. Individual regions and the galaxy as a whole fall within the observed scatter of recent high-resolution studies of the radio-far IR correlation in nearby spiral galaxies; this is likely the result of depleted radio and far-IR emission strengths in the ISM of this dwarf galaxy.

We calculate the scattering of X-rays by interstellar dust, for a dust model that reproduces the observed wavelength-dependent extinction and polarization of starlight. On interstellar sightlines that produce appreciable starlight polarization, we predict that the dust-scattered X-ray halo around point sources will have measurable azimuthal asymmetry due to scattering by partially-aligned nonspherical grains. We calculate the expected halo asymmetry. X-ray halo asymmetry provides a new test of interstellar dust models.

We explore the nature of variations in dust emission within an individual galaxy using 3.6 - 160 micron Spitzer Space Telescope observations and 450 and 850 micron James Clerk Maxwell Telescope observations of the edge-on Sd spiral galaxy NGC 4631 with the goals of understanding the relation between polycyclic aromatic hydrocarbons (PAH) and dust emission, studying the variations in the colors of the dust emission, and searching for possible excess submillimeter emission compared to what is anticipated based on dust models applied to the mid- and far-infrared data. PAH emission at 8 micron is found to correlate best with hot dust emission at 24 micron on kiloparsec scales, although the relation breaks down on scales equal to hundreds of parsecs, possibly because of differences in the mean free paths between the photons that excite the PAHs and heat the dust and possibly because the PAHs are destroyed by the hard radiation fields in the centers of some star formation regions. The ratio of 8 micron PAH emission to 160 micron cool dust emission appears to be a function of radius. The 70/160 and 160/450 micron flux density ratios are remarkably constant even though the surface brightness varies by a factor of 25 in each wave band, which suggests that the emission is from dust heated by a nearly-uniform radiation field. Globally, we find an excess of 850 micron emission relative to what would be predicted by dust models. The 850 micron excess is highest in regions with low 160 micron surface brightness, although the strength and statistical significance of this result depends on the model fit to the data. We rule out variable emissivity functions or ~4 K dust as the possible origins of this 850 micron emission, but we do discuss the other possible mechanisms that could produce the emission.

The recent star formation (SF) in the early-type spiral galaxy M81 is characterized using imaging observations from the far-ultraviolet (UV) to the far-infrared (IR). We compare these data with models of the stellar, gas, and dust emission for sub-galactic regions. Our results suggest the existence of a diffuse dust emission not directly linked to the recent SF. We find a radial decrease of the dust temperature and dust mass density, and in the attenuation of the stellar light. The IR emission in M81 can be modeled with three components: 1) cold dust with a temperature <T_c>=18+-2 K, concentrated near the HII regions but also presenting a diffuse distribution; 2) warm dust with T_w=53+-7 K, directly linked with the HII regions; and 3) aromatic molecules, with diffuse morphology peaking around the HII regions. We derive several relationships to obtain total IR luminosities from IR monochromatic fluxes, and we compare five different star formation rate (SFR) estimators for HII regions in M81 and M51: the UV, Halpha, and three estimators based on Spitzer data. We find that the Halpha luminosity absorbed by dust correlates tightly with the 24 microns emission. The correlation with the total IR luminosity is not as good. Important variations from galaxy to galaxy are found when estimating the total SFR with the 24 microns or the total IR emission alone. The most reliable estimations of the total SFRs are obtained by combining the Halpha emission (or the UV) and an IR luminosity (especially the 24 microns emission), which probe the unobscured and obscured SF, respectively. For the entire M81 galaxy, about 50% of the total SF is obscured by dust. The percentage of obscured SF ranges from 60% in the inner regions of the galaxy to 30% in the outer zones.

We present Spitzer observations of the dwarf starburst galaxy NGC 1705 obtained as part of SINGS. The galaxy morphology is very different shortward and longward of ~5 microns: short-wavelength imaging shows an underlying red stellar population, with the central super star cluster (SSC) dominating the luminosity; longer-wavelength data reveals warm dust emission arising from two off-nuclear regions offset by ~250 pc from the SSC. These regions show little extinction at optical wavelengths. The galaxy has a relatively low global dust mass (~2E5 solar masses, implying a global dust-to-gas mass ratio ~2--4 times lower than the Milky Way average). The off-nuclear dust emission appears to be powered by photons from the same stellar population responsible for the excitation of the observed H Alpha emission; these photons are unassociated with the SSC (though a contribution from embedded sources to the IR luminosity of the off-nuclear regions cannot be ruled out). Low-resolution IRS spectroscopy shows moderate-strength PAH emission in the 11.3 micron band in the eastern peak; no PAH emission is detected in the SSC or the western dust emission complex. There is significant diffuse 8 micron emission after scaling and subtracting shorter wavelength data; the spatially variable PAH emission strengths revealed by the IRS data suggest caution in the interpretation of diffuse 8 micron emission as arising from PAH carriers alone. The metallicity of NGC 1705 falls at the transition level of 35% solar found by Engelbracht and collaborators; the fact that a system at this metallicity shows spatially variable PAH emission demonstrates the complexity of interpreting diffuse 8 micron emission. A radio continuum non-detection, NGC 1705 deviates significantly from the canonical far-IR vs. radio correlation. (Abridged)

We analyze extensive data on NGC1377 from the Spitzer Infrared Nearby Galaxies Survey (SINGS). Within the category of nascent starbursts, that we previously selected by their infrared to radio continuum ratios in large excess of the average and their hot dust, NGC1377 has the largest infrared excess yet measured. Optical imaging reveals a morphological distortion suggestive of a recent accretion event. Infrared spectroscopy reveals a compact and opaque source dominated by a hot, self-absorbed continuum (tau ~ 20 in the 10 micron silicate band). We provide physical evidence against non-stellar activity being the heating source. HII regions are detected through the single [NeII] line, probing <1% of the ionizing radiation. Not only is the optical depth very high, but >85% of ionizing photons are suppressed by dust. The only other detected emission features are molecular hydrogen lines, arguably excited mainly by shocks, besides photodissociation regions, and weak aromatic bands. The new observations support our interpretation in terms of an extremely young starburst (<1 Myr). More generally, galaxies deficient in radio synchrotron are likely observed within a few Myr of the onset of a starburst and after a long quiescence, prior to the replenishment of the ISM with cosmic rays. The similar infrared-radio properties of NGC1377 and some infrared-luminous galaxies suggest that NGC1377 constitutes an archetype to better understand starburst evolution. Although rare locally because observed in a brief evolutionary stage, nascent starbursts may represent a non-negligible fraction of merger-induced starbursts dominating deep infrared counts. Since they differ dramatically from usual starburst templates, they bear important consequences for the interpretation of deep surveys.

Mid-infrared diagnostics are presented for a large portion of the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample plus archival data from the Infrared Space Observatory and the Spitzer Space Telescope. Our main result is that these mid-infrared diagnostics effectively constrain a target's dominant power source. The combination of a high ionization line index and PAH strength serves as an efficient discriminant between AGN and star-forming nuclei, confirming progress made with ISO spectroscopy on starbursting and ultraluminous infrared galaxies. The sensitivity of Spitzer allows us to probe fainter nuclei and star-forming regions within galaxy disks. We find that both star-forming nuclei and extranuclear regions stand apart from nuclei that are powered by Seyfert or LINER activity. In fact, we identify areas within four diagnostic diagrams containing >90% Seyfert/LINER nuclei or >90% HII regions/HII nuclei. We also find that, compared to starbursting nuclei, extranuclear regions typically separate even further from AGN, especially for low-metallicity extranuclear environments. In addition, instead of the traditional mid-infrared approach to differentiating between AGN and star-forming sources that utilizes relatively weak high-ionization lines, we show that strong low-ionization cooling lines of X-ray dominated regions like [SiII] 34.82 micron can alternatively be used as excellent discrimants. Finally, the typical target in this sample shows relatively modest interstellar electron density (~400 cm^{-3}) and obscuration (A_V ~ 1.0 mag for a foreground screen), consistent with a lack of dense clumps of highly obscured gas and dust residing in the emitting regions. [Abridged]

We present new images (groundbased optical and mid-infrared (MIR) from the Spitzer Space Telescope) and spectra (from Spitzer) of the archetypal starburst galaxy M 82. The Spitzer data show that the MIR emission extends at least 6 kpc along the minor axis of the galaxy. We use the optical and infrared data to demonstrate that the extended emission is dominated by emission from dust. The colors of the MIR emission and the spectra indicate that there is a strong component of aromatic feature emission (the MIR features commonly attributed to polycyclic aromatic hydrocarbons). The dust continuum and aromatic feature emission are both strong in the well-known superwind region of this galaxy; clearly the carrier of the aromatic features can survive in close proximity to the wind, far from the plane of the galaxy. We also see significant emission by dust well outside the superwind region, providing the clearest picture to date of the dust distribution in the halo of this galaxy.

We present Spitzer 3.6-160 micron images, Spitzer mid-infrared spectra, and JCMT SCUBA 850 micron images of the Sombrero Galaxy (NGC 4594), an Sa galaxy with a 10^9 M_solar low luminosity active galactic nucleus (AGN). The brightest infrared sources in the galaxy are the nucleus and the dust ring. The spectral energy distribution of the AGN demonstrates that, while the environment around the AGN is a prominent source of mid-infrared emission, it is a relatively weak source of far-infrared emission, as had been inferred for AGN in previous research. The weak nuclear 160 micron emission and the negligible polycyclic aromatic hydrocarbon emission from the nucleus also implies that the nucleus is a site of only weak star formation activity and the nucleus contains relatively little cool interstellar gas needed to fuel such activity. We propose that this galaxy may be representative of a subset of low ionization nuclear emission region galaxies that are in a quiescent AGN phase because of the lack of gas needed to fuel circumnuclear star formation and Seyfert-like AGN activity. Surprisingly, the AGN is the predominant source of 850 micron emission. We examine the possible emission mechanisms that could give rise to the 850 micron emission and find that neither thermal dust emission, CO line emission, bremsstrahlung emission, nor the synchrotron emission observed at radio wavelengths can adequately explain the measured 850 micron flux density by themselves. The remaining possibilities for the source of the 850 micron emission include a combination of known emission mechanisms, synchrotron emission that is self-absorbed at wavelengths longer than 850 microns, or unidentified spectral lines in the 850 micron band.

Photoelectric emission from dust plays an important role in grain charging and gas heating. To date, detailed models of these processes have focused primarily on grains exposed to soft radiation fields. We provide new estimates of the photoelectric yield for neutral and charged carbonaceous and silicate grains, for photon energies exceeding 20 eV. We include the ejection of electrons from both the band structure of the material and the inner shells of the constituent atoms, as well as Auger and secondary electron emission. We apply the model to estimate gas heating rates in planetary nebulae and grain charges in the outflows of broad absorption line quasars. For these applications, secondary emission can be neglected; the combined effect of inner shell and Auger emission is small, though not always negligible. Finally, we investigate the survivability of dust entrained in quasar outflows. The lack of nuclear reddening in broad absorption line quasars may be explained by sputtering of grains in the outflows.

(Abridged) We present an initial look at the far infrared-radio correlation within the star-forming disks of four nearby, nearly face-on galaxies (NGC 2403, NGC 3031, NGC 5194, and NGC 6946). Using Spitzer MIPS imaging and WSRT radio continuum data, observed as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS), we are able to probe variations in the logarithmic 24mu/22cm (q_24) and 70mu/22cm (q_70) surface brightness ratios across each disk at sub-kpc scales. We find general trends of decreasing q_24 and q_70 with declining surface brightness and with increasing radius. The residual dispersion around the trend of q_24 and q_70 versus surface brightness is smaller than the residual dispersion around the trend of q_24 and q_70 versus radius, on average by ~0.1 dex, indicating that the distribution of star formation sites is more important in determining the infrared/radio disk appearance than the exponential profiles of disks. We have also performed preliminary phenomenological modeling of cosmic ray electron (CRe^-) diffusion using an image-smearing technique, and find that smoothing the infrared maps improves their correlation with the radio maps. Exponential kernels tend to work better than Gaussian kernels which suggests that additional processes besides simple random-walk diffusion in three dimensions must affect the evolution of CRe^-s. The best fit smoothing kernels for the two less active star-forming galaxies (NGC 2403 and NGC 3031) have much larger scale-lengths than those of the more active star-forming galaxies (NGC 5194 and NGC 6946). This difference may be due to the relative deficit of recent CRe^- injection into the interstellar medium (ISM) for the galaxies having largely quiescent disks.

We present spatially resolved Spitzer imaging of the supergiant shell region of the M81 group dwarf galaxy IC 2574 obtained as part of the Spitzer Infrared Nearby Galaxies Survey. This region harbors one of the best nearby examples of a kinematically distinct HI shell, with an associated remnant stellar cluster; the shell is initiating sequential star formation as it interacts with the surrounding interstellar medium. This region dominates the infrared luminosity of IC 2574 and is spatially resolved in all Spitzer imaging bands. We study the differences in dust temperature as a function of local environment and compare local star formation rates as inferred from H Alpha and total infrared luminosities. We find that the strong H Alpha sources are associated with regions of warm dust; however, the most luminous infrared and H Alpha sources are not necessarilyco-spatial. The coolest dust is found in the regions farthest from the rim of the shell; these regions show the best agreement between star formation rates derived from H Alpha and from total infrared luminosities (although discrepancies at the factor of 3-4 level still exist). There is considerable variation in the radio-far infrared correlation in different regions surrounding the shell. The low dust content of the region may influence the scatter seen in these relations; these data demonstrate that the expanding shell is dramatically affecting its surroundings by triggering star formation and altering the dust temperature.

The Spitzer Infrared Nearby Galaxies Survey (SINGS) is carrying out a comprehensive multi-wavelength survey on a sample of 75 nearby galaxies. The 1-850um spectral energy distributions are presented using broadband imaging data from Spitzer, 2MASS, ISO, IRAS, and SCUBA. The infrared colors derived from the globally-integrated Spitzer data are generally consistent with the previous generation of models that were developed based on global data for normal star-forming galaxies, though significant deviations are observed. Spitzer's excellent sensitivity and resolution also allow a detailed investigation of the infrared spectral energy distributions for various locations within the three large, nearby galaxies NGC3031 (M81), NGC5194 (M51), and NGC7331. Strong correlations exist between the local star formation rate and the infrared colors f_nu(70um)/f_nu(160um) and f_nu(24um)/f_nu(160um), suggesting that the 24 and 70um emission are useful tracers of the local star formation activity level. Preliminary evidence indicates that variations in the 24um emission, and not variations in the emission from polycyclic aromatic hydrocarbons at 8um, drive the variations in the f_nu(8.0um)/f_nu(24um) colors within NGC3031, NGC5194, and NGC7331. If the galaxy-to-galaxy variations in spectral energy distributions seen in our sample are representative of the range present at high redshift then extrapolations of total infrared luminosities and star formation rates from the observed 24um flux will be uncertain at the factor-of-five level (total range). The corresponding uncertainties using the redshifted 8.0um flux (e.g. observed 24um flux for a z=2 source) are factors of 10-20. Considerable caution should be used when interpreting such extrapolated infrared luminosities.

(Abridged) Far ultraviolet to far infrared images of the nearby galaxy NGC5194, from Spitzer, GALEX, Hubble Space Telescope and ground--based data, are used to investigate local and global star formation, and the impact of dust extinction in HII-emitting knots. In the IR/UV-UV color plane, the NGC5194 HII knots show the same trend observed for normal star-forming galaxies, having a much larger dispersion than starburst galaxies. We identify the dispersion as due to the UV emission predominantly tracing the evolved, non-ionizing stellar population, up to ages 50-100 Myr. While in starbursts the UV light traces the current SFR, in NGC5194 it traces a combination of current and recent-past SFR. Unlike the UV emission, the monochromatic 24 micron luminosity is an accurate local SFR tracer for the HII knots in NGC5194; this suggests that the 24 micron emission carriers are mainly heated by the young, ionizing stars. However, preliminary results show that the ratio of the 24 micron emission to the SFR varies by a factor of a few from galaxy to galaxy. While also correlated with star formation, the 8 micron emission is not directly proportional to the number of ionizing photons. This confirms earlier suggestions that the carriers of the 8 micron emission are heated by more than one mechanism.

Solid particles with the composition of interstellar dust and power-law size distribution dn/da propto a^{-p} for a < a_max with a_{max} > 3 lambda and 3 < p < 4 will have submm opacity spectral index beta(lambda) = dln(kappa)/dln(nu) approx (p-3) beta_{ism}, where beta_{ism} approx 1.7 is the opacity spectral index of interstellar dust material in the Rayleigh limit. For the power-law index p approx 3.5 that characterizes interstellar dust, and that appears likely for particles growing by agglomeration in protoplanetary disks, grain growth to sizes a > 3 mm will result in beta(1 mm) < ~1. Grain growth can naturally account for beta approx 1 observed for protoplanetary disks, provided that a_{max} > ~ 3 lambda.

Using the Texas Echelon Cross Echelle Spectrograph (TEXES) we mapped emission in the H_2 v = 0-0 S(1) and S(2) lines toward the Orion Bar PDR at 2" resolution. We also observed H_2 v = 0-0 S(4) at selected points toward the front of the PDR. Our maps cover a 12" by 40" region of the bar where H_2 ro-vibrational lines are bright. The distributions of H_2 0-0 S(1), 0-0 S(2), and 1-0 S(1) line emission agree in remarkable detail. The high spatial resolution (0.002 pc) of our observations allows us to probe the distribution of warm gas in the Orion Bar to a distance approaching the scale length for FUV photon absorption. We use these new observational results to set parameters for the PDR models described in a companion paper (Draine et al. 2005, in prep). The best-fit model can account for the separation of the H_2 emission from the ionization front and the intensities of the ground state rotational lines as well as the 1-0 S(1) and 2-1 S(1) lines. This model requires significant adjustments to the commonly used values for the dust UV attenuation cross section and the photoelectric heating rate.

We present an extinction analysis of 9 paths through the LMC and SMC based on FUSE observations. To date, just two LMC sight lines have probed dust grain composition and size distributions in the Clouds using spectra including wavelengths as short as 950 A. We supplement these with results from 4 regions distinguished by their IR through UV extinction curves and grouped as LMCAvg, LMC2, SMC bar and SMC wing. Despite the distinct characters of Milky Way and Magellanic Cloud extinction, our results are generally analogous to those found for Galactic curves in that the FUSE portions of each extinction curve are described reasonably well by FM curves fitted only to longer wavelength data and lack any dramatic new extinction features, and any deviations from the CCM formalism continue into FUV wavelengths. An MEM analysis of these curves suggests that LMCAvg and SMC wing sight lines require more silicon and/or carbon in dust than current abundance measurements would allow, while the requirements for LMC2 and SMC bar sight lines do not fully tax the available reservoirs. An intermediate product of this analysis is the measurement of new H_2 abundances in the Magellanic Clouds. Collectively considering Cloud sight lines that possess significant H_2 column densities, E(B-V)/N(HI) ratios are reduced by significant factors relative to the Galactic mean, whereas the corresponding E(B-V)/N(H_2) values more closely resemble their Galactic counterpart. These trends reflect the fact that among these sight lines f(H_2)-values are lower than those common in the Milky Way for paths with similar degrees of reddening.

We have used the Gemini Multi-object Spectrograph (GMOS) on the Gemini South 8m telescope in exceptional conditions (0.6" FWHM seeing) to observe the outer stellar disk of the Sculptor group galaxy NGC 300 at two locations. At our point source detection threshold of r' = 27.0 (3-sigma) mag, we trace the stellar disk out to a radius of 24', or 2.2 R_25 where R_25 is the 25 mag/arcsec**2 isophotal radius. This corresponds to about 10 scale lengths in this low-luminosity spiral (M_B = -18.6), or about 14.4 kpc at a cepheid distance of 2.0 +/- 0.07 Mpc. The background galaxy counts are derived in the outermost field, and these are within 10% of the mean survey counts from both Hubble Deep Fields. The luminosity profile is well described by a nucleus plus a simple exponential profile out to 10 optical scale lengths. We reach an effective surface brightness of 30.5 mag/arcsec**2 (2-sigma) at 55% completeness which doubles the known radial extent of the optical disk. These levels are exceedingly faint in the sense that the equivalent surface brightness in B or V is about 32 mag/arcsec**2. We find no evidence for truncation of the stellar disk. Only star counts can be used to reliably trace the disk to such faint levels, since surface photometry is ultimately limited by nonstellar sources of radiation. In the Appendix, we derive the expected surface brightness of one such source: dust scattering of starlight in the outer disk.

We present extinction curves that include data down to far ultraviolet wavelengths (FUV; 1050 - 1200 A) for nine Galactic sight lines. The FUV extinction was measured using data from the Far Ultraviolet Spectroscopic Explorer. The sight lines were chosen for their unusual extinction properties in the infrared through the ultraviolet; that they probe a wide range of dust environments is evidenced by the large spread in their measured ratios of total-to-selective extinction, R_V = 2.43 - 3.81. We find that extrapolation of the Fitzpatrick & Massa relationship from the ultraviolet appears to be a good predictor of the FUV extinction behavior. We find that predictions of the FUV extinction based upon the Cardelli, Clayton & Mathis (CCM) dependence on R_V give mixed results. For the seven extinction curves well represented by CCM in the infrared through ultraviolet, the FUV extinction is well predicted in three sight lines, over-predicted in two sight lines, and under-predicted in 2 sight lines. A Maximum Entropy Method analysis using a simple three component grain model shows that seven of the nine sight lines in the study require a larger fraction of grain materials to be in dust when FUV extinction is included in the models. Most of the added grain material is in the form of small (radii < 200 A) grains.

The D/H ratio in interstellar gas varies on scales of a few hundred pc in the local Milky Way, with D/H values ranging from ~7 ppm to ~22 ppm. The reduction in D/H relative to the primordial value of D/H (~ 26 ppm) is usually attributed to "astration" -- conversion of D into other elements by nuclear fusion in stars. However, it is shown here that astration has difficulty accounting for the observations because the expected associated variations in O/H are not seen. The lower D/H values are instead likely due to "depletion" of the D onto dust grains. Polycyclic aromatic hydrocarbons (PAHs) are a possible repository for the missing D, and it appears possible for gas-grain reactions to achieve extreme deuteration of carbonaceous grain material. Grain destruction will release D from the grains; the gas phase abundance of D should therefore be positively correlated with the gas phase abundances of other elements that exhibit strong depletions, such as Mg, Si, Ti, and Fe, which also will be returned to the gas by grain destruction.

We calculate the relative grain-grain motions arising from interstellar magnetohydrodynamic (MHD) turbulence. The MHD turbulence includes both fluid motions and magnetic fluctuations. While the fluid motions accelerate grains through hydro-drag, the electromagnetic fluctuations accelerate grains through resonant interactions. We consider both incompressive (Alfv\'{e}n) and compressive (fast and slow) MHD modes and use descriptions of MHD turbulence obtained in Cho & Lazarian (2002). Calculations of grain relative motion are made for realistic grain charging and interstellar turbulence that is consistent with the velocity dispersions observed in diffuse gas, including cutoff of the turbulence from various damping processes. We show that fast modes dominate grain acceleration, and can drive grains to supersonic velocities. Grains are also scattered by gyroresonance interactions, but the scattering is less important than acceleration for grains moving with sub-Alfv\'{e}nic velocities. Since the grains are preferentially accelerated with large pitch angles, the supersonic grains will be aligned with long axes perpendicular to the magnetic field. We compare grain velocities arising from MHD turbulence with those arising from photoelectric emission, radiation pressure and H$_{2}$ thrust. We show that for typical interstellar conditions turbulence should prevent these mechanisms from segregating small and large grains. Finally, gyroresonant acceleration is bound to preaccelerate grains that are further accelerated in shocks. Grain-grain collisions in the shock may then contribute to the overabundance of refractory elements in the composition of galactic cosmic rays.

We examine the feasibility of using dust-scattered X-rays for direct determination of distances to nearby galaxies with bright background AGNs, QSOs, or GRBs. We show how the Chandra X-Ray Observatory could be used to determine the distance to M31 to an unprecedented absolute accuracy of 1%.

IRAC images of M81 show three distinct morphological constituents: a smooth distribution of evolved stars with bulge, disk, and spiral arm components; a clumpy distribution of dust emission tracing the spiral arms; and a pointlike nuclear source. The bulge stellar colors are consistent with M-type giants, and the disk colors are consistent with a slightly younger population. The dust emission generally follows the blue and ultraviolet emission, but there are large areas that have dust emission without ultraviolet and smaller areas with ultraviolet but little dust emission. The former are presumably caused by extinction, and the latter may be due to cavities in the gas and dust created by supernova explosions. The nucleus appears fainter at 8 um than expected from ground-based 10 um observations made four years ago.

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

We present Multiband Imaging Photometer for Spitzer (MIPS) observations of M81 at 24, 70, and 160 micron. The grand design nature of M81 is clearly seen, showing two well resolved spiral arms containing numerous bright star forming regions. The MIPS images reveal a significant amount of cold dust associated with the spiral arms. We investigate the variation of the ultraviolet (UV), H-alpha, and infrared (IR) luminosities and star formation rate (SFR) indicators across the face of M81 using the MIPS images and archival UV and H-alpha images. For regions in M81, we find that UV and H-alpha SFRs (uncorrected for dust attenuation) are always lower than the IR SFR. The cause of this behavior is dust attenuation and/or using SFR calibrations appropriate for entire galaxies, not regions in galaxies. The characteristics of the dust attenuation for the regions indicate the dust grains and/or geometry are different from those in starburst galaxies. The behavior of the infrared-radio correlation in M81 is seen to vary from the global average, with variations correlated with the morphology of M81.

We discuss the propagation of electromagnetic waves on a rectangular lattice of polarizable point dipoles. For wavelengths long compared to the lattice spacing, we obtain the dispersion relation in terms of the lattice spacing and the dipole polarizabilities. We also obtain the polarizabilities required for the lattice to have the same dispersion relation as a continuum medium of given refractive index m.; our result differs slightly from previous work by Draine & Goodman (1993). Our new prescription can be used to assign dipole polarizabilities when the discrete dipole approximation is used to study scattering by finite targets. Results are shown for selected cases.

We consider the absorption and scattering of X-rays observed from the Galactic center. One objective is to characterize the intrinsic X-ray emission from the central black hole, SgrA*, in its quiescent and flaring states. We correct the fluxes observed by the Chandra and XMM telescopes for absorption and scattering, using realistic models for the dust grain size distribution, gas and dust abundances and spatial distributions. However, there are still large uncertainties. Since much of the intervening dust is close to SgrA*, the scattered halo of X-ray photons is very concentrated: its intensity can dominate the Chandra PSF inside 1", and so affects estimates of the point source flux. It also broadens the radial intensity profiles of Galactic center sources - observations of this effect can help constrain models of the spatial distribution of the dust. We estimate that the combined scattering halos from observed Galactic center sources within 10" of SgrA* contribute up to ~10% of the observed diffuse emission in this region. Dust-scattered photons suffer a time delay relative to the photons that arrive directly. For dust that is 100 pc towards us from the Galactic center, this delay is about 1000 s at angles of 1", and 100 ks at 10". We illustrate how the evolution of the scattering halo following X-ray flares from SgrA* or other sources can also help to constrain the dust's line of sight distribution. We discuss the implications of X-ray scattering halos for the intensity of diffuse emission that has been reported within a few arcseconds of SgrA*: in the most extreme, yet viable, model we consider, ~1/3 of it is due to dust scattering of an unresolved source. The remainder results from an extended source of emission.

In the interstellar medium of the Milky Way, certain elements -- e.g., Mg, Si, Al, Ca, Ti, Fe -- reside predominantly in interstellar dust grains. These grains absorb, scatter, and emit electromagnetic radiation, heat the interstellar medium by photoelectric emission, play a role in the ionization balance of the gas, and catalyze the formation of molecules, particularly H2. I review the state of our knowledge of the composition and sizes of interstellar grains, including what we can learn from spectral features, luminescence, scattering, infrared emission, and observed gas-phase depletions. The total grain volume in dust models which reproduce interstellar extinction is significantly greater than estimated from observed depletions. Dust grains might reduce the gas-phase D/H ratio, providing an alternative mechanism to explain observed variations in the gas-phase D/H ratio in the local interstellar medium. Transport in dust grains could cause elemental abundances in newly-formed stars to differ from interstellar abundances.

The discrete-dipole approximation (DDA) is a powerful method for calculating absorption and scattering by targets that have sizes smaller than or comparable to the wavelength of the incident radiation. We present a new prescription -- the Surface-Corrected Lattice Dispersion Relation (SCLDR) -- for assigning the dipole polarizabilities that takes into account both target geometry and finite wavelength. We test the SCLDR in DDA calculations using spherical and ellipsoidal targets and show that for a fixed number of dipoles, the SCLDR prescription results in increased accuracy in the calculated cross sections for absorption and scattering. We discuss extension of the SCLDR prescription to irregular targets.

DDSCAT.6.0 is a freely available software package (http://www.astro.princeton.edu/~draine/DDSCAT.6.0.html) which applies the "discrete dipole approximation" (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. DDSCAT.6.0 allows accurate calculations of electromagnetic scattering from targets with ``size parameters'' 2*pi*a/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 1). DDSCAT.6.0 includes the option of using the FFTW (Fastest Fourier Transform in the West) package. DDSCAT.6.0 also includes MPI support, permitting parallel calculations on multiprocessor systems. DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to ``import'' arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. This User Guide explains how to use DDSCAT.6.0 to carry out electromagnetic scattering calculations. CPU and memory requirements are described.

Scattering and absorption of X-rays by interstellar dust is calculated for a model consisting of carbonaceous grains and amorphous silicate grains. The calculations employ realistic dielectric functions with structure near X-ray absorption edges, with resulting features in absorption, scattering, and extinction. Differential scattering cross sections are calculated for energies between 0.3 and 10 keV. The median scattering angle is given as a function of energy, and simple but accurate approximations are found for the X-ray scattering properties of the dust mixture, as well as for the angular distribution of the scattered X-ray halo for dust with simple spatial distributions. Observational estimates of the X-ray scattering optical depth are compared to model predictions. Observations of X-ray halos to test interstellar dust grain models are best carried out using extragalactic point sources.

Scattering and absorption properties at optical and ultraviolet wavelengths are calculated for an interstellar dust model consisting of carbonaceous grains and amorphous silicate grains. Polarization as a function of scattering angle is calculated for selected wavelengths from the IR to the vacuum UV. The widely-used Henyey-Greenstein phase function provides a good approximation for the scattering phase function between ~0.4 and 1 micron, but fails to fit the calculated phase functions at shorter wavelengths and longer wavelengths. A new analytic phase function is presented. It is exact at long wavelengths, and provides a good fit to the numerically-calculated phase function for lambda > 0.27um. Observational determinations of the scattering albedo and <cos(theta)> show considerable disagreement, especially in the ultraviolet. Possible reasons for this are discussed.

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 multi-wavelength library of ancillary and complementary observations, including radio continuum, HI, CO, submillimeter, BVRIJHK, H-alpha, Paschen-alpha, 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.

This review surveys the observed properties of interstellar dust grains: the wavelength-dependent extinction of starlight, including absorption features, from UV to IR; optical luminescence; IR 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 PAHs is compared with observed emission and scattering. Some issues concerning evolution of interstellar dust are discussed.

Nine lectures reviewing the astrophysics of dust in interstellar clouds. Topics include: (1) Summary of observational evidence concerning interstellar dust: broadband extinction, scattering of starlight, polarization of starlight, spectroscopy of dust, IR and FIR emission, and depletions of grain-forming elements. (2) Optics of interstellar dust grains: dielectric functions of nonconducting and conducting materials, calculational techniques, formulae valid in the Rayleigh limit, Kramers-Kronig relations, microwave emission mechanisms, and X-ray scattering. (3) IR and FIR emission: heating of interstellar dust, including single-photon heating, and resulting IR emission spectrum. (4) Charging of dust grains: collisional charging, photoelectric emission, and resulting charge distribution functions. (5) Dynamics: gas drag, Lorentz force, forces due to anisotropic radiation, and resulting drift velocities. (6) Rotational dynamics: brownian rotation, suprathermal rotation, and effects of starlight torques. (7) Alignment of interstellar dust: observations and theories. (8) Evolution of the grain population: dust formation in outflows, grain growth in the ISM, photodesorption, and grain destruction in shock waves. (9) Effects of dust grains: photoelectric heating, H2 formation, ion recombination, coupling of gas to magnetic fields, and dust grains as indicators of magnetic field direction.

We use published ROSAT observations of the X-ray Nova V1974 Cygni 1992 to test a model for interstellar dust, consisting of a mixture of carbonaceous grains and silicate grains. The time-dependent X-ray emission from the nova is modelled, and X-ray scattering is calculated for a dust mixture with a realistic size distribution. Model results are compared with the scattered X-ray halos measured by ROSAT at 9 epochs. The observed X-ray halos are consistent with the halos calculated for the size distribution of Weingartner & Draine which reproduces the Milky Way extinction with R_V=3.1, provided that the reddening E(B-V) approx 0.20, consistent with E(B-V) = 0.19 inferred from the late-time Balmer decrement. The time delay of the scattered halo relative to the direct flux from the nova is clearly detected. Models with smoothly-distributed dust give good overall agreement with the observed scattering halo, but tend to produce somewhat more scattering than observed at 200-300", and insufficient scattering at 50-100". While an additional population of large grains can increase the scattered intensity at 50-100", this could also be achieved by having ~30% of the dust in a cloud at a distance from us equal to ~95% of the distance to the nova. Such a model also improves agreement with the data at larger angles, and illustrates the sensitivity of X-ray scattering halos to the location of the dust. The observations therefore do not require a population of micron-sized dust grains.

In the previous papers in this series, we found that radiative torques can play a major role in the alignment of grains with the interstellar magnetic field. Since the radiative torques can drive the grains to suprathermal rotational speeds, in previous work we made the simplifying assumption that the grain principal axis of greatest moment of inertia is always parallel to the grain angular momentum. This enabled us to describe many of the features of the grain dynamics. However, this assumption fails when the grains enter periods of thermal rotation, which occur naturally in the radiative torque alignment scenario. In the present paper, we relax this assumption and explore the consequences for the grain dynamics. We develop a treatment to follow the grain dynamics including thermal fluctuations and ``thermal flipping'', and show results for one illustrative example. By comparing with a treatment without thermal fluctuations, we see that inclusion of thermal fluctuations can lead to qualitative changes in the grain dynamics. In a future installment in this series, we will use the more complete dynamical treatment developed here to perform a systematic study of grain alignment by radiative torques.

DDSCAT.5a is a freely available software package which applies the "discrete dipole approximation" (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The DDA approximates the target by an array of polarizable points. DDSCAT.5a requires that these polarizable points be located on a cubic lattice. DDSCAT.5a10 allows accurate calculations of electromagnetic scattering from targets with "size parameters" 2 pi a/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 1). The DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. This User Guide explains how to use DDSCAT.5a10 to carry out EM scattering calculations. CPU and memory requirements are described.

The infrared (IR) emission from interstellar dust in the Small Magellanic Cloud (SMC) is modelled using a mixture of amorphous silicate and carbonaceous grains, including a population of polycyclic aromatic hydrocarbon (PAH) molecules. (1) It is shown that this dust model is able to reproduce the spectral energy distribution from near-IR to far-IR for the entire SMC Bar region, provided the PAH abundance in the SMC Bar region is very low. (2) The IR spectrum of the SMCB1#1 molecular cloud can also be reproduced by our dust model provided the PAH abundance is increased relative to the overall SMC Bar. The PAHs in SMCB1#1 incorporate ~3% of the SMC C abundance, compared to <0.4% in the SMC Bar. (3) The spectrum of SMCB1#1 is best reproduced if the PAH mixture has intrinsic IR band strengths which differ from the band strengths which best fit Milky Way PAH mixtures. The variation in the PAH IR band strengths would imply different PAH mixtures, presumably a consequence of differing metallicity or environmental conditions. Other possibilities such as super-hydrogenation of PAHs and softening of the starlight spectrum are also discussed.

One of the major challenges to identification of the 3.3, 6.2, 7.7, 8.6, and 11.3um interstellar IR emission bands with polycyclic aromatic hydrocarbon (PAH) molecules has been the recent detection of these bands in regions with little UV illumination, since small, neutral PAH molecules have little or no absorption at visible wavelengths and therefore require UV photons for excitation. We show here that our "astronomical" PAH model, incorporating the experimental result that the visual absorption edge shifts to longer wavelength upon ionization and/or as the PAH size increases, can closely reproduce the observed infrared emission bands of vdB 133, a UV-poor reflection nebula. It is also shown that single-photon heating of ``astronomical'' PAHs in reflection nebulae near stars as cool as T_eff=3000K can result in observable emission at 6.2, 7.7, 8.6, and 11.3um. Illustrative mid-IR emission spectra are also calculated for reflection nebulae illuminated by cool stars with T_eff=3500, 4500, 5000K. These will allow comparison with future Space Infrared Telescope Facility (SIRTF) observations of vdB 135 (T_eff=3600K), vdB 47 (T_eff=4500K), and vdB 101 (T_eff=5000K). It is also shown that the dependence of the 12um IRAS emission relative to the total far-IR emission on the effective temperature of the exciting star is consistent with the PAH model expectation for 3000K < T_eff < 30000K.

A gamma ray burst with strong optical-UV emission occuring in a molecular cloud will photodissociate H2, photoionize H2, H, and He, and destroy dust grains. We model these processes, including time-dependent radiative transfer in both continuum radiation and the resonance lines of H2. The UV will pump H2 into vibrationally-excited levels. We calculate the absorption spectrum imprinted on radiation from the GRB at various times. In addition to the strong absorption lines of H2(v=0) at lambda < 1110 Angstrom, due to cold ambient gas, we find that radiation reaching us from the GRB and its afterglow will show strong absorption lines due to vibrationally-excited H2 at 1110 < lambda < 1705 Angstrom. These absorption lines, if observed, would provide unequivocal evidence for association of the GRB with molecular gas. Low-resolution spectra will exhibit conspicuous features due to clustering of individual lines; a list of the strongest such absorption features is given for spectral resolution R approx 350 characteristic of the grism on the Swift UV-Optical Telescope.

Silicon nanoparticles (SNPs) have been proposed as the source of the observed "extended red emission" (ERE) from interstellar dust. We calculate the thermal emission expected from such particles in a reflection nebula such as NGC 2023, and in the diffuse interstellar medium (ISM). Pure neutral Si SNPs would emit at 16.4 micron (um), while Si/SiO2 SNPs (both neutral and charged) produce a feature at 20 um. Observational upper limits on the 16.4 um and 20 um features in NGC 2023 impose upper limits of < 1.5 ppm in pure Si SNPs, and or < 0.2 ppm in Si/SiO2 SNPs. The observed ERE intensity from NGC 2023 then gives a lower bound on the required photoluminescence efficiency eta. For foreground extinction A(0.68um) = 1.2, we find eta > 5% for Si SNPs, or eta > 24% for Si/SiO2 SNPs in NGC 2023. Measurement of the R band extinction toward the ERE-emitting region could strengthen these lower limits. The ERE emissivity of the diffuse ISM appears to require > 42% (> 33%) of solar Si abundance in Si/SiO2 (Si) SNPs. We predict IR emission spectra and show that DIRBE photometry appears to rule out such high abundances of free-flying SNPs in the diffuse ISM. If the ERE is due to SNPs, they must be either in clusters or attached to larger grains. Future observations by SIRTF will be even more sensitive to the presence of free-flying SNPs.

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 um in diameter. We assume that the carbonaceous grains have polycyclic aromatic hydrocarbon (PAH)-like properties 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" due to 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.3um. The model is compared with the observed emission at high Galactic latitudes as well as in the Galactic plane, as measured by COBE and IRTS. We calculate infrared emission spectra for our dust model heated by a range of starlight intensities, and we provide tabulated dust opacities (extended tables available at http://www.astro.princeton.edu/~draine/dust/dustmix.html)

With the high-resolution spectroscopy now available in the optical and satellite UV, it is possible to determine the neutral/ionized column density ratios for several different elements in a single cloud. Assuming ionization equilibrium for each element, one can make several independent determinations of the electron density. For the clouds for which such an analysis has been carried out, these different estimates disagree by large factors, suggesting that some process (or processes) besides photoionization and radiative recombination might play an important role in the ionization balance. One candidate process is collisions of ions with dust grains. Making use of recent work quantifying the abundances of polycyclic aromatic hydrocarbon molecules and other grains in the interstellar medium, as well as recent models for grain charging, we estimate the grain-assisted ion recombination rates for several astrophysically important elements. We find that these rates are comparable to the rates for radiative recombination for conditions typical of the cold neutral medium. Including grain-assisted ion recombination in the ionization equilibrium analysis leads to increased consistency in the various electron density estimates for the gas along the line of sight to 23 Orionis. However, not all of the discrepancies can be eliminated in this way; we speculate on some other processes that might play a role. We also note that grain-assisted recombination of H+ and He+ leads to significantly lower electron fractions than usually assumed for the cold neutral medium.

We model the polarized thermal dust emission from protostellar cores that are assembled by super-sonic turbulent flows in molecular clouds. Self-gravitating cores are selected from a three dimensional simulation of super-sonic and super-Alfvenic magneto-hydrodynamic (MHD) turbulence. The polarization is computed in two ways. In model A it is assumed that dust properties and grain alignment efficiency are uniform; in model B it is assumed that grains are not aligned at visual extinction larger than 3 mag. The main results of this work are: i) Values of the degree of polarization P between 1 and 10% are typical, despite the super-Alfvenic nature of the turbulence; ii) A steep decrease of P with increasing values of the sub-mm dust continuum intensity I is always found in self--gravitating cores selected from the MHD simulations, if grains are not aligned above a certain value of visual extinction (model B); iii) The same behavior is hard to reproduce if grains are aligned independently of visual extinction (model A); iv) The Chandrasekhar-Fermi formula, corrected by a factor f=0.4, provides an approximate estimate of the average magnetic field strength in the cores. Sub-mm dust continuum polarization maps of quiescent protostellar cores and Bok globules always show a decrease in P with increasing value of I consistent with the predictions of our model B. We therefore conclude that sub-mm polarization maps of quiescent cores do not map the magnetic field inside the cores at visual extinction larger than approximately 3 mag. There is no inconsistency between the results from optical and near-IR polarized absorption of background stars, and the observed polarization of sub-mm dust continuum from quiescent cores. In both cases, grains at large visual extinction appear to be virtually unaligned.

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.

It has been proposed that the Galaxy might contain a population of cold clouds in numbers sufficient to account for a substantial fraction of the total mass of the Galaxy. These clouds would have masses of the order of 10^{-3} Solar mass and sizes of the order of 10 AU. We consider here the lensing effects of such clouds on the light from background stars. A semianalytical formalism for calculation of the magnification event rate produced by such gaseous lensing is developed, taking into account the spatial distribution of the dark matter in the Galaxy, the velocity distribution of the lensing clouds and source stars, and motion of the observer. Event rates are calculated for the case of gaseous lensing of stars in the Large Magellanic Cloud and results are directly compared with the results of the search for gravitational microlensing events undertaken by the MACHO collaboration. The MACHO experiment strongly constrains the properties of the proposed molecular clouds, but does not completely rule them out. Future monitoring programs will either detect or more strongly constrain this proposed population.

The abundance of both amorphous and crystalline silicates in very small grains is limited by the fact that the 10 micron silicate emission feature is not detected in the diffuse ISM. On the basis of the observed IR emission spectrum for the diffuse ISM, the observed ultraviolet extinction curve, and the 10 micron silicate absorption profile, we obtain upper limits on the abundances of ultrasmall (a < 15 Angstrom) amorphous and crystalline silicate grains. Contrary to previous work, as much as ~20% of interstellar Si could be in a < 15 Angstrom silicate grains without violating observational constraints. Not more than ~5% of the Si can be in crystalline silicates (of any size).

We present a method for calculating the infrared emission from a population of dust grains heated by starlight, including very small grains for which stochastic heating by starlight photons results in high temperature transients. Because state-to-state transition rates are generally unavailable for complex molecules, we consider model PAH, graphitic, and silicate grains with realistic vibrational mode spectra and realistic radiative properties. The vibrational density of states is used in a statistical-mechanical description of the emission process. Unlike previous treatments, our approach fully incorporates multiphoton heating effects, important for large grains or strong radiation fields. We discuss how the "temperature" of the grain is related to its vibrational energy. By comparing with an "exact" statistical calculation of the emission process, we determine the conditions under which the "thermal" and the "continuous cooling" approximations can be used to calculate the emission spectrum. We present results for the infrared emission spectra of PAH grains of various sizes heated by starlight. We show how the relative strengths of the 6.2, 7.7, and 11.3um features depend on grain size, starlight spectrum and intensity, and grain charging conditions. We show results for grains in the "cold neutral medium", "warm ionized medium", and representative conditions in photodissociation regions. Our model results are compared to observed ratios of emission features for reflection nebulae and photodissociation regions, the Milky Way, normal spiral galaxies, and starburst galaxies.

Grains exposed to anisotropic radiation fields are subjected to forces due to the asymmetric photon-stimulated ejection of particles. These forces act in addition to the ``radiation pressure'' due to absorption and scattering. Here we model the forces due to photoelectron emission and the photodesorption of adatoms. The ``photoelectric'' force depends on the ambient conditions relevant to grain charging. We find that it is comparable to the radiation pressure when the grain potential is relatively low and the radiation spectrum is relatively hard. The calculation of the ``photodesorption'' force is highly uncertain, since the surface physics and chemsitry of grain materials are poorly understood at present. For our simple yet plausible model, the photodesorption force dominates the radiation pressure for grains with size >~0.1 micron exposed to starlight from OB stars. We find that the anisotropy of the interstellar radiation field is ~10% in the visible and ultraviolet. We estimate size-dependent drift speeds for grains in the cold and warm neutral media and find that micron-sized grains could potentially be moved across a diffuse cloud during its lifetime.

Photodissociation regions (PDRs) and shocks give rise to conspicuous emission from rotationally and vibrationally excited molecular hydrogen. This line emission has now been studied with ISO and from the ground in great detail. A remarkable discovery has been that toward the Orion outflow and other shock-excited regions, the H_2 level populations show a very high excitation component. We suggest that these high-excitation populations may arise from non-thermal pumping processes, such as H_2 formation and high-velocity ion-molecule collision in partially dissociative shocks. In PDRs such as NGC 7023 however, formation pumping is always less important than fluorescent pumping. We furthermore present two HD emission line detections toward Orion Peak 1. This enables the first comparison of the H_2 and the HD excitation, which surprisingly turn out to be identical.

We model the photoelectric emission from and charging of interstellar dust and obtain photoelectric gas heating efficiencies as a function of grain size and the relevant ambient conditions. Using realistic grain size distributions, we evaluate the net gas heating rate for various interstellar environments, and find less heating for dense regions characterized by R_V=5.5 than for diffuse regions with R_V=3.1. We provide fitting functions which reproduce our numerical results for photoelectric heating and recombination cooling for a wide range of interstellar conditions. In a separate paper we will examine the implications of these results for the thermal structure of the interstellar medium. Finally, we investigate the potential importance of photoelectric heating in H II regions, including the warm ionized medium. We find that photoelectric heating could be comparable to or exceed heating due to photoionization of H for high ratios of the radiation intensity to the gas density. We also find that photoelectric heating by dust can account for the observed variation of temperature with distance from the galactic midplane in the warm ionized medium.

Observations of H2 line emission have revealed higher-than-expected gas temperatures in a number of photodissociation fronts. We discuss the heating and cooling processes in photodissociation regions. Observations of NGC 2023 are compared to a theoretical model in which there is substantial gas at temperatures T = 500-1000K heated by photoelectric emission and collisional deexcitation of H2. In general the model successfully reproduces the observed H2 line emission from a wide range of energy levels. The observed [SiII]34.8um emission appears to indicate substantial depletion of Si in NGC 2023.

Our current understanding of interstellar dust is summarized at an introductory level. Submicron-sized interstellar dust grains absorb and scatter light, and reradiate the absorbed energy in the infrared. The grain population spans a range of sizes, from molecules containing only tens of atoms, to particles containing 10^{10} atoms. Most of the grain mass appears to be due to two types of solid, in approximately equal amounts: amorphous silicate mineral, and carbonaceous material. Approximately 2/3 of the interstellar carbon in diffuse clouds is in solid form. Interstellar grains are important for many reasons: they obscure our view (in a wavelength-dependent manner); the infrared emission from dust provides a valuable probe of dense regions; and because dust grains play a central role in interstellar chemistry, gas dynamics, and heating and cooling of interstellar gas.

We show that the energy-level splitting arising from grain rotation ensures that paramagnetic dissipation acts at its maximum rate, i.e., the conditions for paramagnetic resonance are automatically fulfilled. We refer to this process as ``resonance relaxation''. The differences between the predictions of classical Davis-Greenstein relaxation and resonance relaxation are most pronounced for grains rotating faster than 1 GHz, i.e., in the domain where classical paramagnetic relaxation is suppressed. This mechanism can partially align even very small grains, resulting in linearly polarized microwave emission which could interfere with efforts to measure the polarization of the cosmic microwave background.

The prompt optical flash recently detected accompanying GRB990123 suggests that, for at least some GRBs, gamma-ray emission is accompanied by prompt optical-UV emission with luminosity L(1-7.5eV)=10^{49}(\Delta\Omega/4\pi)erg/s, where \Delta\Omega is the solid angle into which gamma-ray and optical-UV emission is beamed. Such an optical-UV flash can destroy dust in the beam by sublimation out to an appreciable distance, approximately 10 pc, and may clear the dust out of as much as 10^7(\Delta\Omega/4\pi)M_sun of molecular cloud material on an apparent time scale of 10 seconds. Detection of time dependent extinction on this time scale would therefore provide strong constraints on the GRB source environment. Dust destruction implies that existing, or future, observations of not-heavily-reddened fireballs are not inconsistent with GRBs being associated with star forming regions. In this case, however, if gamma-ray emission is highly beamed, the expanding fireball would become reddened on a 1 week time scale. If the optical depth due to dust beyond approximately 8 pc from the GRB is 0.2<\tau_V<2, most of the UV flash energy is converted to infra-red, \lambda \sim 1 micron, radiation with luminosity \sim 10^{41} erg/s extending over an apparent duration of \sim 20(1+z)(\Delta\Omega/0.01) day. Dust infra-red emission may already have been observed in GRB970228 and GRB980326, and may possibly explain their unusual late time behavior.

If a gamma ray burst with strong UV emission occurs in a molecular cloud, there will be observable consequences resulting from excitation of the surrounding H2. The UV pulse from the GRB will pump H2 into vibrationally-excited levels which produce strong absorption at wavelengths < 1650 A. As a result, both the prompt flash and later afterglow will exhibit strong absorption shortward of 1650 A, with specific spectroscopic features. Such a cutoff in the emission from GRB 980329 may already have been observed by Fruchter et al.; if so, GRB 980329 was at redshift 3.0 < z < 4.4 . BVRI photometry of GRB 990510 could also be explained by H2 absorption if GRB 990510 is at redshift 1.6 < z < 2.3. The fluorescence accompanying the UV pumping of the H2 will result in UV emission from the GRB which can extend over days or months, depending on parameters of the ambient medium and beaming of the GRB flash. The 7.5-13.6 eV fluorescent luminosity is \sim 10^{41.7} erg/s for standard estimates of the parameters of the GRB and the ambient medium. Spectroscopy can distinguish this fluorescent emission from other possible sources of transient optical emission, such as a supernova.

We identify a new mechanism of internal dissipation of rotational kinetic energy in spinning dust grains, arising from the reorientation of nuclear angular momentum, e.g., spins of protons. Grain rotation induces magnetization of the nuclear spin system, with net alignment of nuclear spins parallel to the grain angular velocity. When the grain does not rotate around a principal axis, the nuclear magnetization vector precesses in grain body coordinates, resulting in dissipation of energy. The analogous process involving electron spins was discovered by Purcell and termed ``Barnett relaxation''. We revisit the physics of the ``Barnett relaxation'' process and correct the estimate for the Barnett relaxation rate. We show that nuclear relaxation can be orders of magnitude more important than Barnett relaxation. This finding implies that the processes of ``thermal flipping'' and ``thermal trapping'' are important for a broad range of grain sizes.

ISO has provided us with a new perspective on gas temperatures in photodissociation regions through measurements of line emission from rotationally-excited levels of the H_2 molecule. The H_2 rotational level populations provide a thermal probe, showing that gas temperatures T_gas approx 500-1000 K prevail in a portion of the PDR where significant quantities of H_2 are present. Such high gas temperatures were unexpected. Theoretical models for the S140 PDR are presented. Possible mechanisms for heating the gas to such high temperatures are discussed.

Observations of the cosmic microwave background have revealed a component of 10-60 GHz emission from the Galaxy which correlates with 100-140um emission from interstellar dust but has an intensity much greater than expected for the low-frequency tail of the "electric dipole vibrational" emission peaking at \~130um. This "anomalous emission" is more than can be accounted for by dust-correlated free-free emission. The anomalous emission could be due in part to magnetic dipole emission from thermal fluctuations of the magnetization within interstellar dust grains, but only if a substantial fraction of the Fe in interstellar dust resides in magnetic materials such as metallic iron or magnetite. The observed anomalous emission is probably due primarily to electric dipole radiation from spinning ultrasmall interstellar dust grains. This rotational emission is expected to be partially polarized. From the standpoint of minimizing confusion with non-CBR foregrounds, 60-120 GHz appears to be the optimal frequency window.

Since the classical work by Purcell (1979) it has been generally accepted that most interstellar grains rotate suprathermally. Suprathermally rotating grains would be nearly perfectly aligned with the magnetic field by paramagnetic dissipation if not for ``crossovers'', intervals of low angular velocity resulting from reversals of the torques responsible for suprathermal rotation; during crossovers grains are susceptible to disalignment by random impulses. Lazarian and Draine (1997) identified thermal fluctuations within grain material as an important component of crossover dynamics. For grains of size less than 0.1 micron, these fluctuations ensure good correlation of angular momentum before and after crossover resulting in good alignment, in accord with observations of starlight polarization. In the present paper we discuss two new processes which are important for the dynamics of grains with a<0.1 micron. The first -- ``thermal flipping'' -- offers a way for small grains to bypass the period of greatly reduced angular momentum which would otherwise take place during a crossover, thereby enhancing the alignment of small grains. The second effect -- ``thermal trapping'' -- arises when thermal flipping becomes rapid enough to prevent the systematic torques from driving the grain to suprathermal rotation. This effect acts to reduce the alignment of small grains. The observed variation of grain alignment with grain size would then result from a combination of the thermal flipping process -- which suppresses suprathermal rotation of small grains -- and due to molecular hydrogen formation and starlight -- which drive large grains to suprathermal rotation rates.

We consider the depletion of elements from the interstellar gas onto a population of very small dust grains. Adopting a grain model in which of order 4% of the cosmic C abundance is in grains with radii <= 10 Angstroms, we find that the rate of accretion onto these grains is adequately fast to account for the observed large depletions of elements like Ti, without invoking unreasonably high rates of mass transfer between interstellar phases or low grain destruction rates. If these grains are composed of arene rings, then only a limited number of metal atoms can be locked up in them. The depletion would be quenched when this limit is reached, unless there is a mechanism for transferring the metals to larger grains and refreshing the very small grain population, for example by grain coagulation and shattering in the diffuse ISM. If Fe depletes onto the very small grains, then for reasonable coagulation rates there is at least one metal atom per five C atoms in the very small grain population. Furthermore, approximately 60% of the cosmic Fe is associated with the carbonaceous grain population. It is unclear whether this scenario is compatible with observations. However, if there is another population of very small grains, with a large capacity for holding Fe atoms, it might be the sink for the most heavily depleted elements.

We discuss the forces on grains exposed to anisotropic radiation fields, including the usual ``radiation pressure'' force and also recoil forces due to photoemitted electrons, photodesorbed hydrogen atoms, and hydrogen molecules which form on the grain surface. We show that these forces can lead to grain dynamics in photodissociation regions which result in enhanced dust-to-gas ratios. Since the gas heating is probably dominated by photoelectric emission from dust, this might explain the unusually high gas temperatures inferred from ISO observations of molecular hydrogen in photodissociation regions.

This paper presents new near-infrared observations of the planetary nebula NGC 2346. The data include a broad K-band image, an image in the H2 vibrationally excited 1-0S(1) line and K band slit spectra at three positions in the nebula. In the H2 1-0S(1) line, the nebula is characterized by a central, bright torus, surrounded by weaker emission with a typical butterfly shape, as seen in Halpha and CO lines. The K band spectra show 11 H2 lines with excitation energies from 6150 to 12552 K. The H2 data have been compared to the predictions of models which follow the evolution with time of the H2 emission in PNe of different core mass and shell properties (Natta & Hollenbach 1998). These models compute the emission originating in the photodissociation region (PDR) created at the inner edge of the neutral shell by the UV radiation of the central core, as well as the emission in the shock associated to the expansion of the shell inside the precursor red-giant wind. In NGC 2346, a PDR origin of the H2 emission in a low-density molecular shell (n<~10^4 cm^-3) is indicated. At these low densities, time-dependent H2 chemistry and X-ray heating of the neutral gas enhance the predicted PDR H2 line intensity by large factors.

If the Galaxy contains ~10^{11}M_sol in cold gas clouds of ~Jovian mass and \~AU size, these clouds will act as converging lenses for optical light, magnifying background stars at a detectable rate. The resulting light curves can resemble those due to gravitational lensing by a point mass, raising the possibility that some of the events attributed to gravitational microlensing might in fact be due to ``gaseous lensing''. During a lensing event, the lens would impose narrow infrared and far-red H_2 absorption lines on the stellar spectrum. Existing programs to observe gravitational microlensing, supplemented by spectroscopy, can therefore be used to either detect such events or place limits on the number of such gas clouds present in the Galaxy.

Thermal fluctuations in the magnetization of interstellar grains will produce magnetic dipole emission at frequencies below ~100 GHz. We show how to calculate absorption and emission from small particles composed of magnetic materials. The Kramers-Kronig relations for a dusty medium are generalized to include the possibility of magnetic grains. The frequency-dependent magnetic permeability is discussed for candidate grain materials, including iron and magnetite. We calculate emission spectra for various interstellar grain candidates. While paramagnetic grains or magnetite grains cannot account for the observed "anomalous" emission from dust in the 14-90 GHz range, stronger magnetic dipole emission will result if a fraction of the grain material is ferromagnetic, as could be the case given the high Fe content of interstellar dust. The observed emission from dust near 90 GHz implies that not more than 5% of interstellar Fe is in the form of metallic iron grains or inclusions (e.g., in "GEMS"). However, we show that if most interstellar Fe is in a moderately ferromagnetic material, it could contribute a substantial fraction of the observed 14-90 GHz emission, perhaps comparable to the contribution from spinning ultrasmall dust grains. The two emission mechanisms can be distinguished by measuring the emission from dark clouds. The expected polarization of magnetic dipole emission is discussed

We derive an analytic expression for the intensity of resonance-line radiation ``trapped'' in a semi-infinite medium. Given a source function and destruction probability per scattering, the radiation pressure due to trapped photons can be calculated by numerically integrating over analytic functions. We apply this formalism to a plane-parallel model stellar atmosphere to calculate the radiation pressure due to Lyman-alpha photons produced following absorption of UV and X-rays from an AGN. For low surface gravity stars near the AGN (g~10 cm/sec^2, r~0.25 pc), we find that the pressure due to Lyman-alpha photons becomes an appreciable fraction of that required for hydrostatic support. If the broad emission line emitting gas in AGNs and QSOs consists of stellar outflows, it may be driven, in part, by Lyman-alpha pressure.

We discuss the rotational excitation of small interstellar grains and the resulting electric dipole radiation from spinning dust. Attention is given to excitation and damping of rotation by: collisions with neutrals; collisions with ions; plasma drag; emission of infrared radiation; emission of microwave radiation; photoelectric emission; and formation of H_2 on the grain surface. We introduce dimensionless functions F and G which allow direct comparison of the contributions of different mechanisms to rotational drag and excitation. Emissivities are estimated for dust in different phases of the interstellar medium, including diffuse HI, warm HI, low-density photoionized gas, and cold molecular gas. Spinning dust grains can explain much, and perhaps all, of the 14-50 GHz background component recently observed in CBR studies. It should be possible to detect rotational emission from small grains by ground-based observations of molecular clouds.

Spinning interstellar dust grains produce detectable rotational emission in the 10-100 GHz frequency range. We calculate the emission spectrum, and show that this emission can account for the ``anomalous'' Galactic background component which correlates with 100um thermal emission from dust. Implications for cosmic background studies are discussed.

We discuss the dynamics of dust grains subjected to uncompensated torques arising from H_2 formation. In particular, we discuss grain dynamics when a grain spins down and goes through a ``crossover''. As first pointed out by Spitzer & McGlynn (1979), the grain angular momentum before and after a crossover event are correlated, and the degree of this correlation affects the alignment of dust grains by paramagnetic dissipation. We calculate the correlation including the important effects of thermal fluctuations within the grain material. These fluctuations limit the degree to which the grain angular momentum J is coupled with the grain principal axis a_1 of maximal inertia. We show that this imperfect coupling of a_1 with J plays a critical role during crossovers and can substantially increase the efficiency of paramagnetic alignment for grains larger than 10^{-5} cm. As a result, we show that for reasonable choices of parameters, the observed alignment of $a\gtrsim 10^{-5}$ cm grains could be produced by paramagnetic dissipation in suprathermally rotating grains, if radiative torques due to starlight were not present. We also show that the efficiency of mechanical alignment in the limit of long alignment times is not altered by the thermal fluctuations in the grain material.

Radiative torques on irregular dust grains, in addition to producing superthermal rotation, play a direct dynamical role in the alignment of interstellar dust with the local magnetic field. The equations governing the orientation of spinning, precessing grains are derived; H_2 formation torques and paramagnetic dissipation are included in the dynamics. Stationary solutions (constant alignment angle and spin rate) are found; these solutions may be stable ("attractors") or unstable ("repellors"). The equations of motion are numerically integrated for three exemplary irregular grain geometries, exposed to anisotropic radiation with the spectrum of interstellar starlight. The resulting "trajectory maps" are classified as "noncyclic", "semicyclic", or "cyclic", with examples of each given. We find that radiative torques result in rapid grain alignment, even in the absence of paramagnetic dissipation. It appears that radiative torques due to starlight can account for the observed alignment of interstellar grains with the Galactic magnetic field.

Irregular dust grains are subject to radiative torques when irradiated by interstellar starlight. It is shown how these radiative torques may be calculated using the discrete dipole approximation. Calculations are carried out for one irregular grain geometry, and three different grain sizes. It is shown that radiative torques can play an important dynamical role in spinup of interstellar dust grains, resulting in rotation rates which may exceed even those expected from H_2 formation on the grain surface. Because the radiative torque on an interstellar grain is determined by the overall grain geometry rather than merely the state of the grain surface, the resulting superthermal rotation is expected to be long-lived. By itself, long-lived superthermal rotation would permit grain alignment by normal paramagnetic dissipation on the "Davis-Greenstein" timescale. However, radiative torques arising from anisotropy of the starlight background can act directly to alter the grain alignment on much shorter timescales, and are therefore central to the process of interstellar grain alignment. Radiative torques depend strongly on the grain size, measured by a_eff, the radius of a sphere of equal volume. In diffuse clouds, radiative torques dominate the torques due to H2 formation for a_eff=0.2micron grains, but are relatively unimportant for a_eff<0.05micron grains. This may provide a natural explanation for the observation that $a_eff>0.1 micron grains in diffuse clouds are aligned, while there is little alignment of a_eff < 0.05 micron grains. We show that radiative torques are ineffective at producing superthermal rotation within quiescent dark clouds, but can be very effective in star-forming regions such as the M17 molecular cloud.

The optical and magnetic properties of dust grains are reviewed, as they relate to the problem of interstellar grain alignment. Grain geometry plays an important role in determining the optical properties, and scattering and absorption of starlight will produce radiative torques which may drive grains to suprathermal rates of rotation in interstellar clouds; these radiative torques appear likely to play an active role in the alignment process. The likely magnetic properties of grains are discussed, with particular attention to the imaginary part of the magnetic susceptibility.

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<3000 and n_H<10^4 cm^{-3}; for stronger radiation fields and higher densities radiative and collisional depopulation of vibrationally-excited levels interferes with the radiative cascade. The emission spectrum from the PDR is essentially independent of the color temperature $T_{color}$ of the incident UV for T_{color}>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 $v$=1 and $v$=2 levels are computed as functions of the temperature and density, for different values of chi and n_H. We apply our models to the reflection nebula NGC 2023. We are best able to reproduce the observations with models having chi=5000, n_H=10^5 cm^{-3}.

Interstellar shock waves can erode and destroy grains present in the shocked gas, primarily as the result of sputtering and grain-grain collisions. Uncertainties in current estimates of sputtering yields are reviewed. Results are presented for the simple case of sputtering of fast grains being stopped in cold gas. An upper limit is derived for sputtering of refractory grains in C-type MHD shocks: shock speeds $v_s \gtrsim 50 \kms$ are required for return of more than 30\% of the silicate to the gas phase. Sputtering can also be important for removing molecular ice mantles from grains in two-fluid MHD shock waves in molecular gas. Recent estimates of refractory grain lifetimes against destruction in shock waves are summarized, and the implications of these short lifetimes are discussed.

We discuss the theory of coupled ionization--dissociation fronts produced when molecular clouds are exposed to $\lambda < 1110$\AA\ radiation from hot stars. A steady, composite structure is developed, which generally includes an ionized outflow away from the cloud, an ionization front, a layer of photodissociated gas, a photodissociation front, and a shock wave preceding the photodissociation front. We show that the properties of the structure are determined by two dimensionless parameters, $\psi$ and $\delta$, and by the Alfv\'en speed in the preshock gas. For a broad range of parameters of interest, the ionization front and the hydrogen photodissociation front do not separate, the H$_2$ photodissociation and photoionization take place together, and a classical hydrogen ``photodissociation region'' (PDR) does not exist. We also show that even when a distinct photodissociation region exists, in many cases the dissociation front propagates too rapidly for the usual stationary models of PDRs to be applicable. We discuss several famous PDRs, e.g., in M17 and Orion and conclude that they cannot be described by equilibrium PDR models.