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

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

The James Webb Space Telescope (JWST) offers unprecedented sensitivity, stability, and wavelength coverage for transiting exoplanet studies, opening up new avenues for measuring atmospheric abundances, structure, and temperature profiles. Taking full advantage of JWST spectroscopy of planets from 0.6um to 28um, however, will require many observations with a combination of the NIRISS, NIRCam, NIRSpec, and MIRI instruments. In this white paper, we discuss a new NIRCam mode (not yet approved or implemented) that can reduce the number of necessary observations to cover the 1.0um to 5.0um wavelength range. Even though NIRCam was designed primarily as an imager, it also includes several grisms for phasing and aligning JWST's 18 hexagonal mirror segments. NIRCam's long-wavelength channel includes grisms that cover 2.4um to 5.0um with a resolving power of R = 1200 - 1550 using two separate configurations. The long-wavelength grisms have already been approved for science operations, including wide field and single object (time series) slitless spectroscopy. We propose a new mode that will simultaneously measure spectra for science targets in the 1.0um to 2.0um range using NIRCam's short-wavelength channel. This mode, if approved, would take advantage of NIRCam's Dispersed Hartmann Sensor (DHS), which produces 10 spatially separated spectra per source at R ~ 300. We discuss the added benefit of the DHS in constraining abundances in exoplanet atmospheres as well as its ability to observe the brightest systems. The DHS essentially comes for free (at no time cost) with any NIRCam long-wavelength grism observation, but the detector integration parameters have to be selected to ensure that the long-wavelength grism observations do not saturate and that JWST data volume downlink constraints are not violated.

An eruptive nova-like event took place in 1994 in the stellar-merger candidate V4332 Sgr. Following the eruption, dust consisting of refractory silicate rich dust grains containing a significant component of AlO bonding was formed sometime between 1998 and 2003. Observations using Spitzer between 2005 and 2009 show significant changes in the 10 micron silicate stretch feature. There is a deepening of the 10 micron silicate stretch as well as the development of a feature between about 13 and 20 microns consistent with a blend of the MgO and FeO stretching features and the O-Si-O bending mode of increasingly ordered silicate dust. Near-infrared observations show the presence of AlO and water vapor in the outflow in 2003, 2004 and 2005: the AlO has significantly decreased in spectra obtained in 2014 while the water vapor remains largely unchanged. An attempt is made to correlate these observations and understand the significance of these changes using DUSTY modeling. The observations appear consistent with the kinetically-controlled, condensation of highly under-oxidized SiO/AlO/Fe/Mg dust grains in the outflow followed by the continuous evolution of the initial condensate due to thermal annealing and oxidation of the dust via reaction with ambient O, OH and H2O in the expanding, cooling shell. Periodic monitoring of this dust shell over the mid-infrared spectral range could yield useful information on the evolution of under-oxidized silicate condensates exposed to hot water vapor in more conventional circumstellar environments.

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

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

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

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 present a study of the infrared/submm emission of the LMC star forming complex N158-N159-N160. Combining observations from the Spitzer Space Telescope (3.6-70um), the Herschel Space Observatory (100-500um) and LABOCA (870um) allows us to work at the best angular resolution available now for an extragalactic source. We observe a remarkably good correlation between SPIRE and LABOCA emission and resolve the low surface brightnesses emission. We use the Spitzer and Herschel data to perform a resolved Spectral Energy Distribution (SED) modelling of the complex. Using MBB, we derive a global emissivity index beta_c of 1.47. If beta cold is fixed to 1.5, we find an average temperature of 27K. We also apply the Galliano et al. (2011) modelling technique (and amorphous carbon to model carbon dust) to derive maps of the star formation rate, the mean starlight intensity, the fraction of PAHs or the dust mass surface density of the region. We observe that the PAH fraction strongly decreases in the HII regions. This decrease coincides with peaks in the mean radiation field intensity map. The dust surface densities follow the FIR distribution, with a total dust mass of 2.1x10^4 Msolar (2.8 times less than when using graphite grains) in the resolved elements we model. We find a non-negligible amount of dust in the molecular cloud N159 South (showing no massive SF). We also investigate the drivers of the Herschel/PACS and SPIRE submm colours as well as the variations in the gas-to-dust mass ratio (G/D) and the XCO conversion factor in the region N159. We finally model individual regions to analyse variations in the SED shape across the complex and the 870um emission in more details. No measurable submm excess emission at 870um seems to be detected in these regions.

Vega and Fomalhaut, are similar in terms of mass, ages, and global debris disk properties; therefore, they are often referred as "debris disk twins". We present Spitzer 10-35 um spectroscopic data centered at both stars, and identify warm, unresolved excess emission in the close vicinity of Vega for the first time. The properties of the warm excess in Vega are further characterized with ancillary photometry in the mid infrared and resolved images in the far-infrared and submillimeter wavelengths. The Vega warm excess shares many similar properties with the one found around Fomalhaut. The emission shortward of ~30 um from both warm components is well described as a blackbody emission of ~170 K. Interestingly, two other systems, eps Eri and HR 8799, also show such an unresolved warm dust using the same approach. These warm components may be analogous to the solar system's zodiacal dust cloud, but of far greater. The dust temperature and tentative detections in the submillimeter suggest the warm excess arises from dust associated with a planetesimal ring located near the water-frost line and presumably created by processes occurring at similar locations in other debris systems as well. We also review the properties of the 2 um hot excess around Vega and Fomalhaut, showing that the dust responsible for the hot excess is not spatially associated with the dust we detected in the warm belt. We suggest it may arise from hot nano grains trapped in the magnetic field of the star. Finally, the separation between the warm and cold belt is rather large with an orbital ratio >~10 in all four systems. In light of the current upper limits on the masses of planetary objects and the large gap, we discuss the possible implications for their underlying planetary architecture, and suggest that multiple, low-mass planets likely reside between the two belts in Vega and Fomalhaut.

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.

A multi-epoch H-alpha survey of the late-type spiral galaxy NGC 2403 has been completed in order to determine its nova rate. A total of 9 nova candidates were discovered in 48 nights of observation with two different telescopes over the period from February 2001 to April 2012. After making corrections for temporal coverage and spatial completeness, a nova rate of 2.0 (+0.5,-0.3) per year was determined. This rate corresponds to a luminosity-specific nova rate of 2.5+/-0.7 novae per year per 10^{10} solar luminosities in K. This value is consistent with that of the similar Hubble type galaxy, M33, and is typical of those of other galaxies with measured nova rates, which range from 1-3 novae per year per 10^{10} solar luminosities in K.

Aims: In this paper, we perform detailed modelling of the Spitzer and Herschel observations of the LMC, in order to: (i) systematically study the uncertainties and biases affecting dust mass estimates; and to (ii) explore the peculiar ISM properties of the LMC. Methods: To achieve these goals, we have modelled the spatially resolved SEDs with two alternate grain compositions, to study the impact of different submillimetre opacities on the dust mass. We have rigorously propagated the observational errors (noise and calibration) through the entire fitting process, in order to derive consistent parameter uncertainties. Results: First, we show that using the integrated SED leads to underestimating the dust mass by ~50 % compared to the value obtained with sufficient spatial resolution, for the region we studied. This might be the case, in general, for unresolved galaxies. Second, we show that Milky Way type grains produce higher gas-to-dust mass ratios than what seems possible according to the element abundances in the LMC. A spatial analysis shows that this dilemma is the result of an exceptional property: the grains of the LMC have on average a larger intrinsic submm opacity (emissivity index beta~1.7 and opacity kappa_abs(160 microns)=1.6 m2/kg) than those of the Galaxy. By studying the spatial distribution of the gas-to-dust mass ratio, we are able to constrain the fraction of unseen gas mass between ~10, and ~100 % and show that it is not sufficient to explain the gas-to-dust mass ratio obtained with Milky Way type grains. Finally, we confirm the detection of a 500 microns extended emission excess with an average relative amplitude of ~15 %, varying up to 40 %. This excess anticorrelates well with the dust mass surface density. Although we do not know the origin of this excess, we show that it is unlikely the result of very cold dust, or CMB fluctuations.

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

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

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

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

We present 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 report the results of a multi-year spectroscopic and photometric survey of novae in M31 that resulted in a total of 53 spectra of 48 individual nova candidates. Two of these, M31N 1995-11e and M31N 2007-11g, were revealed to be long-period Mira variables, not novae. These data double the number of spectra extant for novae in M31 through the end of 2009 and bring to 91 the number of M31 novae with known spectroscopic classifications. We find that 75 novae (82%) are confirmed or likely members of the Fe II spectroscopic class, with the remaining 16 novae (18%) belonging to the He/N (and related) classes. These numbers are consistent with those found for Galactic novae. We find no compelling evidence that spectroscopic class depends sensitively on spatial position or population within M31 (i.e., bulge vs. disk), although the distribution for He/N systems appears slightly more extended than that for the Fe II class. We confirm the existence of a correlation between speed class and ejection velocity (based on line width), as in the case of Galactic novae. Follow-up photometry allowed us to determine light-curve parameters for a total of 47 of the 91 novae with known spectroscopic class. We confirm that more luminous novae generally fade the fastest, and that He/N novae are typically faster and brighter than their Fe II counterparts. In addition, we find a weak dependence of nova speed class on position in M31, with the spatial distribution of the fastest novae being slightly more extended than that of slower novae.

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 report the results of the first infrared survey of novae in the nearby spiral galaxy, M31. Both photometric and spectroscopic observations of a sample of 10 novae (M31N 2006-09c, 2006-10a, 2006-10b, 2006-11a, 2007-07f, 2007-08a, 2007-08d, 2007-10a, 2007-11d, and 2007-11e) were obtained with the Spitzer Space Telescope. Eight of the novae were observed with the IRAC (all but M31N 2007-11d and 2007-11e) and eight with the IRS (all but 2007-07f and 2007-08a), resulting in six in common between the two instruments. The observations, which were obtained between ~3 and ~7 months after discovery, revealed evidence for dust formation in two of the novae: M31N 2006-10a and (possibly) 2007-07f, and [Ne II] 12.8 micron line emission in a third (2007-11e). The Spitzer observations were supplemented with ground-based optical photometric and spectroscopic data that were used to determine the speed classes and spectroscopic types of the novae in our survey. After including data for dust-forming Galactic novae, we show that dust formation timescales are correlated with nova speed class in that dust typically forms earlier in faster novae. We conclude that our failure to detect the signature of dust formation in most of our M31 sample is likely a result of the relatively long delay between nova eruption and our Spitzer observations. Indeed, the two novae for which we found evidence of dust formation were the two "slowest" novae in our sample. Finally, as expected, we found that the majority of the novae in our sample belong to the Fe II spectroscopic class, with only one clear example of the He/N class (M31N 2006-10b). Typical of an He/N system, M31N 2006-10b was the fastest nova in our sample, not detected with the IRS, and just barely detected in three of the IRAC bands when it was observed ~4 months after eruption.

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

The LMC star, SSTISAGE1C J050756.44-703453.9, was first noticed during a survey of EROS-2 lightcurves for stars with large irregular brightness variations typical of the R Coronae Borealis (RCB) class. However, the visible spectrum showing emission lines including the Balmer and Paschen series as well as many Fe II lines is emphatically not that of an RCB star. This star has all of the characteristics of a typical UX Ori star. It has a spectral type of approximately A2 and has excited an H II region in its vicinity. However, if it is an LMC member, then it is very luminous for a Herbig Ae/Be star. It shows irregular drops in brightness of up to 2 mag, and displays the reddening and "blueing" typical of this class of stars. Its spectrum, showing a combination of emission and absorption lines, is typical of a UX Ori star that is in a decline caused by obscuration from the circumstellar dust. SSTISAGE1C J050756.44-703453.9 has a strong IR excess and significant emission is present out to 500 micron. Monte Carlo radiative transfer modeling of the SED requires that SSTISAGE1C J050756.44-703453.9 has both a dusty disk as well as a large extended diffuse envelope to fit both the mid- and far-IR dust emission. This star is a new member of the UX Ori subclass of the Herbig Ae/Be stars and only the second such star to be discovered in the LMC.

A multi-epoch H$\alpha$ survey of the early-type spiral galaxy M94 (NGC 4736) has been completed as part of a program to establish the galaxy's nova rate. A total of four nova candidates were discovered in seven epochs of observation during the period from 2005 to 2007. After making corrections for temporal coverage and spatial completeness, a global nova rate of 5.0$^{+1.8}_{-1.4}$ yr$^{-1}$ was determined. This rate corresponds to a specific-luminosity nova rate of 1.4 $\pm$ 0.5 novae per year per 10$^{10} L_{\odot,K}$ when the $K$ luminosity is determined from the $B-K$ color, or 1.5 $\pm$ 0.4 novae per year per 10$^{10} L_{\odot,K}$ when the $K$ luminosity is derived from the Two Micron All Sky Survey. These values are slightly lower than that of other galaxies with measured nova rates, which typically lie in the range of $2-3$ novae per year per 10$^{10} L_{\odot}$ in the $K$ band.

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.

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

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

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

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

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

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

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

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

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

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 present new mid- to far-infrared images of the two dwarf compact elliptical galaxies that are satellites of M31, NGC 185 and NGC 147, obtained with the Spitzer Space Telescope. Spitzer's high sensitivity and spatial resolution enable us for the first time to look directly into the detailed spatial structure and properties of the dust in these systems. The images of NGC 185 at 8 and 24 micron display a mixed morphology characterized by a shell-like diffuse emission region surrounding a central concentration of more intense infrared emission. The lower resolution images at longer wavelengths show the same spatial distribution within the central 50" but beyond this radius, the 160 micron emission is more extended than that at 24 and 70micron. On the other hand, the dwarf galaxy NGC 147 located only a small distance away from NGC 185 shows no significant infrared emission beyond 24 micron and therefore its diffuse infrared emission is mainly stellar in origin. For NGC 185, the derived dust mass based on the best fit to the spectral energy distribution is 1.9e3 Msol, implying a gas mass of ~3.0e5 Msol. These values are in agreement with those previously estimated from infrared as well as from CO and HI observations and are consistent with the predicted mass return from dying stars based on the last burst of star formation ~1.0e9 yr ago. Based on the 70 to 160micron flux density ratio, we estimate a temperature for the dust of ~17K. For NGC 147, we obtain an upper limit for the dust mass of 4.5e2 Msol at 160 micron (assuming a temperature of ~20K), a value consistent with the previous upper limit derived using ISO observations of this galaxy. In the case of NGC 185, we also present full 5-38 micron low-resolution (R~100) spectra of the main emission regions.

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 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 have obtained a full suite of Spitzer observations to characterize the debris disk around HR 8799 and to explore how its properties are related to the recently discovered set of three massive planets orbiting the star. We distinguish three components to the debris system: (1) warm dust (T ~150 K) orbiting within the innermost planet; (2) a broad zone of cold dust (T ~45 K) with a sharp inner edge, orbiting just outside the outermost planet and presumably sculpted by it; and (3) a dramatic halo of small grains originating in the cold dust component. The high level of dynamical activity implied by this halo may arise due to enhanced gravitational stirring by the massive planets. The relatively young age of HR 8799 places it in an important early stage of development and may provide some help in understanding the interaction of planets and planetary debris, an important process in the evolution of our own solar system.

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

(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 bright Type II-plateau supernova (SN) 2004dj occurred within the young, massive stellar cluster Sandage-96 in a spiral arm of NGC 2403. New multi-wavelength observations obtained with several ground-based and space-based telescopes are combined to study the radiation from Sandage-96 after SN 2004dj faded away. Sandage-96 started to dominate the flux in the optical bands starting September 2006 (~800 d after explosion). The optical fluxes are equal to the pre-explosion ones within the observational uncertainties. An optical Keck spectrum obtained ~900 d after explosion shows the dominant blue continuum from the cluster stars shortward of 6000 \AA as well as strong SN nebular emission lines redward. The integrated spectral energy distribution (SED) of the cluster has been extended into the ultraviolet region by archival XMM-Newton and new Swift observations, and compared with theoretical models. The outer parts of the cluster have been resolved by the Hubble Space Telescope, allowing the construction of a color-magnitude diagram. The fitting of the cluster SED with theoretical isochrones results in cluster ages between 10--40 Myr, depending on metallicity and the model family. The isochrone fitting indicates that the resolved part of the cluster has a bimodal age distribution: a younger population at ~10--16 Myr, and an older one at ~32--100 Myr which is similar to the age distribution of the nearby field stars. These stars may have been captured from the field during the cluster formation. The young age of Sandage-96 suggest 12 < M_prog < 20 M_\odot as the most probable mass range for the progenitor of SN 2004dj. This is consistent with, but perhaps slightly higher than, most of the other Type II-plateau SN progenitor masses determined so far.

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 Spitzer infrared observations of the starless core L429. The IR images of this core show an absorption feature, caused by the dense core material, at wavelengths <= 70 micron. The core has a steep density profile, and reaches A_V > 35 mag near the center. We show that L429 is either collapsing or in a near-collapse state.

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

We are currently involved in a multifaceted campaign to study extragalactic classical novae in the Local Group and beyond. Here we report on-going results from the exploitation of the POINT-AGAPE M31 dataset; initial results from our Local Group imaging, and spectroscopic CNe follow-up campaign and introduce the Liverpool Extragalactic Nova Survey.

We report a series of extensive photometric and spectroscopic observations of the luminous M31 nova M31N 2007-11d. Our photometric observations coupled with previous measurements show that the nova took at least four days to reach peak brightness at R~14.9 on 20 Nov 2007 UT. After reaching maximum, the time for the nova to decline 2 and 3 magnitudes from maximum light (t_2 and t_3) was ~9.5 and ~13 days, respectively, establishing that M31N 2007-11d was a moderately fast declining nova. During the nova's evolution a total of three spectra were obtained. The first spectrum was obtained one day after maximum light (5 days post-discovery), followed by two additional spectra taken on the decline at two and three weeks post-maximum. The initial spectrum reveals narrow Balmer and Fe II emission with P Cygni profiles superimposed on a blue continuum. These data along with the spectra obtained on the subsequent decline clearly establish that M31N 2007-11d belongs to the Fe II spectroscopic class. The properties of M31N 2007-11d are discussed within the context of other luminous novae in M31, the Galaxy, and the LMC. Overall, M31N 2007-11d appears to be remarkably similar to Nova LMC 1991, which was another bright, slowly-rising, Fe II nova. A comparison of the available data for luminous extragalactic novae suggest that the >~4 day rise to maximum light seen in M31N 2007-11d may not be unusual, and that the rise times of luminous Galactic novae, usually assumed to be <~2 days, may have been underestimated.

A new multi-epoch H-alpha imaging study of M101 (NGC 5457) has been carried out as part of a larger campaign to study the rate and stellar population of extragalactic novae. The survey yielded a total of 13 nova detections from 10 epochs of M101 observations spanning a three year period. After correcting for the temporal coverage and survey completeness, a global nova rate of 11.7^{+1.9}_{-1.5} per year is found. This value corresponds to a luminosity-specific nova rate of 1.23\pm0.27 novae per year per 10^{10} L_{\sun,K} when the K luminosity is derived from the B-K color, or 1.94\pm0.42 novae per year per 10^{10} L_{\sun,K} when the K magnitude from the Two Micron All Sky Survey is used. These values are consistent with previous estimates by Shafter et al. that were based on more limited data. The spatial distribution of the combined nova sample from the present survey and from the earlier Shafter et al. survey shows that the specific frequency of novae closely follows the integrated background light of the galaxy.

Cassiopeia A is one of the youngest supernova remnants known in the Milky Way and a unique laboratory for supernova physics. We present an optical spectrum of the Cassiopeia A supernova near maximum brightness, obtained from observations of a scattered light echo - more than three centuries after the direct light of the explosion swept past Earth. The spectrum shows that Cassiopeia A was a type IIb supernova and originated from the collapse of the helium core of a red supergiant that had lost most of its hydrogen envelope prior to exploding. Our finding concludes a longstanding debate on the Cassiopeia A progenitor and provides new insight into supernova physics by linking the properties of the explosion to the wealth of knowledge about its remnant.

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

We present a three-year series of observations at 24 microns with the Spitzer Space Telescope of the interstellar material in a 200 x 200 arcmin square area centered on Cassiopeia A. Interstellar dust heated by the outward light pulse from the supernova explosion emits in the form of compact, moving features. Their sequential outward movements allow us to study the complicated three-dimensional structure of the interstellar medium (ISM) behind and near Cassiopeia A. The ISM consists of sheets and filaments, with many structures on a scale of a parsec or less. The spatial power spectrum of the ISM appears to be similar to that of fractals with a spectral index of 3.5. The filling factor for the small structures above the spatial wavenumber k ~ 0.5 cycles/pc is only ~ 0.4%.

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 broad-band 24, 70 and 160 micron photometry, 5-35 micron and 55-90 micron spectra of the eruptive variable V4332 Sgr from Spitzer observations. The distinguishing feature of the 5-35 micron spectrum is an unusually broad absorption feature near 10 micron at the position generally associated with silicate-rich dust. Through radiative transfer modeling, we show that this broad feature cannot arise from silicates alone but requires the inclusion of alumina (Al2O3) as a dust condensate. The case for including Al2O3 is strengthened further by the presence of the AlO radical, a potentially important molecule in forming Al2O3. The present detection indicates that porous alumina manifests itself through a broadening of the 9.7 micron silicate feature and additionally displays, on the shoulder of the silicate feature, a component at ~11.5 micron. We discuss how further observations of V4332 Sgr may have the potential of verifying some general predictions of the dust condensation process.

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

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

We examine the infrared properties of five low surface brightness galaxies (LSBGs) and compare them with related but higher surface brightness galaxies, using Spitzer Space Telescope images and spectra. All the LSBGs are detected in the 3.6 and 4.5um bands, representing the stellar population. All but one are detected at 5.8 and 8.0um, revealing emission from hot dust and aromatic molecules, though many are faint or point-like at these wavelengths. Detections of LSBGs at the far-infrared wavelengths, 24, 70, and 160um, are varied in morphology and brightness, with only two detections at 160um, resulting in highly varied spectral energy distributions. Consistent with previous expectations for these galaxies, we find that detectable dust components exist for only some LSBGs, with the strength of dust emission dependent on the existence of bright star forming regions. However, the far-infrared emission may be relatively weak compared with normal star-forming galaxies.

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

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

Multicolor (BVRI) light curves of the eclipsing classical nova V Per are presented, and a total of twelve new eclipse timings are measured for the system. When combined with previous eclipse timings from the literature, these timings yield a revised ephemeris for the times of mid-eclipse given by HJD = 2,447,442.8260(1) + 0.107123474(3) E. The eclipse profiles are analyzed with a parameter-fitting model that assumes four sources of luminosity: a white dwarf primary star, a main-sequence secondary star, a flared accretion disk with a rim, and a bright spot at the intersection of the mass-transfer stream and the disk periphery. A matrix of model solutions are computed, covering an extensive range of plausible parameter values. The solution matrix is then explored to determine the optimum values for the fitting parameters and their associated errors. For models that treat the accretion disk as a flat structure without a rim, optimum fits require that the disk have a flat temperature profile. Although models with a truncated inner disk (R_in >> R_wd) result in a steeper temperature profile, steady-state models with a temperature profile characterized by T(r) \propto r^{-3/4} are found only for models with a significant disk rim. A comparison of the observed brightness and color at mid-eclipse with the photometric properties of the best-fitting model suggests that V Per lies at a distance of ~ 1 kpc.

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

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

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

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

We present Spitzer observations of the unusual variable V838 Monocerotis. Extended emission is detected around the object at 24, 70 and 160um. The extended infrared emission is strongly correlated spatially with the HST optical light echo images taken at a similar epoch. We attribute this diffuse nebulosity to be from an infrared light echo caused by reprocessed thermal emission from dust heated by the outward-propagating radiation from the 2002 eruption. The detection of an IR light echo provides an opportunity to estimate the mass in dust of the echo material and hence constrain its origin. We estimate the dust mass of the light echo to be on the order of a solar mass - thereby implying the total gas plus dust mass to be considerably more - too massive for the echo material to be the ejecta from previous outburst/mass-losing events. This is therefore suggestive that a significant fraction of the matter seen through the light echo is interstellar in origin. Unresolved emission at 24 and 70um is also seen at the position of the central star possibly indicating the presence of hot dust freshly condensed in the outburst ejecta.

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

We are currently undertaking a Spitzer GTO program to image ~30 massive lensing clusters at moderate redshift with both IRAC and MIPS. By taking advantage of the gravitatinoal lensing power of these clusters, we will study the population of faint galaxies that are below the nominal Spitzer detection limits. Here, we present a few examples of our science programs.

We have studied the infrared properties of the brightest cluster galaxies (BCGs) located in the cores of X-ray-luminous clusters at 0.15 < z < 0.35. The majority of the BCGs are not particularly infrared-luminous compared with other massive early-type galaxies, suggesting that the cluster environment has little influence on the infrared luminosities of the BCGs. The exceptions, however, are the BCGs in the three X-ray-brightest clusters in the sample, A1835, Z3146, and A2390. These BCGs have a prominent far-infrared peak in their spectral energy distributions (SEDs), and two of them (those in A1835 and Z3146) can be classified as luminous infrared galaxies (LIRGs: L_{IR} > 10^{11} L_{sun}). Although radio AGNs are found to be prevalent among the BCGs, the infrared luminosities of these three BCGs, judged from the infrared SED signatures, are likely to be powered by star formation. Considering the overall trend that clusters with shorter radiative gas cooling times harbor more infrared-luminous BCGs, the enhanced star formation may be caused by the cooling cluster gas accreting onto the BCGs.

The dust in the Small Magellanic Cloud (SMC), an ideal analog of primordial galaxies at high redshifts, differs markedly from that in the Milky Way by exhibiting a steeply rising far-ultraviolet extinction curve, an absence of the 2175 Angstrom extinction feature, and a local minimum at ~12 micron in its infrared emission spectrum, suggesting the lack of ultrasmall carbonaceous grains (i.e. polycyclic aromatic hydrocarbon molecules) which are ubiquitously seen in the Milky Way. While current models for the SMC dust all rely heavily on silicates, recent observations of the SMC sightline toward Sk 155 indicated that Si and Mg are essentially undepleted and the depletions of Fe range from mild to severe, suggesting that metallic grains and/or iron oxides, instead of silicates, may dominate the SMC dust. However, in this Letter we apply the Kramers-Kronig relation to demonstrate that neither metallic grains nor iron oxides are capable of accounting for the observed extinction; silicates remain as an important contributor to the extinction, consistent with current models for the SMC dust.

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

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

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

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

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

Spitzer photometry and spectroscopy of the star HD69830 reveal an excess of emission relative to the stellar photosphere between 8 and 35 $\mu$m dominated by strong features attributable to crystalline silicates with an emitting surface area more than 1,000 times that of our zodiacal cloud. The spectrum closely resembles that of the comet C/1995 O1 (Hale-Bopp). Since no excess is detected at 70 $\mu$m, the emitting material must be quite warm, be confined within a few AU of the star, and originate in grains with low long-wavelength emissivity, i.e. grains much smaller than $70 \mu$m$/2\pi\sim 10 $ $\mu$m. The strong mineralogical features are evidence for even smaller, possibly sub-micron sized grains. This small grain size is in direct contrast to the 10-100 $\mu$m grains that dominate the relatively featureless spectra of our zodiacal dust cloud and most other main sequence stars with excesses. The upper limit at 70 $\mu$m also implies that any Kuiper Belt analog must either be very cold or be less massive than $\sim$ 5 times our own Kuiper Belt. With collisional and Poynting-Robertson drag times of about a thousand years for small grains, the emitting material must either: a) be created through continual grinding down of material in a dense asteroid belt; or b) originate in cometary debris arising from either a single ``super-comet'' or a very large number of individual comets arriving from a distant reservoir. In the case of a cometary origin for the emission, the mass requirements for continuous generation by many individual comets are unreasonable and we favor the capture of a single ``super comet'' into a 0.5-1 AU orbit where it can evolve a large number of small grains over a 2 Myr period.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

We present time-resolved HST UV spectroscopy and ground-based optical photometry of the dwarf nova IP Peg in a quiescent state. The observations were obtained prior to an eclipse, when the bright spot caused by the impact of the accretion stream with the edge of the disk dominates the light output. The optical light curve is strongly correlated with the UV spectrophotometric flux curve. An emission-like feature near 1820 A in the UV spectrum is likely to be a manifestation of the ``Fe II curtain.'' Composite spectra constructed from the peaks and troughs of flickers in the light curve show substantial differences. The spectrum of the flickers is not adequately modelled by a simple blackbody, suggesting that a more sophisticated model is appropriate. We perform a cross-correlation analysis of the variability in spectrophotometric flux curves of the UV continuum and prominent UV emission lines (C II 1335, Si IV 1400, C IV 1550). The continuum and lines are not correlated, suggesting that they are produced separately. The C II and Si IV lines are moderately correlated with each other, but neither line is correlated with C IV, suggesting that the latter forms in a different region than the former. We briefly discuss a qualitative model for the geometry of the emission regions in IP Peg that is consistent with the observed behavior of the UV lines and continuum.