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

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 study the coherence of the near-infrared and X-ray background fluctuations and the X-ray spectral properties of the sources producing it. We use data from multiple Spitzer and Chandra surveys, including the UDS/SXDF surveys, the Hubble Deep Field North, the EGS/AEGIS field, the Chandra Deep Field South and the COSMOS surveys, comprising $\sim$2275 Spitzer/IRAC hours and $\sim$~16 Ms of Chandra data collected over a total area of $\sim$~1~deg$^2$. We report an overall $\sim$5$\sigma$ detection of a cross-power signal on large angular scales $>$ 20$''$ between the 3.6 and 4.5\mum\ and the X-ray bands, with the IR vs [1-2] keV signal detected at 5.2$\sigma$. The [0.5-1] and [2-4] keV bands are correlated with the infrared wavelengths at a $\sim$1$-$3$\sigma$ significance level. The hardest X-ray band ([4-7] keV) alone is not significantly correlated with any infrared wavelengths due to poor photon and sampling statistics. We study the X-ray SED of the cross-power signal. We find that its shape is consistent with a variety of source populations of accreting compact objects, such as local unabsorbed AGNs or high-z absorbed sources. We cannot exclude that the excess fluctuations are produced by more than one population. Because of poor statistics, the current relatively broad photometric bands employed here do not allow distinguishing the exact nature of these compact objects or if a fraction of the fluctuations have instead a local origin.

Astrophysical measurements away from the 1 AU orbit of Earth can enable several astrophysical science cases that are challenging or impossible to perform from Earthbound platforms, including: building a detailed understanding of the extragalactic background light throughout the electromagnetic spectrum; measurements of the properties of dust and ice in the inner and outer solar system; determinations of the mass of planets and stellar remnants far from luminous stars using gravitational microlensing; and stable time-domain astronomy. Though potentially transformative for astrophysics, opportunities to fly instrumentation capable of these measurements are rare, and a mission to the distant solar system that includes instrumentation expressly designed to perform astrophysical science, or even one primarily for a different purpose but capable of precise astronomical investigation, has not yet been flown. In this White Paper, we describe the science motivations for this kind of measurement, and advocate for future flight opportunities that permit intersectional collaboration and cooperation to make these science investigations a reality.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes).

We describe the results of principal component analysis (PCA) of up-the-ramp sampled IR array data from the HST WFC3 IR, JWST NIRSpec, and prototype WFIRST WFI detectors. These systems use respectively Teledyne H1R, H2RG, and H4RG-10 near-IR detector arrays with a variety of IR array controllers. The PCA shows that the Legendre polynomials approximate the principal components of these systems (i.e. they roughly diagonalize the covariance matrix). In contrast to the monomial basis that is widely used for polynomial fitting and linearization today, the Legendre polynomials are an orthonormal basis. They provide a quantifiable, compact, and (nearly) linearly uncorrelated representation of the information content of the data. By fitting a few Legendre polynomials, nearly all of the meaningful information in representative WFC3 astronomical datacubes can be condensed from 15 up-the-ramp samples down to 6 compressible Legendre coefficients per pixel. The higher order coefficients contain time domain information that is lost when one projects up-the-ramp sampled datacubes onto 2-dimensional images by fitting a straight line, even if the data are linearized before fitting the line. Going forward, we believe that this time domain information is potentially important for disentangling the various non-linearities that can affect IR array observations, i.e. inherent pixel non-linearity, persistence, burn in, brighter-fatter effect, (potentially) non-linear inter-pixel capacitance (IPC), and perhaps others.

The source-subtracted cosmic infrared background (CIB) fluctuations uncovered in deep Spitzer data cannot be explained by known galaxy populations and appear strongly coherent with unresolved cosmic X-ray background (CXB). This suggests that the source-subtracted CIB contains emissions from significantly abundant accreting black holes (BHs). We show that theoretically such populations would have the angular power spectrum which is largely independent of the epochs occupied by these sources, provided they are at z>~ 4, offering an important test of the origin of the new populations. Using the current measurements we reconstruct the underlying soft X-ray CXB from the new sources and show that its fluctuations, while consistent with a high-z origin, have an amplitude that cannot be reached in direct measurements with the foreseeable X-ray space missions. This necessitates application of the methods developed by the authors to future IR and X-ray datasets, which must cover large areas of the sky in order to measure the signal with high precision. The LIBRAE project within ESA's Euclid mission will probe source-subtracted CIB over ~1/2 the sky at three near-IR bands, and its cross-power with unresolved CXB can be measured then from the concurrent eROSITA mission covering the same areas of the sky. We discuss the required methodology for this measurement and evaluate its projected S/N to show the unique potential of this experimental configuration to accurately probe the CXB from the new BH sources and help identify their epochs.

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 HIgh-Resolution Mid-infrarEd Spectrometer (HIRMES) is the 3rd Generation Instrument for the Stratospheric Observatory For Infrared Astronomy (SOFIA), currently in development at the NASA Goddard Space Flight Center (GSFC), and due for commissioning in 2019. By combining direct-detection Transition Edge Sensor (TES) bolometer arrays, grating-dispersive spectroscopy, and a host of Fabry-Perot tunable filters, HIRMES will provide the ability for High Resolution (R~100,000), Mid-Resolution (R~10,000), and Low-Resolution (R~600) slit-spectroscopy, and 2D Spectral Imaging (R~2000 at selected wavelengths) over the 25 - 122 {\mu}m mid-far infrared waveband. The driving science application is the evolution of proto-planetary systems via measurements of water-vapor, water-ice, deuterated hydrogen (HD), and neutral oxygen lines. However, HIRMES has been designed to be as flexible as possible to cover a wide range of science cases that fall within its phase-space, all whilst reaching sensitivities and observing powers not yet seen thus far on SOFIA, providing unique observing capabilities which will remain unmatched for decades.

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 outer solar system provides a unique, quiet vantage point from which to observe the universe around us, where measurements could enable several niche astrophysical science cases that are too difficult to perform near Earth. NASA's New Horizons mission comprises an instrument package that provides imaging capability from ultraviolet (UV) to near-infrared (near-IR) wavelengths with moderate spectral resolution located beyond the orbit of Pluto. A carefully designed survey with New Horizons can optimize the use of expendable propellant and the limited data telemetry bandwidth to allow several measurements, including a detailed understanding of the cosmic extragalactic background light; studies of the local and extragalactic UV background; measurements of the properties of dust and ice in the outer solar system; confirmation and characterization of transiting exoplanets; determinations of the mass of dark objects using gravitational microlensing; and rapid follow-up of transient events. New Horizons is currently in an extended mission designed to focus on Kuiper Belt science that will conclude in 2021. The astrophysics community has a unique, generational opportunity to use this mission for astronomical observation at heliocentric distances beyond 50 au in the next decade. In this paper, we discuss the potential science cases for such an extended mission, and provide an initial assessment of the most important operational requirements and observation strategies it would require. We conclude that New Horizons is capable of transformative science, and that it would make a valuable and unique asset for astrophysical science that is unlikely to be replicated in the near future.

We study the source-subtracted near-infrared and X-ray background fluctuations of the COSMOS field using data from the Spitzer SPLASH program ($\sim$1272 hours) and Chandra COSMOS Legacy Survey (4.6 Ms). The new auto power spectra of the cosmic infrared and X-ray background fluctuations reach maximum angular scales of $\sim$ 3000$''$ and $\sim$ 5000$''$, respectively. We measure the cross power spectra between each infrared and X-ray band and calculate the mean power above 20$''$. We find that the soft X-ray band is correlated with 3.6 and 4.5$\mu$m at $\sim$ 4 $\sigma$ significance level. The significance between hard X-ray and the 3.6$\mu$m (4.5$\mu$m) band is $\sim$2.2 $\sigma$ ($\sim$3.8$\sigma$). The combined infrared (3.6 + 4.5$\mu$m) data are correlated with the X-ray data in soft ([0.5-2] keV), hard ([2-7] keV) and broad ([0.5-7] keV) bands at $\sim$5.6$\sigma$, $\sim$4.4$\sigma$ and $\sim$6.6$\sigma$ level, respectively. We compare the new measurements with existing models for the contributions from known populations at $z$$<$7, which are not subtracted. The model predictions are consistent with the measurements but we cannot rule out contributions from other components, such as Direct Collapse Black Holes (DCBH). However, the stacked cross-power spectra, combining other available data, show excess fluctuations about an order of magnitude on average at $\sim$4$\sigma$ confidence at scales within $\sim$300$''$. By studying the X-ray SED of the cross-power signal, assuming no significant variation from the infrared, we find that its shape is consistent with DCBHs.

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.

The cosmic infrared background (CIB) contains emissions accumulated over the entire history of the Universe, including from objects inaccessible to individual telescopic studies. The near-IR (~1-10 mic) part of the CIB, and its fluctuations, reflects emissions from nucleosynthetic sources and gravitationally accreting black holes (BHs). If known galaxies are removed to sufficient depths the source-subtracted CIB fluctuations at near-IR can reveal sources present in the first-stars-era and possibly new stellar populations at more recent times. This review discusses the recent progress in this newly emerging field which identified, with new data and methodology, significant source-subtracted CIB fluctuations substantially in excess of what can be produced by remaining known galaxies. The CIB fluctuations further appear coherent with unresolved cosmic X-ray background (CXB) indicating a very high fraction of BHs among the new sources producing the CIB fluctuations. These observations have led to intensive theoretical efforts to explain the measurements and their properties. While current experimental configurations have limitations in decisively probing these theories, their potentially remarkable implications will be tested in the upcoming CIB measurements with the ESA's Euclid dark energy mission. We describe the goals and methodologies of LIBRAE (Looking at Infrared Background Radiation with Euclid), a NASA-selected project for CIB science with Euclid, which has the potential for transforming the field into a new area of precision cosmology.

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.

We present new measurements of the large scale clustering component of the cross-power spectra of the source-subtracted Spitzer-IRAC Cosmic Infrared Background (CIB) and Chandra-ACIS Cosmic X-ray Background (CXB) surface brightness fluctuations. Our investigation uses data from the Chandra Deep Field South (CDFS), Hubble Deep Field North (HDFN), EGS/AEGIS field and UDS/SXDF surveys, comprising 1160 Spitzer hours and $\sim$ 12 Ms of Chandra data collected over a total area of 0.3 deg$^2$. We report the first ($>$5$\sigma$) detection of a cross-power signal on large angular scales >20 arcsec between [0.5-2] keV and the 3.6 and 4.5$\mu$m bands, at $\sim$5$\sigma$ and 6.3$\sigma$ significance, respectively. The correlation with harder X-ray bands is marginally significant. Comparing the new observations with existing models for the contribution of the known unmasked source population at $z<$7, we find an excess of about an order of magnitude at 5$\sigma$ confidence. We discuss possible interpretations for the origin of this excess in terms of the contribution from accreting early black holes, including both direct collapse black holes and primordial black holes, as well as from scattering in the interstellar medium and intra-halo light.

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

Near-infrared array detectors, like the \JWST NIRSpec's Teledyne's H2RGs, often provide reference pixels and a reference output. These are used to remove correlated noise. Improved Reference Sampling and Subtraction (\IRSSquare, pronounced "IRS-square") is a statistical technique for using this reference information optimally in a least squares sense. Compared to "traditional" H2RG readout, \IRSSquare uses a different clocking pattern to interleave many more reference pixels into the data than is otherwise possible. Compared to standard reference correction techniques, \IRSSquare subtracts the reference pixels and reference output using a statistically optimized set of frequency dependent weights. The benefits include somewhat lower noise variance and much less obvious correlated noise. NIRSpec's \IRSSquare images are cosmetically clean, with less $1/f$ banding than in traditional data from the same system. This article describes the \IRSSquare clocking pattern and presents the equations that are needed to use \IRSSquare in systems other than NIRSpec. For NIRSpec, applying these equations is already an option in the calibration pipeline. As an aid to instrument builders, we provide our prototype \IRSSquare calibration software and sample \JWST NIRSpec data. The same techniques are applicable to other detector systems, including those based on Teledyne's H4RG arrays. The H4RG's "interleaved reference pixel readout" mode is effectively one \IRSSquare pattern.

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.

Stellar streams have become central to studies of the interaction histories of nearby galaxies. To characterize the most prominent parts of the stellar stream around the well-known nearby (d = 17 Mpc) edge-on disk galaxy NGC 5907, we have obtained and analyzed new, deep gri Subaru/Suprime-Cam and 3.6 micron Spitzer/Infrared Array Camera (IRAC) observations. Combining the near-infrared 3.6 micron data with visible-light images allows us to use a long wavelength baseline to estimate the metallicity and age of the stellar population along a ~60 kpc long segment of the stream. We have fitted the stellar spectral energy distribution (SED) with a single-burst stellar population synthesis model and we use it to distinguish between the proposed satellite accretion and minor/major merger formation models of the stellar stream around this galaxy. We conclude that a massive minor merger (stellar mass ratio of at least 1:8) can best account for the metallicity of -0.3 inferred along the brightest parts of the stream.

We have performed a specific observational test to measure the effect that the zodiacal light can have on measurements of the spatial fluctuations of the near-IR background. Previous estimates of possible fluctuations caused by zodiacal light have often been extrapolated from observations of the thermal emission at longer wavelengths and low angular resolution, or from IRAC observations of high latitude fields where zodiacal light is faint and not strongly varying with time. The new observations analyzed here target the COSMOS field, at low ecliptic latitude where the zodiacal light intensity varies by factors of $\sim2$ over the range of solar elongations at which the field can be observed. We find that the white noise component of the spatial power spectrum of the background is correlated with the modeled zodiacal light intensity. Roughly half of the measured white noise is correlated with the zodiacal light, but a more detailed interpretation of the white noise is hampered by systematic uncertainties that are evident in the zodiacal light model. At large angular scales ($\gtrsim100"$) where excess power above the white noise is observed, we find no correlation of the power with the modeled intensity of the zodiacal light. This test clearly indicates that the large scale power in the infrared background is not being caused by the zodiacal light.

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.

We present high-resolution multiwavelength radio continuum images of the region within 150$"$ of Sgr A*, revealing a number of new extended features and stellar sources in this region. First, we detect a continuous 2" east-west ridge of radio emission, linking Sgr A* and a cluster of stars associated with IRS 13N and IRS 13E. The ridge suggests that an outflow of east-west blob-like structures is emerging from Sgr A*. We also find arc-like features within the ridge with morphologies suggestive of photoevaporative protoplanetary disks. We use near-IR fluxes to show that the emission has similar characteristics to those of a protoplanetary disk irradiated by the intense radiation field at the Galactic center. This suggests that star formation has taken place within the S cluster 2$"$ from Sgr A*. We suggest that the diffuse X-ray emission associated with Sgr A* is due to an expanding hot wind produced by the mass loss from B-type main sequence stars, and/or the disks of photoevaporation of low mass YSOs at a rate ~10^{-6} solar mass per year. The proposed model naturally reduces the inferred accretion rate and is an alternative to the inflow-outflow style models to explain the underluminous nature of Sgr A*. Second, we detect new cometary radio and nea-IR sources and a striking tower of radio emission suggesting that they are tracing interaction sites of a mildly relativistic jet from Sgr A* with the atmosphere of stars and the nonthermal Sgr A East shell at a PA$\sim50-60^\circ$ with ~10^{-7} solar mass per year, and opening angle 10 degrees.

The Euclid space mission, designed to probe evolution of the Dark Energy, will map a large area of the sky at three adjacent near-IR filters, Y, J and H. This coverage will also enable mapping source-subtracted cosmic infrared background (CIB) fluctuations with unprecedented accuracy on sub-degree angular scales. Here we propose methodology, using the Lyman-break tomography applied to the Euclid-based CIB maps, to accurately isolate the history of CIB emissions as a function of redshift from 10 < z < 20, and to identify the baryonic acoustic oscillations (BAOs) at those epochs. To identify the BAO signature, we would assemble individual CIB maps over conservatively large contiguous areas of >~ 400 sq deg. The method can isolate the CIB spatial spectrum by z to sub-percent statistical accuracy. We illustrate this with a specific model of CIB production at high z normalized to reproduce the measured Spitzer-based CIB fluctuation. We show that even if the latter contain only a small component from high-z sources, the amplitude of that component can be accurately isolated with the methodology proposed here and the BAO signatures at z>~ 10 are recovered well from the CIB fluctuation spatial spectrum. Probing the BAO at those redshifts will be an important test of the underlying cosmological paradigm, and would narrow the overall uncertainties on the evolution of cosmological parameters, including the Dark Energy. Similar methodology is applicable to the planned WFIRST mission, where we show that a possible fourth near-IR channel at > 2 micron would be beneficial.

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.

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.

We present radio and infrared observations indicating on-going star formation activity inside the $\sim2-5$ pc circumnuclear ring at the Galactic center. Collectively these measurements suggest a continued disk-based mode of on-going star formation has taken place near Sgr A* over the last few million years. First, VLA observations with spatial resolution 2.17$"\times0.81"$ reveal 13 water masers, several of which have multiple velocity components. The presence of interstellar water masers suggests gas densities that are sufficient for self-gravity to overcome the tidal shear of the 4$\times10^6$ \msol\, black hole. Second, SED modeling of stellar sources indicate massive YSO candidates interior to the molecular ring, supporting in-situ star formation near Sgr A* and appear to show a distribution similar to that of the counter-rotating disks of $\sim$100 OB stars orbiting Sgr A*. Some YSO candidates (e.g., IRS~5) have bow shock structures suggesting that they have have gaseous disks that are phototoevaporated and photoionized by the strong radiation field. Third, we detect clumps of SiO (2-1) and (5-4) line emission in the ring based on CARMA and SMA observations. The FWHM and luminosity of the SiO emission is consistent with shocked protostellar outflows. Fourth, two linear ionized features with an extent of $\sim0.8$ pc show blue and redshifted velocities between $+50$ and $-40$ \kms, suggesting protostellar jet driven outflows with mass loss rates of $\sim5\times10^{-5}$ solar mass yr$^{-1}$. Finally, we present the imprint of radio dark clouds at 44 GHz, representing a reservoir of molecular gas that feeds star formation activity close to Sgr A*.

We present new methodology to use cosmic infrared background (CIB) fluctuations to probe sources at 10<z<30 from a JWST/NIRCam configuration that will isolate known galaxies to 28 AB mag at 0.5--5 micron. At present significant mutually consistent source-subtracted CIB fluctuations have been identified in the Spitzer and Akari data at 2--5 micron, but we demonstrate internal inconsistencies at shorter wavelengths in the recent CIBER data. We evaluate CIB contributions from remaining galaxies and show that the bulk of the high-z sources will be in the confusion noise of the NIRCam beam, requiring CIB studies. The accurate measurement of the angular spectrum of the fluctuations and probing the dependence of its clustering component on the remaining shot noise power would discriminate between the various currently proposed models for their origin and probe the flux distribution of its sources. We show that the contribution to CIB fluctuations from remaining galaxies is large at visible wavelengths for the current instruments precluding probing the putative Lyman-break of the CIB fluctuations. We demonstrate that with the proposed JWST configuration such measurements will enable probing the Lyman break. We develop a Lyman-break tomography method to use the NIRCam wavelength coverage to identify or constrain, via the adjacent two-band subtraction, the history of emissions over 10<z<30 as the Universe comes out of the 'Dark Ages'. We apply the proposed tomography to the current Spitzer/IRAC measurements at 3.6 and 4.5 micron, to find that it already leads to interestingly low upper limit on emissions at z>30.

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.

Re-examination of the COBE DIRBE data reveals the thermal emission of several comet dust trails. The dust trails of 1P/Halley, 169P/NEAT, and 3200 Phaethon have not been previously reported. The known trails of 2P/Encke, and 73P/Schwassmann-Wachmann 3 are also seen. The dust trails have 12 and 25 micron surface brightnesses of <0.1 and <0.15 MJy/sr, respectively, which is <1% of the zodiacal light intensity. The trails are very difficult to see in any single daily image of the sky, but are evident as rapidly moving linear features in movies of the DIRBE data. Some trails are clearest when crossing through the orbital plane of the parent comet, but others are best seen at high ecliptic latitudes as the Earth passes over or under the dust trail. All these comets have known associations with meteor showers. This re-examination also reveals one additional comet and 13 additional asteroids that had not previously been recognized in the DIRBE data.

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

In order to understand the nature of the sources producing the recently uncovered cosmic infrared background (CIB) fluctuations, we study cross-correlations between the fluctuations in the source-subtracted CIB from Spitzer/IRAC data and the unresolved cosmic X-ray background from deep Chandra observations. Our study uses data from the EGS/AEGIS field, where both data sets cover an 8'x45' region of the sky. Our measurement is the cross-power spectrum between the IR and X-ray data. The cross-power signal between the IRAC maps at 3.6 um and 4.5 um and the Chandra [0.5-2] keV data has been detected, at angular scales 20'', with an overall significance of 3.8 sigma and 5.6 sigma, respectively. At the same time we find no evidence of significant cross-correlations at the harder Chandra bands. The cross-correlation signal is produced by individual IR sources with 3.6 um and 4.5 um magnitudes m_AB 25-26 and [0.5-2] keV X-ray fluxes 7e-17 cgs.We determine that at least 15%-25% of the large scale power of the CIB fluctuations is correlated with the spatial power spectrum of the X-ray fluctuations. If this correlation is attributed to emission from accretion processes at both IR and X-ray wavelengths, this implies a much higher fraction of accreting black holes than among the known populations. We discuss the various possible origins for the cross-power signal and show that neither local foregrounds nor the known remaining normal galaxies and active galactic nuclei can reproduce the measurements. These observational results are an important new constraint on theoretical modeling of the near-IR CIB fluctuations.

ALMA observations of the Galactic center with spatial resolution $2.61"\times0.97"$ resulted in the detection of 11 SiO (5-4) clumps of molecular gas within 0.6pc (15$"$) of Sgr A*, interior to the 2-pc circumnuclear molecular ring. The three SiO (5-4) clumps closest to Sgr A* show the largest central velocities, $\sim150$ \kms, and broadest asymmetric linewidths with full width zero intensity (FWZI) $\sim110-147$ \kms. The remaining clumps, distributed mainly to the NE of the ionized mini-spiral, have narrow FWZI ($\sim18-56$ \kms). Using CARMA SiO (2-1) data, LVG modeling of the the SiO line ratios for the broad velocity clumps, constrains the column density N(SiO) $\sim10^{14}$ cm$^{-2}$, and the H$_2$ gas density n$_{\rm H_2}=(3-9)\times10^5$ cm$^{-3}$ for an assumed kinetic temperature 100-200K. The SiO clumps are interpreted as highly embedded protostellar outflows, signifying an early stage of massive star formation near Sgr A* in the last $10^4-10^5$ years. Support for this interpretation is provided by the SiO (5-4) line luminosities and velocity widths which lie in the range measured for protostellar outflows in star forming regions in the Galaxy. Furthermore, SED modeling of stellar sources shows two YSO candidates near SiO clumps, supporting in-situ star formation near Sgr A*. We discuss the nature of star formation where the gravitational potential of the black hole dominates. In particular, we suggest that external radiative pressure exerted on self-shielded molecular clouds enhances the gas density, before the gas cloud become gravitationally unstable near Sgr A*. Alternatively, collisions between clumps in the ring may trigger gravitational collapse.

The compact radio source Sgr A* is coincident with a 4 million solar mass black hole at the dynamical center of the Galaxy and is surrounded by dense orbiting ionized and molecular gas. We present high resolution radio continuum images of the central 3' and report a faint continuous linear structure centered on Sgr A* with a PA~60 degrees. The extension of this feature appears to be terminated symmetrically by two linearly polarized structures at 8.4 GHz, ~75" from Sgr A*. A number of weak blobs of radio emission with X-ray counterparts are detected along the axis of the linear structure. The linear structure is best characterized by a mildly relativistic jet from Sgr A* with an outflow rate 10^-6 solar mass per year. The near and far-sides of the jet are interacting with orbiting ionized and molecular gas over the last 1-3 hundred years and are responsible for a 2" hole, the "minicavity", characterized by disturbed kinematics, enhanced FeII/III line emission, and diffuse X-ray gas. The estimated kinetic luminosity of the outflow is ~1.2x10^{41} erg/s, so the interaction with the bar may be responsible for the Galactic center X-ray flash inferred to be responsible for much of the fluorescent Fe Kalpha line emission from the inner 100pc of the Galaxy.

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.

We extend previous measurements of cosmic infrared background (CIB) fluctuations to ~ 1 deg using new data from the Spitzer Extended Deep Survey. Two fields, with depths of ~12 hr/pixel over 3 epochs, are analyzed at 3.6 and 4.5 mic. Maps of the fields were assembled using a self-calibration method uniquely suitable for probing faint diffuse backgrounds. Resolved sources were removed from the maps to a magnitude limit of AB mag ~ 25, as indicated by the level of the remaining shot noise. The maps were then Fourier-transformed and their power spectra were evaluated. Instrumental noise was estimated from the time-differenced data, and subtracting this isolates the spatial fluctuations of the actual sky. The power spectra of the source-subtracted fields remain identical (within the observational uncertainties) for the three epochs indicating that zodiacal light contributes negligibly to the fluctuations. Comparing to 8 mic power spectra shows that Galactic cirrus cannot account for the fluctuations. The signal appears isotropically distributed on the sky as required for an extragalactic origin. The CIB fluctuations continue to diverge to > 10 times those of known galaxy populations on angular scales out to < 1 deg. The low shot noise levels remaining in the diffuse maps indicate that the large scale fluctuations arise from the spatial clustering of faint sources well below the confusion noise. The spatial spectrum of these fluctuations is in reasonable agreement with an origin in populations clustered according to the standard cosmological model (LCDM) at epochs coinciding with the first stars era.

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.

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.

We present results from our spectroscopic study, using the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, designed to identify massive young stellar objects (YSOs) in the Galactic Center (GC). Our sample of 107 YSO candidates was selected based on IRAC colors from the high spatial resolution, high sensitivity Spitzer/IRAC images in the Central Molecular Zone (CMZ), which spans the central ~300 pc region of the Milky Way Galaxy. We obtained IRS spectra over 5um to 35um using both high- and low-resolution IRS modules. We spectroscopically identify massive YSOs by the presence of a 15.4um shoulder on the absorption profile of 15um CO2 ice, suggestive of CO2 ice mixed with CH3OH ice on grains. This 15.4um shoulder is clearly observed in 16 sources and possibly observed in an additional 19 sources. We show that 9 massive YSOs also reveal molecular gas-phase absorption from CO2, C2H2, and/or HCN, which traces warm and dense gas in YSOs. Our results provide the first spectroscopic census of the massive YSO population in the GC. We fit YSO models to the observed spectral energy distributions and find YSO masses of 8 - 23 Msun, which generally agree with the masses derived from observed radio continuum emission. We find that about 50% of photometrically identified YSOs are confirmed with our spectroscopic study. This implies a preliminary star formation rate of ~0.07 Msun/yr at the GC.

We present results of 2 mm observations of the Crab Nebula, obtained using the Goddard-IRAM Superconducting 2 Millimeter Observer (GISMO) bolometer camera on the IRAM 30 m telescope. Additional 3.3 mm observations with the MUSTANG bolometer array on the Green Bank Telescope are also presented. The integrated 2 mm flux density of the Crab Nebula provides no evidence for the emergence of a second synchrotron component that has been proposed. It is consistent with the radio power law spectrum, extrapolated up to a break frequency of log(nu_{b} [GHz]) = 2.84 +/- 0.29 or nu_{b} = 695^{+651}_{-336} GHz. The Crab Nebula is well-resolved by the ~16.7" beam (FWHM) of GISMO. Comparison to radio data at comparable spatial resolution enables us to confirm significant spatial variation of the spectral index between 21 cm and 2 mm. The main effect is a spectral flattening in the inner region of the Crab Nebula, correlated with the toroidal structure at the center of the nebula that is prominent in the near-IR through X-ray regime.

A group of four compact HII regions associated with the well-known 50 km/s molecular cloud is the closest site of on-going star formation to the dynamical center of the Galaxy, at a projected distance of ~6 pc. We present a study of ionized gas based on the [NeII] (12.8 micron) line, as well as multi-frequency radio continuum, HST Pa alpha and Spitzer IRAC observations of the most compact member of the HII group, Sgr A East HII D. The radio continuum image at 6cm shows that this source breaks up into two equally bright ionized features, D1 and D2. The SED of the D source is consistent with it being due to a 25\pm3 solar mass, star with a luminosity of 8\pm3x10^4 solar luminosity. The inferred mass, effective temperature of the UV source and the ionization rate are compatible with a young O9-B0 star. The ionized features D1 and D2 are considered to be ionized by UV radiation collimated by an accretion disk. We consider that the central massive star photoevaporates its circumstellar disk on a timescale of 3x10^4 years giving a mass flux ~3x10^{-5} solar mass per yr and producing the ionized material in D1 and D2 expanding in an inhomogeneous medium. The ionized gas kinematics, as traced by the [Ne II] emission, is difficult to interpret, but it could be explained by the interaction of a bipolar jet with surrounding gas along with what appears to to be a conical wall of lower velocity gas. The other HII regions, Sgr A East A-C, have morphologies and kinematics that more closely resemble cometary flows seen in other compact HII regions, where gas moves along a paraboloidal surface formed by the interaction of a stellar wind with a molecular cloud.

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.

This paper provides a detailed description of the data reduction and analysis procedures that have been employed in our previous studies of spatial fluctuation of the cosmic infrared background (CIB) using deep Spitzer IRAC observations. The self-calibration we apply removes a strong instrumental signal from the fluctuations which would otherwise corrupt our results. The procedures and results for masking bright sources, and modeling faint sources down to levels set by the instrumental noise are presented. Various tests are performed to demonstrate that the resulting power spectra of these fields are not dominated by instrumental or procedural effects. These tests indicate that the large scale (>~30') fluctuations that remain in the deepest fields are not directly related to the galaxies that are bright enough to be individually detected. We provide the parameterization of these power spectra in terms of separate instrument noise, shot noise, and power law components. Our measurements of spatial fluctuations of the CIB intensity indicate the mean emission from the objects producing the fluctuations is quite low (>~1 nW m-2 sr-1 at 3-5 micron), and thus consistent with current gamma-ray absorption constraints. The source of the fluctuations may be high-z Population III objects, or a more local component of very low luminosity objects with clustering properties that differ from the resolved galaxies. Finally, we discuss the prospects of the upcoming space-based surveys to directly measure the epochs inhabited by the populations producing these source-subtracted CIB fluctuations, and to isolate the individual fluxes of these populations.

We report the detection of several molecular gas-phase and ice absorption features in three photometrically-selected young stellar object (YSO) candidates in the central 280 pc of the Milky Way. Our spectra, obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, reveal gas-phase absorption from CO2 (15.0um), C2H2 (13.7um) and HCN (14.0um). We attribute this absorption to warm, dense gas in massive YSOs. We also detect strong and broad 15um CO2 ice absorption features, with a remarkable double-peaked structure. The prominent long-wavelength peak is due to CH3OH-rich ice grains, and is similar to those found in other known massive YSOs. Our IRS observations demonstrate the youth of these objects, and provide the first spectroscopic identification of massive YSOs in the Galactic Center.

The central kpc of the Milky Way might be expected to differ significantly from the rest of the Galaxy with regard to gas dynamics and the formation of YSOs. We probe this possibility with mid-infrared observations obtained with IRAC and MIPS on Spitzer and with MSX. We use color-color diagrams and SED fits to explore the nature of YSO candidates (including objects with 4.5 micron excesses possibly due to molecular emission). There is an asymmetry in the distribution of the candidate YSOs, which tend to be found at negative Galactic longitudes; this behavior contrasts with that of the molecular gas, approximately 2/3 of which is at positive longitudes. The small scale height of these objects suggests that they are within the Galactic center region and are dynamically young. They lie between two layers of infrared dark clouds and may have originated from these clouds. We identify new sites for this recent star formation. The methanol masers appear to be associated with young, embedded YSOs characterized by 4.5 micron excesses. We use the SEDs of these sources to estimate their physical characteristics. Within the central 400x50 pc (|l|<1.3\degr and |b|<10') the star formation rate based on the identification of Stage I evolutionary phase of YSO candidates is about 0.14 solar mass/yr. We suggest that a recent burst of star formation took place within the last 10^5 years. This suggestion is also consistent with estimates of star formation rates within the last ~10^7 years showing a peak around 10^5 years ago. Lastly, we find that the Schmidt-Kennicutt Law applies well in the central 400 pc of the Galaxy. This implies that star formation does not appear to be dramatically affected by the extreme physical conditions in the Galactic center region.

We have studied the correlation between 2357 Chandra X-ray point sources in a 40 x 40 parsec field and ~20,000 infrared sources we observed in the corresponding subset of our 2 x 1.4 degree Spitzer/IRAC Galactic Center Survey at 3.6-8.0 um, using various spatial and X-ray hardness thresholds. The correlation was determined for source separations of less than 0.5", 1" or 2". Only the soft X-ray sources show any correlation with infrared point sources on these scales, and that correlation is very weak. The upper limit on hard X-ray sources that have infrared counterparts is <1.7% (3 sigma). However, because of the confusion limit of the IR catalog, we only detect IR sources with absolute magnitudes < ~1. As a result, a stronger correlation with fainter sources cannot be ruled out. Only one compact infrared source, IRS 13, coincides with any of the dozen prominent X-ray emission features in the 3 x 3 parsec region centered on Sgr A*, and the diffuse X-ray and infrared emission around Sgr A* seems to be anti-correlated on a few-arcsecond scale. We compare our results with previous identifications of near-infrared companions to Chandra X-ray sources.

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

A mid-infrared (3.6-8 um) survey of the Galactic Center has been carried out with the IRAC instrument on the Spitzer Space Telescope. This survey covers the central 2x1.4 degree (~280x200 pc) of the Galaxy. At 3.6 and 4.5 um the emission is dominated by stellar sources, the fainter ones merging into an unresolved background. At 5.8 and 8 um the stellar sources are fainter, and large-scale diffuse emission from the ISM of the Galaxy's central molecular zone becomes prominent. The survey reveals that the 8 to 5.8 um color of the ISM emission is highly uniform across the surveyed region. This uniform color is consistent with a flat extinction law and emission from polycyclic aromatic hydrocarbons (PAHs). Models indicate that this broadband color should not be expected to change if the incident radiation field heating the dust and PAHs is <10^4 times that of the solar neighborhood. The few regions with unusually red emission are areas where the PAHs are underabundant and the radiation field is locally strong enough to heat large dust grains to produce significant 8 um emission. These red regions include compact H II regions, Sgr B1, and wider regions around the Arches and Quintuplet Clusters. In these regions the radiation field is >10^4 times that of the solar neighborhood. Other regions of very red emission indicate cases where thick dust clouds obscure deeply embedded objects or very early stages of star formation.

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.

We have obtained Spitzer/IRAC observations of the central 2.0 x 1.4 degrees (~280 x 200 pc) of the Galaxy at 3.6-8.0 microns. A point source catalog of 1,065,565 objects is presented. The catalog includes magnitudes for the point sources at 3.6, 4.5, 5.8, and 8.0 microns, as well as JHK photometry from 2MASS. The point source catalog is confusion limited with average limits of 12.4, 12.1, 11.7, and 11.2 magnitudes for [3.6], [4.5], [5.8], and [8.0], respectively. We find that the confusion limits are spatially variable because of stellar surface density, background surface brightness level, and extinction variations across the survey region. The overall distribution of point source density with Galactic latitude and longitude is essentially constant, but structure does appear when sources of different magnitude ranges are selected. Bright stars show a steep decreasing gradient with Galactic latitude, and a slow decreasing gradient with Galactic longitude, with a peak at the position of the Galactic center. From IRAC color-magnitude and color-color diagrams, we conclude that most of the point sources in our catalog have IRAC magnitudes and colors characteristic of red giant and AGB stars.

We study the possible contribution of optical galaxies detected with the {\it Hubble} ACS instrument to the near-IR cosmic infrared (CIB) fluctuations in deep {\it Spitzer} images. The {\it Spitzer} data used in this analysis are obtained in the course of the GOODS project from which we select four independent $10^\prime\times10^\prime$ regions observed at both 3.6 and 4.5 \um. ACS source catalogs for all of these areas are used to construct maps containing only their emissions in the ACS $B, V, i, z$-bands. We find that deep Spitzer data exhibit CIB fluctuations remaining after removal of foreground galaxies of a very different clustering pattern at both 3.6 and 4.5 \um than the ACS galaxies could contribute. We also find that there are very good correlations between the ACS galaxies and the {\it removed} galaxies in the Spitzer maps, but practically no correlations remain with the residual Spitzer maps used to identify the CIB fluctuations. These contributions become negligible on larger scales used to probe the CIB fluctuations arising from clustering. This means that the ACS galaxies cannot contribute to the large-scale CIB fluctuations found in the residual Spitzer data. The absence of their contributions also means that the CIB fluctuations arise at $z\gsim 7.5$ as the Lyman break of their sources must be redshifted past the longest ACS band, or the fluctuations have to originate in the more local but extremely low luminosity galaxies.

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.

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.

OH(1720 MHz) and methanol masers are now recognized to be excellent probes of the interactions of supernova remnants with molecular clouds and tracers of massive star formation, respectively. To better understand the nature of star formation activity in the central region of the Galaxy, we have used these two classes of masers combined with the IRAC and MIPS data to study prominent sites of ongoing star formation in the nuclear disk. The nuclear disk is characterized by massive GMCs with elevated gas temperatures, compared to their dust temperatures. We note an association between methanol masers and a class of mid-infrared ``green sources''. These highly embedded YSOs show enhanced 4.5micron emission due to excited molecular lines. The distribution of methanol masers and supernova remnants suggest a low efficiency of star formation (with the exception of Sgr B2), which we believe is due to an enhanced flux of cosmic ray electrons impacting molecular clouds in the nuclear disk. We also highlight the importance of cosmic rays in their ability to heat molecular clouds, and thus increase the gas temperature.

We discuss interpretation of the cosmic infrared background (CIB) anisotropies detected by us recently in the Spitzer IRAC based measurements. The fluctuations are approximately isotropic on the sky consistent with their cosmological origin. They remain after removal of fairly faint intervening sources and must arise from a population which has a strong CIB clustering component with only a small shot-noise level. We discuss the constraints the data place on the luminosities, epochs and mass-to-light ratios of the indvidual sources producing them. Assuming the concordance LambdaCDM cosmology the measurements imply that the luminous sources producing them lie at cosmic times < 1 Gyr and were individually much brighter per unit mass than the present stellar populations.

Cosmic infrared background fluctuations may contain measurable contribution from objects inaccessible to current telescopic studies, such as the first stars and other luminous objects in the first Gyr of the Universe's evolution. In an attempt to uncover this contribution we have analyzed the GOODS data obtained with the Spitzer IRAC instrument, which are deeper and cover larger scales than the Spitzer data we have previously analyzed. Here we report these new measurements of the cosmic infrared background (CIB) fluctuations remaining after removing cosmic sources to fainter levels than before. The remaining anisotropies on scales > 0.5 arcmin have a significant clustering component with a low shot-noise contribution. We show that these fluctuations cannot be accounted for by instrumental effects, nor by the Solar system and Galactic foreground emissions and must arise from extragalactic sources.

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.

The deepest space and ground-based observations find metal-enriched galaxies at cosmic times when the Universe was <1 Gyr old. These stellar populations had to be preceded by the metal-free first stars, Population III. Recent cosmic microwave background polarization measurements indicate that stars started forming early when the Universe was $\lsim200$ Myr old. Theoretically it is now thought that Population III stars were significantly more massive than the present metal-rich stellar populations. While such sources will not be individually detectable by existing or planned telescopes, they would have produced significant cosmic infrared background radiation in the near- infrared, whose fluctuations reflect the conditions in the primordial density field. Here we report a measurement of diffuse flux fluctuations after removing foreground stars and galaxies. The anisotropies exceed the instrument noise and the more local foregrounds and can be attributed to emission from massive Population III stars, providing observational evidence of an era dominated by these objects.

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)

The IRAC shallow survey covers 8.5 square degrees in the NOAO Deep Wide-Field Survey in Bootes with 3 or more 30 second exposures per position. An overview of the survey design, reduction, calibration, star-galaxy separation, and initial results is provided. The survey includes approximately 370,000, 280,000, 38,000, and 34,000 sources brighter than the 5 sigma limits of 6.4, 8.8, 51, and 50 microJy at 3.6, 4.5, 5.8, and 8 microns respectively, including some with unusual spectral energy distributions.

As part of the In-Orbit Checkout activities for the Spitzer Space Telescope, the IRAC team carried out a deep observation (average integration time ~8 hours) of a field surrounding the bright QSO HS 1700+6416. This field contains several hundred z~3 Lyman-break galaxy candidates, and we report here on their mid-infrared properties, including the IRAC detection rate, flux densities and colors, and the results of fitting population synthesis models to the optical, near-infrared, and IRAC magnitudes. The results of the model-fitting show that previous optical/near-infrared studies of LBGs were not missing large, hidden old stellar populations. The LBG candidates' properties are consistent with those of massive, star-forming galaxies at z~3. Other IRAC sources in the same field have similar properties, so IRAC selection may prove a promising method of finding additional high-redshift galaxies.

Infrared source counts at wavelengths 3 < lambda < 10 um cover more than 10 magnitudes in source brightness, four orders of magnitude in surface density, and reach an integrated surface density of 10^5 sources/deg^2. At m<14 mag, most of the sources are Galactic stars, in agreement with models. After removal of Galactic stars, galaxy counts are consistent with what few measurements exist at nearby wavelengths. At 3.6 and 4.5 um, the galaxy counts follow the expectations of a Euclidean world model down to ~16 mag and drop below the Euclidean curve for fainter magnitudes. Counts at these wavelengths begin to show decreasing completeness around magnitude 19.5. At 5.8 and 8 um, the counts relative to a Euclidean world model show a large excess at bright magnitudes. This is probably because local galaxies emit strongly in the aromatic dust (``PAH'') features. The counts at 3.6 um resolve <50% of the Cosmic Infrared Background at that wavelength.

Following the successful mapping of the last scattering surface by WMAP and balloon experiments, the epoch of the first stars, when Population III stars formed, is emerging as the next cosmological frontier. It is not clear what these stars' properties were, when they formed or how long their era lasted before leading to the stars and galaxies we see today. We show that these questions can be answered with the current and future measurements of the near-IR cosmic infrared background (CIB). Theoretical arguments suggest that Population III stars were very massive and short-lived stars that formed at $z\sim 10-20$ at rare peaks of the density field in the cold-dark-matter Universe. Because Population III stars probably formed individually in small mini-halos, they are not directly accessible to current telescopic studies. We show that these stars left a strong and measurable signature via their contribution to the CIB anisotropies for a wide range of their formation scenarios. The excess in the recently measured near-IR CIB anisotropies over that from normal galaxies can be explained by contribution from early Population III stars. These results imply that Population III were indeed very massive stars and their epoch started at $z\sim 20$ and lasted past $z\lsim 13$. We show the importance of accurately measuring the CIB anisotropies produced by Population III with future space-based missions.

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.

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.

With high sensitivity imaging arrays, accurate calibration is essential to achieve the limits of detection of space observatories. One can simultaneously extract information about the scene being observed and the calibration properties of the detector and imaging system from redundant dithered images of a scene. There are large differences in the effectiveness of dithering strategies for allowing the separation of detector properties from sky brightness measurements. In this paper, we quantify these differences by developing a figure of merit (FOM) for dithering procedures based on their usefulness for allowing calibration on all spatial scales. The figure of merit measures how well the gain characteristics of the detector are encoded in the measurements, and is independent of the techniques used to analyze the data. Patterns similar to the antenna arrangements of radio interferometers with good u-v plane coverage, are found to have good figures of merit. We present patterns for both deep surveys of limited sky areas and for shallow surveys. By choosing a strategy that encodes the calibration in the observations in an easily extractable way, we enhance our ability to calibrate our detector systems and to reach the ultimate limits of sensitivity which are required to achieve the promise of many missions.

The development of sensitive large format imaging arrays for the infrared promises to provide revolutionary capabilities for space astronomy. For example, the Infrared Array Camera (IRAC) on SIRTF will use four 256 x 256 arrays to provide background limited high spatial resolution images of the sky in the 3 to 8 micron spectral region. In order to reach the performance limits possible with this generation of sensitive detectors, calibration procedures must be developed so that uncertainties in detector calibration will always be dominated by photon statistics from the dark sky as a major system noise source. In the near infrared, where the faint extragalactic sky is observed through the scattered and reemitted zodiacal light from our solar system, calibration is particularly important. Faint sources must be detected on this brighter local foreground. We present a procedure for calibrating imaging systems and analyzing such data. In our approach, by proper choice of observing strategy, information about detector parameters is encoded in the sky measurements. Proper analysis allows us to simultaneously solve for sky brightness and detector parameters, and provides accurate formal error estimates. This approach allows us to extract the calibration from the observations themselves; little or no additional information is necessary to allow full interpretation of the data. Further, this approach allows refinement and verification of detector parameters during the mission, and thus does not depend on a priori knowledge of the system or ground calibration for interpretation of images.

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.

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.

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.

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