The stellar metallicity distribution of disc galaxies and bulges in cosmological simulations

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.Comment: 20 pages, 14 figures, MNRAS, accepte

Combinatorial identities for binary necklaces from exact ray-splitting trace formulae

Based on an exact trace formula for a one-dimensional ray-splitting system, we derive novel combinatorial identities for cyclic binary sequences (P\'olya necklaces).Comment: 15 page

The Bispectrum as a Signature of Gravitational Instability in Redshift-Space

The bispectrum provides a characteristic signature of gravitational instability that can be used to probe the Gaussianity of the initial conditions and the bias of the galaxy distribution. We study how this signature is affected by redshift distortions using perturbation theory and high-resolution numerical simulations. We obtain perturbative results for the multipole expansion of the redshift-space bispectrum which provide a natural way to break the degeneracy between bias and $\Omega$ present in measurements of the redshift-space power spectrum. We propose a phenomenological model that incorporates the perturbative results and also describes the bispectrum in the transition to the non-linear regime. We stress the importance of non-linear effects and show that inaccurate treatment of these can lead to significant discrepancies in the determination of bias from galaxy redshift surveys. At small scales we find that the bispectrum monopole exhibits a strong configuration dependence that reflects the velocity dispersion of clusters. Therefore, the hierarchical model for the three-point function does not hold in redshift-space.Comment: 19 pages, 4 figures. Revised version accepted for publication in Ap

First Structure Formation: I. Primordial Star Forming Regions in hierarchical models

We investigate the possibility of very early formation of primordial star clusters from high-\sigma perturbations in cold dark matter dominated structure formation scenarios. For this we have developed a powerful 2-level hierarchical cosmological code with a realistic and robust treatment of multi-species primordial gas chemistry, paying special attention to the formation and destruction of hydrogen molecules, non-equilibrium ionization, and cooling processes. We performed 3-D simulations at small scales and at high redshifts and find that, analogous to simulations of large scale structure, a complex system of filaments, sheets, and spherical knots at the intersections of filaments form. On the mass scales covered by our simulations (5x10^5 - 1x10^9\Ms) that collapse at redshifts z>25, we find that only at the spherical knots can enough H2 be formed (n_{H_2}/n_H > 5x10^-4) to cool the gas appreciably. Quantities such as the time dependence of the formation of H2 molecules, the final H2 fraction, and central densities from the simulations are compared to the theoretical predictions of Abel (1995) and Tegmark et al. (1997) and found to agree remarkably well. Comparing the 3-D results to an isobaric collapse model we further discuss the possible implications of the extensive merging of small structure that is inherent in hierarchical models. Typically only 5-8% percent of the total baryonic mass in the collapsing structures is found to cool significanlty. Assuming the Padoan (1995) model for star formation our results would predict the first stellar systems to be as small as ~30\Ms. Some implications for primordial globular cluster formation scenarios are also discussed.Comment: 22 pages, 13 Figures. Submitted to ApJ. Laboratory for Computational Astrophysics at the National Center for Supercomputing Application

Modeling scale-dependent bias on the baryonic acoustic scale with the statistics of peaks of Gaussian random fields

Models of galaxy and halo clustering commonly assume that the tracers can be treated as a continuous field locally biased with respect to the underlying mass distribution. In the peak model pioneered by BBKS, one considers instead density maxima of the initial, Gaussian mass density field as an approximation to the formation site of virialized objects. In this paper, the peak model is extended in two ways to improve its predictive accuracy. Firstly, we derive the two-point correlation function of initial density peaks up to second order and demonstrate that a peak-background split approach can be applied to obtain the k-independent and k-dependent peak bias factors at all orders. Secondly, we explore the gravitational evolution of the peak correlation function within the Zel'dovich approximation. We show that the local (Lagrangian) bias approach emerges as a special case of the peak model, in which all bias parameters are scale-independent and there is no statistical velocity bias. We apply our formulae to study how the Lagrangian peak biasing, the diffusion due to large scale flows and the mode-coupling due to nonlocal interactions affect the scale dependence of bias from small separations up to the baryon acoustic oscillation (BAO) scale. For 2-sigma density peaks collapsing at z=0.3, our model predicts a ~ 5% residual scale-dependent bias around the acoustic scale that arises mostly from first-order Lagrangian peak biasing (as opposed to second-order gravity mode-coupling). We also search for a scale dependence of bias in the large scale auto-correlation of massive halos extracted from a very large N-body simulation provided by the MICE collaboration. For halos with mass M>10^{14}Msun/h, our measurements demonstrate a scale-dependent bias across the BAO feature which is very well reproduced by a prediction based on the peak model.Comment: (v1): 23 pages text, 8 figures + appendix (v2): typos fixed, references added, accepted for publication in PR

The Universe Was Reionized Twice

We show the universe was reionized twice, first at z~15-16 and second at z~6. Such an outcome appears inevitable, when normalizing to two well determined observational measurements, namely, the epoch of the final cosmological reionization at z~6 and the density fluctuations at z~6, which in turn are tight ly constrained by Lyman alpha forest observations at z~3. These two observations most importantly fix the product of star formation efficiency and ionizing photon escape fraction from galaxies at high redshift. To the extent that the relative star formation efficiencies in gaseous minihalos with H2 cooling and large halos with atomic cooling at high redshift are still unknown, the primary source for the first reionization could be Pop III stars either in minihalos or in large halos. We show that gas in minihalos can be cooled efficiently by H2 molecules and star formation can continue to take place largely unimpeded throughout the first reionization period, thanks to two new mechanisms for generating a high X-ray background during the Pop III era, put forth here. Moreover, an important process for producing a large number of H2 molecules in relic HII regions of Pop III galaxies, first pointed out by Ricotti, Gnedin, & Shull, is quantified here. It is shown that the Lyman-Werner background may never build up during the Pop III era. The long cosmological reionization and reheating history is complex. We discuss a wide range of implications and possible tests for this new reionization picture. In particular, Thomson scattering optical depth is increased to 0.10 +- 0.03, compared to 0.027 for the case of only one rapid reionization at z=6. Upcoming Microwave Anisotropy Probe observation of the polarization of the cosmic microwave background should be able to distinguish between these two scenarios.Comment: submitted to ApJ, 69 pages, substantial revision made and conclusions strengthene

The mass function of the Las Campanas loose groups of galaxies

We have determined the mass function of loose groups of galaxies in the Las Campanas Redshift Survey. Loose groups of galaxies in the LCRS range in mass from M \sim 10^{12} {\rm M}_{\sun} to 10^{15} {\rm M}_{\sun}. We find that the sample is almost complete for masses in the interval 5\cdot 10^{13}-8\cdot 10^{14} {\rm M}_{\sun}. Comparison of the observed mass function with theoretical mass functions obtained from N-body simulations shows good agreement with a CDM model with the parameters $\Omega_m = 0.3$, $\Omega_{\Lambda} = 0.7$ and the amplitude of perturbations about $\sigma_8=0.78-0.87$. For smaller masses the mass function of LCRS loose groups flattens out, differing considerably from the group mass function found by Girardi and Giuricin (2000) and from mass functions obtained by numerical simulations.Comment: 9 pages, 7 figures, AA accepte

Observational Constraints on Open Inflation Models

We discuss observational constraints on models of open inflation. Current data from large-scale structure and the cosmic microwave background prefer models with blue spectra and/or Omega_0 >= 0.3--0.5. Models with minimal anisotropy at large angles are strongly preferred.Comment: 4 pages, RevTeX, with 2 postscript figures included. Second Figure correcte

Natural Inflation: Particle Physics Models, Power Law Spectra for Large Scale Structure, and Constraints from COBE

A pseudo-Nambu-Goldstone boson, with a potential of the form $V(\phi) = \Lambda^4[1 \pm \cos(\phi/f)], naturally gives rise to inflation if$f \sim M_{Pl}$and$\Lambda \sim M_{GUT}$. We show how this can arise in technicolor-like and superstring models, and work out an explicit string example in the context of multiple gaugino condensation models. We study the cosmology of this model in detail, and find that sufficient reheating to ensure that baryogenesis can take place requires$f > 0.3 M_{Pl}$. The primordial density fluctuation spectrum generated is a non-scale-invariant power law,$P(k) \propto k^{n_s}$, with$n_s \simeq 1 - (M^2_{Pl}/8\pi f^2)$, leading to more power on large length scales than the$n_s = 1$Harrison-Zeldovich spectrum. The standard CDM model with$0 \la n_s \la 0.6-0.7$could in principle explain the large-scale clustering observed in the APM and IRAS galaxy surveys as well as large-scale flows, but the COBE microwave anisotropy implies such low amplitudes (or high bias factors,$b>2$) for these CDM models that galaxy formation occurs too late to be viable; combining COBE with sufficiently early galaxy formation or the large-scale flows leads to$n_s >0.6$, or$f > 0.3 M_{Pl}$as well. For extended and power law inflation models, this constraint is even tighter,$n_s > 0.7$; combined with other bounds on large bubbles in extended inflation, this leaves little room for most extended models.Comment: 42 pages, (12 figures not included but available from the authors

Simulation techniques for cosmological simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.Comment: 42 pages, 16 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 12; work done by an international team at the International Space Science Institute (ISSI), Bern, organised by J.S. Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke