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

In this work we incorporate the newest ISO results on the mid-infrared spectral-energy-distributions (MIR SEDs) of galaxies into models for the number counts and redshift distributions of MIR surveys. A three-component model, with empirically determined MIR SED templates of (1) a cirrus/PDR component (2) a starburst component and (3) an AGN component, is developed for infrared (3--120\micron) SEDs of galaxies. The model includes a complete IRAS 25\micron selected sample of 1406 local galaxies ($z \leq 0.1$; Shupe et al. 1998a). Results based on these 1406 spectra show that the MIR emission features cause significant effects on the redshift dependence of the K-corrections for fluxes in the WIRE 25\micron band and ISOCAM 15\micron band. This in turn will affect deep counts and redshift distributions in these two bands, as shown by the predictions of two evolution models (a luminosity evolution model with $L\propto (1+z)^3$ and a density evolution model with $\rho\propto (1+z)^4$). The dips-and-bumps on curves of MIR number counts, caused by the emission features, should be useful indicators of evolution mode. The strong emission features at $\sim 6$--8\micron will help the detections of relatively high redshift ($z\sim 2$) galaxies in MIR surveys. On the other hand, determinations of the evolutionary rate based on the slope of source counts, and studies on the large scale structures using the redshift distribution of MIR sources, will have to treat the effects of the MIR emission features carefully. We have also estimated a 15\micron local luminosity function from the predicted 15\micron fluxes of the 1406 galaxies using the bivariate (15\micron vs. 25\micron luminosities) method. This luminosity function will improve our understanding of the ISOCAM 15\micron surveys.

We have established a model to systematically estimate the contribution of the mid-infrared emission features between 3 and 11.6 micron to the IRAS in-band fluxes, using the results of ISO PHT-S observation of 16 galaxies by Lu et al. (1997). The model is used to estimate more properly the k-corrections for calculating the restframe 12 and 25 micron fluxes and luminosities of IRAS galaxies. We have studied the 12-25 micron color-luminosity relation for a sample of galaxies selected at 25 microns. The color is found to correlate well with the 25 micron luminosity, the far-infrared and the blue luminosities. The relations with the mid-infrared luminosities are more sensitive to different populations of galaxies, while a single relation of the 12-25 micron color vs. the ratio of the far-infrared and the blue luminosities applies equally well to these different populations. The luminous and ultraluminous infrared galaxies have redder 12-25 micron colors than those of the quasars. These relations provide powerful tools to differentiate different populations of galaxies. We have also selected a sample of 668 galaxies from the IRAS Faint Source Survey flux-density limited at 200 mJy at 12 microns. A 12 micron local luminosity function is derived and, for the first time in the literature, effects of density variation in the local universe are considered and corrected in the calculation of the 12 micron luminosity function. It is also found that the 12 micron-selected sample are dominated by quasars and active galaxies, which therefore strongly affect the 12 micron luminosity function at high luminosities. The ultraluminous infrared galaxies are relatively rare at 12 micron comparing with a 25 micron sample.

The local luminosity function at 25 $\mu$m provides the basis for interpreting the results of deep mid-infrared surveys planned or in progress with space astrophysics missions including ISO, WIRE and SIRTF. We have selected a sample of 1458 galaxies from the IRAS Faint Source Survey with a flux density limit of 250 mJy at 25 $\mu$m. The local luminosity function is derived using both parametric and non-parametric maximum-likelihood techniques, and the classical $1/V_{max}$ estimator. Comparison of these results shows that the $1/V_{max}$ estimate of the luminosity function is significantly affected by the Local Supercluster. A maximum-likelihood fit to the radial density shows no systematic increase that would be caused by density evolution of the galaxy population. The density fit is used to correct the $1/V_{max}$ estimate. We also demonstrate the high quality and completeness of our sample by a variety of methods. The luminosity function derived from this sample is compared to previously published estimates, showing the prior estimates to have been strongly affected by the Local Supercluster. Our new luminosity function leads to lower estimates of mid-infrared backgrounds and number counts.