Search

We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5-18 $\mu$m as metallicity decreases continuously from solar ($Z_\odot$) to 0.2 $Z_\odot$. The total PAH to dust luminosity ratio remains relatively constant until reaching a threshold of $\sim$$\frac{2}{3}$$Z_\odot$, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest wavelength PAH features is expected. In contrast, we find a steep decline in long wavelength power below $Z_\odot$, especially in the 17 $\mu$m feature, with the shorter wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies., Comment: 23 pages, 13 figures, accepted for publication in ApJ

We present the first mid-IR detection of the linear polarization toward the star CygOB2-12, a luminous blue hypergiant that, with Av of ~10 mag of foreground extinction, is a benchmark in the study of the properties of dust in the diffuse interstellar medium. The 8-13 micrometer spectropolarimetry, obtained with the CanariCam multi-mode camera at the Gran Telescopio Canarias (GTC), shows clear trends with wavelength characteristic of silicate grains aligned in the interstellar magnetic field. The maximum polarization, detected with 7.8 statistical significance near 10.2 micrometers, is (1.24 +/- 0.28) % with position angle (126 +/- 8) deg. We comment on these measurements in the context of recent models for the dust composition in the diffuse interstellar medium., Comment: Accepted Astrophysical Journal Letters, October 9, 2022

We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 A, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-two increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemisphere. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 A, as evidence for luminescence in the form of extended red emission (ERE). The observed strength of the 4000 A break indicates that at most ~7% of the dust-correlated light at 4000 A can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence independent of measurements of extinction., Comment: 16 pages, 8 figures. Revised version after reviewer comments

We present a comparative study of four physical dust models and two single-temperature modified blackbody models by fitting them to the resolved WISE, Spitzer, and Herschel photometry of M101 (NGC 5457). Using identical data and a grid-based fitting technique, we compare the resulting dust and radiation field properties derived from the models. We find that the dust mass yielded by the different models can vary by up to factor of 3 (factor of 1.4 between physical models only), although the fits have similar quality. Despite differences in their definition of the carriers of the mid-IR aromatic features, all physical models show the same spatial variations for the abundance of that grain population. Using the well determined metallicity gradient in M101 and resolved gas maps, we calculate an approximate upper limit on the dust mass as a function of radius. All physical dust models are found to exceed this maximum estimate over some range of galactocentric radii. We show that renormalizing the models to match the same Milky Way high latitude cirrus spectrum and abundance constraints can reduce the dust mass differences between models and bring the total dust mass below the maximum estimate at all radii., Comment: Accepted for publication in ApJ -- 28 pp (12 figures) + appendices

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., Comment: 15 pages, 10 figures, published in ApJ

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., Comment: 19 pages, 8 figures, 8 tables

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., Comment: Astro2020 decadal survey science white paper submitted to the National Academy of Sciences

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., Comment: Astro2020 decadal survey science white paper submitted to the National Academy of Sciences

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., Comment: Accepted for publication in ApJ

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., Comment: 8 pages, accepted for publication in ApJ

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., Comment: 19 pages, 16 figures, accepted for publication in ApJ

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., Comment: 17 pages, 5 figures, accepted for publication in ApJ

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

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., Comment: 32 pages, 15 figures

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 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., Comment: Accepted for publication in the Astrophysical Journal, 2015 January 12: 13 pages, 8 figures, 4 tables

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., Comment: 11 pages, 4 pages, to be published in Proceedings of Science: "Advancing Astrophysics with the Square Kilometre Array"

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., Comment: Final version with a corrected large-scale fit

'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., Comment: 13 pages, 4 figures, accepted for publication in Nature

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., Comment: 14 pages, 11 figures, accepted in ApJ

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., Comment: 95 pages. Expanded and updated to conform to latest release, DDSCAT 7.2.1. Complete software package available at http://code.google.com/p/ddscat/

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%., Comment: 14 pages, 12 figures, ApJ accepted with modest corrections and clarifications. Replaced with accepted version

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., Comment: accepted for publication in ApJ

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 (=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., Comment: accepted for publication in ApJ, abstract abridged

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., Comment: MNRAS, in press

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., Comment: 33 pages, 4 figures, 4 tables, accepted for publication in AJ

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]., Comment: accepted for publication in the Astrophysical Journal

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., Comment: 10 pages (emulateapj), 8 figures, 1 table. Submitted to ApJ

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., Comment: 9 pages, 4 figures, 1 table, submitted to ApJ

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., Comment: 25 pages, 7 figures, 2 tables, accepted to ApJ

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., Comment: 5 pages, ApJL accepted

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., Comment: 17 pages, 13 figures; submitted to ApJ

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., Comment: 12 pages, 8 figures, 1 table

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., Comment: ApJ accepted

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., Comment: 46 pages, 10 figures, accepted for publication in The Astrophysical Journal

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., Comment: 19 pages, 9 figures

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., Comment: ApJ, in press; please retrieve full-resolution version from http://www.astro.wesleyan.edu/~cannon/pubs.html

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., Comment: Accepted to the Astrophysical Journal

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 =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., Comment: Accepted for publication in ApJ, 22 pages, 11 figures, 4 tables, the complete data tables can be found in the source

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), Comment: ApJ, in press; please retrieve full-resolution version from http://www.astro.wesleyan.edu/~cannon/pubs.html

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., Comment: ApJ Letters, in press; please retrieve full-resolution version from http://www.mpia.de/homes/cannon/pubs.html

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., Comment: ApJ accepted -- 30 pages, 13 figures -- see ftp://www.aao.gov.au/pub/local/jbh/astro-ph/N300 for full resolution figures and preprint

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., Comment: 59 pages. Updated to describe new target option CYLCAP, to correct error in formula for linear polarization P of scattered light, and provide more explicit instructions for use of multiple dielectric functions

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., Comment: 21 pages, color figures, submitted to ApJ

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., Comment: 36 pages, 9 figures, 1 table, typos corrected, minor changes

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., Comment: LaTeX (4 pages, 4 eps figures, uses iso98.sty and twocolumn.sty macros); to appear in "The Universe as Seen by ISO", Eds. P. Cox and M.F. Kessler, ESA Special Publications Series (SP-427)

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., Comment: LaTeX (21 pages, 4 eps figures, uses aaspp4.sty); ApJ, accepted