The Stratospheric Observatory for Infrared Astronomy (SOFIA) has recently concluded a set of engineering flights for Observatory performance evaluation. These in-flight opportunities are viewed as a first comprehensive assessment of the Observatory's performance and are used to guide future development activities, as well as to identify additional Observatory upgrades. Pointing stability was evaluated, including the image motion due to rigid-body and flexible-body telescope modes as well as possible aero-optical image motion. We report on recent improvements in pointing stability by using an active mass damper system installed on the telescope. Measurements and characterization of the shear layer and cavity seeing, as well as image quality evaluation as a function of wavelength have also been performed. Additional tests targeted basic Observatory capabilities and requirements, including pointing accuracy, chopper evaluation and imager sensitivity. This paper reports on the data collected during these flights and presents current SOFIA Observatory performance and characterization.

SOFIA is a 2.5 meter airborne infrared telescope, mounted in a Boeing 747SP aircraft. Due to the large size of the telescope, only a few degrees of azimuth are available at the telescope bearing. This means the heading of the aircraft is fundamentally associated with the telescope's observation targets, and the ground track necessary to enable a given mission is highly complex and dependent on the coordinates, duration, and order of observations to be performed. We have designed and implemented a Flight Management Infrastructure (FMI) product in order to plan and execute such missions in the presence of a large number of external constraints (e.g. restricted airspace, international boundaries, elevation limits of the telescope, aircraft performance, winds at altitude, and ambient temperatures). We present an overview of the FMI, including the process, constraints and basic algorithms used to plan and execute SOFIA missions.

SOFIA is an airborne, gyroscopically stabilized 2.5m infrared telescope, mounted to a spherical bearing. Unlike its predecessors, SOFIA will work in absolute coordinates, despite its continually changing position and attitude. In order to manage this, SOFIA must relate equatorial and telescope coordinates using a combination of avionics data and star identification, manage field rotation and track sky images. We describe the algorithms and systems required to acquire and maintain the equatorial reference frame, relate it to tracking imagers and the science instrument, set up the oscillating secondary mirror, and aggregate pointings into relocatable nods and dithers.

Cosmological dark matter in the form of neutrinos with masses of up to a few electron volts is known as hot dark matter. After an historical review of the subject, this article considers constraints on hot dark matter from current data on neutrino oscillations and on cosmology. The atmospheric neutrino oscillation data imply a lower limit on the HDM contribution to the cosmological density $\Omega_\nu \gsim 0.001$. The possible improvement of low-$\Omega_m$ flat ($\Lambda$CDM) cosmological models with the addition of light neutrinos appears to be rather limited, but$\Lambda$CDM models with $\Omega_\nu \lsim 0.1$ may be consistent with presently available data. Data expected soon may permit detection of such a hot dark matter contribution, or alternatively provide stronger upper limits on $\Omega_\nu$ and neutrino masses.

We review our analysis of the clustering properties of ``Lyman-break'' galaxies (LBGs) at redshift z~3, previously discussed in Wechsler et al (1998). We examine the likelihood of spikes found by Steidel et al (1998) in the redshift distribution of LBGs, within a suite of models for the evolution of structure in the Universe. Using high-resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys from these models in the same way that they are actually observed, identifying LBGs with the most massive dark matter halos. We find that all the models (with SCDM as a marginal exception) have a substantial probability of producing spikes similar to those observed, because the massive halos are much more clumped than the underlying matter -- i.e., they are biased. Therefore, the likelihood of such a spike is not a good discriminator among these models. The LBG correlation functions are less steep than galaxies today (gamma~1.4), but show similar or slightly longer correlation lengths. We have extened this analysis and include a preliminary comparison to the new data presented in Adelberger et al (1998). We also discuss work in progress, in which we use semi-analytic models to identify Lyman-break galaxies within dark-matter halos.

The recent atmospheric neutrino data from Super-Kamiokande provide strong evidence of neutrino oscillations and therefore of non-zero neutrino mass. These data imply a lower limit on the hot dark matter (i.e., light neutrino) contribution to the cosmological density $\Omega_\nu >~ 0.001$ --- almost as much as that of all the stars in the universe --- and permit higher $\Omega_\nu$. The ``standard'' COBE-normalized critical-matter-density (i.e., $\Omega_m=1$) Cold Dark Matter (CDM) model has too much power on small scales. But adding to CDM neutrinos with mass of about 5 eV, corresponding to $\Omega_\nu \approx 0.2$, results in a much improved fit to data on the nearby galaxy and cluster distribution. Indeed, the resulting Cold + Hot Dark Matter (CHDM) cosmological model is arguably the most successful $\Omega_m=1$ model for structure formation. However, other recent data has begun to make the case for $\Omega_m <~ 0.6$ fairly convincing. In light of all this new data, we reconsider whether cosmology still provides evidence favoring neutrino mass of a few eV in flat models with cosmological constant $\Omega_\Lambda = 1 - \Omega_m$. We find that the possible improvement of the low-$\Omega_m$ flat (LCDM) cosmological models with the addition of light neutrinos appears to be rather limited.

We present two matched sets of five simulations each, covering five presently favored simple modifications to the standard cold dark matter (CDM) scenario. One simulation suite, with a linear box size of 75 Mpc/h, is designed for high resolution and good statistics on the group/poor cluster scale, and the other, with a box size of 300 Mpc/h, is designed for good rich cluster statistics. All runs had 57 million cold particles, and models with massive neutrinos had an additional 113 million hot particles. We consider separately models with massive neutrinos, tilt, curvature, and a nonzero cosmological constant in addition to the standard CDM model. We find that our tilted Omega+Omega_Lambda=1 (TLCDM) model produces too much small-scale power by a factor of ~3, and our open Lambda=0 (OCDM) model also exceeds observed small-scale power by a factor of 2. In addition, we take advantage of the large dynamic range in detectable halo masses our simulations allow to check the shape of the Press-Schechter approximation. We find good fits at cluster masses for delta_c=1.27--1.35 for a Gaussian filter and delta_c=1.57--1.73 for a tophat filter. However, Press-Schechter overpredicts the number density of halos compared to the simulations in the high resolution suite by a weakly cosmology-dependent factor of 1.5--2 at galaxy and group masses, which cannot be fixed by adjusting delta_c within reasonable bounds. An appendix generalizes the spherical collapse model to any isotropic cosmology.

We address the high peaks found by Steidel et al (1997) in the redshift distribution of ``Lyman-break'' objects (LBOs) at redshift z~3. The highest spike represents a relative overdensity of 2.6 in the distribution of LBOs in pixels of comoving size ~10Mpc/h. We examine the likelihood of such a spike in the redshift distribution within a suite of models for the evolution of structure in the Universe, including models with Omega=1 (SCDM and CHDM) and with Omega=0.4-0.5 (LCDM and OCDM). Using high-resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys from these models in the same way that they are actually observed, identifying LBOs with the most massive dark matter halos. We find that all the models (with SCDM as a marginal exception) have a substantial probability of producing spikes similar to those observed, because the massive halos are much more clumped than the underlying matter -- i.e., they are biased. Therefore, the likelihood of such a spike is not a good discriminator among these models. We find in these models that the mean biasing parameter b of LBOs with respect to dark matter varies within a range b ~2-5 on a scale of ~10Mpc/h. We also compute the two-body correlation functions of LBOs predicted in these models. The LBO correlation functions are less steep than galaxies today (gamma ~1.4), but show similar or slightly longer correlation lengths.

The combination of measurements of the local abundance of rich clusters of galaxies and its evolution to higher redshift offers the possibility of a direct measurement of $\Omega_0$ with little contribution from other cosmological parameters. We investigate the significance of recent claims that this evolution indicates that $\Omega_0$ must be small. The most recent cluster velocity dispersion function from a compilation including the ESO Northern Abell Cluster Survey (ENACS) results in a significantly higher normalization for models, corresponding to $\sigma_8\approx 0.6$ for $\Omega_0=1$, compared to the Eke, Cole, & Frenk result of $\sigma_8=0.52\pm 0.04$. Using the ENACS data for a $z=0$ calibration results in strong evolution in the abundance of clusters, and we find that the velocity dispersion function is consistent with $\Omega_0=1$. The results are dependent upon the choice and analysis of low-redshift and high-redshift data, so at present, the data is not good enough to determine $\Omega_0$ unambiguously.

Reconstruction of the linear power spectrum from observational data provides a way to compare cosmological models to a large amount of data, as Peacock & Dodds (1994, 1996) have shown. By applying the appropriate corrections to the observational power spectrum it is possible to recover the underlying linear power spectrum for any cosmological model. Using extensive N-body simulations we demonstrate that the method is applicable to a wide range of cosmological models. However, we find that the recovery of the linear power spectrum from observations following PD94 is misleading because the corrections are model- dependent. When we apply the proper corrections for a given model to the observational power spectrum, we find that no model in our test group recovers the linear power spectrum well for the bias suggested by PD94 between Abell, Radio, Optical, and IRAS catalogs 4.5:1.9:1.3:1, with b_IRAS=1. When we allow b_IRAS to vary we find that: (i)CHDM models give very good fits to observations if optically-selected galaxies are slightly biased b_Opt=1.1 (ii) Most LCDM models give worse but acceptable fits if blue galaxies are considerably antibiased: 0.6<b_Opt<0.9 and fail if optical galaxies are biased. (iii)There is a universal shape of the recovered linear power spectrum of all LCDM models over their entire range of explored wavenumbers,0.01<k<0.6h\Mpc. Recovered spectra of CDM and CHDM models are nearly the same as that of LCDM in the region 0.01<k<0.2h/Mpc but diverge from this spectrum at higher k.

The year 1996 will mark the initiation of a number of new space missions to the planet Mars from which we expect to obtain a rich set of data, including spacecraft radio tracking data. Anticipating these events, we have analyzed the feasibility of testing a violation of the strong equivalence principle (SEP) with Earth-Mars ranging. Using analytic and numerical methods, we have demonstrated that ranging data can provide a useful estimate of the SEP parameter $\eta$. Two estimates of the predicted accuracy are quoted, one based on conventional covariance analysis, and the other based on ``modified worst case'' analysis, which assumes that systematic errors dominate the experiment. If future Mars missions provide ranging measurements with an accuracy of $\sigma$ meters, after ten years of ranging the expected accuracy for the parameter $\eta$ will be of order $\sigma_\eta\approx (1-12)\times 10^{-4}\sigma $. In addition, these ranging measurements will provide a significantly improved determination of the mass of the Jupiter system, independent of the test of the {\small SEP} polarization effect.

The Earth, Mars, Sun, Jupiter system allows for a sensitive test of the strong equivalence principle (SEP) which is qualitatively different from that provided by Lunar Laser Ranging. Using analytic and numerical methods we demonstrate that Earth-Mars ranging can provide a useful estimate of the SEP parameter $\eta$. Two estimates of the predicted accuracy are derived and quoted, one based on conventional covariance analysis, and another (called ``modified worst case'' analysis) which assumes that systematic errors dominate the experiment. If future Mars missions provide ranging measurements with an accuracy of $\sigma$ meters, after ten years of ranging the expected accuracy for the SEP parameter $\eta$ will be of order $(1-12)\times 10^{-4}\sigma$. These ranging measurements will also provide the most accurate determination of the mass of Jupiter, independent of the SEP effect test.