Source author record

Freeke van de Voort

Freeke van de Voort appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

20 published item(s)

preprint2022arXiv

Formation and fate of low metallicity stars in IllustrisTNG50

Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in particular to early enrichment and to later accretion and mixing of lower metallicity gas. Because low metallicity stars are difficult to observe directly, cosmological simulations are crucial for understanding their formation. Here we quantify the rates and locations of low metallicity star formation using the high-resolution TNG50 magnetohydrodynamical cosmological simulation, and we examine where low metallicity stars end up at $z=0$. We find that $20\%$ of stars with $Z_*<0.1\,\mathrm{Z_\odot}$ form after $z=2$, and that such stars are still forming in galaxies of all masses at $z=0$ today. Moreover, most low-metallicity stars at $z=0$ reside in massive galaxies. We analyse the radial distribution of low metallicity star formation, and discuss the curious case of seven galaxies in TNG50 that form stars from primordial gas even at $z=0$.

preprint2022arXiv

Quenching in Cosmic Sheets: Tracing the Impact of Large Scale Structure Collapse on the Evolution of Dwarf Galaxies

Dwarf galaxies are thought to quench primarily due to environmental processes most typically occurring in galaxy groups and clusters or around single, massive galaxies. However, at earlier epochs, ($5 < z < 2$), the collapse of large scale structure (forming Zel'dovich sheets and subsequently filaments of the cosmic web) can produce volume-filling accretion shocks which elevate large swaths of the intergalactic medium (IGM) in these structures to a hot ($T>10^6$ K) phase. We study the impact of such an event on the evolution of central dwarf galaxies ($5.5 < \log M_* < 8.5$) in the field using a spatially large, high resolution cosmological zoom simulation which covers the cosmic web environment between two protoclusters. We find that the shock-heated sheet acts as an environmental quencher much like clusters and filaments at lower redshift, creating a population of quenched, central dwarf galaxies. Even massive dwarfs which do not quench are affected by the shock, with reductions to their sSFR and gas accretion. This process can potentially explain the presence of isolated quenched dwarf galaxies, and represents an avenue of pre-processing, via which quenched satellites of bound systems quench before infall.

preprint2022arXiv

The impact of natal kicks on galactic r-process enrichment by neutron star mergers

We study galactic enrichment with rapid neutron capture (r-process) elements in cosmological, magnetohydrodynamical simulations of a Milky Way-mass galaxy. We include a variety of enrichment models, based on either neutron star mergers or a rare class of core-collapse supernova as sole r-process sources. For the first time in cosmological simulations, we implement neutron star natal kicks on-the-fly to study their impact. With kicks, neutron star mergers are more likely to occur outside the galaxy disc, but how far the binaries travel before merging also depends on the kick velocity distribution and shape of the delay time distribution for neutron star mergers. In our fiducial model, the median r-process abundance ratio is somewhat lower and the trend with metallicity is slightly steeper when kicks are included. In a model 'optimized' to better match observations, with a higher rate of early neutron star mergers, the median r-process abundances are fairly unaffected by kicks. In both models, the scatter in r-process abundances is much larger with natal kicks, especially at low metallicity, giving rise to more r-process enhanced stars. We experimented with a range of kick velocities and find that with lower velocities, the scatter is reduced, but still larger than without natal kicks. We discuss the possibility that the observed scatter in r-process abundances is predominantly caused by natal kicks removing the r-process sources far from their birth sites, making enrichment more inhomogeneous, rather than the usual interpretation that the scatter is set by the rarity of its production source.

preprint2020arXiv

AlFoCS + Fornax3D: resolved star formation in the Fornax cluster with ALMA and MUSE

We combine data from ALMA and MUSE to study the resolved (~300 pc scale) star formation relation (star formation rate vs. molecular gas surface density) in cluster galaxies. Our sample consists of 9 Fornax cluster galaxies, including spirals, ellipticals, and dwarfs, covering a stellar mass range of ~10^8.8 - 10^11 M_Sun. CO(1-0) and extinction corrected Halpha were used as tracers for the molecular gas mass and star formation rate, respectively. We compare our results with Kennicutt (1998) and Bigiel et al. (2008). Furthermore, we create depletion time maps to reveal small-scale variations in individual galaxies. We explore these further in FCC290, using the 'uncertainty principle for star formation' (Kruijssen & Longmore, 2014a) to estimate molecular cloud lifetimes, which we find to be short (<10 Myr) in this galaxy. Galaxy-averaged depletion times are compared with other parameters such as stellar mass and cluster-centric distance. We find that the star formation relation in the Fornax cluster is close to those from Kennicutt (1998) and Bigiel et al. (2008}), but overlaps mostly with the shortest depletion times predicted by Bigiel et al. (2008). This slight decrease in depletion time is mostly driven by dwarf galaxies with disturbed molecular gas reservoirs close to the virial radius. In FCC90, a dwarf galaxy with a molecular gas tail, we find that depletion times are a factor >~10 higher in its tail than in its stellar body.

preprint2020arXiv

Magnetising the circumgalactic medium of disk galaxies

The circumgalactic medium (CGM) is one of the frontiers of galaxy formation and intimately connected to the galaxy via accretion of gas on to the galaxy and gaseous outflows from the galaxy. Here we analyse the magnetic field in the CGM of the Milky Way-like galaxies simulated as part of the \textsc{Auriga} project that constitutes a set of high resolution cosmological magnetohydrodynamical zoom simulations. We show that before $z=1$ the CGM becomes magnetised via galactic outflows that transport magnetised gas from the disk into the halo. At this time the magnetisation of the CGM closely follows its metal enrichment. We then show that at low redshift an in-situ turbulent dynamo that operates on a timescale of Gigayears further amplifies the magnetic field in the CGM and saturates before $z=0$. The magnetic field strength reaches a typical value of $0.1\,μG$ at the virial radius at $z=0$ and becomes mostly uniform within the virial radius. Its Faraday rotation signal is in excellent agreement with recent observations. For most of its evolution the magnetic field in the CGM is an unordered small scale field. Only strong coherent outflows at low redshift are able to order the magnetic field in parts of the CGM that are directly displaced by these outflows.

preprint2020arXiv

Neutron star mergers and rare core-collapse supernovae as sources of r-process enrichment in simulated galaxies

We use cosmological, magnetohydrodynamical simulations of Milky Way-mass galaxies from the Auriga project to study their enrichment with rapid neutron capture (r-process) elements. We implement a variety of enrichment models from both binary neutron star mergers and rare core-collapse supernovae. We focus on the abundances of (extremely) metal-poor stars, most of which were formed during the first ~Gyr of the Universe in external galaxies and later accreted onto the main galaxy. We find that the majority of metal-poor stars are r-process enriched in all our enrichment models. Neutron star merger models result in a median r-process abundance ratio which increases with metallicity, whereas the median trend in rare core-collapse supernova models is approximately flat. The scatter in r-process abundance increases for models with longer delay times or lower rates of r-process producing events. Our results are nearly perfectly converged, in part due to the mixing of gas between mesh cells in the simulations. Additionally, different Milky Way-mass galaxies show only small variation in their respective r-process abundance ratios. Current (sparse and potentially biased) observations of metal-poor stars in the Milky Way seem to prefer rare core-collapse supernovae over neutron star mergers as the dominant source of r-process elements at low metallicity, but we discuss possible caveats to our models. Dwarf galaxies which experience a single r-process event early in their history show highly enhanced r-process abundances at low metallicity, which is seen both in observations and in our simulations. We also find that the elements produced in a single event are mixed with ~10^8 Msun of gas relatively quickly, distributing the r-process elements over a large region.

preprint2020arXiv

The effect of magnetic fields on properties of the circumgalactic medium

We study the effect of magnetic fields on a simulated galaxy and its surrounding gaseous halo, or circumgalactic medium (CGM), within cosmological 'zoom-in' simulations of a Milky Way-mass galaxy as part of the 'Simulating the Universe with Refined Galaxy Environments' (SURGE) project. We use three different galaxy formation models, each with and without magnetic fields, and include additional spatial refinement in the CGM to improve its resolution. The central galaxy's star formation rate and stellar mass are not strongly affected by the presence of magnetic fields, but the galaxy is more disc-dominated and its central black hole is more massive when $B>0$. The physical properties of the CGM change significantly. With magnetic fields, the circumgalactic gas flows are slower, the atomic hydrogen-dominated extended discs around the galaxy are more massive and the densities in the inner CGM are therefore higher, the temperatures in the outer CGM are higher, and the pressure in the halo is higher and smoother. The total gas fraction and metal mass fraction in the halo are also higher when magnetic fields are included, because less gas escapes the halo. Additionally, we find that the CGM properties depend on azimuthal angle and that magnetic fields reduce the scatter in radial velocity, whilst enhancing the scatter in metallicity at fixed azimuthal angle. The metals are thus less well-mixed throughout the halo, resulting in more metal-poor halo gas. These results together show that magnetic fields in the CGM change the flow of gas in galaxy haloes, making it more difficult for metal-rich outflows to mix with the metal-poor CGM and to escape the halo, and therefore should be included in simulations of galaxy formation.

preprint2016arXiv

A massive, quiescent, population II galaxy at a redshift of 2.1

Unlike spiral galaxies such as the Milky Way, the majority of the stars in massive elliptical galaxies were formed in a short period early in the history of the Universe. The duration of this formation period can be measured using the ratio of magnesium to iron abundance ([Mg/Fe]), which reflects the relative enrichment by core-collapse and type Ia supernovae. For local galaxies, [Mg/Fe] probes the combined formation history of all stars currently in the galaxy, including younger and metal-poor stars that were added during late-time mergers. Therefore, to directly constrain the initial star-formation period, we must study galaxies at earlier epochs. The most distant galaxy for which [Mg/Fe] had previously been measured is at z~1.4, with [Mg/Fe]=0.45(+0.05,-0.19). A slightly earlier epoch (z~1.6) was probed by stacking the spectra of 24 massive quiescent galaxies, yielding an average [Mg/Fe] of 0.31+/-0.12. However, the relatively low S/N of the data and the use of index analysis techniques for both studies resulted in measurement errors that are too large to allow us to form strong conclusions. Deeper spectra at even earlier epochs in combination with analysis techniques based on full spectral fitting are required to precisely measure the abundance pattern shortly after the major star-forming phase (z>2). Here we report a measurement of [Mg/Fe] for a massive quiescent galaxy at z=2.1. With [Mg/Fe]=0.59+/-0.11, this galaxy is the most Mg-enhanced massive galaxy found so far, having twice the Mg enhancement of similar-mass galaxies today. The abundance pattern of the galaxy is consistent with enrichment exclusively by core-collapse supernovae and with a star-formation timescale of 0.1-0.5 Gyr - characteristics that are similar to population II stars in the Milky Way. With an average past SFR of 600-3000 Msol/yr, this galaxy was among the most vigorous star-forming galaxies in the Universe.

preprint2016arXiv

Cold-Mode Accretion: Driving the Fundamental Mass-Metallicity Relation at z~2

We investigate the star formation rate (SFR) dependence on the stellar mass and gas-phase metallicity relation at z=2 with MOSFIRE/Keck as part of the ZFIRE survey. We have identified 117 galaxies (1.98 < z < 2.56), with $8.9\leq$log(M/M$_{\odot}$)$\leq11.0$, for which we can measure gas-phase metallicities. For the first time, we show discernible difference between the mass-metallicity relation, using individual galaxies, when deviding the sample by low ($<10$~M$_{\odot}$yr$^{-1}$) and high ($>10$~M$_{\odot}$yr$^{-1}$) SFRs. At fixed mass, low star-forming galaxies tend to have higher metallicity than high star-forming galaxies. Using a few basic assumptions, we further show that the gas masses and metallicities required to produce the fundamental mass--metallicity relation, and its intrinsic scatter, are consistent with cold-mode accretion predictions obtained from the OWLS hydrodynamical simulations. Our results from both simulations and observations are suggestive that cold-mode accretion is responsible for the fundamental mass-metallicity relation at $z=2$ and demonstrates the direct relationship between cosmological accretion and the fundamental properties of galaxies.

preprint2016arXiv

Galaxy growth from redshift 5 to 0 at fixed comoving number density

Studying the average properties of galaxies at a fixed comoving number density over a wide redshift range has become a popular observational method, because it may trace the evolution of galaxies statistically. We test this method by comparing the evolution of galaxies at fixed number density and by following individual galaxies through cosmic time (z=0-5) in cosmological, hydrodynamical simulations from OWLS. Comparing progenitors, descendants, and galaxies selected at fixed number density at each redshift, we find differences of up to a factor of three for galaxy and interstellar medium (ISM) masses. The difference is somewhat larger for black hole masses. The scatter in ISM mass increases significantly towards low redshift with all selection techniques. We use the fixed number density technique to study the assembly of dark matter, gas, stars, and black holes and the evolution in accretion and star formation rates. We find three different regimes for massive galaxies, consistent with observations: at high redshift the gas accretion rate dominates, at intermediate redshifts the star formation rate is the highest, and at low redshift galaxies grow mostly through mergers. Quiescent galaxies have much lower ISM masses (by definition) and much higher black hole masses, but the stellar and halo masses are fairly similar. Without active galactic nucleus (AGN) feedback, massive galaxies are dominated by star formation down to z=0 and most of their stellar mass growth occurs in the centre. With AGN feedback, stellar mass is only added to the outskirts of galaxies by mergers and they grow inside-out.

preprint2016arXiv

The impact of stellar feedback on hot gas in galaxy haloes: the Sunyaev-Zel'dovich effect and soft X-ray emission

The thermal Sunyaev-Zel'dovich (SZ) effect and soft X-ray emission are routinely observed around massive galaxies and in galaxy groups and clusters. We study these observational diagnostics of galaxy haloes for a suite of cosmological `zoom-in' simulations from the `Feedback In Realistic Environments' project, which spans a large range in halo mass 10^10-10^13 Msun). We explore the effect of stellar feedback on the hot gas observables. The properties of our simulated groups, such as baryon fractions, SZ flux, and X-ray luminosities (L_X), are broadly consistent with existing observations, even though feedback from active galactic nuclei is not included. We make predictions for future observations of lower-mass objects for both SZ and diffuse X-ray measurements, finding that they are not just scaled-down versions of massive galaxies, but more strongly affected by galactic winds driven by star formation. Low-mass haloes (<~10^11 Msun) retain a low fraction of their baryons, which results in a strong suppression of the SZ signal. Our simulations therefore predict a scaling with halo mass that is steeper than self-similar for haloes less massive than 10^13 Msun. For halo masses <~10^12 Msun, L_X is time-variable and correlated primarily with the star formation rate (SFR). For these objects, the diffuse X-ray emission is powered mostly by galactic winds and the gas dominating the X-ray emission is flowing out with radial velocities close to the halo's circular velocity. For halo masses >~10^13 Msun, on the other hand, L_X is much less variable and not correlated with the SFR, because the emission originates from the quasi-hydrostatic, virialized halo gas.

preprint2015arXiv

The creation and persistence of a misaligned gas disc in a simulated early-type galaxy

Massive early-type galaxies commonly have gas discs which are kinematically misaligned with the stellar component. These discs feel a torque from the stars and the angular momentum vectors are expected to align quickly. We present results on the evolution of a misaligned gas disc in a cosmological simulation of a massive early-type galaxy from the Feedback In Realistic Environments project. This galaxy experiences a merger which, together with a strong galactic wind, removes most of the original gas disc. The galaxy subsequently reforms a gas disc through accretion of cold gas, but it is initially 120 degrees misaligned with the stellar rotation axis. This misalignment persists for about 2 Gyr before the gas-star misalignment angle drops below 20 degrees. The time it takes for the gaseous and stellar components to align is much longer than previously thought, because the gas disc is accreting a significant amount of mass for about 1.5 Gyr after the merger, during which the angular momentum change induced by accreted gas dominates over that induced by stellar torques. Once the gas accretion rate has decreased sufficiently, the gas disc decouples from the surrounding halo gas and realigns with the stellar component in about 6 dynamical times. During the late evolution of the misaligned gas disc, the centre aligns faster than the outskirts, resulting in a warped disc. We discuss the observational consequences of the long survival of our misaligned gas disc and how our results can be used to calibrate merger rate estimates from observed gas misalignments.

preprint2014arXiv

Characterizing simulated galaxy stellar mass histories

Cosmological galaxy formation simulations can now produce rich and diverse ensembles of galaxy histories. These simulated galaxy histories, taken all together, provide an answer to the question "how do galaxies form?" for the models used to construct them. We characterize such galaxy history ensembles both to understand their properties and to identify points of comparison for histories within a given galaxy formation model or between different galaxy formation models and simulations. We focus primarily on stellar mass histories of galaxies with the same final stellar mass, for six final stellar mass values and for three different simulated galaxy formation models (a semi-analytic model built upon the dark matter Millennium simulation and two models from the hydrodynamical OverWhelmingly Large Simulations project). Using principal component analysis (PCA) to classify scatter around the average stellar mass history, we find that one fluctuation dominates for all sets of histories we consider, although its shape and contribution can vary between samples. We correlate the PCA characterization with several $z=0$ galaxy properties (to connect with survey observables) and also compare it to some other galaxy history properties. We then explore separating galaxy stellar mass histories into classes, using the largest PCA contribution, k-means clustering, and simple Gaussian mixture models. For three component models, these different methods often gave similar results. These history classification methods provide a succinct and often quick way to characterize changes in the full ensemble of histories of a simulated population as physical assumptions are varied, to compare histories of different simulated populations to each other, and to assess the relation of simulated histories to fixed time observations.

preprint2014arXiv

Galactic r-process enrichment by neutron star mergers in cosmological simulations of a Milky Way-mass galaxy

We quantify the stellar abundances of neutron-rich r-process nuclei in cosmological zoom-in simulations of a Milky Way-mass galaxy from the Feedback In Realistic Environments project. The galaxy is enriched with r-process elements by binary neutron star (NS) mergers and with iron and other metals by supernovae. These calculations include key hydrodynamic mixing processes not present in standard semi-analytic chemical evolution models, such as galactic winds and hydrodynamic flows associated with structure formation. We explore a range of models for the rate and delay time of NS mergers, intended to roughly bracket the wide range of models consistent with current observational constraints. We show that NS mergers can produce [r-process/Fe] abundance ratios and scatter that appear reasonably consistent with observational constraints. At low metallicity, [Fe/H]<-2, we predict there is a wide range of stellar r-process abundance ratios, with both supersolar and subsolar abundances. Low-metallicity stars or stars that are outliers in their r-process abundance ratios are, on average, formed at high redshift and located at large galactocentric radius. Because NS mergers are rare, our results are not fully converged with respect to resolution, particularly at low metallicity. However, the uncertain rate and delay time distribution of NS mergers introduces an uncertainty in the r-process abundances comparable to that due to finite numerical resolution. Overall, our results are consistent with NS mergers being the source of most of the r-process nuclei in the Universe.

preprint2013arXiv

Soft X-ray and ultra-violet metal-line emission from the gas around galaxies

(Abridged) A large fraction of the gas in galactic haloes has temperatures between 10^4.5 and 10^7 K. At these temperatures, cooling is dominated by metal-line emission if the metallicity Z>~0.1 Zsun. We explore the detectability of several lines using large cosmological, hydrodynamical simulations. We stack surface brightness maps centred on galaxies to calculate the expected mean surface brightness profiles for different halo masses. Assuming a detection limit of 10^-1 photon s^-1 cm^-2 sr^-1, proposed X-ray telescopes can detect O VIII emission from z=0.125 out to 80% of the virial radius (Rvir) of groups and clusters and out to 0.4Rvir for haloes with masses Mhalo=10^12-13 Msun. Emission lines from C VI, N VII, O VII, and Ne X can be detected out to smaller radii, 0.1-0.5Rvir. With a detection limit of 10^-20 erg s^-1 cm^-2 arcsec^-2, future UV telescopes can detect C III emission out to 0.2-0.6Rvir at z=0.25. C IV, O VI, Si III, and Si IV can be seen out to 0.1-0.2Rvir for Mhalo>10^12 Msun. Optical HI H-alpha emission is comparable in strength to C III emission. At z=3 it may be possible to observe C III out to 0.2-0.3Rvir and other rest-frame UV lines out to ~0.1Rvir for Mhalo>10^11 Msun with upcoming optical instruments. Metal-line emission is typically biased towards high density and metallicity and towards the temperature at which the emissivity curve of the corresponding metal line peaks. The bias is similar for the different soft X-ray lines considered, whereas it varies strongly between different UV lines. Active galactic nucleus (AGN) feedback can change the inner surface brightness profiles significantly, but it generally does not change the radius out to which the emission can be observed. Metal-line emission is a promising probe of the warm and hot, enriched gas around galaxies and provides a unique window into the interactions between galaxies and their gaseous haloes.

preprint2012arXiv

Cold accretion flows and the nature of high column density HI absorption at redshift 3

Simulations predict that galaxies grow primarily through the accretion of gas that has not gone through an accretion shock near the virial radius and that this cold gas flows towards the central galaxy along dense filaments and streams. There is, however, little observational evidence for the existence of these cold flows. We use a large, cosmological, hydrodynamical simulation that has been post-processed with radiative transfer to study the contribution of cold flows to the observed z=3 column density distribution of neutral hydrogen, which our simulation reproduces. We find that nearly all of the HI absorption arises in gas that has remained colder than 10^5.5 K, at least while it was extragalactic. In addition, the majority of the HI is rapidly falling towards a nearby galaxy, with non-negligible contributions from outflowing and static gas. Above a column density of N_HI = 10^17 cm^-2, most of the absorbers reside inside haloes, but the interstellar medium only dominates for N_HI > 10^21 cm^-2. Haloes with total mass below 10^10 Msun dominate the absorption for 10^17<N_HI < 10^21 cm^-2, but the average halo mass increases sharply for higher column densities. Although very little of the HI in absorbers with N_HI <~ 10^20 cm^-2 resides inside galaxies, systems with N_HI > 10^17 cm^-2 are closely related to star formation: most of their HI either will become part of the interstellar medium before z=2 or has been ejected from a galaxy at z>3. Cold accretion flows are critical for the success of our simulation in reproducing the observed rate of incidence of damped Lyman-alpha and particularly that of Lyman limit systems. We therefore conclude that cold accretion flows exist and have already been detected in the form of high column density HI absorbers.

preprint2012arXiv

Properties of gas in and around galaxy haloes

We study the properties of gas inside and around galaxy haloes as a function of radius and halo mass, focussing mostly on z=2, but also showing some results for z=0. For this purpose, we use a suite of large cosmological, hydrodynamical simulations from the OverWhelmingly Large Simulations project. The properties of cold- and hot-mode gas, which we separate depending on whether the temperature has been higher than 10^5.5 K while it was extragalactic, are clearly distinguishable in the outer parts of massive haloes (virial temperatures >> 10^5 K. The differences between cold- and hot-mode gas resemble those between inflowing and outflowing gas. The cold-mode gas is mostly confined to clumpy filaments that are approximately in pressure equilibrium with the diffuse, hot-mode gas. Besides being colder and denser, cold-mode gas typically has a much lower metallicity and is much more likely to be infalling. However, the spread in the properties of the gas is large, even for a given mode and a fixed radius and halo mass, which makes it impossible to make strong statements about individual gas clouds. Metal-line cooling causes a strong cooling flow near the central galaxy, which makes it hard to distinguish gas accreted through the cold and hot modes in the inner halo. Stronger feedback results in larger outflow velocities and pushes hot-mode gas to larger radii. The gas properties evolve as expected from virial arguments, which can also account for the dependence of many gas properties on halo mass. We argue that cold streams penetrating hot haloes are observable as high-column density HI Lyman-alpha absorption systems in sightlines near massive foreground galaxies.

preprint2011arXiv

The drop in the cosmic star formation rate below redshift 2 is caused by a change in the mode of gas accretion and by AGN feedback

The cosmic star formation rate is observed to drop sharply after redshift z=2. We use a large, cosmological, smoothed particle hydrodynamics simulation to investigate how this decline is related to the evolution of gas accretion and to outflows driven by active galactic nuclei (AGN). We find that the drop in the star formation rate follows a corresponding decline in the global cold-mode accretion rate density onto haloes, but with a delay of order the gas consumption time scale in the interstellar medium. Here we define cold-mode (hot-mode) accretion as gas that is accreted and whose temperature has never exceeded (did exceed) 10^5.5 K. In contrast to cold-mode accretion, which peaks at z~3, the hot mode continues to increase to z~1 and remains roughly constant thereafter. By the present time, the hot mode strongly dominates the global accretion rate onto haloes. Star formation does not track hot-mode halo accretion because most of the hot halo gas never accretes onto galaxies. AGN feedback plays a crucial role by preferentially preventing gas that entered haloes in the hot mode from accreting onto their central galaxies. Consequently, in the absence of AGN feedback, gas accreted in the hot mode would become the dominant source of fuel for star formation and the drop off in the cosmic star formation rate would be much less steep.

preprint2011arXiv

The rates and modes of gas accretion on to galaxies and their gaseous haloes

(Abridged) We study the rate at which gas accretes onto galaxies and haloes and investigate whether the accreted gas was shocked to high temperatures before reaching a galaxy. For this purpose we use a suite of large cosmological, hydrodynamical simulations from the OWLS project. We improve on previous work by considering a wider range of halo masses and redshifts, by distinguishing accretion onto haloes and galaxies, by including important feedback processes, and by comparing simulations with different physics. The specific rate of gas accretion onto haloes is, like that for dark matter, only weakly dependent on halo mass. For halo masses Mhalo>>10^11 Msun it is relatively insensitive to feedback processes. In contrast, accretion rates onto galaxies are determined by radiative cooling and by outflows driven by supernovae and active galactic nuclei. Galactic winds increase the halo mass at which the central galaxies grow the fastest by about two orders of magnitude to Mhalo~10^12 Msun. Gas accretion is bimodal, with maximum past temperatures either of order the virial temperature or <~10^5 K. The fraction of gas accreted on to haloes in the hot mode is insensitive to feedback and metal-line cooling. It increases with decreasing redshift, but is mostly determined by halo mass, increasing gradually from less than 10% for ~10^11 Msun to greater than 90% at 10^13 Msun. In contrast, for accretion onto galaxies the cold mode is always significant and the relative contributions of the two accretion modes are more sensitive to feedback and metal-line cooling. The majority of stars present in any mass halo at any redshift were formed from gas accreted in the cold mode, although the hot mode contributes typically over 10% for Mhalo>~10^11 Msun. Galaxies, but not necessarily their gaseous haloes, are predominantly fed by gas that did not experience an accretion shock when it entered the host halo.

preprint2009arXiv

The physics driving the cosmic star formation history

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.