Source author record

Romeel Dave

Romeel Dave appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph astro-ph.IM Computer Vision Machine Learning

Catalog footprint

What is connected

52works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments.

preprint2022arXiv

Co-evolution of massive black holes and their host galaxies at high redshift: discrepancies from six cosmological simulations and the key role of JWST

The James Webb Space Telescope will have the power to characterize high-redshift quasars at z>6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e., with bolometric luminosity L_bol> 10^46 erg/s) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with L_bol= 10^45-10^46 erg/s opens a path to understand the build-up of more normal BHs at z>6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z>4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the M_BH-M_star scaling relation at z=0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z>4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with L_bol= 10^45-10^46 erg/s observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z=6, with the same diversity of sub-grid physics, in order to gain statistics on the most extreme objects at high redshift.

preprint2022arXiv

Effects of Active Galactic Nucleus Feedback on Cold Gas Depletion and Quenching of Central Galaxies

We investigate the influence of active galactic nucleus (AGN) feedback on the galaxy cold gas content and its connection to galaxy quenching in three hydrodynamical simulations of Illustris, IllustrisTNG and SIMBA. By comparing to the observed atomic and molecular neutral hydrogen measurements for central galaxies, we find that Illustris over-predicts the cold gas masses in star-forming galaxies and significantly under-predicts them for quenched galaxies. IllustrisTNG performs better in this comparison than Illustris, but quenched galaxies retain too much cold gas compared with observations. SIMBA shows good agreement with observations, by depleting the global cold gas reservoir for quenched galaxies. We find that the discrepancies in IllustrisTNG are caused by its weak kinetic AGN feedback that only redistributes the cold gas from the inner disks to the outer regions and reduces the inner cold gas densities. It agrees with observations much better when only the cold gas within the stellar disk is considered to infer the star formation rates. From dependences of cold gas reservoir on the black hole mass and Eddington ratio, we find that the cumulative energy release during the black hole growth is the dominant reason for the cold gas depletion and thus the galaxy quenching. We further measure the central stellar surface density within 1 kpc ($Σ_1$) for the high-resolution run of IllustrisTNG and find a tight correlation between $Σ_1$ and black hole mass. It suggests that the observed decreasing trend of cold gas mass with $Σ_1$ is also a reflection of the black hole growth.

preprint2022arXiv

How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-Dust Geometry Model in SED Fitting

One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the Simba cosmological simulation with the Prospector SED fitting code. We consider the impact of the commonly-assumed uniform screen model and introduce a new non-uniform screen model parameterized by the fraction of unobscured stellar light. This non-uniform screen model allows for a non-zero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a non-uniform screen model, with the median offset in UV attenuation decreasing from $-0.30$ dex with a uniform screen model to $-0.17$ dex with the non-uniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the non-uniform screen model, decreasing the SFR offset on average by $0.12$ dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.

preprint2022arXiv

Rapidly quenched galaxies in the Simba cosmological simulation and observations

A wide range of mechanisms have been put forward to explain the quenching of star formation in galaxies with cosmic time, however, the true balance of responsible mechanisms remains unknown. The identification and study of galaxies that have shut down their star formation on different timescales might elucidate which mechanisms dominate at different epochs and masses. Here we study the population of rapidly quenched galaxies (RQGs) in the SIMBA cosmological hydrodynamic simulation at $0.5<z<2$, comparing directly to observational post-starburst galaxies in the UKIDSS Ultra Deep Survey via their colour distributions and mass functions. We find that the fraction of quiescent galaxies that are rapidly quenched in SIMBA is 59% (or 48% in terms of stellar mass), which is higher than observed. A similar "downsizing" of RQGs is observed in both SIMBA and the UDS, with RQGs at higher redshift having a higher average mass. However, SIMBA produces too many RQGs at $1<z_q<1.5$ and too few low mass RQGs at $0.5<z_q<1$. The precise colour distribution of SIMBA galaxies compared to the observations also indicates various inconsistencies in star formation and chemical enrichment histories, including an absence of short, intense starbursts. Our results will help inform the next generation of galaxy evolution models, particularly with respect to the quenching mechanisms employed.

preprint2022arXiv

The black hole population in low-mass galaxies in large-scale cosmological simulations

Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that are too massive, and that the BH occupation fraction at z=0 is not inherited from the simulation seeding modeling. The ability of BHs and their host galaxies to power an AGN depends on BH and galaxy subgrid modeling. The fraction of AGN in low-mass galaxies is not used to calibrate the simulations, and thus can be used to differentiate galaxy formation models. AGN fractions at z=0 span two orders of magnitude at fixed galaxy stellar mass in simulations, similarly to observational constraints, but uncertainties and degeneracies affect both observations and simulations. The agreement is difficult to interpret due to differences in the masses of simulated and observed BHs, BH occupation fraction affected by numerical choices, and an unknown fraction of obscured AGN. Our work advocates for more thorough comparisons with observations to improve the modeling of cosmological simulations, and our understanding of BH and galaxy physics in the low-mass regime. The mass of BHs, their ability to efficiently accrete gas, and the AGN fraction in low-mass galaxies have important implications for the build-up of the entire BH and galaxy populations with time.

preprint2022arXiv

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope.

preprint2022arXiv

The Dwarf Galaxy Population at $z\sim 0.7$: A Catalog of Emission Lines and Redshifts from Deep Keck Observations

We present a catalog of spectroscopically measured redshifts over $0 < z < 2$ and emission line fluxes for 1440 galaxies. The majority ($\sim$65\%) of the galaxies come from the HALO7D survey, with the remainder from the DEEPwinds program. This catalog includes redshifts for 646 dwarf galaxies with $\log(M_{\star}/M_{\odot}) < 9.5$. 810 catalog galaxies did not have previously published spectroscopic redshifts, including 454 dwarf galaxies. HALO7D used the DEIMOS spectrograph on the Keck II telescope to take very deep (up to 32 hours exposure, with a median of $\sim$7 hours) optical spectroscopy in the COSMOS, EGS, GOODS-North, and GOODS-South CANDELS fields, and in some areas outside CANDELS. We compare our redshift results to existing spectroscopic and photometric redshifts in these fields, finding only a 1\% rate of discrepancy with other spectroscopic redshifts. We measure a small increase in median photometric redshift error (from 1.0\% to 1.3\%) and catastrophic outlier rate (from 3.5\% to 8\%) with decreasing stellar mass. We obtained successful redshift fits for 75\% of massive galaxies, and demonstrate a similar 70-75\% successful redshift measurement rate in $8.5 < \log(M_{\star}/M_{\odot}) < 9.5$ galaxies, suggesting similar survey sensitivity in this low-mass range. We describe the redshift, mass, and color-magnitude distributions of the catalog galaxies, finding HALO7D galaxies representative of CANDELS galaxies up to \textit{i}-band magnitudes of 25. The catalogs presented will enable studies of star formation (SF), the mass-metallicity relation, SF-morphology relations, and other properties of the $z\sim0.7$ dwarf galaxy population.

preprint2022arXiv

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.

preprint2021arXiv

Detection of a Multiphase Intragroup Medium: Results from the COS-IGrM Survey

We present the results of the COS Intragroup Medium (COS-IGrM) Survey that used the Cosmic Origins Spectrograph on the Hubble Space Telescope to observe a sample of 18 UV bright quasars, each probing the intragroup medium (IGrM) of a galaxy group. We detect Ly$α$, C II, N V, Si II, Si III, and O VI in multiple sightlines. The highest ionization species detected in our data is O VI, which was detected in 8 out of 18 quasar sightlines. The wide range of ionization states observed provide evidence that the IGrM is patchy and multiphase. We find that the O VI detections generally align with radiatively cooling gas between $10^{5.8}$ and $10^6$ K. The lack of O VI detections in 10 of the 18 groups illustrates that O VI may not be the ideal tracer of the volume filling component of the IGrM. Instead, it either exists at trace levels in a hot IGrM or is generated in the boundary between the hotter IGrM and cooler gas.

preprint2021arXiv

Finding universal relations in subhalo properties with artificial intelligence

We use a generic formalism designed to search for relations in high-dimensional spaces to determine if the total mass of a subhalo can be predicted from other internal properties such as velocity dispersion, radius, or star-formation rate. We train neural networks using data from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and show that the model can predict the total mass of a subhalo with high accuracy: more than 99% of the subhalos have a predicted mass within 0.2 dex of their true value. The networks exhibit surprising extrapolation properties, being able to accurately predict the total mass of any type of subhalo containing any kind of galaxy at any redshift from simulations with different cosmologies, astrophysics models, subgrid physics, volumes, and resolutions, indicating that the network may have found a universal relation. We then use different methods to find equations that approximate the relation found by the networks and derive new analytic expressions that predict the total mass of a subhalo from its radius, velocity dispersion, and maximum circular velocity. We show that in some regimes, the analytic expressions are more accurate than the neural networks. We interpret the relation found by the neural network and approximated by the analytic equation as being connected to the virial theorem.

preprint2021arXiv

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.

preprint2021arXiv

The Three Hundred: Cluster Dynamical States and Relaxation Time Scale

Using the galaxy clusters from The Three Hundred Project, we define a new parameter: $λ_{DS}$ to describe the dynamical state of clusters, which assumes a double-Gaussian distribution in logarithm scale for our mass-complete cluster sample at $z=0$ from the dark-matter-only (DMO) run. Therefore, the threshold for distinguishing relaxed and unrelaxed clusters is naturally determined by the crossing point of the double-Gaussian fitting which has a value of $λ_{DS} = 3.424$. By applying $λ_{DS}$ with the same parameters from the DMO run to the hydro-dynamically simulated clusters (Gadget-X run and GIZMO-SIMBA run), we investigate the effect of baryons on the cluster dynamical state. We find a weak baryon-model dependence for the $λ_{DS}$ parameter. Finally, we study the evolution of $λ_{DS}$ along with clusters mass accretion history. We notice an upper limit of halo mass change $\frac{ΔM_{200}}{M_{200}} \sim 0.12$ that don't alter the cluster dynamical state, i.e. from relaxed to unrelaxed. We define relaxation period (from the most relaxed state to disturb and relaxed again) which reflects how long the dynamical state of a cluster restores its relaxation state, and propose a correlation between this relaxation period and the strength of halo mass change $\frac{ΔM_{200}}{M_{200}}$. With the proposed fitting to such correlation, we verify the relaxation period can be estimated from $\frac{ΔM_{200}}{M_{200}}$ (including multi mass change peaks) with considerably small error.

preprint2020arXiv

Predictions of the L$_{\rm[CII]}$-SFR and [C$_{\rm II}$] Luminosity Function at the Epoch of Reionization

We present the first predictions for the $L_{\rm [CII]}$ - SFR relation and [CII] luminosity function (LF) in the Epoch of Reionization (EoR) based on cosmological hydrodynamics simulations using the SIMBA suite plus radiative transfer calculations via SÍGAME. The sample consists of 11,137 galaxies covering halo mass $\log M_{\rm halo}\in$[9, 12.4] $M_\odot$, star formation rate SFR$\in$[0.01, 330] $M_\odot$ yr$^{-1}$, and metallicity $<Z_{\rm gas}>_{\rm SFR}\in$[0.1, 0.9] $Z_\odot$. The simulated $L_{\rm [CII]}$-SFR relation is consistent with the range observed, but with a spread of $\simeq$0.3 dex at the high end of SFR ($>$100 $M_\odot$ yr$^{-1}$) and $\simeq$0.6 dex at the lower end, and there is tension between our predictions and the values of $L_{\rm [CII]}$ above 10$^{8.5}$ $L_\odot$ observed in some galaxies reported in the literature. The scatter in the $L_{\rm [CII]}$-SFR relation is mostly driven by galaxy properties, such that at a given SFR, galaxies with higher molecular gas mass and metallicity have higher $L_{\rm [CII]}$. The [CII] LF predicted by SIMBA is consistent with the upper limits placed by the only existing untargeted flux-limited [CII] survey at the EoR (ASPECS) and those predicted by semi-analytic models. We compare our results with existing models and discuss differences responsible for the discrepant slopes in the $L_{\rm [CII]}$-SFR relatiion.

preprint2020arXiv

Protoclusters at z=5.7: A view from the MultiDark galaxies

Protoclusters, which will yield galaxy clusters at lower redshift, can provide valuable information on the formation of galaxy clusters. However, identifying progenitors of galaxy clusters in observations is not an easy task, especially at high redshift. Different priors have been used to estimate the overdense regions that are thought to mark the locations of protoclusters. In this paper, we use mimicked Ly$α$-emitting galaxies at $z=5.7$ to identify protoclusters in the MultiDark galaxies, which are populated by applying three different semi-analytic models to the 1 $Gpc h^{-1}$ MultiDark Planck2 simulation. To compare with observational results, we extend the criterion 1 (a Ly$α$ luminosity limited sample), to criterion 2 (a match to the observed mean galaxy number density). To further statistically study the finding efficiency of this method, we enlarge the identified protocluster sample (criterion 3) to about 3500 at $z=5.7$ and study their final mass distribution. The number of overdense regions and their selection probability depends on the semi-analytic models and strongly on the three selection criteria (partly by design). The protoclusters identified with criterion 1 are associated with a typical final cluster mass of $2.82\pm0.92 \times 10^{15} M_\odot$ which is in agreement with the prediction (within $\pm 1 σ$) of an observed massive protocluster at $z=5.7$. Identifying more protoclusters allows us to investigate the efficiency of this method, which is more suitable for identifying the most massive clusters: completeness ($\mathbb{C}$) drops rapidly with decreasing halo mass. We further find that it is hard to have a high purity ($\mathbb{P}$) and completeness simultaneously.

preprint2020arXiv

The Three Hundred Project: the stellar and gas profiles

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

preprint2019arXiv

And yet it flips: connecting galactic spin and the cosmic web

We study the spin alignment of galaxies and halos with respect to filaments and walls of the cosmic web, identified with DisPerSE, using the SIMBA simulation from z=0-2. Massive halos' spins are oriented perpendicularly to their closest filament's axis and walls, while low mass halos tend to have their spins parallel to filaments and in the plane of walls. A similar mass-dependent spin flip is found for galaxies, albeit with a weaker signal particularly at low mass and low-z, suggesting that galaxies' spins retain memory of their larger-scale environment. Low-z star-forming and rotation-dominated galaxies tend to have spins parallel to nearby filaments, while quiescent and dispersion-dominated galaxies show preferentially perpendicular orientation; the star formation trend can be fully explained by the stellar mass correlation, but the morphology trend cannot. There is a strong dependence on HI mass, such that high-HI galaxies tend to have parallel spins while low-HI galaxies are perpendicular, which persists even when matching samples in stellar mass, suggesting that HI content traces anisotropic infall more faithfully than the stellar component. Finally, at fixed stellar mass, the strength of spin alignments correlates with the filament's density, with parallel alignment for galaxies in high density environments. These findings are consistent with conditional tidal torque theory, and highlight a significant correlation between galactic spin and the larger scale tides that are important e.g. for interpreting weak lensing studies. SIMBA allows us to rule out numerical grid locking as the cause of previously-seen low mass alignment.

preprint2019arXiv

Cosmological baryon transfer in the SIMBA simulations

We present a framework for characterizing the large scale movement of baryons relative to dark matter in cosmological simulations, requiring only the initial conditions and final state of the simulation. This is performed using the spread metric which quantifies the distance in the final conditions between initially neighbouring particles, and by analysing the baryonic content of final haloes relative to that of the initial Lagrangian regions defined by their dark matter component. Applying this framework to the SIMBA cosmological simulations, we show that 40% (10%) of cosmological baryons have moved $> 1h^{-1}~ {\rm Mpc}^{-1}$ ($3h^{-1}~ {\rm Mpc}^{-1}$) by $z=0$, due primarily to entrainment of gas by jets powered by AGN, with baryons moving up to $12h^{-1}~ {\rm Mpc}^{-1}$ away in extreme cases. Baryons decouple from the dynamics of the dark matter component due to hydrodynamic forces, radiative cooling, and feedback processes. As a result, only 60% of the gas content in a given halo at $z=0$ originates from its Lagrangian region, roughly independent of halo mass. A typical halo in the mass range $M_{\rm vir} = 10^{12}$--$10^{13}{\rm M}_\odot$ only retains 20% of the gas originally contained in its Lagrangian region. We show that up to 20% of the gas content in a typical Milky Way mass halo may originate in the region defined by the dark matter of another halo. This inter-Lagrangian baryon transfer may have important implications for the origin of gas and metals in the circumgalactic medium of galaxies, as well as for semi-analytic models of galaxy formation and "zoom-in" simulations.

preprint2016arXiv

The Properties of the Circumgalactic Medium in Red and Blue Galaxies: Results from the COS-GASS+COS-Halos Surveys

We use the combined data from the COS-GASS and COS-Halos surveys to characterize the Circum-Galactic Medium (CGM) surrounding typical low-redshift galaxies in the mass range $\rm~M_*\sim~10^{9.5-11.5}~M_{\odot} $, and over a range of impact parameters extending to just beyond the halo virial radius ($\rm~R_{vir}$). We find the radial scale length of the distributions of the equivalent widths of the Lyman~$α$ and Si III absorbers to be 0.9 and 0.4 $\rm~R_{vir}$, respectively. The radial distribution of equivalent widths is relatively uniform for the blue galaxies, but highly patchy (low covering fraction) for the red galaxies. We also find that the Lyman~$α$ and Si III equivalent widths show significant positive correlations with the specific star-formation rate (sSFR) of the galaxy. We find a surprising lack of correlations between the halo mass (virial velocity) and either the velocity dispersions or velocity offsets of the Lyman~$α$ lines. The ratio of the velocity offset to the velocity dispersion for the Lyman~$α$ absorbers has a mean value of $\sim$ 4, suggesting that a given the line-of-sight is intersecting a dynamically coherent structure in the CGM rather than a sea of orbiting clouds. The kinematic properties of the CGM are similar in the blue and red galaxies, although we find that a significantly larger fraction of the blue galaxies have large Lyman~$α$ velocity offsets (>200 km s$^{-1}$). We show that - if the CGM clouds represent future fuel for star-formation - our new results could imply a large drop in the specific star-formation rate across the galaxy mass-range we probe.

preprint2015arXiv

Baryon Cycling in the Low-Redshift Circumgalactic Medium: A Comparison of Simulations to the COS-Halos Survey

We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass.

preprint2015arXiv

Connection Between the Circumgalactic Medium and the Interstellar Medium of Galaxies: Results from the COS-GASS Survey

We present a study exploring the nature and properties of the Circum-Galactic Medium (CGM) and its connection to the atomic gas content in the interstellar medium (ISM) of galaxies as traced by the HI 21cm line. Our sample includes 45 low-z (0.026-0.049) galaxies from the GALEX Arecibo SDSS Survey. Their CGM was probed via absorption in the spectra of background Quasi-Stellar Objects at impact parameters of 63 to 231kpc. The spectra were obtained with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. We detected neutral hydrogen (Ly$α$ absorption-lines) in the CGM of 92% of the galaxies. We find the radial profile of the CGM as traced by the Ly$α$ equivalent width can be fit as an exponential with a scale length of roughly the virial radius of the dark matter halo. We found no correlation between the orientation of sightline relative to the galaxy major axis and the Ly$α$ equivalent width. The velocity spread of the circumgalactic gas is consistent with that seen in the atomic gas in the interstellar medium. We find a strong correlation (99.8% confidence) between the gas fraction (M(HI)/M*) and the impact-parameter-corrected Ly$α$ equivalent width. This is stronger than the analogous correlation between corrected Ly$α$ equivalent width and SFR/M* (97.5% confidence). These results imply a physical connection between the HI disk and the CGM, which is on scales an order-of-magnitude larger. This is consistent with the picture in which the HI disk is nourished by accretion of gas from the CGM.

preprint2015arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at $z\ge6$. III. Stellar Populations from SED Modeling with Secure Ly$α$ Emission and Redshifts

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$α$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$α$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$α$ emission are indistinguishable in terms of age, stellar mass, and SFR.

preprint2015arXiv

The MOSFIRE Deep Evolution Field (MOSDEF) Survey: Rest-Frame Optical Spectroscopy for ~1500 H-Selected Galaxies at 1.37 < z < 3.8

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

preprint2014arXiv

Constraining Very High Mass Population III Stars through He II Emission in Galaxy BDF-521 at z = 7.01

Numerous theoretical models have long proposed that a strong He II 1640 emission line is the most prominent and unique feature of massive Population III (Pop III) stars in high redshift galaxies. The He II 1640 line strength can constrain the mass and IMF of Pop III stars. We use F132N narrowband filter on the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) to look for strong He II lambda 1640 emission in the galaxy BDF-521 at z=7.01, one of the most distant spectroscopically-confirmed galaxies to date. Using deep F132N narrowband imaging, together with our broadband imaging with F125W and F160W filters, we do not detect He II emission from this galaxy, but place a 2-sigma upper limit on the flux of 5.3x10^-19 ergs s^-1 cm^-2. This measurement corresponds to a 2-sigma upper limit on the Pop III star formation rate (SFR_PopIII) of ~ 0.2 M_solar yr^-1, assuming a Salpeter IMF with 50< M/M_solar < 1000. From the high signal-to-noise broadband measurements in F125W and F160W, we fit the UV continuum for BDF-521. The spectral flux density is ~ 3.6x 10^-11 lambda^-2.32 ergs s^-1 cm^-2 A^-1, which corresponds to an overall unobscured SFR of ~ 5 M_solar yr^-1. Our upper limit on SFR_PopIII suggests that massive Pop III stars represent < 4% of the total star formation. Further, the HST high resolution imaging suggests that BDF-521 is an extremely compact galaxy, with a half-light radius of 0.6 kpc.

preprint2014arXiv

Optical-faint, Far-infrared-bright Herschel Sources in the CANDELS Fields: Ultra-Luminous Infrared Galaxies at z>1 and the Effect of Source Blending

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

preprint2014arXiv

Parametrising Star Formation Histories

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

preprint2014arXiv

Properties of Submillimeter Galaxies in the CANDELS GOODS-S Field

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

preprint2013arXiv

Empirical Constraints for the Magnitude and Composition of Galactic Winds

Galactic winds are a key physical mechanism for understanding galaxy formation and evolution, yet empirical and theoretical constraints for the character of winds are limited and discrepant. Recent empirical models find that local star-forming galaxies have a deficit of oxygen that scales with galaxy stellar mass. The oxygen deficit provides unique empirical constraints on the magnitude of mass loss, composition of outflowing material and metal reaccretion onto galaxies. We formulate the oxygen deficit constraints so they may be easily implemented into theoretical models of galaxy evolution. We parameterize an effective metal loading factor which combines the uncertainties of metal outflows and metal reaccretion into a single function of galaxy virial velocity. We determine the effective metal loading factor by forward-fitting the oxygen deficit. The effective metal loading factor we derive has important implications for the implementation of mass loss in models of galaxy evolution.

preprint2013arXiv

On the connection between the intergalactic medium and galaxies: The HI-galaxy cross-correlation at z < 1

We present a new optical spectroscopic survey of 1777 'star-forming' ('SF') and 366 'non-star-forming' ('non-SF') galaxies at redshifts z < 1 (2143 in total), 22 AGN and 423 stars, observed by instruments such as DEIMOS, VIMOS and GMOS, in 3 fields containing 5 quasi-stellar objects (QSOs) with HST UV spectroscopy. We also present a new spectroscopic survey of 165 'strong' (10^14 < NHI < 10^17 cm^-2), and 489 'weak' (10^13 < NHI < 10^14 cm^-2) intervening HI absorption line systems at z < 1 (654 in total), observed in the spectra of 8 QSOs by COS and FOS on the HST. Combining these new data with previously published galaxy catalogs such as VVDS and GDDS, we have gathered a sample of 654 HI absorption systems and 17509 galaxies at transverse scales < 50 Mpc. We present observational results on the HI-galaxy and galaxy-galaxy correlations at transverse scales r < 10 Mpc, and the HI-HI auto-correlation at transverse scales r < 2 Mpc. The two-point correlation functions are measured both along and transverse to the line-of-sight. We constrain the HI-galaxy statistical connection, as a function of both HI column density and galaxy star-forming activity. Our results are consistent with the following conclusions: (1) the bulk of HI systems on Mpc scales have little velocity dispersion (<120 km/s) with respect to the bulk of galaxies; (2) the vast majority of strong HI systems and SF galaxies are distributed in the same locations, together with 75+-15% of non-SF galaxies, all of which typically reside in dark matter haloes of similar masses; (3) 25+-15% of non-SF galaxies reside in galaxy clusters and are not correlated with strong HI systems at scales < 2 Mpc; and (4) 50% of weak HI systems reside within galaxy voids (hence not correlated with galaxies), and are confined in dark matter haloes of masses smaller than those hosting... [abridged]

preprint2013arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. I. Basic Characteristics of the Rest-Frame UV Continuum and Lyman-alpha emission

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 Å. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 Å) that could be the dominant contribution to the total ionizing flux at z >= 6.

preprint2013arXiv

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. II. Morphology of the Rest-Frame UV Continuum and Lyman-alpha Emission

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

preprint2013arXiv

The COS-Halos Survey: Rationale, Design, and A Census of Circumgalactic Neutral Hydrogen

We present the design and methods of the COS-Halos survey, a systematic investigation of the gaseous halos of 44 z = 0.15-0.35 galaxies using background QSOs observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. This survey has yielded 39 spectra of z_em ~ 0.5 QSOs with S/N ~ 10-15 per resolution element. The QSO sightlines pass within 150 physical kpc of the galaxies, which span early and late types over stellar mass log M* / Msun= 9.5 - 11.5. We find that the CGM exhibits strong HI, averaging 1 Ang in Lya equivalent width out to 150 kpc, with 100% covering fraction for star-forming galaxies and 75% covering for passive galaxies. We find good agreement in column densities between this survey and previous studies over similar range of impact parameter. There is weak evidence for a difference between early- and late-type galaxies in the strength and distribution of HI. Kinematics indicate that the detected material is bound to the host galaxy, such that >~90% of the detected column density is confined within +/-200 km s^-1 of the galaxies. This material generally exists well below the halo virial temperatures at T<~ 10^5 K. We evaluate a number of possible origin scenarios for the detected material, and in the end favor a simple model in which the bulk of the detected HI arises in a bound, cool, low-density photoionized diffuse medium that is generic to all L* galaxies and may harbor a total gaseous mass comparable to galactic stellar masses.

preprint2013arXiv

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

preprint2013arXiv

The High-Ion Content and Kinematics of Low-Redshift Lyman Limit Systems

We study the high-ionization phase and kinematics of the circumgalactic medium around low-redshift galaxies using a sample of 23 Lyman Limit Systems (LLSs) at 0.08<z<0.93 observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. In Lehner et al. (2013), we recently showed that low-z LLSs have a bimodal metallicity distribution. Here we extend that analysis to search for differences between the high-ion and kinematic properties of the metal-poor and metal-rich branches. We find that metal-rich LLSs tend to show higher O VI columns and broader O VI profiles than metal-poor LLSs. The total H I line width (dv90 statistic) in LLSs is not correlated with metallicity, indicating that the H I kinematics alone cannot be used to distinguish inflow from outflow and gas recycling. Among the 17 LLSs with O VI detections, all but two show evidence of kinematic sub-structure, in the form of O VI-H I centroid offsets, multiple components, or both. Using various scenarios for how the metallicity in the high-ion and low-ion phases of each LLS compare, we constrain the ionized hydrogen column in the O VI phase to lie in the range log N(H II)~17.6-20. The O VI phase of LLSs is a substantial baryon reservoir, with M(high-ion)~10^{8.5-10.9}(r/150 kpc)^2 solar masses, similar to the mass in the low-ion phase. Accounting for the O VI phase approximately doubles the contribution of low-z LLSs to the cosmic baryon budget.

preprint2013arXiv

The LBT Bootes Field Survey: I. The Rest-frame UV and Luminosity Functions and Clustering of Bright Lyman Break Galaxies at z~3

We present a deep LBT/LBC U-band imaging survey (9 deg2) covering the NOAO Bootes field. A total of 14,485 Lyman Break Galaxies (LBGs) at z~3 are selected, which are used to measure the rest-frame UV luminosity function (LF). The large sample size and survey area reduce the LF uncertainties due to Poisson statistics and cosmic variance by >3 compared to previous studies. At the bright end, the LF shows excess power compared to the best-fit Schechter function, which can be attributed to the contribution of $z\sim3$ quasars. We compute the rest-frame near-infrared LF and stellar mass function (SMF) of z~3 LBGs based on the R-band and IRAC [4.5 micro m]-band flux relation. We investigate the evolution of the UV LFs and SMFs between z~7 and z~3, which supports a rising star formation history in the LBGs. We study the spatial correlation function of two bright LBG samples and estimate their average host halo mass. We find a tight relation between the host halo mass and the galaxy star formation rate (SFR),which follows the trend predicted by the baryonic accretion rate onto the halo, suggesting that the star formation in LBGs is fueled by baryonic accretion through the cosmic web. By comparing the SFRs with the total baryonic accretion rates, we find that cosmic star formation efficiency is about 5%-20% and it does not evolve significantly with redshift, halo mass, or galaxy luminosity.

preprint2013arXiv

The Redshift and Mass Dependence on the Formation of The Hubble Sequence at z>1 from CANDELS/UDS

In this paper we present a detailed study of the structures and morphologies of a sample of 1188 massive galaxies with Mstar>10^10Msun between redshifts z=1-3 within the Ultra Deep Survey (UDS) region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. Using this sample we determine how galaxy structure and morphology evolve with time. We visually classify our sample into disks, ellipticals and peculiar systems and correct for redshift effects on our classifications through simulations. We find evolution in the fractions of galaxies at a given visual classification as a function of redshift. The peculiar population is dominant at z>2 with a substantial spheroid population, and a negligible disk population. We compute the transition redshift, ztrans, where the combined fraction of spheroids and disks is equal to that of peculiar galaxies, as ztrans=1.86+/-0.62 for galaxies in our stellar mass range. We find that this changes as a function of stellar mass, with Hubble-type systems becoming dominant at higher redshifts for higher mass galaxies (ztrans=2.22+/-0.82), than for the lower mass galaxies (ztrans=1.73+/-0.57). Higher mass galaxies become morphologically settled before their lower mass counterparts, a form of morphological downsizing. We furthermore compare our visual classifications with Sersic index, the concentration, asymmetry and clumpiness (CAS) parameters, star formation rate and rest frame U-B colour. We find links between the colour of a galaxy, its star formation rate and how extended or peculiar it appears. Finally, we discuss the negligible z>2 disk fraction based on visual morphologies and speculate that this is an effect of forming disks appearing peculiar through processes such as violent disk instabilities or mergers. We conclude that to properly define high redshift morphology a new and more exact classification scheme is needed.

preprint2013arXiv

The Slow Flow Model of Dust Efflux in Local Star-Forming Galaxies

We develop a dust efflux model of radiation pressure acting on dust grains which successfully reproduces the relation between stellar mass, dust opacity and star formation rate observed in local star-forming galaxies. The dust content of local star-forming galaxies is set by the competition between the physical processes of dust production and dust loss in our model. The dust loss rate is proportional to the dust opacity and star formation rate. Observations of the relation between stellar mass and star formation rate at several epochs imply that the majority of local star-forming galaxies are best characterized as having continuous star formation histories. Dust loss is a consequence of sustained interaction of dust with the radiation field generated by continuous star formation. Dust efflux driven by radiation pressure rather than dust destruction offers a more consistent physical interpretation of the dust loss mechanism. By comparing our model results with the observed relation between stellar mass, dust extinction and star formation rate in local star-forming galaxies we are able to constrain the timescale and magnitude of dust loss. The timescale of dust loss is long and therefore dust is effluxed in a "Slow Flow". Dust loss is modest in low mass galaxies but massive galaxies may lose up to 70~80% of their dust over their lifetime. Our Slow Flow model shows that mass loss driven by dust opacity and star formation may be an important physical process for understanding normal star-forming galaxy evolution.

preprint2013arXiv

Theoretical Evolution of Optical Strong Lines across Cosmic Time

We use the chemical evolution predictions of cosmological hydrodynamic simulations with our latest theoretical stellar population synthesis, photoionization and shock models to predict the strong line evolution of ensembles of galaxies from z=3 to the present day. In this paper, we focus on the brightest optical emission-line ratios, [NII]/H-alpha and [OIII]/H-beta. We use the optical diagnostic Baldwin-Phillips-Terlevich (BPT) diagram as a tool for investigating the spectral properties of ensembles of active galaxies. We use four redshift windows chosen to exploit new near-infrared multi-object spectrographs. We predict how the BPT diagram will appear in these four redshift windows given different sets of assumptions. We show that the position of star-forming galaxies on the BPT diagram traces the ISM conditions and radiation field in galaxies at a given redshift. Galaxies containing AGN form a mixing sequence with purely star-forming galaxies. This mixing sequence may change dramatically with cosmic time, due to the metallicity sensitivity of the optical emission-lines. Furthermore, the position of the mixing sequence may probe metallicity gradients in galaxies as a function of redshift, depending on the size of the AGN narrow line region. We apply our latest slow shock models for gas shocked by galactic-scale winds. We show that at high redshift, galactic wind shocks are clearly separated from AGN in line ratio space. Instead, shocks from galactic winds mimic high metallicity starburst galaxies. We discuss our models in the context of future large near-infrared spectroscopic surveys.

preprint2012arXiv

A Census of Oxygen in Star-Forming Galaxies: An Empirical Model Linking Metallicities, Star Formation Rates and Outflows

In this contribution we present the first census of oxygen in star-forming galaxies in the local universe. We examine three samples of galaxies with metallicities and star formation rates at z = 0.07, 0.8 and 2.26, including the SDSS and DEEP2 surveys. We infer the total mass of oxygen produced and mass of oxygen found in the gas-phase from our local SDSS sample. The star formation history is determined by requiring that galaxies evolve along the relation between stellar mass and star formation rate observed in our three samples. We show that the observed relation between stellar mass and star formation rate for our three samples is consistent with other samples in the literature. The mass-metallicity relation is well established for our three samples and from this we empirically determine the chemical evolution of star-forming galaxies. Thus, we are able to simultaneously constrain the star formation rates and metallicities of galaxies over cosmic time allowing us to estimate the mass of oxygen locked up in stars. Combining this work with independent measurements reported in the literature we conclude that the loss of oxygen from the interstellar medium of local star-forming galaxies is likely to be a ubiquitous process with the oxygen mass loss scaling (almost) linearly with stellar mass. We estimate the total baryonic mass loss and argue that only a small fraction of the baryons inferred from cosmological observations accrete onto galaxies.

preprint2012arXiv

Galaxy Gas Fractions at High-Redshift: The Tension between Observations and Cosmological Simulations

CO measurements of z~1-4 galaxies have found that their baryonic gas fractions are significantly higher than galaxies at z=0, with values ranging from 20-80 %. Here, we suggest that the gas fractions inferred from observations of star-forming galaxies at high-z are overestimated, owing to the adoption of locally-calibrated CO-H2 conversion factors (Xco). Evidence from both observations and numerical models suggest that Xco varies smoothly with the physical properties of galaxies, and that Xco can be parameterised simply as a function of both gas phase metallicity and observed CO surface brightness. When applying this functional form, we find fgas ~10-40 % in galaxies with M*=10^10-10^12 Msun at high-z. Moreover, the scatter in the observed fgas-M* relation is lowered by a factor of two. The lower inferred gas fractions arise physically because the interstellar media of high-z galaxies have higher velocity dispersions and gas temperatures than their local counterparts, which results in an Xco that is lower than the z=0 value for both quiescent discs and starbursts. We further compare these gas fractions to those predicted by cosmological galaxy formation models. We show that while the canonically inferred gas fractions from observations are a factor of 2-3 larger at a given stellar mass than predicted by models, our rederived Xco values for z=1-4 galaxies results in revised gas fractions that agree significantly better with the simulations.

preprint2012arXiv

Luminous and High Stellar Mass Candidate Galaxies at z ~ 8 Discovered in CANDELS

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

preprint2012arXiv

Not Dead Yet: Cool Circumgalactic Gas in the Halos of Early Type Galaxies

We report new observations of circumgalactic gas in the halos of early type galaxies obtained by the COS-Halos Survey with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We find that detections of HI surrounding early type galaxies are typically as common and strong as around star-forming galaxies, implying that the total mass of circumgalactic material is comparable in the two populations. For early type galaxies, the covering fraction for HI absorption above 10^16 cm^2 is ~40-50% within ~150 kpc. Line widths and kinematics of the detected material show it to be cold (T ~< 10^5 K) in comparison to the virial temperature of the host halos. The implied masses of cool, photoionized CGM baryons may be up to 10^9 --- 10^11 Msun. Contrary to some theoretical expectations, strong halo HI absorbers do not disappear as part of the quenching of star-formation. Even passive galaxies retain significant reservoirs of halo baryons which could replenish the interstellar gas reservoir and eventually form stars. This halo gas may feed the diffuse and molecular gas that is frequently observed inside ETGs.

preprint2012arXiv

Testing Subhalo Abundance Matching in Cosmological Smoothed Particle Hydrodynamics Simulations

Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearly monotonic relation between stellar mass and halo mass provided that, for satellite galaxies, we use the mass of the subhalo at the epoch when it became a satellite. In each simulation, the median relation for central and satellite galaxies is nearly identical, though a somewhat larger fraction of satellites are outliers. SHAM-assigned masses (at z=0-2), luminosities (R-band at z=0), or star formation rates (at z=2) have a 68% scatter of 0.09-0.15 dex relative to the true simulation values. When we apply SHAM to the subhalo population of collisionless N-body simulation with the same initial conditions as the SPH runs, we find generally good agreement for the halo occupation distributions and halo radial profiles of galaxy samples defined by thresholds in stellar mass. However, because a small fraction of SPH galaxies suffer severe stellar mass loss after becoming satellites, SHAM slightly overpopulates high mass halos; this effect is more significant for the wind simulation, which produces galaxies that are less massive and more fragile. SHAM recovers the two-point correlation function of the SPH galaxies in the no-wind simulation to better than 10% at scales 0.1 < r < 10 Mpc/h. For the wind simulation, agreement is better than 15% at r > 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales.

preprint2011arXiv

Probing Population III Stars in Galaxy IOK-1 at z = 6.96 through He II Emission

The He II λ1640 emission line has been suggested as a direct probe of Population III (Pop III) stars at high-redshift, since it can arise from highly energetic ionizing photons associated with hot, metal free stars. We use the HST WFC3/F130N IR narrowband filter to probe He II λ1640 emission in galaxy IOK-1 at z=6.96. The sensitivity of this measurement is >5x deeper than for previous measurements. From this deep narrowband imaging, combined with broadband observations in the F125W and F160W filters, we find the He II flux to be 1.2+/- 1.0x 10^-18 ergs/s/cm^2, corresponding to a 1σupper limit on the Pop III star formation rate (SFR) of ~ 0.5 M_sun/yr for the case of a Salpeter IMF with 50-500M_sun and mass loss. Given that the broadband measurements can be fit with a UV continuum spectral flux density of ~ 4.85x 10^-10x λ^-2.46 ergs/s/cm^2/A, which corresponds to an overall SFR of ~16+/-2.6 M_sun/yr, massive Pop III stars represent < 6% of the total star formation. This measurement places the strongest limit yet on metal-free star formation at high redshift, although the exact conversion from He II luminosity to Pop III SFR is highly uncertain due to the unknown IMF, stellar evolution, and photoionization effects. Although we have not detected He II λ1640 at more than the 1.2σlevel, our work suggests that a > 3σlevel detection is possible with JWST.

preprint2011arXiv

The Large, Oxygen-Rich Halos of Star-Forming Galaxies Are A Major Reservoir of Galactic Metals

The circumgalactic medium (CGM) is fed by galaxy outflows and accretion of intergalactic gas, but its mass, heavy element enrichment, and relation to galaxy properties are poorly constrained by observations. In a survey of the outskirts of 42 galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, we detected ubiquitous, large (150 kiloparsec) halos of ionized oxygen surrounding star-forming galaxies, but we find much less ionized oxygen around galaxies with little or no star formation. This ionized CGM contains a substantial mass of heavy elements and gas, perhaps far exceeding the reservoirs of gas in the galaxies themselves. It is a basic component of nearly all star-forming galaxies that is removed or transformed during the quenching of star formation and the transition to passive evolution.

preprint2010arXiv

Imaging the molecular gas in a submm galaxy at z = 4.05: cold mode accretion or a major merger?

We present a high resolution (down to 0.18"), multi-transition imaging study of the molecular gas in the z = 4.05 submillimeter galaxy GN20. GN20 is one of the most luminous starburst galaxy known at z > 4, and is a member of a rich proto-cluster of galaxies at z = 4.05 in GOODS-North. We have observed the CO 1-0 and 2-1 emission with the VLA, the CO 6-5 emission with the PdBI Interferometer, and the 5-4 emission with CARMA. The H_2 mass derived from the CO 1-0 emission is 1.3 \times 10^{11} (α/0.8) Mo. High resolution imaging of CO 2-1 shows emission distributed over a large area, appearing as partial ring, or disk, of ~ 10kpc diameter. The integrated CO excitation is higher than found in the inner disk of the Milky Way, but lower than that seen in high redshift quasar host galaxies and low redshift starburst nuclei. The VLA CO 2-1 image at 0.2" resolution shows resolved, clumpy structure, with a few brighter clumps with intrinsic sizes ~ 2 kpc. The velocity field determined from the CO 6-5 emission is consistent with a rotating disk with a rotation velocity of ~ 570 km s^{-1} (using an inclination angle of 45^o), from which we derive a dynamical mass of 3 \times 10^{11} \msun within about 4 kpc radius. The star formation distribution, as derived from imaging of the radio synchrotron and dust continuum, is on a similar scale as the molecular gas distribution. The molecular gas and star formation are offset by ~ 1" from the HST I-band emission, implying that the regions of most intense star formation are highly dust-obscured on a scale of ~ 10 kpc. The large spatial extent and ordered rotation of this object suggests that this is not a major merger, but rather a clumpy disk accreting gas rapidly in minor mergers or smoothly from the proto-intracluster medium. ABSTRACT TRUNCATED

preprint2009arXiv

Galaxies in a simulated $Λ$CDM Universe I: cold mode and hot cores

We study the formation of galaxies in a (50 Mpc/h)^3 cosmological simulation (2x288^3 particles), evolved using the entropy conserving SPH code Gadget-2. Most of the baryonic mass in galaxies of all masses is originally acquired through filamentary "cold mode" accretion of gas that was never shock heated to its halo virial temperature, confirming the key feature of our earlier results obtained with a different SPH code (Keres et al. 2005). Atmospheres of hot, virialized gas develop in halos above ~2.5e11 Msun, a transition mass that is nearly constant from z=3 to z=0. Cold accretion persists in halos above the transition mass, especially at z>=2. It dominates the growth of galaxies in low mass halos at all times, and it is the main driver of the cosmic star formation history. Satellite galaxies have accretion rates similar to central galaxies of the same baryonic mass at high redshifts, but they have less accretion than comparable central galaxies at low redshift. Relative to our earlier results, the Gadget-2 simulations predict much lower rates of "hot mode" accretion from the virialized gas component of massive halos. At z<=1, typical hot accretion rates in halos above 5e12 Msun are below 1 Msun/yr, even though our simulation does not include AGN heating or other forms of "preventive" feedback. The inner density profiles of hot gas in these halos are shallow, with long associated cooling times. The cooling recipes typically used in semi-analytic models can overestimate the accretion rates in these halos by orders of magnitude, so such models may overemphasize the role of preventive feedback in producing observed galaxy masses and colors. A fraction of the massive halos develop cuspy profiles and significant cooling rates between z=1 and z=0, a redshift trend similar to the observed trend in the frequency of cooling flow clusters.

preprint2009arXiv

Galaxies in a Simulated $Λ$CDM Universe II: Observable Properties and Constraints on Feedback

We compare the properties of galaxies that form in a cosmological simulation without strong feedback to observations at z=0. We confirm previous findings that models without strong feedback overproduce the observed galaxy baryonic mass function, especially at the low and high mass extremes. Through post-processing we investigate what kinds of feedback would be required to reproduce observed galaxy masses and star formation rates. To mimic an extreme form of "preventive" feedback (e.g., AGN radio mode) we remove all baryonic mass that was originally accreted via "hot mode" from shock-heated gas. This does not bring the high mass end of the galaxy mass function into agreement with observations because much of the stellar mass in these systems formed at high redshift from baryons that originally accreted via "cold mode" onto lower mass progenitors. An efficient "ejective" feedback mechanism, such as supernova driven winds, must reduce the masses of these progenitors. Feedback must also reduce the masses of lower mass z=0 galaxies, which assemble at lower redshifts and have much lower star formation rates. If we monotonically re-map galaxy masses to reproduce the observed mass function, but retain the simulation's predicted star formation rates, we obtain fairly good agreement with the observed sequence of star-forming galaxies but fail to recover the observed population of passive, low star formation rate galaxies. Suppressing all hot mode accretion improves agreement for high mass galaxies but worsens the agreement at intermediate masses. Reproducing these z=0 observations requires a feedback mechanism that dramatically suppresses star formation in a fraction of galaxies, increasing with mass, while leaving star formation rates of other galaxies essentially unchanged.

preprint2009arXiv

The Simulated HI Sky at low redshift

Observations of intergalactic neutral hydrogen can provide a wealth of information about structure and galaxy formation, potentially tracing accretion and feedback processes on Mpc scales. Below a column density of NHI ~ 10^19 cm-2, the "edge" or typical observational limit for HI emission from galaxies, simulations predict a cosmic web of extended emission and filamentary structures. We study the distribution of neutral hydrogen and its 21cm emission properties in a cosmological hydrodynamic simulation, to gain more insights into the distribution of HI below NHI ~ 10^19 cm-2. Such Lyman Limit systems are expected to trace out the cosmic web, and are relatively unexplored. Beginning with a 32 h^-1 Mpc simulation, we extract the neutral hydrogen component by determining the neutral fraction, including a post-processed correction for self-shielding based on the thermal pressure. We take into account molecular hydrogen, assuming an average density ratio Omega_H2 / Omega_HI = 0.3 at z = 0. The statistical properties of the HI emission are compared with observations, to assess the reliability of the simulation. The simulated HI distribution robustly describes the full column density range between NHI ~ 10^14 and NHI ~ 10^21 cm-2 and agrees very well with available measurements from observations. Furthermore there is good correspondence in the statistics when looking at the two-point correlation function and the HI mass function. The reconstructed maps are used to simulate observations of existing and future telescopes by adding noise and taking account of the sensitivity of the telescopes. The general agreement in statistical properties of HI suggests that neutral hydrogen as modeled in this hydrodynamic simulation is a fair representation of that in the Universe. (abridged)

preprint2008arXiv

The Flux Auto- and Cross-Correlation of the Lyman-alpha Forest. II. Modelling Anisotropies with Cosmological Hydrodynamic Simulations

The isotropy of the Lyman-alpha forest in real-space uniquely provides a measurement of cosmic geometry at z > 2. The angular diameter distance for which the correlation function along the line of sight and in the transverse direction agree corresponds to the correct cosmological model. However, the Lyman-alpha forest is observed in redshift-space where distortions due to Hubble expansion, bulk flows, and thermal broadening introduce anisotropy. Similarly, a spectrograph's line spread function affects the autocorrelation and cross-correlation differently. In this the second paper of a series on using the Lyman-alpha forest observed in pairs of QSOs for a new application of the Alcock-Paczynski (AP) test, these anisotropies and related sources of potential systematic error are investigated with cosmological hydrodynamic simulations. Three prescriptions for galactic outflow were compared and found to have only a marginal effect on the Lyman-alpha flux correlation (which changed by at most 7% with use of the currently favored variable-momentum wind model vs. no winds at all). An approximate solution for obtaining the zero-lag cross-correlation corresponding to arbitrary spectral resolution directly from the zero-lag cross-correlation computed at full-resolution (good to within 2% at the scales of interest) is presented. Uncertainty in the observationally determined mean flux decrement of the Lyman-alpha forest was found to be the dominant source of systematic error; however, this is reduced significantly when considering correlation ratios. We describe a simple scheme for implementing our results, while mitigating systematic errors, in the context of a future application of the AP test.

preprint2008arXiv

The Growth of Central and Satellite Galaxies in Cosmological Smoothed Particle Hydrodynamics Simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7e9 M[Sun]. Most friends-of-friends halos in the simulation have a distinct central galaxy, typically two to five times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive halos form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Since z=1, 27% of central galaxies (above 3e10 M[Sun]) and 22% of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain "central" objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with the observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.

preprint2002arXiv

How Do Galaxies Get Their Gas?

Not the way one might have thought. In hydrodynamic simulations of galaxy formation, some gas follows the traditionally envisioned route, shock heating to the halo virial temperature before cooling to the much lower temperature of the neutral ISM. But most gas enters galaxies without ever heating close to the virial temperature, gaining thermal energy from weak shocks and adiabatic compression, and radiating it just as quickly. This ``cold mode'' accretion is channeled along filaments, while the conventional, ``hot mode'' accretion is quasi-spherical. Cold mode accretion dominates high redshift growth by a substantial factor, while at z<1 the overall accretion rate declines and hot mode accretion has greater relative importance. The decline of the cosmic star formation rate at low z is driven largely by geometry, as the typical cross section of filaments begins to exceed that of the galaxies at their intersections.

Romeel Dave

What is connected

Connect this record

See the researcher in context

Building this map preview

52 published item(s)

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

Co-evolution of massive black holes and their host galaxies at high redshift: discrepancies from six cosmological simulations and the key role of JWST

Effects of Active Galactic Nucleus Feedback on Cold Gas Depletion and Quenching of Central Galaxies

How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-Dust Geometry Model in SED Fitting

Rapidly quenched galaxies in the Simba cosmological simulation and observations

The black hole population in low-mass galaxies in large-scale cosmological simulations

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

The Dwarf Galaxy Population at $z\sim 0.7$: A Catalog of Emission Lines and Redshifts from Deep Keck Observations

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

Detection of a Multiphase Intragroup Medium: Results from the COS-IGrM Survey

Finding universal relations in subhalo properties with artificial intelligence

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

The Three Hundred: Cluster Dynamical States and Relaxation Time Scale

Predictions of the L$_{\rm[CII]}$-SFR and [C$_{\rm II}$] Luminosity Function at the Epoch of Reionization

Protoclusters at z=5.7: A view from the MultiDark galaxies

The Three Hundred Project: the stellar and gas profiles

And yet it flips: connecting galactic spin and the cosmic web

Cosmological baryon transfer in the SIMBA simulations

The Properties of the Circumgalactic Medium in Red and Blue Galaxies: Results from the COS-GASS+COS-Halos Surveys

Baryon Cycling in the Low-Redshift Circumgalactic Medium: A Comparison of Simulations to the COS-Halos Survey

Connection Between the Circumgalactic Medium and the Interstellar Medium of Galaxies: Results from the COS-GASS Survey

Physical Properties of Spectroscopically-Confirmed Galaxies at $z\ge6$. III. Stellar Populations from SED Modeling with Secure Ly$α$ Emission and Redshifts

The MOSFIRE Deep Evolution Field (MOSDEF) Survey: Rest-Frame Optical Spectroscopy for ~1500 H-Selected Galaxies at 1.37 < z < 3.8

Constraining Very High Mass Population III Stars through He II Emission in Galaxy BDF-521 at z = 7.01

Optical-faint, Far-infrared-bright Herschel Sources in the CANDELS Fields: Ultra-Luminous Infrared Galaxies at z>1 and the Effect of Source Blending

Parametrising Star Formation Histories

Properties of Submillimeter Galaxies in the CANDELS GOODS-S Field

Empirical Constraints for the Magnitude and Composition of Galactic Winds

On the connection between the intergalactic medium and galaxies: The HI-galaxy cross-correlation at z < 1

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. I. Basic Characteristics of the Rest-Frame UV Continuum and Lyman-alpha emission

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. II. Morphology of the Rest-Frame UV Continuum and Lyman-alpha Emission

The COS-Halos Survey: Rationale, Design, and A Census of Circumgalactic Neutral Hydrogen

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

The High-Ion Content and Kinematics of Low-Redshift Lyman Limit Systems

The LBT Bootes Field Survey: I. The Rest-frame UV and Luminosity Functions and Clustering of Bright Lyman Break Galaxies at z~3

The Redshift and Mass Dependence on the Formation of The Hubble Sequence at z>1 from CANDELS/UDS

The Slow Flow Model of Dust Efflux in Local Star-Forming Galaxies

Theoretical Evolution of Optical Strong Lines across Cosmic Time

A Census of Oxygen in Star-Forming Galaxies: An Empirical Model Linking Metallicities, Star Formation Rates and Outflows

Galaxy Gas Fractions at High-Redshift: The Tension between Observations and Cosmological Simulations

Luminous and High Stellar Mass Candidate Galaxies at z ~ 8 Discovered in CANDELS

Not Dead Yet: Cool Circumgalactic Gas in the Halos of Early Type Galaxies

Testing Subhalo Abundance Matching in Cosmological Smoothed Particle Hydrodynamics Simulations

Probing Population III Stars in Galaxy IOK-1 at z = 6.96 through He II Emission

The Large, Oxygen-Rich Halos of Star-Forming Galaxies Are A Major Reservoir of Galactic Metals

Imaging the molecular gas in a submm galaxy at z = 4.05: cold mode accretion or a major merger?

Galaxies in a simulated $Λ$CDM Universe I: cold mode and hot cores

Galaxies in a Simulated $Λ$CDM Universe II: Observable Properties and Constraints on Feedback

The Simulated HI Sky at low redshift

The Flux Auto- and Cross-Correlation of the Lyman-alpha Forest. II. Modelling Anisotropies with Cosmological Hydrodynamic Simulations

The Growth of Central and Satellite Galaxies in Cosmological Smoothed Particle Hydrodynamics Simulations

How Do Galaxies Get Their Gas?