The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne telescope mission to search for inflationary gravitational waves from the early universe. PIPER employs two 32x40 arrays of superconducting transition-edge sensors, which operate at 100 mK. An open bucket dewar of liquid helium maintains the receiver and telescope optics at 1.7 K. We describe the thermal design of the receiver and sub-kelvin cooling with a continuous adiabatic demagnetization refrigerator (CADR). The CADR operates between 70-130 mK and provides ~10 uW cooling power at 100 mK, nearly five times the loading of the two detector assemblies. We describe electronics and software to robustly control the CADR, overall CADR performance in flight-like integrated receiver testing, and practical considerations for implementation in the balloon float environment.

Astrophotonics is the application of versatile photonic technologies to channel, manipulate, and disperse guided light from one or more telescopes to achieve scientific objectives in astronomy in an efficient and cost-effective way. The developments and demands from the telecommunication industry have driven a major boost in photonic technology and vice versa in the last 40 years. The photonic platform of guided light in fibers and waveguides has opened the doors to next-generation instrumentation for both ground- and space-based telescopes in optical and near/mid-IR bands, particularly for the upcoming extremely large telescopes (ELTs). The large telescopes are pushing the limits of adaptive optics to reach close to a near-diffraction-limited performance. The photonic devices are ideally suited for capturing this AO-corrected light and enabling new and exciting science such as characterizing exoplanet atmospheres. The purpose of this white paper is to summarize the current landscape of astrophotonic devices and their scientific impact, highlight the key issues, and outline specific technological and organizational approaches to address these issues in the coming decade and thereby enable new discoveries as we embark on the era of extremely large telescopes.

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.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.

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 Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

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.

We describe feedhorn-coupled polarization-sensitive detector arrays that utilize monocrystalline silicon as the dielectric substrate material. Monocrystalline silicon has a low-loss tangent and repeatable dielectric constant, characteristics that are critical for realizing efficient and uniform superconducting microwave circuits. An additional advantage of this material is its low specific heat. In a detector pixel, two Transition-Edge Sensor (TES) bolometers are antenna-coupled to in-band radiation via a symmetric planar orthomode transducer (OMT). Each orthogonal linear polarization is coupled to a separate superconducting microstrip transmission line circuit. On-chip filtering is employed to both reject out-of-band radiation from the upper band edge to the gap frequency of the niobium superconductor, and to flexibly define the bandwidth for each TES to meet the requirements of the application. The microwave circuit is compatible with multi-chroic operation. Metalized silicon platelets are used to define the backshort for the waveguide probes. This micro-machined structure is also used to mitigate the coupling of out-of-band radiation to the microwave circuit. At 40 GHz, the detectors have a measured efficiency of 90%. In this paper, we describe the development of the 90 GHz detector arrays that will be demonstrated using the Cosmology Large Angular Scale Surveyor (CLASS) ground-based telescope.

The Cosmology Large Angular Scale Surveyor (CLASS) is a four telescope array designed to characterize relic primordial gravitational waves from inflation and the optical depth to reionization through a measurement of the polarized cosmic microwave background (CMB) on the largest angular scales. The frequencies of the four CLASS telescopes, one at 38 GHz, two at 93 GHz, and one dichroic system at 145/217 GHz, are chosen to avoid spectral regions of high atmospheric emission and span the minimum of the polarized Galactic foregrounds: synchrotron emission at lower frequencies and dust emission at higher frequencies. Low-noise transition edge sensor detectors and a rapid front-end polarization modulator provide a unique combination of high sensitivity, stability, and control of systematics. The CLASS site, at 5200 m in the Chilean Atacama desert, allows for daily mapping of up to 70\% of the sky and enables the characterization of CMB polarization at the largest angular scales. Using this combination of a broad frequency range, large sky coverage, control over systematics, and high sensitivity, CLASS will observe the reionization and recombination peaks of the CMB E- and B-mode power spectra. CLASS will make a cosmic variance limited measurement of the optical depth to reionization and will measure or place upper limits on the tensor-to-scalar ratio, $r$, down to a level of 0.01 (95\% C.L.).

Direct spectroscopic biosignature characterization (hereafter "biosignature characterization") will be a major focus for future space observatories equipped with coronagraphs or starshades. Our aim in this article is to provide an introduction to potential detector and cooling technologies for biosignature characterization. We begin by reviewing the needs. These include nearly noiseless photon detection at flux levels as low as $<0.001~\textrm{photons}~s^{-1}~\textrm{pixel}^{-1}$ in the visible and near-IR. We then discuss potential areas for further testing and/or development to meet these needs using non-cryogenic detectors (EMCCD, HgCdTe array, HgCdTe APD array), and cryogenic single photon detectors (MKID arrays and TES microcalorimeter arrays). Non-cryogenic detectors are compatible with the passive cooling that is strongly preferred by coronagraphic missions, but would add non-negligible noise. Cryogenic detectors would require active cooling, but in return deliver nearly quantum limited performance. Based on the flight dynamics of past NASA missions, we discuss reasonable vibration expectations for a large UV-Optical-IR space telescope (LUVOIR) and preliminary cooling concepts that could potentially fit into a vibration budget without being the largest element. We believe that a cooler that meets the stringent vibration needs of a LUVOIR is also likely to meet those of a starshade-based Habitable Exoplanet Imaging Mission.

The Primordial Inflation Polarization ExploreR (PIPER) is a balloon-borne telescope designed to measure the polarization of the Cosmic Microwave Background on large angular scales. PIPER will map 85% of the sky at 200, 270, 350, and 600 GHz over a series of 8 conventional balloon flights from the northern and southern hemispheres. The first science flight will use two 32x40 arrays of backshort-under-grid transition edge sensors, multiplexed in the time domain, and maintained at 100 mK by a Continuous Adiabatic Demagnetization Refrigerator. Front-end cryogenic Variable-delay Polarization Modulators provide systematic control by rotating linear to circular polarization at 3 Hz. Twin telescopes allow PIPER to measure Stokes I, Q, U, and V simultaneously. The telescope is maintained at 1.5 K in an LHe bucket dewar. Cold optics and the lack of a warm window permit sensitivity at the sky-background limit. The ultimate science target is a limit on the tensor-to-scalar ratio of r ~ 0.007, from the reionization bump to l ~ 300. PIPER's first flight will be from the Northern hemisphere, and overlap with the CLASS survey at lower frequencies. We describe the current status of the PIPER instrument.

We have used the unique far-UV imaging capability offered by a sounding rocket borne instrument to acquire observations of C/2012 S1 (ISON) when its angular separation with respect to the sun was 26.3deg, on 2013 November 20.49. At the time of observation the comet's heliocentric distance and velocity relative to the sun were rh = 0.43 AU and rh_dot = -62.7 km s^-1. Images dominated by C I 1657 A and H I 1216 A were acquired over a 1e6 x 1e6 km^2 region. The water production rate implied by the Lyman alpha observations is constrained to be Q_H2O approximately 8e29 s^-1 while the neutral carbon production rate was Q_C approximately 4e28 s^-1. The radial profile of C I was consistent with it being a dissociation product of a parent molecule with a lifetime approximately 5e4 seconds, favoring a parent other than CO. We constrain the Q_CO production rate to 5(+1.5, -7.5)e28 s^-1 with 1sigma errors derived from photon statistics. The upper limit on the Q_CO/Q_H2O < 6%.

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.

Variable-delay Polarization Modulators (VPMs) are currently being implemented in experiments designed to measure the polarization of the cosmic microwave background on large angular scales because of their capability for providing rapid, front-end polarization modulation and control over systematic errors. Despite the advantages provided by the VPM, it is important to identify and mitigate any time-varying effects that leak into the synchronously modulated component of the signal. In this paper, the effect of emission from a $300$ K VPM on the system performance is considered and addressed. Though instrument design can greatly reduce the influence of modulated VPM emission, some residual modulated signal is expected. VPM emission is treated in the presence of rotational misalignments and temperature variation. Simulations of time-ordered data are used to evaluate the effect of these residual errors on the power spectrum. The analysis and modeling in this paper guides experimentalists on the critical aspects of observations using VPMs as front-end modulators. By implementing the characterizations and controls as described, front-end VPM modulation can be very powerful for mitigating $1/f$ noise in large angular scale polarimetric surveys. None of the systematic errors studied fundamentally limit the detection and characterization of B-modes on large scales for a tensor-to-scalar ratio of $r=0.01$. Indeed, $r<0.01$ is achievable with commensurately improved characterizations and controls.

The Cosmology Large Angular Scale Surveyor (CLASS) will measure the polarization of the Cosmic Microwave Background to search for and characterize the polarized signature of inflation. CLASS will operate from the Atacama Desert and observe $\sim$70% of the sky. A variable-delay polarization modulator (VPM) modulates the polarization at $\sim$10 Hz to suppress the 1/f noise of the atmosphere and enable the measurement of the large angular scale polarization modes. The measurement of the inflationary signal across angular scales that span both the recombination and reionization features allows a test of the predicted shape of the polarized angular power spectra in addition to a measurement of the energy scale of inflation. CLASS is an array of telescopes covering frequencies of 38, 93, 148, and 217 GHz. These frequencies straddle the foreground minimum and thus allow the extraction of foregrounds from the primordial signal. Each focal plane contains feedhorn-coupled transition-edge sensors that simultaneously detect two orthogonal linear polarizations. The use of single-crystal silicon as the dielectric for the on-chip transmission lines enables both high efficiency and uniformity in fabrication. Integrated band definition has been implemented that both controls the bandpass of the single mode transmission on the chip and prevents stray light from coupling to the detectors.

Are we alone? Answering this ageless question will be a major focus for astrophysics in coming decades. Our tools will include unprecedentedly large UV-Optical-IR space telescopes working with advanced coronagraphs and starshades. Yet, these facilities will not live up to their full potential without better detectors than we have today. To inform detector development, this paper provides an overview of visible and near-IR (VISIR; $\lambda=0.4-1.8~\mu\textrm{m}$) detector needs for the Advanced Technology Large Aperture Space Telescope (ATLAST), specifically for spectroscopic characterization of atmospheric biosignature gasses. We also provide a brief status update on some promising detector technologies for meeting these needs in the context of a passively cooled ATLAST.

The far-infrared and submillimeter portions of the electromagnetic spectrum provide a unique view of the astrophysical processes present in the early universe. Micro-Spec ($\mu$-Spec), a high-efficiency direct-detection spectrometer concept working in the 450-1000-$\mu$m wavelength range, will enable a wide range of spaceflight missions that would otherwise be challenging due to the large size of current instruments and the required spectral resolution and sensitivity. This paper focuses on the $\mu$-Spec two-dimensional multimode region, where the light of different wavelengths diffracts and converges onto a set of detectors. A two-step optimization process is used to generate geometrical configurations given specific requirements on spectrometer size, operating spectral range, and performance. The canonically employed focal-plane constraints for the power combiner were removed to probe the design space in its entirety. A new four-stigmatic-point optical design solution is identified and explored for use in far-infrared and submillimeter spectroscopy.

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

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

A method is proposed and experimentally explored for in-situ calibration of complex transmission data for superconducting microwave resonators. This cryogenic calibration method accounts for the instrumental transmission response between the vector network analyzer reference plane and the device calibration plane. Once calibrated, the observed resonator response is analyzed in detail by two approaches. The first, a phenomenological model based on physically realizable rational functions, enables the extraction of multiple resonance frequencies and widths for coupled resonators without explicit specification of the circuit network. In the second, an ABCD-matrix representation for the distributed transmission line circuit is used to model the observed response from the characteristic impedance and propagation constant. When used in conjunction with electromagnetic simulations, the kinetic inductance fraction can be determined with this method with an accuracy of 2%. Datasets for superconducting microstrip and coplanar-waveguide resonator devices were investigated and a recovery within 1% of the observed complex transmission amplitude was achieved with both analysis approaches. The experimental configuration used in microwave characterization of the devices and self-consistent constraints for the electromagnetic constitutive relations for parameter extraction are also presented.

We report on the status and development of polarization-sensitive detectors for millimeter-wave applications. The detectors are fabricated on single-crystal silicon, which functions as a low-loss dielectric substrate for the microwave circuitry as well as the supporting membrane for the Transition-Edge Sensor (TES) bolometers. The orthomode transducer (OMT) is realized as a symmetric structure and on-chip filters are employed to define the detection bandwidth. A hybridized integrated enclosure reduces the high-frequency THz mode set that can couple to the TES bolometers. An implementation of the detector architecture at Q-band achieves 90% efficiency in each polarization. The design is scalable in both frequency coverage, 30-300 GHz, and in number of detectors with uniform characteristics. Hence, the detectors are desirable for ground-based or space-borne instruments that require large arrays of efficient background-limited cryogenic detectors.

The Cosmology Large Angular Scale Surveyor (CLASS) is an experiment to measure the signature of a gravita-tional-wave background from inflation in the polarization of the cosmic microwave background (CMB). CLASS is a multi-frequency array of four telescopes operating from a high-altitude site in the Atacama Desert in Chile. CLASS will survey 70\% of the sky in four frequency bands centered at 38, 93, 148, and 217 GHz, which are chosen to straddle the Galactic-foreground minimum while avoiding strong atmospheric emission lines. This broad frequency coverage ensures that CLASS can distinguish Galactic emission from the CMB. The sky fraction of the CLASS survey will allow the full shape of the primordial B-mode power spectrum to be characterized, including the signal from reionization at low $\ell$. Its unique combination of large sky coverage, control of systematic errors, and high sensitivity will allow CLASS to measure or place upper limits on the tensor-to-scalar ratio at a level of $r=0.01$ and make a cosmic-variance-limited measurement of the optical depth to the surface of last scattering, $\tau$.

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne cosmic microwave background (CMB) polarimeter designed to search for evidence of inflation by measuring the large-angular scale CMB polarization signal. BICEP2 recently reported a detection of B-mode power corresponding to the tensor-to-scalar ratio r = 0.2 on ~2 degree scales. If the BICEP2 signal is caused by inflationary gravitational waves (IGWs), then there should be a corresponding increase in B-mode power on angular scales larger than 18 degrees. PIPER is currently the only suborbital instrument capable of fully testing and extending the BICEP2 results by measuring the B-mode power spectrum on angular scales $\theta$ = ~0.6 deg to 90 deg, covering both the reionization bump and recombination peak, with sensitivity to measure the tensor-to-scalar ratio down to r = 0.007, and four frequency bands to distinguish foregrounds. PIPER will accomplish this by mapping 85% of the sky in four frequency bands (200, 270, 350, 600 GHz) over a series of 8 conventional balloon flights from the northern and southern hemispheres. The instrument has background-limited sensitivity provided by fully cryogenic (1.5 K) optics focusing the sky signal onto four 32x40-pixel arrays of time-domain multiplexed Transition-Edge Sensor (TES) bolometers held at 140 mK. Polarization sensitivity and systematic control are provided by front-end Variable-delay Polarization Modulators (VPMs), which rapidly modulate only the polarized sky signal at 3 Hz and allow PIPER to instantaneously measure the full Stokes vector (I, Q, U, V) for each pointing. We describe the PIPER instrument and progress towards its first flight.

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.

We report on the diffusive-ballistic thermal conductance of multi-moded single-crystal silicon beams measured below 1 K. It is shown that the phonon mean-free-path $\ell$ is a strong function of the surface roughness characteristics of the beams. This effect is enhanced in diffuse beams with lengths much larger than $\ell$, even when the surface is fairly smooth, 5-10 nm rms, and the peak thermal wavelength is 0.6 $\mu$m. Resonant phonon scattering has been observed in beams with a pitted surface morphology and characteristic pit depth of 30 nm. Hence, if the surface roughness is not adequately controlled, the thermal conductance can vary significantly for diffuse beams fabricated across a wafer. In contrast, when the beam length is of order $\ell$, the conductance is dominated by ballistic transport and is effectively set by the beam area. We have demonstrated a uniformity of $\pm$8% in fractional deviation for ballistic beams, and this deviation is largely set by the thermal conductance of diffuse beams that support the micro-electro-mechanical device and electrical leads. In addition, we have found no evidence for excess specific heat in single-crystal silicon membranes. This allows for the precise control of the device heat capacity with normal metal films. We discuss the results in the context of the design and fabrication of large-format arrays of far-infrared and millimeter wavelength cryogenic detectors.

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.

Image registration, or alignment of two or more images covering the same scenes or objects, is of great interest in many disciplines such as remote sensing, medical imaging, astronomy, and computer vision. In this paper, we introduce a new application of image registration algorithms. We demonstrate how through a wavelet based image registration algorithm, engineers can evaluate stability of Micro-Electro-Mechanical Systems (MEMS). In particular, we applied image registration algorithms to assess alignment stability of the MicroShutters Subsystem (MSS) of the Near Infrared Spectrograph (NIRSpec) instrument of the James Webb Space Telescope (JWST). This work introduces a new methodology for evaluating stability of MEMS devices to engineers as well as a new application of image registration algorithms to computer scientists.

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 investigate the polarization modulation properties of a variable-delay polarization modulator (VPM). The VPM modulates polarization via a variable separation between a polarizing grid and a parallel mirror. We find that in the limit where the wavelength is much larger than the diameter of the metal wires that comprise the grid, the phase delay derived from the geometric separation between the mirror and the grid is sufficient to characterize the device. However, outside of this range, additional parameters describing the polarizing grid geometry must be included to fully characterize the modulator response. In this paper, we report test results of a VPM at wavelengths of 350 microns and 3 mm. Electromagnetic simulations of wire grid polarizers were performed and are summarized using a simple circuit model that incorporates the loss and polarization properties of the device.

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

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

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.

[abridged] We present constraints on the nature of the first galaxies selected at 350 microns. The sample includes galaxies discovered in the deepest blank-field survey at 350 microns (in the Bootes Deep Field) and also later serendipitous detections in the Lockman Hole. Spectral energy distribution templates are fit to identified counterparts, and the sample is found to comprise IR-luminous galaxies at 1<z<3 predominantly powered by star formation. The first spectrum of a 350-micron-selected galaxy provides an additional confirmation, showing prominent dust grain features typically associated with star-forming galaxies. Compared to submillimeter galaxies selected at 850 and 1100 microns, galaxies selected at 350 microns have a similar range of far-infrared color temperatures. However, no 350-micron-selected sources are reliably detected at 850 or 1100 microns. Galaxies in our sample with redshifts 1<z<2 show a tight correlation between the far- and mid-infrared flux densities, but galaxies at higher redshifts show a large dispersion in their mid- to far-infrared colors. The 350-micron flux densities (15<S(350)<40 mJy) place these objects near the Herschel/SPIRE 350-micron confusion threshold, with the lower limit on the star formation rate density suggesting the bulk of the 350-micron contribution will come from less luminous infrared sources and normal galaxies.

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.

MUSTANG is a 90 GHz bolometer camera built for use as a facility instrument on the 100 m Robert C. Byrd Green Bank radio telescope (GBT). MUSTANG has an 8 by 8 focal plane array of transition edge sensor bolometers read out using time-domain multiplexed SQUID electronics. As a continuum instrument on a large single dish MUSTANG has a combination of high resolution (8") and good sensitivity to extended emission which make it very competitive for a wide range of galactic and extragalactic science. Commissioning finished in January 2008 and some of the first science data have been collected.

The first structures in the Universe formed at z>7, at higher redshift than all currently known galaxies. Since GRBs are brighter than other cosmological sources at high redshift and exhibit simple power-law afterglow spectra that is ideal for absorption studies, they serve as powerful tools for studying the early universe. New facilities planned for the coming decade will be able to obtain a large sample of high-redshift GRBs. Such a sample would constrain the nature of the first stars, galaxies, and the reionization history of the Universe.

We explore the design of a space mission, Project Lyman, which has the goal of quantifying the ionization history of the universe from the present epoch to a redshift of z ~ 3. Observations from WMAP and SDSS show that before a redshift of z >~ 6 the first collapsed objects, possibly dwarf galaxies, emitted Lyman continuum (LyC) radiation shortward of 912 A, reionizing most of the universe. How LyC escapes from galactic environments, whether it induces positive or negative feedback on the local and global collapse of structures, and the role played by clumping, molecules, metallicity and dust are major unanswered theoretical questions, requiring observational constraint. Numerous intervening Lyman limit systems, which frustrate the detection of LyC from high z objects, thin below z ~ 3 where there are a few objects with apparently very high fesc. At low z there are only controversial detections and a handful of upper limits. A wide-field multi-object spectroscopic survey with moderate spectral and spatial resolution can quantify fesc within diverse spatially resolved galactic environments over redshifts with significant evolution in galaxy assemblage and quasar activity. It can also calibrate LyC escape against Ly-alpha escape, providing an essential tool to JWST for probing the beginnings of reionization. We present calculations showing the evolution of the characteristic apparent magnitude of star-forming galaxy luminosity functions at 900 A, as a function of redshift and assumed escape fraction to determine the required aperture for detecting LyC. We review our efforts to build a pathfinding dual order multi-object spectro/telescope with a (0.5deg)^2 field-of-view, using a GSFC microshutter array, and crossed delay-line micro-channel plate detector.

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.

We present first results of Hertz/VPM, the first submillimeter polarimeter employing the dual Variable-delay Polarization Modulator (dual-VPM). This device differs from previously used polarization modulators in that it operates in translation rather than mechanical rotation. We discuss the basic theory behind this device, and its potential advantages over the commonly used half wave plate (HWP). The dual-VPM was tested both at the Submillimeter Telescope Observatory (SMTO) and in the lab. In each case we present a detailed description of the setup. We discovered nonideal behavior in the system. This is at least in part due to properties of the VPM wire grids (diameter, spacing) employed in our experiment. Despite this, we found that the dual-VPM system is robust, operating with high efficiency and low instrumental polarization. This device is well suited for air and space-borne applications.

The simulated Doppler shifts of the solar Mg I Fraunhofer line produced by scattering on the solar light by asteroidal, cometary, and trans-Neptunian dust particles are compared with the shifts obtained by Wisconsin H-Alpha Mapper (WHAM) spectrometer. The simulated spectra are based on the results of integrations of the orbital evolution of particles. The deviation of the derived spectral parameters for various sources of dust used in the model reached maximum at the elongation (measured eastward from the Sun) between 90 deg and 120 deg. For the future zodiacal light Doppler shifts measurements, it is important to pay a particular attention to observing at this elongation range. At the elongations of the fields observed by WHAM, the model-predicted Doppler shifts were close to each other for several scattering functions considered. Therefore the main conclusions of our paper don't depend on a scattering function and mass distribution of particles if they are reasonable. A comparison of the dependencies of the Doppler shifts on solar elongation and the mean width of the Mg I line modeled for different sources of dust with those obtained from the WHAM observations shows that the fraction of cometary particles in zodiacal dust is significant and can be dominant. Cometary particles originating inside Jupiter's orbit and particles originating beyond Jupiter's orbit (including trans-Neptunian dust particles) can contribute to zodiacal dust about 1/3 each, with a possible deviation from 1/3 up to 0.1-0.2. The fraction of asteroidal dust is estimated to be about 0.3-0.5. The mean eccentricities of zodiacal particles located at 1-2 AU from the Sun that better fit the WHAM observations are between 0.2 and 0.5, with a more probable value of about 0.3.

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.

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

Filled arrays of bolometers are currently being employed for use in astronomy from the far-infrared through millimeter parts of the electromagnetic spectrum. Because of the large range of wavelengths for which such detectors are applicable, the number of modes supported by a pixel will vary according to the specific application of a given available technology. We study the dependence of image fidelity and induced polarization on the size of the pixel by employing a formalism in which diffraction due to the pixel boundary is treated by propagating the second-order statistical correlations of the radiation field through a model optical system. We construct polarized beam pattern images of square pixels for various ratios of p/\lambda where p is the pixel size and \lambda is the wavelength of the radiation under consideration. For the limit in which few modes are supported by the pixel (p/\lambda<1), we find that the diffraction due to the pixel edges is non-negligible and hence must be considered along with the telescope diffraction pattern in modeling the ultimate spatial resolution of an imaging system. For the case in which the pixel is over-moded (p/\lambda>1), the geometric limit is approached as expected. This technique gives a quantitative approach to optimize the imaging properties of arrays of planar detectors in the few-mode limit.

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

The Primordial Anisotropy Polarization Pathfinder Array (PAPPA) is a balloon-based instrument to measure the polarization of the cosmic microwave background and search for the signal from gravity waves excited during an inflationary epoch in the early universe. PAPPA will survey a 20 x 20 deg patch at the North Celestial Pole using 32 pixels in 3 passbands centered at 89, 212, and 302 GHz. Each pixel uses MEMS switches in a superconducting microstrip transmission line to combine the phase modulation techniques used in radio astronomy with the sensitivity of transition-edge superconducting bolometers. Each switched circuit modulates the incident polarization on a single detector, allowing nearly instantaneous characterization of the Stokes I, Q, and U parameters. We describe the instrument design and status.

We present simulated observations of the Doppler shifts of the solar Mg I Fraunhofer line scattered by asteroidal, cometary, and trans-Neptunian dust particles. The studies are based on the results of integrations of orbital evolution of particles under the gravitational influence of planets, the Poynting-Robertson drag, radiation pressure, and solar wind drag. The derived shifts in the centroid and profile of the line with solar elongation are different for different sources of dust. A comparison of the velocities of zodiacal dust particles based on these numerical integrations with the velocities obtained from WHAM observations shows that the fraction of cometary dust particles among zodiacal dust particles is significant and can be dominant. A considerable fraction of trans-Neptunian dust particles among zodiacal dust particles also fits different observations. The mean eccentricity of zodiacal dust particles is estimated to be about 0.5.

We report on new observations of the motion of zodiacal dust using optical absorption line spectroscopy of zodiacal light. We have measured the change in the profile shape of the scattered solar Mg I 5184 line toward several lines of sight in the ecliptic plane as well as the ecliptic pole. The variation in line centroid and line width as a function of helio-ecliptic longitude show a clear prograde signature and suggest that significant fraction of the dust follows non-circular orbits that are not confined to the ecliptic plane. When combined with dynamical models, the data suggest that the zodiacal dust is largely cometary, rather than asteroidal, in origin.

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 Single Aperture Far-InfraRed (SAFIR) Observatory's science goals are driven by the fact that the earliest stages of almost all phenomena in the universe are shrouded in absorption by and emission from cool dust and gas that emits strongly in the far-infrared and submillimeter. Over the past several years, there has been an increasing recognition of the critical importance of this spectral region to addressing fundamental astrophysical problems, ranging from cosmological questions to understanding how our own Solar System came into being. The development of large, far-infrared telescopes in space has become more feasible with the combination of developments for the James Webb Space Telescope and of enabling breakthroughs in detector technology. We have developed a preliminary but comprehensive mission concept for SAFIR, as a 10 m-class far-infrared and submillimeter observatory that would begin development later in this decade to meet the needs outlined above. Its operating temperature (<4K) and instrument complement would be optimized to reach the natural sky confusion limit in the far-infrared with diffraction-limited peformance down to at least 40 microns. This would provide a point source sensitivity improvement of several orders of magnitude over that of Spitzer or Herschel, with finer angular resolution, enabling imaging and spectroscopic studies of individual galaxies in the early universe. We have considered many aspects of the SAFIR mission, including the telescope technology, detector needs and technologies, cooling method and required technology developments, attitude and pointing, power systems, launch vehicle, and mission operations. The most challenging requirements for this mission are operating temperature and aperture size of the telescope, and the development of detector arrays.

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

We have resolved the ultraluminous infrared galaxy (ULIRG), IRAS 09111-1007, with the new 350 micron-optimised Second Generation Submillimeter High Angular Resolution Camera (SHARC II) and present the first submillimetre fluxes and images for the system. IRAS 09111-1007 comprises two interacting luminous infrared galaxies (LIRGs) with a projected nuclear separation of 39 kpc. The Western galaxy is roughly four times more luminous in the submillimetre than its Eastern counterpart. It is an extremely bright LIRG with an AGN. The classification of the Eastern source is uncertain: it could be a Seyfert 2 galaxy or a LINER. We highlight IRAS 09111-1007 as a system that necessitates further study: a double AGN ULIRG whose molecular gas content differs from other widely separated pairs and whose ULIRG phase might not be explained by current multiple merger and/or final stage ULIRG scenarios.

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.

Using the Wisconsin H-Alpha Mapper (WHAM), we have measured at high spectral resolution and high signal-to-noise the profile of the scattered solar Mg I 5184 absorption line in the zodiacal light. The observations were carried out toward 49 directions that sampled the ecliptic equator from solar elongations of 48\dg (evening sky) to 334\dg (morning sky) plus observations near +47\dg and +90\dg ecliptic latitude. The spectra show a clear prograde kinematic signature that is inconsistent with dust confined to the ecliptic plane and in circular orbits influenced only by the sun's gravity. In particular, the broadened widths of the profiles, together with large amplitude variations in the centroid velocity with elongation angle, indicate that a significant population of dust is on eccentric orbits. In addition, the wide, flat-bottomed line profile toward the ecliptic pole indicates a broad distribution of orbital inclinations extending up to about 30\dg - 40\dg with respect to the ecliptic plane. The absence of pronounced asymmetries in the shape of the profiles limits the retrograde population to less than 10% of the prograde population and also places constraints on the scattering phase function of the particles. These results do not show the radial outflow or evening--morning velocity amplitude asymmetry reported in some earlier investigations. The reduction of the spectra included the discovery and removal of extremely faint, unidentified terrestrial emission lines that contaminate and distort the underlying Mg I profile. This atmospheric emission is too weak to have been seen in earlier, lower signal-to-noise observations, but it probably affected the line centroid measurements of previous investigations.

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.

To understand the cosmic history of element synthesis it will be important to obtain extinction-free measures of the heavy element contents of high-redshift objects and to chart two monumental events: the collapse of the first metal-free clouds to form stars, and the initial seeding of the universe with dust. The information needed to achieve these objectives is uniquely available in the far-infrared/submillimeter (FIR/SMM) spectral region. Following the Decadal Report and anticipating the development of the Single Aperture Far-IR (SAFIR) telescope and FIR/SMM interferometry, we estimate the measurement capabilities of a large-aperture, background-limited FIR/SMM observatory and an interferometer on a boom, and discuss how such instruments could be used to measure the element synthesis history of the universe.

A high spectral resolution observation of the diffuse X-ray background in the 60 - 1000 eV energy range has been made using an array of thirty-six 1 mm^2 micro-calorimeters flown on a sounding rocket. Detector energy resolution ranged from 5-12 eV FWHM, and a composite spectrum of ~ 1 steradian of the background centered at l = 90, b = +60 was obtained with a net resolution of ~ 9 eV. The target area includes bright 1/4 keV regions, but avoids Loop I and the North Polar Spur. Lines of C VI, O VII, and O VIII are clearly detected with intensities of 5.4 +/- 2.3, 4.8 +/- 0.8, and 1.6 +/- 0.4 photons cm^-2 s^-1 sr^-1, respectively. The oxygen lines alone account for a majority of the diffuse background observed in the ROSAT R4 band that is not due to resolved extragalactic discrete sources. We also have a positive detection of the Fe-M line complex near 70 eV at an intensity consistent with previous upper limits that indicate substantial gas phase depletion of iron. We include a detailed description of the instrument and its detectors.

The question "How did we get here and what will the future bring?" captures the human imagination and the attention of the National Academy of Science's Astronomy and Astrophysics Survey Commitee (AASC). Fulfillment of this "fundamental goal" requires astronomers to have sensitive, high angular and spectral resolution observations in the far-infrared/submillimeter (far-IR/sub-mm) spectral region. With half the luminosity of the universe and vital information about galaxy, star and planet formation, observations in this spectral region require capabilities similar to those currently available or planned at shorter wavelengths. In this paper we summarize the scientific motivation, some mission concepts and technology requirements for far-IR/sub-mm space interferometers that can be developed in the 2010-2020 timeframe.

We discuss the utility of a low resolution prism as a component of a Multi-Object Spectrometer for NASA's proposed Next Generation Space Telescope (NGST). Low resolution prism spectroscopy permits simultaneous observation of the 0.6-5micron wavelength regime at R~50. To such data we can apply the modern techniques in spectral energy distribution (SED) fitting to determine source redshifts, sometimes called ``photometric redshifts''. Low resolution prism observations of galaxy SED's provide a significant advantage over multi-filter observations for any realistic observing strategy, and we have performed extensive simulations to quantify it. The advantage of the prism is largest in the regime of faint (K_AB>30) objects at high redshift (z>4), the primary discovery space of NGST. Many important objects will be too faint for follow up at higher spectral resolution, so prism observations are the optimal technique to study them. Prism observations also reduce the contamination of high redshift samples by lower redshift interlopers.

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.

The advent of submillimeter wavelength array cameras operating on large ground-based telescopes is revolutionizing imaging at these wavelengths, enabling high-resolution submillimeter surveys of dust emission in star-forming regions and galaxies. Here we present a recent 350 micron image of the edge-on galaxy NGC 891, which was obtained with the Submillimeter High Angular Resolution Camera (SHARC) at the Caltech Submillimeter Observatory (CSO). We find that high resolution submillimeter data is a vital complement to shorter wavelength satellite data, which enables a reliable separation of the cold dust component seen at millimeter wavelengths from the warmer component which dominates the far-infrared (FIR) luminosity.

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.

The FIRAS instrument on the COBE satellite has conducted an unbiased survey of the far-infrared emission from our Galaxy. The first results of this survey were reported by Wright et al. (1991). We report the results of new analyses of this spectral survey, which includes emission lines from 158 um C+, 122 um and 205 um N+, 370 um and 609 um C, and CO J=2-1 through 5-4. We report the morphological distribution along the galactic plane (b=0) of the spectral line emission, and the high galactic latitude intensities of the C+ and 205 um N+ emission. The high galactic latitude intensity cosecant of the 158 um fine structure transition from C+ is presented, and C+ is seen to decrease more rapidly than the far infrared intensity with increasing galactic latitude. C+ and H I emission are closely correlated with a C+ cooling rate of (2.65 +/- 0.15)x10^{-26} erg/s/H-atom. We conclude that this emission arises almost entirely from the Cold Neutral Medium. The high galactic latitude intensity of the 205 um fine structure transition from N+ arises entirely from the Warm Ionized Medium, and its cosecant is presented. We estimate the total ionizing photon rate in the Galaxy to be 3.5x10^{53} ionizing photons per second, based on the 205 um N+ transition.