The formation of the first stars is an exciting frontier area in astronomy. Early redshifts z ~ 20 have become observationally promising as a result of a recently recognized effect of a supersonic relative velocity between the dark matter and gas. This effect produces prominent structure on 100 comoving Mpc scales, which makes it much more feasible to detect 21-cm fluctuations from the epoch of first heating. We use semi-numerical hybrid methods to follow for the first time the joint evolution of the X-ray and Lyman-Werner radiative backgrounds, including the effect of the supersonic streaming velocity on the cosmic distribution of stars. We incorporate self-consistently the negative feedback on star formation induced by the Lyman-Werner radiation, which dissociates molecular hydrogen and thus suppresses gas cooling. We find that the feedback delays the X-ray heating transition by a Delta z ~ 2, but leaves a promisingly large fluctuation signal over a broad redshift range. The large-scale power spectrum is predicted to reach a maximal signal-to-noise ratio of S/N ~ 3-4 at z ~ 18 (for a projected first-generation instrument), with S/N > 1 out to z ~ 22-23. We hope to stimulate additional numerical simulations as well as observational efforts focused on the epoch prior to cosmic reionization.

Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.

In astrophysical regimes where the collisional excitation of hydrogen atoms is relevant, the cross sections for the interactions of hydrogen atoms with electrons and protons are necessary for calculating line profiles and intensities. In particular, at relative velocities exceeding ~1000 km/s, collisional excitation by protons dominates over that by electrons. Surprisingly, the hydrogen-proton cross sections at these velocities do not exist for atomic levels of n >= 4, forcing researchers to utilize extrapolation via inaccurate scaling laws. In this study, we present a faster and improved algorithm for computing cross sections for the hydrogen-proton collisional system, including excitation and charge transfer to the n >= 2 levels of the hydrogen atom. We develop a code named BDSCx which directly solves the Schrodinger equation with variable (but non-adaptive) resolution and utilizes a hybrid spatial-Fourier grid. Our novel hybrid grid reduces the number of grid points needed from ~4000 n^6 (for a "brute force", Cartesian grid) to ~2000 n^4 and speeds up the computation by a factor ~50 for calculations going up to n = 4 . We present (l,m)-resolved results for charge-transfer and excitation final states for n = 2--4 and for projectile energies of 5--80 keV, as well as fitting functions for the cross sections. The ability to accurately compute proton-hydrogen cross sections to n = 4 allows us to calculate the Balmer decrement, the ratio of Balmer alpha to Balmer beta line intensities. We find that the Balmer decrement starts to increase beyond its largely constant value of 2--3 below 10 keV, reaching values of 4--5 at 5 keV, thus complicating its use as a diagnostic of dust extinction when fast (~1000$ km/s) shocks are impinging upon the ambient interstellar medium.

Recently the initial supersonic relative velocity between the dark matter and baryons was shown to have an important effect on galaxy formation at high redshift. We study the impact of this relative motion on the distribution of the star-forming halos and on the formation redshift of the very first star. We include a new aspect of the relative velocity effect found in recent simulations by fitting their results to obtain the spatially-varying minimum halo mass needed for molecular cooling. Thus, the relative velocities have three separate effects: suppression of the halo abundance, suppression of the gas content within each halo, and boosting of the minimum cooling mass. We show that the two suppressions (of gas content and of halo abundance) are the primary effects on the small minihalos that cannot form stars, while the cooling mass boost combines with the abundance suppression to produce order unity fluctuations in stellar density. We quantify the large-scale inhomogeneity of galaxies, finding that 68% of the star formation (averaged on a 3 Mpc scale) is confined to 35% of the volume at z=20 (and just 18% at z=40). In addition, we estimate the redshift of the first star to be z ~ 65, which includes a delay of Dz ~ 5 due to the relative velocity.

We study the effect of the relative velocity of dark matter and baryonic fluids after the epoch of recombination on the evolution of the first bound objects in the early universe. Recent work has shown that, although relative motion of the two fluids is formally a second order effect in density, it has a dramatic impact on the formation and distribution of the first cosmic structures. Focusing on the gas content, we analyze the effect of relative velocity on the properties of halos over a wide range of halo masses and redshifts. We calculate accurately the linear evolution of the baryon and dark matter fluctuations, and quantify the resulting effect on halos based on an analytical formalism that has been carefully checked with simulations in the case with no relative velocity. We estimate the effect on the abundance of and gas fraction in early halos. We find that the relative velocity effect causes several changes: (i) the characteristic mass that divides gas-rich and gas-poor halos is increased by roughly an order of magnitude, from 2 10^4 Msun to about 2 10^5 Msun; (ii) this characteristic mass has a large scatter (full width at half maximum is ~ 1.5 10^5 Msun at z=20); (iii) the fraction of baryons in star-less gas minihalos is suppressed by a factor of 4 at z=20; (iv) the fraction of baryons in halos that can cool and form stars is suppressed by a factor of 1.5 at z=20; and (v) there are enhanced spatial variations of these various fractions.

Single field inflationary models predict nearly Gaussian initial conditions and hence a detection of non-Gaussianity would be a signature of the more complex inflationary scenarios. In this paper we study the effect on the cosmic microwave background and on large scale structure from primordial non-Gaussianity in a two-field inflationary model in which both the inflaton and curvaton contribute to the density perturbations. We show that in addition to the previously described enhancement of the galaxy bias on large scales, this setup results in large-scale stochasticity. We provide joint constraints on the local non-Gaussianity parameter $\tilde f_{\rm NL}$ and the ratio $\xi$ of the amplitude of primordial perturbations due to the inflaton and curvaton using WMAP and SDSS data.

At the time of recombination, baryons and photons decoupled and the sound speed in the baryonic fluid dropped from relativistic to the thermal velocities of the hydrogen atoms. This is less than the relative velocities of baryons and dark matter computed via linear perturbation theory, so we infer that there are supersonic coherent flows of the baryons relative to the underlying potential wells created by the dark matter. As a result, the advection of small-scale perturbations (near the baryonic Jeans scale) by large-scale velocity flows is important for the formation of the first baryonic structures. This effect involves a quadratic term in the cosmological perturbation theory equations and hence has not been included in studies based on linear perturbation theory. We show that the relative motion suppresses the abundance of the first bound objects, even if one only investigates dark matter haloes, and leads to qualitative changes in their spatial distribution, such as introducing scale-dependent bias and stochasticity. We discuss the possible observable implications for high-redshift galaxy clustering and reionization.