Cameron B. Hummels

ENZO: AN ADAPTIVE MESH REFINEMENT CODE FOR ASTROPHYSICS

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in one, two, and three dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the codes parallel performance, and discuss the Enzo collaborations code development methodology.

The GALEX Arecibo SDSS Survey - I. Gas fraction scaling relations of massive galaxies and first data release

We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large programme that is gathering high quality H i-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^(10) M_⊙ and redshifts 0.025 < z < 0.05 , selected from the Sloan Digital Sky Survey (SDSS) spectroscopic and Galaxy Evolution Explorer (GALEX) imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5–5 per cent) is reached. This paper presents the first Data Release, consisting of ~20 per cent of the final GASS sample. We use this data set to explore the main scaling relations of the H i gas fraction with galaxy structure and NUV−r colour. A large fraction (~60 per cent) of the galaxies in our sample are detected in H i. Even at stellar masses above 10^(11) M_⊙, the detected fraction does not fall below ~40 per cent. We find that the atomic gas fraction M_(HI)/M★ decreases strongly with stellar mass, stellar surface mass density and NUV−r colour, but is only weakly correlated with the galaxy bulge-to-disc ratio (as measured by the concentration index of the r-band light). We also find that the fraction of galaxies with significant (more than a few per cent) H I decreases sharply above a characteristic stellar surface mass density of 10^(8.5) M_⊙ kpc^(−2). The fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of the key goals of GASS is to identify and quantify the incidence of galaxies that are transitioning between the blue, star-forming cloud and the red sequence of passively evolving galaxies. Likely transition candidates can be identified as outliers from the mean scaling relations between M_(HI)/M★ and other galaxy properties. We have fitted a plane to the two-dimensional relation between the H I mass fraction, stellar surface mass density and NUV−r colour. Interesting outliers from this plane include gas-rich red sequence galaxies that may be in the process of regrowing their discs, as well as blue, but gas-poor spirals.

The GALEX Arecibo SDSS Survey – VIII. Final data release. The effect of group environment on the gas content of massive galaxies

We present the final data release from the GALEX Arecibo SDSS Survey (GASS), a large Arecibo program that measured the Hi properties for an unbiased sample of �800 galaxies with stellar masses greater than 10 10 M⊙ and redshifts 0.025 < z < 0.05. This release includes new Arecibo observations for 250 galaxies. We use the full GASS sample to investigate environmental effects on the cold gas content of massive galaxies at fixed stellar mass. The environment is characterized in terms of dark matter halo mass, obtained by cross-matching our sample with the SDSS group catalog of Yang et al. Our analysis provides, for the first time, clear statistical evidence that massive galaxies located in halos with masses of 10 13 10 14 M⊙ have at least 0.4 dex less Hi than objects in lower density environments. The process responsible for the suppression of gas in group galaxies most likely drives the observed quenching of the star formation in these systems. Our findings strongly support the importance of the group environment for galaxy evolution, and have profound implications for semi-analytic models of galaxy formation, which currently do not allow for stripping of the cold interstellar medium in galaxy groups.

The GALEX Arecibo SDSS Survey II: The Star Formation Efficiency of Massive Galaxies

We use measurements of the H I content, stellar mass and star formation rates (SFRs) in ~190 massive galaxies with M_★ > 10^(10) M_⊙, obtained from the GALEX (Galaxy Evolution Explorer) Arecibo SDSS (Sloan Digital Sky Survey) survey described in Paper I to explore the global scaling relations associated with the bin-averaged ratio of the SFR over the H I mass (i.e. ΣSFR/ΣM_(HI)), which we call the H I-based star formation efficiency (SFE). Unlike the mean specific star formation rate (sSFR), which decreases with stellar mass and stellar mass surface density, the SFE remains relatively constant across the sample with a value close to SFE = 10^(−9.5) yr^(−1) (or an equivalent gas consumption time-scale of ~3 × 10^9 yr). Specifically, we find little variation in SFE with stellar mass, stellar mass surface density, NUV −r colour and concentration (R_(90)/R_(50)). We interpret these results as an indication that external processes or feedback mechanisms that control the gas supply are important for regulating star formation in massive galaxies. An investigation into the detailed distribution of SFEs reveals that approximately 5 per cent of the sample shows high efficiencies with SFE > 10^(−9) yr^(−1), and we suggest that this is very likely due to a deficiency of cold gas rather than an excess SFR. Conversely, we also find a similar fraction of galaxies that appear to be gas-rich for their given sSFR, although these galaxies show both a higher than average gas fraction and lower than average sSFR. Both of these populations are plausible candidates for ‘transition’ galaxies, showing potential for a change (either decrease or increase) in their sSFR in the near future. We also find that 36 ± 5 per cent of the total H I mass density and 47 ± 5 per cent of the total SFR density are found in galaxies with M_★ > 10^(10) M_⊙.

Constraints on hydrodynamical subgrid models from quasar absorption line studies of the simulated circumgalactic medium

Cosmological hydrodynamical simulations of galaxy evolution are increasingly able to produce realistic galaxies, but the largest hurdle remaining is in constructing subgrid models that accurately describe the behavior of stellar feedback. As an alternate way to test and calibrate such models, we propose to focus on the circumgalactic medium. To do so, we generate a suite of adaptive-mesh refinement (AMR) simulations for a Milky-Way-massed galaxy run to z=0, systematically varying the feedback implementation. We then post-process the simulation data to compute the absorbing column density for a wide range of common atomic absorbers throughout the galactic halo, including H I, Mg II, Si II, Si III, Si IV, C IV, N V, O VI, and O VII. The radial profiles of these atomic column densities are compared against several quasar absorption line studies, to determine if one feedback prescription is favored. We find that although our models match some of the observations (specifically those ions with lower ionization strengths), it is particularly difficult to match O VI observations. There is some indication that the models with increased feedback intensity are better matches. We demonstrate that sufficient metals exist in these halos to reproduce the observed column density distribution in principle, but the simulated circumgalactic medium lacks significant multiphase substructure and is generally too hot. Furthermore, we demonstrate the failings of inflow-only models (without energetic feedback) at populating the CGM with adequate metals to match observations even in the presence of multiphase structure. Additionally, we briefly investigate the evolution of the CGM from z=3 to present. Overall, we find that quasar absorption line observations of the gas around galaxies provide a new and important constraint on feedback models.

FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (“FIRE-1”) for consistency. Motivated by the development of more accurate numerics – including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms – and exploration of new physics (e.g. magnetic fields), we introduce “FIRE-2”, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (∼kpc) mass concentrations in massive (>L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

The GALEX Arecibo SDSS Survey - VI. Second data release and updated gas fraction scaling relations

We present the second data release from the GALEX Arecibo SDSS Survey (GASS), an ongoing large Arecibo program to measure the Hi properties for an unbiased sample of∼1000 galaxies with stellar masses greater than 10 10 M⊙ and redshifts 0.025< z< 0.05. GASS targets are selected from the Sloan Digital Sky Survey (SDSS) spectroscopic and Galaxy Evolution Explorer (GALEX)imaging surveys, and are observed until detected or until a gas mass fraction limit of a few per cent is reached. This second data installment includes new Arecibo observations of 240 galaxies, and marks the 50% of the complete survey. We present catalogs of the Hi, optical and ultraviolet parameters for these galaxies, and their Hi-line profiles. Having more than doubled the size of the sampl e since the first data release, we also revisit the main scaling relations of t he Hi mass fraction with galaxy stellar mass, stellar mass surfac e density, concentration index, and NUV−r color, as well as the gas fraction plane introduced in our ear lier work.

Adaptive Mesh Refinement Simulations of Galaxy Formation: Exploring Numerical and Physical Parameters

We carry out adaptive mesh refinement cosmological simulations of Milky Way mass halos in order to investigate the formation of disk-like galaxies in a Λ-dominated cold dark matter model. We evolve a suite of five halos to z = 0 and find a gas disk formation in each; however, in agreement with previous smoothed particle hydrodynamics simulations (that did not include a subgrid feedback model), the rotation curves of all halos are centrally peaked due to a massive spheroidal component. Our standard model includes radiative cooling and star formation, but no feedback. We further investigate this angular momentum problem by systematically modifying various simulation parameters including: (1) spatial resolution, ranging from 1700 to 212 pc; (2) an additional pressure component to ensure that the Jeans length is always resolved; (3) low star formation efficiency, going down to 0.1%; (4) fixed physical resolution as opposed to comoving resolution; (5) a supernova feedback model that injects thermal energy to the local cell; and (6) a subgrid feedback model which suppresses cooling in the immediate vicinity of a star formation event. Of all of these, we find that only the last (cooling suppression) has any impact on the massive spheroidal component. In particular, a simulation with cooling suppression and feedback results in a rotation curve that, while still peaked, is considerably reduced from our standard runs.

The GALEX Arecibo SDSS Survey IV: baryonic mass-velocity-size relations of massive galaxies

We present dynamical scaling relations for a homogeneous and representative sample of �500 massive galaxies, selected only by stellar mass (> 10 10 M⊙) and redshift (0.025 < z < 0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate how galaxies scatter around the best fits obtained for pruned subsets of disk-dominated and bulge-dominated systems. The BFJ relation is significantly less scattered than the BTF when the relations are applied to their maximum samples (for the BTF, only galaxies with Hi detections), and is not affected by the inclination problems that plague the BTF. Disk-dominated, gas-rich galaxies systematically deviate from the BFJ relation defined by the spheroids. We demonstrate that by applying a simple correction to the stellar velocity dispersions that depends only on the concentration index of the galaxy, we are able to bring disks and spheroids onto the same dynamical relation — in other words, we obtain a generalized BFJ relation that holds for all the galaxies in our sample, regardless of morphology, inclination or gas content, and has a scatter smaller than 0.1 dex. We compare the velocity-size relation for the three dynamical indicators used in this work, i.e., rotational velocity, observed and concentration-corrected stellar dispersion. We find that disks and spheroids are offset in the stellar dispersion-size relation, and that the offset is removed when corrected dispersions are used instead. The generalized BFJ relation represents a fundamental correlation between the global dark matter and baryonic content of galaxies, which is obeyed by all (massive) systems regardless of morphology.

COSMOLOGICAL SIMULATIONS OF GALAXY FORMATION WITH COSMIC RAYS

We investigate the dynamical impact of cosmic rays in cosmological simulations of galaxy formation using adaptive-mesh refinement simulations of a