New HST/STIS optical spectra were obtained for a sample of early type stars with existing IUE UV spectra. These data were used to construct optical extinction curves whose general properties are discussed elsewhere. In this paper, we identify extinction features in the curves that are wider than diffuse interstellar bands (DIBs) but narrower than the well known broad band variability. This intermediate scale structure, or ISS, contains distinct features whose peaks can contribute a few percent to 20% of the total extinction. Most of the ISS variation can be captured by three principal components. We model the ISS with three Drude profiles and show that their strengths and widths vary from one sight line to another, but their central positions are stable, near 4370, 4870 and 6300A. The Very Broad Structure, VBS, in optical curves appears to be a minimum between the 4870 and 6300A absorption peaks. We find relations among the fit parameters and provide a physical interpretation of them in terms of a simplistic grain model. Finally, we note that the strengths of the 4370 and 4870A features are correlated to the strength of the 2175A UV bump, but that the 6300A feature is not, and that none of the ISS features are related to R(V). However, we verify that the broad band curvature of the continuous optical extinction is strongly related to R(V).

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.

Photodissociation regions (PDRs) are parts of the ISM consisting of predominantly neutral gas, located at the interface between H II regions and molecular clouds. The physical conditions within these regions show variations on very short spatial scales, and therefore PDRs constitute ideal laboratories for investigating the properties and evolution of dust grains. We have mapped IC 63 at high resolution from the UV to the NIR (275 nm to 1.6 $\mu$m), using the Hubble Space Telescope WFC3. Using a Bayesian SED fitting tool, we simultaneously derive a set of stellar ($T_\text{eff}$, $\log(g)$, distance) and extinction ($A_V$, $R_V$) parameters for 520 background stars. We present maps of $A_V$ and $R_V$ with a resolution of 25 arcsec based on these results. The extinction properties vary across the PDR, with values for $A_V$ between 0.5 and 1.4 mag, and a decreasing trend in $R_V$, going from 3.7 at the front of the nebula to values as low as 2.5 further in. This provides evidence for evolution of the dust optical properties. We fit two modified blackbodies to the MIR and FIR SED, obtained by combining the $A_V$ map with data from Spitzer and Herschel. We derive effective temperatures (30 K and 227 K) and the ratio of opacities at 160 $\mu$m to V band $\kappa_{160} / \kappa_V$ ($7.0 \times 10^{-4}$ and $2.9 \times 10^{-9}$) for the two dust populations. Similar fits to individual pixels show spatial variations of $\kappa_{160} / \kappa_{V}$. The analysis of our HST data, combined with these Spitzer and Herschel data, provides the first panchromatic view of dust within a PDR.

We produce a set of 72 NIR through UV extinction curves by combining new HST/STIS optical spectrophotometry with existing IUE spectrophotometry (yielding gapless coverage from 1150 to 10000 Angstroms) and NIR photometry. These curves are used to determine a new, internally consistent, NIR through UV Milky Way Mean Curve and to characterize how the shapes of the extinction curves depend on R(V). We emphasize that while this dependence captures much of the curve variability, there remains considerable variation which is independent of R(V). We use the optical spectrophotometry to verify the presence of structure at intermediate wavelength scales in the curves. The fact that the optical through UV portions of the curves are sampled at relatively high resolution makes them very useful for determining how extinction affects different broad band systems, and we provide several examples. Finally, we compare our results to previous investigations.

Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.

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.

Interstellar dust plays a central role in shaping the detailed structure of the interstellar medium, thus strongly influencing star formation and galaxy evolution. Dust extinction provides one of the main pillars of our understanding of interstellar dust while also often being one of the limiting factors when interpreting observations of distant objects, including resolved and unresolved galaxies. The ultraviolet (UV) and mid-infrared (MIR) wavelength regimes exhibit features of the main components of dust, carbonaceous and silicate materials, and therefore provide the most fruitful avenue for detailed extinction curve studies. Our current picture of extinction curves is strongly biased to nearby regions in the Milky Way. The small number of UV extinction curves measured in the Local Group (mainly Magellanic Clouds) clearly indicates that the range of dust properties is significantly broader than those inferred from the UV extinction characteristics of local regions of the Milky Way. Obtaining statistically significant samples of UV and MIR extinction measurements for all the dusty Local Group galaxies will provide, for the first time, a basis for understanding dust grains over a wide range of environments. Obtaining such observations requires sensitive medium-band UV, blue-optical, and mid-IR imaging and followup R ~ 1000 spectroscopy of thousands of sources. Such a census will revolutionize our understanding of the dependence of dust properties on local environment providing both an empirical description of the effects of dust on observations as well as strong constraints on dust grain and evolution models.

The comparison of chemical abundances in the neutral gas of galaxies to photospheric abundances of old and young stars, ionized gas abundances, and abundances in galactic halos can trace the chemical enrichment of the universe through cosmic times. In particular, our understanding of chemical enrichment through spectroscopic observations of damped Lyman alpha systems (DLAs) relies on corrections for depletion of metals from the gas to the dust phase. These corrections must be determined in the nearby universe, where both gas-phase abundances and photospheric abundances of young stars recently formed out of the interstellar medium can be measured. Multi-object high-resolution (R>50,000) ultraviolet (970-2400 A) and optical (300-600 nm) spectroscopy toward massive stars in local volume galaxies (D < 15 Mpc) covering a wide range of metallicities (a few % solar to solar) and morphological types will provide the abundance and depletion measurements needed to obtain a detailed and comprehensive characterization of the lifecycle of metals in neutral gas and dust in galaxies, thereby observationally addressing important questions about chemical enrichment and galaxy evolution.

We present the severe, systematic uncertainties currently facing our understanding of dust emission, which stymie our ability to truly exploit dust as a tool for studying galaxy evolution. We propose a program of study to tackle these uncertainties, describe the necessary facilities, and discuss the potential science gains that will result. This white paper was submitted to the US National Academies' Astro2020 Decadal Survey on Astronomy and Astrophysics.

Metal Evolution, Transport, and Abundance in the LMC (METAL) is a large Cycle 24 program on the Hubble Space Telescope aimed at measuring dust extinction properties and interstellar depletions in the Large Magellanic Cloud (LMC) at half-solar metallicity. The 101-orbit program is comprised of COS and STIS spectroscopy toward 33 LMC massive stars between 1150 A and 3180 A, and parallel WFC3 imaging in 7 NUV-NIR filters. The fraction of silicon in the gas-phase (depletion) obtained from the spectroscopy decreases with increasing hydrogen column density. Depletion patterns for silicon differ between the Milky Way, LMC, and SMC, with the silicon depletion level offsetting almost exactly the metallicity differences, leading to constant gas-phase abundances in those galaxies for a given hydrogen column density. The silicon depletion correlates linearly with the absolute-to-selective extinction, R$_V$, indicating a link between gas depletion and dust grain size. Extinction maps are derived from the resolved stellar photometry in the parallel imaging, which can be compared to FIR images from Herschel and Spitzer to estimate the emissivity of dust at LMC metallicity. The full METAL sample of depletions, UV extinction curves, and extinction maps will inform the abundance, size, composition, and optical properties of dust grains in the LMC, comprehensively improve our understanding of dust properties, and the accuracy with which dust-based gas masses, star formation rates and histories in nearby and high-redshift galaxies are estimated. This overview paper describes the goals, design, data reduction, and initial results of the METAL survey.

We present interferometric observations of six O-type stars that were made with the Precision Astronomical Visible Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array. The observations include multiple brackets for three targets, $\lambda$~Ori~A, $\zeta$~Oph, and 10~Lac, but there are only preliminary, single observations of the other three stars, $\xi$~Per, $\alpha$~Cam, and $\zeta$~Ori~A. The stellar angular diameters range from 0.55 milliarcsec for $\zeta$~Ori~A down to 0.11 mas for 10~Lac, the smallest star yet resolved with the CHARA Array. The rotational oblateness of the rapidly rotating star $\zeta$ Oph is directly measured for the first time. We assembled ultraviolet to infrared flux measurements for these stars, and then derived angular diameters and reddening estimates using model atmospheres and an effective temperature set by published results from analysis of the line spectrum. The model-based angular diameters are in good agreement with observed angular diameters. We also present estimates for the effective temperatures of these stars derived by setting the interferometric angular size and fitting the spectrophotometry.

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.

At intermediate and high redshifts, measurements of galaxy star-formation rates are usually based on rest-frame ultraviolet (UV) data. A correction for dust attenuation, A_UV, is needed for these measurements. This correction is typically inferred from UV spectral slopes (beta) using an equation known as "Meurer's Relation." In this paper, we study this relation at a redshift of 1.5 using images and photometric measurements in the rest-frame UV (HST) through mid-infrared (Spitzer). It is shown that massive star-forming galaxies (above 10^10 Msun) have dust corrections that are dependent on their inclination to the line-of-sight. Edge-on galaxies have higher A_UV and infrared excess (IRX=L(IR)/L(UV)) than face-on galaxies at a given beta. Interestingly, dust corrections for low-mass star-forming galaxies do not depend on inclination. This is likely because more massive galaxies have more disk-like shapes/kinematics, while low-mass galaxies are more prolate and have more disturbed kinematics. To account for an inclination-dependent dust correction, a modified Meurer's Relation is derived: A_UV=4.43+1.99 beta - 1.73 (b/a-0.67), where b/a is the galaxy axis ratio. This inclination-dependence of A_UV can be explained by a two-component model of the dust distribution inside galaxies. In such a model, the dust attenuation of edge-on galaxies has a higher contribution from a "mixture" component (dust uniformly mixed with stars in the diffuse interstellar medium), and a lower contribution from a "birth cloud" component (near-spherical dust shells surrounding young stars in H II regions) than that of face-on galaxies. The difference is caused by the larger path-lengths through disks at higher inclinations.

We use Hubble Space Telescope (HST) observations of red clump stars taken as part of the Small Magellanic Cloud Investigation of Dust and Gas Evolution (SMIDGE) program to measure the average dust extinction curve in a ~ 200 pc x 100 pc region in the southwest bar of the Small Magellanic Cloud (SMC). The rich information provided by our 8-band ultra-violet through near-infrared photometry allows us to model the color-magnitude diagram of the red clump accounting for the extinction curve shape, a log-normal distribution of $A_{V}$, and the depth of the stellar distribution along the line of sight. We measure an extinction curve with $R_{475} = A_{475}/(A_{475}-A_{814})$ = 2.65 $\pm$ 0.11. This measurement is significantly larger than the equivalent values of published Milky Way $R_{V}$ = 3.1 ($R_{475} = 1.83$) and SMC Bar $R_{V}$ = 2.74 ($R_{475} = 1.86$) extinction curves. Similar extinction curve offsets in the Large Magellanic Cloud (LMC) have been interpreted as the effect of large dust grains. We demonstrate that the line-of-sight depth of the SMC (and LMC) introduces an apparent "gray" contribution to the extinction curve inferred from the morphology of the red clump. We show that no gray dust component is needed to explain extinction curve measurements when a full-width half-max depth of 10 $\pm$ 2 kpc in the stellar distribution of the SMC (5 $\pm$ 1 kpc for the LMC) is considered, which agrees with recent studies of Magellanic Cloud stellar structure. The results of our work demonstrate the power of broad-band HST imaging for simultaneously constraining dust and galactic structure outside the Milky Way.

The most powerful tests of stellar models come from the brightest stars in the sky, for which complementary techniques, such as astrometry, asteroseismology, spectroscopy, and interferometry can be combined. The K2 Mission is providing a unique opportunity to obtain high-precision photometric time series for bright stars along the ecliptic. However, bright targets require a large number of pixels to capture the entirety of the stellar flux, and bandwidth restrictions limit the number and brightness of stars that can be observed. To overcome this, we have developed a new photometric technique, that we call halo photometry, to observe very bright stars using a limited number of pixels. Halo photometry is simple, fast and does not require extensive pixel allocation, and will allow us to use K2 and other photometric missions, such as TESS, to observe very bright stars for asteroseismology and to search for transiting exoplanets. We apply this method to the seven brightest stars in the Pleiades open cluster. Each star exhibits variability; six of the stars show what are most-likely slowly pulsating B-star (SPB) pulsations, with amplitudes ranging from 20 to 2000 ppm. For the star Maia, we demonstrate the utility of combining K2 photometry with spectroscopy and interferometry to show that it is not a 'Maia variable', and to establish that its variability is caused by rotational modulation of a large chemical spot on a 10 d time scale.

We use the Spitzer SAGE survey of the Magellanic Clouds to evaluate the relationship between the 8-mic PAH emission, 24-mic hot dust emission, and HII region radiative transfer. We confirm that in the higher-metallicity Large Magellanic Cloud, PAH destruction is sensitive to optically thin conditions in the nebular Lyman continuum: objects identified as optically thin candidates based on nebular ionization structure show 6 times lower median 8-mic surface brightness (0.18 mJy arcsec^-2) than their optically thick counterparts (1.2 mJy arcsec^-2). The 24-mic surface brightness also shows a factor of 3 offset between the two classes of objects (0.13 vs 0.44 mJy arcsec^-2, respectively), which is driven by the association between the very small dust grains and higher density gas found at higher nebular optical depths. In contrast, PAH and dust formation in the low-metallicity Small Magellanic Cloud is strongly inhibited such that we find no variation in either 8-mic or 24-mic emission between our optically thick and thin samples. This is attributable to extremely low PAH and dust production together with high, corrosive UV photon fluxes in this low-metallicity environment. The dust mass surface densities and gas-to-dust ratios determined from dust maps using Herschel HERITAGE survey data support this interpretation.

Observations and modeling suggest that the dust abundance (gas-to-dust ratio, G/D) depends on (surface) density. The variations of the G/D provide constraints on the timescales for the different processes involved in the lifecycle of metals in galaxies. Recent G/D measurements based on Herschel data suggest a factor 5---10 decrease in the dust abundance between the dense and diffuse interstellar medium (ISM) in the Magellanic Clouds. However, the relative nature of the Herschel measurements precludes definitive conclusions on the magnitude of those variations. We investigate the variations of the dust abundance in the LMC and SMC using all-sky far-infrared surveys, which do not suffer from the limitations of Herschel on their zero-point calibration. We stack the dust spectral energy distribution (SED) at 100, 350, 550, and 850 microns from IRAS and Planck in intervals of gas surface density, model the stacked SEDs to derive the dust surface density, and constrain the relation between G/D and gas surface density in the range 10---100 \Msu pc$^{-2}$ on $\sim$ 80 pc scales. We find that G/D decreases by factors of 3 (from 1500 to 500) in the LMC and 7 (from 1.5$\times 10^4$ to 2000) in the SMC between the diffuse and dense ISM. The surface density dependence of G/D is consistent with elemental depletions and with simple modeling of the accretion of gas-phase metals onto dust grains. This result has important implications for the sub-grid modeling of galaxy evolution, and for the calibration of dust-based gas mass estimates, both locally and at high-redshift.

The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.

The Space Telescope Imaging Spectrograph (STIS) has measured the spectral energy distributions (SEDs) for several stars of types O, B, A, F, and G. These absolute fluxes from the CALSPEC database are fit with a new spectral grid computed from the ATLAS-APOGEE ATLAS9 model atmosphere database using a chi-square minimization technique in four parameters. The quality of the fits are compared for complete LTE grids by Castelli & Kurucz (CK04) and our new comprehensive LTE grid (BOSZ). For the cooler stars, the fits with the MARCS LTE grid are also evaluated, while the hottest stars are also fit with the NLTE Lanz & Hubeny OB star grids. Unfortunately, these NLTE models do not transition smoothly in the infrared to agree with our new BOSZ LTE grid at the NLTE lower limit of Teff =15,000K. The new BOSZ grid is available via the Space Telescope Institute MAST archive and has a much finer sampled IR wavelength scale than CK04, which will facilitate the modeling of stars observed by the James Webb Space Telescope (JWST). Our result for the angular diameter of Sirius agrees with the ground- based interferometric value.

We present an update to the ultraviolet-to-radio database of global broadband photometry for the 79 nearby galaxies that comprise the union of the KINGFISH (Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel) and SINGS (Spitzer Infrared Nearby Galaxies Survey) samples. The 34-band dataset presented here includes contributions from observational work carried out with a variety of facilities including GALEX, SDSS, PS, NOAO, 2MASS, WISE, Spitzer, Herschel, Planck, JCMT, and the VLA. Improvements of note include recalibrations of previously-published SINGS BVRcIc and KINGFISH far-infrared/submillimeter photometry. Similar to previous results in the literature, an excess of submillimeter emission above model predictions is seen primarily for low-metallicity dwarf/irregular galaxies. This 34-band photometric dataset for the combined KINGFISH$+$SINGS sample serves as an important multi-wavelength reference for the variety of galaxies observed at low redshift. A thorough analysis of the observed spectral energy distributions is carried out in a companion paper.

Dust modeling is crucial to infer dust properties and budget for galaxy studies. However, there are systematic disparities between dust grain models that result in corresponding systematic differences in the inferred dust properties of galaxies. Quantifying these systematics requires a consistent fitting analysis. We compare the output dust parameters and assess the differences between two dust grain models, Compiegne et al (2011), and THEMIS (Jones et al 2013, Kohler et al 2015). In this study, we use a single fitting method applied to all the models to extract a coherent and unique statistical analysis. We fit the models to the dust emission seen by Spitzer and Herschel in the Small and Large Magellanic Clouds (SMC and LMC). The observations cover the infrared (IR) spectrum from a few microns to the sub-millimeter range. For each fitted pixel, we calculate the full n-D likelihood, based on the method described in Gordon et al (2014). The free parameters are both environmental ($U$, the interstellar radiation field strength; $\alpha_\mathrm{ISRF}$, power-law coefficient for a multi-U environment; $\Omega^*$, the starlight strength) and intrinsic to the model ($Y_\mathrm{i}$: abundances of the grain species $i$; $\alpha_\mathrm{sCM20}$, coefficient in the small carbon grain size distribution). Fractional residuals of 5 different sets of parameters show that fitting THEMIS brings a more accurate reproduction of the observations than the Compiegne model. However, independent variations of the dust species show strong model-dependencies. We find that the abundance of silicates can only be constrained to an upper-limit and the silicate/carbon ratio is different than that seen in our Galaxy. In the LMC, our fits result in dust masses slightly lower than those found in literature , by a factor lower than 2. In the SMC, we find dust masses in agreement with previous studies.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

The Mid-Infrared Instrument (MIRI) Medium Resolution Spectrometer (MRS) is the only mid-IR Integral Field Spectrometer on board James Webb Space Telescope. The complexity of the MRS requires a very specialized pipeline, with some specific steps not present in other pipelines of JWST instruments, such as fringe corrections and wavelength offsets, with different algorithms for point source or extended source data. The MRS pipeline has also two different variants: the baseline pipeline, optimized for most foreseen science cases, and the optimal pipeline, where extra steps will be needed for specific science cases. This paper provides a comprehensive description of the MRS Calibration Pipeline from uncalibrated slope images to final scientific products, with brief descriptions of its algorithms, input and output data, and the accessory data and calibration data products necessary to run the pipeline.

We report on interferometric observations with the CHARA Array of two classical Wolf-Rayet stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H-band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems ($f_{\rm WR, 137} = 0.59\pm0.04$; $f_{\rm WR, 138} = 0.67\pm0.01$), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modeling of each system. Our spectral modeling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edge-on. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane.

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys.

The Magellanic Clouds provide the only laboratory to study the effect of metallicity and galaxy mass on molecular gas and star formation at high (~20 pc) resolution. We use the dust emission from HERITAGE Herschel data to map the molecular gas in the Magellanic Clouds, avoiding the known biases of CO emission as a tracer of H$_{2}$. Using our dust-based molecular gas estimates, we find molecular gas depletion times of ~0.4 Gyr in the LMC and ~0.6 SMC at 1 kpc scales. These depletion times fall within the range found for normal disk galaxies, but are shorter than the average value, which could be due to recent bursts in star formation. We find no evidence for a strong intrinsic dependence of the molecular gas depletion time on metallicity. We study the relationship between gas and star formation rate across a range in size scales from 20 pc to ~1 kpc, including how the scatter in molecular gas depletion time changes with size scale, and discuss the physical mechanisms driving the relationships. We compare the metallicity-dependent star formation models of Ostriker, McKee, and Leroy (2010) and Krumholz (2013) to our observations and find that they both predict the trend in the data, suggesting that the inclusion of a diffuse neutral medium is important at lower metallicity.

Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.

We present and describe the astro-photometric catalog of more than 800,000 sources found in the Hubble Tarantula Treasury Project (HTTP). HTTP is a Hubble Space Telescope (HST) Treasury program designed to image the entire 30 Doradus region down to the sub-solar (~0.5 solar masses) mass regime using the Wide Field Camera 3 (WFC3) and the Advanced Camera for Surveys (ACS). We observed 30 Doradus in the near ultraviolet (F275W, F336W), optical (F555W, F658N, F775W), and near infrared (F110W, F160W) wavelengths. The stellar photometry was measured using point-spread function (PSF) fitting across all the bands simultaneously. The relative astrometric accuracy of the catalog is 0.4 mas. The astro-photometric catalog, results from artificial star experiments and the mosaics for all the filters are available for download. Color-magnitude diagrams are presented showing the spatial distributions and ages of stars within 30 Dor as well as in the surrounding fields. HTTP provides the first rich and statistically significant sample of intermediate and low mass pre-main sequence candidates and allows us to trace how star formation has been developing through the region. The depth and high spatial resolution of our analysis highlight the dual role of stellar feedback in quenching and triggering star formation on the giant HII region scale. Our results are consistent with stellar sub-clustering in a partially filled gaseous nebula that is offset towards our side of the Large Magellanic Cloud.

New low-resolution UV spectra of a sample of reddened OB stars in M31 were obtained with HST/STIS to study the wavelength dependence of interstellar extinction and the nature of the underlying dust grain populations. Extinction curves were constructed for four reddened sightlines in M31 paired with closely matching stellar atmosphere models. The new curves have a much higher S/N than previous studies. Direct measurements of N(H I) were made using the Ly$\alpha$ absorption lines enabling gas-to-dust ratios to be calculated. The sightlines have a range in galactocentric distance of 5 to 14 kpc and represent dust from regions of different metallicities and gas-to-dust ratios. The metallicities sampled range from Solar to 1.5 Solar. The measured curves show similarity to those seen in the Milky Way and the Large Magellanic Cloud. The Maximum Entropy Method was used to investigate the dust composition and size distribution for the sightlines observed in this program finding that the extinction curves can be produced with the available carbon and silicon abundances if the metallicity is super-Solar.

MWC 314 is a bright candidate luminous blue variable that resides in a fairly close binary system, with an orbital period of 60.753$\pm$0.003 d. We observed MWC 314 with a combination of optical spectroscopy, broad-band ground- and space-based photometry, as well as with long baseline, near-infrared interferometry. We have revised the single-lined spectroscopic orbit and explored the photometric variability. The orbital light curve displays two minima each orbit that can be partially explained in terms of the tidal distortion of the primary that occurs around the time of periastron. The emission lines in the system are often double-peaked and stationary in their kinematics, indicative of a circumbinary disc. We find that the stellar wind or circumbinary disc is partially resolved in the K\prime-band with the longest baselines of the CHARA Array. From this analysis, we provide a simple, qualitative model in an attempt to explain the observations. From the assumption of Roche Lobe overflow and tidal synchronisation at periastron, we estimate the component masses to be M1 $\approx 5$ M$_\odot$ and M2$\approx 15$ M$_\odot$, which indicates a mass of the LBV that is extremely low. In addition to the orbital modulation, we discovered two pulsational modes with the MOST satellite. These modes are easily supported by a low-mass hydrogen-poor star, but cannot be easily supported by a star with the parameters of an LBV. The combination of these results provides evidence that the primary star was likely never a normal LBV, but rather is the product of binary interactions. As such, this system presents opportunities for studying mass-transfer and binary evolution with many observational techniques.

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

We present Spitzer/Infrared Spectrograph (IRS) 5-21 micron spectroscopic maps towards 12 regions in the Andromeda galaxy (M31). These regions include the nucleus, bulge, an active region in the star-forming ring, and 9 other regions chosen to cover a range of mid-to-far-infrared colours. In line with previous results, PAH feature ratios (6.2 micron and 7.7 micron features compared to the 11.2 micron feature) measured from our extracted M31 spectra, except the nucleus, strongly correlate. The equivalent widths of the main PAH features, as a function of metallicity and radiation hardness, are consistent with those observed for other nearby spiral and starburst galaxies. Reprocessed data from the ISOCAM instrument on the Infrared Space Observatory agree with the IRS data; early reports of suppressed 6-8 micron features and enhanced 11.3 micron feature intensity and FWHM apparently resulted from background-subtraction problems. The nucleus does not show any PAH emission but does show strong silicate emission at 9.7 micron. Furthermore, different spectral features (11.3 micron PAH emission, silicate emission and [NeIII] 15.5 micron line emission) have distinct spatial distributions in the nuclear region: the silicate emission is strongest towards the stellar nucleus, while the PAH emission peaks 15 arcsec north of the nucleus. The PAH feature ratios at this position are atypical with strong emission at 11.2 microns and 15-20 microns but weak emission at 6--8 microns. The nucleus itself is dominated by stellar light giving rise to a strong blue continuum and silicate emission.

In this article, we describe the MIRI Imager module (MIRIM), which provides broad-band imaging in the 5 - 27 microns wavelength range for the James Webb Space Telescope. The imager has a 0"11 pixel scale and a total unobstructed view of 74"x113". The remainder of its nominal 113"x113" field is occupied by the coronagraphs and the low resolution spectrometer. We present the instrument optical and mechanical design. We show that the test data, as measured during the test campaigns undertaken at CEA-Saclay, at the Rutherford Appleton Laboratory, and at the NASA Goddard Space Flight Center, indicate that the instrument complies with its design requirements and goals. We also discuss the operational requirements (multiple dithers and exposures) needed for optimal scientific utilization of the MIRIM.

We describe the operations concept and data reduction plan for the Mid- Infrared Instrument (MIRI) for the James Webb Space Telescope (JWST). The overall JWST operations concept is to use Observation Templates (OTs) to provide a straightforward and intuitive way for users to specify observations. MIRI has four OTs that correspond to the four observing modes: 1.) Imaging, 2.) Coronagraphy, 3.) Low Resolution Spectroscopy, and 4.) Medium Resolution Spectroscopy. We outline the user choices and expansion of these choices into detailed instrument operations. The data reduction plans for MIRI are split into three stages, where the specificity of the reduction steps to the observation type increases with stage. The reduction starts with integration ramps: stage 1 yields uncalibrated slope images; stage 2 calibrates the slope images; and then stage 3 combines multiple calibrated slope images into high level data products (e.g. mosaics, spectral cubes, and extracted source information). Finally, we give examples of the data and data products that will be derived from each of the four different OTs.

MIRI (the Mid-Infrared Instrument for the James Webb Space Telescope (JWST)) operates from 5 to 28.5 microns and combines over this range: 1.) unprecedented sensitivity levels; 2.) sub-arcsec angular resolution; 3.) freedom from atmospheric interference; 4.) the inherent stability of observing in space; and 5.) a suite of versatile capabilities including imaging, low and medium resolution spectroscopy (with an integral field unit), and coronagraphy. We illustrate the potential uses of this unique combination of capabilities with various science examples: 1.) imaging exoplanets; 2.) transit and eclipse spectroscopy of exoplanets; 3.) probing the first stages of star and planet formation, including identifying bioactive molecules; 4.) determining star formation rates and mass growth as galaxies are assembled; and 5.) characterizing the youngest massive galaxies. This paper is the introduction to a series of ten covering all aspects of the instrument.

We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.

We present a new method for deriving the stellar birth function (SBF) of resolved stellar populations. The SBF (stars born per unit mass, time, and metallicity) is the combination of the initial mass function (IMF), the star-formation history (SFH), and the metallicity distribution function (MDF). The framework of our analysis is that of Poisson Point Processes (PPPs), a class of statistical models suitable when dealing with points (stars) in a multidimensional space (the measurement space of multiple photometric bands). The theory of PPPs easily accommodates the modeling of measurement errors as well as that of incompleteness. Compared to most of the tools used to study resolved stellar populations, our method avoids binning stars in the color-magnitude diagram and uses the entirety of the information (i.e., the whole likelihood function) for each data point; the proper combination of the individual likelihoods allows the computation of the posterior probability for the global population parameters. This includes unknowns such as the IMF slope and combination of SFH and MDF, which are rarely solved for simultaneously in the literature, however entangled and correlated they might be. Our method also allows proper inclusion of nuisance parameters, such as distance and extinction distributions. The aim of this paper, is to assess the validity of this new approach under a range of assumptions, using only simulated data. Forthcoming work will show applications to real data. Although it has a broad scope of possible applications, we have developed this method to study multi-band HST observations of the Milky Way Bulge. Therefore we will focus on simulations with characteristics similar to those of the Galactic Bulge.

We present the analysis of supernova remnants (SNRs) in the Large Magellanic Cloud (LMC) and their influence on the environment at far-infrared (FIR) and submillimeter wavelengths. We use new observations obtained with the {\it Herschel} Space Observatory and archival data obtained with the {\it Spitzer} Space Telescope, to make the first FIR atlas of these objects. The SNRs are not clearly discernible at FIR wavelengths, however their influence becomes apparent in maps of dust mass and dust temperature, which we constructed by fitting a modified black-body to the observed spectral energy distribution in each sightline. Most of the dust that is seen is pre-existing interstellar dust in which SNRs leave imprints. The temperature maps clearly reveal SNRs heating surrounding dust, while the mass maps indicate the removal of 3.7$^{+7.5}_{-2.5}$ M$_{\odot}$ of dust per SNR. This agrees with the calculations by others that significant amounts of dust are sputtered by SNRs. Under the assumption that dust is sputtered and not merely pushed away, we estimate a dust destruction rate in the LMC of $0.037^{+0.075}_{-0.025}$ M$_\odot$ yr$^{-1}$ due to SNRs, yielding an average lifetime for interstellar dust of $2^{+4.0}_{-1.3}\times10^7$ yr. We conclude that sputtering of dust by SNRs may be an important ingredient in models of galactic evolution, that supernovae may destroy more dust than they produce, and that they therefore may not be net producers of long lived dust in galaxies.

The rate at which interstellar gas is converted into stars, and its dependence on environment, is one of the pillars on which our understanding of the visible Universe is build. We present a comparison of the surface density of young stars (Sigma_*) and dust surface density (Sigma_d) across NGC346 (N66) in 115 independent pixels of 6x6 pc^2. We find a correlation between Sigma_* and Sigma_d with a considerable scatter. A power law fit to the data yields a steep relation with an exponent of 2.6+-0.2. We convert Sigma_d to gas surface density (Sigma_g) and Sigma_* to star formation rate (SFR) surface densities (Sigma_SFR), using simple assumptions for the gas-to-dust mass ratio and the duration of star formation. The derived total SFR (4+-1 10^-3 M_sun/yr) is consistent with SFR estimated from the Ha emission integrated over the Ha nebula. On small scales the Sigma_SFR derived using Ha systematically underestimates the count-based Sigma_SFR, by up to a factor of 10. This is due to ionizing photons escaping the area, where the stars are counted. We find that individual 36 pc^2 pixels fall systematically above integrated disc-galaxies in the Schmidt-Kennicutt diagram by on average a factor of ~7. The NGC346 average SFR over a larger area (90 pc radius) lies closer to the relation but remains high by a factor of ~3. The fraction of the total mass (gas plus young stars) locked in young stars is systematically high (~10 per cent) within the central 15 pc and systematically lower outside (2 per cent), which we interpret as variations in star formation efficiency. The inner 15 pc is dominated by young stars belonging to a centrally condensed cluster, while the outer parts are dominated by a dispersed population. Therefore, the observed trend could reflect a change of star formation efficiency between clustered and non-clustered star-formation.

The DUSTiNGS survey (DUST in Nearby Galaxies with Spitzer) is a 3.6 and 4.5 micron imaging survey of 50 nearby dwarf galaxies designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. Using two epochs, spaced approximately six months apart, we identify a total of 526 dusty variable AGB stars (sometimes called "extreme" or x-AGB stars; [3.6]-[4.5]>0.1 mag). Of these, 111 are in galaxies with [Fe/H] < -1.5 and 12 are in galaxies with [Fe/H] < -2.0, making them the most metal-poor dust-producing AGB stars known. We compare these identifications to those in the literature and find that most are newly discovered large-amplitude variables, with the exception of approximately 30 stars in NGC 185 and NGC 147, one star in IC 1613, and one star in Phoenix. The chemical abundances of the x-AGB variables are unknown, but the low metallicities suggest that they are more likely to be carbon-rich than oxygen-rich and comparisons with existing optical and near-IR photometry confirms that 70 of the x-AGB variables are confirmed or likely carbon stars. We see an increase in the pulsation amplitude with increased dust production, supporting previous studies suggesting that dust production and pulsation are linked. We find no strong evidence linking dust production with metallicity, indicating that dust can form in very metal-poor environments.

The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.

Nearby resolved dwarf galaxies provide excellent opportunities for studying the dust-producing late stages of stellar evolution over a wide range of metallicity (-2.7 < [Fe/H] < -1.0). Here, we describe DUSTiNGS (DUST in Nearby Galaxies with Spitzer): a 3.6 and 4.5 micron post-cryogen Spitzer Space Telescope imaging survey of 50 dwarf galaxies within 1.5 Mpc that is designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. The survey includes 37 dwarf spheroidal, 8 dwarf irregular, and 5 transition-type galaxies. This near-complete sample allows for the building of statistics on these rare phases of stellar evolution over the full metallicity range. The photometry is >75% complete at the tip of the Red Giant Branch for all targeted galaxies, with the exception of the crowded inner regions of IC 10, NGC 185, and NGC 147. This photometric depth ensures that the majority of the dust-producing stars, including the thermally-pulsing AGB stars, are detected in each galaxy. The images map each galaxy to at least twice the half-light radius to ensure that the entire evolved star population is included and to facilitate the statistical subtraction of background and foreground contamination, which is severe at these wavelengths. In this overview, we describe the survey, the data products, and preliminary results. We show evidence for the presence of dust-producing AGB stars in 8 of the targeted galaxies, with metallicities as low as [Fe/H] = -1.9, suggesting that dust production occurs even at low metallicity.

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm \mu m}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$\sigma$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm \mu m}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions.

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.

We have measured stellar photometry with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) in near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) bands for 117 million resolved stars in M31. As part of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, we measured photometry with simultaneous point spread function fitting across all bands and at all source positions after precise astrometric image alignment (<5-10 milliarcsecond accuracy). In the outer disk, the photometry reaches a completeness-limited depth of F475W~28, while in the crowded, high surface brightness bulge, the photometry reaches F475W~25. We find that simultaneous photometry and optimized measurement parameters significantly increase the detection limit of the lowest resolution filters (WFC3/IR) providing color-magnitude diagrams that are up to 2.5 magnitudes deeper when compared with color-magnitude diagrams from WFC3/IR photometry alone. We present extensive analysis of the data quality including comparisons of luminosity functions and repeat measurements, and we use artificial star tests to quantify photometric completeness, uncertainties and biases. We find that largest sources of systematic error in the photometry are due to spatial variations in the point spread function models and charge transfer efficiency corrections. This stellar catalog is the largest ever produced for equidistant sources, and is publicly available for download by the community.

We model multi-wavelength broadband UBVIJHKs and Spitzer IRAC and MIPS photometry and IRS spectra from the SAGE and SAGE-Spec observing programs of two oxygen-rich asymptotic giant branch (O-rich AGB) stars in the Large Magellanic Cloud (LMC) using radiative transfer models of dust shells around stars. We chose a star from each of the bright and faint O-rich AGB populations found by earlier studies of the SAGE sample in order to derive a baseline set of dust properties to be used in the construction of an extensive grid of radiative transfer models of the O-rich AGB stars found in the SAGE surveys. From the bright O-rich AGB population we chose HV 5715, and from the faint O-rich AGB population we chose SSTISAGE1C J052206.92-715017.6 (SSTSAGE052206). We found the complex indices of refraction of oxygen-deficient silicates from Ossenkopf et al (1992) and a "KMH"-like grain size distribution with gamma of -3.5, a_min of 0.01 microns, and a_0 of 0.1 microns to be reasonable dust properties for these models. There is a slight indication that the dust around the faint O-rich AGB may be more silica-rich than that around the bright O-rich AGB. Simple models of gas emission suggest a relatively extended gas envelope for the faint O-rich AGB star modeled, consistent with the relatively large dust shell inner radius for the same model. The models compute the dust mass-loss rates for the two stars to be 2.0x10^-9 M_Sun/yr and 2.3x10^-9 M_Sun/yr, respectively. When a dust-to-gas mass ratio of 0.002 is assumed for SSTSAGE052206 and HV 5715, the dust mass-loss rates imply total mass-loss rates of 1.0x10^-6 M_Sun/yr and 1.2x10^-6 M_Sun/yr, respectively. These properties of the dust shells and stars, as inferred from our models of the two stars, are found to be consistent with properties observed or assumed by detailed studies of other O-rich AGB stars in the LMC and elsewhere.

We search Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and Wide Field Camera 3 (WFC3) broadband imaging data from the Panchromatic Hubble Andromeda Treasury (PHAT) survey to identify detections of cataloged planetary nebulae (PNe). Of the 711 PNe currently in the literature within the PHAT footprint, we find 467 detected in the broadband. For these 467 we are able to refine their astrometric accuracy from ~0."3 to 0."05. Using the resolution of HST, we are able to show that 152 objects currently in the catalogs are definitively not PNe, and we show that 32 objects thought to be extended in ground-based images are actually point-like and therefore good PN candidates. We also find one PN candidate that is marginally resolved. If this is a PN, it is up to 0.7 pc in diameter. With our new photometric data, we develop a method of measuring the level of excitation in individual PNe by comparing broadband and narrowband imaging and describe the effects of excitation on a PN's photometric signature. Using the photometric properties of the known PNe in the PHAT catalogs, we search for more PN, but do not find any new candidates, suggesting that ground-based emission-line surveys are complete in the PHAT footprint to F475W $\simeq$ 24.

In this paper we present a catalogue of Giant Molecular Clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Hershel Exploitation of Local Galaxy Andromeda (HELGA) dataset. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalogue and characterise the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al (2012). 236 GMCs in the mass range 10^4-10^7 M_sol are identified, their cumulative mass distribution is found to be proportional to M^-1.45 in agreement with earlier studies. The GMCs appear to follow the same cloud mass to L_CO correlation observed in the Milky Way. However, comparison between this catalogue and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit suggesting that we are observing associations of GMCs. Following Gordon et al. (2006), we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.5 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.9 deg to the GMCs not associated with either ring. Lastly, we comment upon the effects of deprojection on our results and investigate the effect different models for M31's inclination will have upon the projection of an unperturbed spiral arm system.

We investigate the clustering of early-type stars younger than 300 Myr on galactic scales in M31. Based on the stellar photometric catalogs of the Panchromatic Hubble Andromeda Treasury program that also provides stellar parameters derived from the individual energy distributions, our analysis is focused on the young stars in three star-forming regions, located at galactocentric distances of about 5, 10, and 15 kpc, corresponding to the inner spiral arms, the ring structure, and the outer arm, respectively. We apply the two-point correlation function to our selected sample to investigate the clustering behavior of these stars across different time- and length-scales. We find that young stellar structure survives across the whole extent of M31 longer than 300 Myr. Stellar distribution in all regions appears to be self-similar, with younger stars being systematically more strongly clustered than the older, which are more dispersed. The observed clustering is interpreted as being induced by turbulence, the driving source for which is probably gravitational instabilities driven by the spiral arms, which are stronger closer to the galactic centre.

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

The measurement of precise absolute fluxes for stellar sources has been pursued with increased vigor since the discovery of the dark energy and the realization that its detailed understanding requires accurate spectral energy distributions (SEDs) of redshifted Ia supernovae in the rest frame. The flux distributions of spectrophotometric standard stars were initially derived from the comparison of stars to laboratory sources of known flux but are now mostly based on calculated model atmospheres. For example, pure hydrogen white dwarf (WD) models provide the basis for the HST CALSPEC archive of flux standards. The basic equations for quantitative spectrophotometry and photometry are explained in detail. Several historical lab based flux calibrations are reviewed; and the SEDs of stars in the major on-line astronomical databases are compared to the CALSPEC reference standard spectrophotometry. There is good evidence that relative fluxes from the visible to the near-IR wavelength of ~2.5 micron are currently accurate to 1% for the primary reference standards; and new comparisons with lab flux standards show promise for improving that precision.

We investigate the far-infrared (IR) dust emission for 20 local star forming galaxies from the Key Insights on Nearby Galaxies: A Far-IR Survey with Herschel (KINGFISH) sample. We model the far-IR/submillimeter spectral energy distribution (SED) using images from Spitzer Space Telescope and Herschel Space Observatory. We calculate the cold dust temperature (T(cold)) and emissivity (beta) on a pixel by pixel basis (where each pixel ranges from 0.1-3 kpc^2) using a two temperature modified blackbody fitting routine. Our fitting method allows us to investigate the resolved nature of temperature and emissivity variations by modeling from the galaxy centers to the outskirts (physical scales of ~15-50 kpc, depending on the size of the galaxy). We fit each SED in two ways: (1) fit T(cold) and beta simultaneously, (2) hold beta constant and fit T(cold). We compare T(cold) and beta with star formation rates (calculated from L(Halpha) and L(24)), the luminosity of the old stellar population (traced through L(3.6), and the dust mass surface density (traced by 500 micron luminosity, L(500)). We find a significant trend between SFR/L(500) and T(cold), implying that the flux of hard UV photons relative to the amount of dust is significantly contributing to the heating of the cold, or diffuse, dust component. We also see a trend between L(3.6)/L(500) and beta, indicating that the old stellar population contributes to the heating at far-IR/submillimeter wavelengths. Finally, we find that when beta is held constant, T(cold) exhibits a strongly decreasing radial trend, illustrating that the shape of the far-IR SED is changing radially through a galaxy, thus confirming on a sample almost double in size the trends observed in Galametz et al. (2012).

We measure the extinction curve in the central 200 pc of M31 at mid-ultraviolet to near-infrared wavelengths (from 1928A to 1.5{\mu}m), using Swift/UVOT and HST WFC3/ACS observations in thirteen bands. Taking advantage of the high angular resolution of the HST WFC3 and ACS detectors, we develop a method to simultaneously determine the relative extinction and the fraction of obscured starlight for five dusty complexes located in the circumnuclear region. The extinction curves of these clumps (RV =2.4-2.5) are steeper than the average Galactic one (RV =3.1), but are similar to optical and near-infrared curves recently measured toward the Galactic Bulge (RV~2.5). This similarity suggests that steep extinction curves may be common in the inner bulge of galaxies. In the ultraviolet, the extinction curves of these clumps are also unusual. We find that one dusty clump (size <2 pc) exhibits a strong UV bump (extinction at 2175A), more than three standard deviation higher than that predicted by common models. Although the high stellar metallicity of the M31 bulge indicates that there are sufficient carbon and silicon to produce large dust grains, the grains may have been destroyed by supernova explosions or past activity of the central super-massive black hole, resulting in the observed steepened extinction curve.

The Magellanic Bridge is the nearest low-metallicity, tidally stripped environment, offering a unique high-resolution view of physical conditions in merging and forming galaxies. In this paper we present analysis of candidate massive young stellar objects (YSOs), i.e., {\it in situ, current} massive star formation (MSF) in the Bridge using {\it Spitzer} mid-IR and complementary optical and near-IR photometry. While we definitely find YSOs in the Bridge, the most massive are $\sim10 M_\odot$, $\ll45 M_\odot$ found in the Large Magellanic Cloud (LMC). The intensity of MSF in the Bridge also appears decreasing, as the most massive YSOs are less massive than those formed in the past. To investigate environmental effects on MSF, we have compared properties of massive YSOs in the Bridge to those in the LMC. First, YSOs in the Bridge are apparently less embedded than in the LMC: 81% of Bridge YSOs show optical counterparts, compared to only 56% of LMC sources with the same range of mass, circumstellar dust mass, and line-of-sight extinction. Circumstellar envelopes are evidently more porous or clumpy in the Bridge's low-metallicity environment. Second, we have used whole samples of YSOs in the LMC and the Bridge to estimate the probability of finding YSOs at a given \hi\ column density, N(HI). We found that the LMC has $\sim3\times$ higher probability than the Bridge for N(HI) $>10\times10^{20}$ cm$^{-2}$, but the trend reverses at lower N(HI). Investigating whether this lower efficiency relative to HI is due to less efficient molecular cloud formation, or less efficient cloud collapse, or both, will require sensitive molecular gas observations.

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.

'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.

Motived by recent ISM studies, we present high quality MMT and Gemini spectroscopic observations of H II regions in the nearby spiral galaxies NGC 628 and NGC 2403 in order to measure their chemical abundance gradients. Using long-slit and multi-object mask optical spectroscopy, we obtained measurements of the temperature sensitive auroral lines [O III] {\lambda}4363 and/or [N II] {\lambda}5755 at a strength of 4{\sigma} or greater in 11 H II regions in NGC 628 and 7 regions in NGC 2403. These observations allow us, for the first time, to derive an oxygen abundance gradient in NGC 628 based solely on "direct" oxygen abundances of H II regions: 12 + log(O/H) = (8.43+/-0.03) + (-0.017+/-0.002) x Rg (dex/kpc), with a dispersion in log(O/H) of {\sigma} = 0.10 dex, from 14 regions with a radial coverage of ~2-19 kpc. This is a significantly shallower slope than found by previous "strong-line" abundance studies. In NGC 2403, we derive an oxygen abundance gradient of 12 + log(O/H) = (8.48+/-0.04) + (-0.032+/-0.007) x Rg (dex/kpc), with a dispersion in log(O/H) of {\sigma} = 0.07 dex, from 7 H II with a radial coverage of ~1-10 kpc. Additionally, we measure the N, S, Ne, and Ar abundances. We find the N/O ratio decreases with increasing radius for the inner disk, but reaches a plateau past R25 in NGC 628. NGC 2403 also has a negative N/O gradient with radius, but we do not sample the outer disk of the galaxy past R25 and so do not see evidence for a plateau. This bi-modal pattern measured for NGC 628 indicates dominant contributions from secondary nitrogen inside of the R25 transition and dominantly primary nitrogen farther out. As expected for {\alpha}-process elements, S/O, Ne/O, and Ar/O are consistent with constant values over a range in oxygen abundance.

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.

Cosmic dust is present in many astrophysical objects, and recent observations across the electromagnetic spectrum show that the dust distribution is often strongly three-dimensional (3D). Dust grains are effective in absorbing and scattering ultraviolet (UV)/optical radiation, and they re-emit the absorbed energy at infrared wavelengths. Understanding the intrinsic properties of these objects, including the dust itself, therefore requires 3D dust radiative transfer (RT) calculations. Unfortunately, the 3D dust RT problem is nonlocal and nonlinear, which makes it one of the hardest challenges in computational astrophysics. Nevertheless, significant progress has been made in the past decade, with an increasing number of codes capable of dealing with the complete 3D dust RT problem. We discuss the complexity of this problem, the two most successful solution techniques [ray-tracing (RayT) and Monte Carlo (MC)], and the state of the art in modeling observational data using 3D dust RT codes. We end with an outlook on the bright future of this field.

Context. Edge-on spiral galaxies with prominent dust lanes provide us with an excellent opportunity to study the distribution and properties of the dust within them. The HEROES project was set up to observe a sample of seven large edge-on galaxies across various wavelengths for this investigation. Aims. Within this first paper, we present the Herschel observations and perform a qualitative and quantitative analysis on them, and we derive some global properties of the far infrared and submillimetre emission. Methods. We determine horizontal and vertical profiles from the Herschel observations of the galaxies in the sample and describe the morphology. Modified black-body fits to the global fluxes, measured using aperture photometry, result in dust temperatures and dust masses. The latter values are compared to those that are derived from radiative transfer models taken from the literature. Results. On the whole, our Herschel flux measurements agree well with archival values. We find that the exponential horizontal dust distribution model often used in the literature generally provides a good description of the observed horizontal profiles. Three out of the seven galaxies show signatures of extended vertical emission at 100 and 160 {\mu}m at the 5{\sigma} level, but in two of these it is probably due to deviations from an exactly edge-on orientation. Only for NGC 4013, a galaxy in which vertically extended dust has already been detected in optical images, we can detect vertically extended dust, and the derived scaleheight agrees with the value estimated through radiative transfer modelling. Our analysis hints at a correlation between the dust scaleheight and its degree of clumpiness, which we infer from the difference between the dust masses as calculated from modelling of optical data and from fitting the spectral energy distribution of Herschel datapoints.

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.

The Hubble Tarantula Treasury Project (HTTP) is an ongoing panchromatic imaging survey of stellar populations in the Tarantula Nebula in the Large Magellanic Cloud that reaches into the sub-solar mass regime (< 0.5 Mo). HTTP utilizes the capability of HST to operate the Advanced Camera for Surveys (ACS) and the Wide Field Camera 3 (WFC3) in parallel to study this remarkable region in the near-ultraviolet, optical, and near-infrared spectral regions, including narrow band H$\alpha$ images. The combination of all these bands provides a unique multi-band view. The resulting maps of the stellar content of the Tarantula Nebula within its main body provide the basis for investigations of star formation in an environment resembling the extreme conditions found in starburst galaxies and in the early Universe. Access to detailed properties of individual stars allows us to begin to reconstruct the evolution of the stellar skeleton of the Tarantula Nebula over space and time with parcsec-scale resolution. In this first paper we describe the observing strategy, the photometric techniques, and the upcoming data products from this survey and present preliminary results obtained from the analysis of the initial set of near-infrared observations.

A crucial aspect of 3D Monte Carlo radiative transfer is the choice of the spatial grid used to partition the dusty medium. We critically investigate the use of octree grids in Monte Carlo dust radiative transfer, with two different octree construction algorithms (regular and barycentric subdivision) and three different octree traversal algorithms (top-down, neighbour list, and the bookkeeping method). In general, regular octree grids need higher levels of subdivision compared to the barycentric grids for a fixed maximum cell mass threshold criterion. The total number of grid cells, however, depends on the geometry of the model. Surprisingly, regular octree grid simulations turn out to be 10 to 20% more efficient in run time than the barycentric grid simulations, even for those cases where the latter contain fewer grid cells than the former. Furthermore, we find that storing neighbour lists for each cell in an octree, ordered according to decreasing overlap area, is worth the additional memory and implementation overhead: using neighbour lists can cut down the grid traversal by 20% compared to the traditional top-down method. In conclusion, the combination of a regular node subdivision and the neighbour list method results in the most efficient octree structure for Monte Carlo radiative transfer simulations.

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 a detailed study of how the Star Formation Rate (SFR) relates to the interstellar medium (ISM) of M31 at ~140pc scales. The SFR is calculated using the far-ultraviolet and 24um emission, corrected for the old stellar population in M31. We find a global value for the SFR of 0.25+/-0.05Msun/yr and compare this with the SFR found using the total far-infrared (FIR) luminosity. There is general agreement in regions where young stars dominate the dust heating. Atomic hydrogen (HI) and molecular gas (traced by carbon monoxide, CO) or the dust mass is used to trace the total gas in the ISM. We show that the global surface densities of SFR and gas mass place M31 amongst a set of low-SFR galaxies in the plot of Kennicutt (1998b). The relationship between SFR and gas surface density is tested in six radial annuli across M31, assuming a power law relationship with index, N. The star formation law using total gas traced by HI and CO gives a global index of N=2.03+/-0.04, with a significant variation with radius; the highest values are observed in the 10kpc ring. We suggest that this slope is due to HI turning molecular at ~10Msun/pc2. When looking at H2 regions, we measure a higher mean SFR suggesting a better spatial correlation between H2 and SF. We find N~0.6 with consistent results throughout the disk - this is at the low end of values found in previous work and argues against a superlinear SF law on small scales.

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.

We present a probabilistic approach for inferring the parameters of the present day power-law stellar mass function (MF) of a resolved young star cluster. This technique (a) fully exploits the information content of a given dataset; (b) accounts for observational uncertainties in a straightforward way; (c) assigns meaningful uncertainties to the inferred parameters; (d) avoids the pitfalls associated with binning data; and (e) is applicable to virtually any resolved young cluster, laying the groundwork for a systematic study of the high mass stellar MF (M > 1 Msun). Using simulated clusters and Markov chain Monte Carlo sampling of the probability distribution functions, we show that estimates of the MF slope, {\alpha}, are unbiased and that the uncertainty, {\Delta}{\alpha}, depends primarily on the number of observed stars and stellar mass range they span, assuming that the uncertainties on individual masses and the completeness are well-characterized. Using idealized mock data, we compute the lower limit precision on {\alpha} and provide an analytic approximation for {\Delta}{\alpha} as a function of the observed number of stars and mass range. We find that ~ 3/4 of quoted literature uncertainties are smaller than the theoretical lower limit. By correcting these uncertainties to the theoretical lower limits, we find the literature studies yield <{\alpha}>=2.46 with a 1-{\sigma} dispersion of 0.35 dex. We verify that it is impossible for a power-law MF to obtain meaningful constraints on the upper mass limit of the IMF. We show that avoiding substantial biases in the MF slope requires: (1) including the MF as a prior when deriving individual stellar mass estimates; (2) modeling the uncertainties in the individual stellar masses; and (3) fully characterizing and then explicitly modeling the completeness for stars of a given mass. (abridged)

We study the emission by dust and stars in the Large and Small Magellanic Clouds, a pair of low-metallicity nearby galaxies, as traced by their spatially resolved spectral energy distributions (SEDs). This project combines Herschel Space Observatory PACS and SPIRE far-infrared photometry with other data at infrared and optical wavelengths. We build maps of dust and stellar luminosity and mass of both Magellanic Clouds, and analyze the spatial distribution of dust/stellar luminosity and mass ratios. These ratios vary considerably throughout the galaxies, generally between the range $0.01\leq L_{\rm dust}/L_\ast\leq 0.6$ and $10^{-4}\leq M_{\rm dust}/M_\ast\leq 4\times10^{-3}$. We observe that the dust/stellar ratios depend on the interstellar medium (ISM) environment, such as the distance from currently or previously star-forming regions, and on the intensity of the interstellar radiation field (ISRF). In addition, we construct star formation rate (SFR) maps, and find that the SFR is correlated with the dust/stellar luminosity and dust temperature in both galaxies, demonstrating the relation between star formation, dust emission and heating, though these correlations exhibit substantial scatter.

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

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.

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local Universe (12 + log(O/H) $=$ 7.17). In this work we study I Zw 18 using data from {\it Spitzer}, {\it Herschel Space Telescope} and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive maps of I Zw 18, to date, in both, the far infrared and the CO $J=1\rightarrow0$ transition. We use dust emission models to derive a dust mass upper limit of only M$_{dust}\leq1.1\times10^4$ M$_{\odot}$ ($3\sigma$ limit). This upper limit is driven by the non-detection at 160 $\mu$m, and it is a factor of 4-10 times smaller than previous estimates (depending upon the model used). We also estimate an upper limit to the total dust-to-gas mass ratio of M$_{Dust}$/M$_{gas}\leq5.0\times10^{-5}$. If a linear correlation between the dust-to-gas mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9$\times10^{-4}$. We also show that the infrared SED is similar to that of starbursting systems.

We present Herschel/SPIRE images at 250, 350, and 500 {\mu}m of NGC 4244, a typical low-mass, disk-only and edge-on spiral galaxy. The dust disk is clumpy and shows signs of truncation at the break radius of the stellar disk. This disk coincides with the densest part of the Hi disk. We compare the Spectral Energy Distribution, including the new SPIRE fluxes, to 3D radiative transfer models; a smooth model disk and a clumpy model with embedded heating. Each model requires a very high value for the dust scale-length (h(dust) = 2 - 5 h(stars)), higher dust masses than previous models of NGC 4244 (Md = 0.47 - 1.39 \times 10e7 Msun) and a face-on optical depth of {\tau}(V) = 0.4 - 1.12, in agreement with previous disk opacity studies. The vertical scales of stars and dust are similar. The clumpy model much better mimics the general morphology in the submm images and the general SED. The inferred gas-to-dust mass ratio is compatible with those of similar low-mass disks. The relatively large radial scale-length of the dust disk points to radial mixing of the dusty ISM within the stellar disk. The large vertical dust scale and the clumpy dust distribution of our SED model are both consistent with a scenario in which the vertical structure of the ISM is dictated by the balance of turbulence and self-gravity.

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 Panchromatic Hubble Andromeda Treasury (PHAT) is an on-going HST Multicycle Treasury program to image ~1/3 of M31's star forming disk in 6 filters, from the UV to the NIR. The full survey will resolve the galaxy into more than 100 million stars with projected radii from 0-20 kpc over a contiguous 0.5 square degree area in 828 orbits, producing imaging in the F275W and F336W filters with WFC3/UVIS, F475W and F814W with ACS/WFC, and F110W and F160W with WFC3/IR. The resulting wavelength coverage gives excellent constraints on stellar temperature, bolometric luminosity, and extinction for most spectral types. The photometry reaches SNR=4 at F275W=25.1, F336W=24.9, F475W=27.9, F814W=27.1, F110W=25.5, and F160W=24.6 for single pointings in the uncrowded outer disk; however, the optical and NIR data are crowding limited, and the deepest reliable magnitudes are up to 5 magnitudes brighter in the inner bulge. All pointings are dithered and produce Nyquist-sampled images in F475W, F814W, and F160W. We describe the observing strategy, photometry, astrometry, and data products, along with extensive tests of photometric stability, crowding errors, spatially-dependent photometric biases, and telescope pointing control. We report on initial fits to the structure of M31's disk, derived from the density of RGB stars, in a way that is independent of the assumed M/L and is robust to variations in dust extinction. These fits also show that the 10 kpc ring is not just a region of enhanced recent star formation, but is instead a dynamical structure containing a significant overdensity of stars with ages >1 Gyr. (Abridged)

Gamma-ray burst (GRB) host galaxies have been studied extensively in optical photometry and spectroscopy. Here we present the first mid-infrared spectrum of a GRB host, HG031203. It is one of the nearest GRB hosts at z=0.1055, allowing both low and high-resolution spectroscopy with Spitzer-IRS. Medium resolution UV-to-K-band spectroscopy with the X-shooter spectrograph on the VLT is also presented, along with Spitzer IRAC and MIPS photometry, as well as radio and sub-mm observations. These data allow us to construct a UV-to-radio spectral energy distribution with almost complete spectroscopic coverage from 0.3-35 micron of a GRB host galaxy for the first time, potentially valuable as a template for future model comparisons. The IRS spectra show strong, high-ionisation fine structure line emission indicative of a hard radiation field in the galaxy, suggestive of strong ongoing star-formation and a very young stellar population. The selection of HG031203 via the presence of a GRB suggests that it might be a useful analogue of very young star-forming galaxies in the early universe, and hints that local BCDs may be used as more reliable analogues of star-formation in the early universe than typical local starbursts. We look at the current debate on the ages of the dominant stellar populations in z~7 and z~8 galaxies in this context. The nebular line emission is so strong in HG031203, that at z~7, it can reproduce the spectral energy distributions of z-band dropout galaxies with elevated IRAC 3.6 and 4.5 micron fluxes without the need to invoke a 4000A break.

We estimate the total dust input from the cool evolved stars in the Small Magellanic Cloud (SMC), using the 8 micron excess emission as a proxy for the dust-production rate. We find that Asymptotic Giant Branch (AGB) and red supergiant (RSG) stars produce (8.6-9.5) x 10^7 solar masses per year of dust, depending on the fraction of far-infrared sources that belong to the evolved star population (with 10%-50% uncertainty in individual dust-production rates). RSGs contribute the least (<4%), while carbon-rich AGB stars (especially the so-called "extreme" AGB stars) account for 87%-89% of the total dust input from cool evolved stars. We also estimate the dust input from hot stars and supernovae (SNe), and find that if SNe produce 10^-3 solar masses of dust each, then the total SN dust input and AGB input are roughly equivalent. We consider several scenarios of SNe dust production and destruction and find that the interstellar medium (ISM) dust can be accounted for solely by stellar sources if all SNe produce dust in the quantities seen around the dustiest examples and if most SNe explode in dense regions where much of the ISM dust is shielded from the shocks. We find that AGB stars contribute only 2.1% of the ISM dust. Without a net positive contribution from SNe to the dust budget, this suggests that dust must grow in the ISM or be formed by another unknown mechanism.

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 present a characterization of the close environment of GRB980425 based on 5-160mic spectro-imaging obtained with Spitzer. The Gamma-Ray Burst GRB980425 occurred in a nearby (z=0.0085) SBc-type dwarf galaxy, at a projected distance of 900pc from an HII region with strong signatures of Wolf-Rayet (WR) stars. While this "WR region" produces less than 5% of the B-band emission of the host, we find that it is responsible for 45+/-10% of the total infrared luminosity, with a maximum contribution reaching 75% at 25-30mic. This atypical property is rarely observed among morphologically-relaxed dwarves, suggesting a strong causal link with the GRB event. The luminosity of the WR region (L_8-1000mic=4.6x10^8 Lsol), the peak of its spectral energy distribution at <~100mic and the presence of highly-ionized emission lines (e.g., [NeIII]) also reveal extremely young (<5Myr) star-forming activity, with a typical time-scale of only 47Myr to double the stellar mass already built. Finally, the mid-IR over B-band luminosity ratio in this region is substantially higher than in star-forming galaxies with similar L_IR, but it is lower than in young dust-enshrouded stellar clusters. Considering the modest obscuration measured from the silicate features (tau_9.7mic ~ 0.015), this suggests that the WR region is dominated by one or several star clusters that have either partly escaped or cleared out their parent molecular cloud. Combined with the properties characterizing the whole population of GRB hosts, our results reinforce the idea that long GRBs mostly happen within or in the vicinity of relatively unobscured galactic regions harboring very recent star formation.

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger.

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.

Edge-on spiral galaxies offer a unique perspective on the vertical structure of spiral disks, both stars and the iconic dark dustlanes. The thickness of these dustlanes can now be resolved for the first time with Herschel in far-infrared and sub-mm emission. We present NHEMESES, an ongoing project that targets 12 edge-on spiral galaxies with the PACS and SPIRE instruments on Herschel. These vertically resolved observations of edge-on spirals will impact on several current topics. First and foremost, these Herschel observations will settle whether or not there is a phase change in the vertical structure of the ISM with disk mass. Previously, a dramatic change in dustlane morphology was observed as in massive disks the dust collapses into a thin lane. If this is the case, the vertical balance between turbulence and gravity dictates the ISM structure and consequently star-formation and related phenomena (spiral arms, bars etc.). We specifically target lower mass nearby edge-ons to complement existing Herschel observations of high-mass edge-on spirals (the HEROES project). Secondly, the combined data-set, together with existing Spitzer observations, will drive a new generation of spiral disk Spectral Energy Distribution models. These model how dust reprocesses starlight to thermal emission but the dust geometry remains the critical unknown. And thirdly, the observations will provide an accurate and unbiased census of the cold dusty structures occasionally seen extending out of the plane of the disk, when backlit by the stellar disk. To illustrate the NHEMESES project, we present early results on NGC 4244 and NGC 891, two well studies examples of a low and high-mass edge-on spiral.

Cosmic rays are a known problem in astronomy, causing both loss of data and data inaccuracy. The problem becomes even more extreme when considering data from a high-radiation environment, such as in orbit around Earth or outside the Earth's magnetic field altogether, unprotected, as will be the case for the James Webb Space Telescope (JWST). For JWST, all the instruments employ nondestructive readout schemes. The most common of these will be "up the ramp" sampling, where the detector is read out regularly during the ramp. We study three methods to correct for cosmic rays in these ramps: a two-point difference method, a deviation from the fit method, and a y-intercept method. We apply these methods to simulated nondestructive read ramps with single-sample groups and varying combinations of flux, number of samples, number of cosmic rays, cosmic-ray location in the exposure, and cosmic-ray strength. We show that the y-intercept method is the optimal detection method in the read-noise-dominated regime, while both the y-intercept method and the two-point difference method are best in the photon-noise-dominated regime, with the latter requiring fewer computations.

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

The SAGE-Spectroscopy Spitzer legacy program (SAGE-Spec; PI: F. Kemper, PID: 40159; Kemper et al. 2010) is the IRS spectroscopic follow-up to the successful SAGE-LMC legacy program (Meixner et al. 2006; PI: M. Meixner, PID: 20203) that mapped the Large Magellanic Cloud (LMC) with all bands of the IRAC and MIPS instruments on board the Spitzer Space Telescope. This technical document gives details of the data-reduction procedure and the various data products that are publicly available through IRSA: http://irsa.ipac.caltech.edu/data/SPITZER/SAGE/

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers.

We report far-infrared and submillimeter observations of Supernova 1987A, the star that exploded on February 23, 1987 in the Large Magellanic Cloud, a galaxy located 160,000 lightyears away. The observations reveal the presence of a population of cold dust grains radiating with a temperature of ~17-23 K at a rate of about 220 solar luminosity. The intensity and spectral energy distribution of the emission suggests a dust mass of ~0.4-0.7 solar mass. The radiation must originate from the SN ejecta and requires the efficient precipitation of all refractory material into dust. Our observations imply that supernovae can produce the large dust masses detected in young galaxies at very high redshifts.

We compare atomic gas, molecular gas, and the recent star formation rate (SFR) inferred from H-alpha in the Small Magellanic Cloud (SMC). By using infrared dust emission and local dust-to-gas ratios, we construct a map of molecular gas that is independent of CO emission. This allows us to disentangle conversion factor effects from the impact of metallicity on the formation and star formation efficiency of molecular gas. On scales of 200 pc to 1 kpc we find a characteristic molecular gas depletion time of ~1.6 Gyr, similar to that observed in the molecule-rich parts of large spiral galaxies on similar spatial scales. This depletion time shortens on much larger scales to ~0.6 Gyr because of the presence of a diffuse H-alpha component, and lengthens on much smaller scales to ~7.5 Gyr because the H-alpha and H2 distributions differ in detail. We estimate the systematic uncertainties in our measurement to be a factor of 2-3. We suggest that the impact of metallicity on the physics of star formation in molecular gas has at most this magnitude. The relation between SFR and neutral (H2+HI) gas surface density is steep, with a power-law index ~2.2+/-0.1, similar to that observed in the outer disks of large spiral galaxies. At a fixed total gas surface density the SMC has a 5-10 times lower molecular gas fraction (and star formation rate) than large spiral galaxies. We explore the ability of the recent models by Krumholz et al. (2009) and Ostriker et al. (2010) to reproduce our observations. We find that to explain our data at all spatial scales requires a low fraction of cold, gravitationally-bound gas in the SMC. We explore a combined model that incorporates both large scale thermal and dynamical equilibrium and cloud-scale photodissociation region structure and find that it reproduces our data well, as well as predicting a fraction of cold atomic gas very similar to that observed in the SMC.

We investigate the infrared (IR) properties of cool, evolved stars in the Small Magellanic Cloud (SMC), including the red giant branch (RGB) stars and the dust-producing red supergiant (RSG) and asymptotic giant branch (AGB) stars using observations from the Spitzer Space Telescope Legacy program entitled: "Surveying the Agents of Galaxy Evolution in the Tidally-stripped, Low Metallicity SMC", or SAGE-SMC. The survey includes, for the first time, full spatial coverage of the SMC bar, wing, and tail regions at infrared (IR) wavelengths (3.6 - 160 microns). We identify evolved stars using a combination of near-IR and mid-IR photometry and point out a new feature in the mid-IR color-magnitude diagram that may be due to particularly dusty O-rich AGB stars. We find that the RSG and AGB stars each contribute ~20% of the global SMC flux (extended + point-source) at 3.6 microns, which emphasizes the importance of both stellar types to the integrated flux of distant metal-poor galaxies. The equivalent SAGE survey of the higher-metallicity Large Magellanic Cloud (SAGE-LMC) allows us to explore the influence of metallicity on dust production. We find that the SMC RSG stars are less likely to produce a large amount of dust (as indicated by the [3.6]-[8] color). There is a higher fraction of carbon-rich stars in the SMC, and these stars appear to able to reach colors as red as their LMC counterparts, indicating that C-rich dust forms efficiently in both galaxies. A preliminary estimate of the dust production in AGB and RSG stars reveals that the extreme C-rich AGB stars dominate the dust input in both galaxies, and that the O-rich stars may play a larger role in the LMC than in the SMC.

In star forming galaxies, dust plays a significant role in shaping the ultraviolet (UV) through infrared (IR) spectrum. Dust attenuates the radiation from stars, and re-radiates the energy through equilibrium and non-equilibrium emission. Polycyclic aromatic hydrocarbons (PAH), graphite, and silicates contribute to different features in the spectral energy distribution; however, they are all highly opaque in the same spectral region - the UV. Compared to old stellar populations, young populations release a higher fraction of their total luminosity in the UV, making them a good source of the energetic UV photons that can power dust emission. However, given their relative abundance, the question of whether young or old stellar populations provide most of these photons that power the infrared emission is an interesting question. Using three samples of galaxies observed with the Spitzer Space Telescope and our dusty radiative transfer model, we find that young stellar populations (on the order of 100 million years old) dominate the dust heating in star forming galaxies, and old stellar populations (13 billion years old) generally contribute less than 20% of the far-IR luminosity.

We present an analysis of the stellar kinematics of the Large Magellanic Cloud based on ~5900 new and existing velocities of massive red supergiants, oxygen-rich and carbon-rich AGB stars, and other giants. After correcting the line-of-sight velocities for the LMC's space motion and accounting for asymmetric drift in the AGB population, we derive a rotation curve that is consistent with all of the tracers used, as well as that of published HI data. The amplitude of the rotation curve is v_0=87+/-5 km s^-1 beyond a radius R_0=2.4+/-0.1 kpc, and has a position angle of the kinematic line of nodes of theta=142 degrees +/-5 degrees. By examining the outliers from our fits, we identify a population of 376 stars, or >~5% of our sample, that have line-of-sight velocities that apparently oppose the sense of rotation of the LMC disk. We find that these kinematically distinct stars are either counter-rotating in a plane closely aligned with the LMC disk, or rotating in the same sense as the LMC disk, but in a plane that is inclined by 54 degrees +/- 2 degrees to the LMC. Their kinematics clearly link them to two known HI arms, which have previously been interpreted as being pulled out from the LMC. We measure metallicities from the Ca triplet lines of ~1000 LMC field stars and 30 stars in the kinematically distinct population. For the LMC field, we find a median [Fe/H]=-0.56 +/- 0.02 with dispersion of 0.5 dex, while for the kinematically distinct stars the median [Fe/H] is -1.25 +/- 0.13 with a dispersion of 0.7 dex. The metallicity differences provide strong evidence that the kinematically distinct population originated in the SMC. This interpretation has the consequence that the HI arms kinematically associated with the stars are likely falling into the LMC, instead of being pulled out.

We present a galaxy (SAGE1CJ053634.78-722658.5) at a redshift of 0.14 of which the IR is entirely dominated by emission associated with the AGN. We present the 5-37 um Spitzer/IRS spectrum and broad wavelength SED of SAGE1CJ053634, an IR point-source detected by Spitzer/SAGE (Meixner et al 2006). The source was observed in the SAGE-Spec program (Kemper et al., 2010) and was included to determine the nature of sources with deviant IR colours. The spectrum shows a redshifted (z=0.14+-0.005) silicate emission feature with an exceptionally high feature-to-continuum ratio and weak polycyclic aromatic hydrocarbon (PAH) bands. We compare the source with models of emission from dusty tori around AGNs from Nenkova et al. (2008). We present a diagnostic diagram that will help to identify similar sources based on Spitzer/MIPS and Herschel/PACS photometry. The SED of SAGE1CJ053634 is peculiar because it lacks far-IR emission and a clear stellar counterpart. We find that the SED and the IR spectrum can be understood as emission originating from the inner ~10 pc around an accreting black hole. There is no need to invoke emission from the host galaxy, either from the stars or from the interstellar medium, although a possible early-type host galaxy cannot be excluded based on the SED analysis. The hot dust around the accretion disk gives rise to a continuum, which peaks at 4 um, whereas the strong silicate features may arise from optically thin emission of dusty clouds within ~10 pc around the black hole. The weak PAH emission does not appear to be linked to star formation, as star formation templates strongly over-predict the measured far-IR flux levels. The SED of SAGE1CJ053634 is rare in the local universe but may be more common in the more distant universe. The conspicuous absence of host-galaxy IR emission places limits on the far-IR emission arising from the dusty torus alone.

The absolute flux calibration of the James Webb Space Telescope will be based on a set of stars observed by the Hubble and Spitzer Space Telescopes. In order to cross-calibrate the two facilities, several A, G, and white dwarf (WD) stars are observed with both Spitzer and Hubble and are the prototypes for a set of JWST calibration standards. The flux calibration constants for the four Spitzer IRAC bands 1-4 are derived from these stars and are 2.3, 1.9, 2.0, and 0.5% lower than the official cold-mission IRAC calibration of Reach et al. (2005), i.e. in agreement within their estimated errors of ~2%. The causes of these differences lie primarily in the IRAC data reduction and secondarily in the SEDs of our standard stars. The independent IRAC 8 micron band-4 fluxes of Rieke et al. (2008) are about 1.5 +/- 2% higher than those of Reach et al. and are also in agreement with our 8 micron result.

We present a study of the dust properties of 12 classical and superbubble HII regions in the Large Magellanic Cloud. We use infrared photometry from Spitzer (8, 24, 70, and 160 \mum bands), obtained as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) program, along with archival spectroscopic classifications of the ionizing stars to examine the role of stellar sources on dust heating and processing. Our infrared observations show surprisingly little correlation between the emission properties of the dust and the effective temperatures or bolometric magnitudes of stars in the HII regions, suggesting that the HII region evolutionary timescale is not on the order of the dust processing timescale. We find that the infrared emission of superbubbles and classical HII regions shows little differentiation between the two classes, despite the significant differences in age and morphology. We do detect a correlation of the 24 \mum emission from hot dust with the ratio of 70 to 160 \mum flux. This correlation can be modeled as a trend in the temperature of a minority hot dust component, while a majority of the dust remains significantly cooler.

We estimate the conversion factor relating CO emission to H2 mass, alpha_CO, in five Local Group galaxies that span approximately an order of magnitude in metallicity - M31, M 33, the Large Magellanic Cloud (LMC), NGC 6822, and the Small Magellanic Cloud (SMC). We model the dust mass along the line of sight from infrared (IR) emission and then solve for the alpha_CO that best allows a single gas-to-dust ratio (delta_GDR) to describe each system. This approach remains sensitive to CO-dark envelopes of H2 surrounding molecular clouds. In M 31, M 33, and the LMC we find alpha_CO \approx 3-9 M_sun pc^-2 (K km s^-1)^-1, consistent with the Milky Way value within the uncertainties. The two lowest metallicity galaxies in our sample, NGC 6822 and the SMC (12 + log(O/H) \approx 8.2 and 8.0), exhibit a much higher alpha_CO. Our best estimates are \alpha_NGC6822 \approx 30 M_sun/pc^-2 (K km s^-1)^-1 and \alpha_SMC \approx 70 M_sun/pc^-2 (K km s-1)-1. These results are consistent with the conversion factor becoming CO a strong function of metallicity around 12 + log(O/H) \sim 8.4 - 8.2. We favor an interpretation where decreased dust-shielding leads to the dominance of CO-free envelopes around molecular clouds below this metallicity.

Fundamental parameters and time-evolution of mass loss are investigated for post-main-sequence stars in the Galactic globular cluster 47 Tucanae (NGC 104). This is accomplished by fitting spectral energy distributions (SEDs) to existing optical and infrared photometry and spectroscopy, to produce a true Hertzsprung--Russell diagram. We confirm the cluster's distance as 4611 (+213, -200) pc and age as 12 +/- 1 Gyr. Horizontal branch models appear to confirm that no more RGB mass loss occurs in 47 Tuc than in the more-metal-poor omega Centauri, though difficulties arise due to inconsistencies between the models. Using our SEDs, we identify those stars which exhibit infrared excess, finding excess only among the brightest giants: dusty mass loss begins at a luminosity of ~ 1000 Lsun, becoming ubiquitous above 2000 Lsun. Recent claims of dust production around lower-luminosity giants cannot be reproduced, despite using the same archival Spitzer imagery.

We present the photometric catalogs for the star-forming cluster NGC 602 in the wing of the Small Magellanic Cloud covering a range of wavelengths from optical HST/ACS (F555W, F814W) and SMARTS/ANDICAM (V, I) to infrared (Spitzer/IRAC 3.6, 4.5, 5.8, and 8 micron and MIPS 24 micron). Combining this with IRSF (InfraRed Survey Facility) near-infrared photometry (J, H, Ks), we compare the young main sequence (MS) and pre-main sequence (PMS) populations prominent in the optical with the current young stellar object (YSO) populations revealed by the infrared (IR). We analyze the MS and PMS population with isochrones in color-magnitude diagrams to derive ages and masses. The optical data reveal ~565 PMS candidates, low mass Stage III YSOs. We characterize ~40 YSOs by fitting their spectral energy distributions (SEDs) to a grid of models (Robitaille et al. 2007) to derive luminosities, masses and evolutionary phase (Stage I-III). The higher resolution HST images reveal that ~70% of the YSO candidates are either multiples or protoclusters. For YSOs and PMS sources found in common, we find a consistency in the masses derived. We use the YSO mass function to derive a present-day star-formation rate of ~0.2-1.0 Msun/yr/kpc^2, similar to the rate derived from the optical star formation history suggesting a constant star formation rate for this region. We demonstrate a progression of star formation from the optical star cluster center to the edge of the star forming dust cloud. We derive lifetimes of a few 10^5 years for the YSO Stages I and II.

We present the largest sample to-date of intermediate-resolution blue-to-red optical spectra of B-type supergiants in M31 and undertake the first survey of diffuse interstellar bands (DIBs) in this galaxy. Spectral classifications, radial velocities and interstellar reddenings are presented for 34 stars in three regions of M31. Radial velocities and equivalent widths are given for the 5780 and 6283 DIBs towards 11 stars. Equivalent widths are also presented for the following DIBs detected in three sightlines in M31: 4428, 5705, 5780, 5797, 6203, 6269, 6283, 6379, 6613, 6660, and 6993. All of these M31 DIB carriers reside in clouds at radial velocities matching those of interstellar Na I and/or H I. The relationships between DIB equivalent widths and reddening (E(B-V)) are consistent with those observed in the local ISM of the Milky Way. Many of the observed sightlines show DIB strengths (per unit reddening) which lie at the upper end of the range of Galactic values. DIB strengths per unit reddening are found (with 68% confidence), to correlate with the interstellar UV radiation field strength. The strongest DIBs are observed where the interstellar UV flux is lowest. The mean Spitzer 8/24 micron emission ratio in our three fields is slightly lower than that measured in the Milky Way, but we identify no correlation between this ratio and the DIB strengths in M31. Interstellar oxygen abundances derived from the spectra of three M31 H II regions in one of the fields indicate that the average metallicity of the ISM in that region is 12 + log[O/H] = 8.54 +- 0.18, which is approximately equal to the value in the solar neighbourhood.

We present a study of the infrared properties of 4922 spectroscopically confirmed massive stars in the Large and Small Magellanic Clouds, focusing on the active OB star population. Besides OB stars, our sample includes yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables (LBVs) and supergiant B[e] stars. We detect a distinct Be star sequence, displaced to the red, and find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We also find photometric variability among the active OB population and evidence for transitions of Be stars to B stars and vice versa. We furthermore confirm the presence of dust around all the supergiant B[e] stars in our sample, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs.

We present results of our study of the infrared properties of massive stars in the Large and Small Magellanic Clouds, which are based on the Spitzer SAGE surveys of these galaxies. We have compiled catalogs of spectroscopically confirmed massive stars in each galaxy, as well as photometric catalogs for a subset of these stars that have infrared counterparts in the SAGE database, with uniform photometry from 0.3 to 24 microns in the UBVIJHKs+IRAC+MIPS24 bands. These catalogs enable a comparative study of infrared excesses of OB stars, classical Be stars, yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables and supergiant B[e] stars, as a function of metallicity, and provide the first roadmaps for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths.

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

We present a 2Dust model for the dust shell around a LMC long-period variable (LPV) previously studied as part of the OGLE survey. OGLE LMC LPV 28579 (SAGE J051306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have photometry and spectra from the Spitzer SAGE and SAGE-Spec programs along with UBVIJHK_s photometry. By modeling this source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce its spectral energy distribution using a mixture of AmC and SiC (15% by mass). The grain sizes are distributed according to the KMH model. The best-fit model has an optical depth of 0.28 for the shell at the peak of the SiC feature, with R_in~1430 R_sun or 4.4 R_star. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s^-1, we obtain a dust mass-loss rate of 2.5x10^-9 M_sun yr-1. We calculate a 15% variation in this rate by testing the fit sensitivity against variation in input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 13.7 \mu m feature seen in the spectrum. We find that a combination of CO and C_2H_2 gas at an excitation temperature of about 1000 K and column densities of 3x10^21 cm^-2 and 10^19 cm^-2 respectively are able to reproduce the observations. Given that the excitation temperature is close to T_dust(R_in), most of the molecular contribution probably arises from the inner shell region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L_sun. For an effective temperature of about 3000 K, this implies a stellar mass of 1.5-2 M_sun and an age of 1-2.5 Gyr. For a gas:dust ratio of 200, we obtain a gas mass-loss rate of 5.0x10^-7 M_sun yr^-1, consistent with the gas mass-loss rates estimated from the period, color and 8 \mu m flux of the source.

The amount of molecular gas is a key for understanding the future star formation in a galaxy. Because H2 is difficult to observe directly in dense and cold clouds, tracers like CO are used. However, at low metallicities especially, CO only traces the shielded interiors of the clouds. mm dust emission can be used as a tracer to unveil the total dense gas masses. The comparison of masses deduced from the continuum SIMBA 1.2 mm emission and virial masses in a sample of giant molecular clouds (GMCs), in the SW region of the Small Magellanic Cloud (SMC), showed a discrepancy that is in need of an explanation. This study aims at better assessing possible uncertainties on the dust emission observed in the sample of GMCs from the SMC and focuses on the densest parts of the GMCs where CO is detected. New observations were obtained with the LABOCA camera on the APEX telescope. All GMCs previously observed in CO are detected and their emission at 870microns is compared to ancillary data. The different contributions to the sub-mmm emission are estimated, as well as dust properties, in order to deduce molecular cloud masses precisely. The (sub-)mm emission observed in the GMCs in the SW region of the SMC is dominated by dust emission and masses are deduced for the part of each cloud where CO is detected and compared to the virial masses. The mass discrepancy between both methods is confirmed at 870microns with the LABOCA observations: the virial masses are on average 4 times smaller than the masses of dense gas deduced from dust emission, contrary to what is observed for equivalent clouds in our Galaxy. The origin of this mass discrepancy in the SMC remains unkown. The direct interpretation of this effect is that the CO linewidth used to compute virial masses do not measure the full velocity distribution of the gas. Geometrical effects and uncertainties on the dust properties are also discussed.

We present a catalog of 5324 massive stars in the Small Magellanic Cloud (SMC), with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 3654 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer, SAGE-SMC survey database, for which we present uniform photometry from 0.3-24 um in the UBVIJHKs+IRAC+MIPS24 bands. We compare the color magnitude diagrams and color-color diagrams to those of the Large Magellanic Cloud (LMC), finding that the brightest infrared sources in the SMC are also the red supergiants, supergiant B[e] (sgB[e]) stars, luminous blue variables, and Wolf-Rayet stars, with the latter exhibiting less infrared excess, the red supergiants being less dusty and the sgB[e] stars being on average less luminous. Among the objects detected at 24 um are a few very luminous hypergiants, 4 B-type stars with peculiar, flat spectral energy distributions, and all 3 known luminous blue variables. We detect a distinct Be star sequence, displaced to the red, and suggest a novel method of confirming Be star candidates photometrically. We find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We estimate mass-loss rates for the red supergiants, confirming the correlation with luminosity even at the metallicity of the SMC. Finally, we confirm the new class of stars displaying composite A & F type spectra, the sgB[e] nature of 2dFS1804 and find the F0 supergiant 2dFS3528 to be a candidate luminous blue variable with cold dust.

Recently, we have discovered an error in our Monte-Carlo spectral fitting routine, more specifically where the errors on the fluxes were rescaled to get a reduced chi2 of 1. The rescaled errors were too big, resulting in too wide a range of good fits in our 100 step Monte-Carlo routine. This problem affects Figs. 7-9 and Tables A.1, A.2 in Gielen et al. (2008), Table 3 in Gielen et al. (2009a), and Table 4 in Gielen et al. (2009b). We corrected for this error and present the new values and errors in the tables below. The new values and errors nearly all fall within the old error range. Our best chi2 values and overall former scientific results are not affected. With these new errors some possible new trends in the dust parameters might be observed. These will be discussed in an upcoming paper where we extend the sample presented in Gielen et al. (2008) with newly obtained SPITZER-IRS data.

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 {\mu}m and SPIRE 250, 350 and 500 {\mu}m. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 {\mu}m emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

Massive evolved stars can produce large amounts of dust, and far-infrared (IR) data are essential for determining the contribution of cold dust to the total dust mass. Using Herschel, we search for cold dust in three very dusty massive evolved stars in the Large Magellanic Cloud: R71 is a Luminous Blue Variable, HD36402 is a Wolf-Rayet triple system, and IRAS05280-6910 is a red supergiant. We model the spectral energy distributions using radiative transfer codes and find that these three stars have mass-loss rates up to 10^-3 solar masses/year, suggesting that high-mass stars are important contributors to the life-cycle of dust. We found far-IR excesses in two objects, but these excesses appear to be associated with ISM and star-forming regions. Cold dust (T < 100 K) may thus not be an important contributor to the dust masses of evolved stars.

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.

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

We investigate the dust associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC) as observed with the Herschel Space Observatory. N49 is unusually bright because of an interaction with a molecular cloud along its eastern edge. We have used PACS and SPIRE to measure the far IR flux densities of the entire SNR and of a bright region on the eastern edge of the SNR where the SNR shock is encountering the molecular cloud. Using these fluxes supplemented with archival data at shorter wavelengths, we estimate the dust mass associated with N49 to be about 10 Msun. The bulk of the dust in our simple two-component model has a temperature of 20-30 K, similar to that of nearby molecular clouds. Unfortunately, as a result of the limited angular resolution of Herschel at the wavelengths sampled with SPIRE, the uncertainties are fairly large. Assuming this estimate of the dust mass associated with the SNR is approximately correct, it is probable that most of the dust in the SNR arises from regions where the shock speed is too low to produce significant X-ray emission. The total amount of warm 50-60 K dust is ~0.1 or 0.4 Msun, depending on whether the dust is modeled in terms of carbonaceous or silicate grains. This provides a firm lower limit to the amount of shock heated dust in N49.

We demonstrate the unique capabilities of Herschel to study very young luminous extragalactic young stellar objects (YSOs) by analyzing a central strip of the Large Magellanic Cloud obtained through the HERITAGE Science Demonstration Program. We combine PACS 100 and 160, and SPIRE 250, 350, and 500 microns photometry with 2MASS (1.25-2.17 microns) and Spitzer IRAC and MIPS (3.6-70 microns) to construct complete spectral energy distributions (SEDs) of compact sources. From these, we identify 207 candidate embedded YSOs in the observed region, ~40% never-before identified. We discuss their position in far-infrared color-magnitude space, comparing with previously studied, spectroscopically confirmed YSOs and maser emission. All have red colors indicating massive cool envelopes and great youth. We analyze four example YSOs, determining their physical properties by fitting their SEDs with radiative transfer models. Fitting full SEDs including the Herschel data requires us to increase the size and mass of envelopes included in the models. This implies higher accretion rates (greater than or equal to 0.0001 M_sun/yr), in agreement with previous outflow studies of high-mass protostars. Our results show that Herschel provides reliable longwave SEDs of large samples of high-mass YSOs; discovers the youngest YSOs whose SEDs peak in Herschel bands; and constrains the physical properties and evolutionary stages of YSOs more precisely than was previously possible.

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.

Spitzer and IRAS observations of the LMC suggest an excess of FIR emission with respect to the gas surface density traced by 12CO and HI 21 cm emission lines. This "FIR excess" is noticeable near molecular clouds in the LMC, and has usually been interpreted as the presence of a self-shielded H2 component not traced by CO molecular clouds' envelopes. Based on Herschel observations, we examine the correlation between gas and dust at higher resolution than previously achieved. We consider three additional causes for the FIR excess: X factor, FIR dust emissivity, and gas-to-dust ratio variations between the diffuse and dense phases of the ISM. We examine the structure of NT80 and NT71, two molecular clouds detected in the NANTEN 12CO survey of the LMC. Dust surface density maps are derived from the HERITAGE data. The gas phase is traced by MAGMA 12CO and ATCA HI 21 cm observations of the LMC. The dust emissivity, gas-to-dust ratio, and X factor required to match the dust and gas surface densities are derived, and their correlations with the dust surface density are examined. The dust surface density is spatially correlated with the atomic and molecular gas phases. The dust temperature is consistently lower in the dense phase of the ISM than in the diffuse phase. We confirm variations in the ratio of FIR emission to gas surface density derived from HI and CO observations. There is an excess of FIR emission, spatially correlated with regions of intermediate HI and dust surface densities (Av = 1-2), and little or no CO. While there is no significant trend in the dust emissivity or gas-to-dust ratio with dust surface density, the X factor is enhanced at Av = 1-2. We conclude that H2 envelopes not traced by CO and X factor variations close to the CO boundary are more likely to cause these deviations between FIR emission and gas surface density than gas-to-dust ratio or emissivity variations.

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.

We study the structure of the medium surrounding sites of high-mass star formation to determine the interrelation between the HII regions and the environment from which they were formed. The density distribution of the surroundings is key in determining how the radiation of the newly formed stars interacts with the surrounds in a way that allows it to be used as a star formation tracer. We present new Herschel/SPIRE 250, 350 and 500 mum data of LHA 120-N44 and LHA 120-N63 in the LMC. We construct average spectral energy distributions (SEDs) for annuli centered on the IR bright part of the star formation sites. The annuli cover ~10-~100 pc. We use a phenomenological dust model to fit these SEDs to derive the dust column densities, characterise the incident radiation field and the abundance of polycyclic aromatic hydrocarbon molecules. We see a factor 5 decrease in the radiation field energy density as a function of radial distance around N63. N44 does not show a systematic trend. We construct a simple geometrical model to derive the 3-D density profile of the surroundings of these two regions. Herschel/SPIRE data have proven very efficient in deriving the dust mass distribution. We find that the radiation field in the two sources behaves very differently. N63 is more or less spherically symmetric and the average radiation field drops with distance. N44 shows no systematic decrease of the radiation intensity which is probably due to the inhomogeneity of the surrounding molecular material and to the complex distribution of several star forming clusters in the region.

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.

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}$.

HII regions are the birth places of stars, and as such they provide the best measure of current star formation rates (SFRs) in galaxies. The close proximity of the Magellanic Clouds allows us to probe the nature of these star forming regions at small spatial scales. We aim to determine the monochromatic IR band that most accurately traces the bolometric IR flux (TIR), which can then be used to estimate an obscured SFR. We present the spatial analysis, via aperture/annulus photometry, of 16 LMC and 16 SMC HII region complexes using the Spitzer IRAC and MIPS bands. UV rocket data and SHASSA H-alpha data are also included. We find that nearly all of the LMC and SMC HII region SEDs peak around 70um, from ~10 to ~400 pc from the central sources. As a result, the sizes of HII regions as probed by 70um is approximately equal to the sizes as probed by TIR (about 70 pc in radius); the radial profile of the 70um flux, normalized by TIR, is constant at all radii (70um ~ 0.45 TIR); the 1-sigma standard deviation of the 70um fluxes, normalized by TIR, is a lower fraction of the mean (0.05 to 0.12 out to ~220 pc) than the normalized 8, 24, and 160um normalized fluxes (0.12 to 0.52); and these results are invariant between the LMC and SMC. From these results, we argue that 70um is the most suitable IR band to use as a monochromatic obscured star formation indicator because it most accurately reproduces the TIR of HII regions in the LMC and SMC and over large spatial scales. We also explore the general trends of the 8, 24, 70, and 160um bands in the LMC and SMC HII region SEDs, radial surface brightness profiles, sizes, and normalized (by TIR) radial flux profiles. We derive an obscured SFR equation that is modified from the literature to use 70um luminosity, SFR(Mo/yr) = 9.7(0.7)x10^{-44} L(70)(ergs/s), which is applicable from 10 to 300 pc distance from the center of an HII region.

Using infrared photometry from the Spitzer Space Telescope, we perform the first inventory of aromatic feature emission (AFE, but also commonly referred to as PAH emission) for a statistically complete sample of star-forming galaxies in the local volume. The photometric methodology involved is calibrated and demonstrated to recover the aromatic fraction of the IRAC 8 micron flux with a standard deviation of 6% for a training set of 40 SINGS galaxies (ranging from stellar to dust dominated) with both suitable mid-infrared Spitzer IRS spectra and equivalent photometry. A potential factor of two improvement could be realized with suitable 5.5 and 10 micron photometry, such as what may be provided in the future by JWST. The resulting technique is then applied to mid-infrared photometry for the 258 galaxies from the Local Volume Legacy (LVL) survey, a large sample dominated in number by low-luminosity dwarf galaxies for which obtaining comparable mid-infrared spectroscopy is not feasible. We find the total LVL luminosity due to five strong aromatic features in the 8 micron complex to be 2.47E10 solar luminosities with a mean volume density of 8.8E6 solar luminosities per cubic Megaparsec. Twenty-four of the LVL galaxies, corresponding to a luminosity cut at M = -18.22 in the B band, account for 90% of the aromatic luminosity. Using oxygen abundances compiled from the literature for 129 of the 258 LVL galaxies, we find a correlation between metallicity and the aromatic to total infrared emission ratio but not the aromatic to total 8 micron dust emission ratio. A possible explanation is that metallicity plays a role in the abundance of aromatic molecules relative to the total dust content, but other factors such as star formation and/or the local radiation field affect the excitation of those molecules.

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

(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.

Using Spitzer IRAC observations from the SAGE-SMC Legacy program and archived Spitzer IRAC data, we investigate dust production in 47 Tuc, a nearby massive Galactic globular cluster. A previous study detected infrared excess, indicative of circumstellar dust, in a large population of stars in 47 Tuc, spanning the entire Red Giant Branch (RGB). We show that those results suffered from effects caused by stellar blending and imaging artifacts and that it is likely that no stars below about 1 mag from the tip of the RGB are producing dust. The only stars that appear to harbor dust are variable stars, which are also the coolest and most luminous stars in the cluster.

We present 52-93 micron spectra, obtained with the Spitzer Space Telescope, of luminous compact far-IR sources in the SMC. These comprise 9 Young Stellar Objects (YSOs), the compact HII region N81 and a similar object within N84, and two red supergiants (RSGs). The spectra of the sources in N81 (of which we also show the ISO-LWS spectrum between 50-170 micron) and N84 both display strong [OI] 63-micron and [OIII] 88-micron fine-structure line emission. We attribute these lines to strong shocks and photo-ionized gas, respectively, in a ``champagne flow'' scenario. The nitrogen content of these two HII regions is very low, definitely N/O<0.04 but possibly as low as N/O<0.01. Overall, the oxygen lines and dust continuum are weaker in star-forming objects in the SMC than in the LMC. We attribute this to the lower metallicity of the SMC compared to that of the LMC. Whilst the dust mass differs in proportion to metallicity, the oxygen mass differs less; both observations can be reconciled with higher densities inside star-forming cloud cores in the SMC than in the LMC. The dust in the YSOs in the SMC is warmer (37-51 K) than in comparable objects in the LMC (32-44 K). We attribute this to the reduced shielding and reduced cooling at the low metallicity of the SMC. On the other hand, the efficiency of the photo-electric effect to heat the gas is found to be indistinguishable to that measured in the same manner in the LMC, 0.1-0.3%. This may result from higher cloud-core densities, or smaller grains, in the SMC. The dust associated with the two RSGs in our SMC sample is cool, and we argue that it is swept-up interstellar dust, or formed (or grew) within the bow-shock, rather than dust produced in these metal-poor RSGs themselves. Strong emission from crystalline water ice is detected in at least one YSO. (abridged)

We present spectroscopic observations of a sample of 15 embedded young stellar objects (YSOs) in the Large Magellanic Cloud (LMC). These observations were obtained with the Spitzer Infrared Spectrograph (IRS) as part of the SAGE-Spec Legacy program. We analyze the two prominent ice bands in the IRS spectral range: the bending mode of CO_2 ice at 15.2 micron and the ice band between 5 and 7 micron that includes contributions from the bending mode of water ice at 6 micron amongst other ice species. The 5-7 micron band is difficult to identify in our LMC sample due to the conspicuous presence of PAH emission superimposed onto the ice spectra. We identify water ice in the spectra of two sources; the spectrum of one of those sources also exhibits the 6.8 micron ice feature attributed to ammonium and methanol. We model the CO_2 band in detail, using the combination of laboratory ice profiles available in the literature. We find that a significant fraction (> 50%) of CO_2 ice is locked in a water-rich component, consistent with what is observed for Galactic sources. The majority of the sources in the LMC also require a pure-CO_2 contribution to the ice profile, evidence of thermal processing. There is a suggestion that CO_2 production might be enhanced in the LMC, but the size of the available sample precludes firmer conclusions. We place our results in the context of the star formation environment in the LMC.

Aims: By studying the photospheric abundances of 4 RV Tauri stars in the LMC, we test whether the depletion pattern of refractory elements, seen in similar Galactic sources, is also common for extragalactic sources. Since this depletion process probably only occurs through interaction with a stable disc, we investigate the circumstellar environment of these sources. Methods: A detailed photospheric abundance study was performed using high-resolution UVES optical spectra. To study the circumstellar environment we use photometric data to construct the spectral energy distributions of the stars, and determine the geometry of the circumstellar environment, whereas low-resolution Spitzer-IRS infrared spectra are used to trace its mineralogy. Results: Our results show that, also in the LMC, the photospheres of RV Tauri stars are commonly affected by the depletion process, although it can differ significantly in strength from source to source. From our detailed disc modelling and mineralogy study, we find that this process, as in the Galaxy, appears closely related to the presence of a stable Keplerian disc. The newly studied extragalactic objects have similar observational characteristics as Galactic post-AGB binaries surrounded by a dusty disc, and are therefore also believed to be part of a binary system. One source shows a very small infrared excess, atypical for a disc source, but still has evidence for depletion. We speculate this could point to the presence of a very evolved disc, similar to debris discs seen around young stellar objects.

[abridged] We present 52-93 micron spectra obtained with Spitzer in the MIPS-SED mode, of a representative sample of luminous compact far-IR sources in the LMC. These include carbon stars, OH/IR AGB stars, post-AGB objects and PNe, RCrB-type star HV2671, OH/IR red supergiants WOHG064 and IRAS05280-6910, B[e] stars IRAS04530-6916, R66 and R126, Wolf-Rayet star Brey3a, Luminous Blue Variable R71, supernova remnant N49, a large number of young stellar objects, compact HII regions and molecular cores, and a background galaxy (z~0.175). We use the spectra to constrain the presence and temperature of cold dust and the excitation conditions and shocks within the neutral and ionized gas, in the circumstellar environments and interfaces with the surrounding ISM. Evolved stars, including LBV R71, lack cold dust except in some cases where we argue that this is swept-up ISM. This leads to an estimate of the duration of the prolific dust-producing phase ("superwind") of several thousand years for both RSGs and massive AGB stars, with a similar fractional mass loss experienced despite the different masses. We tentatively detect line emission from neutral oxygen in the extreme RSG WOHG064, with implications for the wind driving. In N49, the shock between the supernova ejecta and ISM is revealed by its strong [OI] 63-micron emission and possibly water vapour; we estimate that 0.2 Msun of ISM dust was swept up. Some of the compact HII regions display pronounced [OIII] 88-micron emission. The efficiency of photo-electric heating in the interfaces of ionized gas and molecular clouds is estimated at 0.1-0.3%. We confirm earlier indications of a low nitrogen content in the LMC. Evidence for solid state emission features is found in both young and evolved object; some of the YSOs are found to contain crystalline water ice.

The Spectrograph for Photometric Imaging with Numeric Reconstruction (SPINR) sounding rocket experiment was launched on 2000 August 4 to record far-ultraviolet (912-1450 A) spectral and spatial information for the giant reflection nebula in the Upper Scorpius region. The data were divided into three arbitrary bandpasses (912-1029 A, 1030-1200 A, and 1235-1450 A) for which stellar and nebular flux levels were derived. These flux measurements were used to constrain a radiative transfer model and to determine the dust albedo for the Upper Scorpius region. The resulting albedos were 0.28+/-0.07 for the 912-1029 A bandpass, 0.33+/-0.07 for the 1030-1200 A bandpass, and 0.77+/-0.13 for the 1235-1450 A bandpass.

Compared to starburst galaxies, normal star forming galaxies have been shown to display a much larger dispersion of the dust attenuation at fixed reddening through studies of the IRX-beta diagram (the IR/UV ratio "IRX" versus the UV color "beta"). To investigate the causes of this larger dispersion and attempt to isolate second parameters, we have used GALEX UV, ground-based optical, and Spitzer infrared imaging of 8 nearby galaxies, and examined the properties of individual UV and 24 micron selected star forming regions. We concentrated on star-forming regions, in order to isolate simpler star formation histories than those that characterize whole galaxies. We find that 1) the dispersion is not correlated with the mean age of the stellar populations, 2) a range of dust geometries and dust extinction curves are the most likely causes for the observed dispersion in the IRX-beta diagram 3) together with some potential dilution of the most recent star-forming population by older unrelated bursts, at least in the case of star-forming regions within galaxies, 4) we also recover some general characteristics of the regions, including a tight positive correlation between the amount of dust attenuation and the metal content. Although generalizing our results to whole galaxies may not be immediate, the possibility of a range of dust extinction laws and geometries should be accounted for in the latter systems as well.

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

We investigate dust production and stellar mass loss in the Galactic globular cluster NGC 362. Due to its close proximity to the Small Magellanic Cloud (SMC), NGC 362 was imaged with the IRAC and MIPS cameras onboard the Spitzer Space Telescope as part of the Surveying the Agents of Galaxy Evolution (SAGE-SMC) Spitzer Legacy program. We detect several cluster members near the tip of the Red Giant Branch that exhibit infrared excesses indicative of circumstellar dust and find that dust is not present in measurable quantities in stars below the tip of the Red Giant Branch. We modeled the spectral energy distribution (SED) of the stars with the strongest IR excess and find a total cluster dust mass-loss rate of 3.0(+2.0/-1.2) x 10^-9 solar masses per year, corresponding to a gas mass-loss rate of 8.6(+5.6/-3.4) x 10^-6 solar masses per year, assuming [Fe/H] = -1.16. This mass loss is in addition to any dust-less mass loss that is certainly occurring within the cluster. The two most extreme stars, variables V2 and V16, contribute up to 45% of the total cluster dust-traced mass loss. The SEDs of the more moderate stars indicate the presence of silicate dust, as expected for low-mass, low-metallicity stars. Surprisingly, the SED shapes of the stars with the strongest mass-loss rates appear to require the presence of amorphous carbon dust, possibly in combination with silicate dust, despite their oxygen-rich nature. These results corroborate our previous findings in omega Centauri.

We present a sample of 75 extinction curves derived from FUSE far-ultraviolet spectra supplemented by existing IUE spectra. The extinction curves were created using the standard pair method based on a new set of dereddened FUSE+IUE comparison stars. Molecular hydrogen absorption features were removed using individualized H_2 models for each sightline. The general shape of the FUSE extinction (8.4 micron^-1 < lambda^-1 < 11 micron^-1) was found to be broadly consistent with extrapolations from the IUE extinction (3.3 micron-1 < lambda^-1 < 8.6 micron^-1) curve. Significant differences were seen in the strength of the far-UV rise and the width of the 2175 A bump. All the FUSE+IUE extinction curves had positive far-UV slopes giving no indication that the far-UV rise was turning over at the shortest wavelengths. The dependence of A(lambda)/A(V) versus R(V)^-1 in the far-UV using the sightlines in our sample was found to be stronger than tentatively indicated by previous work. We present an updated R(V) dependent relationship for the full UV wavelength range (3.3 micron^-1 <= lambda^-1 <= 11 micron^-1). Finally, we searched for discrete absorption features in the far-ultraviolet. We found a 3 sigma upper limit of ~0.12 A(V) on features with a resolution of 250 (~4 A width) and 3 sigma upper limits of ~0.15 A(V) for lambda^-1 < 9.6 micron^-1 and ~0.68 A(V) for lambda^-1 > 9.6 micron^-1 on features with a resolution of 10^4 (~0.1 A width).

The survey description and the near-, mid-, and far-infrared flux properties are presented for the 258 galaxies in the Local Volume Legacy (LVL). LVL is a Spitzer Space Telescope legacy program that surveys the local universe out to 11 Mpc, built upon a foundation of ultraviolet, H-alpha, and HST imaging from 11HUGS (11 Mpc H-alpha and Ultraviolet Galaxy Survey) and ANGST (ACS Nearby Galaxy Survey Treasury). LVL covers an unbiased, representative, and statistically robust sample of nearby star-forming galaxies, exploiting the highest extragalactic spatial resolution achievable with Spitzer. As a result of its approximately volume-limited nature, LVL augments previous Spitzer observations of present-day galaxies with improved sampling of the low-luminosity galaxy population. The collection of LVL galaxies shows a large spread in mid-infrared colors, likely due to the conspicuous deficiency of 8um PAH emission from low-metallicity, low-luminosity galaxies. Conversely, the far-infrared emission tightly tracks the total infrared emission, with a dispersion in their flux ratio of only 0.1 dex. In terms of the relation between infrared-to-ultraviolet ratio and ultraviolet spectral slope, the LVL sample shows redder colors and/or lower infrared-to-ultraviolet ratios than starburst galaxies, suggesting that reprocessing by dust is less important in the lower mass systems that dominate the LVL sample. Comparisons with theoretical models suggest that the amplitude of deviations from the relation found for starburst galaxies correlates with the age of the stellar populations that dominate the ultraviolet/optical luminosities.

We present a catalog of 1750 massive stars in the Large Magellanic Cloud, with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 1268 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer SAGE survey database, for which we present uniform photometry from 0.3-24 microns in the UBVIJHKs+IRAC+MIPS24 bands. The resulting infrared color-magnitude diagrams illustrate that the supergiant B[e], red supergiant and luminous blue variable (LBV) stars are among the brightest infrared point sources in the Large Magellanic Cloud, due to their intrinsic brightness, and at longer wavelengths, due to dust. We detect infrared excesses due to free-free emission among ~900 OB stars, which correlate with luminosity class. We confirm the presence of dust around 10 supergiant B[e] stars, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs. The similar luminosities of B[e] supergiants (log L/Lo>=4) and the rare, dusty progenitors of the new class of optical transients (e.g. SN 2008S and NGC 300 OT), plus the fact that dust is present in both types of objects, suggests a common origin for them. We find the infrared colors for Wolf-Rayet stars to be independent of spectral type and their SEDs to be flatter than what models predict. The results of this study provide the first comprehensive roadmap for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths.

To understand the impact of low metallicities on giant molecular cloud (GMC) structure, we compare far infrared dust emission, CO emission, and dynamics in the star-forming complex N83 in the Wing of the Small Magellanic Cloud. Dust emission (measured by Spitzer as part of the S3MC and SAGE-SMC surveys) probes the total gas column independent of molecular line emission and traces shielding from photodissociating radiation. We calibrate a method to estimate the dust column using only the high-resolution Spitzer data and verify that dust traces the ISM in the HI-dominated region around N83. This allows us to resolve the relative structures of H2, dust, and CO within a giant molecular cloud complex, one of the first times such a measurement has been made in a low-metallicity galaxy. Our results support the hypothesis that CO is photodissociated while H2 self-shields in the outer parts of low-metallicity GMCs, so that dust/self shielding is the primary factor determining the distribution of CO emission. Four pieces of evidence support this view. First, the CO-to-H2 conversion factor averaged over the whole cloud is very high 4-11 \times 10^21 cm^-2/(K km/s), or 20-55 times the Galactic value. Second, the CO-to-H2 conversion factor varies across the complex, with its lowest (most nearly Galactic) values near the CO peaks. Third, bright CO emission is largely confined to regions of relatively high line-of-sight extinction, A_V >~ 2 mag, in agreement with PDR models and Galactic observations. Fourth, a simple model in which CO emerges from a smaller sphere nested inside a larger cloud can roughly relate the H2 masses measured from CO kinematics and dust.

Using the data obtained with the Spitzer Space telescope as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) legacy survey, we have studied the variations of the dust composition and abundance across the Large Magellanic Cloud (LMC). Such variations are expected, as the explosive events which have lead to the formation of the many HI shells observed should have affected the dust properties. Using a model and comparing with a reference spectral energy distribution from our Galaxy, we deduce the relative abundance variations of small dust grains across the LMC. We examined the infrared color ratios as well as the relative abundances of very small grains (VSGs) and polycyclic aromatic hydrocarbons (PAHs) relative to the big grain (BG) abundance. Results show that each dust component could have different origins or evolution in the interstellar medium (ISM). The VSG abundance traces the star formation activity and could result from shattering of larger grains, whereas the PAH abundance increases around molecular clouds as well as in the stellar bar, where they could have been injected into the ISM during mass loss from old stars.

We combine new Spitzer Space Telescope observations in the mid- and far-infrared with SCUBA 850 micron observations to improve the measurement of dust temperatures, masses and luminosities for 11 galaxies of the SCUBA Local Universe Galaxy Survey (SLUGS). By fitting dust models we measure typical dust masses of 10E7.9 M_sol and dust luminosities of ~ 10E10 L_sol, for galaxies with modest star formation rates. The data presented in this paper combined with previous observations show that cold dust is present in all types of spiral galaxies and is a major contributor to their total luminosity. Because of the lower dust temperature of the SCUBA sources measured in this paper, they have flatter Far-IR nu F_nu(160um)/nu F_nu(850um) slopes than the larger Spitzer Nearby Galaxies Survey (SINGS), the sample that provides the best measurements of the dust properties of galaxies in the nearby universe. The new data presented here added to SINGS extend the parameter space that is well covered by local galaxies, providing a comprehensive set of templates that can be used to interpret the observations of nearby and distant galaxies.

We present empirical relations describing excess emission from evolved stars in the Large Magellanic Cloud (LMC) using data from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey which includes the IRAC 3.6, 4.5, 5.8 and 8.0 \mu m and MIPS 24, 70 and 160 \mu m bands. We combine the SAGE data with the Two Micron All Sky Survey (2MASS; J, H and Ks) and the optical Magellanic Cloud Photometric Survey (MCPS; U, B, V and I) point source catalogs to create complete spectral energy distributions (SEDs) of the asymptotic giant branch (AGB) star candidates in the LMC. AGB star outflows are among the main producers of dust in a galaxy, and this mass loss results in an excess in the fluxes observed in the 8 and 24 \mic m bands. We identify oxygen-rich, carbon-rich and extreme AGB star populations in our sample based on their 2MASS and IRAC colors. We calculate excesses fluxes in the mid-IR bands by comparison of the SEDs with model photospheres. We find about 16,000 O-rich, 6300 C-rich and 1000 extreme sources with reliable 8 \mu m excesses, and about 4500 O-rich, 5300 C-rich and 960 extreme sources with reliable 24 \mic m excesses. The excesses are in the range 0.1 mJy--5 Jy. The 8 and 24 \mic m excesses for all three types of AGB candidates show a general increasing trend with luminosity. The dust color temperature derived from the ratio of the 8 and 24 \mic m excesses decreases with an increase in excess, while the 24 \mic m optical depth increases with excess. The extreme AGB candidates are the major contributors to the mass loss, and we estimate the total AGB mass-loss return to the LMC to be (5.9--13)\times 10^{-3} Msun yr^{-1}.

We present a mid-infrared spectroscopic data cube of the central part of 30 Doradus, observed with Spitzer's IRS and MIPS/SED mode. Aromatic dust emission features and emission lines from molecular and atomic hydrogen are detected but not particularly strong. The dominant spectral features are emission lines from moderately ionized species of argon, neon, and sulphur, which are used to determine the physical conditions in the ionized gas. The ionized gas excitation shows strong variations on parsec scales, some of which can plausibly be associated with individual hot stars. We fit the ionic line strengths with photoionization and shock models, and find that photoionization dominates in the region. The ionization parameter U traces the rim of the central bubble, as well as highlighting isolated sources of ionization, and at least one quiescent clump. The hardness of the ionizing radiation field T_rad reveals several "hot spots" that are either the result of individual very hot stars or trace the propagation of the diffuse ionizing field through the surrounding neutral cloud. Consistent with other measurements of giant molecular hydrogen regions, log(U) ranges between -3 and -0.75, and T_rad between 30000 and 85000K.

The HII complex N44 in the Large Magellanic Cloud (LMC) provides an excellent site to perform a detailed study of star formation in a mild starburst, as it hosts three regions of star formation at different evolutionary stages and it is not as complicated and confusing as the 30 Doradus giant HII region. We have obtained Spitzer Space Telescope observations and complementary ground-based 4m uBVIJK observations of N44 to identify candidate massive young stellar objects (YSOs). We further classify the YSOs into Types I, II, and III, according to their spectral energy distributions (SEDs). In our sample of 60 YSO candidates, ~65% of them are resolved into multiple components or extended sources in high-resolution ground-based images. We have modeled the SEDs of 36 YSOs that appear single or dominant within a group. We find good fits for Types I and I/II YSOs,but Types II and II/III YSOs show deviations between their observed SEDs and models that do not include PAH emission. We have also found that some Type III YSOs have central holes in their disk components. YSO counterparts are found in four ultracompact HII regions and their stellar masses determined from SED model fits agree well with those estimated from the ionization requirements of the HII regions. The distribution of YSOs is compared with those of the underlying stellar population and interstellar gas conditions to illustrate a correlation between the current formation of O-type stars and previous formation of massive stars. Evidence of triggered star formation is also presented.

(abridged) We study photoelectric heating throughout the Large Magellanic Cloud. We quantify the importance of the [CII] cooling line and the photoelectric heating process of various environments in the LMC and investigate which parameters control the extent of photoelectric heating. We use the BICE [CII] map and the Spitzer/SAGE infrared maps. We examine the spatial variations in the efficiency of photoelectric heating: photoelectric heating rate over power absorbed by grains. We correlate the photoelectric heating efficiency and the emission from various dust constituents and study the variations as a function of H\alpha emission, dust temperatures, and the total infrared luminosity. From this we estimate radiation field, gas temperature, and electron density. We find systematic variations in photoelectric efficiency. The highest efficiencies are found in the diffuse medium, while the lowest coincide with bright star-forming regions (~1.4 times lower). The [CII] line emission constitutes 1.32% of the far infrared luminosity across the whole of the LMC. We find correlations between the [CII] emission and ratios of the mid infrared and far infrared bands, which comprise various dust constituents. The correlations are interpreted in light of the spatial variations of the dust abundance and by the local environmental conditions that affect the dust emission properties. As a function of the total infrared surface brightness, S_{TIR}, the [CII] surface brightness can be described as: S_{[CII]}=1.25 S_{TIR}^{0.69} [10^{-3} erg s^{-1} cm^{-2} sr^{-1}]. The [CII] emission is well-correlation with the 8 micrometer emission, suggesting that the polycyclic aromatic hydrocarbons play a dominant role in the photoelectric heating process.

High resolution mid-infrared spectra are presented for 155 nuclear and extranuclear regions from the Spitzer Infrared Nearby Galaxies Survey (SINGS). The fluxes for nine atomic forbidden and three molecular hydrogen mid-infrared emission lines are also provided, along with upper limits in key lines for infrared-faint targets. The SINGS sample shows a wide range in the ratio of [SIII]18.71um/[SIII]33.48um, but the average ratio of the ensemble indicates a typical interstellar electron density of 300-400 cm^{-3} on ~23"x15" scales and 500-600 cm^{-3} using ~11"x9" apertures, independent of whether the region probed is a star-forming nuclear, a star-forming extranuclear, or an AGN environment. Evidence is provided that variations in gas-phase metallicity play an important role in driving variations in radiation field hardness, as indicated by [NeIII]15.56um/[NeII]12.81um, for regions powered by star formation. Conversely, the radiation hardness for galaxy nuclei powered by accretion around a massive black hole is independent of metal abundance. Furthermore, for metal-rich environments AGN are distinguishable from star-forming regions by significantly larger [NeIII]15.56um/[NeII]12.81um ratios. Finally, [FeII]25.99um/[NeII]12.81um versus [SiII]34.82um/[SIII]33.48um also provides an empirical method for discerning AGN from normal star-forming sources. However, similar to [NeIII]15.56um/[NeII]12.81um, these mid-infrared line ratios lose their AGN/star-formation diagnostic powers for very low metallicity star-forming systems with hard radiation fields.

The Tail region of the Small Magellanic Cloud (SMC) was imaged using the MIPS instrument on the Spitzer Space Telescope as part of the SAGE-SMC Spitzer Legacy. Diffuse infrared emission from dust was detected in all the MIPS bands. The Tail gas-to-dust ratio was measured to be 1200 +/- 350 using the MIPS observations combined with existing IRAS and HI observations. This gas-to-dust ratio is higher than the expected 500-800 from the known Tail metallicity indicating possible destruction of dust grains. Two cluster regions in the Tail were resolved into multiple sources in the MIPS observations and local gas-to-dust ratios were measured to be ~440 and ~250 suggests dust formation and/or significant amounts of ionized gas in these regions. These results support the interpretation that the SMC Tail is a tidal tail recently stripped from the SMC that includes gas, dust, and young stars.

We present initial results and source lists of variable sources in the Large Magellanic Cloud (LMC) for which we detect thermal infrared variability from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey, which had 2 epochs of photometry separated by three months. The SAGE survey mapped a 7 degree by 7 degree region of the LMC using the IRAC and the MIPS instruments on board Spitzer. Variable sources are identified using a combination of the IRAC 3.6, 4.5, 5.8, 8.0 \micron bands and the MIPS 24 \micron bands. An error-weighted flux difference between the two epochs is used to assess the variability. Of the ~ 3 million sources detected at both epochs we find ~ 2,000 variable sources for which we provide electronic catalogs. Most of the variable sources can be classified as asymptotic giant branch (AGB) stars. A large fraction (> 66%) of the extreme AGB stars are variable and only smaller fractions of carbon-rich (6.1%) and oxygen-rich (2.0%) stars are detected as variable. We also detect a population of variable young stellar object candidates.

We explore the stellar population of M31 in a Spitzer Space Telescope survey utilizing IRAC and MIPS observations. Red supergiants are the brightest objects seen in the infrared; they are a prominent evolutionary phase. Due to their circumstellar envelopes, many of these radiate the bulk of their luminosity at IRAC wavelengths and do not stand out in the near infrared or optically. Going fainter, we see large numbers of luminous asymptotic giant branch (AGB) stars, many of which are known long period variables. Relative to M33, the AGB carbon star population of M31 appears sparse, but this needs to be spectroscopically confirmed.

We determine an absolute calibration for the MIPS 24 microns band and recommend adjustments to the published calibrations for 2MASS, IRAC, and IRAS photometry to put them on the same scale. We show that consistent results are obtained by basing the calibration on either an average A0V star spectral energy distribution (SED), or by using the absolutely calibrated SED of the sun in comparison with solar-type stellar photometry (the solar analog method). After the rejection of a small number of stars with anomalous SEDs (or bad measurements), upper limits of ~ 1.5% (rms) are placed on the intrinsic infrared SED variations in both A dwarf and solar-type stars. These types of stars are therefore suitable as general-purpose standard stars in the infrared. We provide absolutely calibrated SEDs for a standard zero magnitude A star and for the sun to allow extending this work to any other infrared photometric system. They allow the recommended calibration to be applied from 1 to 25 microns with an accuracy of ~2 %, and with even higher accuracy at specific wavelengths such as 2.2, 10.6, and 24 microns, near which there are direct measurements. However, we confirm earlier indications that Vega does not behave as a typical A0V star between the visible and the infrared, making it problematic as the defining star for photometric systems. The integration of measurements of the sun with those of solar-type stars also provides an accurate estimate of the solar SED from 1 through 30 microns, which we show agrees with theoretical models.

We present Spitzer observations of IC 443 obtained with MIPS and IRS as part of our GTO program on the astrophysics of ejecta from evolved stars. We find that the overall morphology at mid/far IR wavelengths resembles even more closely a loop or a shell than the ground based optical and/or near IR images.The dust temperature map, based on the 70/160micron ratio, shows a range from 18 to 30 K degrees. The IRS spectra confirm the findings from previous near+mid IR spectroscopic observations of a collisionally excited gas, atomic and molecular, rich in fine structure atomic and pure H2 rotational emission lines, respectively. The spectroscopic shock indicator, [Ne II] 12.8micron, suggests shock velocities ranging from 60-90 km/s, consistent with the values derived from other indicators.

The aromatic features in M101 were studied spectroscopically and photometrically using observations from all three instruments on the Spitzer Space Telescope. The global SED of M101 shows strong aromatic feature (commonly called PAH feature) emission. The spatially resolved spectral and photometric measurements of the aromatic feature emission show strong variations with significantly weaker emission at larger radii. We compare these variations with changes in the ionization index (as measured by [NeIII]/[NeII] and [SIV/SIII], which we probe over the ranges 0.03-20 and 0.044-15 respectively) and metallicity (expressed as log(O/H)+12, which ranges from 8.1 to 8.8). Over these ranges, the spectroscopic equivalent widths of the aromatic features from 7 HII regions and the nucleus were found to correlate better with ionization index than with metallicity. This implies that the weakening of the aromatic emission in massive star forming regions is due primarily to processing of the dust grains in these environments, not to differences in how they form (although formation could still be important on a secondary basis). The behavior of the correlation between the aromatic feature equivalent widths and ionization index can be described as a constant equivalent width until a threshold in ionization index is reached ([NeIII]/[NeII] ~ 1), above which the equivalent widths decrease with a power law dependence. This behavior for M101 HII regions is also seen for the sample of starburst galaxies presented in the companion study of Engelbracht et al. (2008) which expands the range of [NeIII]/[NeII] ratios to 0.03-25 and log(O/H)+12 values to 7.1-8.8. The form of the correlation explains seemingly contradictory results present in the literature. [abridged]

The Spectral Energy Distribution (SED) mode of the Multiband Imaging Photometer for Spitzer (MIPS) Space Telescope provides low-spectral resolution (R ~ 15-25) spectroscopy in the far infrared using the MIPS 70 um detector. A reflective grating provides a dispersion of 1.7 um per pixel, and an effective wavelength coverage of 52.8--98.7 um over detector rows 1-27. The final 5 detector rows are contaminated by second-order diffracted light and are left uncalibrated. The flux calibration is based on observations of MIPS calibration stars with 70 um flux densities of 0.5--15 Jy. The point-source flux calibration accuracy is estimated to be 10% or better down to about 0.5 Jy at the blue end of the spectrum and to 2 Jy near the red end. With additional uncertainties from the illumination and aperture corrections included, the surface brightness calibration of extended sources is accurate to ~15%. Repeatability of better than 5% is found for the SED mode through multiple measurements of several calibration stars.

We determine the variation in the nonthermal radio spectral index in the nearby spiral galaxy M33. We separate the thermal and nonthermal components of the radio continuum emission without the assumption of a constant nonthermal spectral index. Using the Spitzer FIR data at 70 and 160 $\mu$m and a standard dust model, we de-redden the H$\alpha$ emission. The extinction-corrected H$\alpha$ emission serves as a template for the thermal free-free radio emission. Subtracting from the observed 3.6 and 20cm emission (Effelsberg and the VLA) this free-free emission, we obtain maps of the nonthermal intensity and spectral index.

The origin of the tight radio--IR correlation in galaxies has not been fully understood. One reason is the uncertainty about which heating sources (stars or diffuse interstellar radiation field)provide the energy that is absorbed by dust and re-radiated in IR. Another problem is caused by comparing the IR emission with the thermal and nonthermal components of the radio continuum emission separated by simplistically assuming a constant nonthermal spectral index. We use the data at the $Spitzer$ MIPS wavelengths of 24, 70, and 160 $\mu$m, as well as recent radio continuum map at 3.6 cm observed with the 100--m Effelsberg telescope. Using the wavelet transformation, we separate diffuse emission components from compact sources and study the radio-IR correlation at various scales. We also investigate the IR correlations with the thermal and nonthermal radio emissions separated by our developed method. A H$\alpha$ map serves as a tracer of star forming regions.

We present infrared observations of 66 starburst galaxies over a wide range of oxygen abundances, to measure how metallicity affects their dust properties. The data include imaging and spectroscopy from the Spitzer Space Telescope, supplemented by groundbased near-infrared imaging. We confirm a strong correlation of aromatic emission with metallicity, with a threshold at a metallicity [12+log(O/H)]~8. The large scatter in both the metallicity and radiation hardness dependence of this behavior implies that it is not due to a single effect, but to some combination. We show that the far-infrared color temperature of the large dust grains increases towards lower metallicity, peaking at a metallicity of 8 before turning over. We compute dust masses and compare them to HI masses from the literature to derive the gas to dust ratio, which increases by nearly 3 orders of magnitude between solar metallicity and a metallicity of 8, below which it flattens out. The abrupt change in aromatic emission at mid-infrared wavelengths thus appears to be reflected in the far-infrared properties, indicating that metallicity changes affect the composition of the full range of dust grain sizes that dominate the infrared emission. In addition, we find that the ratio L(8 micron)/L(TIR), important for calibrating 24 micron measurements of high redshift galaxies, increases slightly as the metallicity decreases from ~solar to ~50% of solar, and then decreases by an order of magnitude with further decreases in metallicity. Although the great majority of galaxies show similar patterns of behavior as described above, there are three exceptions, SBS 0335-052E, Haro 11, and SHOC 391. Their infrared SEDs are dominated energetically by the mid-IR near 24 micron rather than by the 60 - 200 micron region. (Abridged)

We present IRAC and MIPS images and photometry of a sample of previously known planetary nebulae (PNe) from the SAGE survey of the Large Magellanic Cloud (LMC) performed with the Spitzer Space Telescope. Of the 233 known PNe in the survey field, 185 objects were detected in at least two of the IRAC bands, and 161 detected in the MIPS 24 micron images. Color-color and color-magnitude diagrams are presented using several combinations of IRAC, MIPS, and 2MASS magnitudes. The location of an individual PN in the color-color diagrams is seen to depend on the relative contributions of the spectral components which include molecular hydrogen, polycyclic aromatic hydrocarbons (PAHs), infrared forbidden line emission from the ionized gas, warm dust continuum, and emission directly from the central star. The sample of LMC PNe is compared to a number of Galactic PNe and found to not significantly differ in their position in color-color space. We also explore the potential value of IR PNe luminosity functions (LFs) in the LMC. IRAC LFs appear to follow the same functional form as the well-established [O III] LFs although there are several PNe with observed IR magnitudes brighter than the cut-offs in these LFs.

We determine the variation in the nonthermal radio spectral index in the nearby spiral galaxy M33 at a linear resolution of 360 pc. We separate the thermal and nonthermal components of the radio continuum emission without the assumption of a constant nonthermal spectral index. Using the Spitzer FIR data at 70 and 160 $\mu$m and a standard dust model, we deredden the H$\alpha$ emission. The extinction corrected H$\alpha$ emission serves as a template for the thermal free-free radio emission. Subtracting from the observed 3.6 cm and 20 cm emission (Effelsberg and the VLA) this free-free emission, we obtain the nonthermal maps. A constant electron temperature used to obtain the thermal radio intensity seems appropriate for M~33 which, unlike the Milky Way, has a shallow metallicity gradient. For the first time, we derive the distribution of the nonthermal spectral index across a galaxy, M33. We detect strong nonthermal emission from the spiral arms and star-forming regions. Wavelet analysis shows that at 3.6 cm the nonthermal emission is dominated by contributions from star-forming regions, while it is smoothly distributed at 20 cm. For the whole galaxy, we obtain thermal fractions of 51% and 18% at 3.6 cm and 20 cm, respectively. The thermal emission is slightly stronger in the southern than in the northern half of the galaxy. We find a clear radial gradient of mean extinction in the galactic plane. The nonthermal spectral index map indicates that the relativistic electrons suffer energy-loss when diffusing from their origin in star-forming regions towards interarm regions and the outer parts of the galaxy. We also conclude that the radio emission is mostly nonthermal at R $>$ 5 kpc in M33.

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.

Dust in galaxies can be mapped by either the FIR/sub-mm emission, the optical or infrared reddening of starlight, or the extinction of a known background source. We compare two dust extinction measurements for a set of fifteen sections in thirteen nearby galaxies, to determine the scale of the dusty ISM responsible for disk opacity: one using stellar reddening and the other a known background source. In our earlier papers, we presented extinction measurements of 29 galaxies, based on calibrated counts of distant background objects identified though foreground disks in HST/WFPC2 images. For the 13 galaxies that overlap with the Spitzer Infrared Nearby Galaxies Survey (SINGS), we now compare these results with those obtained from an I-L color map. Our goal is to determine whether or not a detected distant galaxy indicates a gap in the dusty ISM, and hence to better understand the nature and geometry of the disk extinction. We find that distant galaxies are predominantly in low-extinction sections marked by the color maps, indicating that their number depends both on the cloud cover of {\it Spitzer}-resolved dust structures --mostly the spiral arms--and a diffuse, unresolved underlying disk. We note that our infrared color map (E[I-L]) underestimates the overall dust presence in these disks severely, because it implicitly assumes the presence of a dust screen in front of the stellar distribution.

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.

We describe the absolute calibration of the Multiband Imaging Photometer for Spitzer (MIPS) 160 micron channel. After the on-orbit discovery of a near-IR ghost image that dominates the signal for sources hotter than about 2000 K, we adopted a strategy utilizing asteroids to transfer the absolute calibrations of the MIPS 24 and 70 micron channels to the 160 micron channel. Near-simultaneous observations at all three wavelengths are taken, and photometry at the two shorter wavelengths is fit using the Standard Thermal Model. The 160 micron flux density is predicted from those fits and compared with the observed 160 micron signal to derive the conversion from instrumental units to surface brightness. The calibration factor we derive is 41.7 MJy/sr/MIPS160 (MIPS160 being the instrumental units). The scatter in the individual measurements of the calibration factor, as well as an assesment of the external uncertainties inherent in the calibration, lead us to adopt an uncertainty of 5.0 MJy/sr/MIPS160 (12%) for the absolute uncertainty on the 160 micron flux density of a particular source as determined from a single measurement. For sources brighter than about 2 Jy, non-linearity in the response of the 160 micron detectors produces an under-estimate of the flux density: for objects as bright as 4 Jy, measured flux densities are likely to be ~20% too low. This calibration has been checked against that of ISO (using ULIRGS) and IRAS (using IRAS-derived diameters), and is consistent with those at the 5% level.

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.

Using the new capabilities of the Spitzer Space Telescope and extensive multiwavelength data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), it is now possible to study the infrared properties of star formation in nearby galaxies down to scales equivalent to large HII regions. We are therefore able to determine what fraction of large, infrared-selected star-forming regions in normal galaxies are highly obscured and address how much of the star formation we miss by relying solely on the optical portion of the spectrum. Employing a new empirical method for deriving attenuations of infrared-selected star-forming regions we investigate the statistics of obscured star formation on 500pc scales in a sample of 38 nearby galaxies. We find that the median attenuation is 1.4 magnitudes in H-alpha and that there is no evidence for a substantial sub-population of uniformly highly-obscured star-forming regions. The regions in the highly-obscured tail of the attenuation distribution (A_H-alpha > 3) make up only ~4% of the sample of nearly 1800 regions, though very embedded infrared sources on the much smaller scales and lower luminosities of compact and ultracompact HII regions are almost certainly present in greater numbers. The highly-obscured cases in our sample are generally the bright, central regions of galaxies with high overall attenuation but are not otherwise remarkable. We also find that a majority of the galaxies show decreasing radial trends in H-alpha attenuation. The small fraction of highly-obscured regions seen in this sample of normal, star-forming galaxies suggests that on 500pc scales the timescale for significant dispersal or break up of nearby, optically-thick dust clouds is short relative to the lifetime of a typical star-forming region.

We present Spitzer observations of the dark globule CB190 (L771). We observe a roughly circular 24 micron shadow with a 70 arcsec radius. The extinction profile of this shadow matches the profile derived from 2MASS photometry at the outer edges of the globule and reaches a maximum of ~32 visual magnitudes at the center. The corresponding mass of CB190 is ~10 Msun. Our 12CO and 13CO J = 2-1 data over a 10 arcmin X 10 arcmin region centered on the shadow show a temperature ~10 K. The thermal continuum indicates a similar temperature for the dust. The molecular data also show evidence of freezeout onto dust grains. We estimate a distance to CB190 of 400 pc using the spectroscopic parallax of a star associated with the globule. Bonnor-Ebert fits to the density profile, in conjunction with this distance, yield xi_max = 7.2, indicating that CB190 may be unstable. The high temperature (56 K) of the best fit Bonnor-Ebert model is in contradiction with the CO and thermal continuum data, leading to the conclusion that the thermal pressure is not enough to prevent free-fall collapse. We also find that the turbulence in the cloud is inadequate to support it. However, the cloud may be supported by the magnetic field, if this field is at the average level for dark globules. Since the magnetic field will eventually leak out through ambipolar diffusion, it is likely that CB190 is collapsing or in a late pre-collapse stage.

The absolute calibration and characterization of the Multiband Imaging Photometer for Spitzer (MIPS) 70 micron coarse- and fine-scale imaging modes are presented based on over 2.5 years of observations. Accurate photometry (especially for faint sources) requires two simple processing steps beyond the standard data reduction to remove long-term detector transients. Point spread function (PSF) fitting photometry is found to give more accurate flux densities than aperture photometry. Based on the PSF fitting photometry, the calibration factor shows no strong trend with flux density, background, spectral type, exposure time, or time since anneals. The coarse-scale calibration sample includes observations of stars with flux densities from 22 mJy to 17 Jy, on backgrounds from 4 to 26 MJy sr^-1, and with spectral types from B to M. The coarse-scale calibration is 702 +/- 35 MJy sr^-1 MIPS70^-1 (5% uncertainty) and is based on measurements of 66 stars. The instrumental units of the MIPS 70 micron coarse- and fine-scale imaging modes are called MIPS70 and MIPS70F, respectively. The photometric repeatability is calculated to be 4.5% from two stars measured during every MIPS campaign and includes variations on all time scales probed. The preliminary fine-scale calibration factor is 2894 +/- 294 MJy sr^-1 MIPS70F^-1 (10% uncertainty) based on 10 stars. The uncertainty in the coarse- and fine-scale calibration factors are dominated by the 4.5% photometric repeatability and the small sample size, respectively. The 5-sigma, 500 s sensitivity of the coarse-scale observations is 6-8 mJy. This work shows that the MIPS 70 micron array produces accurate, well calibrated photometry and validates the MIPS 70 micron operating strategy, especially the use of frequent stimulator flashes to track the changing responsivities of the Ge:Ga detectors.

We present the stellar calibrator sample and the conversion from instrumental to physical units for the 24 micron channel of the Multiband Imaging Photometer for Spitzer (MIPS). The primary calibrators are A stars, and the calibration factor based on those stars is 4.54*10^{-2} MJy sr^{-1} (DN/s)^{-1}, with a nominal uncertainty of 2%. We discuss the data-reduction procedures required to attain this accuracy; without these procdures, the calibration factor obtained using the automated pipeline at the Spitzer Science Center is 1.6% +/- 0.6% lower. We extend this work to predict 24 micron flux densities for a sample of 238 stars which covers a larger range of flux densities and spectral types. We present a total of 348 measurements of 141 stars at 24 micron. This sample covers a factor of ~460 in 24 micron flux density, from 8.6 mJy up to 4.0 Jy. We show that the calibration is linear over that range with respect to target flux and background level. The calibration is based on observations made using 3-second exposures; a preliminary analysis shows that the calibration factor may be 1% and 2% lower for 10- and 30-second exposures, respectively. We also demonstrate that the calibration is very stable: over the course of the mission, repeated measurements of our routine calibrator, HD 159330, show a root-mean-square scatter of only 0.4%. Finally, we show that the point spread function (PSF) is well measured and allows us to calibrate extended sources accurately; Infrared Astronomy Satellite (IRAS) and MIPS measurements of a sample of nearby galaxies are identical within the uncertainties.

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.

We investigate the energy sources of the infrared (IR) emission and their relation to the radio continuum emission at various spatial scales within the Scd galaxy M33. We use the wavelet transform to analyze IR data at the Spitzer wavelengths of 24, 70, and 160$\mu$m, as well as recent radio continuum data at 3.6cm and 20cm. An H$\alpha$ map serves as a tracer of the star forming regions and as an indicator of the thermal radio emission. We find that the dominant scale of the 70$\mu$m emission is larger than that of the 24$\mu$m emission, while the 160$\mu$m emission shows a smooth wavelet spectrum. The radio and H$\alpha$ maps are well correlated with all 3 MIPS maps, although their correlations with the 160$\mu$m map are weaker. After subtracting the bright HII regions, the 24 and 70$\mu$m maps show weaker correlations with the 20cm map than with the 3.6cm map at most scales. We also find a strong correlation between the 3.6cm and H$\alpha$ emission at all scales. Comparing the results with and without the bright HII regions, we conclude that the IR emission is influenced by young, massive stars increasingly with decreasing wavelength from 160 to 24$\mu$m. The radio-IR correlations indicate that the warm dust-thermal radio correlation is stronger than the cold dust-nonthermal radio correlation at scales smaller than 4kpc. A perfect 3.6cm-H$\alpha$ correlation implies that extinction has no significant effect on H$\alpha$ emitting structures.

Mid-infrared observations of the Andromeda galaxy, M31, obtained with the Infrared Array Camera on board the Spitzer Space Telescope, are presented. The image mosaics cover areas of approximate 3.7deg x 1.6deg and include the satellite galaxies M32 and NGC 205. The appearance of M31 varies dramatically in the different mid-infrared bands, from the smooth bulge and disk of the old stellar population seen at 3.6um to the well-known '10 kpc ring' dominating the 8um image. The similarity of the 3.6um and optical isophotes and nearly constant optical-mid-infrared color over the inner 400 arcsec confirms that there is no significant extinction at optical wavelengths in M31's bulge. The nuclear colors indicate the presence of dust but not an infrared-bright nucleus. The integrated 8um non-stellar luminosity implies a star formation rate of 0.4 Msun/yr, consistent with other indicators that show M31 to be a quiescent galaxy.

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.

The ultraviolet-to-radio continuum spectral energy distributions are presented for all 75 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). A principal component analysis of the sample shows that most of the sample's spectral variations stem from two underlying components, one representative of a galaxy with a low infrared-to-ultraviolet ratio and one representative of a galaxy with a high infrared-to-ultraviolet ratio. The influence of several parameters on the infrared-to-ultraviolet ratio is studied (e.g., optical morphology, disk inclination, far-infrared color, ultraviolet spectral slope, and star formation history). Consistent with our understanding of normal star-forming galaxies, the SINGS sample of galaxies in comparison to more actively star-forming galaxies exhibits a larger dispersion in the infrared-to-ultraviolet versus ultraviolet spectral slope correlation. Early type galaxies, exhibiting low star formation rates and high optical surface brightnesses, have the most discrepant infrared-to-ultraviolet correlation. These results suggest that the star formation history may be the dominant regulator of the broadband spectral variations between galaxies. Finally, a new discovery shows that the 24 micron morphology can be a useful tool for parametrizing the global dust temperature and ultraviolet extinction in nearby galaxies. The dust emission in dwarf/irregular galaxies is clumpy and warm accompanied by low ultraviolet extinction, while in spiral galaxies there is typically a much larger diffuse component of cooler dust and average ultraviolet extinction. For galaxies with nuclear 24 micron emission, the dust temperature and ultraviolet extinction are relatively high compared to disk galaxies.

We discuss images of M87 from 3.6 to 160 um obtained with Spitzer. As found previously, there is an excess in the far infrared over a simple power law interpolation from the radio to the resolved nonthermal features in the mid-infrared and optical. We show that this excess is most likely warm dust in the galaxy itself, and that the properties of this emission component are similar to the far infrared emission of normal giant elliptical galaxies. The new observations let us determine the spectrum of the jet and surrounding lobes of nonthermal emission. We find that even in the lobes the synchrotron break frequency is in the optical, probably requiring in situ particle acceleration not only in the jet but in the lobes as well.

We describe Spitzer images of a sample of dwarf and low surface brightness galaxies, using the high sensitivity and spatial resolution to explore the morphologies of dust in these galaxies. For the starbursting dwarf UGC 10445, we present a complete infrared spectral energy distribution and modeling of its individual dust components. We find that its diffuse cold (T~19K) dust component extends beyond its near-infrared disk and speculate that the most plausible source of heating is ultraviolet photons from starforming complexes. We find that the mass of T~19K dust in UGC 10445 is surprisingly large, with a lower limit of 3 x 10^6 M_solar. We explore the implications of having such a high dust content on the nature and evolution of the galaxy.

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.

Color-magnitude diagrams (CMDs) are presented for the Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) survey of the Large Magellanic Cloud (LMC). IRAC and MIPS 24 um epoch one data are presented. These data represent the deepest, widest mid-infrared CMDs of their kind ever produced in the LMC. Combined with the 2MASS survey, the diagrams are used to delineate the evolved stellar populations in the Large Magellanic Cloud as well as Galactic foreground and extragalactic background populations. Some 32000 evolved stars brighter than the tip of the red giant branch are identified. Of these, approximately 17500 are classified as oxygen-rich, 7000 carbon-rich, and another 1200 as ``extreme'' asymptotic giant branch (AGB) stars. Brighter members of the latter group have been called ``obscured'' AGB stars in the literature owing to their dusty circumstellar envelopes. A large number (1200) of luminous oxygen--rich AGB stars/M supergiants are also identified. Finally, there is strong evidence from the 24 um MIPS channel that previously unexplored, lower luminosity oxygen-rich AGB stars contribute significantly to the mass loss budget of the LMC (1200 such sources are identified).

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.

We present Spitzer Space Telescope images of the isolated dwarf galaxy UGC 10445. The galaxy is detected at all photometric bands (3.6-160um) as well as in the Multiband Imaging Photometer for Spitzer (MIPS) spectral energy distribution mode (55-95um). We derive a star formation rate of 0.25 M_sun/yr based on H-alpha and infrared flux densities. There is over 10^6 solar masses of cold dust (T~18K) in the galaxy, represented by 160um emission, that extends to a larger radius than the ultraviolet (UV), optical and near-infrared light. Such extended emission has been seen previously only in dwarf galaxies in cluster environments. We suggest the source of heating for this dust is UV light originating in star forming complexes. To produce the large quantity of dust requires a higher rate of star formation in the past than is observed currently.

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

We present Spitzer Space Telescope IRAC and MIPS observations of the galactic globular cluster M15 (NGC 7078), one of the most metal-poor clusters with a [Fe/H] = -2.4. Our Spitzer images reveal a population of dusty red giants near the cluster center, a previously detected planetary nebula (PN) designated K648, and a possible detection of the intra-cluster medium (ICM) arising from mass loss episodes from the evolved stellar population. Our analysis suggests 9 (+/-2) x 10^-4 solar masses of dust is present in the core of M15, and this material has accumulated over a period of approximately 10^6 years, a timescale ten times shorter than the last galactic plane crossing event. We also present Spitzer IRS follow up observations of K648, including the detection of the [NeII] 12.81 micron line, and discuss abundances derived from infrared fine structure lines.

We present late-time optical and mid-infrared observations of the Type-II supernova 2003gd in NGC 628. Mid-infrared excesses consistent with cooling dust in the ejecta are observed 499-678 days after outburst, and are accompanied by increasing optical extinction and growing asymmetries in the emission-line profiles. Radiative-transfer models show that up to 0.02 solar masses of dust has formed within the ejecta, beginning as early as 250 days after outburst. These observations show that dust formation in supernova ejecta can be efficient and that massive-star supernovae can be major dust producers throughout the history of the Universe.

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.

The behavior of the aromatic features as a function of metallicity/radiation field hardness in M101 HII regions and starburst galaxies is investigated using Spitzer/IRS spectra. These spectra reveal that the 6.2, 7.8+8.6, and 11.3 micron aromatics have different dependencies on metallicity/radiation field hardness. Specifically, the 6.2 and 7.8+8.6 micron aromatics are weak or absent at a metallicity, 12+log(O/H), of approx. 8.0 while the 11.3 micron feature remains relatively strong. These results apply to both the M101 HII regions and starbursts showing that HII regions can indeed be thought of as ``mini-starbursts.'' Comparison of this work with various candidate materials results in a good match with the annealing behavior in Quenched Carbonaceous Composite (QCC) which is an amorphous, hydrogenated, carbonaceous solid.

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.

We present observations at 3.6, 4.5, 5.8, 8, 24, 70 & 160um of NGC 205, the dwarf elliptical companion of M31, obtained with the Spitzer Space Telescope. The point sources subtracted images at 8 and 24um display a complex and fragmented infrared emission coming from both very small dust particles and larger grains. The extended dust emission is spatially concentrated in three main emission regions, seen at all wavelengths from 8 to 160um. These regions lie approximately along NGC 205's semi-major axis and range from ~100 to 300 pc in size. Based on our mid-to-far infrared flux density measurements alone, we derive a total dust mass estimate of the order of 3.2e4 solar masses, mainly at a temperature of ~20K. The gas mass associated with this component matches the predicted mass returned by the dying stars from the last burst of star formation in NGC 205 (~0.5 Gyr ago). Analysis of the Spitzer data combined with previous 1.1mm observations over a small central region or "Core'' (18" diameter), suggest the presence of very cold (~12K) dust and a dust mass about sixteen times higher than is estimated from the Spitzer measurements alone. Assuming a gas to dust mass ratio of 100, these two datasets, i.e.with and without the millimeter observations, suggest a total gas mass range of 3.2e6 to 5e7 solar masses.

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.

New images of M31 at 24, 70, and 160 micron taken with the Multiband Imaging Photometer for Spitzer (MIPS) reveal the morphology of the dust in this galaxy. This morphology is well represented by a composite of two logarithmic spiral arms and a circular ring (radius ~10 kpc) of star formation offset from the nucleus. The two spiral arms appear to start at the ends of a bar in the nuclear region and extend beyond the star forming ring. As has been found in previous work, the spiral arms are not continuous but composed of spiral segments. The star forming ring is very circular except for a region near M32 where it splits. The lack of well defined spiral arms and the prominence of the nearly circular ring suggests that M31 has been distorted by interactions with its satellite galaxies. Using new dynamical simulations of M31 interacting with M32 and NGC 205 we find that, qualitatively, such interactions can produce an offset, split ring like that seen in the MIPS images.

Variable A in M33 is a member of a rare class of highly luminous, evolved stars near the upper luminosity boundary that show sudden and dramatic shifts in apparent temperature due to the formation of optically thick winds in high mass loss episodes. Recent optical and infrared spectroscopy and imaging reveal that its ``eruption'' begun in $\sim$1950 has ended, {\it lasting $\approx$ 45 yrs}. Our current observations show major changes in its wind from a cool, dense envelope to a much warmer state surrounded by low density gas with rare emission lines of Ca II, [Ca II] and K I. Its spectral energy distribution has unexpectedly changed, especially at the long wavelengths, with a significant decrease in its apparent flux, while the star remains optically obscured. We conclude that much of its radiation is now escaping out of our line of sight. We attribute this to the changing structure and distribution of its circumstellar ejecta corresponding to the altered state of its wind as the star recovers from a high mass loss event.

(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.

Using Spitzer Space Telescope and Hubble Space Telescope observations of the superbubble N51D, we have identified three young stellar objects (YSOs) in dust globules, and made the first detection of a Herbig-Haro object outside the Galaxy. The spectral energy distributions of these YSOs suggest young massive stars with disk, envelope, and outflow cavities. The interstellar conditions are used to assess whether the star formation was spontaneous or induced by external pressure.

We report the discovery of a debris system associated with the $\sim 30$ Myr old G3/5V star HD 12039 using {\it Spitzer Space Telescope} observations from 3.6 -- 160$\mu$m. An observed infrared excess (L$_{\rm IR}$/L$_{\ast} = 1\times10^{-4}$) above the expected photosphere for $\lambda \gtrsim 14\mu$m is fit by thermally emitting material with a color temperature of T$\sim 110$ K, warmer than the majority of debris disks identified to date around Sun-like stars. The object is not detected at 70$\mu$m with a 3$\sigma$ upper limit 6 times the expected photospheric flux. The spectrum of the infrared excess can be explained by warm, optically thin material comprised of blackbody-like grains of size $\gtrsim 7 \mu$m that reside in a belt orbiting the star at 4--6 AU. An alternate model dominated by smaller grains, near the blow-out size $a\sim 0.5\mu$m, located at 30-40AU is also possible, but requires the dust to have been produced recently since such small grains will be expelled from the system by radiation pressure in $\sim$ few $\times 10^{2}$yrs.

We present ultraviolet, optical, near-infrared, Spitzer mid-infrared, and radio images of 14 radio-selected objects in M 33. These objects are thought to represent the youngest phase of star cluster formation. We have detected the majority of cluster candidates in M 33 at all wavelengths. From the near-IR images, we derived ages 2-10 Myr, K_S-band extinctions (A_K_S) of 0-1 mag, and stellar masses of 10^3-10^4 M_solar. We have generated spectral energy distributions (SEDs) of each cluster from 0.1 micron to 160 microns. From these SEDs, we have modeled the dust emission around these star clusters to determine the dust masses (1-10^3 M_solar) and temperatures (40-90 K) of the clusters' local interstellar medium. Extinctions derived from the JHK_S, Halpha, and UV images are similar to within a factor of 2 or 3. These results suggest that eleven of the fourteen radio-selected objects are optically-visible young star clusters with a surrounding H II region, that two are background objects, possibly AGN, and that one is a Wolf-Rayet star with a surrounding H II region.

We have identified five new R Coronae Borealis (RCB) stars in the Galactic bulge using the MACHO Project photometry database, raising the total number of known Galactic RCB stars to about 40. We have obtained spectra to confirm the identifications. The fact that four out of the five newly identified RCB stars are ``cool'' (T(eff) < 6000 K) rather than ``warm'' (T(eff) > 6000 K) suggests that the preponderance of warm RCB stars among the existing sample is a selection bias. These cool RCB stars are redder and fainter than their warm counterparts and may have been missed in surveys done with blue plates. Based on the number of new RCB stars discovered in the MACHO bulge fields, there may be ~250 RCB stars in the reddened "exclusion" zone toward the bulge.

We present Spitzer 70um and 160um observations of the Spitzer extragalactic First Look Survey (xFLS). The data reduction techniques and the methods for producing co-added mosaics and source catalogs are discussed. Currently, 26% of the 70um sample and 49% of the 160um-selected sources have redshifts. The majority of sources with redshifts are star-forming galaxies at z<0.5, while about 5% have infrared colors consistent with AGN. The observed infrared colors agree with the spectral energy distribution (SEDs) of local galaxies previously determined from IRAS and ISO data. The average 160um/70um color temperature for the dust is Td~= 30+/-5 K, and the average 70um/24um spectral index is alpha~= 2.4+/-0.4. The observed infrared to radio correlation varies with redshift as expected out to z~1 based on the SEDs of local galaxies. The xFLS number counts at 70um and 160um are consistent within uncertainties with the models of galaxy evolution, but there are indications that the current models may require slight modifications. Deeper 70um observations are needed to constrain the models, and redshifts for the faint sources are required to measure the evolution of the infrared luminosity function.

The carrier of the dust-associated photoluminescence process causing the extended red emission (ERE) in many dusty interstellar environments remains unidentified. Several competing models are more or less able to match the observed broad, unstructured ERE band. We now constrain the character of the ERE carrier further by determining the wavelengths of the radiation that initiates the ERE. Using the imaging capabilities of the Hubble Space Telescope, we have resolved the width of narrow ERE filaments appearing on the surfaces of externally illuminated molecular clouds in the bright reflection nebula NGC 7023 and compared them with the depth of penetration of radiation of known wavelengths into the same cloud surfaces. We identify photons with wavelengths shortward of 118 nm as the source of ERE initiation, not to be confused with ERE excitation, however. There are strong indications from the well-studied ERE in the Red Rectangle nebula and in the high-|b| Galactic cirrus that the photon flux with wavelengths shortward of 118 nm is too small to actually excite the observed ERE, even with 100% quantum efficiency. We conclude, therefore, that ERE excitation results from a two-step process. While none of the previously proposed ERE models can match these new constraints, we note that under interstellar conditions most polycyclic aromatic hydrocarbon (PAH) molecules are ionized to the di-cation stage by photons with E > 10.5 eV and that the electronic energy level structure of PAH di-cations is consistent with fluorescence in the wavelength band of the ERE. Therefore, PAH di-cations deserve further study as potential carriers of the ERE. (abridged)

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.

We investigate which conditions of dust attenuation and stellar populations allow models of dusty, continuously star-forming, bulge-less disk galaxies at 0.8<z<3.2 to meet the different colour selection criteria of high-z ``red'' galaxies (e.g. Rc-K>5.3, Ic-K>4, J-K>2.3). As a main novelty, we use stellar population models that include the thermally pulsating Asymptotic Giant Branch (TP-AGB) phase of stellar evolution. The star formation rate of the models declines exponentially as a function of time, the e-folding time being longer than 3 Gyr. In addition, we use calculations of radiative transfer of the stellar and scattered radiation through different dusty interstellar media in order to explore the wide parameter space of dust attenuation. We find that synthetic disks can exhibit red optical/near-infrared colours because of reddening by dust, but only if they have been forming stars for at least about 1 Gyr. Extremely few models barely exhibit Rc-K>5.3, if the inclination i=90 deg and if the opacity 2*tauV>6. Hence, Rc-K-selected galaxies at 1<z<2 most probably are either systems with an old, passively evolving bulge or starbursts. Synthetic disks at 1<z<2 exhibit 4<Ic-K<4.8, if they are seen edge on (i.e. at i about 90 deg) and if 2*tauV>0.5. This explains the large fraction of observed, edge-on disk-like galaxies with Ks<19.5 and F814W-Ks>4. Finally, models with 2<z<3.2 exhibit 2.3<J-K<3, with no bias towards i about 90 deg and for a large range in opacity (e.g. 2*tauV>1 for i about 70 deg). In conclusion, red disk-like galaxies at 0.8<z<3.2 may not necessarily be dustier than nearby disk galaxies (with 0.5<2*tauV<2) and/or much older than about 1 Gyr. This result is due both to a realistic description of dust attenuation and to the emission contribution by TP-AGB stars... (Abridged)

Blue Luminescence (BL) was first discovered in a proto-planetary nebula, the Red Rectangle (RR) surrounding the post-AGB star HD 44179. BL has been attributed to fluorescence by small, 3-4 ringed neutral polycyclic aromatic hydrocarbon (PAH) molecules, and was thought to be unique to the RR environment where such small molecules are actively being produced and shielded from the harsh interstellar radiation by a dense circumstellar disk. In this paper we present the BL spectrum detected in several ordinary reflection nebulae illuminated by stars having temperatures between 10,000 -- 23,000 K. All these nebulae are known to also exhibit the infrared emission features called aromatic emission features (AEFs) attributed to large PAHs. We present the spatial distribution of the BL in these nebulae. In the case of Ced~112, the BL is spatially correlated with mid-IR emission structures attributed to AEFs. These observations provide evidence for grain processing and possibly for in-situ formation of small grains and large molecules from larger aggregates. Most importantly, the detection of BL in these ordinary reflection nebulae suggests that the BL carrier is an ubiquitous component of the ISM and is not restricted to the particular environment of the RR.

Dust emission from the Type II supernova SN 2002hh in NGC 6946 has been detected at mid-infrared wavelengths by the Spitzer Space Telescope from 590 to 758 days after outburst and confirmed by higher angular resolution Gemini-N mid-IR observations. The day-600 5.8-24 micron emission can be fit by a 290-K blackbody having a luminosity of 1.6E+7 L_Sol. The minimum emitting radius of 1.1E+17 cm is too large for the emitting dust to have been formed in the supernova ejecta. Using radiative transfer models and realistic dust grain parameters, fits to the observed flux distribution could be obtained with an optically thick dust shell having a mass of 0.10-0.15 M_Sol, corresponding to a total dust+gas mass in excess of 10 M_Sol, suggesting a massive M supergiant or luminous blue variable precursor to this self-obscured object.

We analyze a sample of ~2600 MIPS/Spitzer 24mic sources brighter than ~80muJy and located in the Chandra Deep Field South to characterize the evolution of the comoving infrared (IR) energy density of the Universe up to z~1. Using published ancillary optical data we first obtain a nearly complete redshift determination for the 24mic objects associated with R<24 counterparts at z<1. We find that the 24mic population at 0.5<z<1 is dominated by ``Luminous Infrared Galaxies'' (i.e., 10^11 L_sol < L_IR < 10^12 L_sol), the counterparts of which appear to be also luminous at optical wavelengths and tend to be more massive than the majority of optically-selected galaxies. We finally derive 15mic and total IR luminosity functions (LFs) up to z~1. In agreement with the previous results from ISO and SCUBA and as expected from the MIPS source number counts, we find very strong evolution of the contribution of the IR-selected population with lookback time. Pure evolution in density is firmly excluded by the data, but we find considerable degeneracy between strict evolution in luminosity and a combination of increases in both density and luminosity (L*_IR prop. to (1+z)^{3.2_{-0.2}^{+0.7}}, Phi*_IR prop. to (1+z)^{0.7_{-0.6}^{+0.2}}). Our results imply that the comoving IR energy density of the Universe evolves as (1+z)^(3.9+/-0.4) up to z~1 and that galaxies luminous in the infrared (i.e., L_IR > 10^11 L_IR) are responsible for 70+/-15% of this energy density at z~1. Taking into account the contribution of the UV luminosity evolving as (1+z)^~2.5, we infer that these IR-luminous sources dominate the star-forming activity beyond z~0.7. The uncertainties affecting these conclusions are largely dominated by the errors in the k-corrections used to convert 24mic fluxes into luminosities.

We present the discovery of debris systems around three solar mass stars based upon observations performed with the Spitzer Space Telescope as part of a Legacy Science Program, ``the Formation and Evolution of Planetary Systems'' (FEPS). We also confirm the presence of debris around two other stars. All the stars exhibit infrared emission in excess of the expected photospheres in the 70 micron band, but are consistent with photospheric emission at <= 33 micron. This restricts the maximum temperature of debris in equilibrium with the stellar radiation to T < 70 K. We find that these sources are relatively old in the FEPS sample, in the age range 0.7 - 3 Gyr. Based on models of the spectral energy distributions, we suggest that these debris systems represent materials generated by collisions of planetesimal belts. We speculate on the nature of these systems through comparisons to our own Kuiper Belt, and on the likely planet(s) responsible for stirring the system and ultimately releasing dust through collisions. We further report observations of a nearby star HD 13974 (d =11 pc) that is indistinguishable from a bare photosphere at both 24 micron and 70 micron. The observations place strong upper limits on the presence of any cold dust in this nearby system (L_IR/L_* < 10^{-5.2}).

We examine colors from 3.6 micron to 24 micron as a function of metallicity (O/H) for a sample of 34 galaxies. The galaxies range over 2 orders of magnitude in metallicity. They display an abrupt shift in the 8 micron to 24 micron color between metallicities 1/3 to 1/5 of the solar value. The mean 8 micron to 24 micron flux density ratio below and above 12 + log (O/H) = 8.2 is 0.08 +/- 0.04 and 0.70 +/- 0.53, respectively. We use mid-infrared colors and spectroscopy to demonstrate that the shift is primarily due to a decrease in the 8 micron flux density as opposed to an increase in the 24 micron flux density. This result is most simply interpreted as due to a weakening at low metallicity of the mid-infrared emission bands usually attributed to PAHs (polycyclic aromatic hydrocarbons) relative to the small-grain dust emission. However, existing empirical spectral energy distribution models cannot account for the observed short-wavelength (i.e., below 8 micron) colors of the low-metallicity galaxies merely by reducing the strength of the PAH features; some other emission source (e.g., hot dust) is required.

Two images of Cassiopeia A obtained at 24 micrometer with the Spitzer Space Telescope over a one year time interval show moving structures outside the shell of the supernova remnant to a distance of more than 20 arcmin. Individual features exhibit apparent motions of 10 to 20 arcsec per year, independently confirmed by near-infrared observations. The observed tangential velocities are at roughly the speed of light. It is likely that the moving structures are infrared echoes, in which interstellar dust is heated by the explosion and by flares from the compact object near the center of the remnant.

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 report MIPS photometry of 20 radio-loud quasars and galaxies at 24 and 70 um (and of five at 160 um). We combine this sample with additional sources detected in the far infrared by IRAS and ISO for a total of 47 objects, including 23 steep spectrum Type I AGNs: radio-loud quasars and broad line radio galaxies; and 24 Type II AGNs: narrow line and weak line FR II radio galaxies. Of this sample, the far infrared emission of all but 3C 380 appears to be dominated by emission by dust heated by the AGN and by star formation. The AGN appears to contribute more than 50% of the far infrared luminosity in most of sources. It is also expected that the material around the nucleus is optically thin in the far infrared. Thus, the measurements at these wavelengths can be used to test the orientation-dependent unification model. As predicted by the model, the behavior of the sources is consistent with the presence of an obscuring circumnuclear torus; in fact, we find it may still have significant optical depth at 24 um.In addition, as expected for the radio-loud quasars, there is a significant correlation between the low frequency radio (178 MHz) and the 70 um emission, two presumably isotropic indicators of nuclear activity. This result is consistent with the simple unified scheme. However, there is a population of radio galaxies that are underluminous at 70 um compared with the radio-loud quasars and hence are a challenge to the simple unified model.

We present high spatial resolution mid- and far-infrared images of the Vega debris disk obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The disk is well resolved and its angular size is much larger than found previously. The radius of the disk is at least 43" (330 AU), 70"(543 AU), and 105" (815 AU) in extent at 24, 70 and 160 um, respectively. The disk images are circular, smooth and without clumpiness at all three wavelengths. The radial surface brightness profiles imply an inner boundary at a radius of 11"+/-2" (86 AU). Assuming an amalgam of amorphous silicate and carbonaceous grains, the disk can be modeled as an axially symmetric and geometrically thin disk, viewed face-on, with the surface particle number density following an r^-1 power law. The disk radiometric properties are consistent with a range of models using grains of sizes ~1 to ~50 um. We find that a ring, containing grains larger than 180 um and at radii of 86-200 AU from the star, can reproduce the observed 850 um flux, while its emission does not violate the observed MIPS profiles. This ring could be associated with a population of larger asteroidal bodies analogous to our own Kuiper Belt. Cascades of collisions starting with encounters amongthese large bodies in the ring produce the small debris that is blown outward by radiation pressure to much larger distances where we detect its thermal emission. The dust production rate is >~10^15 g/s based on the MIPS results. This rate would require a very massive asteroidal reservoir for the dust to be produced in a steady state throughout Vega's life. Instead, we suggest that the disk we imaged is ephemeral and that we are witnessing the aftermath of a large and relatively recent collisional event, and subsequent collisional cascade.

R Coronae Borealis (RCB) stars are hydrogen-deficient, carbon-rich pulsating post-AGB stars that experience massive irregular declines in brightness caused by circumstellar dust formation. The mechanism of dust formation around RCB stars is not well understood. It has been proposed that CO molecules play an important role in cooling the circumstellar gas so that dust may form. We report on a survey for CO in a sample of RCB stars. We obtained H- and K-band spectra including the first and second overtone CO bands for eight RCB stars, the RCB-like star, DY Per and the final-helium-flash star, FG Sge. The first and second overtone CO bands were detected in the cooler (T(eff)<6000 K) RCB stars, Z Umi, ES Aql, SV Sge and DY Per. The bands are not present in the warmer (T(eff)>6000 K) RCB stars, R CrB, RY Sgr, SU Tau, XX Cam. In addition, first overtone bands are seen in FG Sge, a final-helium-flash star that is in an RCB-like phase at present. Effective temperatures of the eight RCB stars range from 4000 to 7250 K. The observed photospheric CO absorption bands were compared to line-blanketed model spectra of RCB stars. As predicted by the models, the CO bands are strongest in the coolest RCB stars and not present in the warmest. No correlation was found between the presence or strength of the CO bands and dust formation activity in the stars.

We report the discovery of a uniquely large excess of 18O in the hydrogen-deficient carbon (HdC) star, HD 137613, based on a spectrum of the first overtone bands of CO at 2.3-2.4 micron in which three strong absorption bands of 12C18O are clearly present. Bands of 12C16O also are present but no bands of 13C16O or 12C17O are seen. We estimate an isotopic ratio 16O/18O ~1. The Solar value of this ratio is ~500. Neither He-core burning nor He-shell flash burning can produce the isotopic ratios of oxygen and carbon observed in HD 137613. However, a remarkable similarity exists between the observed abundances and those found in the outer layers of the broad He-shell of early-AGB stars, soon after the end of He-core burning. It is not known how the outer envelope down to the He-shell could be lost but some mechanism of enhanced mass loss must be involved. HD 137613 may be a post-early-AGB star with the outer layers of the former He-burning shell as its photosphere. The unusual elemental abundances of the HdC stars resemble those of the R Coronae Borealis (RCB) stars, but HdC stars do not produce clouds of dust that produce declines in brightness. None of the other RCB or HdC stars observed shows significant 18O.

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.

We describe the data reduction algorithms for the Multiband Imaging Photometer for Spitzer (MIPS) instrument. These algorithms were based on extensive preflight testing and modeling of the Si:As (24 micron) and Ge:Ga (70 and 160 micron) arrays in MIPS and have been refined based on initial flight data. The behaviors we describe are typical of state-of-the-art infrared focal planes operated in the low backgrounds of space. The Ge arrays are bulk photoconductors and therefore show a variety of artifacts that must be removed to calibrate the data. The Si array, while better behaved than the Ge arrays, does show a handful of artifacts that also must be removed to calibrate the data. The data reduction to remove these effects is divided into three parts. The first part converts the non-destructively read data ramps into slopes while removing artifacts with time constants of the order of the exposure time. The second part calibrates the slope measurements while removing artifacts with time constants longer than the exposure time. The third part uses the redundancy inherit in the MIPS observing modes to improve the artifact removal iteratively. For each of these steps, we illustrate the relevant laboratory experiments or theoretical arguments along with the mathematical approaches taken to calibrate the data. Finally, we describe how these preflight algorithms have performed on actual flight data.

The first six months of flight data from the Multiband Imaging Photometer for Spitzer (MIPS) were used to test MIPS reduction algorithms based on extensive preflight laboratory data and modeling. The underlying approach for the preflight algorithms has been found to be sound, but some modifications have improved the performance. The main changes are scan mirror dependent flat fields at 24 micron, hand processing to remove the time dependent stim flash latents and fast/slow response variations at 70 micron, and the use of asteroids and other sources instead of stars for flux calibration at 160 micron due to a blue ``leak.'' The photometric accuracy of flux measurpements is currently 5%, 10%, and 20% at 24, 70, and 160 micron, respectively. These numbers are expected to improve as more flight data are analyzed and data reduction algorithms refined.

This paper summarizes the evidence of the ultraviolet properties of dust grains found in starburst galaxies. Observations of starburst galaxies clearly show that the 2175 A feature is weak or absent. This can be the result of radiative transfer effects (mixing the dust and stars) or due to dust grains which do not have this feature. Spherical DIRTY radiative transfer models imply that it is not radiative transfer effects, but other radiative transfer models with disk/bulge geometries have found cases where it could be radiative transfer effects. Recent work on the extinction curves in the Magellanic Clouds and Milky Way has revealed that the traditional explanation of low metallicity for the absence of the 2175 A feature in the Small Magellanic Cloud is likely incorrect. The SMC has one sightline with a 2175 A feature and the Milky Way has sightlines without this feature. In addition, where the 2175 A feature is found to be weak or absent in both Magellanic Clouds and the Milky Way, there is evidence for recent star formation. Taking the sum of the radiative transfer modeling of starburst galaxies and the behavior of Local Group extinction curves, it is likely that the dust grains in starburst galaxies intrinsically lack the 2175 A feature.

We critique the method of constructing extinction curves of lensing galaxies using multiply imaged QSOs. If one of the two QSO images is lightly reddened or if the dust along both sightlines has the same properties then the method works well and produces an extinction curve for the lensing galaxy. These cases are likely rare and hard to confirm. However, if the dust along each sightline has different properties then the resulting curve is no longer a measurement of extinction. Instead, it is a measurement of the difference between two extinction curves. This "lens difference curve'' does contain information about the dust properties, but extracting a meaningful extinction curve is not possible without additional, currently unknown information. As a quantitative example, we show that the combination of two Cardelli, Clayton, & Mathis (CCM) type extinction curves having different values of R(V) will produce a CCM extinction curve with a value of R(V) which is dependent on the individual R(V) values and the ratio of V band extinctions. The resulting lens difference curve is not an average of the dust along the two sightlines. We find that lens difference curves with any value of R(V), even negative values, can be produced by a combination of two reddened sightlines with different CCM extinction curves with R(V) values consistent with Milky Way dust (2.1 < R(V) < 5.6). This may explain extreme values of R(V) inferred by this method in previous studies. But lens difference curves with more normal values of R(V) are just as likely to be composed of two dust extinction curves with R(V) values different than that of the lens difference curve. While it is not possible to determine the individual extinction curves making up a lens difference curve, there is information about a galaxy's dust contained in the lens difference curves.

We have obtained Spitzer Space Telescope MIPS observations of 40 F- and G-type common proper motion members of the Scorpius-Centaurus OB Association with ages between 5 and 20 Myr at 24 micron and 70 micron. We report the detection of fourteen objects which possess 24 micron fluxes >30% larger than their predicted photospheres, tentatively corresponding to a disk fraction of >35%, including seven objects which also possess 70 micron excesses >100 times larger than their predicted photospheres. The 24 micron plus 70 micron excess sources possess high fractional infrared luminosities, LIR/L* = 7.0e-4 - 3.0e-3; either they possess optically thin, dusty beta Pictoris-like disks or compact, opaque HD 98800-like disks.

We present Spitzer Space Telescope observations of V1647 Ori, the outbursting source lighting McNeil's nebula, taken near the optical peak of the outburst in early March 2004. The source is easily detected in all Spitzer imaging bands from 3.6 - 70 microns. The fluxes at all wavelengths are roughly a factor of 15 brighter than pre-outburst levels; we measure a bolometric luminosity of 44 Lsun. We posit that this event is due to an increase in the accretion luminosity of the source. Simple models of an accretion disk plus tenuous envelope can qualitatively explain the observed pre- and post-outburst spectral energy distributions. The accretion activity implied by our results indicates that the outburst may be intermediate between FUor and EXor-type events. We also report the discovery of a previously unknown mid-infrared counterpart to the nearby Herbig-Haro object HH 22.

As part of a large Spitzer/MIPS GTO program, we have searched for infrared excesses due to debris disks toward 26 FGK field stars known from radial velocity (RV) studies to have one or more planets. While none of these stars show excesses at 24um, we have detected 70um excesses around 6 stars at the 3-sigma confidence level. The excesses are produced by cool material (< 100 K) located beyond 10 AU, well outside the "habitable zones" of these systems and consistent with the presence of Kuiper Belt analogues with ~100 times more emitting surface area than in our own planetary system. These planet-bearing stars are, by selection for RV studies, typically older than 1 Gyr, and the stars identified here with excesses have a median age of 4 Gyr. We find a preliminary correlation of both the frequency and the magnitude of dust emission with the presence of known planets. These are the first stars outside the solar system identified as having both well-confirmed planetary systems and well-confirmed IR excesses.

A large amount (about three solar masses) of cold (18 K) dust in the prototypical type II supernova remnant Cassiopeia A was recently reported. It was concluded that dust production in type II supernovae can explain how the large quantities (10^8 solar masses) of dust observed in the most distant quasars could have been produced within only 700 million years after the Big Bang. Foreground clouds of interstellar material, however, complicate the interpretation of the earlier submillimetre observations of Cas A. Here we report far-infrared and molecular line observations that demonstrate that most of the detected submillimetre emission originates from interstellar dust in a molecular cloud complex located in the line of sight between the Earth and Cas A, and is therefore not associated with the remnant. The argument that type II supernovae produce copious amounts of dust is not supported by the case of Cas A, which previously appeared to provide the best evidence for this possibility.

Following our initial discovery of blue luminescence in the spectrum of the Red Rectangle (RR) and its identification as fluorescence by small three- to four-ringed polycyclic aromatic hydrocarbon (PAH) molecules, we report on the spatial correlation between the blue luminescence and the 3.3 micron emission, commonly attributed to small, neutral PAH molecules, and on the newly-derived UV/optical attenuation curve for the central source of the RR, HD 44179. Both results provide strong additional evidence for the presence of small PAH molecules with masses of less than 250 a.m.u. in the RR, which supports the attribution of the blue luminescence to fluorescence by the same molecules. We contrast the excellent spatial correlation of the two former emissions with the distinctly different spatial distribution of the extended red emission (ERE) and of the dust-scattered light within the RR. The UV/optical attenuation curve of the central star is unlike any interstellar extinction curve and is interpreted as resulting from circumstellar opacity alone. Major contributions to this opacity are absorptions in broad bands in the mid-UV, contributing to the electronic excitation of the luminescing PAH molecules, and a sharp ionization discontinuity near 7.5 eV in the far-UV, which places a sharp upper limit on the masses of the PAH molecules that are responsible for this absorption. The strength of the far-UV absorption leads to an abundance of the PAH molecules of 10^{-5} relative to hydrogen in the RR. Such small PAHs are perhaps unique to the environment in the RR, where they are shielded from harsh radiation by the dense circmstellar disk.

We present results of new Monte Carlo calculations made with the DIRTY code of radiative transfer of stellar and scattered radiation for a dusty giant late-type galaxy like the Milky Way, which illustrate the effect of the attenuation of stellar light by internal dust on the integrated photometry of the individual bulge and disk components. Here we focus on the behavior of the attenuation function, the color excess, and the fraction of light scattered or directly transmitted towards the outside observer as a function of the total amount of dust and the inclination of the galaxy, and the structure of the dusty interstellar medium (ISM) of the disk. We confirm that dust attenuation produces qualitatively and quantitatively different effects on the integrated photometry of bulge and disk, whatever the wavelength. In addition, we find that the structure of the dusty ISM affects more sensitively the observed magnitudes than the observed colors of both bulge and disk. Finally, we show that the contribution of the scattered radiation to the total monochromatic light received by the outside observer is significant, particularly at UV wavelengths, even for a two-phase, clumpy, dusty ISM. Thus understanding dust scattering properties is fundamental for the interpretation of extragalactic observations in the rest-frame UV.

We have assembled a homogeneous database of 417 ultraviolet (UV) extinction curves for reddened sightlines having International Ultraviolet Explorer (IUE) spectra. We have combined these with optical and 2MASS photometry allowing estimates of the ratio of total-to-selective extinction, R(V), for the entire sample. Fitzpatrick-Massa (FM) parameters have also been found for the entire sample. This is the largest study of parameterized UV extinction curves yet published and it covers a wide range of environments, from dense molecular clouds to the diffuse interstellar medium (ISM), with extinctions A(V) ranging from 0.50 to 4.80. It is the first to extend far beyond the solar neighborhood and into the Galaxy at large, with 30 sightlines having distances > 5 kpc. Previously, the longest sightlines with FM parameters and R(V) extended ~ 1 kpc. We find that (1.) the CCM extinction law applies for 93% of the sightlines, implying that dust processing in the Galaxy is efficient and systematic; (2.) the central wavelength of the 2175 A bump is constant; (3.) the 2175 A bump width is dependent on environment. Only four sightlines show systematic deviations from CCM, HD 29647, 62542, 204827, and 210121. These sightlines all sample dense, molecule-rich clouds. The new extinction curves and values of R(V) allow us to revise the CCM law.

The Spitzer Space Telescope was used to study the mid-infrared to far-infrared properties of NGC 300, and to compare dust emission to Halpha to elucidate the heating of the ISM and the star formation cycle at scales < 100 pc. The new data allow us to discern clear differences in the spatial distribution of 8 micron dust emission with respect to 24 micron dust and to HII regions traced by the Halpha light. The 8 micron emission highlights the rims of HII regions, and the 24 micron emission is more strongly peaked in star forming regions than at 8 microns. We confirm the existence and approximate amplitude of interstellar dust emission at 4.5 microns, detected statistically in Infrared Space Observatory (ISO) data, and conclude it arises in star forming regions. When averaging over regions larger than ~ 1 kpc, the ratio of Halpha to Aromatic Feature emission in NGC 300 is consistent with the values observed in disks of spiral galaxies. The mid-to-far-infrared spectral energy distribution of dust emission is generally consistent with pre-Spitzer models.

We have observed the reflection nebula NGC 7023, with the Short-High module and long-slit Short-Low and Long-Low modules of the Infrared Spectrograph on the Spitzer Space Telescope. We also present Infrared Array Camera and Multiband Imaging Photometer for Spitzer images of NGC 7023 at 3.6, 4.5, 8.0, and 24 microns. We observe the aromatic emission features (AEFs) at 6.2, 7.7, 8.6, 11.3, and 12.7 microns, plus a wealth of weaker features. We find new unidentified interstellar emission features at 6.7, 10.1, 15.8, 17.4, and 19.0 microns. Possible identifications include aromatic hydrocarbons or nanoparticles of unknown mineralogy. We see variations in relative feature strengths, central wavelengths, and feature widths, in the AEFs and weaker emission features, depending on both distance from the star and nebular position (southeast vs. northwest).

We derive galaxy source counts at 70 and 160 microns using the Multiband Imaging Photometer for Spitzer (MIPS) to map the Chandra Deep Field South (CDFS) and other fields. At 70 microns, our observations extend upwards about 2 orders of magnitude in flux density from a threshold of 15 mJy, and at 160 microns they extend about an order of magnitude upward from 50 mJy. The counts are consistent with previous observations on the bright end. Significant evolution is detected at the faint end of the counts in both bands, by factors of 2-3 over no-evolution models. This evolution agrees well with models that indicate most ofthe faint galaxies lie at redshifts between 0.7 and 0.9. The new Spitzer data already resolve about 23% of the Cosmic Far Infrared Background at 70 microns and about 7% at 160 microns.

We present Spitzer IRAC (3.6, 4.5, 5.8, 8.0 um) and MIPS (24 um) observations of the SCUBA submillimeter sources and uJy VLA radio sources in a 5'x5' area in the Lockman Hole East region. Out of the ~40 SCUBA/VLA sources in the field, Spitzer counterparts were detected for nearly all except for the few low-weight SCUBA detections. We show that the majority (80-90%) of the detected sources are cold (i.e., starburst-like) infrared-luminous galaxies (L_{IR} > 10^{11} Lsun) at redshift 0.5<z<3.5, whose star-formation rate density (SFRD) is comparable to that of the optically-selected star-forming galaxies.

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 3-160 micron photometry obtained with the IRAC and MIPS instruments for the first five targets from the Spitzer Legacy Science Program "Formation and Evolution of Planetary Systems" and 4-35 micron spectro-photometry obtained with the IRS for two sources. We discuss in detail our observations of the debris disks surrounding HD 105 (G0V, 30 +- 10 Myr) and HD 150706 (G3V, ~ 700 +- 300 Myr). For HD 105, possible interpretations include large bodies clearing the dust inside of 45 AU or a reservoir of gas capable of sculpting the dust distribution. The disk surrounding HD 150706 also exhibits evidence of a large inner hole in its dust distribution. Of the four survey targets without previously detected IR excess, spanning ages 30 Myr to 3 Gyr, the new detection of excess in just one system of intermediate age suggests a variety of initial conditions or divergent evolutionary paths for debris disk systems orbiting solar-type stars.

We investigate the luminous X-ray sources in the Lockman Hole (LH) and the Extended Groth Strip (EGS) detected at 24microns using MIPS and also with IRAC on board Spitzer. We assemble optical/infrared spectral energy distributions (SEDs) for 45 X-ray/24micron sources in the EGS and LH. Only about 1/4 of the hard X-ray/24micron sources show pure type 1 AGN SEDs. More than half of the X-ray/24micron sources have stellar-emission-dominated or obscured SEDs, similar to those of local type 2 AGN and spiral/starburst galaxies. One-third of the sources detected in hard X-rays do not have a 24micron counterpart. Two such sources in the LH have SEDs resembling those of S0/elliptical galaxies. The broad variety of SEDs in the optical-to-Spitzer bands of X-ray selected AGN means that AGN selected according to the behavior in the optical/infrared will have to be supplemented by other kinds of data (e.g., X-ray) to produce unbiased samples of AGN.

We present images of the galaxy NGC 55 at 24, 70, and 160 micron obtained with the Multiband Imaging Photometer for Spitzer (MIPS) instrument aboard the Spitzer Space Telescope. The new images display the far infrared emission in unprecedented detail and demonstrate that the infrared morphology differs dramatically from that at shorter wavelengths. The most luminous emission region in the galaxy is marginally resolved at 24 micron and has a projected separation of nearly 520 pc from the peak emission in the optical and near infrared. This region is responsible for ~9% of the total emission at 24 micron and is likely a young star formation region. We show that this and other compact sources account for more than 1/3 of the total 24 micron emission. We compute a total infrared luminosity for NGC 55 of 1.2*10^9 L_sun. The star formation rate implied by our measurements is 0.22 M_sun/yr. We demonstrate that the cold dust is more extended than the warm dust in NGC 55--the minor-axis scale heights are 0.32, 0.43, and 0.49 kpc at 24, 70 and 160 micron, respectively. The dust temperature map shows a range of temperatures that are well-correlated with the 24 micron surface brightness, from 20 K in low-surface-brightness regions to 26 K in high-surface-brightness regions.

We present Spitzer observations in five wavebands between 3.6 and 24um of an unbiased sample of 9 luminous, dusty galaxies selected at 1200um by the MAMBO camera on the IRAM 30-m telescope, a population akin to the well-known submm or `SCUBA' galaxies (hereafter SMGs). Owing to the coarse resolution of submm/mm instrumentation, SMGs have traditionally been difficult to identify at other wavelengths. We compare our multi-wavelength catalogs to show that the overlap between 24 and 1200um must be close to complete at these flux levels. We find that all (4/4) of the most secure >=4sigma SMGs have robust >=4sigma counterparts at 1.4GHz, while the fraction drops to 7/9 using all >=3sigma SMGs. We show that combining mid-IR and marginal (>=3sigma) radio detections provides plausible identifications in the remaining cases, enabling us to identify the complete sample. Accretion onto an obscured central engine is betrayed by the shape of the mid-IR continuum emission for several sources, confirming Spitzer's potential to weed out active galaxies. We demonstrate the power of a S(24um)/S(8um) vs S(8um)/S(4.5um) color-color plot as a diagnostic for this purpose. However, we conclude that the majority (~75%) of SMGs have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. Sensitive 24-um observations are clearly a useful route to identify and characterize reliable counterparts to high-redshift far-IR-bright galaxies, complementing what is possible via deep radio imaging.

We present preliminary results on 24micron detections of luminous infrared galaxies at z>1 with the Multiband Imaging Photometer for Spitzer (MIPS). Observations were performed in the Lockman Hole and the Extended Groth Strip (EGS), and were supplemented by data obtained with the Infrared Array Camera (IRAC) between 3 and 9microns. The positional accuracy of ~2arcsec for most MIPS/IRAC detections provides unambiguous identifications of their optical counterparts. Using spectroscopic redshifts from the Deep Extragalactic Evolutionary Probe survey, we identify 24micron sources at z>1 in the EGS, while the combination of the MIPS/IRAC observations with $BVRIJHK$ ancillary data in the Lockman Hole also shows very clear cases of galaxies with photometric redshifts at 1<z<2.5. The observed 24micron fluxes indicate infrared luminosities greater than 10^11 L_sol, while the data at shorter wavelengths reveal rather red and probably massive (M>=M*) galaxy counterparts. It is the first time that this population of luminous objects is detected up to z~2.5 in the infrared. Our work demonstrates the ability of the MIPS instrument to probe the dusty Universe at very high redshift, and illustrates how the forthcoming Spitzer deep surveys will offer a unique opportunity to illuminate a dark side of cosmic history not explored by previous infrared experiments.

We present observations of the star formation region NGC 7129 taken with the Multiband Imaging Photometer for Spitzer (MIPS). A significant population of sources, likely pre-main sequence members of the young stellar cluster, is revealed outside the central photoionization region. Combining with Infrared Array Camera (IRAC) and ground-based near-infrared images, we have obtained colors and spectral energy distributions for some 60 objects. The [3.6]-[4.5] vs. [8]-[24] color-color plane shows sources clustered at several different loci, which roughly correspond to the archetypal evolutionary sequence Class 0, I, II, and III. We obtain preliminary classifications for 36 objects, and find significant numbers of both Class I and II objects. Most of the pre-main sequence candidates are associated with the densest part of the molecular cloud surrounding the photoionization region, indicating active star formation over a broad area outside the central cluster. We discuss three Class II candidates that exhibit evidence of inner disk clearing, which would be some of the youngest known examples of a transition from accretion to optically thin quiescent disks.

IRAC imaging of a 4'7x4'7 area in the Lockman Hole detected over 400 galaxies in the IRAC 3.6 micron and 4.5 micron bands, 120 in the 5.8 micron, and 80 in the 8 micron bandin 30 minutes of observing time. Color-color diagrams suggest that about half of these galaxies are at redshifts 0.6<z<1.3 with about a quarter at higher redshifts (z>1.3). We also detect IRAC counterparts for 6 of the 7 SCUBA sources and all 9 XMM sources in this area. The detection of the counterparts of the SCUBA sources and galaxies at z>1.3 demonstrates the ability of IRAC to probe the universe at very high redshifts.

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.

Sixty three members of the 100 Myr old open cluster M47 (NGC 2422) have been detected at 24 micron with Spitzer. The Be star V 378 Pup shows an excess both in the near-infrared and at 24 micron (K-[24] = 2.4 mag), probably due to free-free emission from the gaseous envelope. Seven other early-type stars show smaller excesses, K-[24] = 0.6-0.9. Among late-type stars, two show large excesses: P922 - a K1V star with K-[24] = 1.08 pm 0.11 and P1121 - an F9V star with K-[24] = 3.72 pm 0.02. P1121 is the first known main-sequence star showing an excess comparable to that of beta Pic, which may indicate the presence of an exceptionally massive debris disk. It is possible that a major planetesimal collision has occurred in this system, consistent with the few hundred Myr time scales estimated for the clearing of the Solar System.

We present an updated phenomenological galaxy evolution model to fit the Spitzer 24, 70 and 160 microns number counts as well as all the previous mid and far infrared observations. Only a minor change of the co-moving luminosity density distribution in the previous model (Lagache, Dole, Puget 2003), combined with a slight modification of the starburst template spectra mainly between 12 and 30 microns, are required to fit all the data available. We show that the peak in the MIPS 24 micron counts is dominated by galaxies with redshift between 1 and 2, with a non negligible contribution from the z>2 galaxies (~30% at S=0.2 mJy). The very close agreement between the model and number counts at 15 and 24 microns strikingly implies that (1) the PAHs (Policyclic Aromatic Hydrocarbons) features remain prominent in the redshift band 0.5 to 2.5 and (2) the IR energy output has to be dominated by ~3 10^11 Lo to ~3 10^12 Lo galaxies from redshift 0.5 to 2.5. Combining Spitzer with the Infrared Space Observatory (ISO) deep cosmological surveys gives for the first time an unbiased view of the infrared Universe from z=0 to z=2.5.

We use the source counts measured with the Multiband Imaging Photometer for Spitzer (MIPS) at 24, 70, and 160 microns to determine the 5-sigma confusion limits due to extragalactic sources: 56 micro-Jy, 3.2 and 40 mJy at 24, 70 and 160 microns, respectively. We also make predictions for confusion limits for a number of proposed far infrared missions of larger aperture (3.5 to 10m diameter).

We examine the 24 micron to X-ray color of 157 X-ray--selected AGN as a function of X-ray obscuration and optical classification in the Chandra Deep Field South. The sample consists of the Chandra hard--band detections with hard flux >10^-15 erg/s/cm^2. A deep 24 micron mosaic obtained with Spitzer provides mid--infrared fluxes for the sample. Since obscured AGN locally have higher 24 micron / 2-8 keV ratios than unobscured AGN, and since X-ray background models predict a large population of obscured AGN, we expect to find many X-ray--hard, IR--bright AGN. Instead, we find that the 24 micron to X-ray flux ratio does not depend on X-ray hardness in the full sample, nor does it differ between narrow--line and broad--line AGN. We identify 5 nearly Compton--thick AGN, and find they have similar 24 micron to X-ray flux ratios compared to the full sample. We consider AGN in the narrow redshift spikes at z~0.7; for these AGN, there is some evidence that the flux ratio increases with X-ray hardness. The redshift slice also shows an odd trend, which is also prominent in the full sample: a group of X-ray--hard AGN with very low 24 micron to X-ray flux ratios.

We present sub-millimeter statistical detections of galaxies discovered in the 5'x5' Spitzer Early Release Observations (to 4-15 microJy 5 sigma at 3.6-8 microns, 170 microJy at 24 microns) through a stacking analysis of our reanalysed SCUBA 8mJy survey maps, and a Spitzer identification of a new sub-millimeter point source in the 8mJy survey region. For sources detected at 5.8 or 8 microns (154 and 111 sources respectively), we detect positive skews in the sub-millimeter flux distributions at 99.2-99.8% confidence using Kolmogorov-Smirnov tests, at both 850 microns and 450 microns. We also marginally detect the Spitzer 24 micron galaxies at 850 microns at 97% confidence, and place limits on the mean sub-millimeter fluxes of the 3.6 and 4.5 micron sources. Integrating the sub-millimeter fluxes of the Spitzer populations, we find the 5.8 micron galaxies contribute 0.12 +/- 0.05 nW/m^2/sr to the 850 micron background, and 2.4 +/- 0.7 nW/m^2/sr to the 450 micron background; similar contributions are made by the 8 micron-selected sample. We infer that the populations dominating the 5.8 and 8 micron extragalactic background light also contribute around a quarter of the 850 micron background and the majority of the 450 micron background.

Galaxy source counts in the infrared provide strong constraints on the evolution of the bolometric energy output from distant galaxy populations. We present the results from deep 24 micron imaging from Spitzer surveys, which include approximately 50,000 sources to an 80% completeness of 60 uJy. The 24 micron counts rapidly rise at near-Euclidean rates down to 5 mJy, increase with a super-Euclidean rate between 0.4 - 4 mJy, and converge below 0.3 mJy. The 24 micron counts exceed expectations from non-evolving models by a factor >10 at 0.1 mJy. The peak in the differential number counts corresponds to a population of faint sources that is not expected from predictions based on 15 micron counts from ISO. We argue that this implies the existence of a previously undetected population of infrared-luminous galaxies at z ~ 1-3. Integrating the counts to 60 uJy, we derive a lower limit on the 24 micron background intensity of 1.9 +/- 0.6 nW m-2 sr-1 of which the majority (~ 60%) stems from sources fainter than 0.4 mJy. Extrapolating to fainter flux densities, sources below 60 uJy contribute 0.8 {+0.9/-0.4} nW m-2 sr-1 to the background, which provides an estimate of the total 24 micron background of 2.7 {+1.1/-0.7} nW m-2 sr-1.

We investigate Extremely Red Objects (EROs) using near- and mid-infrared observations in five passbands (3.6 to 24 micron) obtained from the Spitzer Space Telescope, and deep ground-based R and K imaging. The great sensitivity of the IRAC camera allows us to detect 64 EROs in only 12 minutes of IRAC exposure time, by means of an R-[3.6] color cut (analogous to the traditional red R-K cut). A pure infrared K-[3.6] red cut detects a somewhat different population and may be more effective at selecting z > 1.3 EROs. We find 17% of all galaxies detected by IRAC at 3.6 or 4.5 micron to be EROs. These percentages rise to about 40% at 5.8 micron, and about 60% at 8.0 micron. We utilize the spectral bump at 1.6 micron to divide the EROs into broad redshift slices using only near-infrared colors (2.2/3.6/4.5 micron). We conclude that two-thirds of all EROs lie at redshift z > 1.3. Detections at 24 micron imply that at least 11% of 0.6 < z < 1.3 EROs and at least 22% of z > 1.3 EROs are dusty star-forming galaxies.

Here we report our discovery of a band of blue luminescence (BL) in the Red Rectangle (RR) nebula. This enigmatic proto-planetary nebula is also one of the brightest known sources of extended red emission as well as of unidentified infra-red (UIR) band emissions. The spectrum of this newly discovered BL is most likely fluorescence from small neutral polycyclic aromatic hydrocarbon (PAH) molecules. PAH molecules are thought to be widely present in many interstellar and circumstellar environments in our galaxy as well as in other galaxies, and are considered likely carriers of the UIR-band emission. However, no specific PAH molecule has yet been identified in a source outside the solar system, as the set of mid-infra-red emission features attributed to these molecules between the wavelengths of 3.3 micron and 16.4 micron is largely insensitive to molecular sizes. In contrast, near-UV/blue fluorescence of PAHs is more specific as to size, structure, and charge state of a PAH molecule. If the carriers of this near-UV/blue fluorescence are PAHs, they are most likely neutral PAH molecules consisting of 3-4 aromatic rings such as anthracene (C14H10) and pyrene (C16H10). These small PAHs would then be the largest molecules specifically identified in the interstellar medium.

The wavelength dependence of interstellar polarization due to dust in M31 has been observed along four sightlines. Only one sightline had been measured previously.The globular clusters, S78, S150, S233 and Baade 327 were used as point sources to probe the interstellar dust in M31. The Serkowski law produces good fits for all the sightlines although the relationship between K and lambda(max) may be different from that found in the Galaxy. The results of this study imply that the slope K/lambda(max) may be significantly larger in M31. The Serkowski curves are significantly narrower than those of the same lambda(max) in the Galaxy and may require extreme modifications to the size distributions of silicate particles. The fits for the four sightlines reveal values of lambda(max) ranging from 4800 to 5500 A. These are consistent with average values of lambda(max) measured in the Galaxy and the Magellanic Clouds. The range measured for M31 corresponds to R(V) values of 2.7 to 3.1. The range in R(V) seen in the Galaxy is 2.5 to 5.5 implying, for this small sample, that the average size ofinterstellar grains in M31 is typically smaller than that seen for Galactic grains if the nature of the grains is the same. Also, the polarization efficiency for these sightlines is large although some bias is expected since sightlines known to have significant interstellar polarization were selected for the sample.

Studies of dust scattering properties in astrophysical objects with Milky Way interstellar dust are reviewed. Such objects are reflection nebulae, dark clouds, and the Diffuse Galactic Light (DGL). To ensure their basic quality, studies had to satisfy four basic criteria to be included in this review. These four criteria significantly reduced the scatter in dust properties measurements, especially in the case of the DGL. Determinations of dust scattering properties were found to be internally consistent for each object type as well as consistent between object types. The 2175 A bump is seen as an absorption feature. Comparisons with dust grain models find general agreement with significant disagreements at particular wavelengths (especially in the far-ultraviolet). Finally, unanswered questions and future directions are enumerated.

It is well known that the sightline toward HD 204827 in the cluster Trumpler 37 shows an UV extinction curve that does not follow CCM. However, when a dust component, foreground to the cluster, is removed, the residual extinction curve is identical to that found in the SMC within the uncertainties. The curve is very steep and has little or no 2175 A bump. The position of HD 204827 in the sky is projected onto the edge of the Cepheus IRAS Bubble. In addition, HD 204827 has an IRAS bowshock indicating that it may be embedded in dust swept up by the supernova that created the IRAS Bubble. Shocks due to the supernova may have led to substantial processing of this dust. The HD 204827 cloud is dense and rich in carbon molecules. The 3.4 $\mu$m feature indicating a C-H grain mantle is present in the dust toward HD 204827. The environment of the HD 204827 cloud dust may be similar to the dust associated with HD 62542 which lies on the edge of a stellar wind bubble and is also dense and rich in molecules. This sightline may be a Rosetta Stone if its environment can be related to those in the SMC having similar dust.

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.

We present an exhaustive, quantitative comparison of all of the known extinction curves in the Small and Large Magellanic Clouds (SMC and LMC) with our understanding of the general behavior of Milky Way extinction curves. The R_V dependent CCM relationship and the sample of extinction curves used to derive this relationship is used to describe the general behavior of Milky Way extinction curves. The ultraviolet portion of the SMC and LMC extinction curves are derived from archival IUE data, except for one new SMC extinction curve which was measured using HST/STIS observations. The optical extinction curves are derived from new (for the SMC) and literature UBVRI photometry (for the LMC). The near-infrared extinction curves are calculated mainly from 2MASS photometry supplemented with DENIS and new JHK photometry. For each extinction curve, we give R_V = A(V)/E(B-V) and N(HI) values which probe the same dust column as the extinction curve. We compare the properties of the SMC and LMC extinction curves with the CCM relationship three different ways: each curve by itself, the behavior of extinction at different wavelengths with R_V, and behavior of the extinction curve FM fit parameters with R_V. As has been found previously, we find that a small number of LMC extinction curves are consistent with the CCM relationship, but majority of the LMC and all of the SMC curves do not follow the CCM relationship. For the first time, we find that the CCM relationship seems to form a bound on the properties of all of the LMC and SMC extinction curves. This result strengthens the picture of dust extinction curves exhibit a continuum of properties between those found in the Milky Way and the SMC Bar. (abridged)

We have obtained new STIS/HST spectra to search for structure in the ultraviolet interstellar extinction curve, with particular emphasis on a search for absorption features produced by polycyclic aromatic hydrocarbons (PAHs). The presence of these molecules in the interstellar medium has been postulated to explain the infrared emission features seen in the 3-13 $\mu$m spectra of numerous sources. UV spectra are uniquely capable of identifying specific PAH molecules. We obtained high S/N UV spectra of stars which are significantly more reddened than those observed in previous studies. These data put limits on the role of small (30-50 carbon atoms) PAHs in UV extinction and call for further observations to probe the role of larger PAHs. PAHs are of importance because of their ubiquity and high abundance inferred from the infrared data and also because they may link the molecular and dust phases of the interstellar medium. A presence or absence of ultraviolet absorption bands due to PAHs could be a definitive test of this hypothesis. We should be able to detect a 20 \AA wide feature down to a 3$\sigma$ limit of $\sim$0.02 A$_V$. No such absorption features are seen other than the well-known 2175 \AA bump.

Employing the Maximum Entropy Method algorithm, we fit interstellar extinction measurements which span the wavelength range 0.125-3 micron. We present a uniform set of MEM model fits, all using the same grain materials, optical constants and abundance constraints. In addition, we are taking advantage of improved UV and IR data and better estimates of the gas-to-dust ratio. The model fits cover the entire range of extinction properties that have been seen in the Galaxy and the Magellanic Clouds. The grain models employed for this presentation are the simplistic homogeneous spheres models (i.e., Mathis, Rumpl, & Nordsieck 1977) with two (graphite, silicate) or three (graphite, silicate, amorphous carbon) components. Though such usage is only a first step, the results do provide interesting insight into the use of grain size as a diagnostic of dust environment. We find that the SMC Bar extinction curve cannot be fit using carbon grains alone. This is a challenge to the recent observational result indicating little silicon depletion in the SMC.

We present our analysis of UV attenuation by internal dust of a large sample (N=906 galaxies) of Lyman Break Galaxies (LBGs). Using spectral energy distributions (SEDs) from the P\'EGASE galaxy spectral evolution model we apply dust attenuation corrections to the G-R colors using the Witt & Gordon (2000) models for radiative transfer in dusty galactic environments to arrive at the UV attenuation factors. We show that the dust in the LBGs exhibit SMC-like characteristics rather than MW-like, and that the dust geometry in these systems is most likely to be represented by a clumpy shell configuration. We show that the attenuation factor exhibits a pronounced dependence on the luminosity of the LBG, $a_{1600}\propto (L/L_\sun)^\alpha$, where $0.5\leq\alpha\leq1.5$. The exponent $\alpha$ depends on the initial parameters of the stellar population chosen to model the galaxies and the dust properties. We find that the luminosity weighted average attenuation factor is likely to be in the range from $5.7-18.5$, which is consistent with the upper limits to the star formation rate at $2<z<4$ set by the FIR background. This implies that the current UV/optical surveys do detect the bulk of the star formation during the epoch $2<z<4$, but require substantial correction for internal dust attenuation.

We have analyzed observations of the Acidalia hemisphere of Mars taken by the Hubble Space Telescope's Near-Infrared Camera Multi-Object Spectrograph (HST/NICMOS) during July of 1997 (Ls = 152 degrees, northern Martian summer). The data consist of images at ~60 km/pixel resolution, using both narrow- and medium-band filters specifically selected to allow us to study the hydration state of the Martian surface. Calibration was performed by comparison to Phobos-2 ISM observations of overlapping regions, and atmospheric gas correction was performed by modeling the atmosphere for each pixel using a line-by-line radiative transfer code coupled with the MOLA altimetry data. Our results indicate the presence of at least three spectrally different large-scale (>1000 km diameter) terrains corresponding to the dark regions of northern Acidalia, the southern hemisphere classical dark terrain, and the classical intermediate terrain adjacent to southern Acidalia. We also identified two other spectrally unique terrains, corresponding to the northern polar ice cap, and to the southern winter polar hood. Comparisons with mineral spectra indicate the possibility of different H2O- or OH-bearing (i.e., hydroxides and/or hydrates) minerals existing both in northern Acidalia and in the nearby intermediate albedo terrain. Hydrated minerals do not appear to be spectrally important components of the southern hemisphere dark terrains imaged by HST in 1997.

Space Telescope Imaging Spectrograph (STIS) and archival International Ultraviolet Explorer (IUE) data were used to search for structure in the ultraviolet interstellar extinction curve. Between 1250 and 3200 A, in addition to the well-known 2175 A bump, there is only one fairly small (10%) absorption feature at 2700 A. The feature is very broad (>400 A). It coincides extremely well in position and shape with laboratory measurements of presolar diamond grains found in meteorites. As this feature appears in both the new STIS and archival IUE data, it is unlikely to be an instrumental effect. This feature appears to be ubiquitous in the ISM as it is seen along sightlines toward various clouds with reddenings in the range of E(B-V) = 0.35 to 1.73 mag. The identification with diamonds remains tentative since it is based on a single, weak feature. If the feature is due to interstellar diamonds then they require an order of magnitude more absorption strength per gram/cm^2 than the presolar meteoritic diamonds measured in the lab.

A stellar wind module has been developed for the PHOENIX stellar atmosphere code for the purpose of computing non-LTE, line-blanketed, expanding atmospheric structures and detailed synthetic spectra of hot luminous stars with winds. We apply the code to observations of Deneb, for which we report the first positive detections of mm and cm emission (obtained using the SCUBA and the VLA), as well a strong upper limit on the 850 micron flux (using the HHT). The slope of the radio spectrum shows that the stellar wind is partially ionized. We report a uniform-disk angular diameter measurement, 2.40 +/- 0.06 mas, from the Navy Prototype Optical Interferometer (NPOI). The measured bolometric flux and corrected NPOI angular diameter yield an effective temperature of 8600 +/- 500 K. Least-squares comparisons of synthetic spectral energy distributions from 1220 A to 3.6 cm with the observations provide estimates for the effective temperature and the mass-loss rate of 8400 +/- 100 K and 8 +/- 3 E-7 M_sun/yr, respectively. This range of mass-loss rates is consistent with that derived from high dispersion UV spectra when non-LTE metal-line blanketing is considered. We are unable achieve a reasonable fit to a typical Halpha P-Cygni profile with any model parameters over a reasonable range. This is troubling because the \ha profile is the observational basis for Wind Momentum-Luminosity Relationship.

We investigate the effects of dust on Large Magellanic Cloud (LMC) HII region spectral energy distributions using arcminute-resolution far-ultraviolet (FUV), H-alpha, far-infrared (FIR), and radio images. Widely-used indicators of the amount of light lost to dust (attenuation) at H-alpha and in the FUV correlate with each other, although often with substantial scatter. There are two interesting systematic discrepancies. First, H-alpha attenuations estimated from the Balmer decrement are lower than those estimated from the H-alpha-to-thermal radio luminosity ratio. Our data, at this stage, cannot unambiguously identify the source of this discrepancy. Second, the attenuation at 1500 angstroms and UV spectral slope, beta, correlate, although the slope and scatter are substantially different from the correlation first derived for starbursting galaxies by Calzetti et al. Combining our result with those of Meurer et al. for ultra-luminous infrared galaxies and Calzetti et al. for starbursting galaxies, we conclude that no single relation between beta and 1500 angstrom attenuation is applicable to all star-forming systems.

We present UV and optical spectropolarimetry of two highly polarized IRAS-selected QSOs, IRAS 13349+2438 and the BALQSO IRAS 14026+4341. The polarization in both objects rises rapidly toward the blue, peaks near 3000 A in the rest frame and remains nearly constant for shorter wavelengths. The rest frame optical polarized flux density spectra also increase rapidly towards the blue, but then decrease dramatically below 3000 A. This distinctive wavelength dependence of polarized flux shows that the polarization is produced by dust scattering. As for many Seyfert, radio and Hyperluminous Infrared Galaxies (HIGs), the lower polarization of the weak [O III] 4959,5007 lines in IRAS 13349+2438 suggests that the scattering grains lie interior to, or mixed with the narrow line gas. We construct full radiative transfer models of these systems consisting of a dusty sphere of modest optical depth illuminated axisymmetrically from within by a powerlaw QSO spectrum. We show that this simple model successfully reproduces the qualitative polarization properties of the objects. Despite similarities to other IRAS-selected BALQSOs, our FOS spectropolarimetry of IRAS 13349+2438 does not reveal broad absorption lines. IRAS 14026+4341 has an Al III BAL in both scattered and total flux density. We discuss these two objects in terms of both orientation and evolutionary unified schemes for QSOs, BALQSOs and HIGs.

We report the discovery of eight new R Coronae Borealis (RCB) stars in the Large Magellanic Cloud (LMC) using the MACHO project photometry database. The discovery of these new stars increases the number of known RCB stars in the LMC to thirteen. We have also discovered four stars similar to the Galactic variable DY Per. These stars decline much more slowly and are cooler than the RCB stars. The absolute luminosities of the Galactic RCB stars are unknown since there is no direct measurement of the distance to any Galactic RCB star. Hence, the importance of the LMC RCB stars. We find a much larger range of absolute magnitudes (M(V) = -2.5 to -5 mag) than inferred from the small pre-MACHO sample of LMC RCB stars. It is likely that there is a temperature - M(V)relationship with the cooler stars being intrinsically fainter. Cool (~5000 K) RCB stars are much more common than previously thought based on the Galactic RCB star sample. Using the fairly complete sample of RCB stars discovered in the MACHO fields, we have estimated the likely number of RCB stars in the Galaxy to be ~3,200. The SMC MACHO fields were also searched for RCB stars but none were found.

We present the DIRTY radiative transfer model in this paper and a companion paper. This model computes the polarized radiative transfer of photons from arbitrary distributions of stars through arbitrary distributions of dust using Monte Carlo techniques. The dust re-emission is done self-consistently with the dust absorption and scattering and includes all three important emission paths: equilibrium thermal emission, non-equilibrium thermal emission, and the aromatic features emission. The algorithm used for the radiative transfer allows for the efficient computation of the appearance of a model system as seen from any viewing direction. We present a simple method for computing an upper limit on the output quantity uncertainties for Monte Carlo radiative transfer models which use the weighted photon approach.

In this paper and a companion paper we present the DIRTY model, a Monte Carlo radiative transfer code, self-consistently including dust heating and emission, and accounting for the effects of the transient heating of small grains. The code is completely general; the density structure of the dust, the number and type of heating sources, and their geometric configurations can be specified arbitrarily within the model space. Source photons are tracked through the scattering and absorbing medium using Monte Carlo techniques and the effects of multiple scattering are included. The dust scattering, absorbing, and emitting properties are calculated from realistic dust models derived by fitting observed extinction curves in Local Group galaxies including the Magellanic Clouds and the Milky Way. The dust temperature and the emitted dust spectrum are calculated self consistently from the absorbed energy including the effects of temperature fluctuations in small grains. Dust self-absorption is also accounted for, allowing the treatment of high optical depths, by treating photons emitted by the dust as an additional heating source and adopting an iterative radiative transfer scheme. As an illustrative case, we apply the DIRTY radiative transfer code to starburst galaxies wherein the heating sources are derived from stellar evolutionary synthesis models. Within the context of the starburst model, we examine the dependence of the UV to FIR SED, dust temperatures, and dust masses predicted by DIRTY on variations of the input parameters.

We present 0.35 to 2.5 micron spectra of the south and northwest filaments in the reflection nebula NGC 7023. These spectra were used to test the theory of Seahra & Duley that carbon nanoparticles are responsible for Extended Red Emission (ERE). Our spectra fail to show their predicted second emission band at 1.0 micron even though both filaments exhibit strong emission in the familiar 0.7 micron ERE band. The northwest filament spectrum does show one, and possibly two, new dust emission features in the near-infrared. We clearly detect a strong emission band at 1.5 micron which we tentatively attribute to beta-FeSi_2 grains. We tentatively detect a weaker emission band at 1.15 micron which coincides with the location expected for transitions from the conduction band to mid-gap defect states of silicon nanoparticles. This is added evidence that silicon nanoparticles are responsible for ERE as they already can explain the observed behavior of the main visible ERE band.

We have studied the UV extinction properties along 30 Galactic sightlines using data from the International Ultraviolet Explorer (IUE) archive that have never been previously examined. These distant (d > 1 kpc) sightlines were selected to investigate the distribution and physical conditions of gas located in low density regions of the Galactic disk and halo. The average densities along these sightlines are extremely low. It is likely that they are dominated by the warm intercloud medium and have little contribution from the cold cloud medium. We find that a subsample of these sightlines has extinction curves with weak bumps and very steep far-UV extinction reminiscent of the Magellanic clouds. These sightlines all lie in the region bounded by 325 < l < 0 and -5 > b > -11. The gas along these sightlines shows forbidden velocities which may indicate that the dust has been subject to shocks. This type of low density sightline may mimic the environments found in the Magellanic Clouds. Large values of N(Ca II)/N(Na I) indicating low depletion are associated with steep far-UV extinction. A possible correlation exists between decreasing bump strength and increasing far-UV steepness for extinction curves in the Galaxy and the Magellanic Clouds.

The presence of dust in starburst galaxies complicates the study of their stellar populations as the dust's effects are similar to those associated with changes in the galaxies' stellar age and metallicity. This degeneracy can be overcome for starburst galaxies if UV/optical/near-infrared observations are combined with far-infrared observations. We present the calibration of the flux ratio method for calculating the dust attenuation at a particular wavelength, Att(\lambda), based on the measurement of F(IR)/F(\lambda) flux ratio. Our calibration is based on spectral energy distributions (SEDs) from the PEGASE stellar evolutionary synthesis model and the effects of dust (absorption and scattering) as calculated from our Monte Carlo radiative transfer model. We tested the attenuations predicted from this method for the Balmer emission lines of a sample starburst galaxies against those calculated using radio observations and found good agreement. The UV attenuation curves for a handful of starburst galaxies were calculated using the flux ratio method, and they compare favorably with past work. The relationship between Att(\lambda) and F(IR)/F(\lambda) is almost completely independent of the assumed dust properties (grain type, distribution, and clumpiness). For the UV, the relationship is also independent of the assumed stellar properties (age, metallicity, etc) accept for the case of very old burst populations. However at longer wavelengths, the relationship is dependent on the assumed stellar properties.

We present and discuss the results of new multiple-scattering radiative transfer calculations for three representative types of galactic environments, filled with either homogeneous or two-phase clumpy dust distributions. Extinction and scattering properties for two types of interstellar dust, similar to those found in the average diffuse medium of the Milky Way Galaxy (MW) and the Bar of the Small Magellanic Cloud (SMC), are considered. The wavelength coverage extends from 1000 A to 30,000 A, with particular emphasis on the rest-frame UV. This makes these models especially applicable to starburst galaxies and Lyman-break galaxy samples. The examination of the models concentrates on the study of UV/visual/near-IR reddening effects, the wavelength dependence of attenuation, and on the changes that arise from the transition from homogeneous to clumpy dust distributions in different star/dust geometries. Embedded dust, especially when clumpy, leads to saturation at fairly low reddening values with correspondingly gray attenuation functions. This makes the assessment of the attenuation of the far-UV flux from starburst galaxies difficult, if only UV/visual/near-IR data are available. Existing data for UV-selected starburst galaxies indicate a range of UV attenuation factors of 0-150. Our models reproduce the "Calzetti Attenuation Law", provided one adopts SMC-type dust and a clumpy shell-type dust distribution surrounding the starbursts. The average far-UV attenuation factor for the Calzetti sample is found to be 7.4. The only relatively reliable measure for the UV attenuation factor for single galaxies was found in the ratio of the integrated far-IR flux to the far-UV flux, measured near 1600 A, requiring the measurement of the entire spectral energy distribution of galaxies.

We present new images of the reflection nebula surrounding the R Coronae Borealis Star, UW Cen. This nebula, first detected in 1990, has changed its appearance significantly. At the estimated distance of UW Cen, this nebula is approximately 0.6 ly in radius so the nebula cannot have physically altered in only 8 years. Instead, the morphology of the nebula appears to change as different parts are illuminated by light from the central star modulated by shifting thick dust clouds near its surface. These dust clouds form and dissipate at irregular intervals causing the well-known declines in the R Coronae Borealis (RCB) stars. In this way, the central star acts like a lighthouse shining through holes in the dust clouds and lighting up different portions of the nebula. The existence of this nebula provides clues to the evolutionary history of RCB stars possibly linking them to the Planetary Nebulae and the final helium shell flash stars.

We have thoroughly characterized the ultraviolet to near-infrared (0.15 - 2.2 micron) spectral energy distribution (SED) of the central parsec of the M33 nucleus through new infrared photometry and optical/near-infrared spectroscopy, combined with ultraviolet/optical observations from the literature and the HST archive. The SED shows evidence for a significant level of attenuation, which we model through a Monte Carlo radiative transfer code as a shell of clumpy Milky Way-type dust (tau_V ~ 2 +/- 1). The discovery of Milky Way-type dust (with a strong 2175 A bump) internal to the M33 nucleus is different from previous work which has found SMC-like dust (no bump) near starburst regions. The amount by which dust can be processed may be related to the mass and age of the starburst as well as the extent to which the dust can shield itself. Our starburst models include the effects of this dust and can fit the SED if the nucleus was the site of a moderate (~10^8 L_sun at 10 Myrs) episode of coeval star formation about 70 Myrs ago. This result is quite different from previous studies which resorted to multiple stellar populations (between 2 and 7) attenuated by either no or very little internal dust. The M33 nuclear starburst is remarkably similar to an older version (70 Myr versus 10 Myr) of the ultra-compact starburst in the center of the Milky Way.

We have reanalyzed the Large Magellanic Cloud's (LMC) ultraviolet (UV) extinction using data from the IUE final archive. Our new analysis takes advantage of the improved signal--to--noise of the IUE NEWSIPS reduction, the exclusion of stars with very low reddening, the careful selection of well matched comparison stars, and an analysis of the effects of Galactic foreground dust. Differences between the average extinction curves of the 30 Dor region and the rest of the LMC are reduced compared to previous studies. We find that there is a group of stars with very weak 2175 Ang. bumps that lie in or near the region occupied by the supergiant shell, LMC 2, on the southeast side of 30 Dor. The average extinction curves inside and outside LMC 2 show a very significant difference in 2175 Ang. bump strength, but their far--UV extinctions are similar. While it is unclear whether or not the extinction outside the LMC 2 region can be fit with the relation of Cardelli, Clayton and Mathis (CCM), sightlines near LMC 2 cannot be fit with CCM due to their weak 2175 Ang. bumps. While the extinction properties seen in the LMC lie within the range of properties seen in the Galaxy, the correlations of UV extinction properties with environment seen in the Galaxy do not appear to hold in the LMC.

We have reviewed the characteristics of the extended red emission (ERE) as observed in many dusty astronomical environments, in particular, the diffuse interstellar medium of the Galaxy. The spectral nature and the photon conversion efficiency of the ERE identify the underlying process as highly efficient photoluminescence by an abundant component of interstellar dust. We have compared the photoluminescence properties of a variety of carbon- and silicon-based materials proposed as sources for the ERE with the observationally established constraints. We found that silicon nanoparticles provide the best match to the spectrum and the efficiency requirement of the ERE. If present in interstellar space with an abundance sufficient to explain the intensity of the ERE, silicon nanoparticles will also contribute to the interstellar 9.7 micron Si-O stretch feature in absorption, to the near- and mid-IR nonequilibrium thermal background radiation, and to the continuum extinction in the near- and far-UV. About 36% of the interstellar silicon depleted into the dust phase would be needed in the form of silicon nanoparticles, amounting to less than 5% of the interstellar dust mass. We propose that silicon nanoparticles form through the nucleation of SiO in oxygen-rich stellar mass outflows and that they represent an important small-grain component of the interstellar dust spectrum.

We have used spectra taken between 1992 and 1997 to derive the spectroscopic orbit of the eclipsing double-lined spectroscopic binary HD 197770. This binary has a period of 99.69 +/- 0.02 days and K amplitudes of 31.2 +/- 0.8 and 47.1 +/- 0.4 km s^{-1} for components A & B, respectively. The msin^{3}i values for A & B are 2.9 and 1.9, respectively, and are close to the actual masses due to the eclipsing nature of this binary. Both components of HD 197770 have spectral types near B2 III. This means both components are undermassive by about a factor of five and, thus, evolved stars. Additional evidence of the evolved nature of HD 197770 is found in 25, 60, and 100 micron IRAS images of HD 197770. These images show 2 apparent shells centered on HD 197770; a bright 60 micron shell with a 14' diameter and a larger (1.2 degeree diameter) bubble-like feature. At least one of the components of HD 197770 is likely to be a post-AGB star.

The recent discovery that the UV dust extinction in starburst galaxies is similar to that found in the Small Magellanic Cloud (SMC) motivated us to re-investigate the ultraviolet (UV) extinction found in the SMC. We have been able to improve significantly on previous studies by carefully choosing pairs of well matched reddened and unreddened stars. In addition, we benefited from the improved S/N of the NEWSIPS IUE data and the larger sample of SMC stars now available. Searching the IUE Final Archive, we found only four suitable early-type stars that were significantly reddened and had well matched comparison stars. The extinction for three of these stars is remarkably similar. The curves are roughly linear with 1/lambda and have no measurable 2175 A bump. The fourth star has an extinction curve with a significant 2175 A bump and weaker far-UV extinction. The dust along all four sightlines is thought to be local to the SMC. There is no significant Galactic foreground component. The first three stars lie in the SMC Bar and the line-of-sight for each of them passes through regions of recent star formation. The fourth star belongs to the SMC Wing and its line-of-sight passes though a much more quiescent region. Thus, the behavior of the dust extinction in the SMC supports a dependence of dust properties on star formation activity. However, other environmental factors (such as galactic metallicity) must also be important. Dust in the 30 Dor region of the LMC, where much more active star formation is present, does not share the extreme extinction properties seen in SMC dust.

Extended Red Emission (ERE) has been detected in many dusty astrophysical objects and this raises the question: Is ERE present only in discrete objects or is it an observational feature of all dust, i.e. present in the diffuse interstellar medium? In order to answer this question, we determined the blue and red intensities of the radiation from the diffuse interstellar medium (ISM) and examined the red intensity for the presence of an excess above that expected for scattered light. The diffuse ISM blue and red intensities were obtained by subtracting the integrated star and galaxy intensities from the blue and red measurements made by the Imaging Photopolarimeter (IPP) aboard the Pioneer 10 and 11 spacecraft. The color of the diffuse ISM was found to be redder than the Pioneer intensities, contrary to that expected for scattered light (Diffuse Galactic Light [DGL]). The red DGL was calculated using the blue diffuse ISM intensities and the approximately invariant color of the DGL calculated with the DGL model. Subtracting the calculated red DGL from the red diffuse ISM intensities resulted in the detection of an excess red intensity. This represents the likely detection of ERE in the diffuse ISM since H-alpha emission cannot account for the strength of this excess and the only other known emission process applicable to the diffuse ISM is ERE. Thus, ERE appears to be a general characteristic of dust. From the correlation between N_HI and ERE intensity the ERE photon conversion efficiency was estimated at 10 +/- 3%.

To investigate the nature of starbursts' dust, we constructed a model of the stars and dust in starburst galaxies and applied it to 30 observed starburst spectral energy distributions (SEDs). The starburst model was constructed by combining two stellar evolutionary synthesis models with a model describing the radiative transfer of stellar photons through dust. The stellar evolutionary synthesis models were used to compute the dust-free SEDs for stellar populations with ages between 1x10^6 and 15x10^9 years. Using a Monte Carlo radiative transfer model, the effects of dust were computed for average Milky Way (MW) and Small Magellanic Cloud (SMC) dust, two different star/dust geometries, and locally homogeneous or clumpy dust. Using color-color plots, the starburst model was used to interpret the behavior of 30 starbursts with aperture-matched UV and optical SEDs (and IR for 19 of the 30) from previous studies. From the color-color plots, it was evident that the dust in starbursts has an extinction curve lacking a 2175 A bump, like the SMC curve, and a steep far-UV rise, intermediate between the MW and SMC curves. The star/dust geometry which is able to explain the distribution of the 30 starbursts in various color-color plots has an inner dust-free sphere of stars surrounded by an outer star-free shell of clumpy dust. When combined with other work from the literature on the Orion region and the 30 Dor region of the Large Magellanic Cloud, this work implies a trend in dust properties with star formation intensity.

Using the 1-m telescope at Ritter Observatory, we took 36 observations of zeta^2 Coronae Borealis with a fiber-fed echelle spectrograph. From these observations, zeta^2 CrB was found to be a triple system and a new spectroscopic orbit was calculated. This orbit has two periods, a 1.72357 day period for the inner binary composed of zeta^2 CrB A & B and a 251 day period for the outer binary composed of zeta^2 CrB AB & C. The inner binary is a double-lined spectroscopic binary composed of two B7 V stars. The inner binary's center of mass (zeta^2 CrB AB) describes a long-period single-lined variation identified with the outer binary. The inner binary period is significantly shorter than the 12.5842 day period previously calculated by Abhyankar & Sarma (1966). The inner binary possesses an essentially circular orbit (e = 0.01) while the outer binary has an eccentric orbit (e = 0.48). From the widths of their Si II 6371 A lines, the vsini's were calculated to be 46 +/- 7 km s^-1 for zeta^2 CrB A and 7.5 +/- 2 km s^-1 for zeta^2 CrB B. As zeta^2 CrB A & B have similar masses, their different rotational velocities make this system a sensitive test of synchronization theories.