The evolution of galaxy metallicity scaling relations in cosmological hydrodynamical simulations

The evolution of the metal content of galaxies and its relations to other global properties [such as total stellar mass (M*), circular velocity, star formation rate (SFR), halo mass, etc.] provides important constraints on models of galaxy formation. Here we examine the evolution of metallicity scaling relations of simulated galaxies in the Galaxies-Intergalactic Medium Interaction Calculation suite of cosmological simulations. We make comparisons to observations of the correlation of gas-phase abundances with M* (the mass-metallicity relation, MZR), as well as with both M* and SFR or gas mass fraction (the so-called 3D fundamental metallicity relations, FMRs). The simulated galaxies follow the observed local MZR and FMRs over an order of magnitude in M*, but overpredict the metallicity of massive galaxies (log M* > 10.5), plausibly due to inefficient feedback in this regime. We discuss the origin of the MZR and FMRs in the context of galactic outflows and gas accretion. We examine the evolution of mass-metallicity relations defined using different elements that probe the three enrichment channels (SNII, SNIa, and AGB stars). Relations based on elements produced mainly by SNII evolve weakly, whereas those based on elements produced preferentially in SNIa/AGB exhibit stronger evolution, due to the longer timescales associated with these channels. Finally, we compare the relations of central and satellite galaxies, finding systematically higher metallicities for satellites, as observed. We show this is due to the removal of the metal poor gas reservoir that normally surrounds galaxies and acts to dilute their gas-phase metallicity (via cooling/accretion onto the disk), but is lost due to ram pressure stripping for satellites

The accretion of galaxies into groups and clusters

We use the galaxy stellar mass and halo merger tree information from the semi-analyticmodel galaxy catalogue of Font et al. (2008) to examine the accretion of galaxies into a large sample of groups and clusters, covering a wide range in halo mass (1012.9 to 1015.3 h−1 M⊙), and selected from each of four redshift epochs (z=0, 0.5, 1.0 and 1.5). We find that clusters at all examined redshifts have accreted a significant fraction of their final galaxy populations through galaxy groups. A 1014.5 h−1 M⊙ mass cluster at z=0 has, on average, accreted_ 40% of its galaxies (Mstellar > 109 h−1 M⊙) from halos with masses greater than 1013 h−1 M⊙. Further, the galaxies which are accreted through groups are more massive, on average, than galaxies accreted through smaller halos or from the field population. We find that at a given epoch, the fraction of galaxies accreted from isolated environments is independent of the final cluster or group mass. In contrast, we find that observing a cluster of the same halo mass at each redshift epoch implies different accretion rates of isolated galaxies, from 5-6 % per Gyr at z=0 to 15% per Gyr at z=1.5. We find that combining the existence of a Butcher Oemler effect at z=0.5 and the observations that galaxies within groups display significant environmental effects with galaxy accretion histories justifies striking conclusions. Namely, that the dominant environmental process must begin to occur in halos of 1012 – 1013 h−1 M⊙, and act over timescales of > 2 Gyrs. This argues in favor of a mechanism like “strangulation”, in which the hot halo of a galaxy is stripped upon infalling into a more massive halo . This simple model predicts that by z=1.5 galaxy groups and clusters will display little to no environmental effects. This conclusion may limit the effectiveness of red sequence cluster finding methods at high redshift

Cosmological simulations of the formation of the stellar haloes around disc galaxies

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) suite of cosmological hydrodynamical simulations to study the formation of stellar spheroids of Milky Way mass disc galaxies. The simulations contain accurate treatments of metal-dependent radiative cooling, star formation, supernova feedback and chemodynamics, and the large volumes that have been simulated yield an unprecedentedly large sample of ≈400 simulated ∼L* disc galaxies. The simulated galaxies are surrounded by low-mass, low surface brightness stellar haloes that extend out to ∼100 kpc and beyond. The diffuse stellar distributions bear a remarkable resemblance to those observed around the Milky Way, M31 and other nearby galaxies, in terms of mass density, surface brightness and metallicity profiles. We show that in situ star formation typically dominates the stellar spheroids by mass at radii of r≲ 30 kpc, whereas accretion of stars dominates at larger radii and this change in origin induces a change in the slope of the surface brightness and metallicity profiles, which is also present in the observational data. The system-to-system scatter in the in situ mass fractions of the spheroid, however, is large and spans over a factor of 4. Consequently, there is a large degree of scatter in the shape and normalization of the spheroid density profile within r≲ 30 kpc (e.g. when fitted by a spherical power-law profile, the indices range from −2.6 to −3.4). We show that the in situ mass fraction of the spheroid is linked to the formation epoch of the system. Dynamically, older systems have, on average, larger contributions from in situ star formation, although there is significant system-to-system scatter in this relationship. Thus, in situ star formation likely represents the solution to the long-standing failure of pure accretion-based models to reproduce the observed properties of the inner spheroid

The BAHAMAS project: the CMB--large-scale structure tension and the roles of massive neutrinos and galaxy formation

Recent studies have presented evidence for tension between the constraints on Omega_m and sigma_8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully-calibrated cosmological hydrodynamical simulations, which thus capture the back reaction of the baryons on the total matter distribution. We have extended the BAHAMAS simulations to include a treatment of massive neutrinos, which currently represents the best motivated extension to the standard model. We make synthetic thermal Sunyaev-Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: i) in general there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and iii) invoking a non-minimal neutrino mass, typically of 0.2-0.4 eV (depending on the priors on the other relevant cosmological parameters and the datasets being modelled), can resolve the tension. This solution is fully consistent with separate constraints on the summed neutrino mass from the primary CMB and baryon acoustic oscillations, given the internal tensions in the Planck primary CMB dataset

The Age of the Milky Way Inner Halo

The Milky Way galaxy is observed to have multiple components with distinct properties, such as the bulge, disk, and halo. Unraveling the assembly history of these populations provides a powerful test to the theory of galaxy formation and evolution, but is often restricted due to difficulties in measuring accurate stellar ages for low mass, hydrogen-burning stars. Unlike these progenitors, the "cinders" of stellar evolution, white dwarf stars, are remarkably simple objects and their fundamental properties can be measured with little ambiguity from spectroscopy. Here I report observations and analysis of newly formed white dwarf stars in the halo of the Milky Way, and a comparison to published analysis of white dwarfs in the well-studied 12.5 billion-year-old globular cluster Messier 4. From this, I measure the mass distribution of the remnants and invert the stellar evolution process to develop a new relation that links this final stellar mass to the mass of their immediate progenitors, and therefore to the age of the parent population. By applying this technique to a small sample of four nearby and kinematically-confirmed halo white dwarfs, I measure the age of local field halo stars to be 11.4 +/- 0.7 billion years. This age is directly tied to the globular cluster age scale, on which the oldest clusters formed 13.5 billion years ago. Future (spectroscopic) observations of newly formed white dwarfs in the Milky Way halo can be used to reduce the present uncertainty, and to probe relative differences between the formation time of the last clusters and the inner halo.Comment: Published in Nature, 2012, 486, 90. Second version corrects a missing reference (#10) in the third paragraph and Figure 1 captio

Kinematic Clues to Bar Evolution for Galaxies in the Local Universe: Why the Fastest Rotating Bars are Rotating Most Slowly

We have used Spitzer images of a sample of 68 barred spiral galaxies in the local universe to make systematic measurements of bar length and bar strength. We combine these with precise determinations of the corotation radii associated with the bars, taken from our previous study, which used the phase change from radial inflow to radial outflow of gas at corotation, based on high-resolution two-dimensional velocity fields in Hα taken with a Fabry-Pérot spectrometer. After presenting the histograms of the derived bar parameters, we study their dependence on the galaxy morphological type and on the total stellar mass of the host galaxy, and then produce a set of parametric plots. These include the bar pattern speed versus bar length, the pattern speed normalized with the characteristic pattern speed of the outer disk versus the bar strength, and the normalized pattern speed versus R, the ratio of corotation radius to bar length. To provide guidelines for our interpretation, we used recently published simulations, including disk and dark matter halo components. Our most striking conclusion is that bars with values of R < 1.4, previously considered dynamically fast rotators, can be among the slowest rotators both in absolute terms and when their pattern speeds are normalized. The simulations confirm that this is because as the bars are braked, they can grow longer more quickly than the outward drift of the corotation radius. We conclude that dark matter halos have indeed slowed down the rotation of bars on Gyr timescales. © 2017. The American Astronomical Society. All rights reserved.

Duckietown: An Innovative Way to Teach Autonomy

Teaching robotics is challenging because it is a multidisciplinary, rapidly evolving and experimental discipline that integrates cutting-edge hardware and software. This paper describes the course design and first implementation of Duckietown, a vehicle autonomy class that experiments with teaching innovations in addition to leveraging modern educational theory for improving student learning. We provide a robot to every student, thanks to a minimalist platform design, to maximize active learning; and introduce a role-play aspect to increase team spirit, by modeling the entire class as a fictional start-up (Duckietown Engineering Co.). The course formulation leverages backward design by formalizing intended learning outcomes (ILOs) enabling students to appreciate the challenges of: (a) heterogeneous disciplines converging in the design of a minimal self-driving car, (b) integrating subsystems to create complex system behaviors, and (c) allocating constrained computational resources. Students learn how to assemble, program, test and operate a self-driving car (Duckiebot) in a model urban environment (Duckietown), as well as how to implement and document new features in the system. Traditional course assessment tools are complemented by a full scale demonstration to the general public. The “duckie” theme was chosen to give a gender-neutral, friendly identity to the robots so as to improve student involvement and outreach possibilities. All of the teaching materials and code is released online in the hope that other institutions will adopt the platform and continue to evolve and improve it, so to keep pace with the fast evolution of the field.National Science Foundation (U.S.) (Award IIS #1318392)National Science Foundation (U.S.) (Award #1405259

Linking dwarf galaxies to halo building blocks with the most metal-poor star in Sculptor

Current cosmological models indicate that the Milky Way's stellar halo was assembled from many smaller systems. Based on the apparent absence of the most metal-poor stars in present-day dwarf galaxies, recent studies claimed that the true Galactic building blocks must have been vastly different from the surviving dwarfs. The discovery of an extremely iron-poor star (S1020549) in the Sculptor dwarf galaxy based on a medium-resolution spectrum cast some doubt on this conclusion. However, verification of the iron-deficiency and measurements of additional elements, such as the alpha-element Mg, are mandatory for demonstrating that the same type of stars produced the metals found in dwarf galaxies and the Galactic halo. Only then can dwarf galaxy stars be conclusively linked to early stellar halo assembly. Here we report high-resolution spectroscopic abundances for 11 elements in S1020549, confirming the iron abundance of less than 1/4000th that of the Sun, and showing that the overall abundance pattern mirrors that seen in low-metallicity halo stars, including the alpha-elements. Such chemical similarity indicates that the systems destroyed to form the halo billions of years ago were not fundamentally different from the progenitors of present-day dwarfs, and suggests that the early chemical enrichment of all galaxies may be nearly identical.Comment: 16 pages, including 2 figures. Accepted for publication in Nature. It is embargoed for discussion in the press until formal publication in Natur

Environment from cross-correlations: connecting hot gas and the quenching of galaxies

The observable properties of galaxies are known to depend on both internal processes and the external environment. In terms of the environmental role, we still do not have a clear picture of the processes driving the transformation of galaxies. This may be due to the fact that these environmental processes depend on local physical conditions (e.g., local tidal force or hot gas density), whereas observations typically probe only broad-brush proxies for these conditions (e.g., host halo mass, distance to the N^th nearest neighbour, etc.). Here we propose a new method that directly links galaxies to their local environments, by using spatial cross-correlations of galaxy catalogues with maps from large-scale structure surveys (e.g., thermal Sunyaev-Zel'dovich [tSZ] effect, diffuse X-ray emission, weak lensing of galaxies or the CMB). We focus here on the quenching of galaxies and its link to local hot gas properties. Maps of galaxy overdensity and quenched fraction excess are constructed from volume-limited SDSS catalogs, which are cross-correlated with maps of tSZ effect from Planck and X-ray emission from ROSAT. Strong signals out to Mpc scales are detected for all cross-correlations and are compared to predictions from cosmological hydrodynamical simulations (the EAGLE and BAHAMAS simulations). The simulations successfully reproduce many, but not all, of the observed power spectra, with an indication that environmental quenching may be too efficient in the simulations. We demonstrate that the cross-correlations are sensitive to both the internal and external processes responsible for quenching. The methods outlined in this paper can be easily adapted to other observables and, with upcoming surveys, will provide a stringent direct test of physical models for environmental transformation