The history of stellar metallicity in a simulated disc galaxy

We explore the chemical distribution of stars in a simulated galaxy. Using simulations of the same initial conditions but with two different feedback schemes (McMaster Unbiased Galaxy Simulations – MUGS – and Making Galaxies in a Cosmological Context – MaGICC), we examine the features of the age–metallicity relation (AMR), and the three-dimensional age– [Fe/H]–[O/Fe] distribution, both for the galaxy as a whole and decomposed into disc, bulge, halo and satellites. The MUGS simulation, which uses traditional supernova feedback, is replete with chemical substructure. This substructure is absent from the MaGICC simulation, which includes early feedback from stellar winds, a modified initial mass function and more efficient feedback. The reduced amount of substructure is due to the almost complete lack of satellites in MaGICC. We identify a significant separation between the bulge and disc AMRs, where the bulge is considerably more metal-rich with a smaller spread in metallicity at any given time than the disc. Our results suggest, however, that identifying the substructure in observations will require exquisite age resolution, of the order of 0.25 Gyr. Certain satellites show exotic features in the AMR, even forming a ‘sawtooth’ shape of increasing metallicity followed by sharp declines which correspond to pericentric passages. This fact, along with the large spread in stellar age at a given metallicity, compromises the use of metallicity as an age indicator, although alpha abundance provides a more robust clock at early times. This may also impact algorithms that are used to reconstruct star formation histories from resolved stellar populations, which frequently assume a monotonically increasing AMR

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

Formation of Sub-galactic Clouds under UV Background Radiation

The effects of the UV background radiation on the formation of sub-galactic clouds are studied by means of one-dimensional hydrodynamical simulations. The radiative transfer of the ionizing photons due to the absorption by HI, HeI and HeII, neglecting the emission, is explicitly taken into account. We find that the complete suppression of collapse occurs for the clouds with circular velocities typically in the range V_c \sim 15-40 km/s and the 50% reduction in the cooled gas mass with V_c \sim 20-55 km/s. These values depend most sensitively on the collapse epoch of the cloud, the shape of the UV spectrum, and the evolution of the UV intensity. Compared to the optically thin case, previously investigated by Thoul & Weinberg (1996), the absorption of the external UV photon by the intervening medium systematically lowers the above threshold values by \Delta V_c \sim 5 km/s. Whether the gas can contract or keeps expanding is roughly determined by the balance between the gravitational force and the thermal pressure gradient when it is maximally exposed to the external UV flux. Based on our simulation results, we discuss a number of implications on galaxy formation, cosmic star formation history, and the observations of quasar absorption lines. In Appendix, we derive analytical formulae for the photoionization coefficients and heating rates, which incorporate the frequency/direction-dependent transfer of external photons.Comment: 38 pages, 16 figures, accepted for publication in Ap

Reviews

Reviews of International and comparative industrial relations, Tatau Tatau - one big union altogether, Remedy for present evils: a history of the New Zealand Public Service Association from 1890, Sexual harassment in the workplace, Employee selection, Legislating for workplace hazards in New Zealand: overseas experience and our present and future needs, People and enterprises - human behaviour in New Zealand organisations and From school to unemployment? The labour market for young peopl

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

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

The Role of Heating and Enrichment in Galaxy Formation

We show that the winds identified with high-redshift low-mass galaxies may strongly affect the formation of stars in more massive galaxies that form later. With 3D realizations of a simple linear growth model we track gas shocking, metal enrichment, and cooling, together with dark halo formation. We show that outflows typically strip baryonic material out of collapsing intermediate mass halos, suppressing star formation. More massive halos can trap the heated gas but collapse late, leading to a broad bimodal redshift distribution, with a larger characteristic mass associated with the lower redshift peak. This scenario accounts for the observed bell-shaped luminosity function of early type galaxies, explains the small number of Milky Way satellite galaxies relative to Cold Dark Matter models predictions, and provides a possible explanation for the lack of metal poor G-dwarfs in the solar neighborhood and the more general lack of low-metallicity stars in massive galaxies relative to ``closed-box'' models of chemical enrichment. Intergalactic medium heating from outflows should produce spectral distortions in the cosmic microwave background that will be measurable with the next generation of experiments.Comment: 19 pages, 12 figures, accepted to ApJ, models refined and minor revisions mad

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

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

Search for Yukawa Production of a Light Neutral Higgs Boson at LEP

Within a Two-Higgs-Doublet Model (2HDM) a search for a light Higgs boson in the mass range of 4-12 GeV has been performed in the Yukawa process e+e- -> b bbar A/h -> b bbar tau+tau-, using the data collected by the OPAL detector at LEP between 1992 and 1995 in e+e- collisions at about 91 GeV centre-of-mass energy. A likelihood selection is applied to separate background and signal. The number of observed events is in good agreement with the expected background. Within a CP-conserving 2HDM type II model the cross-section for Yukawa production depends on xiAd = |tan beta| and xihd = |sin alpha/cos beta| for the production of the CP-odd A and the CP-even h, respectively, where tan beta is the ratio of the vacuum expectation values of the Higgs doublets and alpha is the mixing angle between the neutral CP-even Higgs bosons. From our data 95% C.L. upper limits are derived for xiAd within the range of 8.5 to 13.6 and for xihd between 8.2 to 13.7, depending on the mass of the Higgs boson, assuming a branching fraction into tau+tau- of 100%. An interpretation of the limits within a 2HDM type II model with Standard Model particle content is given. These results impose constraints on several models that have been proposed to explain the recent BNL measurement of the muon anomalous magnetic moment.Comment: 24 pages, 9 figures, Submitted to Euro. Phys. J.