We describe the operations concept and data reduction plan for the Mid- Infrared Instrument (MIRI) for the James Webb Space Telescope (JWST). The overall JWST operations concept is to use Observation Templates (OTs) to provide a straightforward and intuitive way for users to specify observations. MIRI has four OTs that correspond to the four observing modes: 1.) Imaging, 2.) Coronagraphy, 3.) Low Resolution Spectroscopy, and 4.) Medium Resolution Spectroscopy. We outline the user choices and expansion of these choices into detailed instrument operations. The data reduction plans for MIRI are split into three stages, where the specificity of the reduction steps to the observation type increases with stage. The reduction starts with integration ramps: stage 1 yields uncalibrated slope images; stage 2 calibrates the slope images; and then stage 3 combines multiple calibrated slope images into high level data products (e.g. mosaics, spectral cubes, and extracted source information). Finally, we give examples of the data and data products that will be derived from each of the four different OTs.

The MIRI Si:As IBC detector arrays extend the heritage technology from the Spitzer IRAC arrays to a 1024 x 1024 pixel format. We provide a short discussion of the principles of operation, design, and performance of the individual MIRI detectors, in support of a description of their operation in arrays provided in an accompanying paper (Ressler et al. (2015)). We then describe modeling of their response. We find that electron diffusion is an important component of their performance, although it was omitted in previous models. Our new model will let us optimize the bias voltage while avoiding avalanche gain. It also predicts the fraction of the IR-active layer that is depleted (and thus contributes to the quantum efficiency) as signal is accumulated on the array amplifier. Another set of models accurately predicts the nonlinearity of the detector-amplifier unit and has guided determination of the corrections for nonlinearity. Finally, we discuss how diffraction at the interpixel gaps and total internal reflection can produce the extended cross-like artifacts around images with these arrays at short wavelengths, ~ 5 microns. The modeling of the behavior of these devices is helping optimize how we operate them and also providing inputs to the development of the data pipeline.

We present deep HST/NICMOS Camera 3 F110W and F160W imaging of a 26'x33', corresponding to 3.1pcx3.8pc, non-contiguous field towards the Orion Nebula Cluster (ONC). The main aim is to determine the ratio of low--mass stars to brown dwarfs for the cluster as a function of radius out to a radial distance of 1.5pc. The sensitivity of the data outside the nebulous central region is F160W=21.0 mag, significantly deeper than previous studies of the region over a comparable area. We create an extinction limited sample and determine the ratio of low-mass stars (0.08-1Msun) to brown dwarfs (0.02-0.08Msun and 0.03-0.08Msun) for the cluster as a whole and for several annuli. The ratio found for the cluster within a radius of 1.5pc is R(02)=N(0.08-1Msun)/N(0.02-0.08Msun)=1.7+-0.2, and R(03)=N(0.08-1Msun)/N(0.03-0.08Msun)=2.4+-0.2, after correcting for field stars. The ratio for the central 0.3pcx0.3pc region down to 0.03Msun was previously found to be R(03)=3.3+0.8-0.7, suggesting the low-mass content of the cluster is mass segregated. We discuss the implications of a gradient in the ratio of stars to brown dwarfs in the ONC in the context of previous measurements of the cluster and for other nearby star forming regions. We further discuss the current evidence for variations in the low-mass IMF and primordial mass segregation.

[Abridged] We detect three (plus one less certain) z-dropout sources in two separate fields of our UDF05 HST NICMOS images. These z~7 Lyman-Break Galaxy (LBG) candidates allow us to constrain the Luminosity Function (LF) of the star forming galaxy population at those epochs. By assuming a change in only M* and adopting a linear evolution in redshift, anchored to the measured values at z~6, the best fit evolution coefficient is found to be 0.43+-0.19 mag per unit redshift (0.36+-0.18, if including all four candidates), which provides a value of M*(z=7.2)=-19.7+-0.3. This implies a steady evolution for the LBG LF out to z~7, at the same rate that is observed throughout the z~3 to 6 period. This puts a strong constraint on the star-formation histories of z~6 galaxies, whose ensemble star-formation rate density must be lower by a factor 2 at ~170 Myr before the epoch at which they are observed. In particular, a large fraction of stars in the z~6 LBG population must form at redshifts well above z~7. Extrapolating this steady evolution of the LF out to higher redshifts, we estimate that galaxies would be able to reionize the universe by z~6, provided that the faint-end slope of the z>7 LF steepens to alpha~-1.9, and that faint galaxies, with luminosities below the current detection limits, contribute a substantial fraction of the required ionizing photons. This scenario gives however an integrated optical depth to electron scattering that is ~2sigma below the WMAP-5 measurement. Therefore, altogether, our results indicate that, should galaxies be the primary contributors to reionization, either the currently detected evolution of the galaxy population slows down at z>7, or the LF evolution must be compensated by a decrease in metallicity and a corresponding increase in ionization efficiency at these early epochs.

We present the UDF05 project, a HST Large Program of deep ACS (F606W, F775W, F850LP) and NICMOS (F110W, F160W) imaging of three fields, two of which coincide with the NICP1-4 NICMOS parallel observations of the Hubble Ultra Deep Field (HUDF). In this first paper we use the ACS data for the NICP12 field, as well as the original HUDF ACS data, to measure the UV Luminosity Function (LF) of z~5 Lyman Break Galaxies (LBGs) down to very faint levels. Specifically, based on a V-i, i-z selection criterion, we identify a sample of 101 and 133 candidate z~5 galaxies down to z_{850}=28.5 and 29.25 magnitudes in the NICP12 and in the HUDF fields, respectively. Using an extensive set of Monte Carlo simulations we derive corrections for observational biases and selection effects, and construct the rest-frame 1400 A LBG LF over the range M_{1400}=[-21.4, -17.1], i.e. down to ~0.04 L* at z~5, and complement it with data from the Subaru Deep Field (SDF) from Yoshida et al. (2006) to extend it to the brighter end (M_{1400}>-22.2). We show that: (i) Different assumptions regarding the SED distribution of the LBG population, dust properties and intergalactic absorption result in a 25% variation in the number density of LBGs at z~5; (ii) Under consistent assumptions for dust properties and intergalactic absorption, the HUDF is about 30% under-dense in z~5 LBGs relative to the NICP12 field, a variation which is well explained by cosmic variance; (iii) The faint-end slope of the LF is independent of the specific assumptions for the input physical parameters, and has a value of alpha ~ -1.6, similar to the faint-end slope of the LF that has been measured for LBGs at z~3 and z~6.

This paper presents the Hubble Ultra Deep Field (HUDF), a one million second exposure of an 11 square minute-of-arc region in the southern sky with the Hubble Space Telescope. The exposure time was divided among four filters, F435W (B435), F606W (V606), F775W (i775), and F850LP (z850), to give approximately uniform limiting magnitudes mAB~29 for point sources. The image contains at least 10,000 objects presented here as a catalog. Few if any galaxies at redshifts greater than ~4 resemble present day spiral or elliptical galaxies. Using the Lyman break dropout method, we find 504 B-dropouts, 204 V-dropouts, and 54 i-dropouts. Using these samples that are at different redshifts but derived from the same data, we find no evidence for a change in the characteristic luminosity of galaxies but some evidence for a decrease in their number densities between redshifts of 4 and 7. The ultraviolet luminosity density of these samples is dominated by galaxies fainter than the characteristic luminosity, and the HUDF reveals considerably more luminosity than shallower surveys. The apparent ultraviolet luminosity density of galaxies appears to decrease from redshifts of a few to redshifts greater than 6. The highest redshift samples show that star formation was already vigorous at the earliest epochs that galaxies have been observed, less than one billion years after the Big Bang.

We provide an overview of the most important calibration aspects of the NICMOS instrument on board of HST. We describe the performance of the instrument after the installation of the NICMOS Cooling System, and show that the behavior of the instrument has become very stable and predictable. We detail the improvements made to the NICMOS pipeline and outline plans for future developments. The derivation of the absolute photometric zero-point calibration is described in detail. Finally, we describe and quantify a newly discovered count-rate dependent non-linearity in the NICMOS cameras. This new non-linearity is distinctly different from the total count dependent non-linearity that is well known for near-infrared detectors. We show that the non-linearity has a power law behavior, with pixels with high system, or vice versa, pixels with low count rate detecting slightly less than expected. The effect has a wavelength dependence with observations at the shortest wavelengths being the most affected (~0.05-0.1 mag per dex flux change at ~1 micron, 0.03 mag per dex at 1.6 micron).

We present grism spectra taken with the Advanced Camera for Surveys to identify 29 red sources with (i'-z') > 0.9 in the Hubble Ultra Deep Field (HUDF). Of these 23 are found to be galaxies at redshifts between z=5.4 and 6.7, identified by the break at 1216 angstrom due to IGM absorption; two are late type dwarf stars with red colors; and four are galaxies with colors and spectral shape similar to dust reddened or old galaxies at redshifts z~1-2. This constitutes the largest uniform, flux-limited sample of spectroscopically confirmed galaxies at such faint fluxes z'<27.5. Many are also among the most distant spectroscopically confirmed galaxies (at redshifts up to z=6.7). We find a significant overdensity of galaxies at redshifts z=5.9 +/- 0.2. Nearly two thirds of the galaxies in our sample (15/23) belong to this peak. Taking into account the selection function and the redshift sensitivity of the survey, we get a conservative overdensity of at least a factor of two along the line-of-sight. The galaxies found in this redshift peak are also localized in the plane of the sky in a non-random manner, occupying about half of the ACS chip. Thus the volume overdensity is a factor of four. The star-formation rate derived from detected sources in this overdense region is sufficient to reionize the local IGM.

We have discovered an extended Lyman alpha plume associated with a compact source at redshift 5.4 in slitless spectroscopic data from the Grism ACS Program for Extragalactic Science (GRAPES) project. The spatial extent of the emission is about 6 x 1.5 kpc (1 x 0.25 arcsec). Combining our grism data and the broadband images from the Hubble UltraDeep Field (UDF) images, we find a Lyman alpha line flux of 2e-17 erg/cm2/s and surface brightness 7e-17 erg/cm2/s/arcsec2. The UDF images show diffuse continuum emission associated with the Lyman alpha plume (hereafter UDF 5225), with three embedded knots. The morphology of UDF 5225 is highly suggestive of a galaxy in assembly. It is moreover possible that the prominent Lyman alpha emission from this object is due to an active nucleus, and that we are seeing the simultaneous growth through accretion of a galaxy and its central black hole. Followup observations at higher spectral resolution could test this hypothesis.

We present a data set of images of the gravitationally lensed quasar Q2237+0305, that was obtained at the Apache Point Observatory (APO) between June 1995 and January 1998. Although the images were taken under variable, often poor seeing conditions and with coarse pixel sampling, photometry is possible for the two brighter quasar images A and B with the help of exact quasar image positions from HST observations. We obtain a light curve with 73 data points for each of the images A and B. There is evidence for a long (>~ 100 day) brightness peak in image A in 1996 with an amplitude of about 0.4 to 0.5 mag (relative to 1995), which indicates that microlensing has been taking place in the lensing galaxy. Image B does not vary much over the course of the observation period. The long, smooth variation of the light curve is similar to the results from the OGLE monitoring of the system (Wozniak et al. 2000a).

We present near-IR imaging polarimetry observations of the nucleus of Cygnus A, taken with the NICMOS camera of the HST at a wavelength of 2.0 microns.These maps reveal a highly collimated region of polarized emission straddling the nucleus and extending to a radius of 1.2 arcseconds. Remarkably, this feature coincides with one, but only one, limb of the edge-brightened bicone structure seen in the total intensity image. The high degree (P_k ~ 25%) and orientation of the extended polarization feature are consistent with a scattering origin. Most plausibly, the detection of polarization along only one limb of the bicone is a consequence of intrinsic anisotropy of the near-IR continuum within the radiation cones, with the direction of maximum intensity of the near-IR radiation field significantly displaced from the direction of the radio axis. The unresolved nuclear core source is also highly polarized (P_k > 28%), with a position angle close to the perpendicular to the radio axis. Given that this high degree of nuclear polarization can only be explained in terms of dichroic extinction if the dichroic mechanism is unusually efficient in Cygnus A, it is more likely that the nuclear polarization is due to the scattering of nuclear light in an unresolved scattering region close to the AGN. In this case, the flux of the core source in the K-band is dominated by scattered rather than transmitted quasar light, and previous extinction estimates based on K-band photometry of the core substantially underestimate the true nuclear extinction.

We report on observations of the fading optical counterpart of the gamma-ray burst GRB 970228, made with the Hubble Space Telescope STIS CCD approximately six months after outburst and with the HST/NICMOS and Keck/NIRC approximately one year after outburst. The unresolved counterpart is detected by STIS at V=28.0 +/- 0.25, consistent with a continued power-law decline with exponent -1.14 +/- 0.05. The counterpart is located within, but near the edge of, a faint extended source with diameter ~0."8 and integrated magnitude V=25.8 +/- 0.25. A reanalysis of HST and NTT observations performed shortly after the burst shows no evidence of proper motion of the point source or fading of the extended emission. Only the extended source is visible in the NICMOS images with a magnitude of H=23.3 +/- 0.1. The Keck observations find K = 22.8 +/- 0.3. Several distinct and independent means of deriving the foreground extinction in the direction of GRB 970228 all agree with A_V = 0.75 +/- 0.2. After adjusting for Galactic extinction, we find that the size of the observed extended emission is consistent with that of galaxies of comparable magnitude found in the Hubble Deep Field (HDF) and other deep HST images. Only 2% of the sky is covered by galaxies of similar or greater surface brightness; therefore the extended source is almost certainly the host galaxy. Additionally, we find that the extinction-corrected V - H and V - K colors of the host are as blue as any galaxy of comparable or brighter magnitude in the HDF. Taken in concert with recent observations of GRB 970508, GRB 971214, and GRB 980703 our work suggests that all four GRBs with spectroscopic identification or deep multicolor broad-band imaging of the host lie in rapidly star-forming galaxies.

We report on Hubble Space Telescope observations of the optical transient discovered in the error box of GRB970508. The object was imaged on 1997 June 2 with the Space Telescope Imaging Spectrograph and Near-Infrared Camera and Multi-Object Spectrometer. The observations reveal a point-like source with R = 23.1 +- 0.2 and H = 20.6 +- 0.3, in agreement with the power-law temporal decay seen in previous ground-based monitoring. Unlike the case of GRB970228, no nebulosity is detected surrounding the OT of GRB970508, although Mg I absorption and [O II] emission seen in Keck spectra at a redshift of z = 0.835 suggest the presence of a dense, star-forming medium. The HST observations set very conservative upper limits of R ~24.5 and H ~22.2 on the brightness of any underlying extended source. If this subtends a substantial fraction of an arcsecond, then the R band limit is ~25.5. Subsequent photometry suggests a flattening of the light curve at later epochs. Assuming the OT decline follows a pure power-law and ascribing the flattening to the presence of an underlying component of constant flux, we find that this must have R = 25.4, consistent with the upper limits determined by HST. At z = 0.8, this would correspond to an absolute magnitude in the U band of ~-18, similar to that of the Large Magellanic Cloud. We propose a scenario in which the host galaxy of the GRB is of Magellanic type, possibly being a ``satellite'' of one of the bright galaxies located at few arcseconds from the optical transient.

We report on Hubble Space Telescope (HST) observations of the optical transient (OT) discovered in the error box of the gamma-ray burst GRB970508. The object was imaged on 1997 June 2 with the Space Telescope Imaging Spectrograph (STIS) and Near-Infrared Camera and Multi-Object Spectrometer (NICMOS). The observations reveal a point-like source with R = 23.1 +- 0.2 and H = 20.6 +- 0.3, in agreement with the power-law temporal decay seen in ground-based monitoring. Unlike the case of GRB970228, no nebulosity is detected surrounding the OT of GRB970508. We set very conservative upper limits of R ~ 24.5 and H ~ 22.2 on the brightness of any underlying extended source. If this subtends a substantial fraction of an arcsecond, then the R band limit is ~25.5. In combination with Keck spectra that show Mg I absorption and [O II] emission at a redshift of z = 0.835, our observations suggest that the OT is located in a star-forming galaxy with total luminosity one order of magnitude lower than the knee of the galaxy luminosity function, L*. Such galaxies are now thought to harbor the majority of star formation at z ~ 1; therefore, these observations may provide support for a link between GRBs and star formation.

Photometric monitoring of the gravitational lens system 0957+561A,B in the g and r bands with the Apache Point Observatory (APO) 3.5 m telescope during 1996 shows a sharp g band event in the trailing (B) image light curve at the precise time predicted from the observation of an event during 1995 in the leading (A) image with a delay of 415 days. This success confirms the "short delay," and the lack of any feature at a delay near 540 days rejects the "long delay" for this system, resolving a long-standing controversy. A series of statistical analyses of our light curve data yield a best fit delay of 417 +/- 3 days (95% confidence interval). Recent improvements in the modeling of the lens system (consisting of a galaxy and cluster) allow us to derive a value of the global (at z = 0.36) value of Hubble's constant H_0 using Refsdal's method, a simple and direct distance determination based on securely understood physics and geometry. The result is H_0 = 63 +/- 12 km/s/Mpc (for Omega = 1) where this 95% confidence interval is dominated by remaining lens model uncertainties.