NGC 4151 is among the most well-studied Seyfert galaxies that does not suffer from strong obscuration along the observer's line-of-sight. This allows to probe the central active galactic nucleus (AGN) engine with photometry, spectroscopy, reverberation mapping or interferometry. Yet, the broadband polarization from NGC 4151 has been poorly examined in the past despite the fact that polarimetry gives us a much cleaner view of the AGN physics than photometry or spectroscopy alone. In this paper, we compile the 0.15 -- 89.0 $\mu$m total and polarized fluxes of NGC 4151 from archival and new data in order to examine the physical processes at work in the heart of this AGN. We demonstrate that, from the optical to the near-infrared (IR) band, the polarized spectrum of NGC 4151 shows a much bluer power-law spectral index than that of the total flux, corroborating the presence of an optically thick, locally heated accretion flow, at least in its near-IR emitting radii. Specific signatures from the atmosphere of the accretion structure are tentatively found at the shortest ultraviolet (UV) wavelengths, before the onset of absorption opacity. Otherwise, dust scattering appears to be the dominant contributor from the near-UV to near-IR polarized spectrum, superimposed onto a weaker electron component. We also identify a change in the polarization processes from the near-IR to the mid-IR, most likely associated with the transition from Mie scattering to dichroic absorption from aligned dust grains in the dusty torus or narrow-line region. Finally, we present and dicuss the very first far-infrared polarization measurement of NGC 4151 at 89 $\mu$m.

We have studied the small scale distribution of atomic hydrogen (HI) using 21-cm absorption spectra against multiple-component background radio continuum sources from the 21-SPONGE survey and the Millennium Arecibo Absorption Line Survey. We have found $>5\sigma$ optical depth variations at a level of $\sim0.03-0.5$ between 13 out of 14 adjacent sightlines separated by a few arcseconds to a few arcminutes, suggesting the presence of neutral structures on spatial scales from a few to thousands of AU (which we refer to as tiny scale atomic structure, TSAS). The optical depth variations are strongest in directions where the HI column density and the fraction of HI in the cold neutral medium (CNM) are highest, which tend to be at low Galactic latitudes. By measuring changes in the properties of Gaussian components fitted to the absorption spectra, we find that changes in both the peak optical depth and the linewidth of TSAS absorption features contribute to the observed optical depth variations, while changes in the central velocity do not appear to strongly impact the observed variations. Both thermal and turbulent motions contribute appreciably to the linewidths, but the turbulence does not appear strong enough to confine overpressured TSAS. In a majority of cases, the TSAS column densities are sufficiently high that these structures can radiatively cool fast enough to maintain thermal equilibrium with their surroundings, even if they are overpressured. We also find that a majority of TSAS is associated with the CNM. For TSAS in the direction of the Taurus molecular cloud and the local Leo cold cloud, we estimate densities over an order of magnitude higher than typical CNM densities.

We developed methods for mapping spatial variations of the spatial power spectrum (SPS) and structure function (SF) slopes, with a goal of connecting neutral hydrogen (HI) statistical properties with the turbulent drivers. The new methods were applied on the HI observations of the Small and Large Magellanic Clouds (SMC and LMC). In the case of the SMC, we find highly uniform turbulent properties of HI, with no evidence for local enhancements of turbulence due to stellar feedback. Such properties could be caused by a significant turbulent driving on large-scales. Alternatively, a significant line-of-sight depth of the SMC could be masking out localized regions with a steeper SPS slope caused by stellar feedback. In contrast to the SMC, the LMC HI shows a large diversity in terms of its turbulent properties. Across most of the LMC, the small-scale SPS slope is steeper than the large-scale slope due to the presence of the HI disk. On small spatial scales, we find several areas of localized steepening of the SPS slope around major HII regions, with the 30 Doradus region being the most prominent. This is in agreement with predictions from numerical simulations which suggest steepening of the SPS slope due to stellar feedback eroding and destroying interstellar clouds. We also find localized steepening of the large-scale SPS slope in the outskirts of the LMC. This is likely caused by the flaring of the HI disk, or alternatively ram-pressure stripping of the LMC disk due to the interactions with the surrounding halo gas.

We study a very young star-forming region in the outer Galaxy that is the most concentrated source of outflows in the Spitzer Space Telescope GLIMPSE360 survey. This region, dubbed CMa-l224, is located in the Canis Major OB1 association. CMa-l224 is relatively faint in the mid-infrared, but it shines brightly at the far-infrared wavelengths as revealed by the Herschel Space Observatory data from the Hi-GAL survey. Using the 3.6 and 4.5 $\mu$m data from the Spitzer/GLIMPSE360 survey, combined with the JHK$_s$ 2MASS and the 70-500 $\mu$m Herschel/Hi-GAL data, we develop a young stellar object (YSO) selection criteria based on color-color cuts and fitting of the YSO candidates' spectral energy distributions with YSO 2D radiative transfer models. We identify 293 YSO candidates and estimate physical parameters for 210 sources well-fit with YSO models. We select an additional 47 sources with GLIMPSE360-only photometry as `possible YSO candidates'. The vast majority of these sources are associated with high H$_2$ column density regions and are good targets for follow-up studies. The distribution of YSO candidates at different evolutionary stages with respect to Herschel filaments supports the idea that stars are formed in the filaments and become more dispersed with time. Both the supernova-induced and spontaneous star formation scenarios are plausible in the environmental context of CMa-l224. However, our results indicate that a spontaneous gravitational collapse of filaments is a more likely scenario. The methods developed for CMa-l224 can be used for larger regions in the Galactic plane where the same set of photometry is available.

We present 21-cm Spectral Line Observations of Neutral Gas with the VLA (21-SPONGE), a Karl G. Jansky Very Large Array (VLA) large project (~600 hours) for measuring the physical properties of Galactic neutral hydrogen (HI). 21-SPONGE is distinguished among previous Galactic HI studies as a result of: (1) exceptional optical depth sensitivity ($\sigma_{\tau} < 10^{-3}$ per $0.42\rm\,km\,s^{-1}$ channels over 57 lines of sight); (2) matching 21 cm emission spectra with highest-possible angular resolution (~4') from the Arecibo Observatory; (3) detailed comparisons with numerical simulations for assessing observational biases. We autonomously decompose 21 cm spectra and derive the physical properties (i.e., spin temperature, $T_s$, column density) of the cold neutral medium (CNM; $T_s<250\rm\,K$), thermally unstable medium (UNM; $250< T_s < 1000\rm\,K$) and warm neutral medium (WNM; $T_s > 1000\rm\,K$) simultaneously. We detect 50% of the total HI mass in absorption, the majority of which is CNM (56 +/- 10%, corresponding to 28% of the total HI mass). Although CNM is detected ubiquitously, the CNM fraction along most lines of sight is <50%. We find that 20% of the total HI mass is thermally unstable (41 +/- 10% of HI detected in absorption), with no significant variation with Galactic environment. Finally, although the WNM comprises 52% of the total HI mass, we detect little evidence for WNM absorption with $1000<T_s<4000\rm\,K$. Following spectral modeling, we detect a stacked residual absorption feature corresponding to WNM with $T_s\sim10^4\rm\,K$. We conclude that excitation in excess of collisions likely produces significantly higher WNM $T_s$ than predicted by steady-state models.

We investigate spatial variations of turbulent properties in the Small Magellanic Cloud (SMC) by using neutral hydrogen HI observations. With the goal of testing the importance of stellar feedback on HI turbulence, we define central and outer SMC regions based on the star formation rate (SFR) surface density, as well as the HI integrated intensity. We use the structure function and the Velocity Channel Analysis (VCA) to calculate the power-law index (gamma) for both underlying density and velocity fields in these regions. In all cases, our results show essentially no difference in gamma between the central and outer regions. This suggests that HI turbulent properties are surprisingly homogeneous across the SMC when probed at a resolution of 30 pc. Contrary to recent suggestions from numerical simulations, we do not find a significant change in gamma due to stellar feedback as traced by the SFR surface density. This could be due to the stellar feedback being widespread over the whole of the SMC, but more likely due to a large-scale gravitational driving of turbulence. We show that the lack of difference between central and outer SMC regions can not be explained by the high optical depth HI.

We analyze synthetic neutral hydrogen (HI) absorption and emission spectral lines from a high- resolution, three-dimensional hydrodynamical simulation to quantify how well observational methods recover the physical properties of interstellar gas. We present a new method for uniformly decomposing HI spectral lines and estimating the properties of associated gas using the Autonomous Gaussian Decomposition (AGD) algorithm. We find that Hi spectral lines recover physical structures in the simulation with excellent completeness at high Galactic latitude, and this completeness declines with decreasing latitude due to strong velocity-blending of spectral lines. The temperature and column density inferred from our decomposition and radiative transfer method agree with the simulated values within a factor of < 2 for the majority of gas structures. We next compare synthetic spectra with observations from the 21-SPONGE survey at the Karl G. Jansky Very Large Array using AGD. We find more components per line of sight in 21-SPONGE than in synthetic spectra, which reflects insufficient simulated gas scale heights and the limitations of local box simulations. In addition, we find a significant population of low-optical depth, broad absorption components in the synthetic data which are not seen in 21-SPONGE. This population is not obvious in integrated or per-channel diagnostics, and reflects the benefit of studying velocity-resolved components. The discrepant components correspond to the highest spin temperatures (1000 < Ts < 4000 K), which are not seen in 21-SPONGE despite sufficient observational sensitivity. We demonstrate that our analysis method is a powerful tool for diagnosing neutral ISM conditions, and future work is needed to improve observational statistics and implementation of simulated physics.

We analyze 45 spectropolarimetric observations of the eclipsing, interacting binary star V356 Sgr, obtained over a period of 21 years, to characterize the geometry of the system's circumstellar material. After removing interstellar polarization from these data, we find the system exhibits a large intrinsic polarization signature arising from electron scattering. In addition, the lack of repeatable eclipses in the polarization phase curves indicates the presence of a substantial pool of scatterers not occulted by either star. We suggest that these scatterers form either a circumbinary disk coplanar with the gainer's accretion disk or an elongated structure perpendicular to the orbital plane of V356 Sgr, possibly formed by bipolar outflows. We also observe small-scale, cycle-to-cycle variations in the magnitude of intrinsic polarization at individual phases, which we interpret as evidence of variability in the amount of scattering material present within and around the system. This may indicate a mass transfer or mass loss rate that varies on the time-scale of the system's orbital period. Finally, we compare the basic polarimetric properties of V356 Sgr with those of the well studied beta Lyr system; the significant differences observed between the two systems suggests diversity in the basic circumstellar geometry of Roche-lobe overflow systems

Using high-resolution data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey, we show that linear structure in Galactic neutral hydrogen (HI) correlates with the magnetic field orientation implied by Planck 353 GHz polarized dust emission. The structure of the neutral interstellar medium is more tightly coupled to the magnetic field than previously known. At high Galactic latitudes, where the Planck data are noise-dominated, the HI data provide an independent constraint on the Galactic magnetic field orientation, and hence the local dust polarization angle. We detect strong cross-correlations between template maps constructed from estimates of dust intensity combined with either HI-derived angles, starlight polarization angles, or Planck 353 GHz angles. The HI data thus provide a new tool in the search for inflationary gravitational wave B-mode polarization in the cosmic microwave background, which is currently limited by dust foreground contamination.

We present the first detection of HCO$^+$ absorption in the Magellanic System. Using the Australia Telescope Compact Array (ATCA), we observed 9 extragalactic radio continuum sources behind the Magellanic System and detected HCO$^+$ absorption towards one source located behind the leading edge of the Magellanic Bridge. The detection is located at LSR velocity of $v=214.0 \pm 0.4\rm\,km\,s^{-1}$, with a full width at half maximum of $\Delta v=4.5\pm 1.0\rm\,km\,s^{-1}$ and optical depth of $\tau(\rm HCO^+)=0.10\pm 0.02$. Although there is abundant neutral hydrogen (HI) surrounding the sightline in position-velocity space, at the exact location of the absorber the HI column density is low, $<10^{20}\rm\,cm^{-2}$, and there is little evidence for dust or CO emission from Planck observations. While the origin and survival of molecules in such a diffuse environment remains unclear, dynamical events such as HI flows and cloud collisions in this interacting system likely play an important role.

We present the first results from the SAGE-Var program, a follow on to the Spitzer legacy program Surveying the Agents of Galaxy Evolution (SAGE; Meixner, et al. 2006). We obtained 4 epochs of photometry at 3.6 & 4.5 microns covering the bar of the Large Magellanic Cloud (LMC) and the central region of the Small Magellanic Cloud (SMC) in order to probe the variability of extremely red sources missed by variability surveys conducted at shorter wavelengths, and to provide additional epochs of observation for known variables. Our 6 total epochs of observations allow us to probe infrared variability on 15 different timescales ranging from ~20 days to ~5 years. Out of a full catalog of 1,717,554 (LMC) and 457,760 (SMC) objects, we find 10 (LMC) and 6 (SMC) large amplitude AGB variables without optically measured variability owing to circumstellar dust obscuration. The catalog also contains multiple observations of known AGB variables, type I and II Cepheids, eclipsing variables, R CrB stars and young stellar objects which will be discussed in following papers. Here we present infrared Period-Luminosity (PL) relations for classical Cepheids in the Magellanic Clouds, as well as improved PL relationships for AGB stars pulsating in the fundamental mode using mean magnitudes constructed from 6 epochs of observations.

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.

We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $\sigma_{\tau}=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

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.

We report the results of infrared photometry on 39 extended red objects (EROs) in the Carina Nebula, observed with the Spitzer Space Telescope. Most EROs are identified by bright, extended 8.0 um emission, which ranges from 10'' to 40'' in size, but our sample also includes 4 EROs identified by extended 24 um emission. Of particular interest are nine EROs associated with late O or early B-type stars and characterized by arc-shaped morphology, suggesting dusty, stellar-wind bow shocks. These objects are preferentially oriented towards the central regions of the Carina Nebula, suggesting that these bow shocks are generally produced by the interactions of OB winds with the bulk expansion of the H II region rather than high proper motion. We identify preferred regions of mid-infrared color space occupied by our bow shock candidates, which also contain bow shock candidates in M17 and RCW 49 but are well-separated from polycyclic aromatic hydrocarbon emission or circumstellar discs. Color cuts identify an additional 12 marginally-resolved bow shock candidates, 10 of which are also associated with known late O or early B stars. H II region expansion velocities derived from bow shock candidate standoff distances are ~10 km/s, assuming typical H II region gas densities, comparable to expansion velocities derived from bow shocks in M17 and RCW 49. One candidate bow shock provides direct evidence of physical interaction between the massive stellar winds originating in the Trumpler 15 and Trumpler 14 clusters, supporting the conclusion that both clusters are at similar heliocentric distances.

We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21cm absorption spectra from the 21cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the HI line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes.

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.

We use the Karl G. Jansky Very Large Array (VLA) to conduct a high-sensitivity survey of neutral hydrogen (HI) absorption in the Milky Way. In combination with corresponding HI emission spectra obtained mostly with the Arecibo Observatory, we detect a widespread warm neutral medium (WNM) component with excitation temperature <Ts>= 7200 (+1800,-1200) K (68% confidence). This temperature lies above theoretical predictions based on collisional excitation alone, implying that Ly-{\alpha} scattering, the most probable additional source of excitation, is more important in the interstellar medium (ISM) than previously assumed. Our results demonstrate that HI absorption can be used to constrain the Ly-{\alpha} radiation field, a critical quantity for studying the energy balance in the ISM and intergalactic medium yet notoriously difficult to model because of its complicated radiative transfer, in and around galaxies nearby and at high redshift.

We characterize the completeness of point source lists from Spitzer Space Telescope surveys in the four Infrared Array Camera (IRAC) bandpasses, emphasizing the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) programs (GLIMPSE I, II, 3D, 360; Deep GLIMPSE) and their resulting point source Catalogs and Archives. The analysis separately addresses effects of incompleteness resulting from high diffuse background emission and incompleteness resulting from point source confusion (i.e., crowding). An artificial star addition and extraction analysis demonstrates that completeness is strongly dependent on local background brightness and structure, with high-surface-brightness regions suffering up to five magnitudes of reduced sensitivity to point sources. This effect is most pronounced at the IRAC 5.8 and 8.0 microns bands where UV-excited PAH emission produces bright, complex structures (photodissociation regions; PDRs). With regard to diffuse background effects, we provide the completeness as a function of stellar magnitude and diffuse background level in graphical and tabular formats. These data are suitable for estimating completeness in the low-source-density limit in any of the four IRAC bands in GLIMPSE Catalogs and Archives and some other Spitzer IRAC programs that employ similar observational strategies and are processed by the GLIMPSE pipeline. Point source incompleteness is primarily a consequence of structure in the diffuse background emission rather than photon noise. With regard to source confusion in the high-source-density regions of the Galactic Plane, we provide figures illustrating the 90% completeness levels as a function of point source density at each band. (Slightly abridged)

We present a self-consistent three-dimensional Monte-Carlo radiative transfer model of the stellar and dust emission in the Milky-Way, and have computed synthetic observations of the 3.6 to 100 microns emission in the Galactic mid-plane. In order to compare the model to observations, we use the GLIMPSE, MIPSGAL, and IRAS surveys to construct total emission spectra, as well as longitude and latitude profiles for the emission. The distribution of stars and dust is taken from the SKY model, and the dust emissivities includes an approximation of the emission from polycyclic aromatic hydrocarbons in addition to thermal emission. The model emission is in broad agreement with the observations, but a few modifications are needed to obtain a good fit. Firstly, by adjusting the model to include two major and two minor spiral arms rather than four equal spiral arms, the fit to the longitude profiles for |l|>30 degrees can be improved. Secondly, introducing a deficit in the dust distribution in the inner Galaxy results in a better fit to the shape of the IRAS longitude profiles at 60 and 100 microns. With these modifications, the model fits the observed profiles well, although it systematically under-estimates the 5.8 and 8.0 microns fluxes. One way to resolve this discrepancy is to increase the abundance of PAH molecules by 50% compared to the original model, although we note that changes to the dust distribution or radiation field may provide alternative solutions. Finally, we use the model to quantify which stellar populations contribute the most to the heating of different dust types, and which stellar populations and dust types contribute the most to the emission at different wavelengths.

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.

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.

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

We present new Spitzer/IRS spectra of two hot R Coronae Borealis (RCB) stars, one in the Galaxy,V348 Sgr, and one lying in the LMC, HV 2671. These two objects may constitute a link between the RCB stars and the late Wolf-Rayet ([WCL]) class of central stars of planetary nebula (CSPNe) such as CPD -56 8032 that has little or no hydrogen in their atmospheres. HV 2671 and V348 Sgr are members of a rare subclass that has significantly higher effective temperatures than most RCB stars, but sharing the traits of hydrogen deficiency and dust formation that define the cooler RCB stars. The [WC] CSPNe star, CPD -56 8032, displays evidence for dual-dust chemistry showing both PAHs and crystalline silicates in its mid-IR spectrum. HV 2671 shows strong PAH emission but shows no sign of having crystalline silicates. The spectrum of V348 Sgr is very different from those of CPD -56 8032 and HV 2671. The PAH emission seen strongly in the other two stars is not present. Instead, the spectrum is dominated by a broad emission centered at about 8.2 micron. The mid-IR spectrum of CPD -56 8032 shows emission features that may be associated with C60. The other two stars do not show evidence for C60. HV 2671 has also been detected by Herschel/PACS and SPIRE. V348 Sgr and CPD -56 8032 have been detected by AKARI/FIS. These data were combined with Spitzer, IRAS, 2MASS and other photometry to produce their spectral energy distributions from the visible to the far-IR. Monte Carlo radiative transfer modeling was used to study the circumstellar dust around these stars. HV 2671 and CPD -56 8032 require both a flared inner disk with warm dust and an extended diffuse envelope with cold dust to to fit their SEDs. The SED of V348 Sgr can be fit with a much smaller disk and envelope.

We present a catalog of 1439 young stellar objects (YSOs) spanning the 1.42 deg^2 field surveyed by the Chandra Carina Complex Project (CCCP), which includes the major ionizing clusters and the most active sites of ongoing star formation within the Great Nebula in Carina. Candidate YSOs were identified via infrared (IR) excess emission from dusty circumstellar disks and envelopes, using data from the Spitzer Space Telescope Vela--Carina survey and the Two-Micron All Sky Survey. We model the 1--24 /mu m IR spectral energy distributions of the YSOs to constrain physical properties. Our Pan-Carina YSO Catalog (PCYC) is dominated by intermediate-mass (2 Msun < m < 10 Msun) objects with disks, including Herbig Ae/Be stars and their less evolved progenitors. The PCYC provides a valuable complementary dataset to the CCCP X-ray source catalogs, identifying 1029 YSOs in Carina with no X-ray detection. We also catalog 410 YSOs with X-ray counterparts, including 62 candidate protostars. Candidate protostars with X-ray detections tend to be more evolved than those without. In most cases, X-ray emission apparently originating from intermediate-mass, disk-dominated YSOs is consistent with the presence of low-mass companions, but we also find that X-ray emission correlates with cooler stellar photospheres and higher disk masses. We suggest that intermediate-mass YSOs produce X-rays during their early pre-main sequence evolution, perhaps driven by magnetic dynamo activity during the convective atmosphere phase, but this emission dies off as the stars approach the main sequence. Extrapolating over the stellar initial mass function scaled to the PCYC population, we predict a total population of >2x10^4 YSOs and a present-day star formation rate (SFR) of >0.008 Msun/yr. The global SFR in the Carina Nebula, averaged over the past ~5 Myr, has been approximately constant.

We report the results of a new survey of massive, OB stars throughout the Carina Nebula using the X-ray point source catalog provided by the Chandra Carina Complex Project (CCCP) in conjunction with infrared (IR) photometry from the Two Micron All-Sky Survey and the Spitzer Space Telescope Vela--Carina survey. Mid-IR photometry is relatively unaffected by extinction, hence it provides strong constraints on the luminosities of OB stars, assuming that their association with the Carina Nebula, and hence their distance, is confirmed. We fit model stellar atmospheres to the optical (UBV) and IR spectral energy distributions (SEDs) of 182 OB stars with known spectral types and measure the bolometric luminosity and extinction for each star. We find that the extinction law measured toward the OB stars has two components: Av=1--1.5 mag produced by foreground dust with a ratio of total-to-selective absorption Rv=3.1 plus a contribution from local dust with Rv>4.0 in the Carina molecular clouds that increases as Av increases. Using X-ray emission as a strong indicator of association with Carina, we identify 94 candidate OB stars with Lbol\geq10^4 Lsun by fitting their IR SEDs. If the candidate OB stars are eventually confirmed by follow-up spectroscopic observations, the number of cataloged OB stars in the Carina Nebula will increase by ~50%. Correcting for incompleteness due to OB stars falling below the Lbol cutoff or the CCCP detection limit, these results potentially double the size of the young massive stellar population.

We present a catalog of ~14,000 X-ray sources observed by the ACIS instrument on the Chandra X-ray Observatory within a 1.42 square degree survey of the Great Nebula in Carina, known as the Chandra Carina Complex Project (CCCP). This study appears in a Special Issue of the ApJS devoted to the CCCP. Here, we describe the data reduction and analysis procedures performed on the X-ray observations, including calibration and cleaning of the X-ray event data, point source detection, and source extraction. The catalog appears to be complete across most of the field to an absorption-corrected total-band luminosity of ~10^{30.7} erg/s for a typical low-mass pre-main sequence star. Counterparts to the X-ray sources are identified in a variety of visual, near-infrared, and mid-infrared surveys. The X-ray and infrared source properties presented here form the basis of many CCCP studies of the young stellar populations in Carina.

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

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.

We report the first results of imaging the Carina Nebula with Spitzer/IRAC, providing a catalog of point sources and YSOs based on SED fits. We discuss several aspects of the extended emission, including dust pillars that result when a clumpy molecular cloud is shredded by massive star feedback. There are few "extended green objects" (EGOs) normally taken as signposts of outflow activity, and none of the HH jets detected optically are seen as EGOs. A population of "extended red objects" tends to be found around OB stars, some with clear bow-shocks. These are dusty shocks where stellar winds collide with flows off nearby clouds. Finally, the relative distributions of O stars and subclusters of YSOs as compared to dust pillars shows that while some YSOs are located within pillars, many more stars and YSOs reside just outside pillar heads. We suggest that pillars are transient phenomena, part of a continuous outwardly propagating wave of star formation driven by massive star feedback. As pillars are destroyed, they leave newly formed stars in their wake, which are then subsumed into the young OB association. Altogether, the current generation of YSOs shows no strong deviation from a normal IMF. The number of YSOs suggests a roughly constant star-formation rate over the past 3Myr, implying that star formation in pillars constitutes an important mechanism to construct unbound OB associations. Accelerated pillars may give birth to O-type stars that, after several Myr, could appear to have formed in isolation.

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 combine near-infrared (2MASS) and mid-infrared (Spitzer-IRAC) photometry to characterize the IR extinction law (1.2-8 microns) over nearly 150 degrees of contiguous Milky Way midplane longitude. The relative extinctions in 5 passbands across these wavelength and longitude ranges are derived by calculating color excess ratios for G and K giant red clump stars in contiguous midplane regions and deriving the wavelength dependence of extinction in each one. Strong, monotonic variations in the extinction law shape are found as a function of angle from the Galactic center, symmetric on either side of it. These longitudinal variations persist even when dense interstellar regions, known a priori to have a shallower extinction curve, are removed. The increasingly steep extinction curves towards the outer Galaxy indicate a steady decrease in the absolute-to-selective extinction ratio (R_V) and in the mean dust grain size at greater Galactocentric angles. We note an increasing strength of the 8 micron extinction inflection at high Galactocentric angles and, using theoretical dust models, show that this behavior is consistent with the trend in R_V. Along several lines of sight where the solution is most feasible, A_lambda/A_Ks as a function of Galactic radius is estimated and shown to have a Galactic radial dependence. Our analyses suggest that the observed relationship between extinction curve shape and Galactic longitude is due to an intrinsic dependence of the extinction law on Galactocentric radius.

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

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

M17 is one of the youngest and most massive nearby star-formation regions in the Galaxy. It features a bright H II region erupting as a blister from the side of a giant molecular cloud (GMC). Combining photometry from the Spitzer GLIMPSE survey with complementary infrared (IR) surveys, we identify candidate young stellar objects (YSOs) throughout a 1.5 deg x 1 deg field that includes the M17 complex. The long sightline through the Galaxy behind M17 creates significant contamination in our YSO sample from unassociated sources with similar IR colors. Removing contaminants, we produce a highly-reliable catalog of 96 candidate YSOs with a high probability of association with the M17 complex. We fit model spectral energy distributions to these sources and constrain their physical properties. Extrapolating the mass function of 62 intermediate-mass YSOs (M >3 Msun), we estimate that >1000 stars are in the process of forming in the extended outer regions of M17. From IR survey images from IRAS and GLIMPSE, we find that M17 lies on the rim of a large shell structure ~0.5 deg in diameter (~20 pc at 2.1 kpc). We present new maps of CO and 13CO (J=2-1) emission, which show that the shell is a coherent, kinematic structure associated with M17 at v = 19 km/s. The shell is an extended bubble outlining the photodissociation region of a faint, diffuse H II region several Myr old. We provide evidence that massive star formation has been triggered by the expansion of the bubble. The formation of the massive cluster ionizing the M17 H II region itself may have been similarly triggered. We conclude that the star formation history in the extended environment of M17 has been punctuated by successive waves of massive star formation propagating through a GMC complex.

We present an analysis of late-O/early-B-powered, parsec-sized bubbles and associated star-formation using 2MASS, GLIMPSE, MIPSGAL and MAGPIS surveys. Three bubbles were selected from the Churchwell et al. (2007) catalog. We confirm that the structure identified in Watson et al. (2008) holds in less energetic bubbles, i.e. a PDR, identified by 8 um emission due to PAHs surrounds hot dust, identified by 24 um emission and ionized gas, identified by 20 cm continuum. We estimate the dynamical age of two bubbles by comparing bubble sizes to numerical models of Hosokawa & Inutsuka (2006). We also identify and analyze candidate young stellar objects (YSOs) using SED fitting and identify sites of possible triggered star-formation. Lastly, we identify likely ionizing sources for two sources based on SED fitting.

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

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.

Using images from the Spitzer GLIMPSE Legacy survey, we have identified more than 300 extended 4.5 micron sources (abbreviated EGO, Extended Green Object, for the common coding of the [4.5] band as green in 3-color composite IRAC images). We present a catalog of these EGOs, including integrated flux density measurements at 3.6, 4.5, 5.8, 8.0, and 24 microns from the GLIMPSE and MIPSGAL surveys. The average angular separation between a source in our sample and the nearest IRAS point source is >1 arcminute. The majority of EGOs are associated with infrared dark clouds (IRDCs), and where high-resolution 6.7 GHz methanol maser surveys overlap the GLIMPSE coverage, EGOs and 6.7 GHz methanol masers are strongly correlated. Extended 4.5 micron emission is thought to trace shocked molecular gas in protostellar outflows; the association of EGOs with IRDCs and 6.7 GHz methanol masers suggests that extended 4.5 micron emission may pinpoint outflows specifically from massive protostars. The mid-infrared colors of EGOs lie in regions of color-color space occupied by young protostars still embedded in infalling envelopes.

We present a highly reliable flux-limited census of 18,949 point sources in the Galactic mid-plane that have intrinsically red mid-infrared colors. These sources were selected from the Spitzer Space Telescope GLIMPSE I and II surveys of 274 deg^2 of the Galactic mid-plane, and consist mostly of high- and intermediate- mass young stellar objects (YSOs) and asymptotic giant branch (AGB) stars. The selection criteria were carefully chosen to minimize the effects of position- dependent sensitivity, saturation, and confusion. The distribution of sources on the sky and their location in IRAC and MIPS 24 microns color-magnitude and color-color space are presented. Using this large sample, we find that YSOs and AGB stars can be mostly separated by simple color-magnitude selection criteria into approximately 50-70% of YSOs and 30-50% of AGB stars. Planetary nebulae and background galaxies together represent at most 2-3% of all the red sources. 1,004 red sources in the GLIMPSE II region, mostly AGB stars with high mass- loss rates, show significant (>0.3 mag) variability at 4.5 and/or 8.0 microns. With over 11,000 likely YSOs and over 7,000 likely AGB stars, this is to date the largest uniform census of AGB stars and high- and intermediate mass YSOs in the Milky-Way Galaxy.

We report the discovery of six infrared stellar-wind bowshocks in the Galactic massive star formation regions M17 and RCW49 from Spitzer GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire) images. The InfraRed Array Camera (IRAC) on the Spitzer Space Telescope clearly resolves the arc-shaped emission produced by the bowshocks. We combine Two Micron All-Sky Survey (2MASS), Spitzer, MSX, and IRAS observations to obtain the spectral energy distributions (SEDs) of the bowshocks and their individual driving stars. We use the stellar SEDs to estimate the spectral types of the three newly-identified O stars in RCW49 and one previously undiscovered O star in M17. One of the bowshocks in RCW49 reveals the presence of a large-scale flow of gas escaping the H II region at a few 10^2 km/s. Radiation-transfer modeling of the steep rise in the SED of this bowshock toward longer mid-infrared wavelengths indicates that the emission is coming principally from dust heated by the star driving the shock. The other 5 bowshocks occur where the stellar winds of O stars sweep up dust in the expanding H II regions.

We present an analysis of wind-blown, parsec-sized, mid-infrared bubbles and associated star-formation using GLIMPSE/IRAC, MIPSGAL/MIPS and MAGPIS/VLA surveys. Three bubbles from the Churchwell et al. (2006) catalog were selected. The relative distribution of the ionized gas (based on 20 cm emission), PAH emission (based on 8 um, 5.8 um and lack of 4.5 um emission) and hot dust (24 um emission) are compared. At the center of each bubble there is a region containing ionized gas and hot dust, surrounded by PAHs. We identify the likely source(s) of the stellar wind and ionizing flux producing each bubble based upon SED fitting to numerical hot stellar photosphere models. Candidate YSOs are also identified using SED fitting, including several sites of possible triggered star formation.

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 present a list of 552 sources with suspected variability, based on a comparison of mid-infrared photometry from the GLIMPSE I and MSX surveys, which were carried out nearly a decade apart. We were careful to address issues such as the difference in resolution and sensitivity between the two surveys, as well as the differences in the spectral responses of the instruments. We selected only sources where the IRAC 8.0 and MSX 8.28 micron fluxes differ by more than a factor of two, in order to minimize contamination from sources where the difference in fluxes at 8 micron is due to a strong 10 micron silicate feature. We present a subset of 40 sources for which additional evidence suggests variability, using 2MASS and MIPSGAL data. Based on a comparison with the variability flags in the IRAS and MSX Point-Source Catalogs we estimate that at least a quarter of the 552 sources, and at least half of the 40 sources are truly variable. In addition, we tentatively confirm the variability of one source using multi-epoch IRAS LRS spectra. We suggest that most of the sources in our list are likely to be Asymptotic Giant Branch stars.

We present new Spitzer photometry of the Eagle Nebula (M16, containing the optical cluster NGC 6611) combined with near-infrared photometry from 2MASS. We use dust radiative transfer models, mid-infrared and near-infrared color-color analysis, and mid-infrared spectral indices to analyze point source spectral energy distributions, select candidate young stellar objects (YSOs), and constrain their mass and evolutionary state. Comparison of the different protostellar selection methods shows that mid-infrared methods are consistent, but as has been known for some time, near-infrared-only analysis misses some young objects. We reveal more than 400 protostellar candidates, including one massive young stellar object (YSO) that has not been previously highlighted. The YSO distribution supports a picture of distributed low-level star formation, with no strong evidence of triggered star formation in the ``pillars''. We confirm the youth of NGC 6611 by a large fraction of infrared-excess sources, and reveal a younger cluster of YSOs in the nearby molecular cloud. Analysis of the YSO clustering properties shows a possible imprint of the molecular cloud's Jeans length. Multiwavelength mid-IR imaging thus allows us to analyze the protostellar population, to measure the dust temperature and column density, and to relate these in a consistent picture of star formation in M16.

We compare H-alpha, radio continuum, and Spitzer Space Telescope (SST) images of 58 planetary nebulae (PNe) recently discovered by the Macquarie-AAO-Strasbo- urg H-alpha PN Project (MASH) of the SuperCOSMOS H-alpha Survey. Using InfraRed Array Camera (IRAC) data we define the IR colors of PNe and demonstrate good isolation between these colors and those of many other types of astronomical object. The only substantive contamination of PNe in the color-color plane we illustrate is due to YSOs. However, this ambiguity is readily resolved by the unique optical characteristics of PNe and their environs. We also examine the relationships between optical and MIR morphologies from 3.6 to 8.0um and explore the ratio of mid-infrared (MIR) to radio nebular fluxes, which is a valuable discriminant between thermal and nonthermal emission. MASH emphasizes late evolutionary stages of PNe compared with previous catalogs, enabling study of the changes in MIR and radio flux that attend the aging process. Spatially integrated MIR energy distributions were constructed for all MASH PNe observed by the GLIMPSE Legacy Project, using the H-alpha morphologies to establish the dimensions for the calculations of the Midcourse Space Experiment (MSX), IRAC, and radio continuum (from the Molonglo Observatory Synthesis Telescope and the Very Large Array) flux densities. The ratio of IRAC 8.0-um to MSX 8.3-um flux densities provides a measure of the absolute diffuse calibration of IRAC at 8.0 um. We independently confirm the aperture correction factor to be applied to IRAC at 8.0um to align it with the diffuse calibration of MSX. The result agrees with the recommendations of the Spitzer Science Center and with results from a parallel study of HII regions. These PNe probe the diffuse calibration of IRAC on a spatial scale of 9-77 arcsec.

We have observed the central region of the IR-dark cloud filament associated with IRAS 18507+0121 at millimeter wavelengths in CO(1-0), 13CO(1-0), and C18O(1-0) line emission and with the Spitzer Space Telescope at mid-IR wavelengths. Five massive outflows from two cloud cores were discovered. Three outflows are centered on or near an Ultracompact HII region (G34.4+0.23) while the remaining two outflows originate from the millimeter core G34.4+0.23 MM. Modeling of the SEDs of the mid-IR sources identified 31 young stellar objects in the filament with a combined stellar mass of ~127 +/- 27 Msun. An additional 22 sources were identified as probable cluster members based on the presence of strong 24 micron emission. The total star formation efficiency in the G34.4 cloud filament is estimated to be ~7% while the massive and intermediate mass star formation efficiency in the entire cloud filament is estimated to be roughly 2%. A comparison of the gravitational binding energy with the outflow kinetic energy suggests that the compact core containing G34.4+0.23 MM is being destroyed by its molecular outflows whereas the outflows associated with more massive core surrounding the G34.4 UC HII region are not likely to totally disrupt the cloud. Additionally, a qualitative evaluation of the region appears to suggest that stars in this region may have formed in two stages: first lower mass stars formed and then, a few Myrs later, the more massive stars began to form.

We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.

We report on the discovery of two spiral galaxies located behind the southern Milky Way, within the least explored region of the Great Attractor. They lie at 317, -0.5 deg galactic, where obscuration from Milky Way stars and dust exceeds 13 to 15 mag of visual extinction. The galaxies were the most prominent of a set identified using mid-infrared images of the low-latitude (|b| < 1 deg) Spitzer Legacy program GLIMPSE. Follow-up HI radio observations reveal that both galaxies have redshifts that place them squarely in the Norma Wall of galaxies, which appears to extend diagonally across the Galactic Plane from Norma in the south to Centaurus & Vela in the north. We report on the near-infrared, mid-infrared and radio properties of these newly discovered galaxies, and discuss their context in the larger view of the Great Attractor. The work presented here demonstrates that mid-infrared surveys open up a new window to study galaxies in the Zone of Avoidance.

We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer Space Telescope at 8.0microns using a sample of 43 HII regions with a wide range of morphologies near GLON=312deg. For each region we carefully measure sky-subtracted,point-source- subtracted, areally-integrated IRAC 8.0-micron fluxes and compare these with Midcourse Space eXperiment (MSX) 8.3-micron images at two different spatial resolutions, and with radio continuum maps. We determine an accurate median ratio of IRAC 8.0-micron/MSX\8.3-micron fluxes, of 1.55+/-0.15. From robust spectral energy distributions of these regions we conclude that the present 8.0-micron diffuse calibration of the SST is 36% too high compared with the MSX validated calibration, perhaps due to scattered light inside the camera. This is an independent confirmation of the result derived for the diffuse calibration of IRAC by the Spitzer Science Center (SSC). From regression analyses we find that 843-MHz radio fluxes of HII regions and mid-infrared (MIR) fluxes are linearly related for MSX at 8.3-microns and Spitzer at 8.0 microns, confirming the earlier MSX result by Cohen & Green. The median ratio of MIR/843-MHz diffuse continuum fluxes is 600 times smaller in nonthermal than thermal regions, making it a sharp discriminant. The ratios are largely independent of morphology up to a size of ~24 arcsec. We provide homogeneous radio and MIR morphologies for all sources. MIR morphology is not uniquely related to radio structure. Compact regions may have MIR filaments and/or diffuse haloes, perhaps infrared counter- parts to weakly ionized radio haloes found around compact HII regions. We offer two IRAC colour-colour plots as quantitative diagnostics of diffuse HII regions.

We present the results of a search for new members of the Taurus star-forming region using the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope}. With IRAC images of 29.7 deg^2 of Taurus at 3.6, 4.5, 5.8, and 8.0 um, we have identified sources with red mid-infrared colors indicative of disk-bearing objects and have obtained optical and infrared spectra of 23 of these candidate members. Through this work, we have discovered 13 new members of Taurus, two of which have spectral types later than M6 and thus are likely to be brown dwarfs according to the theoretical evolutionary models of Chabrier and Baraffe. This survey indicates that the previous census of Taurus has a completeness of ~80% for members with disks. The new members that we have found do not significantly modify the previously measured distributions of Taurus members as a function of position, mass, and extinction. For instance, we find no evidence for a population of highly reddened brown dwarfs (A_K~2) that has been missed by previous optical and near-infrared surveys, which suggests that brown dwarf disks are not significantly more flared than disks around stars. In addition to the new members, we also present IRAC photometry for the 149 previously known members that appear within this survey, which includes 27 objects later than M6.

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

Diffuse emission in the mid-infrared shows a wealth of structure, that lends itself to high-resolution structure analysis of the interstellar gas. A large part of the emission comes from polycyclic aromatic hydrocarbons, excited by nearby ultra-violet sources. Can the observed diffuse emission structure be interpreted as column density structure? We discuss this question with the help of a set of model molecular clouds bathed in the radiation field of a nearby O-star. The correlation strength between column density and ``observed'' flux density strongly depends on the absolute volume density range in the region. Shadowing and irradiation effects may completely alter the appearance of an object. Irradiation introduces additional small-scale structure and it can generate structures resembling shells around HII-regions in objects that do not possess any shell-like structures whatsoever. Nevertheless, structural information about the underlying interstellar medium can be retrieved. In the more diffuse regime ($n({HI})\lesssim 100$cm$^{-3}$), flux density maps may be used to trace the 3D density structure of the cloud via density gradients. Thus, while caution definitely is in order, mid-infrared surveys such as GLIMPSE will provide quantitative insight into the turbulent structure of the interstellar medium.

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 the results of long-term spectropolarimetric and spectroscopic monitoring of MWC 314, a candidate Luminous Blue Variable star. We detect the first evidence of H$\alpha$ variability in MWC 314, and find no apparent periodicity in this emission. The total R-band polarization is observed to vary between 2.21% and 3.00% at a position angle consistently around $\sim0^{\circ}$, indicating the presence of a time-variable intrinsic polarization component, hence an asymmetrical circumstellar envelope. We find suggestive evidence that MWC 314's intrinsic polarization exhibits a wavelength-independent magnitude varying between 0.09% and 0.58% at a wavelength-independent position angle covering all four quadrants of the Stokes Q-U plane. Electron scattering off of density clumps in MWC 314's wind is considered as the probable mechanism responsible for these variations.

Using Infrared Array Camera (IRAC) images at 3.6, 4.5, 5.8, and 8 microns from the GLIMPSE Legacy science program on the Spitzer Space Telescope, we searched for infrared counterparts to the 95 known supernova remnants that are located within galactic longitudes 65>|l|>10 degrees and latitudes |b|<1 degree. Eighteen infrared counterparts were detected. Many other supernova remnants could have significant infrared emission but are in portions of the Milky Way too confused to allow separation from bright HII regions and pervasive mid-infrared emission from atomic and molecular clouds along the line of sight. Infrared emission from supernova remnants originates from synchrotron emission, shock-heated dust, atomic fine-structure lines, and molecular lines. The detected remnants are G11.2-0.3, Kes 69, G22.7-0.2, 3C 391, W 44, 3C 396, 3C 397, W 49B, G54.4-0.3, Kes 17, Kes 20A, RCW 103, G344.7-0.1, G346.6-0.2, CTB 37A, G348.5-0.0, and G349.7+0.2. The infrared colors suggest emission from molecular lines (9 remnants), fine-structure lines (3), and PAH (4), or a combination; some remnants feature multiple colors in different regions. None of the remnants are dominated by synchrotron radiation at mid-infrared wavelengths. The IRAC-detected sample emphasizes remnants interacting with relatively dense gas, for which most of the shock cooling occurs through molecular or ionic lines in the mid-infrared.

The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0), where l_0 is insensitive to limiting magnitude, band choice, and side of Galactic center: l_0= 17-30 degrees with a best fit value in the the 4.5 micron band of l_0=24 +/- 4 degrees. Modeling the source distribution as an exponential disk yields a radial scale length of H= 3.9 +/- 0.6 kpc. There is a pronounced north-south asymmetry in source counts for |l| < 30 degrees, with ~ 25% more stars in the north. For l=10-30 degrees, there is a strong enhancement of stars of m= 11.5-13.5 mag. A linear bar passing through the Galactic center with half-length R_bar=4.4 +/- 0.5 kpc, tilted by phi=44 +/- 10 degrees to the Sun-Galactic Center line, provides the simplest interpretation of this data. We examine the possibility that enhanced source counts at l=26-28 degrees, 31.5-34 degrees, and 306-309 degrees are related to Galactic spiral structure. Total source counts are depressed in regions where the counts of red objects (m_K-m_[8.0] >3) peak. In these areas, the counts are reduced by extinction due to molecular gas and/or high diffuse backgrounds associated with star formation.

Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K.

We present a new planetary nebula, first identified in images from the Australia Telescope Compact Array, although not recognized at that time. Recent observations with the Spitzer Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high-resolution radio and infrared images enable the identification of the central star or its wind, the recognition of the radio emission as thermal, and the probable presence of polycylic aromatic hydrocarbons in and around the source. These lead to the conclusion that G313.3+00.3 is a planetary nebula. This object is of particular interest because it was discovered solely through radio and mid-infrared imaging, without any optical (or near-infrared) confirmation, and acts as a proof of concept for the discovery of many more highly extinguished planetary nebulae. G313.3+00.3 is well-resolved by both the instruments with which it was identified, and suffers extreme reddening due to its location in the Scutum-Crux spiral arm.

We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitudes were imaged at J, H, and K bands by 2MASS and at 3--8um by Spitzer for the GLIMPSE Legacy program. We measure the mean values of the color excess ratios (A(lambda)-A(K))/(A(J)-A(K)) directly from the color distributions of observed stars. The extinction ratio between two of the filters, e.g. A(J)/A(K), is required to calculate A(lambda)/A(K) from those measured ratios. We use the apparent JHK magnitudes of giant stars along our two sightlines, and fit the reddening as a function of magnitude (distance) to determine A(J)/kpc, A(K)/kpc, and A(J)/A(K). Our values of A(lambda)/A(K) show a flattening across the 3--8um wavelength range, roughly consistent with the Lutz 1996 extinction measurements derived for the sightline toward the Galactic center.

Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity.

GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions.

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.

GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse.

We report on new data, including four vacuum ultraviolet spectropolarimetric observations by the Wisconsin Ultraviolet Photo-Polarimeter Experiment ("WUPPE") on the Astro-1 (one observation) and Astro-2 (three observations) shuttle missions, and 15 new visible-wavelength observations obtained by the HPOL CCD spectropolarimeter at the Pine Bluff Observatory of the University of Wisconsin. This includes three HPOL observations made within 12 hours of each of the three Astro-2 WUPPE observations, giving essentially simultaneous observations extending from 1500 to 10500 Angstroms. An analysis of these data yields estimates on the properties of wind "clumps" when they are detectable by the polarization of their scattered light. We find that the clumps must be detected near the base of the wind, r/R* = 1.3 - 2.5. At this time the clump density is at least 10^13 cm^-3, 20 times the mean wind, and the temperature is roughly 10,000 K, about 20% cooler than the mean wind. The clumps leading to the largest observed polarization must have electron optical depth of 0.1 - 1 and an angular extent of 0.1 - 1 ster, and account for at most 2% of the wind mass loss. The H and HeI emission lines from the wind are unpolarized, but their P Cygni absorption is enhanced by a factor of 4 in the scattered light. We speculate on the possible relationship of the clumps detected polarimetrically to those seen as Discrete Absorption Features (DACs) and in interferometry.