We present constraints on cosmological and star formation parameters based on combining observations of the Wilkinson Microwave Anisotropy Probe (WMAP) and high-redshift quasars from the Sloan Digital Sky Survey (SDSS). We use a semi-analytic model for reionization (Chiu and Ostriker 2000) that takes into account a number of important physical processes both within collapsing halos and in the intergalactic medium. Assuming that the efficiency of producing UV photons per baryon is constant, we derive a constraint of the form sigma_8 Omega_0^0.5~0.33 in a flat, Lambda-dominated universe with h=0.72, n=0.99, and Omega_b h^2=0.024. However, the calculated optical depth to electron scattering of tau_es~0.06 is well below the value found by WMAP of 0.17+/-(0.04~0.07) (Spergel et al 2003). Since the WMAP constraints on tau_es are somewhat degenerate with the value of the spectral index n, we then permit the primordial spectral index n to float and fix Omega_0 h^2=0.14, while normalizing the power spectrum using WMAP. In addition, we allow the UV-efficiency to have time-dependence. Combining the WMAP constraints with the quasar transmission data, our analysis then favors a model with tau_es=0.11^{+0.02}_{-0.03}, n=0.96^{+0.02}_{-0.03}$, implying sigma_8=0.83^{+0.03}_{-0.05} (95% confidence), and an effective UV-efficiency that was at least ~10x greater at z >> 6. These results indicate that the quasar and WMAP observations are consistent. If future observations confirm an optical depth to electron scattering tau_es~0.1, then it would appear that no more "exotic" sources of UV-photons, such as mini-quasars or AGNs, are necessary; but our analysis indicates that a determination of tau_es>~0.17 would require a more radical solution.

It is well known that the mass function for_halos_ in CDM cosmology is a relatively steep power law for low masses, possibly too steep to be consistent with observations. But how steep is the_galaxy_ mass function? We have analyzed the stellar and gas mass functions of the first massive luminous objects formed in a \Lambda CDM universe, as calculated in the numerical simulation described in Gnedin (2000ab). We found that while the dark matter mass function is steep, the stellar and gas mass functions are flatter for low mass objects. The stellar mass function is consistently flat at the low mass end. Moreover, while the gas mass function follows the dark matter mass function until reionization at z~7, between z=7 and z=4, the gas mass function also flattens considerably at the low mass end. At z=4, the gas and stellar mass functions are fit by a Schechter function with \alpha ~ -1.2 +/- 0.1, significantly shallower than the dark matter halo mass function and consistent with some recent observations. The baryonic mass functions are shallower because (a) the dark matter halo mass function is consistent with the Press-Schechter formulation at low masses n(M) M^-2 and (b) heating/cooling and ionization processes appear to cause baryons to collect in halos with the relationship M_b M_d^4 at low masses. Combining (a) and (b) gives n(M_b) M_b^-5/4, comparable to the simulation results. Thus, the well known observational fact that low mass galaxies are underabundant as compared to expectations from numerical dark matter simulations or Press-Schechter modeling of CDM universes emerges naturally from these results, implying that perhaps no ``new physics'' beyond the standard model is needed.

We present a semi-analytic model for the thermal and ionization history of the universe at 1000 >~ z >~ 3. This model incorporates much of the essential physics included in full-scale hydrodynamical simulations, such as (1) gravitational collapse and virialization; (2) star/quasar formation and subsequent ionizing radiation; (3) heating and cooling; (4) atomic and molecular physics of hydrogen; and (5) the feedback relationships between these processes. In addition, we model the process of reheating and reionization using two separate phases, self-consistently calculating the filling factor of each phase. Thus radiative transfer is treated more accurately than simulations to date have done: we allow to lowest order for both the inhomogeneity of the sources and the sinks of radiation. After calibrating and checking the results of this model against a hydrodynamical simulation, we apply our model to a variety of Gaussian and non-Gaussian CDM-dominated cosmological models. Our major conclusions include: (1) the epoch of reheating depends most strongly on the power spectrum of density fluctuations at small scales; (2) because of the effects of gas clumping, full reionization occurs at z ~ 10 in all models; (3) the CMBR polarization and the stars and quasars to baryons ratio are strong potential discrimants between different assumed power spectra; (4) the formation of galactic spheroids may be regulated by the evolution of reheating through feedback, so that the Jeams mass tracks the non-linear mass scale; and (5) the evolution of the bias of luminous objects can potentially discriminate strongly between Gaussian and non-Gaussian PDFs.

We derive an analytic expression for the intensity of resonance-line radiation ``trapped'' in a semi-infinite medium. Given a source function and destruction probability per scattering, the radiation pressure due to trapped photons can be calculated by numerically integrating over analytic functions. We apply this formalism to a plane-parallel model stellar atmosphere to calculate the radiation pressure due to Lyman-alpha photons produced following absorption of UV and X-rays from an AGN. For low surface gravity stars near the AGN (g~10 cm/sec^2, r~0.25 pc), we find that the pressure due to Lyman-alpha photons becomes an appreciable fraction of that required for hydrostatic support. If the broad emission line emitting gas in AGNs and QSOs consists of stellar outflows, it may be driven, in part, by Lyman-alpha pressure.

(Abridged) By comparing the frequency of typical events with that of unusual events, one can test whether the cosmological density distribution function is consistent with the normally made assumption of Gaussianity. To this end, we compare the consistency of the tail-inferred (from clusters) and measured values (from large-scale flows) of the rms level of mass fluctuations for two distribution functions: a Gaussian, and a texture (positively-skewed) PDF. Averaging the recent large-scale flow measurements, we find that observations of the rms and the tail at the 10 h^-1 Mpc scale disfavor a texture PDF at ~1.5 sigma in all cases. However, taking only the most recent measurement of the rms, that from Willick et al. (1997b), the comparison disfavors textures for low Omega_0=0.3, and disfavors Gaussian models if Omega_0=1 (again at ~1.5 sigma). Predictions for evolution of high temperature clusters can also be made for the models considered, and strongly disfavor Omega_0=1 in Gaussian models and marginally disfavor Omega_0=1 in texture models. Only Omega_0=0.3 Gaussian models are consistent with all the data considered.