-
Simulating the diversity of shapes of the Lyman-$α$ line
Authors:
Jeremy Blaizot,
Thibault Garel,
Anne Verhamme,
Harley Katz,
Taysun Kimm,
Léo Michel-Dansac,
Peter D. Mitchell,
Joakim Rosdahl,
Maxime Trebitsch
Abstract:
The Ly$α$ line is a powerful probe of distant galaxies, which contains information about inflowing/outflowing gas through which Ly$α$ photons scatter. To develop our understanding of this probe, we post-process a zoom-in radiation-hydrodynamics simulation of a low-mass ($M_* \sim 10^9 M_\odot$) galaxy to construct 22500 mock spectra in 300 directions from $z = 3$ to 4. Remarkably, we show that one…
▽ More
The Ly$α$ line is a powerful probe of distant galaxies, which contains information about inflowing/outflowing gas through which Ly$α$ photons scatter. To develop our understanding of this probe, we post-process a zoom-in radiation-hydrodynamics simulation of a low-mass ($M_* \sim 10^9 M_\odot$) galaxy to construct 22500 mock spectra in 300 directions from $z = 3$ to 4. Remarkably, we show that one galaxy can reproduce the variety of a large sample of spectroscopically observed Ly$α$ line profiles. While most mock spectra exhibit double-peak profiles with a dominant red peak, their shapes cover a large parameter space in terms of peak velocities, peak separation and flux ratio. This diversity originates from radiative transfer effects at ISM and CGM scales, and depends on galaxy inclination and evolutionary phase. Red-dominated lines preferentially arise in face-on directions during post-starburst outflows and are bright. Conversely, accretion phases usually yield symmetric double peaks in the edge-on direction and are fainter. While resonant scattering effects at $< 0.2\times R_{\rm vir}$ are responsible for the broadening and velocity shift of the red peak, the extended CGM acts as a screen and impacts the observed peak separation. The ability of simulations to reproduce observed Ly$α$ profiles and link their properties with galaxy physical parameters offers new perspectives to use Ly$α$ to constrain the mechanisms that regulate galaxy formation and evolution. Notably, our study implies that deeper Ly$α$ surveys may unveil a new population of blue-dominated lines tracing inflowing gas.
△ Less
Submitted 17 May, 2023;
originally announced May 2023.
-
Baryonic mass budgets for haloes in the EAGLE simulation, including ejected and prevented gas
Authors:
Peter D. Mitchell,
Joop Schaye
Abstract:
Feedback processes are expected to shape galaxy evolution by ejecting gas from galaxies and their associated dark matter haloes, and also by preventing diffuse gas from ever being accreted. We present predictions from the EAGLE simulation project for the mass budgets associated with "ejected" and "prevented" gas, as well as for ejected metals. We find that most of the baryons that are associated w…
▽ More
Feedback processes are expected to shape galaxy evolution by ejecting gas from galaxies and their associated dark matter haloes, and also by preventing diffuse gas from ever being accreted. We present predictions from the EAGLE simulation project for the mass budgets associated with "ejected" and "prevented" gas, as well as for ejected metals. We find that most of the baryons that are associated with haloes of mass $10^{11} < M_{200} \, /\mathrm{M_\odot} < 10^{13}$ at $z=0$ have been ejected beyond the virial radius after having been accreted. When the gas ejected from satellites (and their progenitors) is accounted for, the combined ejected mass represents half of the total baryon budget even in the most massive simulated galaxy clusters ($M_{200} \approx 10^{14.5} \, \mathrm{M_\odot}$), with the consequence that the total baryon budget exceeds the cosmic average if ejected gas is included. We find that gas is only prevented from being accreted onto haloes for $M_{200} < 10^{12} \, \mathrm{M_\odot}$, and that this component accounts for about half the total baryon budget for $M_{200} < 10^{11} \, \mathrm{M_\odot}$, with ejected gas making up most of the remaining half. For metals, most of the mass that is not locked into stars has been ejected beyond the virial radius, at least for $M_{200} < 10^{13} \, \mathrm{M_\odot}$. Finally, within the virial radius we find that most of the mass in the circum-galactic medium (CGM) has not passed through the ISM of a progenitor galaxy, for all halo masses and redshifts. About half of the CGM within half the virial radius has passed through the ISM in the past, however.
△ Less
Submitted 15 December, 2021;
originally announced December 2021.
-
How gas flows shape the stellar-halo mass relation in the EAGLE simulation
Authors:
Peter D. Mitchell,
Joop Schaye
Abstract:
The difference in shape between the observed galaxy stellar mass function and the predicted dark matter halo mass function is generally explained primarily by feedback processes. Feedback can shape the stellar-halo mass (SHM) relation by driving gas out of galaxies, by modulating the first-time infall of gas onto galaxies (i.e., preventative feedback), and by instigating fountain flows of recycled…
▽ More
The difference in shape between the observed galaxy stellar mass function and the predicted dark matter halo mass function is generally explained primarily by feedback processes. Feedback can shape the stellar-halo mass (SHM) relation by driving gas out of galaxies, by modulating the first-time infall of gas onto galaxies (i.e., preventative feedback), and by instigating fountain flows of recycled wind material. We present and apply a method to disentangle these effects for hydrodynamical simulations of galaxy formation. We build a model of linear coupled differential equations that by construction reproduces the flows of gas onto and out of galaxies and haloes in the EAGLE cosmological simulation. By varying individual terms in this model, we isolate the relative effects of star formation, ejection via outflow, first-time inflow and wind recycling on the SHM relation. We find that for halo masses $M_{200} < 10^{12} \, \mathrm{M_\odot}$ the SHM relation is shaped primarily by a combination of ejection from galaxies and haloes, while for larger $M_{200}$ preventative feedback is also important. The effects of recycling and the efficiency of star formation are small. We show that if, instead of $M_{200}$, we use the cumulative mass of dark matter that fell in for the first time, the evolution of the SHM relation nearly vanishes. This suggests that the evolution is due to the definition of halo mass rather than to an evolving physical efficiency of galaxy formation. Finally, we demonstrate that the mass in the circum-galactic medium is much more sensitive to gas flows, especially recycling, than is the case for stars and the interstellar medium.
△ Less
Submitted 15 December, 2021; v1 submitted 19 March, 2021;
originally announced March 2021.
-
Explaining the scatter in the galaxy mass-metallicity relation with gas flows
Authors:
Maria L. van Loon,
Peter D. Mitchell,
Joop Schaye
Abstract:
The physical origin of the scatter in the relation between galaxy stellar mass and the metallicity of the interstellar medium, i.e. the Mass-Metallicity Relation (MZR), reflects the relative importance of key processes in galaxy evolution. The \eagle cosmological hydrodynamical simulation is used to investigate the correlations between the residuals of the MZR and the residuals of the relations be…
▽ More
The physical origin of the scatter in the relation between galaxy stellar mass and the metallicity of the interstellar medium, i.e. the Mass-Metallicity Relation (MZR), reflects the relative importance of key processes in galaxy evolution. The \eagle cosmological hydrodynamical simulation is used to investigate the correlations between the residuals of the MZR and the residuals of the relations between stellar mass and, respectively, specific inflow, outflow and star formation rate as well as the gas fraction for central galaxies. At low redshift, all these residuals are found to be anti-correlated with the residuals of the MZR for $M_\star/\mathrm{M}_\odot \lesssim 10^{10}$. The correlations between the residuals of the MZR and the residuals of the other relations with mass are interrelated, but we find that gas fraction, specific inflow rate and specific outflow rate all have at least some independent influence on the scatter of the MZR. We find that, while for $M_\star/\mathrm{M}_\odot > 10^{10.4}$ the specific mass of the nuclear black hole is most important, for $M_\star/\mathrm{M}_\odot \lesssim 10^{10.3}$ gas fraction and specific inflow rate are the variables that correlate most strongly with the MZR scatter. The timescales involved in the residual correlations and the time that galaxies stay above the MZR are revealed to be a few Gyr. However, most galaxies that are below the MZR at $z=0$ have been below the MZR throughout their lifetimes.
△ Less
Submitted 29 April, 2021; v1 submitted 26 January, 2021;
originally announced January 2021.
-
The impact of stellar and AGN feedback on halo-scale baryonic and dark matter accretion in the EAGLE simulations
Authors:
Ruby J. Wright,
Claudia del P. Lagos,
Chris Power,
Peter D. Mitchell
Abstract:
We use the EAGLE suite of hydrodynamical simulations to analyse accretion rates (and the breakdown of their constituent channels) onto haloes over cosmic time, comparing the behaviour of baryons and dark matter (DM). We also investigate the influence of sub-grid baryon physics on halo-scale inflow, specifically the consequences of modelling radiative cooling, as well as feedback from stars and act…
▽ More
We use the EAGLE suite of hydrodynamical simulations to analyse accretion rates (and the breakdown of their constituent channels) onto haloes over cosmic time, comparing the behaviour of baryons and dark matter (DM). We also investigate the influence of sub-grid baryon physics on halo-scale inflow, specifically the consequences of modelling radiative cooling, as well as feedback from stars and active galactic nuclei (AGN). We find that variations in halo baryon fractions at fixed mass (particularly their circum-galactic medium gas content) are very well correlated with variations in the baryon fraction of accreting matter, which we show to be heavily suppressed by stellar feedback in low-mass haloes, $M_{\rm halo}\lesssim10^{11.5}M_{\odot}$. Breaking down accretion rates into first infall, recycled, transfer and merger components, we show that baryons are much more likely to be smoothly accreted than to have originated from mergers when compared to DM, finding (averaged across halo mass) a merger contribution of $\approx6\%$ for baryons, and $\approx15\%$ for DM at $z\approx0$. We also show that the breakdown of inflow into different channels is strongly dependent on sub-grid physics, particularly the contribution of recycled accretion (accreting matter that has been previously ejected from progenitor haloes). Our findings highlight the dual role that baryonic feedback plays in regulating the evolution of galaxies and haloes: by (i) directly removing gas from haloes, and (ii) suppressing gas inflow to haloes.
△ Less
Submitted 5 August, 2020; v1 submitted 1 June, 2020;
originally announced June 2020.
-
Galactic inflow and wind recycling rates in the EAGLE simulations
Authors:
Peter D. Mitchell,
Joop Schaye,
Richard G. Bower
Abstract:
The role of galactic wind recycling represents one of the largest unknowns in galaxy evolution, as any contribution of recycling to galaxy growth is largely degenerate with the inflow rates of first-time infalling material, and the rates with which outflowing gas and metals are driven from galaxies. We present measurements of the efficiency of wind recycling from the EAGLE cosmological simulation…
▽ More
The role of galactic wind recycling represents one of the largest unknowns in galaxy evolution, as any contribution of recycling to galaxy growth is largely degenerate with the inflow rates of first-time infalling material, and the rates with which outflowing gas and metals are driven from galaxies. We present measurements of the efficiency of wind recycling from the EAGLE cosmological simulation project, leveraging the statistical power of large-volume simulations that reproduce a realistic galaxy population. We study wind recycling at the halo scale, i.e. gas that has been ejected beyond the halo virial radius, and at the galaxy scale, i.e. gas that has been ejected from the ISM to at least $\approx 10 \, \%$ of the virial radius (thus excluding smaller-scale galactic fountains). Galaxy-scale wind recycling is generally inefficient, with a characteristic return timescale that is comparable or longer than a Hubble time, and with an efficiency that clearly peaks at the characteristic halo mass of $M_{200} = 10^{12} \, \mathrm{M_\odot}$. Correspondingly, the majority of gas being accreted onto galaxies in EAGLE is infalling for the first time. At the halo scale, the efficiency of recycling onto haloes differs by orders of magnitude from values assumed by semi-analytic galaxy formation models. Differences in the efficiency of wind recycling with other hydrodynamical simulations are currently difficult to assess, but are likely smaller. We are able to show that the fractional contribution of wind recycling to galaxy growth is smaller in EAGLE than in some other simulations. We find that cumulative first-time gas accretion rates at the virial radius are reduced relative to the expectation from dark matter accretion for haloes with mass, $M_{200} < 10^{12} \, \mathrm{M_\odot}$, indicating efficient preventative feedback on halo scales.
△ Less
Submitted 7 August, 2020; v1 submitted 20 May, 2020;
originally announced May 2020.
-
Galactic outflow rates in the EAGLE simulations
Authors:
Peter D. Mitchell,
Joop Schaye,
Richard G. Bower,
Robert A. Crain
Abstract:
We present measurements of galactic outflow rates from the EAGLE suite of cosmological simulations. We find that gas is removed from the interstellar medium (ISM) of central galaxies with a dimensionless mass loading factor that scales approximately with circular velocity as $V_{\mathrm{c}}^{-3/2}$ in the low-mass regime where stellar feedback dominates. Feedback from active galactic nuclei (AGN)…
▽ More
We present measurements of galactic outflow rates from the EAGLE suite of cosmological simulations. We find that gas is removed from the interstellar medium (ISM) of central galaxies with a dimensionless mass loading factor that scales approximately with circular velocity as $V_{\mathrm{c}}^{-3/2}$ in the low-mass regime where stellar feedback dominates. Feedback from active galactic nuclei (AGN) causes an upturn in the mass loading for halo masses $> 10^{12} \, \mathrm{M_\odot}$. We find that more gas outflows through the halo virial radius than is removed from the ISM of galaxies, particularly at low redshift, implying substantial mass loading within the circum-galactic medium (CGM). Outflow velocities span a wide range at a given halo mass/redshift, and on average increase positively with redshift and halo mass up to $M_{200} \sim 10^{12} \, \mathrm{M_\odot}$. Outflows exhibit a bimodal flow pattern on circum-galactic scales, aligned with the galactic minor axis. We present a number of like-for-like comparisons to outflow rates from other recent cosmological hydrodynamical simulations, and show that comparing the propagation of galactic winds as a function of radius reveals substantial discrepancies between different models. Relative to some other simulations, EAGLE favours a scenario for stellar feedback where agreement with the galaxy stellar mass function is achieved by removing smaller amounts of gas from the ISM, but with galactic winds that then propagate and entrain ambient gas out to larger radii.
△ Less
Submitted 1 April, 2020; v1 submitted 21 October, 2019;
originally announced October 2019.
-
The effect of gas accretion on the radial gas metallicity profile of simulated galaxies
Authors:
Florencia Collacchioni,
Claudia D. P. Lagos,
Peter D. Mitchell,
Joop Schaye,
Emily Wisnioski,
Sofía A. Cora,
Camila A. Correa
Abstract:
We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where…
▽ More
We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where higher $\dot M_{\rm accr}$ are associated with more negative slopes. The slope of the RMPs depends more strongly on $\dot M_{\rm accr}$ than on stellar mass, star formation rate or gas fraction, suggesting $\dot M_{\rm accr}$ to be a more fundamental driver of the RMP slope of galaxies. We find that eliminating the dependence on stellar mass is essential for pinning down the properties that shape the slope of the RMP. Although $\dot M_{\rm accr}$ is the main property modulating the slope of the RMP, we find that it causes other correlations that are more easily testable observationally: at fixed stellar mass, galaxies with more negative RMP slopes tend to have higher gas fractions and SFRs, while galaxies with lower gas fractions and SFRs tend to have flatter metallicity profiles within an effective radius.
△ Less
Submitted 28 May, 2020; v1 submitted 11 October, 2019;
originally announced October 2019.
-
Fast Outflows Identified in Early Star-Forming Galaxies at $z = 5-6$
Authors:
Yuma Sugahara,
Masami Ouchi,
Yuichi Harikane,
Nicolas Bouché,
Peter D. Mitchell,
Jérémy Blaizot
Abstract:
We present velocities of galactic outflows in seven star-forming galaxies at $z=$5-6 with stellar masses of $M_* \sim10^{10.1}\ \rm{M_\odot}$. Although it is challenging to observationally determine the outflow velocities, we overcome this by using ALMA [CII]158 $μ$m emission lines for systemic velocities and deep Keck spectra with metal absorption lines for velocity profiles available to date. We…
▽ More
We present velocities of galactic outflows in seven star-forming galaxies at $z=$5-6 with stellar masses of $M_* \sim10^{10.1}\ \rm{M_\odot}$. Although it is challenging to observationally determine the outflow velocities, we overcome this by using ALMA [CII]158 $μ$m emission lines for systemic velocities and deep Keck spectra with metal absorption lines for velocity profiles available to date. We construct a composite Keck spectrum of the galaxies at $z=$5-6 with the [CII]-systemic velocities, and fit outflow-line profiles to the SiII1260, CII1335, and SiIV1394,1403 absorption lines in the composite spectrum. We measure the maximum (90%) and central outflow velocities to be $v_\rm{max}=700^{+180}_{-110}\ \rm{km\ s^{-1}}$ and $v_\rm{out}= 400^{+100}_{-150}\ \rm{km\ s^{-1}}$ on average, respectively, showing no significant differences between the outflow velocities derived with the low to high-ionization absorption lines. For $M_* \sim10^{10.1}\ \rm{M_\odot}$, we find that the $v_\rm{max}$ value of our $z=$5-6 galaxies is 3 times higher than those of $z\sim0$ galaxies and comparable to $z\sim2$ galaxies. Estimating the halo circular velocity $v_\rm{cir}$ from the stellar masses and the abundance matching results, we investigate a $v_\rm{max}$-$v_\rm{cir}$ relation. Interestingly, $v_\rm{max}$ for galaxies with $M_*=10^{10.0-10.8}\ \rm{M_\odot}$ shows a clear positive correlation with $v_\rm{cir}$ and/or the galaxy star formation rate over $z=$0-6 with a small scatter of $\simeq \pm 0.1$ dex, which is in good agreement with theoretical predictions. This positive correlation suggests that the outflow velocity is physically related to the halo circular velocity, and that the redshift evolution of $v_\rm{max}$ at fixed $M_*$ is explained by the increase in $v_\rm{cir}$ toward high redshift.
△ Less
Submitted 30 October, 2019; v1 submitted 5 April, 2019;
originally announced April 2019.
-
Shark: introducing an open source, free and flexible semi-analytic model of galaxy formation
Authors:
Claudia del P. Lagos,
Rodrigo J. Tobar,
Aaron S. G. Robotham,
Danail Obreschkow,
Peter D. Mitchell,
Chris Power,
Pascal J. Elahi
Abstract:
We present a new, open source, free semi-analytic model (SAM) of galaxy formation, Shark, designed to be highly flexible and modular, allowing easy exploration of different physical processes and ways of modelling them. We introduce the philosophy behind Shark and provide an overview of the physical processes included in the model. Shark is written in C++11 and has been parallelized with OpenMP. I…
▽ More
We present a new, open source, free semi-analytic model (SAM) of galaxy formation, Shark, designed to be highly flexible and modular, allowing easy exploration of different physical processes and ways of modelling them. We introduce the philosophy behind Shark and provide an overview of the physical processes included in the model. Shark is written in C++11 and has been parallelized with OpenMP. In the released version (v1.1), we implement several different models for gas cooling, active galactic nuclei, stellar and photo-ionisation feedback, and star formation (SF). We demonstrate the basic performance of Shark using the Planck15 cosmology SURFS simulations, by comparing against a large set of observations, including: the stellar mass function (SMF) and stellar-halo mass relation at z=0-4; the cosmic evolution of the star formation rate density (SFRD), stellar mass, atomic and molecular hydrogen; local gas scaling relations; and structural galaxy properties, finding excellent agreement. Significant improvements over previous SAMs are seen in the mass-size relation for disks/bulges, the gas-stellar mass and stellar mass-metallicity relations. To illustrate the power of Shark in exploring the systematic effects of the galaxy formation modelling, we quantify how the scatter of the SF main sequence and the gas scaling relations changes with the adopted SF law, and the effect of the starbursts H$_2$ depletion timescale on the SFRD and $Ω_{\rm H_2}$. We compare Shark with other SAMs and the hydrodynamical simulation EAGLE, and find that SAMs have a much higher halo baryon fractions due to large amounts of intra-halo gas, which in the case of EAGLE is in the intergalactic medium.
△ Less
Submitted 24 September, 2018; v1 submitted 30 July, 2018;
originally announced July 2018.
-
Comparing galaxy formation in semi-analytic models and hydrodynamical simulations
Authors:
Peter D. Mitchell,
Cedric G. Lacey,
Claudia D. P. Lagos,
Carlos S. Frenk,
Richard G. Bower,
Shaun Cole,
John C. Helly,
Matthieu Schaller,
Violeta Gonzalez-Perez,
Tom Theuns
Abstract:
It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We…
▽ More
It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.
△ Less
Submitted 27 September, 2017; v1 submitted 25 September, 2017;
originally announced September 2017.
-
A unified multi-wavelength model of galaxy formation
Authors:
Cedric G. Lacey,
Carlton M. Baugh,
Carlos S. Frenk,
Andrew J. Benson,
Richard G. Bower,
Shaun Cole,
Violeta Gonzalez-Perez,
John C. Helly,
Claudia D. P. Lagos,
Peter D. Mitchell
Abstract:
We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy di…
▽ More
We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged]
△ Less
Submitted 6 August, 2016; v1 submitted 28 September, 2015;
originally announced September 2015.
-
Constraining the properties of AGN host galaxies with Spectral Energy Distribution modeling
Authors:
L. Ciesla,
V. Charmandaris,
A. Georgakakis,
E. Bernhard,
P. D. Mitchell,
V. Buat,
D. Elbaz,
E. Le Floc'h,
C. G. Lacey,
G. E. Magdis,
M. Xilouris
Abstract:
[abridged] We use the latest release of CIGALE, a galaxy SED fitting model relying on energy balance, to study the influence of an AGN in estimating both the SFR and stellar mass in galaxies, as well as the contribution of the AGN to the power output of the host. Using the galaxy formation SAM GALFORM, we create mock galaxy SEDs using realistic star formation histories (SFH) and add an AGN of Type…
▽ More
[abridged] We use the latest release of CIGALE, a galaxy SED fitting model relying on energy balance, to study the influence of an AGN in estimating both the SFR and stellar mass in galaxies, as well as the contribution of the AGN to the power output of the host. Using the galaxy formation SAM GALFORM, we create mock galaxy SEDs using realistic star formation histories (SFH) and add an AGN of Type 1, Type 2, or intermediate type whose contribution to the bolometric luminosity can be variable. We perform an SED fitting of these catalogues with CIGALE assuming three different SFHs: a single- and double-exponentially-decreasing, and a delayed SFH. Constraining thecontribution of an AGN to the LIR (fracAGN) is very challenging for fracAGN<20%, with uncertainties of ~5-30% for higher fractions depending on the AGN type, while FIR and sub-mm are essential. The AGN power has an impact on the estimation of $M_*$ in Type 1 and intermediate type AGNs but has no effect for galaxies hosting Type 2 AGNs. We find that in the absence of AGN emission, the best estimates of $M_*$ are obtained using the double-exponentially-decreasing model but at the expense of realistic ages of the stellar population. The delayed SFH model provides good estimates of $M_*$ and SFR, with a maximum offset of 10% as well as better estimates of the age. Our analysis shows that the underestimation of the SFR increases with fracAGN for Type 1 systems, as well as for low contributions of an intermediate AGN type, but it is quite insensitive to the emission of Type 2 AGNs up to fracAGN~45%. Similarly the UV emission is critical in accurately retrieving the $M_*$ for Type 1 and intermediate type AGN, and the SFR of all of the three AGN types. We show that the presence of AGN emission introduces a scatter to the SFR-$M_*$ main sequence relation derived from SED fitting, which is driven by the uncertainties on $M_*$.
△ Less
Submitted 16 January, 2015; v1 submitted 15 January, 2015;
originally announced January 2015.
-
A new methodology to test galaxy formation models using the dependence of clustering on stellar mass
Authors:
David J. R. Campbell,
Carlton M. Baugh,
Peter D. Mitchell,
John C. Helly,
Violeta Gonzalez-Perez,
Cedric G. Lacey,
Claudia del P. Lagos,
Vimal Simha,
Daniel J. Farrow
Abstract:
We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different…
▽ More
We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger timescales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from GAMA, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from SDSS and VIPERS, depending on the GALFORM model used.
△ Less
Submitted 10 June, 2015; v1 submitted 11 December, 2014;
originally announced December 2014.
-
The evolution of the star forming sequence in hierarchical galaxy formation models
Authors:
Peter D. Mitchell,
Cedric G. Lacey,
Shaun Cole,
Carlton M. Baugh
Abstract:
It has been argued that the specific star formation rates of star forming galaxies inferred from observational data decline more rapidly below z = 2 than is predicted by hierarchical galaxy formation models. We present a detailed analysis of this problem by comparing predictions from the GALFORM semi-analytic model with an extensive compilation of data on the average star formation rates of star-f…
▽ More
It has been argued that the specific star formation rates of star forming galaxies inferred from observational data decline more rapidly below z = 2 than is predicted by hierarchical galaxy formation models. We present a detailed analysis of this problem by comparing predictions from the GALFORM semi-analytic model with an extensive compilation of data on the average star formation rates of star-forming galaxies. We also use this data to infer the form of the stellar mass assembly histories of star forming galaxies. Our analysis reveals that the currently available data favour a scenario where the stellar mass assembly histories of star forming galaxies rise at early times and then fall towards the present day. In contrast, our model predicts stellar mass assembly histories that are almost flat below z = 2 for star forming galaxies, such that the predicted star formation rates can be offset with respect to the observational data by factors of up to 2-3. This disagreement can be explained by the level of coevolution between stellar and halo mass assembly that exists in contemporary galaxy formation models. In turn, this arises because the standard implementations of star formation and supernova feedback used in the models result in the efficiencies of these process remaining approximately constant over the lifetime of a given star forming galaxy. We demonstrate how a modification to the timescale for gas ejected by feedback to be reincorporated into galaxy haloes can help to reconcile the model predictions with the data.
△ Less
Submitted 26 August, 2014; v1 submitted 6 March, 2014;
originally announced March 2014.
-
A new model of galaxy formation: How sensitive are predicted galaxy luminosities to the choice of SPS model?
Authors:
V. Gonzalez-Perez,
C. G. Lacey,
C. M. Baugh,
C. D. P. Lagos,
J. Helly,
D. J. R. Campbell,
P. D. Mitchell
Abstract:
We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the WMAP7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven diff…
▽ More
We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the WMAP7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven different SPS models. In each case we obtain the rest-frame luminosity function in the far-ultra-violet, optical and near-infrared (NIR) wavelength ranges. We find that both the predicted rest-frame ultra-violet and optical luminosity function are insensitive to the choice of SPS model. However, we find that the predicted evolution of the rest-frame NIR luminosity function depends strongly on the treatment of the thermally pulsating asymptotic giant branch (TP-AGB) stellar phase in the SPS models, with differences larger than a factor of 2 for model galaxies brighter than $M_{\rm AB}(K)-5$log$h<-22$ ($\sim$L$_*$ for $0\leq z\leq 1.5$). We have also explored the predicted number counts of galaxies, finding remarkable agreement between the results with different choices of SPS model, except when selecting galaxies with very red optical-NIR colours. The predicted number counts of these extremely red galaxies appear to be more affected by the treatment of star formation in disks than by the treatment of TP-AGB stars in the SPS models.
△ Less
Submitted 6 February, 2014; v1 submitted 26 September, 2013;
originally announced September 2013.
-
How well can we really estimate the stellar masses of galaxies from broadband photometry?
Authors:
P. D. Mitchell,
C. G. Lacey,
C. M. Baugh,
S. Cole
Abstract:
The estimated stellar masses of galaxies are widely used to characterize how the galaxy population evolves over cosmic time. If stellar masses can be estimated in a robust manner, free from any bias, global diagnostics such as the stellar mass function can be used to constrain the physics of galaxy formation. We explore how galaxy stellar masses, estimated by fitting broad-band spectral energy dis…
▽ More
The estimated stellar masses of galaxies are widely used to characterize how the galaxy population evolves over cosmic time. If stellar masses can be estimated in a robust manner, free from any bias, global diagnostics such as the stellar mass function can be used to constrain the physics of galaxy formation. We explore how galaxy stellar masses, estimated by fitting broad-band spectral energy distributions (SEDs) with stellar population models, can be biased as a result of commonly adopted assumptions for the star-formation and chemical enrichment histories, recycled fractions and dust attenuation curves of galaxies. We apply the observational technique of broad-band SED fitting to model galaxy SEDs calculated by the theoretical galaxy formation model GALFORM, isolating the effect of each of these assumptions. We find that, averaged over the entire galaxy population, the common assumption of exponentially declining star-formation histories does not adversely affect stellar mass estimation. We show that fixing the metallicity in SED fitting or using sparsely sampled metallicity grids can introduce mass dependent systematics into stellar mass estimates. We find that the common assumption of a star-dust geometry corresponding to a uniform foreground dust screen can cause the stellar masses of dusty model galaxies to be significantly underestimated. Finally, we show that stellar mass functions recovered by applying SED fitting to model galaxies at high redshift can differ significantly in both shape and normalization from the intrinsic mass functions predicted by a given model. Given these differences, our methodology of using stellar masses estimated from model galaxy SEDs offers a new, self-consistent way to compare model predictions with observations.
△ Less
Submitted 4 September, 2013; v1 submitted 27 March, 2013;
originally announced March 2013.