We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM by star formation, reheating, inflows and outflows. Our model successfully reproduces the observed dust mass function at redshift z = 0 and the observed scaling relations for dust across a wide range of redshifts. We find that the dust mass content in the present Universe is mainly produced via grain growth in the interstellar medium (ISM). By contrast, in the early Universe, the primary production mechanism for dust is the condensation in stellar ejecta. The shift of the significant production channel for dust characterises the scaling relations of dust-to-gas (DTG) and dust-to-metal (DTM) ratios. In galaxies where the grain growth dominates, we find positive correlations for DTG and DTM ratios with both metallicity and stellar mass. On the other hand, in galaxies where dust is produced primarily via condensation, we find negative or no correlation for DTM and DTG ratios with either metallicity or stellar mass. In agreement with observation showing that the circumgalactic medium (CGM) contains more dust than the ISM, our model also shows the same trend for z < 4. Our semi-analytic model is publicly available at https: //github.com/dptriani/dusty-sage.

Spitzer's final Infrared Array Camera (IRAC) observations of SN 1987A show the 3.6 and 4.5 $\mu$m emission from the equatorial ring (ER) continues a period of steady decline. Deconvolution of the images reveals that the emission is dominated by the ring, not the ejecta, and is brightest on the west side. Decomposition of the marginally resolved emission also confirms this, and shows that the west side of the ER has been brightening relative to the other portions of the ER. The infrared (IR) morphological changes resemble those seen in both the soft X-ray emission and the optical emission. The integrated ER light curves at 3.6 and 4.5 $\mu$m are more similar to the optical light curves than the soft X-ray light curve, though differences would be expected if dust is responsible for this emission and its destruction is rapid. Future observations with the James Webb Space Telescope will continue to monitor the ER evolution, and will reveal the true spectrum and nature of the material responsible for the broadband emission at 3.6 and 4.5 $\mu$m. The present observations also serendipitously reveal a nearby variable source, subsequently identified as a Be star, that has gone through a multi-year outburst during the course of these observations.

We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit.

Infrared observations of active galactic nucleus (AGN) reveal emission from the putative dusty circumnuclear 'torus' invoked by AGN unification, that is heated up by radiation from the central accreting black hole (BH). The strong 9.7 and 18 micron silicate features observed in the AGN spectra both in emission and absorption, further indicate the presence of such dusty environments. We present detailed calculations of the chemistry of silicate dust formation in AGN accretion disk winds. The winds considered herein are magnetohydrodynamic (MHD) winds driven off the entire accretion disk domain that extends from the BH vicinity to the radius of BH influence, of order of 1 to 100 pc depending on the mass of the resident BH. Our results indicate that these winds provide conditions conducive to the formation of significant amounts of dust, especially for objects accreting close to their Eddington limit, making AGN a significant source of dust in the universe, especially for luminous quasars. Our models justify the importance of a r to the power -1 density law in the winds for efficient formation and survival of dust grains. The dust production rate scales linearly with the mass of the central BH and varies as a power law of index between 2 to 2.5 with the dimensionless mass accretion rate. The resultant distribution of the dense dusty gas resembles a toroidal shape, with high column density and optical depths along the equatorial viewing angles, in agreement with the AGN unification picture.

The Central Molecular Zone (CMZ), covering the inner ~1$^\circ$ of the Galactic plane has been mapped at 2 mm using the GISMO bolometric camera on the 30 m IRAM telescope. The $21''$ resolution maps show abundant emission from cold molecular clouds, from star forming regions, and from one of the Galactic center nonthermal filaments. In this work we use the Herschel Hi-GAL data to model the dust emission across the Galactic center. We find that a single-temperature fit can describe the 160 -- 500 $\mu$m emission for most lines of sight, if the long-wavelength dust emissivity scales as $\lambda^{-\beta}$ with $\beta \approx 2.25$. This dust model is extrapolated to predict the 2 mm dust emission. Subtraction of the model from the GISMO data provides a clearer look at the 2 mm emission of star-forming regions and the brightest nonthermal filament.

We have used the Goddard IRAM 2-Millimeter Observer (GISMO) with the 30 m IRAM telescope to carry out a 2 mm survey of the Galaxy's central molecular zone (CMZ). These observations detect thermal emission from cold ISM dust, thermal free-free emission from ionized gas, and nonthermal synchrotron emission from relatively flat-spectrum sources. Archival data sets spanning $3.6 \mu$m to 90 cm are used to distinguish different emission mechanisms. After the thermal emission of dust is modeled and subtracted, the remaining 2 mm emission is dominated by free-free emission, with the exception of the brightest nonthermal filament (NTF) that runs though the middle of the bundle of filaments known as the Radio Arc. This is the shortest wavelength at which any NTF has been detected. The GISMO observations clearly trace this NTF over a length of ~0.2$^\circ$, with a mean 2 mm spectral index which is steeper than at longer wavelengths. The 2 mm to 6 cm (or 20 cm) spectral index steepens from $\alpha \approx -0.2$ to $-0.7$ as a function distance from the Sickle H II region, suggesting that this region is directly related to the NTF. A number of unresolved (at $21''$) 2 mm sources are found nearby. One appears to be thermal dust emission from a molecular cloud that is associated with an enigmatic radio point source whose connection to the Radio Arc is still debated. The morphology and colors at shorter IR wavelengths indicate other 2 mm unresolved sources are likely to be compact H II regions.

New determinations are presented of the cosmic infrared background monopole brightness in the Planck HFI bands from 100 GHz to 857 GHz. Planck was not designed to measure the monopole component of sky brightness, so cross-correlation of the 2015 HFI maps with COBE/FIRAS data is used to recalibrate the zero level of the HFI maps. For the HFI 545 and 857 GHz maps, the brightness scale is also recalibrated. Correlation of the recalibrated HFI maps with a linear combination of Galactic H I and H alpha data is used to separate the Galactic foreground emission and determine the cosmic infrared background brightness in each of the HFI bands. We obtain CIB values of 0.007 +- 0.014, 0.010 +- 0.019, 0.060 +- 0.023, 0.149 +- 0.017, 0.371 +- 0.018, and 0.576 +- 0.034 MJy/sr at 100, 143, 217, 353, 545, and 857 GHz, respectively. The estimated uncertainties for the 353 to 857 GHz bands are about 3 to 6 times smaller than those of previous direct CIB determinations at these frequencies. Our results are compared with integrated source brightness results from selected recent submillimeter and millimeter wavelength imaging surveys.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

Infrared (IR) observations of core collapse supernovae (CCSNe) have been used to infer the mass of dust that has formed in their ejecta. A plot of inferred dust masses versus supernova (SN) ages shows a trend of increasing dust mass with time, spanning a few decades of observations. This trend has been interpreted as evidence for the slow and gradual formation of dust in CCSNe. Observationally, the trend exhibits a $ t^2$ behavior, exactly what is expected from an expanding optically thick ejecta. In this case, the observed dust resides in the infrared (IR)-thin "photosphere" of the ejecta, and constitutes only a fraction of the total dust mass. We therefore propose that dust formation proceeds very rapidly, condensing most available refractory elements within two years after the explosion. At early epochs, only a fraction of the dust emission escapes the ejecta accounting for the low observed dust mass. The ejecta's entire dust content is unveiled only a few decades after the explosion, with the gradual decrease in its IR opacity. Corroborating evidence for this picture includes the early depletions of refractory elements in the ejecta of SN1987A and the appearance of a silicate emission band around day 300 in SN~2004et.

The equatorial ring of Supernova (SN) 1987A has been exposed to forward shocks from the SN blast wave, and it has been suggested that these forward shocks have been causing on-going destruction of dust in the ring. We obtained SOFIA FORCAST 11.1, 19.7 and 31.5 micron photometry of SN\,1987A in 2016. Compared with Spitzer measurements 10 years earlier, the 31.5 micron flux has significantly increased. The excess at 31.5 micron appears to be related to the Herschel 70 micron excess, which was detected 5 years earlier. The dust mass needed to account for the the 31.5--70 micron excess is 3--7x10^-4 Msun, more than ten times larger than the ring dust mass (1x10^-5 Msun) estimate from the data 10-years earlier. We argue that dust grains are re-formed or grown in the post-shock regions in the ring after forward shocks have destroyed pre-existing dust grains in the ring and released refractory elements into gas. In the post-shock region, atoms can stick to surviving dust grains, and the dust mass may have increased (grain growth), or dust grains might have condensed directly from the gas. An alternative possibility is that the outer part of the expanding ejecta dust might have been heated by X-ray emission from the circumstellar ring. The future development of this excess could reveal whether grains are reformed in the post-shocked region of the ring or eject dust is heated by X-ray.

The light curves of Type IIn supernovae are dominated by the radiative energy released through the interaction of the supernova shockwaves with their dense circumstellar medium (CSM). The ultraluminous Type IIn supernova SN 2010jl exhibits an infrared emission component that is in excess of the extrapolated UV-optical spectrum as early as a few weeks post-explosion. This emission has been attributed by some as evidence for rapid formation of dust in the cooling postshock CSM. We investigate the physical processes that may inhibit or facilitate the formation of dust in the CSM. When only radiative cooling is considered, the temperature of the dense shocked gas rapidly drops below the dust condensation temperature. However, by accounting for the heating of the postshock gas by the downstream radiation from the shock, we show that dust formation is inhibited until the radiation from the shock weakens, as the shock propagates into the less dense outer regions of the CSM. In SN 2010jl dust formation can therefore only commence after day 380. Only the IR emission since that epoch can be attributed to the newly formed CSM dust. Observations on day 460 and later show that the IR luminosity exceeds the UV-optical luminosity. The post-shock dust cannot extinct the radiation emitted by the expanding SN shock. Therefore, its IR emission must be powered by an interior source, which we identify with the reverse shock propagating through the SN ejecta. IR emission before day 380 must therefore be an IR echo from preexisting CSM dust.

We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determine SFH step intensities and combines this with a grid-based approach to determine three extinction parameters. We apply our procedure to the extensive FUV-to-FIR photometric data of M51, convolved to a common spatial resolution and pixel scale, and make the resulting maps publicly available. We recover, for M51a, a peak star formation rate (SFR) between 0.1 and 5 Gyr ago, with much lower star formation activity over the last 100 Myr. For M51b, we find a declining SFR toward the present day. In the outskirt regions of M51a, which includes regions between M51a and M51b, we recover a SFR peak between 0.1 and 1 Gyr ago, which corresponds to the effects of the interaction between M51a and M51b. We utilize our results to (1) illustrate how UV+IR hybrid SFR laws vary across M51, and (2) provide first-order estimates for how the IR luminosity per unit stellar mass varies as a function of the stellar age. From the latter result, we find that IR emission from dust heated by stars is not always associated with young stars, and that the IR emission from M51b is primarily powered by stars older than 5 Gyr.

A search for the progenitor of SN~2010jl, an unusually luminous core-collapse supernova of Type~IIn, using pre-explosion {\it Hubble}/WFPC2 and {\it Spitzer}/IRAC images of the region, yielded upper limits on the UV and near-infrared (IR) fluxes from any candidate star. These upper limits constrain the luminosity and effective temperature of the progenitor, the mass of any preexisting dust in its surrounding circumstellar medium (CSM), and dust proximity to the star. A {\it lower} limit on the CSM dust mass is required to hide a luminous progenitor from detection by {\it Hubble}. {\it Upper} limits on the CSM dust mass and constraints on its proximity to the star are set by requiring that the absorbed and reradiated IR emission not exceed the IRAC upper limits. Using the combined extinction-IR emission constraints we present viable $M_d-R_1$ combinations, where $M_d$ and $R_1$ are the CSM dust mass and its inner radius. These depend on the CSM outer radius, dust composition and grain size, and the properties of the progenitor. The results constrain the pre-supernova evolution of the progenitor, and the nature and origin of the observed post-explosion IR emission from SN~2010jl. In particular, an $\eta$~Car-type progenitor will require at least 4~mag of visual extinction to avoid detection by the {\it Hubble}. This can be achieved with dust masses $\gtrsim 10^{-3}$~\msun\ (less than the estimated 0.2-0.5~\msun\ around $\eta$~Car) which must be located at distances of $\gtrsim 10^{16}$~cm from the star to avoid detection by {\it Spitzer}.

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

While theoretical dust condensation models predict that most refractory elements produced in core-collapse supernovae (SNe) efficiently condense into dust, a large quantity of dust has so far only been observed in SN 1987A. We present the analysis of Spitzer Space Telescope, Herschel Space Observatory, Stratospheric Observatory for Infrared Astronomy (SOFIA), and AKARI observations of the infrared (IR) shell surrounding the pulsar wind nebula in the supernova remnant G54.1+0.3. We attribute a distinctive spectral feature at 21 $\mu$m to a magnesium silicate grain species that has been invoked in modeling the ejecta-condensed dust in Cas A, which exhibits the same spectral signature. If this species is responsible for producing the observed spectral feature and accounts for a significant fraction of the observed IR continuum, we find that it would be the dominant constituent of the dust in G54.1+0.3, with possible secondary contributions from other compositions, such as carbon, silicate, or alumina grains. The smallest mass of SN-formed dust required by our models is 1.1 $\pm$ 0.8 $\rm M_{\odot}$. We discuss how these results may be affected by varying dust grain properties and self-consistent grain heating models. The spatial distribution of the dust mass and temperature in G54.1+0.3 confirms the scenario in which the SN-formed dust has not yet been processed by the SN reverse shock and is being heated by stars belonging to a cluster in which the SN progenitor exploded. The dust mass and composition suggest a progenitor mass of 16$-$27 $\rm M_{\odot}$ and imply a high dust condensation efficiency, similar to that found for Cas A and SN 1987A. The study provides another example of significant dust formation in a Type IIP SN and sheds light on the properties of pristine SN-condensed dust.

The nature of the progenitor star (or system) for the Type IIn supernova (SN) subclass remains uncertain. While there are direct imaging constraints on the progenitors of at least four Type IIn supernovae, one of them being SN 2010jl, ambiguities remain in the interpretation of the unstable progenitors and the explosive events themselves. A blue source in pre-explosion HST/WFPC2 images falls within the 5 sigma astrometric error circle derived from post-explosion ground-based imaging of SN 2010jl. At the time the ground-based astrometry was published, however, the SN had not faded sufficiently for post-explosion HST follow-up observations to determine a more precise astrometric solution and/or confirm if the pre-explosion source had disappeared, both of which are necessary to ultimately disentangle the possible progenitor scenarios. Here we present HST/WFC3 imaging of the SN 2010jl field obtained in 2014 and 2015, when the SN had faded sufficiently to allow for new constraints on the progenitor. The SN, which is still detected in the new images, is offset by 0."099 +/- 0."008 (24 +/- 2 pc) from the underlying and extended source of emission that contributes at least partially, if not entirely, to the blue source previously suggested as the candidate progenitor in the WFPC2 data. This point alone rules out the possibility that the blue source in the pre-explosion images is the exploding star, but may instead suggest an association with a young (<5-6 Myr) cluster and still argues for a massive (>30 solar masses) progenitor. We obtain new upper limits on the flux from a single star at the SN position in the pre-explosion WFPC2 and Spitzer/IRAC images that may ultimately be used to constrain the progenitor properties.

Updated imaging and photometric results from Chandra observations of SN 1987A, covering the last 16 years, are presented. We find that the 0.5-2 keV light curve has remained constant at ~8x10^-12 erg s^-1 cm^-2 since 9500 days, with the 3-8 keV light curve continuing to increase until at least 10000 days. The expansion rate of the ring is found to be energy dependent, such that after day 6000 the ring expands faster in the 2-10 keV band than it does at energies <2 keV. Images show a reversal of the east-west asymmetry between 7000 and 8000 days after the explosion. The latest images suggest the southeastern side of the equatorial ring is beginning to fade. Consistent with the latest optical and infrared results, our Chandra analysis indicates the blast wave is now leaving the dense equatorial ring, which marks the beginning of a major change in the evolutionary phase of the supernova remnant 1987A.

The origin and depletion of iron differ from all other abundant refractory elements that make up the composition of the interstellar dust. Iron is primarily synthesized in Type Ia supernovae (SNe Ia) and in core collapse supernovae (CCSN), and is present in the outflows from AGB stars. Only the latter two are observed to be sources of interstellar dust, since searches for dust in SN~Ia have provided strong evidence for the absence of any significant mass of dust in their ejecta. Consequently, more than 65% of the iron is injected into the ISM in gaseous form. Yet, ultraviolet and X-ray observations along many lines of sight in the ISM show that iron is severely depleted in the gas phase compared to expected solar abundances. The missing iron, comprising about 90% of the total, is believed to be locked up in interstellar dust. This suggests that most of the missing iron must have precipitated from the ISM gas by cold accretion onto preexisting silicate, carbon, or composite grains. Iron is thus the only element that requires most of its growth to occur outside the traditional stellar condensation sources. This is a robust statement that does not depend on our evolving understanding of the dust destruction efficiency in the ISM. Reconciling the physical, optical, and chemical properties of such composite grains with their many observational manifestations is a major challenge for understanding the nature and origin of interstellar dust.

Core collapse supernovae (CCSNe) are important sources of interstellar dust, potentially capable of producing one solar mass of dust in their explosively expelled ejecta. However, unlike other dust sources, the dust has to survive the passage of the reverse shock, generated by the interaction of the supernova blast wave with its surrounding medium. Knowledge of the net amount of dust produced by CCSNe is crucial for understanding the origin and evolution of dust in the local and high-redshift universe. Our aim is to identify the dust destruction mechanisms in the ejecta, and derive the net amount of dust that survives the passage of the reverse shock. We use analytical models for the evolution of a supernova blast wave and of the reverse shock, with special application to the clumpy ejecta of the remnant of Cassiopeia A. We assume that the dust resides in cool oxygen-rich clumps that are uniformly distributed within the remnant and surrounded by a hot X-ray emitting plasma, and that the dust consists of silicates (MgSiO3) and amorphous carbon grains. The passage of the reverse shock through the clumps gives rise to a relative gas-grain motion and also destroys the clumps. Inside the ejecta clouds, dust is processed via kinetic sputtering, which is terminated either when the grains escape the clumps, or when the clumps are destroyed by the reverse shock. In either case, grain destruction proceeds thereafter by thermal sputtering in the hot ambient gas. We find that 11.8% and 15.9% of, respectively, the silicate and carbon dust survives the passage of the reverse shock by the time the shock has reached the center of the remnant. These fractions depend on the morphology of the ejecta and the medium into which the remnant is expanding, as well as the composition and size distribution of the grains that formed in the ejecta. Results will therefore differ for different types of supernovae.

Spitzer observations of SN 1987A have now spanned more than a decade. Since day ~4,000, mid-infrared (mid-IR) emission has been dominated by that from shock-heated dust in the equatorial ring (ER). From 6,000 to 8,000 days after the explosion, Spitzer observations included broadband photometry at 3.6 - 24 micron, and low and moderate resolution spectroscopy at 5 - 35 micron. Here we present later Spitzer observations, through day 10,377, which include only the broadband measurements at 3.6 and 4.5 micron. These data show that the 3.6 and 4.5 micron brightness has clearly begun to fade after day ~8,500, and no longer tracks the X-ray emission as well as it did at earlier epochs. This can be explained by the destruction of the dust in the ER on time scales shorter than the cooling time for the shocked gas. We find that the evolution of the late time IR emission is also similar to the now fading optical emission. We provide the complete record of the IR emission lines, as seen by Spitzer prior to day 8,000. The past evolution of the gas as seen by the IR emission lines seems largely consistent with the optical emission, although the IR [Fe II] and [Si II] lines show different, peculiar velocity structures.

Two millimeter observations of the MACS J1149.6+2223 cluster have detected a source that was consistent with the location of the lensed MACS1149-JD galaxy at z=9.6. A positive identification would have rendered this galaxy as the youngest dust forming galaxy in the universe. Follow up observation with the AzTEC 1.1 mm camera and the IRAM NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm have not confirmed this association. In this paper we show that the NOEMA observations associate the 2 mm source with [PCB2012] 2882 ([PCB2012] 2882 is the NED-searchable name for this source.), source number 2882 in the Hubble Space Telescope ( HST) Cluster Lensing and Supernova (CLASH) catalog of MACS J1149.6+2223. This source, hereafter referred to as CLASH 2882, is a gravitationally lensed spiral galaxy at z=0.99. We combine the GISMO 2 mm and NOEMA 1.3 mm fluxes with other (rest frame) UV to far-IR observations to construct the full spectral energy distribution (SED) of this galaxy, and derive its star formation history, and stellar and interstellar dust content. The current star formation rate of the galaxy is 54/mu Msun yr-1, and its dust mass is about 5 10^7/mu Msun, where mu is the lensing magnification factor for this source, which has a mean value of 2.7. The inferred dust mass is higher than the maximum dust mass that can be produced by core collapse supernovae (CCSN) and evolved AGB stars. As with many other star forming galaxies, most of the dust mass in CLASH 2882 must have been accreted in the dense phases of the ISM.

SN 2006gy was the most luminous SN ever observed at the time of its discovery and the first of the newly defined class of superluminous supernovae (SLSNe). The extraordinary energetics of SN 2006gy and all SLSNe (>10^51 erg) require either atypically large explosion energies (e.g., pair-instability explosion) or the efficient conversion of kinetic into radiative energy (e.g., shock interaction). The mass-loss characteristics can therefore offer important clues regarding the progenitor system. For the case of SN 2006gy, both a scattered and thermal light echo from circumstellar material (CSM) have been reported at later epochs (day ~800), ruling out the likelihood of a pair-instability event and leading to constraints on the characteristics of the CSM. Owing to the proximity of the SN to the bright host-galaxy nucleus, continued monitoring of the light echo has not been trivial, requiring the high resolution offered by the Hubble Space Telescope (HST) or ground-based adaptive optics (AO). Here we report detections of SN 2006gy using HST and Keck AO at ~3000 days post-explosion and consider the emission mechanism for the very late-time light curve. While the optical light curve and optical spectral energy distribution are consistent with a continued scattered-light echo, a thermal echo is insufficient to power the K'-band emission by day 3000. Instead, we present evidence for late-time infrared emission from dust that is radiatively heated by CSM interaction within an extremely dense dust shell, and we consider the implications on the CSM characteristics and progenitor system.

We present a new analysis of the infrared (IR) emission from the ejecta of SN1987A covering days 615, 775, 1144, 8515, and 9090 after the explosion. We show that the observations are consistent with the rapid formation of about 0.4 Msun of dust, consisting of mostly silicates, near day 615, and evolving to about 0.45 Msun of composite dust grains consisting of ~0.4 Msun of silicates and ~ 0.05 Msun of amorphous carbon after day ~8500. The proposed scenario challenges previous claims that dust in SN ejecta is predominantly carbon, and that it grew from an initial mass of ~1e-3 Msun, to over 0.5 Msun by cold accretion. It alleviates several problems with previous interpretations of the data: (1) it reconciles the abundances of silicon, magnesium, and carbon with the upper limits imposed by nucleosynthesis calculations; (2) it eliminates the requirement that most of the dust observed around day 9000 grew by cold accretion onto the1e-3 Msun of dust previously inferred for days 615 and 775 after the explosion; and (3) establishes the dominance of silicate over carbon dust in the SN ejecta. At early epochs, the IR luminosity of the dust is powered by the radioactive decay of 56Co, and at late times by at least (1.3-1.6)e-4 Msun 44Ti. Even if only a fraction greater than ~10% of the silicate dust survives the injection into the ISM, the observations firmly establish the role of core collapse SNe as the major source of thermally condensed silicate dust in the universe.

Far infrared cooling lines are ubiquitous features in the spectra of star forming galaxies. Surveys of redshifted fine-structure lines provide a promising new tool to study structure formation and galactic evolution at redshifts including the epoch of reionization as well as the peak of star formation. Unlike neutral hydrogen surveys, where the 21 cm line is the only bright line, surveys of red-shifted fine-structure lines suffer from confusion generated by line broadening, spectral overlap of different lines, and the crowding of sources with redshift. We use simulations to investigate the resulting spectral confusion and derive observing parameters to minimize these effects in pencil-beam surveys of red-shifted far-IR line emission. We generate simulated spectra of the 17 brightest far-IR lines in galaxies, covering the 150 to 1300 micron wavelength region corresponding to redshifts 0 < z < 7, and develop a simple iterative algorithm that successfully identifies the 158 micron [CII] line and other lines. Although the [CII] line is a principal coolant for the interstellar medium, the assumption that the brightest observed lines in a given line of sight are always [CII] lines is a poor approximation to the simulated spectra once other lines are included. Blind line identification requires detection of fainter companion lines from the same host galaxies, driving survey sensitivity requirements. The observations require moderate spectral resolution 700 < R < 4000 with angular resolution between 20 arcsec and 10 armin, sufficiently narrow to minimize confusion yet sufficiently large to include a statistically meaningful number of sources.

We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds (LMC and SMC, together -- the MCs) as measured by UV absorption spectroscopy. We have measured P II and Fe II along 85 sightlines toward the MCs using archival FUSE observations. For 16 of those sightlines, we have measured Si II, Cr II, and Zn II from new HST COS observations. We have combined these measurements with H I and H$_2$ column densities and reference stellar abundances from the literature to derive gas-phase abundances, depletions, and gas-to-dust ratios (GDRs). 80 of our 84 P measurements and 13 of our 16 Zn measurements are depleted by more than 0.1 decades, suggesting that P and Zn abundances are not accurate metallicity indicators at and above the metallicity of the SMC. The maximum P and Zn depletions are the same in the MW, LMC, and SMC. Si, Cr, and Fe are systematically less depleted in the SMC than in the MW or LMC. The minimum Si depletion in the SMC is consistent with zero. Our depletion-derived GDRs broadly agree with GDRs from the literature. The GDR varies from location to location within a galaxy by a factor of up to 2 in the LMC and up to 5 in the SMC. This variation is evidence of dust destruction and/or growth in the diffuse neutral phase of the interstellar medium.

Supernova generated shock waves are responsible for most of the destruction of dust grains in the interstellar medium (ISM). Calculations of the dust destruction timescale have so far been carried out using plane parallel steady shocks, however that approximation breaks down when the destruction timescale becomes longer than that for the evolution of the supernova remnant (SNR) shock. In this paper we present new calculations of grain destruction in evolving, radiative SNRs. To facilitate comparison with the previous study by Jones et al. (1996), we adopt the same dust properties as in that paper. We find that the efficiencies of grain destruction are most divergent from those for a steady shock when the thermal history of a shocked gas parcel in the SNR differs significantly from that behind a steady shock. This occurs in shocks with velocities >~ 200 km/s for which the remnant is just beginning to go radiative. Assuming SNRs evolve in a warm phase dominated ISM, we find dust destruction timescales are increased by a factor of ~2 compared to those of Jones et al. (1996), who assumed a hot gas dominated ISM. Recent estimates of supernova rates and ISM mass lead to another factor of ~3 increase in the destruction timescales, resulting in a silicate grain destruction timescale of ~2-3 Gyr. These increases, while not able resolve the problem of the discrepant timescales for silicate grain destruction and creation, are an important step towards understanding the origin, and evolution of dust in the ISM.

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

The nature, composition, abundance, and size distribution of dust in galaxies is determined by the rate at which it is created in the different stellar sources and destroyed by interstellar shocks. Because of their extensive wavelength coverage, proximity, and nearly face-on geometry, the Magellanic Clouds (MCs) provide a unique opportunity to study these processes in great detail. In this paper we use the complete sample of supernova remnants (SNRs) in the MCs to calculate the lifetime and destruction efficiencies of silicate and carbon dust in these galaxies. We find dust lifetimes of 22 +- 13 Myr (30 +- 17 Myr) for silicate (carbon) grains in the LMC, and 54 +- 32 Myr (72 +- 43 Myr) for silicate (carbon) grains in the SMC. The significantly shorter lifetimes in the MCs, as compared to the Milky Way, are explained as the combined effect of their lower total dust mass, and the fact that the dust-destroying isolated SNe in the MCs seem to be preferentially occurring in regions with higher than average dust-to-gas (D2G) mass ratios. We also calculate the supernova rate and the current star formation rate in the MCs, and use them to derive maximum dust injection rates by asymptotic giant branch stars and core collapse supernovae. We find that the injection rates are an order of magnitude lower than the dust destruction rates by the SNRs. This supports the conclusion that, unless the dust destruction rates have been considerably overestimated, most of the dust must be reconstituted from surviving grains in dense molecular clouds. More generally, we also discuss the dependence of the dust destruction rate on the local D2G mass ratio, the ambient gas density and metallicity, as well as the application of our results to other galaxies and dust evolution models.

NGC 6872, hereafter the Condor, is a large spiral galaxy that is interacting with its closest companion, the S0 galaxy IC 4970. The extent of the Condor provides an opportunity for detailed investigation of the impact of the interaction on the current star formation rate and its history across the galaxy, on the age and spatial distribution of its stellar population, and on the mechanism that drive the star formation activity. To address these issues we analyzed the far-ultraviolet (FUV) to near-infrared (near-IR) spectral energy distribution (SED) of 17, 10 kpc diameter, regions across the galaxy, and derived their star formation history, current star formation rate, and stellar population and mass. We find that most of the star formation takes place in the extended arms, with very little star formation in the central 5 kpc of the galaxy, in contrast to what was predicted from previous numerical simulations. There is a trend of increasing star formation activity with distance from the nucleus of the galaxy, and no evidence for a recent increase in the current star formation rate (SFR) due to the interaction. The nucleus itself shows no significant current star formation activity. The extent of the Condor also provides an opportunity to test the applicability of a single standard prescription for conversion of the FUV + IR (22 micron) intensities to a star formation rate for all regions. We find that the conversion factor differs from region to region, arising from regional differences in the stellar populations.

The origin of dust in galaxies is still a mystery. The majority of the refractory elements are produced in supernova explosions but it is unclear how and where dust grains condense and grow, and how they avoid destruction in the harsh environments of star-forming galaxies. The recent detection of 0.1-0.5 solar masses of dust in nearby supernova remnants suggests in situ dust formation, while other observations reveal very little dust in supernovae the first few years after explosion. Observations of the bright SN 2010jl have been interpreted as pre-existing dust, dust formation or no dust at all. Here we report the rapid (40-240 days) formation of dust in its dense circumstellar medium. The wavelength dependent extinction of this dust reveals the presence of very large (> 1 micron) grains, which are resistant to destructive processes. At later times (500-900 days), the near-IR thermal emission shows an accelerated growth in dust mass, marking the transition of the supernova from a circumstellar- to an ejecta-dominated source of dust. This provides the link between the early and late dust mass evolution in supernovae with dense circumstellar media.

We present deep continuum observations using the GISMO camera at a wavelength of 2 mm centered on the Hubble Deep Field (HDF) in the GOODS-N field. These are the first deep field observations ever obtained at this wavelength. The 1 sigma sensitivity in the innermost approx. 4 arcminutes of the 7 utes map is approx. 135 uJy/beam, a factor of three higher in flux/beam sensitivity than the deepest available SCUBA 850 um observations, and almost a factor of four higher in flux / beam sensitivity than the combined MAMBO/AzTEC 1.2 mm observations of this region. Our source extraction algorithm identifies 12 sources directly, and another 3 through correlation with known sources at 1.2 mm and 850 um. Five of the directly detected GISMO sources have counterparts in the MAMBO/AzTEC catalog, and four of those also have SCUBA counterparts. HDF850.1, one of the first blank-field detected submillimeter galaxies, is now detected at 2 mm. The median redshift of all sources with counterparts of known redshifts is med(z) = 2.91 +/- 0.94. Statistically, the detections are most likely real for 5 of the seven 2 mm sources without shorter wavelength counterparts, while the probability for none of them being real is negligible.

The evolution of dust at redshifts z>9, and consequently the dust properties, differs greatly from that in the local universe. In contrast to the local universe, core collapse supernovae (CCSNe) are the only source of thermally-condensed dust. Because of the low initial dust-to-gas mass ratio, grain destruction rates are low, so that CCSNe are net producers of interstellar dust. Galaxies with large initial gas mass or high mass infall rate will therefore have a more rapid net rate of dust production comported to galaxies with lower gas mass, even at the same star formation rate. The dust composition is dominated by silicates, which exhibit a strong rise in the UV opacity near the Lyman break. This "silicate-UV break" may be confused with the Lyman break, resulting in a misidentification of a galaxies' photometric redshift. In this paper we demonstrate these effects by analyzing the spectral energy distribution (SED) of MACS1149-JD, a lensed galaxy at z=9.6. A potential 2mm counterpart of MACS1149-JD has been identified with GISMO. While additional observations are required to corroborate this identification, we use this possible association to illustrate the physical processes and the observational effects of dust in the very high redshift universe.

Infrared continuum observations provide a means of investigating the physical composition of the dust in the ejecta and swept up medium of the Cas A supernova remnant. Using low resolution Spitzer IRS spectra (5-35 $\mu$m), and broad-band Herschel PACS imaging (70, 100, and 160 $\mu$m), we identify characteristic dust spectra, associated with ejecta layers that underwent distinct nuclear burning histories. The most luminous spectrum exhibits strong emission features at $\sim9$ and 21 $\mu$m and is closely associated with ejecta knots with strong Ar emission lines. The dust features can be reproduced by magnesium silicate grains with relatively low Mg to Si ratios. Another dust spectrum is associated with ejecta having strong Ne emission lines. It has no indication of any silicate features, and is best fit by Al$_2$O$_3$ dust. A third characteristic dust spectrum shows features that are best matched by magnesium silicates with a relatively high Mg to Si ratio. This dust is primarily associated with the X-ray emitting shocked ejecta, but it is also evident in regions where shocked interstellar or circumstellar material is expected. However, the identification of dust composition is not unique, and each spectrum includes an additional featureless dust component of unknown composition. Colder dust of indeterminate composition is associated with emission from the interior of the SNR, where the reverse shock has not yet swept up and heated the ejecta. Most of the dust mass in Cas A is associated with this unidentified cold component, which is $\lesssim0.1$ $M_{\odot}$. The mass of warmer dust is only $\sim 0.04$ $M_{\odot}$.

Supernova (SN) explosions are crucial engines driving the evolution of galaxies by shock heating gas, increasing the metallicity, creating dust, and accelerating energetic particles. In 2012 we used the Atacama Large Millimeter/Submillimeter Array to observe SN 1987A, one of the best-observed supernovae since the invention of the telescope. We present spatially resolved images at 450um, 870um, 1.4mm, and 2.8mm, an important transition wavelength range. Longer wavelength emission is dominated by synchrotron radiation from shock-accelerated particles, shorter wavelengths by emission from the largest mass of dust measured in a supernova remnant (>0.2Msun). For the first time we show unambiguously that this dust has formed in the inner ejecta (the cold remnants of the exploded star's core). The dust emission is concentrated to the center of the remnant, so the dust has not yet been affected by the shocks. If a significant fraction survives, and if SN 1987A is typical, supernovae are important cosmological dust producers.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large amounts of newly-condensed dust in their ejecta, comparable to the total mass of available refractory elements. The dust masses derived from these observations assume that all the grains of a given species radiate at the same temperature, regardless of the dust heating mechanism or grain radius. In this paper, we derive the dust mass in the ejecta of the Crab Nebula, using a physical model for the heating and radiation from the dust. We adopt a power-law distribution of grain sizes and two different dust compositions (silicates and amorphous carbon), and calculate the heating rate of each dust grain by the radiation from the pulsar wind nebula (PWN). We find that the grains attain a continuous range of temperatures, depending on their size and composition. The total mass derived from the best-fit models to the observed IR spectrum is 0.019-0.13 solar masses, depending on the assumed grain composition. We find that the power-law size distribution of dust grains is characterized by a power-law index of 3.5-4.0 and a maximum grain size larger than 0.1 microns. The grain sizes and composition are consistent with what is expected for dust grains formed in a Type IIP SN. Our derived dust mass is at least a factor of two less than the mass reported in previous studies of the Crab Nebula that assumed more simplified two-temperature models. The results of this study show that a physical model resulting in a realistic distribution of dust temperatures can constrain the dust properties and affect the derived dust masses. Our study may also have important implications for deriving grain properties and mass estimates in other SNRs and for the ultimate question of whether SNe are major sources of dust in the Galactic interstellar medium (ISM) and in external galaxies.

We report on the present stage of SN 1987A as observed by the Chandra X-ray Observatory. We reanalyze published Chandra observations and add three more epochs of Chandra data to get a consistent picture of the evolution of the X-ray fluxes in several energy bands. We discuss the implications of several calibration issues for Chandra data. Using the most recent Chandra calibration files, we find that the 0.5-2.0 keV band fluxes of SN 1987A have increased by ~6 x 10 ^-13 erg s^-1 cm^-2 per year since 2009. This is in contrast with our previous result that the 0.5-2.0 keV light curve showed a sudden flattening in 2009. Based on our new analysis, we conclude that the forward shock is still in full interaction with the equatorial ring.

The extragalactic background light (EBL) is one of the fundamental observational quantities in cosmology. All energy releases from resolved and unresolved extragalactic sources, and the light from any truly diffuse background, excluding the cosmic microwave background (CMB), contribute to its intensity and spectral energy distribution. It therefore plays a crucial role in cosmological tests for the formation and evolution of stellar objects and galaxies, and for setting limits on exotic energy releases in the universe. The EBL also plays an important role in the propagation of very high energy gamma-rays which are attenuated en route to Earth by pair producing gamma-gamma interactions with the EBL and CMB. The EBL affects the spectrum of the sources, predominantly blazars, in the ~10 GeV to 10 TeV energy regime. Knowledge of the EBL intensity and spectrum will allow the determination of the intrinsic blazar spectrum in a crucial energy regime that can be used to test particle acceleration mechanisms and VHE gamma-ray production models. Conversely, knowledge of the intrinsic gamma-ray spectrum and the detection of blazars at increasingly higher redshifts will set strong limits on the EBL and its evolution. This paper reviews the latest developments in the determination of the EBL and its impact on the current understanding of the origin and production mechanisms of gamma-rays in blazars, and on energy releases in the universe. The review concludes with a summary and future directions in Cherenkov Telescope Array techniques and in infrared ground-based and space observatories that will greatly improve our knowledge of the EBL and the origin and production of very high energy gamma-rays.

Recent infrared (IR) observations of freshly-formed dust in supernova remnants (SNRs) have yielded significantly lower dust masses than predicted by theoretical models and measured from high redshift observations. The Crab Nebula's pulsar wind is thought to be sweeping up freshly-formed supernova (SN) dust along with the ejected gas. The evidence for this dust was found in the form of an IR excess in the integrated spectrum of the Crab and in extinction against the synchrotron nebula that revealed the presence of dust in the filament cores. We present the first spatially resolved emission spectra of dust in the Crab Nebula acquired with the Infrared Spectrograph aboard the Spitzer Space Telescope. The IR spectra are dominated by synchrotron emission and show forbidden line emission from from S, Si, Ne, Ar, O, Fe, and Ni. We derived a synchrotron spectral map from the 3.6 and 4.5 microns images, and subtracted this contribution from our data to produce a map of the residual continuum emission from dust. The dust emission appears to be concentrated along the ejecta filaments and is well described by an amorphous carbon or silicate grain compositions. We find a dust temperature of 55+/- 4 K for silicates and 60 +/- 7 K for carbon grains. The total estimated dust mass is 0.0012-0.012 solar masses, well below the theoretical dust yield predicted for a core-collapse supernova. Our grain heating model implies that the dust grain radii are relatively small, unlike what is expected for dust grains formed in a Type IIP SN.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of ~ 1.5 keV, or with two thermal components with temperatures of 1.5 and 0.2 keV. Previous studies suggest that the hot component may originate from reverse-shocked SN ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from supernova (SN) ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of ~ 140 K by a relatively dense, hot plasma, that also gives rise to the hot X-ray emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-ray emitting gas. The total mass of the warm dust component is at least 0.013 solar masses, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission.

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 present the late-time optical light curve of the ejecta of SN 1987A measured from HST imaging observations spanning the past 17 years. We find that the flux from the ejecta declined up to around year 2001, powered by the radioactive decay of 44Ti. Then the flux started to increase, more than doubling by the end of 2009. We show that the increase is the result of energy deposited by X-rays produced in the interaction with the circumstellar medium. We suggest that the change of the dominant energy input to the ejecta, from internal to external, marks the transition from supernova to supernova remnant. The details of the observations and the modelling are described in the accompanying supplementary information.

Analyses of high-redshift ultraluminous infrared (IR) galaxies traditionally use the observed optical to submillimeter spectral energy distribution (SED) and estimates of the dynamical mass as observational constraints to derive the star formation rate (SFR), the stellar mass, and age of these objects. An important observational constraint neglected in the analysis is the mass of dust giving rise to the IR emission. In this paper we add this constraint to the analysis of AzTEC-3. Adopting an upper limit to the mass of stars and a bolometric luminosity for this object, we construct stellar and chemical evolutionary scenarios, constrained to produce the inferred dust mass and observed luminosity before the associated stellar mass exceeds the observational limit. We find that the model with a Top Heavy IMF provided the most plausible scenario consistent with the observational constraints. In this scenario the dust formed over a period of ~200 Myr, with a SFR of ~500 Msun/yr. These values for the age and SFR in AzTEC-3 are significantly higher and lower, respectively, from those derived without the dust mass constraint. However, this scenario is not unique, and others cannot be completely ruled out because of the prevailing uncertainties in the age of the galaxy, its bolometric luminosity, and its stellar and dust masses. A robust result of our models is that all scenarios require most of the radiating dust mass to have been accreted in molecular clouds. Our new procedure highlights the importance of a multiwavelength approach, and of the use of dust evolution models in constraining the age and the star formation activity and history in galaxies.

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission to measure the gravity-wave signature of primordial inflation through its distinctive imprint on the linear polarization of the cosmic microwave background. The instrument consists of a polarizing Michelson interferometer configured as a nulling polarimeter to measure the difference spectrum between orthogonal linear polarizations from two co-aligned beams. Either input can view the sky or a temperature-controlled absolute reference blackbody calibrator. PIXIE will map the absolute intensity and linear polarization (Stokes I, Q, and U parameters) over the full sky in 400 spectral channels spanning 2.5 decades in frequency from 30 GHz to 6 THz (1 cm to 50 um wavelength). Multi-moded optics provide background-limited sensitivity using only 4 detectors, while the highly symmetric design and multiple signal modulations provide robust rejection of potential systematic errors. The principal science goal is the detection and characterization of linear polarization from an inflationary epoch in the early universe, with tensor-to-scalar ratio r < 10^{-3} at 5 standard deviations. The rich PIXIE data set will also constrain physical processes ranging from Big Bang cosmology to the nature of the first stars to physical conditions within the interstellar medium of the Galaxy.

Direct measurements of the extragalactic background light (EBL) in the near-IR to mid-IR waveband are extremely difficult due to an overwhelming foreground from the zodiacal light that outshines the faint cosmological diffuse radiation field by more than an order of magnitude. Indirect constraints on the EBL are provided by gamma-ray observations of AGN. Using the combination of the Fermi Gamma-Ray Space Telescope together with the current generation of ground-based air Cherenkov telescopes (H.E.S.S., MAGIC, and VERITAS) provides unprecedented sensitivity and spectral coverage for constraining the EBL in the near- to mid-IR. In this paper we present new limits on the EBL based on the analysis of the broad-band spectra of a select set of gamma-ray blazars covering 200 MeV to several TeV. The EBL intensity at 15 microns is constrained to be 1.36 +/- 0.58 nW m^-2 sr^-1. We find that the fast evolution and baseline EBL models of Stecker et al. (2006), as well as the model of Kneiske et al. (2004), predict significantly higher EBL intensities in the mid-IR (15 microns) than is allowed by the constraints derived here. In addition, the model of Franceschini et al. (2008) and the fiducial model of Dominguez et al. (2011) predict near- to mid-IR ratios smaller than that predicted by our analysis. Namely, their intensities in the near-IR are too low while their intensities in the mid-IR are marginally too high. All of the aforementioned models are inconsistent with our analysis at the >3 sigma level.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

Two distinct scenarios for the origin of the ~ 4 \times 10^8 M\odot of dust observed in the high-redshift (z = 6.4) quasar J1148+5251 have been proposed. The first assumes that this galaxy is much younger than the age of the universe at that epoch so that only supernovae (SNe) could have produced this dust. The second scenario assumes a significantly older galactic age, so that the dust could have formed in lower-mass asymptotic giant branch (AGB) stars. Presenting new integral solutions for the chemical evolution of metals and dust in galaxies, we offer a critical evaluation of these two scenarios, and observational consequences that can discriminate between the two. We show that AGB stars can produce the inferred mass of dust in this object, however, the final mass of surviving dust depends on the galaxy's star formation history (SFH). In general supernovae cannot produce the observed amount of dust unless the average SN event creates over ~ 1 M\odot of dust in its ejecta. However, special SFHs can be constructed in which SNe can produce the inferred dust mass with a reasonable average dust yield of ~ 0.15 M\odot. The two scenarios propose different origins for the galaxy's spectral energy distribution, different star formation efficiencies and stellar masses, and consequently different comoving number densities of J1148+5251-type hyperluminous infrared (IR) objects. The detection of diagnostic mid-IR fine structure lines, and more complete surveys determining the comoving number density of these objects can discriminate between the two scenarios.

This paper presents late-time near-infrared and {\it Spitzer} mid-infrared photometric and spectroscopic observations of warm dust in the Type IIn SN 2005ip in NGC 2906. The spectra show evidence for two dust components with different temperatures. Spanning the peak of the thermal emission, these observations provide strong constraints on the dust mass, temperature, and luminosity, which serve as critical diagnostics for disentangling the origin and heating mechanism of each component. The results suggest the warmer dust has a mass of $\sim 5 \times 10^{-4}~$\msolar, originates from newly formed dust in the ejecta, or possibly the cool, dense shell, and is continuously heated by the circumstellar interaction. By contrast, the cooler component likely originates from a circumstellar shock echo that forms from the heating of a large, pre-existing dust shell $\sim 0.01 - 0.05$~\msolar~by the late-time circumstellar interaction. The progenitor wind velocity derived from the blue edge of the He 1 1.083 \micron~P Cygni profile indicates a progenitor eruption likely formed this dust shell $\sim$100 years prior to the supernova explosion, which is consistent with a Luminous Blue Variable (LBV) progenitor star.

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &\mu&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

We have used the Spitzer satellite to monitor the mid-IR evolution of SN 1987A over a 5 year period spanning the epochs between days 6000 and 8000 since the explosion. The supernova (SN) has evolved into a supernova remnant (SNR) and its radiative output is dominated by the interaction of the SN blast wave with the pre-existing equatorial ring (ER). The mid-IR spectrum is dominated by emission from ~180 K silicate dust, collisionally-heated by the hot X-ray emitting gas with a temperature and density of ~5x10^6 K and 3x10^4 cm-3, respectively. The mass of the radiating dust is ~1.2x10^(-6) Msun on day 7554, and scales linearly with IR flux. The infrared to soft-X-ray flux ratio is roughly constant with a value of 2.5. Gas-grain collisions therefore dominate the cooling of the shocked gas. The constancy of of this ratio suggests that very little grain processing or gas cooling have occurred throughout this epoch. The shape of the dust spectrum remained unchanged during the observations while the total flux increased with a time dependence of t^(0.87), t being the time since the first encounter between the blast wave and the ER. These observations are consistent with the transitioning of the blast wave from free expansion to a Sedov phase as it propagates into the main body of the ER.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

Observations with the Hubble Space Telescope (HST), conducted since 1990, now offer an unprecedented glimpse into fast astrophysical shocks in the young remnant of supernova 1987A. Comparing observations taken in 2010 using the refurbished instruments on HST with data taken in 2004, just before the Space Telescope Imaging Spectrograph failed, we find that the Ly-a and H-a lines from shock emission continue to brighten, while their maximum velocities continue to decrease. We observe broad blueshifted Ly-a, which we attribute to resonant scattering of photons emitted from hotspots on the equatorial ring. We also detect NV~\lambda\lambda 1239,1243 A line emission, but only to the red of Ly-A. The profiles of the NV lines differ markedly from that of H-a, suggesting that the N^{4+} ions are scattered and accelerated by turbulent electromagnetic fields that isotropize the ions in the collisionless shock.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

The standard method of mapping the interstellar medium in a galaxy, by observing the molecular gas in the CO 1-0 line and the atomic gas in the 21-cm line, is largely limited with current telescopes to galaxies in the nearby universe. In this letter, we use SPIRE observations of the galaxies M99 and M100 to explore the alternative approach of mapping the interstellar medium using the continuum emission from the dust. We have compared the methods by measuring the relationship between the star-formation rate and the surface density of gas in the galaxies. We find the two methods give relationships with a similar dispersion, confirming that observing the continuum emission from the dust is a promising method of mapping the interstellar medium in galaxies.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ 0.35 solar) nearby dwarf galaxy, NGC 1705, in six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub-mm by LABOCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a significant cold dust component within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 350 microns, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an additional cold dust component.

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 {\mu}m) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

We present the first Herschel PACS and SPIRE images of the low-metallicity galaxy NGC6822 observed from 70 to 500 mu and clearly resolve the HII regions with PACS and SPIRE. We find that the ratio 250/500 is dependent on the 24 mu surface brightness in NGC6822, which would locally link the heating processes of the coldest phases of dust in the ISM to the star formation activity. We model the SEDs of some regions HII regions and less active regions across the galaxy and find that the SEDs of HII regions show warmer ranges of dust temperatures. We derive very high dust masses when graphite is used in our model to describe carbon dust. Using amorphous carbon, instead, requires less dust mass to account for submm emission due to its lower emissivity properties. This indicates that SED models including Herschel constraints may require different dust properties than commonly used.

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

We present Herschel-SPIRE observations at 250-500um of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust lane extending north-south; and a bright emission feature 10kpc to the south-east. We estimate that approximately 10^6 solar masses of dust is spatially coincident with atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC4438 due to recent tidal interactions with M86. The gas-to-dust ratio of the cold gas component ranges from ~20-80. We discuss the different heating mechanisms for the dust features.

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

We present infrared colours (in the 25-500 mic spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated bynthe synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f350/f500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter beta<2 in low metallicity, active systems.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

We present a detailed analysis of the far-IR properties of the bright, lensed, z = 2.3, SMG, SMM J2135-0102, using new observations with Herschel, SCUBA-2 and the VLA. These data allow us to constrain the galaxy's SED and show that it has an intrinsic rest-frame 8-1000um luminosity, L(bol), of (2.3 +/- 0.2) x 10^12 L(sun) and a likely SFR of ~400 M(sun)/yr. The galaxy sits on the far-IR/radio correlation for far-IR-selected galaxies. At ~>70um, the SED can be described adequately by dust components with T(d) ~ 30 and 60K. Using SPIRE's Fourier Transform Spectrometer we report a detection of the [CII] 158um cooling line. If the [CII], CO and far-IR continuum arise in photo-dissociation regions, we derive a characteristic gas density, n ~ 10^3 cm^-3, and a far-UV radiation field, G_0, 10^3x stronger than the Milky Way. L([CII])/L(bol) is significantly higher than in local ULIRGs but similar to the values found in local star-forming galaxies and starburst nuclei. This is consistent with SMM J2135-0102 being powered by starburst clumps distributed across ~2 kpc, evidence that SMGs are not simply scaled-up ULIRGs. Our results show that SPIRE's FTS has the ability to measure the redshifts of distant, obscured galaxies via the blind detection of atomic cooling lines, but it will not be competitive with ground-based CO-line searches. It will, however, allow detailed study of the integrated properties of high-redshift galaxies, as well as the chemistry of their ISM, once more suitably bright candidates have been found.

We study the formation of molecular precursors to dust in the ejecta of Population III supernovae using a chemical kinetic approach. Our work focuses on zero-metallicity 20 Msun and 170 Msun progenitors, and we consider fully-macroscopically mixed and unmixed ejecta. The nucleation stage for small silica, metal oxides and sulphides, pure metal, and carbon clusters is described with a new chemical reaction network. We consider the effect of the pressure dependence of critical nucleation rates, and test the impact of microscopically-mixed He+ on carbon dust formation. The unmixed ejecta of a 170 Msun progenitor supernova synthesizes ~ 5.6 Msun of small clusters, while its 20 Msun counterpart produces ~ 0.103 Msun. Our results point to smaller amounts of dust formed in the ejecta of Pop. III supernovae by a factor ~ 5 compared to values derived by previous studies, and to different dust chemical composition. Such deviations result from some erroneous assumptions made, the inappropriate use of classical nucleation theory to model dust formation, and the omission of the synthethis of molecules in supernova ejecta. Unmixed ejecta of massive Pop. III supernovae chiefly form silica and/or silicates, and pure silicon grains whereas their lower mass counterparts form a dust mixture dominated by silica and/or silicates, pure silicon and iron sulphides. Amorphous carbon can only condense in ejecta where the carbon-rich zone is deprived of He+. The first dust enrichment to the primordial gas in the early universe from Pop. III massive supernova comprises primarily pure silicon, silica and silicates. If carbon dust is present at redshift z> 6, alternative dust sources must be considered.

The Herschel Reference Survey is a guaranteed time Herschel key project and will be a benchmark study of dust in the nearby universe. The survey will complement a number of other Herschel key projects including large cosmological surveys that trace dust in the distant universe. We will use Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350 and 500 micron. The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux limits in the K-band to minimize the selection effects associated with dust and with young high-mass stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological types (ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass and environment, (ii) the connection between the dust content and composition and the other phases of the interstellar medium and (iii) the origin and evolution of dust in galaxies. In this paper, we describe the goals of the survey, the details of the sample and some of the auxiliary observing programs that we have started to collect complementary data. We also use the available multi-frequency data to carry out an analysis of the statistical properties of the sample.

We study the formation and destruction of molecules in the ejecta of Population III supernovae (SNe) using a chemical kinetic approach to follow the evolution of molecular abundances from day 100 to day 1000 after explosion. The chemical species included range from simple di-atomic molecules to more complex dust precursor species. All relevant chemical processes that are unique to the SN environment are considered. Our work focuses on zero-metallicity progenitors with masses of 20, 170, and 270 Msun, and we study the effect of different levels of heavy element mixing and the inward diffusion of hydrogen on the ejecta chemistry. We show that the ejecta chemistry does not reach a steady state within the relevant time-span for molecule formation. The primary species formed are O2, CO, SiS, and SO. The SiO, formed as early as 200 days after explosion, is rapidly depleted by the formation of silica molecular precursors in the ejecta. The rapid conversion of CO to C2 and its thermal fractionation at temperatures above 5000 K allow for the formation of carbon chains in the oxygen-rich zone of the unmixed models, providing an important pathway for the formation of carbon dust in hot environments where the C/O ratio is less than 1. We show that the fully-mixed ejecta of a 170 Msun progenitor synthesizes 11.3 Mun of molecules whereas 20 Msun and 270 Msun progenitors produce 0.78, and 3.2 Msun of molecules, respectively. The admixing of 10 % of hydrogen into the fully-mixed ejecta of the 170 Msun progenitor increases its molecular yield to ~ 47Msun. The unmixed ejecta of a 170 Msun progenitor supernova without hydrogen penetration synthesizes ~37 Msun of molecules, whereas its 20 Msun counterpart produces ~ 1.2 Msun. Finally, we discuss the cosmological implication of molecule formation by Pop. III SNe in the early universe.

The Diffuse Supernova Neutrino Background (DSNB) provides an immediate opportunity to study the emission of MeV thermal neutrinos from core-collapse supernovae. The DSNB is a powerful probe of stellar and neutrino physics, provided that the core-collapse rate is large enough and that its uncertainty is small enough. To assess the important physics enabled by the DSNB, we start with the cosmic star formation history of Hopkins & Beacom (2006) and confirm its normalization and evolution by cross-checks with the supernova rate, extragalactic background light, and stellar mass density. We find a sufficient core-collapse rate with small uncertainties that translate into a variation of +/- 40% in the DSNB event spectrum. Considering thermal neutrino spectra with effective temperatures between 4-6 MeV, the predicted DSNB is within a factor 4-2 below the upper limit obtained by Super-Kamiokande in 2003. Furthermore, detection prospects would be dramatically improved with a gadolinium-enhanced Super-Kamiokande: the backgrounds would be significantly reduced, the fluxes and uncertainties converge at the lower threshold energy, and the predicted event rate is 1.2-5.6 events /yr in the energy range 10-26 MeV. These results demonstrate the imminent detection of the DSNB by Super-Kamiokande and its exciting prospects for studying stellar and neutrino physics.

We present recent advances in theoretical studies of the formation and evolution of dust in primordial supernovae (SNe) that are considered to be the main sources of dust in the early universe. Being combined with the results of calculations of dust formation in the ejecta of Population III SNe, the investigations of the evolution of newly formed dust within supernova remnants (SNRs) show that smaller grains are predominantly destroyed by sputtering in the shocked gas, while larger grains are injected into the ambient medium. The mass of dust grains surviving the destruction in SNRs reaches up to 0.1--15 $M_\odot$, which is high enough to account for the content of dust observed for the host galaxies of quasars at $z > 5$. In addition, the transport of dust formed in the ejecta causes the formation of low-mass stars in the dense shell of primordial SNRs and affects the elemental composition of those stars. We also show that the flat extinction curve is expected in the high-redshift universe where SNe are the possible sources of dust.

Models for the evolution of dust are used to show that the observed trend of the abundance of polycyclic aromatic hydrocarbons (PAHs) with metallicity is the result of the delayed injection of carbon dust that formed in low mass asymptotic giant branch (AGB) stars into the interstellar medium. We also use our dust evolution models to examine the origin of dust at redshifts > 6, when only supernovae and their remnants could have been, respectively, their sources of production and destruction. Unless an average supernova (or its progenitor) produces between 0.1 and 1 Msun of dust, alternative sources will need to be invoked to account for the massive amount of dust observed at these redshifts.

The abundances of gas and dust (solids and complex molecules) in the interstellar medium (ISM) as well as their composition and structures impact practically all of astrophysics. Fundamental processes from star formation to stellar winds to galaxy formation all scale with the number of metals. However, significant uncertainties remain in both absolute and relative abundances, as well as how these vary with environment, e.g., stellar photospheres versus the interstellar medium (ISM). While UV, optical, IR, and radio studies have considerably advanced our understanding of ISM gas and dust, they cannot provide uniform results over the entire range of column densities needed. In contrast, X-rays will penetrate gas and dust in the cold (3K) to hot (100,000,000K) Universe over a wide range of column densities (log NH=20-24 cm^-2), imprinting spectral signatures that reflect the individual atoms which make up the gas, molecule or solid. *X-rays therefore are a powerful and viable resource for delving into a relatively unexplored regime for determining gas abundances and dust properties such as composition, charge state, structure, and quantity via absorption studies, and distribution via scattering halos.*

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

This is a report on the findings of the extragalactic science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses extragalactic science topics including active galactic nuclei, cosmic ray acceleration in galaxies, galaxy clusters and large scale structure formation shocks, and the study of the extragalactic infrared and optical background radiation. The scientific potential of ground based gamma-ray observations of Gamma-Ray Bursts and dark matter annihilation radiation is covered in other sections of the white paper.

We combine the recent estimate of the contribution of galaxies to the 3.6 micron intensity of the extragalactic background light (EBL) with optical and near-infrared (IR) galaxy counts to set new limits on intrinsic spectra of some of the most distant TeV blazars 1ES 0229+200, 1ES 1218+30.4, and 1ES 1101-232, located at redshifts 0.1396, 0.182, and 0.186, respectively. The new lower limit on the 3.6 micron EBL intensity is significantly higher than the previous one set by the cumulative emission from resolved Spitzer galaxies. Correcting for attenuation by the revised EBL, we show that the differential spectral index of the intrinsic spectrum of the three blazars is 1.28 +- 0.20 or harder. These results present blazar emission models with the challenge of producing extremely hard intrinsic spectra in the sub-TeV to multi-TeV regime. These results also question the reliability of recently derived upper limits on the near-IR EBL intensity that are solely based on the assumption that intrinsic blazar spectra should not be harder than 1.5.

We investigate the evolution of dust formed in Population III supernovae (SNe) by considering its transport and processing by sputtering within the SN remnants (SNRs). We find that the fates of dust grains within SNRs heavily depend on their initial radii $a_{\rm ini}$. For Type II SNRs expanding into the ambient medium with density of $n_{\rm H,0} = 1$ cm$^{-3}$, grains of $a_{\rm ini} < 0.05$ $\mu$m are detained in the shocked hot gas and are completely destroyed, while grains of $a_{\rm ini} > 0.2$ $\mu$m are injected into the surrounding medium without being destroyed significantly. Grains with $a_{\rm ini}$ = 0.05-0.2 $\mu$m are finally trapped in the dense shell behind the forward shock. We show that the grains piled up in the dense shell enrich the gas up to 10$^{-6}-10^{-4}$ $Z_\odot$, high enough to form low-mass stars with 0.1-1 $M_\odot$. In addition, [Fe/H] in the dense shell ranges from -6 to -4.5, which is in good agreement with the ultra-metal-poor stars with [Fe/H] < -4. We suggest that newly formed dust in a Population III SN can have great impacts on the stellar mass and elemental composition of Population II.5 stars formed in the shell of the SNR.

(Condensed form) - The serendipitous discovery of infrared echoes around the supernova remnant of Cas A by the Spitzer satellite has provided astronomers with a unique opportunity to study the properties of the echoing material and the history and nature of the outburst that generated these echoes. All the echoes located within a distance of ~15 arcmin from the SN are caused by the delayed arrival of thermal emission from dust located at a distance of 160 lyr (corresponding to half the adopted age of the remnant) directly behind the origin of the explosion. The spectra of the echoes are distinct from that of the general diffuse interstellar medium (ISM) revealing hot silicate grains that are either stochastically heated to temperatures in excess of ~150 K, or radiating at an equilibrium temperature of this value. We show that the optical light curve from the supernova, is not capable of producing such spectra, and could therefore not have given rise to the echoes. Instead, we find that the echoes were generated by an intense and short burst of EUV-UV radiation with a luminosity of ~ 1.5E11 Lsun. The average H-column density of the IR emitting region in the echoing clouds is about 5E17 cm-2. Taking a burst time of ~1 d gives a cloud density of ~200 cm-3, typical of dense IR cirrus.

Multiwavelength observations of supernova remnant (SNR) 1987A show that its morphology and luminosity are rapidly changing at X-ray, optical, infrared, and radio wavelengths as the blast wave from the explosion expands into the circumstellar equatorial ring, produced by mass loss from the progenitor star. The observed infrared (IR) radiation arises from the interaction of dust grains that formed in mass outflow with the soft X-ray emitting plasma component of the shocked gas. Spitzer IRS spectra at 5 - 30 microns taken on day 6190 since the explosion show that the emission arises from ~ 1.1E-6 Msun of silicate grains radiating at a temperature of ~180+20-15 K. Subsequent observations on day 7137 show that the IR flux had increased by a factor of 2 while maintaining an almost identical spectral shape. The observed IR-to-X-ray flux ratio (IRX) is consistent with that of a dusty plasma with standard Large Magellanic Cloud dust abundances. IRX has decreased by a factor of ~ 2 between days 6190 and 7137, providing the first direct observation of the ongoing destruction of dust in an expanding SN blast wave on dynamic time scales. Detailed models consistent with the observed dust temperature, the ionization timescale of the soft X-ray emission component, and the evolution of IRX suggest that the radiating silicate grains are immersed in a 3.5E6 K plasma with a density of (0.3-1)E4 cm^{-3}, and have a size distribution that is confined to a narrow range of radii between 0.023 and 0.22 microns. Smaller grains may have been evaporated by the initial UV flash from the supernova.

Dusty hyperluminous galaxies in the early universe provide unique environments for studying the role of massive stars in the formation and destruction of dust. At redshifts above 6, when the universe was less than 1 Gyr old, dust could have only condensed in the explosive ejecta of Type II supernovae (SNe), since most of the progenitors of the AGB stars, the major alternative source of interstellar dust, did not have time to evolve off the main sequence. We present analytical models for the evolution of the gas, dust, and metals in high redshift galaxies, with a special application to SDSS J1148+5251, a hyperluminous quasar at z = 6.4. We show that an average SN must condense at least 1 Msun of dust to account for the mass of dust in this object, when grain destruction by supernova remnants (SNRs) is taken into account. This required yield is in excess of ~0.05 Msun, the largest mass of dust inferred from infrared observations of Cas A. If the yield of Cas A is typical, then other processes, such as accretion onto preexisting grains in molecular clouds is needed to produce the mass of dust in J1148+5251. For such process to be effective, SNR must significantly increase, presumably by non-evaporative grain-grain collisions during the late stages of their evolution, the number of nucleation centers onto which refractory elements can condense in molecular clouds.

Spectral and photometric observations of nearby galaxies show a correlation between the strength of their mid-IR aromatic features, attributed to PAH molecules, and their metal abundance, leading to a deficiency of these features in low-metallicity galaxies. In this paper, we suggest that the observed correlation represents a trend of PAH abundance with galactic age, reflecting the delayed injection of carbon dust into the ISM by AGB stars in the final post-AGB phase of their evolution. AGB stars are the primary sources of PAHs and carbon dust in galaxies, and recycle their ejecta back to the interstellar medium only after a few hundred million years of evolution on the main sequence. In contrast, more massive stars that explode as Type II supernovae inject their metals and dust almost instantaneously after their formation. We first determined the PAH abundance in galaxies by constructing detailed models of UV-to-radio SED of galaxies that estimate the contribution of dust in PAH-free HII regions, and PAHs and dust from photodissociation regions, to the IR emission. All model components: the galaxies' stellar content, properties of their HII regions, and their ionizing and non-ionizing radiation fields and dust abundances, are constrained by their observed multiwavelength spectrum. After determining the PAH and dust abundances in 35 nearby galaxies using our SED model, we use a chemical evolution model to show that the delayed injection of carbon dust by AGB stars provides a natural explanation to the dependence of the PAH content in galaxies with metallicity. We also show that larger dust particles giving rise to the far-IR emission follow a distinct evolutionary trend closely related to the injection of dust by massive stars into the ISM.

We investigate the evolution of dust that formed at Population III supernova (SN) explosions and its processing through the collisions with the reverse shocks resulting from the interaction of the SN ejecta with the ambient medium. In particular, we investigate the transport of the shocked dust within the SN remnant (SNR), and its effect on the chemical composition, the size distribution, and the total mass of dust surviving in SNR. We find that the evolution of the reverse shock, and hence its effect on the processing of the dust depends on the thickness of the envelope retained by the progenitor star. Furthermore, the transport and survival of the dust grains depend on their initial radius, a_{ini}, and composition: For Type II SNRs expanding into the interstellar medium (ISM) with a density of n_{H,0}=1 cm^{-3}, small grains with a_{ini} < ~ 0.05 micron are completely destroyed by sputtering in the postshock flow, while grains with a_{ini}= 0.05--0.2 micron are trapped into the dense shell behind the forward shock. Very large grains of a_{ini} > ~ 0.2 micron are ejected into the ISM without decreasing their sizes significantly. We find that the total mass fraction of dust that is destroyed by the reverse shock ranges from 0.2 to 1.0, depending on the energy of the explosion and the density of the ambient ISM. The results of our calculations have significant impact on the abundance pattern of subsequent generation of stars that form in the dense shell of primordial SNRs.

Determination of the cosmic infrared background (CIB) at far infrared wavelengths using COBE/DIRBE data is limited by the accuracy to which foreground interplanetary and Galactic dust emission can be modeled and subtracted. Previous determinations of the far infrared CIB (e.g., Hauser et al. 1998) were based on the detection of residual isotropic emission in skymaps from which the emission from interplanetary dust and the neutral interstellar medium were removed. In this paper we use the Wisconsin H-alpha Mapper (WHAM) Northern Sky Survey as a tracer of the ionized medium to examine the effect of this foreground component on determination of the CIB. We decompose the DIRBE far infrared data for five high Galactic latitude regions into H I and H-alpha correlated components and a residual component. We find the H-alpha correlated component to be consistent with zero for each region, and we find that addition of an H-alpha correlated component in modeling the foreground emission has negligible effect on derived CIB results. Our CIB detections and 2 sigma upper limits are essentially the same as those derived by Hauser et al. and are given by nu I_nu (nW m-2 sr-1) < 75, < 32, 25 +- 8, and 13 +- 3 at 60, 100, 140, and 240 microns, respectively. Our residuals have not been subjected to a detailed anisotropy test, so our CIB results do not supersede those of Hauser et al. We derive upper limits on the 100 micron emissivity of the ionized medium that are typically about 40% of the 100 micron emissivity of the neutral atomic medium. This low value may be caused in part by a lower dust-to-gas mass ratio in the ionized medium than in the neutral medium, and in part by a shortcoming of using H-alpha intensity as a tracer of far infrared emission.

Dusty hyperluminous galaxies in the early universe provide unique environments for studying the role of massive stars in the formation and destruction of dust. At redshifts above ~ 6, when the universe was less than ~ 1 Gyr old, dust could have only condensed in the explosive ejecta of Type II supernovae (SNe), since most of the progenitors of the AGB stars, the major alternative source of interstellar dust, did not have time to evolve off the main sequence since the onset of star formation. In this paper we present analytical models for the evolution of the gas, dust, and metals in high redshift galaxies, with a special application to SDSS J1148+5251, a hyperluminous quasar at z = 6.4. We find that an average supernova must condense at least 1 Msun of dust to account for the observed dust mass in this quasar. Observationally, it is in excess of the largest dust yield of ~0.02 Msun found thus far in the ejecta of any SN. If future observations find this to be a typical supernova dust yield, then additional processes, such as accretion onto preexisting grains, or condensation around the AGN will need to be invoked to account for the large amount of dust in this and similar objects. The galaxy's star formation history is still uncertain, and current observations of the gas, metal, and dust contents of J1148 can be reproduced by either an intensive and short burst of star formation (~ 1000 Msun/yr) with a duration of ~ 100 Myr, or a much lower star formation rate (~ 100 Msun/yr) occurring over the lifetime of the galaxy.

Multiwavelength observations of SNR 1987A show that its morphology is rapidly changing at X-ray, radio, and optical wavelengths as the blast wave from the explosion expands into the circumstellar equatorial ring. Infrared emission arises from the interaction of dust grains with the hot X-ray emitting gas. We show that the IR emission provides important complementary information on the interaction of the SN blast wave with the circumstellar equatorial ring that cannot be obtained at any other wavelength.

We have imaged a $\sim$6 arcminute$^2$ region in the Bo\"otes Deep Field using the 350 $\mu$m-optimised second generation Submillimeter High Angular Resolution Camera (SHARC II), achieving a peak 1$\sigma$ sensitivity of $\sim$5 mJy. We detect three sources above 3$\sigma$, and determine a spurious source detection rate of 1.09 in our maps. In the absence of $5\sigma$ detections, we rely on deep 24 $\mu$m and 20 cm imaging to deduce which sources are most likely to be genuine, giving two real sources. From this we derive an integral source count of 0.84$^{+1.39}_{-0.61}$ sources arcmin$^{-2}$ at $S>13$ mJy, which is consistent with 350 $\mu$m source count models that have an IR-luminous galaxy population evolving with redshift. We use these constraints to consider the future for ground-based short-submillimetre surveys.

We present high resolution 11.7 and 18.3um mid-IR images of SN 1987A obtained on day 6526 with T-ReCS attached to the Gemini telescope. The 11.7um flux has increased significantly since our last observations on day 6067. The images clearly show that all the emission arises from the equatorial ring (ER). Spectra obtained with Spitzer, on day 6184 with MIPS at 24um, on day 6130 with IRAC in the 3.6-8um region, and on day 6190 with IRS in the 12-37um region show that the emission consists of thermal emission from silicate dust that condensed out in the red giant wind of the progenitor star. The dust temperature is ~166K, and the emitting dust mass is ~2.6 x 10-6 Msun. Lines of [Ne II]12.82um and [Ne III]15.56um are clearly present, as well as a weak [Si II]34.8um line. We also detect two lines near 26um which we tentatively ascribe to [Fe II]25.99um and [O IV]25.91um. Comparison of the Gemini 11.7um image with X-ray images from Chandra, UV-optical images from HST, and radio synchrotron images obtained by the ATCA show generally good correlation of the images across all wavelengths. Because of the limited resolution of the mid-IR images we cannot uniquely determine the location or heating mechanism of the dust giving rise to the emission. The dust could be collisionally heated by the X-ray emitting plasma, providing a unique diagnostic of plasma conditions. Alternatively, the dust could be radiatively heated in the dense UV-optical knots that are overrun by the advancing supernova blast wave. In either case the dust-to-gas mass ratio in the circumstellar medium around the supernova is significantly lower than that in the general ISM of the LMC, suggesting either a low condensation efficiency in the wind of the progenitor star, or the efficient destruction of the dust by the SN blast wave.

We report the detection of a 3.6sigma 350micron-selected source in the Bootes Deep Field. The source, the first Short-wavelength Submillimeter-selected Galaxy (SSG 1), was discovered as part of a blank field extragalactic survey using the 350micron-optimised Submillimeter High Angular Resolution Camera (SHARC II) at the Caltech Submillimeter Observatory. With multiwavelength photometry from NOAO-NDWFS (R and I band), FLAMEX (J and K_s), Spitzer (IRAC and MIPS) and the Westerbork 1.4GHz Deep Survey (radio upper limit), we are able to constrain the photometric redshift using different methods, all of which suggest a redshift of close to 1. In the absence of long-wavelength submillimeter data we use SED templates to infer that this source is an ultraluminous infrared galaxy (ULIRG) with a dust temperature of 30+/-5 K, occupying a region of luminosity-temperature space shared by modarate redshift ISO-selected ULIRGs (rather than high redshift SCUBA-selected SMGs). SHARC II can thus select SMGs with moderately "warm'' dust that might be missed in submillimeter surveys at longer wavelengths.

The intensity of the diffuse ~ 1 - 4 micron sky emission from which solar system and Galactic foregrounds have been subtracted is in excess of that expected from energy released by galaxies and stars that formed during the z < 5 redshift interval (Arendt & Dwek 2003, Matsumoto et al. 2005). The spectral signature of this excess near-infrared background light (NIRBL) component is almost identical to that of reflected sunlight from the interplanetary dust cloud, and could therefore be the result of the incomplete subtraction of this foreground emission component from the diffuse sky maps. Alternatively, this emission component could be extragalactic. Its spectral signature is consistent with that of redshifted continuum and recombination line emission from HII regions formed by the first generation of very massive stars. In this paper we analyze the implications of this spectral component for the formation rate of these Population III stars, the redshift interval during which they formed, the reionization of the universe and evolution of collapsed halo masses. We find that to reproduce the intensity and spectral shape of the NIRBL requires a peak star formation rate that is higher by about a factor of 4 to 10 compared to those derived from hierarchical models. Furthermore, an extragalactic origin for the NIRBL leads to physically unrealistic absorption-corrected spectra of distant TeV blazars. All these results suggest that Pop III stars contribute only a fraction of the NIRBL intensity with zodiacal light, star forming galaxies, and/or non-nuclear sources giving rise to the remaining fraction.

The 1 - 5 micron diffuse sky emission from which local foreground emission from the solar system and the Galaxy have been subtracted exceeds the brightness that can be attributed to normal star forming galaxies. The nature of this excess near-infrared background light (NIRBL) is still controversial. On one hand, it has been interpreted as a distinct spectral feature created by the redshifted emission from primordial (Population III) stars that have formed at redshifts > 8. On the other hand, the NIRBL spectrum is almost identical to that of the zodiacal cloud, raising the possibility that it is of local origin. Blazars can, in principle, offer a simple test for the nature and origin of the NIRBL. Very high energy gamma-ray photons emitted by these objects are attenuated on route to earth by photon-photon interactions with the extragalactic background light (EBL). This paper examines whether the extragalactic nature of the NIRBL can be determined from the analysis of the TeV spectra of blazars. (Abridged)

This paper reviews our current understanding of interstellar dust models, what constitutes a viable dust model, what observational constraints are essential for deriving such model, and the current viable dust models. Interstellar dust exhibits spatial and temporal variations, and the paper reviews the ingredients in constructing models for the evolution of dust, stressing the current uncertainties in the yield of dust from supernovae and AGB stars. Finally, the paper briefly describes the effect of dust evolution on the spectral energy distribution of galaxies, and, using a very simple criterion, present a rough estimate of the redshift when galaxies first become opaque.

Very high energy (~ TeV) $\gamma$-rays from blazars are attenuated by photons from the extragalactic background light (EBL). Observations of blazars can therefore provide an ideal opportunity for determining the EBL intensity if their intrinsic spectrum is known. Conversely, knowledge of the EBL intensity can be used to determine the intrinsic blazar spectrum. Unfortunately, neither the EBL intensity nor the intrinsic blazar spectrum is known with high enough precision to accurately derive one quantity from the other. In this paper we use the most recent data on the EBL to construct twelve different realizations representing all possible permutations between EBL limits and the detections in the different wavelength regions. We use these realizations to explore the effects of the EBL on the inferred spectra of blazars. In particular, we show that the frequently cited "IR background-TeV gamma-ray crisis" does not exist, and derive the intrinsic spectra and peak energies of the blazars Mrk 421, 501 and H1426+428 for EBL realizations that give rise to physically viable intrinsic blazar spectra. We also show that the intrinsic spectrum of Mrk~421 during a period of intense flaring activity has a peak energy that seems to shift to higher energies at higher flux states. Finally, we also explore the effect of the uncertainties in the absolute calibration of the gamma-ray energies on derived TeV opacities and the intrinsic blazar spectra.

The extinction curve derived from ISO mid-infrared (IR) observations of the Galactic center (GC) exhibits a surprisingly flat behavior in the ~ 3 to 8 micron region, contrary to the deep minimum expected from standard interstellar dust models consisting of bare silicate and graphite dust particles. We show that this extinction is likely caused by the presence of metallic needles in the interstellar medium (ISM) towards the Galactic center. If the needles contribute only to the 3 - 8 micron extinction, they must have a long wavelength cutoff at ~ 8 microns, and therefore a typical length over radius ratio of about 600, smaller than the 3000 aspect ratio determined for the needles in Cas A. Homogeneously distributed throughout the ISM, they comprise only a minor mass fraction of the ISM, with a needle-to-H mass ratio of \~5x10^{-6}, which is equivalent to 0.14% of the silicate dust mass. Their total ISM abundance can then be readily explained by the combined production in SNe and O-rich stellar outflows. It is currently unclear how universal the GC extinction law is. Local 2 - 5 micron extinction measurements seem to be consistent with the standard extinction law, suggesting a non uniform distribution of needles in the ISM. The GC observations show that metallic needles, in spite of their low abundance or non uniform distribution, can be the dominant source of opacity in the 3 - 8 micron wavelength region. However, expelled into the intergalactic medium, their abundance is too low to cause any dimming of cosmological sources, and their length is too short to make them a significant source of submillimeter emission.

Recently, Dunne et al. (2003) obtained 450 and 850 micron SCUBA images of CasA, and reported the detection of 2-4 M_sun of cold, 18K, dust in the remnant. Here we show that their interpretation of the observations faces serious difficulties. Their inferred dust mass is larger than the mass of refractory material in the ejecta of a 10 to 30 M_sun star. The cold dust model faces even more difficulties if the 170 micron observations of the remnant are included in the analysis, decreasing the cold dust temperature to ~ 8K, and increasing its mass to > 20 M_sun. We offer here a more plausible interpretation of their observation, in which the cold dust emission is generated by conducting needles with properties that are completely determined by the combined submillimeter and X-ray observations of the remnant. The needles consist of metallic whiskers with <1% of embedded impurities that may have condensed out of blobs of material that were expelled at high velocities from the inner metal-rich layers of the star in an asymmetric explosion. The needles are collisionally heated by the shocked gas to a temperature of 8K. Taking the destruction of needles into account, a dust mass of only 1E-4 to 1E-3M_sun is needed to account for the observed SCUBA emission. Aligned in the magnetic field, needles may give rise to observable polarized emission. The detection of submillimeter polarization will therefore offer definitive proof for a needle origin for the cold dust emission. Supernovae may still be proven to be important sources of interstellar dust, but the evidence is still inconclusive.

We present new interstellar dust models which have been derived by simultaneously fitting the far-ultraviolet to near-infrared extinction, the diffuse infrared (IR) emission and, unlike previous models, the elemental abundance constraints on the dust for different interstellar medium abundances, including solar, F and G star, and B star abundances. The fitting problem is a typical ill-posed inversion problem, in which the grain size distribution is the unknown, which we solve by using the method of regularization. The dust model contains various components: PAHs, bare silicate, graphite, and amorphous carbon particles, as well as composite particles containing silicate, organic refractory material, water ice, and voids. The optical properties of these components were calculated using physical optical constants. As a special case, we reproduce the Li & Draine (2001) results, however their model requires an excessive amount of silicon, magnesium, and iron to be locked up in dust: about 50 ppm (atoms per million of H atoms), significantly more than the upper limit imposed by solar abundances of these elements, about 34, 35, and 28 ppm, respectively. A major conclusion of this paper is that there is no unique interstellar dust model that simultaneously fits the observed extinction, diffuse IR emission, and abundances constraints.

A viable interstellar dust model - characterized by the composition, morphology, and size distribution of the dust grains and by the abundance of the different elements locked up in the dust - should fit all observational constraints arising primarily from the interactions of the dust with incident radiation or the ambient gas. In spite of the many different manifestations of these interactions, we still lack a comprehensive dust model that is consistent with {\it all} the observational constraints. An important advance towards the construction of such a model was recently made by Zubko, Dwek, and Arendt (2003, ZDA03) who, for the first time, included the average interstellar extinction, the diffuse infrared emission, {\it and} the interstellar abundances as explicit constraints in the construction of models consisting of astrophysical dust particles without any predetermined functional form for their grain size distribution. The results showed the existence of many distinct dust models that satisfy a basic set of observational constraints. X-ray halos, produced primarily by small angle scattering off large dust particles along the line of sight to bright X-ray sources, probe dust properties largely inaccessible at other wavelengths. In this contribution we briefly review the microscopic physics and the macroscopic effects that determine the intensity and spatial profile of X-ray halos. Focusing on the X-ray halo around the bright source GX~13+1, we show that halo observations can be used to discriminate between the different, currently viable, ZDA03 dust models. X-ray halos may therefore be an essential constraint that need to be explicitly included in all future dust models.

Energy spectra of gamma-ray blazars may contain an imprint from the cosmic infrared background radiation due to gamma-ray absorption (pair-production) by soft photons constituting the extragalactic background light (EBL). The signature of this imprint depends on the spectral shape of the EBL. In this work we correct the observed spectra of Mrk 421 and Mrk 501 for absorption using different possible realizations of the EBL, consistent with the most recent detections and limits. We present the intrinsic gamma-ray spectrum of these sources for the different EBL scenarios. These spectra reveal their true peak energy and luminosities, which provide important information on the nature and physical characteristics of the particle acceleration mechanism operating in these sources.

We decompose the COBE/DIRBE observations of the near-IR sky brightness (minus zodiacal light) into Galactic stellar and interstellar medium (ISM) components and an extragalactic background. This empirical procedure allows us to estimate the 4.9 micron cosmic infrared background (CIB) as a function of the CIB intensity at shorter wavelengths. A weak indication of a rising CIB intensity at wavelengths > 3.5 micron hints at interesting astrophysics in the CIB spectrum, or warns that the foreground zodiacal emission may be incompletely subtracted. Subtraction of only the stellar component from the zodiacal-light-subtracted all-sky map reveals the clearest 3.5 micron ISM emission map, which is found to be tightly correlated with the ISM emission at far-IR wavelengths.

We use the radio-infrared (IR) flux correlation bet ween star-forming galaxies in the local universe to derive a simple analytical expression between the intensity of the IR background and the brightness t emperature of the radio background. This relation i s insensitive to the star formation history of the galaxies that produce the cosmic IR background (CIB). We use the observed CIB intensity to constrain t he cosmic star formation history, and the relation between the CIB and the cosmic radio background (CR B) to constrain the relative contribution of star-f orming galaxies to the CRB. Current limits on the C IB intensity predict a 178 MHz brightness temperature of ~18+-9 K, about half of the 37+-8 K inferred for an isotropic radio component. This suggests that s tar-forming galaxies and AGN contribute about equal ly to the CRB intensity at that frequency.

The cosmic infrared background records much of the radiant energy released by processes of structure formation that have occurred since the decoupling of matter and radiation following the Big Bang. In the past few years, data from the Cosmic Background Explorer mission provided the first measurements of this background, with additional constraints coming from studies of the attenuation of TeV gamma-rays. At the same time there has been rapid progress in resolving a significant fraction of this background with the deep galaxy counts at infrared wavelengths from the Infrared Space Observatory instruments and at submillimeter wavelengths from the Submillimeter Common User Bolometer Array instrument. This article reviews the measurements of the infrared background and sources contributing to it, and discusses the implications for past and present cosmic processes.

The extragalactic background light (EBL), exclusive of the cosmic microwave background, consists of the cumulative radiative output from all energy sources in the universe since the epoch of recombination. Most of this energy is released at ultraviolet and optical wavelengths. However, observations show that a significant fraction of the EBL falls in the 10 to 1000 micron wavelength regime. This provides conclusive evidence that we live in a dusty universe, since only dust can efficiently absorbs a significant fraction of the background energy and reemit it at infrared wavelengths. The general role of dust in forming the cosmic infrared background (CIB) is therefore obvious. However, its role in determining the exact spectral shape of the CIB is quite complex. The CIB spectrum depends on the microscopic physical properties of the dust, its composition, abundance, and spatial distribution relative to the emitting sources, and its response to evolutionary processes that can modify all the factors listed above. This paper will present a brief summary of the many ways dust affects the intensity and spectral shape of the cosmic infrared background. In an Appendix we present new limits on the mid-infrared intensity of the CIB using TeV gamma-ray observations of Mrk 501.

Recent observations with dust detectors on board the interplanetary spacecraft Ulysses and Galileo have recorded a substantial flux of large interstellar grains with radii between 0.25 and 2.0 mu entering the solar system from the local interstellar cloud. The most commonly used interstellar grain size distribution is characterized by a a^-3.5 power law in grain radii a, and extends to a maximum grain radius of 0.25 mu. The extension of the interstellar grain size distribution to such large radii will have a major effect on the median grain size, and on the amount of mass needed to be tied up in dust for a given visual optical depth. It is therefore important to investigate whether this population of larger dust particles prevails in the general interstellar medium, or if it is merely a local phenomenon. The presence of large interstellar grains can be mainly inferred from their effect on the intensity and radial profiles of scattering halos around X-ray sources. In this paper we examine the grain size distribution that gives rise to the X-ray halo around Nova Cygni 1992. The results of our study confirm the need to extend the interstellar grain size distribution in the direction of this source to and possibly beyond 2.0 mu. The model that gives the best fit to the halo data is characterized by: (1) a grain size distribution that follows an a^-3.5 power law up to 0.50 mu, followed by an a^-4.0 extension from 0.50 mu to 2.0 mu; and (2) silicate and graphite (carbon) dust-to-gas mass ratios of 0.0044 and 0.0022, respectively, consistent with solar abundances constraints. Additional observations of X-ray halos probing other spatial directions are badly needed to test the general validity of this result.

SAFIRE is a versatile imaging Fabry-Perot spectrograph covering 145 to 655 microns, with spectral resolving powers ranging over 5-10,000. Selected as a "PI" instrument for the airborne Stratospheric Observatory for Infrared Astronomy (SOFIA), SAFIRE will apply two-dimensional pop-up bolometer arrays to provide background-limited imaging spectrometry. Superconducting transition edge bolometers and SQUID multiplexers are being developed for these detectors. SAFIRE is expected to be a "First Light" instrument, useable during the initial SOFIA operations. Although a PI instrument rather than a "Facility Class" science instrument, it will be highly integrated with the standard SOFIA planning, observation, and data analysis tools.

We present analytical approximations for calculating the scattering, absorption and escape of nonionizing photons from a spherically symmetric two-phase clumpy medium, with either a central point source of isotropic radiation, a uniform distribution of isotropic emitters, or uniformly illuminated by external sources. The analytical approximations are based on the mega-grains model of two-phase clumpy media, as proposed by Hobson & Padman, combined with escape and absorption probability formulae for homogeneous media. The accuracy of the approximations is examined by comparison with 3D Monte Carlo simulations of radiative transfer, including multiple scattering. Our studies show that the combined mega-grains and escape/absorption probability formulae provide a good approximation of the escaping and absorbed radiation fractions for a wide range of parameters characterizing the medium. A realistic test is performed by modeling the absorption of a starlike source of radiation by interstellar dust in a clumpy medium, and by calculating the resulting equilibrium dust temperatures and infrared emission spectrum of both the clumps and the interclump medium. In particular, we find that the temperature of dust in clumps is lower than in the interclump medium if clumps are optically thick. Comparison with Monte Carlo simulations of radiative transfer in the same environment shows that the analytic model yields a good approximation of dust temperatures and the emerging UV to FIR spectrum of radiation for all three types of source distributions mentioned above. Our analytical model provides a numerically expedient way to estimate radiative transfer in a variety of interstellar conditions and can be applied to a wide range of astrophysical environments, from star forming regions to starburst galaxies.

A Monte Carlo model of radiative transfer in multi-phase dusty media is applied to the situation of stars and clumpy dust in a sphere or a disk. The distribution of escaping and absorbed photons are shown for various filling factors and densities. Analytical methods of approximating the escaping fraction of radiation, based on the Mega-Grains approach, are discussed. Comparison with the Monte Carlo results shows that the escape probability formulae provide a reasonable approximation of the escaping/absorbed fractions, for a wide range of parameters characterizing a clumpy dusty medium. A possibly more realistic model of the interstellar medium is one in which clouds have a self-similar hierarchical structure of denser and denser clumps within clumps, resulting in a fractal distribution of gas and dust. Monte Carlo simulations of radiative transfer in such multi-phase fractal media are compared with the two-phase clumpy case.

Foreground emission and scattered light from interplanetary dust (IPD) particles and emission from Galactic stellar sources are the greatest obstacles for determining the cosmic infrared background (CIB) from diffuse sky measurements in the ~ 1 to 5 micron range. We describe a new method for deriving the CIB at near infrared wavelengths which reduces the uncertainties associated with the removal of the Galactic stellar emission component from the sky maps. The method produces positive residuals at 3.5 and 4.9 microns, of which only the 3.5 micron residual is nearly isotropic. We consider our result as a tentative detection of the CIB at this wavelength.

Spectroscopic observations at 2.4 - 45 microns of the young supernova remnant Cas A with the Infrared Space Observatory (ISO) Short Wavelength Spectrometer (SWS) reveal strong emission lines of O, Ne, Si, S, and Ar. These lines are observed at high velocities (several thousand km/s), and are therefore associated with the supernova ejecta known as the fast-moving knots (FMKs). Continuum emission from dust is also seen in the Cas A spectrum. The continuum strength is spatially well correlated with the O and Ar line strengths, indicating that the dust emission also arises from the FMKs. The dust continuum has an emission feature at ~22 microns which cannot be fit by typical astronomical silicates, but can be fit with a particular class of silicate minerals. This suggests that the dust in Cas A is silicate material that has freshly condensed from the Cas A ejecta into a mineral form that is uncharacteristic of typical ISM dust grains.

A major goal of modern astrophysics is to understand the processes by which the universe evolved from its initial simplicity, as seen in measurements of the Cosmic Microwave Background, to the universe we see today, with complexity on all scales. While the diffuse background measurements of COBE reveal the importance of the far infrared and submillimeter in early galaxy and star formation, the understanding of the development of complex structure requires high resolution imaging and spectroscopy. We present a concept for a space mission called SPECS, the Submillimeter Probe of the Evolution of Cosmic Structure, which provides these capabilities. SPECS is a cold, spatial and spectral Michelson interferometer with adjustable baselines ranging up to 1 km. It has Hubble sensitivity and angular resolution in the far IR and submillimeter, spectral resolution up to 10,000, and a 14' field of view. SPECS will be able to image thermal dust continuum and infrared cooling and diagnostic line emission over a wide range of redshifts, providing extinction-free astrophysical probes of young galaxies and early cosmic structures, and measures of the luminosity and heavy element formation history of the universe. SPECS would also have the potential to improve vastly our knowledge of protostars, protoplanetary systems, Active Galactic Nuclei and other objects in the local universe. We recommend that a concerted effort be made during the next decade to develop and test certain critical technologies (photon-counting far IR detectors, formation flying spacecraft, cold, lightweight mirrors, and active coolers), so that SPECS can be deployed in about 2015.

Foreground emission and scattered light from interplanetary dust (IPD) particles and emission from Galactic stellar sources are the greatest obstacles for determining the cosmic infrared background (CIB) from diffuse sky measurements in the ~ 1 to 5 micron range. We use ground-based observational limits on the K-band intensity of the CIB in conjunction with skymaps obtained by the Diffuse Infrared Background Experiment (DIRBE) on the COBE satellite to reexamine the limits on the CIB at 1.25, 3.5, and 4.9 microns. Adopting a CIB intensity of 7.4 nW m-2 sr-1 at 2.2 microns, and using the 2.2 micron DIRBE skymap from which the emission from IPD cloud has been subtracted, we create a spatial template of the Galactic stellar contribution to the diffuse infrared sky. This template is then used to subtract the contribution of the diffuse Galactic stellar emission from the IPD-emission-subtracted DIRBE skymaps. The DIRBE 100 micron data are used to estimate the small contribution of emission from interstellar dust at 3.5 and 4.9 microns. Our method significantly reduces the errors associated with the subtraction of Galactic starlight, leaving only the IPD emission component as the primary obstacle for the detection of the CIB at these wavelengths. This analysis leads to a tentative detection of the CIB at 3.5 microns. The cosmological implications of these results are discussed in the paper.

The COBE Diffuse Infrared Background Experiment (DIRBE) was designed to search for the cosmic infrared background (CIB) radiation. Scattered light and thermal emission from the interplanetary dust (IPD) are major contributors to the diffuse sky brightness at most infrared wavelengths. Accurate removal of this zodiacal light foreground is a necessary step toward a direct measurement of the CIB. The zodiacal light foreground contribution in each of the 10 DIRBE wavelength bands ranging from 1.25 to 240 microns is distinguished by its apparent seasonal variation over the whole sky. This contribution has been extracted by fitting the brightness calculated from a parameterized physical model to the time variation of the all-sky DIRBE measurements over 10 months of observations. The model brightness is evaluated as the integral along the line of sight of the product of a source function and a three-dimensional dust density distribution function. The dust density distribution is composed of multiple components: a smooth cloud, three asteroidal dust bands, and a circumsolar ring near 1 A.U. By using a directly measurable quantity which relates only to the IPD cloud, we exclude other contributors to the sky brightness from the IPD model. Using the IPD model described here, high-quality maps of the infrared sky with the zodiacal foreground removed have been generated. Imperfections in the model reveal themselves as low-level systematic artifacts in the residual maps which correlate with components of the IPD. The most evident of these artifacts are located near the ecliptic plane in the mid-infrared, and are less than 2% of the zodiacal foreground brightness. Uncertainties associated with the model are discussed, including implications for the CIB search.

The DIRBE on the COBE spacecraft was designed primarily to conduct systematic search for an isotropic CIB in ten photometric bands from 1.25 to 240 microns. The results of that search are presented here. Conservative limits on the CIB are obtained from the minimum observed brightness in all-sky maps at each wavelength, with the faintest limits in the DIRBE spectral range being at 3.5 microns (\nu I_\nu < 64 nW/m^2/sr, 95% CL) and at 240 microns (\nu I_\nu < 28 nW/m^2/sr, 95% CL). The bright foregrounds from interplanetary dust scattering and emission, stars, and interstellar dust emission are the principal impediments to the DIRBE measurements of the CIB. These foregrounds have been modeled and removed from the sky maps. Assessment of the random and systematic uncertainties in the residuals and tests for isotropy show that only the 140 and 240 microns data provide candidate detections of the CIB. The residuals and their uncertainties provide CIB upper limits more restrictive than the dark sky limits at wavelengths from 1.25 to 100 microns. No plausible solar system or Galactic source of the observed 140 and 240 microns residuals can be identified, leading to the conclusion that the CIB has been detected at levels of \nu I_\nu = 25+-7 and 14+-3 nW/m^2/sr at 140 and 240 microns respectively. The integrated energy from 140 to 240 microns, 10.3 nW/m^2/sr, is about twice the integrated optical light from the galaxies in the Hubble Deep Field, suggesting that star formation might have been heavily enshrouded by dust at high redshift. The detections and upper limits reported here provide new constraints on models of the history of energy-releasing processes and dust production since the decoupling of the cosmic microwave background from matter.

The Cosmic Infrared Background (CIB) is hidden behind veils of foreground emission from our own solar system and Galaxy. This paper describes procedures for removing the Galactic IR emission from the 1.25 - 240 micron COBE DIRBE maps as steps toward the ultimate goal of detecting the CIB. The Galactic emission models are carefully chosen and constructed so that the isotropic CIB is completely retained in the residual sky maps. We start with DIRBE data from which the scattered light and thermal emission of the interplanetary dust (IPD) cloud have already been removed. Locations affected by the emission from bright compact and stellar sources are excluded from the analysis. The unresolved emission of faint stars at near- and mid-IR wavelengths is represented by a model based on Galactic source counts. The 100 micron DIRBE observations are used as the spatial template for the interstellar medium (ISM) emission at high latitudes. Correlation of the 100 micron data with H I column density allows us to isolate the component of the observed emission that is associated with the ISM. Limits are established on the far-IR emissivity of the diffuse ionized medium, which indicate a lower emissivity per H nucleus than in the neutral medium. At 240 micron, we find that adding a second spatial template to the ISM model can greatly improve the accuracy of the model at low latitudes. The crucial product of this analysis is a set of all-sky IR maps from which the Galactic (and IPD) emission has been removed. We discuss systematic uncertainties and potential errors in the foreground subtraction process that may have an impact on studies seeking to detect the CIB in the residual maps.

In this paper we examine the cosmological constraints of the recent DIRBE and FIRAS detection of the extragalactic background light between 125-5000 microns on the metal and star formation histories of the universe.

The COBE FIRAS data contain foreground emission from interplanetary, Galactic interstellar dust and extragalactic background emission. We use three different methods to separate the various emission components, and derive the spectrum of the extragalactic Far InfraRed Background (FIRB). Each method relies on a different set of assumptions, which affect the FIRB spectrum in different ways. Despite this, the FIRB spectra derived by these different methods are remarkably similar. The average spectrum that we derive in the \nu = 5 - 80 cm^{-1} (2000-125 \um) frequency interval is: I(\nu) = (1.3 \pm 0.4)\times10^{-5} (\nu /\nu_0)^{0.64 \pm 0.12} P_{\nu}(18.5 \pm 1.2 K), where \nu_0=100 cm^{-1} (\lambda_0=100 \um) and P is the Planck function. The derived FIRB spectrum is consistent with the DIRBE 140 and 240 \um detections. The total intensity received in the 5 - 80 cm^{-1} frequency interval is 14 nW m^{-2} sr^{-1}, and comprises about 20% of the total intensity expected from the energy release from nucleosynthesis throughout the history of the universe.

Small angle scatterings of X-rays by interstellar dust particles create halos around x-ray sources. By observing these halos with sufficient angular resolution it is possible to derive information about both the spatial distribution of the dust along the line of sight to the source and some of the physical properties of the scattering dust particles. Previous analyses of halos have been used to infer that the scattering dust particles must be fluffy. These fluffy dust particles have also been suggested as a possible solution to the interstellar carbon "crisis." This crisis results from the discrepancy between the available ISM carbon abundance and the amount required to be in dust to account for the average ISM extinction. Dwek (1997), however, pointed out that fluffy dust particles may not solve the carbon crisis. In light of this, we re-examine the evidence for fluffy dust in the ISM from x-ray scattering data. We find that previous studies overestimated the intensity of the halos by using the RG approximation, which fails for energies below 1 keV. Using the exact Mie solution we find that below 1 keV the scattering becomes much less efficient. Applying the Mie calculations to the x-ray halo observed around Nova Cygni 1992, we conclude that the data can be explained by scattering from regular dust, and does not require interstellar dust particles to be fluffy.

We present models for the evolution of the elemental abundances in the gas and dust phases of the interstellar medium (ISM) of our Galaxy by generalizing standard models for its dynamical and chemical evolution. In these models, the stellar birthrate history is determined by the infall rate of primordial gas, and by its functional dependence on the mass surface density of the stars and gas. We adopt a two component model for the Galaxy, consisting of a central bulge and an exponential disk with different infall rates and stellar birthrate histories. Condensation in stellar winds, Type Ia and Type II supernovae, and the accretion of refractory elements onto preexisting grains in dense molecular clouds are the dominant contributors to the abundance of elements locked up in the dust. Grain destruction by sputtering and evaporative grain-grain collisions in supernova remnants are the most important mechanisms that return these elements back to the gas phase. We calculate the dust production rate by the various dust sources, analyze the origin of the elemental depletion pattern, and study the relation between dust abundance and ISM metallicity, and the evolution of the the dust abundance and composition at each Galactocentric radius as a function of time. The derived relation of dust mass with metallicity is compared to the observed Galactic dust abundance gradient, and to the M$_{dust}$ versus log(O/H) relation that is observed in external Dwarf galaxies. The dependence of dust composition on the mass of the progenitor star, and the delayed recycling of newly synthesized dust by low mass stars back to the ISM give rise to variations in the dust composition as a function of time. Our models provide a framework for the self-consistent inclusion of dust in

Interstellar dust models are facing a "carbon crisis", so called because recent observations suggest that the abundance of carbon available for dust in the interstellar medium is less than half of the amount required to be tied up in graphite grains in order to explain the interstellar extinction curve. This paper presents an detailed assessment of a newly-proposed dust model (Mathis 1996), in which the majority of the interstellar carbon is contained in composite and fluffy dust (CFD) grains. Per unit mass, these grains produce more UV extinction, and can therefore account for the interstellar extinction curve with about half the carbon required in traditional dust models. The results of our analysis show that the CFD model falls short in solving the carbon crisis, in providing a fit to the UV-optical interstellar extinction curve. It also predicts a far-infrared emissivity in excess of that observed with the COBE/DIRBE and FIRAS instruments from the diffuse interstellar medium. The failure of the new model highlights the interrelationships between the various dust properties and their observational consequences, and the need to satisfy them all simultaneously in any comprehensive interstellar dust model. In light of these problems, the paper examines other possible solutions to the carbon crisis.

Using data obtained by the DIRBE and FIRAS instruments on the COBE spacecraft we present the mean 3.5-1000 um dust spectrum from the high latitude cirrus clouds. These data represent the most comprehensive wavelength coverage of dust emission, and first strong evidence for the presence of PAHs in cirrus. The COBE data are well fit with a dust model consisting of a mixture of PAH molecules and bare astronomical silicate and graphite grains. From the model we derive the size distribution, abundances relative to the total hydrogen column density, and relative contribution of each dust component to the observed IR emission. The model provides a good fit to the FIRAS spectrum in the 140-500 um wavelength regime, but leaves an excess Galactic emission component at 500-1000 um. The nature of this component is still unresolved. The model requires about 85% of the solar carbon abundance to be locked up in dust, a third of that in PAHs. The remaining ~ 15% is in the gas as inferred from the [C II] 158 um line emission detected by the FIRAS. PAHs emit about 25% of the total IR power, and must therefore be responsible for the same fraction of the total UV-optical extinction.

We derive Galactic continuum spectra from 5-96/cm from COBE/FIRAS observations. The spectra are dominated by warm dust emission, which may be fit with a single temperature in the range 16-21 K (for nu^2 emissivity) along each line of sight. Dust heated by the attenuated radiation field in molecular clouds gives rise to intermediate temperature (10-14 K) emission in the inner Galaxy only. A widespread, very cold component (4-7 K) with optical depth that is spatially correlated with the warm component is also detected. The cold component is unlikely to be due to very cold dust shielded from starlight, because it is present at high latitude. We consider hypotheses that the cold component is due to enhanced submillimeter emissivity of the dust that gives rise to the warm component, or that it may be due to very small, large, or fractal particles. Lack of substantial power above the emission from warm dust places strong constraints on the amount of cold gas in the Galaxy. The microwave sky brightness due to interstellar dust is dominated by the cold component, and its angular variation could limit our ability to discern primordial fluctuations in the cosmic microwave background radiation.