Showing 1–50 of 132 results for author: Thomas, P A

First Light and Reionization Epoch Simulations (FLARES) -- XV: The physical properties of super-massive black holes and their impact on galaxies in the early universe

Authors: Stephen M. Wilkins, Jussi K. Kuusisto, Dimitrios Irodotou, Shihong Liao, Christopher C. Lovell, Sonja Soininen, Sabrina C. Berger, Sophie L. Newman, William J. Roper, Louise T. C. Seeyave, Peter A. Thomas, Aswin P. Vijayan

Abstract: Understanding the co-evolution of super-massive black holes (SMBHs) and their host galaxies remains a key challenge of extragalactic astrophysics, particularly the earliest stages at high-redshift. However, studying SMBHs at high-redshift with cosmological simulations, is challenging due to the large volumes and high-resolution required. Through its innovative simulation strategy, the First Light… ▽ More Understanding the co-evolution of super-massive black holes (SMBHs) and their host galaxies remains a key challenge of extragalactic astrophysics, particularly the earliest stages at high-redshift. However, studying SMBHs at high-redshift with cosmological simulations, is challenging due to the large volumes and high-resolution required. Through its innovative simulation strategy, the First Light And Reionisation Epoch Simulations (FLARES) suite of cosmological hydrodynamical zoom simulations allows us to simulate a much wider range of environments which contain SMBHs with masses extending to $M_{\bullet}>10^{9}\ M_{\odot}$ at $z=5$. In this paper, we use FLARES to study the physical properties of SMBHs and their hosts in the early Universe ($5\le\, z \le10$). FLARES predicts a sharply declining density with increasing redshift, decreasing by a factor of 100 over the range $z=5\to 10$. Comparison between our predicted bolometric luminosity function and pre-\emph{JWST} observations yield a good match. However, recent \emph{JWST} observations appear to suggest a larger contribution of SMBHs than previously observed, or predicted by FLARES. Finally, by using a re-simulation with AGN feedback disabled, we explore the impact of AGN feedback on their host galaxies. This reveals that AGN feedback results in a reduction of star formation activity, even at $z>5$, but only in the most massive galaxies. A deeper analysis reveals that AGN are also the cause of suppressed star formation in passive galaxies but that the presence of an AGN doesn't necessarily result in the suppression of star formation. △ Less

Submitted 9 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 17 pages, 13 figures, to be submitted to the Open Journal of Astrophysics. The first two authors are joint primary authors. Updated with new observational constraints

The impact of binary stars on the dust and metal evolution of galaxies

Authors: Robert M. Yates, David Hendriks, Aswin P. Vijayan, Robert G. Izzard, Peter A. Thomas, Payel Das

Abstract: We present detailed implementations of (a) binary stellar evolution (using binary_c) and (b) dust production and destruction into the cosmological semi-analytic galaxy evolution simulation, L-Galaxies. This new version of L-Galaxies is compared to a version assuming only single stars and to global and spatially-resolved observational data across a range of redshifts ($z$). We find that binaries ha… ▽ More We present detailed implementations of (a) binary stellar evolution (using binary_c) and (b) dust production and destruction into the cosmological semi-analytic galaxy evolution simulation, L-Galaxies. This new version of L-Galaxies is compared to a version assuming only single stars and to global and spatially-resolved observational data across a range of redshifts ($z$). We find that binaries have a negligible impact on the stellar masses, gas masses, and star formation rates of galaxies only if the total mass ejected by massive stars is unchanged. This is because massive stars determine the strength of supernova (SN) feedback, which in turn regulates galaxy growth. Binary effects, such as common envelope ejection and novae, affect carbon and nitrogen enrichment in galaxies, however heavier alpha elements are more affected by the choice of SN and wind yields. Unlike many other simulations, the new L-Galaxies reproduces observed dust-to-metal (DTM) and dust-to-gas (DTG) ratios at $z\sim{}0-4$. This is mainly due to shorter dust accretion timescales in dust-rich environments. However, dust masses are under-predicted at $z>4$, highlighting the need for enhanced dust production at early times in simulations, possibly accompanied by increased star formation. On sub-galactic scales, there is very good agreement between L-Galaxies and observed dust and metal radial profiles at $z=0$. A drop in DTM ratio is also found in diffuse, low-metallicity regions, contradicting the assumption of a universal value. We hope that this work serves as a useful template for binary stellar evolution implementations in other cosmological simulations in future. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 19 pages, 10 figures. Submitted to MNRAS. Comments welcome!

Non-polaritonic effects in cavity-modified photochemistry

Authors: Philip A. Thomas, Wai Jue Tan, Vasyl G. Kravets, Alexander N. Grigorenko, William L. Barnes

Abstract: Strong coupling of molecules to vacuum fields has been widely reported to lead to modified chemical properties such as reaction rates. However, some recent attempts to reproduce infrared strong coupling results have not been successful, suggesting that factors other than strong coupling may sometimes be involved. Here we re-examine the first of these vacuum-modified chemistry experiments in which… ▽ More Strong coupling of molecules to vacuum fields has been widely reported to lead to modified chemical properties such as reaction rates. However, some recent attempts to reproduce infrared strong coupling results have not been successful, suggesting that factors other than strong coupling may sometimes be involved. Here we re-examine the first of these vacuum-modified chemistry experiments in which changes to a molecular photoisomerisation process in the UV-vis spectral range were attributed to strong coupling of the molecules to visible light. We observed significant variations in photoisomerisation rates consistent with the original work; however, we found no evidence that these changes need to be attributed to strong coupling. Instead, we suggest that the photoisomerisation rates involved are most strongly influenced by the absorption of ultraviolet radiation in the cavity. Our results indicate that care must be taken to rule out non-polaritonic effects before invoking strong coupling to explain any changes of chemical properties arising in cavity-based experiments. △ Less

Submitted 14 July, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 33 pages, 17 figures

Journal ref: Adv. Mater. 2023, 2309393

First Light And Reionisation Epoch Simulations (FLARES) XIII: The Lyman-continuum emission of high-redshift galaxies

Authors: Louise T. C. Seeyave, Stephen M. Wilkins, Jussi K. Kuusisto, Christopher C. Lovell, Dimitrios Irodotou, Charlotte Simmonds, Aswin P. Vijayan, Peter A. Thomas, William J. Roper, Conor M. Byrne, Gareth T. Jones, Jack C. Turner, Christopher J. Conselice

Abstract: The history of reionisation is highly dependent on the ionising properties of high-redshift galaxies. It is therefore important to have a solid understanding of how the ionising properties of galaxies are linked to physical and observable quantities. In this paper, we use the First Light and Reionisation Epoch Simulations (FLARES) to study the Lyman-continuum (LyC, i.e. hydrogen-ionising) emission… ▽ More The history of reionisation is highly dependent on the ionising properties of high-redshift galaxies. It is therefore important to have a solid understanding of how the ionising properties of galaxies are linked to physical and observable quantities. In this paper, we use the First Light and Reionisation Epoch Simulations (FLARES) to study the Lyman-continuum (LyC, i.e. hydrogen-ionising) emission of massive ($M_*>10^8\,\mathrm{M_\odot}$) galaxies at redshifts $z=5-10$. We find that the specific ionising emissivity (i.e. intrinsic ionising emissivity per unit stellar mass) decreases as stellar mass increases, due to the combined effects of increasing age and metallicity. FLARES predicts a median ionising photon production efficiency (i.e. intrinsic ionising emissivity per unit intrinsic far-UV luminosity) of $\log_{10}(ξ_{\rm ion}\rm{/erg^{-1}Hz})=25.40^{+0.16}_{-0.17}$, with values spanning the range $\log_{10}(ξ_{\rm ion}\rm{/erg^{-1}Hz})=25-25.75$. This is within the range of many observational estimates, but below some of the extremes observed. We compare the production efficiency with observable properties, and find a weak negative correlation with the UV-continuum slope, and a positive correlation with the OIII equivalent width. We also consider the dust-attenuated production efficiency (i.e. intrinsic ionising emissivity per unit dust-attenuated far-UV luminosity), and find a median of $\log_{10}(ξ_{\rm ion}\rm{/erg^{-1}Hz})\sim25.5$. Within our sample of $M_*>10^8\,\mathrm{M_\odot}$ galaxies, it is the stellar populations in low mass galaxies that contribute the most to the total ionising emissivity. Active galactic nuclei (AGN) emission accounts for $10-20$ % of the total emissivity at a given redshift, and extends the LyC luminosity function by $\sim0.5$ dex. △ Less

Submitted 21 August, 2023; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: 19 pages, 17 figures, accepted for publication in MNRAS

First Light And Reionisation Epoch Simulations (FLARES) XII: The consequences of star-dust geometry on galaxies in the EoR

Authors: Aswin P. Vijayan, Peter A. Thomas, Christopher C. Lovell, Stephen M. Wilkins, Thomas R. Greve, Dimitrios Irodotou, William J. Roper, Louise T. C. Seeyave

Abstract: Using the First Light And Reionisation Epoch Simulations (${\rm F{\small LARES}}$), a suite of hydrodynamical simulations we explore the consequences of a realistic model for star--dust geometry on the observed properties of galaxies. We find that the UV attenuation declines rapidly from the central regions of galaxies, and bright galaxies have spatially extended star formation that suffers less o… ▽ More Using the First Light And Reionisation Epoch Simulations (${\rm F{\small LARES}}$), a suite of hydrodynamical simulations we explore the consequences of a realistic model for star--dust geometry on the observed properties of galaxies. We find that the UV attenuation declines rapidly from the central regions of galaxies, and bright galaxies have spatially extended star formation that suffers less obscuration than their fainter counterparts, demonstrating a non-linear relationship between the UV luminosity and the UV attenuation, giving a double power-law shape to the UVLF. Spatially distinct stellar populations within galaxies experience a wide range of dust attenuation due to variations in the dust optical depth along their line-of-sight; which can range from completely dust obscured to being fully unobscured. The overall attenuation curve of a galaxy is then a complex combination of various lines-of-sight within the galaxy. We explore the manifestation of this effect to study the reliability of line ratios to infer galaxy properties, in particular the Balmer decrement and the BPT diagram. We find the Balmer decrement predicted Balmer line attenuation to be higher (factor of $1$ to $\gtrsim10$) than expected from commonly used attenuation curves. The observed BPT line ratios deviate from their intrinsic values (median difference of 0.08 (0.02) and standard deviation of 0.2 (0.05) for log$_{10}$([N${\rm {\small II}}$]$λ6585/$H$_α$) (log$_{10}$([O${\rm {\small III}}$]$λ5008/$H$_β$)). Finally, we explore the variation in observed properties (UV attenuation, UV slope and Balmer decrement) with viewing angle, finding average differences of $\sim0.3$ magnitudes in the UV attenuation. △ Less

Submitted 11 December, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: 17 pages, 13 figures. Published MNRAS. Data available at https://flaresimulations.github.io

First Light and Reionisation Epoch Simulations (FLARES) X: Environmental Galaxy Bias and Survey Variance at High Redshift

Authors: Peter A. Thomas, Christopher C. Lovell, Maxwell G. A. Maltz, Aswin P. Vijayan, Stephen M. Wilkins, Dimitrios Irodotou, William J. Roper, Louise Seeyave

Abstract: Upcoming deep galaxy surveys with JWST will probe galaxy evolution during the epoch of reionisation (EoR, $5\leq z\leq10$) over relatively compact areas (e.g. $\sim$ 300\,arcmin$^2$ for the JADES GTO survey). It is therefore imperative that we understand the degree of survey variance, to evaluate how representative the galaxy populations in these studies will be. We use the First Light And Reionis… ▽ More Upcoming deep galaxy surveys with JWST will probe galaxy evolution during the epoch of reionisation (EoR, $5\leq z\leq10$) over relatively compact areas (e.g. $\sim$ 300\,arcmin$^2$ for the JADES GTO survey). It is therefore imperative that we understand the degree of survey variance, to evaluate how representative the galaxy populations in these studies will be. We use the First Light And Reionisation Epoch Simulations (FLARES) to measure the galaxy bias of various tracers over an unprecedentedly large range in overdensity for a hydrodynamic simulation, and use these relations to assess the impact of bias and clustering on survey variance in the EoR. Star formation is highly biased relative to the underlying dark matter distribution, with the mean ratio of the stellar to dark matter density varying by a factor of 100 between regions of low and high matter overdensity (smoothed on a scale of 14$\,h^{-1}$cMpc). This is reflected in the galaxy distribution --the most massive galaxies are found solely in regions of high overdensity. As a consequence of the above, galaxies in the EoR are highly clustered, which can lead to large variance in survey number counts. For mean number counts $N\lesssim 100$ (1000), in a unit redshift slice of angular area 300\,arcmin$^2$ (1.4\,deg$^2$), the 2-sigma range in $N$ is roughly a factor of four (two). We present relations between the expected variance and survey area for different survey geometries; these relations will be of use to observers wishing to understand the impact of survey variance on their results. △ Less

Submitted 14 June, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: Accepted version. Thanks to the referee contains more extensive description of testing, comparison to online Cosmic Variance Calculator, and greater range of survey areas than the original submission

FLARES IX: The Physical Mechanisms Driving Compact Galaxy Formation and Evolution

Authors: William J. Roper, Christopher C. Lovell, Aswin P. Vijayan, Dimitrios Irodotou, Jussi K. Kuusisto, Jasleen Matharu, Louise T. C. Seeyave, Peter A. Thomas, Stephen M. Wilkins

Abstract: In the FLARES (First Light And Reionisation Epoch Simulations) suite of hydrodynamical simulations, we find the high redshift ($z>5$) intrinsic size-luminosity relation is, surprisingly, negatively sloped. However, after including the effects of dust attenuation we find a positively sloped UV observed size-luminosity relation in good agreement with other simulated and observational studies. In thi… ▽ More In the FLARES (First Light And Reionisation Epoch Simulations) suite of hydrodynamical simulations, we find the high redshift ($z>5$) intrinsic size-luminosity relation is, surprisingly, negatively sloped. However, after including the effects of dust attenuation we find a positively sloped UV observed size-luminosity relation in good agreement with other simulated and observational studies. In this work, we extend this analysis to probe the underlying physical mechanisms driving the formation and evolution of the compact galaxies driving the negative size-mass/size-luminosity relation. We find the majority of compact galaxies ($R_{1/2, \star}< 1 \mathrm{pkpc}$), which drive the negative slope of the size-mass relation, have transitioned from extended to compact sizes via efficient centralised cooling, resulting in high specific star formation rates in their cores. These compact stellar systems are enshrouded by non-star forming gas distributions as much as $100\times$ larger than their stellar counterparts. By comparing with galaxies from the EAGLE simulation suite, we find that these extended gas distributions `turn on' and begin to form stars between $z=5$ and $z=0$ leading to increasing sizes, and thus the evolution of the size-mass relation from a negative to a positive slope. This explicitly demonstrates the process of inside-out galaxy formation in which compact bulges form earlier than the surrounding discs. △ Less

Submitted 6 February, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Molecular Strong Coupling and Cavity Finesse

Authors: Kishan S. Menghrajani, Adarsh B. Vasista, Wai Jue Tan, Philip A. Thomas, Felipe Herrera, William L. Barnes

Abstract: Molecular strong coupling offers exciting prospects in physics, chemistry and materials science. Whilst attention has been focused on including realistic models for the molecular systems involved, the important role played by the entire mode structure of the optical cavities employed so far has been largely overlooked. We show that the extent and effectiveness of molecular strong coupling is criti… ▽ More Molecular strong coupling offers exciting prospects in physics, chemistry and materials science. Whilst attention has been focused on including realistic models for the molecular systems involved, the important role played by the entire mode structure of the optical cavities employed so far has been largely overlooked. We show that the extent and effectiveness of molecular strong coupling is critically dependent on cavity finesse. Low finesse cavities can show strong coupling as judged from the presence of Rabi splitting in reflectivity measurements, but give photoluminescence signals equivalent to that of bare molecules outside the cavity. The emission of light is shown to involve hybridized light-matter polaritonic states only for cavities with high finesse. By developing an analytical model of cavity photoluminescence in multimode structures, we clarify the role of finite finesse in polariton formation. The detailed nature of the modes supported by a cavity and how these modes interact with the molecular system - whether on resonance or not - will be as important in developing a coherent framework for molecular strong coupling as the inclusion of realistic molecular models. △ Less

Submitted 18 March, 2024; v1 submitted 15 November, 2022; originally announced November 2022.

First Light And Reionisation Epoch Simulations (FLARES) VII: The Star Formation and Metal Enrichment Histories of Galaxies in the early Universe

Authors: Stephen M. Wilkins, Aswin P. Vijayan, Christopher C. Lovell, William J. Roper, Erik Zackrisson, Dimitrios Irodotou, Louise T. C. Seeyave, Jussi K. Kuusisto, Peter A. Thomas, Joseph Caruana, Christopher J. Conselice

Abstract: The star formation and metal enrichment histories of galaxies - at any epoch - constitute one of the key properties of galaxies, and their measurement is a core aim of observational extragalactic astronomy. The lack of deep rest-frame optical coverage at high-redshift has made robust constraints elusive, but this is now changing thanks to the \emph{James Webb Space Telescope (JWST)}. In preparatio… ▽ More The star formation and metal enrichment histories of galaxies - at any epoch - constitute one of the key properties of galaxies, and their measurement is a core aim of observational extragalactic astronomy. The lack of deep rest-frame optical coverage at high-redshift has made robust constraints elusive, but this is now changing thanks to the \emph{James Webb Space Telescope (JWST)}. In preparation for the constraints provided by \emph{JWST} we explore the star formation and metal enrichment histories of galaxies at $z=5-13$ using the First Light And Reionisation Epoch Simulations (FLARES) suite. Built on the EAGLE model, the unique strategy of FLARES allows us to simulate a wide range of stellar masses (and luminosities) and environments. While we predict significant redshift evolution of average ages and specific star formation rates our core result is a mostly flat relationship of age and specific star formation rate with stellar mass. We also find that galaxies in this epoch predominantly have strongly rising star formation histories, albeit with the magnitude dropping with redshift and stellar mass. In terms of chemical enrichment we predict a strong stellar mass - metallicity relation present at $z=10$ and beyond alongside significant $α$-enhancement. Finally, we find no environmental dependence of the relationship between age, specific star formation rate, or metallicity with stellar mass. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: To be submitted to MNRAS. Comments welcome

First Light And Reionisation Epoch Simulations (FLARES) VI: The colour evolution of galaxies $z=5-15$

Authors: Stephen M. Wilkins, Aswin P. Vijayan, Christopher C. Lovell, William J. Roper, Dimitrios Irodotou, Joseph Caruana, Louise T. C. Seeyave, Jussi K. Kuusisto, Peter A. Thomas

Abstract: With its exquisite sensitivity, wavelength coverage, and spatial and spectral resolution, the James Webb Space Telescope is poised to revolutionise our view of the distant, high-redshift ($z>5$) Universe. While Webb's spectroscopic observations will be transformative for the field, photometric observations play a key role in identifying distant objects and providing more comprehensive samples than… ▽ More With its exquisite sensitivity, wavelength coverage, and spatial and spectral resolution, the James Webb Space Telescope is poised to revolutionise our view of the distant, high-redshift ($z>5$) Universe. While Webb's spectroscopic observations will be transformative for the field, photometric observations play a key role in identifying distant objects and providing more comprehensive samples than accessible to spectroscopy alone. In addition to identifying objects, photometric observations can also be used to infer physical properties and thus be used to constrain galaxy formation models. However, inferred physical properties from broadband photometric observations, particularly in the absence of spectroscopic redshifts, often have large uncertainties. With the development of new tools for forward modelling simulations it is now routinely possible to predict observational quantities, enabling a direct comparison with observations. With this in mind, in this work, we make predictions for the colour evolution of galaxies at $z=5-15$ using the FLARES: First Light And Reionisation Epoch Simulations cosmological hydrodynamical simulation suite. We predict a complex evolution, driven predominantly by strong nebular line emission passing through individual bands. These predictions are in good agreement with existing constraints from Hubble and Spitzer as well as some of the first results from Webb. We also contrast our predictions with other models in the literature: while the general trends are similar we find key differences, particularly in the strength of features associated with strong nebular line emission. This suggests photometric observations alone should provide useful discriminating power between different models. △ Less

Submitted 6 September, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 10 pages, 9 figures, accepted for publication in MNRAS

Exploring the effect of baryons on the radial distribution of satellite galaxies with GAMA and IllustrisTNG

Authors: Stephen D. Riggs, Jon Loveday, Peter A. Thomas, Annalisa Pillepich, Dylan Nelson, Benne W. Holwerda

Abstract: We explore the radial distribution of satellite galaxies in groups in the Galaxy and Mass Assembly (GAMA) survey and the IllustrisTNG simulations. Considering groups with masses $12.0 \leq \log_{10} (\mathcal{M}_h / h^{-1} \mathrm{M}_{\odot}) < 14.8$ at $z<0.267$, we find a good agreement between GAMA and a sample of TNG300 groups and galaxies designed to match the GAMA selection. Both display a f… ▽ More We explore the radial distribution of satellite galaxies in groups in the Galaxy and Mass Assembly (GAMA) survey and the IllustrisTNG simulations. Considering groups with masses $12.0 \leq \log_{10} (\mathcal{M}_h / h^{-1} \mathrm{M}_{\odot}) < 14.8$ at $z<0.267$, we find a good agreement between GAMA and a sample of TNG300 groups and galaxies designed to match the GAMA selection. Both display a flat profile in the centre of groups, followed by a decline that becomes steeper towards the group edge, and normalised profiles show no dependence on group mass. Using matched satellites from TNG and dark matter-only TNG-Dark runs we investigate the effect of baryons on satellite radial location. At $z=0$, we find that the matched subhaloes from the TNG-Dark runs display a much flatter radial profile: namely, satellites selected above a minimum stellar mass exhibit both smaller halo-centric distances and longer survival times in the full-physics simulations compared to their dark-matter only analogues. We then divide the TNG satellites into those which possess TNG-Dark counterparts and those which do not, and develop models for the radial positions of each. We find the satellites with TNG-Dark counterparts are displaced towards the halo centre in the full-physics simulations, and this difference has a power-law behaviour with radius. For the `orphan' galaxies without TNG-Dark counterparts, we consider the shape of their radial distribution and provide a model for their motion over time, which can be used to improve the treatment of satellite galaxies in semi-analytic and semi-empirical models of galaxy formation. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: 20 pages, 22 figures. MNRAS, version after referee report

First Light And Reionisation Epoch Simulations (FLARES) V: The redshift frontier

Authors: Stephen M. Wilkins, Aswin P. Vijayan, Christopher C. Lovell, William J. Roper, Dimitrios Irodotou, Joseph Caruana, Louise T. C. Seeyave, Jussi K. Kuusisto, Peter A. Thomas, Shedeur A. K. Parris

Abstract: The James Webb Space Telescope (JWST) is set to transform many areas of astronomy, one of the most exciting is the expansion of the redshift frontier to $z>10$. In its first year alone JWST should discover hundreds of galaxies, dwarfing the handful currently known. To prepare for these powerful observational constraints, we use the First Light And Reionisation Epoch (FLARES) simulations to predict… ▽ More The James Webb Space Telescope (JWST) is set to transform many areas of astronomy, one of the most exciting is the expansion of the redshift frontier to $z>10$. In its first year alone JWST should discover hundreds of galaxies, dwarfing the handful currently known. To prepare for these powerful observational constraints, we use the First Light And Reionisation Epoch (FLARES) simulations to predict the physical and observational properties of the $z>10$ population of galaxies accessible to JWST. This is the first time such predictions have been made using a hydrodynamical model validated at low redshift. Our predictions at $z=10$ are broadly in agreement with current observational constraints on the far-UV luminosity function and UV continuum slope $β$, though the observational uncertainties are large. We note tension with recent constraints $z\sim 13$ from Harikane et al. 2022 - compared to these constraints, FLARES predicts objects with the same space density should have an order of magnitude lower luminosity, though this is mitigated slightly if dust attenuation is negligible in these systems. Our predictions suggest that in JWST's first cycle alone, around $600$ galaxies should be identified at $z>10$, with the first small samples available at $z>13$. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: Submitted to MNRAS 20/04/2022

First Light And Reionisation Epoch Simulations (FLARES) IV: The size evolution of galaxies at $z\geq5$

Authors: William J. Roper, Christopher C. Lovell, Aswin P. Vijayan, Madeline A. Marshall, Dimitrios Irodotou, Jussi K. Kuusisto, Peter A. Thomas, Stephen M. Wilkins

Abstract: We present the intrinsic and observed sizes of galaxies at $z\geq5$ in the First Light And Reionisation Epoch Simulations (FLARES). We employ the large effective volume of FLARES to produce a sizeable sample of high redshift galaxies with intrinsic and observed luminosities and half light radii in a range of rest frame UV and visual photometric bands. This sample contains a significant number of i… ▽ More We present the intrinsic and observed sizes of galaxies at $z\geq5$ in the First Light And Reionisation Epoch Simulations (FLARES). We employ the large effective volume of FLARES to produce a sizeable sample of high redshift galaxies with intrinsic and observed luminosities and half light radii in a range of rest frame UV and visual photometric bands. This sample contains a significant number of intrinsically ultra-compact galaxies in the far-UV (1500 angstrom), leading to a negative intrinsic far-UV size-luminosity relation. However, after the inclusion of the effects of dust these same compact galaxies exhibit observed sizes that are as much as 50 times larger than those measured from the intrinsic emission, and broadly agree with a range of observational samples. This increase in size is driven by the concentration of dust in the core of galaxies, heavily attenuating the intrinsically brightest regions. At fixed luminosity we find a galaxy size redshift evolution with a slope of $m=1.21-1.87$ depending on the luminosity sample in question, and we demonstrate the wavelength dependence of the size-luminosity relation which will soon be probed by the Webb Space Telescope. △ Less

Submitted 23 March, 2022; originally announced March 2022.

All-optical control of phase singularities using strong light-matter coupling

Authors: Philip A. Thomas, Kishan S. Menghrajani, William L. Barnes

Abstract: Strong light-matter coupling occurs when the coupling strength between a confined electromagnetic mode and a molecular resonance exceeds losses to the environment. The study of strong coupling has been motivated by applications such as lasing and the modification of chemical processes. Here we show that strong coupling can be used to create phase singularities. Many nanophotonic structures have be… ▽ More Strong light-matter coupling occurs when the coupling strength between a confined electromagnetic mode and a molecular resonance exceeds losses to the environment. The study of strong coupling has been motivated by applications such as lasing and the modification of chemical processes. Here we show that strong coupling can be used to create phase singularities. Many nanophotonic structures have been designed to generate phase singularities for use in sensing and optoelectronics. We utilise the concept of cavity-free strong coupling, where electromagnetic modes sustained by a material are strong enough to strongly couple to the material's own molecular resonance, to create phase singularities in a simple thin film of organic molecules. We show that the use of photochromic molecules allows for all-optical control of phase singularities. Our results suggest a new application for strong light-matter coupling and a new, simplified, more versatile pathway to singular phase optics. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 18 pages, 6 figures

First Light And Reionisation Epoch Simulations (FLARES) III: The properties of massive dusty galaxies at cosmic dawn

Authors: Aswin P. Vijayan, Stephen M. Wilkins, Christopher C. Lovell, Peter A. Thomas, Peter Camps, Maarten Baes, James Trayford, Jussi Kuusisto, William J. Roper

Abstract: Using the First Light And Reionisation Epoch Simulations (\textsc{Flares}) we explore the dust driven properties of massive high-redshift galaxies at $z\in[5,10]$. By post-processing the galaxy sample using the radiative transfer code \textsc{skirt} we obtain the full spectral energy distribution. We explore the resultant luminosity functions, IRX-$β$ relations as well as the luminosity-weighted d… ▽ More Using the First Light And Reionisation Epoch Simulations (\textsc{Flares}) we explore the dust driven properties of massive high-redshift galaxies at $z\in[5,10]$. By post-processing the galaxy sample using the radiative transfer code \textsc{skirt} we obtain the full spectral energy distribution. We explore the resultant luminosity functions, IRX-$β$ relations as well as the luminosity-weighted dust temperatures in the Epoch of Reionisation (EoR). We find that most of our results are in agreement with the current set of observations, but under-predict the number densities of bright IR galaxies, which are extremely biased towards the most overdense regions. We see that the \textsc{Flares} IRX-$β$ relation (for $5\le z\le8$) predominantly follows the local starburst relation. The IRX shows an increase with stellar mass, plateauing at the high-mass end ($\sim10^{10}$M$_{\odot}$) and shows no evolution in the median normalisation with redshift. We also look at the dependence of the peak dust temperature ($T_{\mathrm{peak}}$) on various galaxy properties including the stellar mass, IR luminosity and sSFR, finding the correlation to be strongest with sSFR. The luminosity-weighted dust temperatures increase towards higher redshifts, with the slope of the $T_{\mathrm{peak}}$ - redshift relation showing a higher slope than the lower redshift relations obtained from previous observational and theoretical works. The results from \textsc{Flares}, which is able to provide a better statistical sample of high-redshift galaxies compared to other simulations, provides a distinct vantage point for the high-redshift Universe. △ Less

Submitted 10 March, 2022; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 19 pages, 16 figures, Published in MNRAS. Data link: https://flaresimulations.github.io/data.html

A machine learning approach to mapping baryons onto dark matter haloes using the EAGLE and C-EAGLE simulations

Authors: Christopher C. Lovell, Stephen M. Wilkins, Peter A. Thomas, Matthieu Schaller, Carlton M. Baugh, Giulio Fabbian, Yannick Bahé

Abstract: High-resolution cosmological hydrodynamic simulations are currently limited to relatively small volumes due to their computational expense. However, much larger volumes are required to probe rare, overdense environments, and measure clustering statistics of the large scale structure. Typically, zoom simulations of individual regions are used to study rare environments, and semi-analytic models and… ▽ More High-resolution cosmological hydrodynamic simulations are currently limited to relatively small volumes due to their computational expense. However, much larger volumes are required to probe rare, overdense environments, and measure clustering statistics of the large scale structure. Typically, zoom simulations of individual regions are used to study rare environments, and semi-analytic models and halo occupation models applied to dark matter only (DMO) simulations are used to study the Universe in the large-volume regime. We propose a new approach, using a machine learning framework to explore the halo-galaxy relationship in the periodic EAGLE simulations, and zoom C-EAGLE simulations of galaxy clusters. We train a tree based machine learning method to predict the baryonic properties of galaxies based on their host dark matter halo properties. The trained model successfully reproduces a number of key distribution functions for an infinitesimal fraction of the computational cost of a full hydrodynamic simulation. By training on both periodic simulations as well as zooms of overdense environments, we learn the bias of galaxy evolution in differing environments. This allows us to apply the trained model to a larger DMO volume than would be possible if we only trained on a periodic simulation. We demonstrate this application using the $(800 \; \mathrm{Mpc})^3$ P-Millennium simulation, and present predictions for key baryonic distribution functions and clustering statistics from the EAGLE model in this large volume. △ Less

Submitted 2 May, 2023; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 17 pages, 16 figures, Accepted to MNRAS. Updated Acknowledgements

L-GALAXIES 2020: The evolution of radial metallicity profiles and global metallicities in disc galaxies

Authors: Robert M. Yates, Bruno M. B. Henriques, Jian Fu, Guinevere Kauffmann, Peter A. Thomas, Qi Guo, Simon D. M. White, Patricia Schady

Abstract: We present a modified version of the L-GALAXIES 2020 semi-analytic model of galaxy evolution, which includes significantly increased direct metal enrichment of the circumgalactic medium (CGM) by supernovae (SNe). These more metal-rich outflows do not require increased mass-loading factors, in contrast to some other galaxy evolution models. This modified L-GALAXIES 2020 model is able to simultaneou… ▽ More We present a modified version of the L-GALAXIES 2020 semi-analytic model of galaxy evolution, which includes significantly increased direct metal enrichment of the circumgalactic medium (CGM) by supernovae (SNe). These more metal-rich outflows do not require increased mass-loading factors, in contrast to some other galaxy evolution models. This modified L-GALAXIES 2020 model is able to simultaneously reproduce the gas-phase metallicity $(Z_{\rm g})$ and stellar metallicity $(Z_{*})$ radial profiles observed in nearby disc galaxies by MaNGA and MUSE, as well as the observed mass - metallicity relations for gas and stars at $z=0$ and their evolution back to $z\sim{}2-3$. A direct CGM enrichment fraction of $\sim{}90\%$ for SNe-II is preferred. We find that massive disc galaxies have slightly flatter $Z_{\rm g}$ profiles than their lower-mass counterparts in L-GALAXIES 2020, due to more efficient enrichment of their outskirts via inside-out growth and metal-rich accretion. Such a weak, positive correlation between stellar mass and $Z_{\rm g}$ profile slope is also seen in our MaNGA-DR15 sample of 571 star-forming disc galaxies. Although, below ${\rm log}(M_{*}/{\rm M}_{\odot})\sim{}10.0$ this observational result is strongly dependent on the metallicity diagnostic and morphological selection chosen. In addition, a lowered maximum SN-II progenitor mass of $25{\rm M}_{\odot}$, reflecting recent theoretical and observational estimates, can also provide a good match to observed metallicity profiles at $z=0$ in L-GALAXIES 2020. However, this model version fails to reproduce an evolution in $Z_{\rm g}$ at fixed mass over cosmic time, or the magnesium abundances observed in the intracluster medium (ICM). △ Less

Submitted 3 June, 2021; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: 22 pages, 19 figures. Accepted by MNRAS. (Missing citation hyperlinks corrected)

Using angular momentum maps to detect kinematically distinct galactic components

Authors: Dimitrios Irodotou, Peter A. Thomas

Abstract: In this work we introduce a physically motivated method of performing disc/spheroid decomposition of simulated galaxies, which we apply to the Eagle sample. We make use of the HEALPix package to create Mollweide projections of the angular momentum map of each galaxy's stellar particles. A number of features arise on the angular momentum space which allows us to decompose galaxies and classify them… ▽ More In this work we introduce a physically motivated method of performing disc/spheroid decomposition of simulated galaxies, which we apply to the Eagle sample. We make use of the HEALPix package to create Mollweide projections of the angular momentum map of each galaxy's stellar particles. A number of features arise on the angular momentum space which allows us to decompose galaxies and classify them into different morphological types. We assign stellar particles with angular separation of less/greater than 30 degrees from the densest grid cell on the angular momentum sphere to the disc/spheroid components, respectively. We analyse the spatial distribution for a subsample of galaxies and show that the surface density profiles of the disc and spheroid closely follow an exponential and a Sersic profile, respectively. In addition discs rotate faster, have smaller velocity dispersions, are younger and are more metal rich than spheroids. Thus our morphological classification reproduces the observed properties of such systems. Finally, we demonstrate that our method is able to identify a significant population of galaxies with counter-rotating discs and provide a more realistic classification of such systems compared to previous methods. △ Less

Submitted 13 January, 2021; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 16 pages, 15 figures, accepted by MNRAS

First Light And Reionisation Epoch Simulations (FLARES) II: The Photometric Properties of High-Redshift Galaxies

Authors: Aswin P. Vijayan, Christopher C. Lovell, Stephen M. Wilkins, Peter A. Thomas, David J. Barnes, Dimitrios Irodotou, Jussi Kuusisto, Will Roper

Abstract: We present the photometric properties of galaxies in the First Light and Reionisation Epoch Simulations (FLARES). The simulations trace the evolution of galaxies in a range of overdensities through the Epoch of Reionistion (EoR). With a novel weighting scheme we combine these overdensities, extending significantly the dynamic range of observed composite distribution functions compared to periodic… ▽ More We present the photometric properties of galaxies in the First Light and Reionisation Epoch Simulations (FLARES). The simulations trace the evolution of galaxies in a range of overdensities through the Epoch of Reionistion (EoR). With a novel weighting scheme we combine these overdensities, extending significantly the dynamic range of observed composite distribution functions compared to periodic simulation boxes. FLARES predicts a significantly larger number of intrinsically bright galaxies, which can be explained through a simple model linking dust-attenuation to the metal content of the interstellar medium, using a line-of-sight (LOS) extinction model. With this model we present the photometric properties of the FLARES galaxies for $z \in [5,10]$. We show that the ultraviolet (UV) luminosity function (LF) matches the observations at all redshifts. The function is fit by Schechter and double power-law forms, with the latter being favoured at these redshifts by the FLARES composite UV LF. We also present predictions for the UV continuum slope as well as the attenuation in the UV. The impact of environment on the UV LF is also explored, with the brightest galaxies forming in the densest environments. We then present the line luminosity and equivalent widths of some prominent nebular emission lines arising from the galaxies, finding rough agreement with available observations. We also look at the relative contribution of obscured and unobscured star formation, finding comparable contributions at these redshifts. △ Less

Submitted 14 December, 2020; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: 21 pages, 22 figures. Data available at https://flaresimulations.github.io/data.html

A new signature for strong light-matter coupling using spectroscopic ellipsometry

Authors: Philip A. Thomas, Wai Jue Tan, Henry A. Fernandez, William L. Barnes

Abstract: Light-matter interactions can occur when an ensemble of molecular resonators is placed in a confined electromagnetic field. In the strong coupling regime the rapid exchange of energy between the molecules and the electromagnetic field results in the emergence of hybrid light-matter states called polaritons. Multiple criteria exist to define the strong coupling regime, usually by comparing the spli… ▽ More Light-matter interactions can occur when an ensemble of molecular resonators is placed in a confined electromagnetic field. In the strong coupling regime the rapid exchange of energy between the molecules and the electromagnetic field results in the emergence of hybrid light-matter states called polaritons. Multiple criteria exist to define the strong coupling regime, usually by comparing the splitting of the polariton bands with the linewidths of the uncoupled modes. Here we highlight the limitations of these criteria and study strong coupling using spectroscopic ellipsometry, a commonly used optical characterisation technique. We identify a new signature of strong coupling in ellipsometric phase spectra. Combining ellipsometric amplitude and phase spectra yields a distinct topological feature that we suggest could serve as a new criterion for strong coupling. Our results introduce the idea of ellipsometric topology and could provide further insight into the transition from the weak to strong coupling regime. △ Less

Submitted 7 August, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: 20 pages, 10 figures

First Light And Reionisation Epoch Simulations (FLARES) I: Environmental Dependence of High-Redshift Galaxy Evolution

Authors: Christopher C. Lovell, Aswin P. Vijayan, Peter A. Thomas, Stephen M. Wilkins, David J. Barnes, Dimitrios Irodotou, Will Roper

Abstract: We introduce the First Light And Reionisation Epoch Simulations (FLARES), a suite of zoom simulations using the EAGLE model. We resimulate a range of overdensities during the Epoch of Reionisation (EoR) in order to build composite distribution functions, as well as explore the environmental dependence of galaxy formation and evolution during this critical period of galaxy assembly. The regions are… ▽ More We introduce the First Light And Reionisation Epoch Simulations (FLARES), a suite of zoom simulations using the EAGLE model. We resimulate a range of overdensities during the Epoch of Reionisation (EoR) in order to build composite distribution functions, as well as explore the environmental dependence of galaxy formation and evolution during this critical period of galaxy assembly. The regions are selected from a large $(3.2 \;\mathrm{cGpc})^{3}$ parent volume, based on their overdensity within a sphere of radius $14\,h^{-1}\;\mathrm{cMpc}$. We then resimulate with full hydrodynamics, and employ a novel weighting scheme that allows the construction of composite distribution functions that are representative of the full parent volume. This significantly extends the dynamic range compared to smaller volume periodic simulations. We present an analysis of the galaxy stellar mass function (GSMF), the star formation rate distribution function (SFRF) and the star forming sequence (SFS) predicted by \flares, and compare to a number of observational and model constraints. We also analyse the environmental dependence over an unprecedented range of overdensity. Both the GSMF and the SFRF exhibit a clear double-Schechter form, up to the highest redshifts ($z = 10$). We also find no environmental dependence of the SFS normalisation. The increased dynamic range probed by FLARES will allow us to make predictions for a number of large area surveys that will probe the EoR in coming years, such as WFIRST and Euclid. △ Less

Submitted 2 September, 2020; v1 submitted 15 April, 2020; originally announced April 2020.

Comments: 21 pages, 19 figures, submitted to MNRAS

MEGA: Merger graphs of structure formation

Authors: William J. Roper, Peter A. Thomas, Chaichalit Srisawat

Abstract: When following the growth of structure in the Universe, we propose replacing merger trees with merger graphs, in which haloes can both merge and split into separate pieces. We show that this leads to smoother mass growth and eliminates catastrophic failures in which massive haloes have no progenitors or descendants. For those who prefer to stick with merger trees, we find that trees derived from o… ▽ More When following the growth of structure in the Universe, we propose replacing merger trees with merger graphs, in which haloes can both merge and split into separate pieces. We show that this leads to smoother mass growth and eliminates catastrophic failures in which massive haloes have no progenitors or descendants. For those who prefer to stick with merger trees, we find that trees derived from our merger graphs have similar mass growth properties to previous methods, but again without catastrophic failures. For future galaxy formation modelling, two different density thresholds can be used to distinguish host haloes (extended galactic haloes, groups and clusters) from higher-density subhaloes: sites of galaxy formation. △ Less

Submitted 7 May, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Detailed dust modelling in the L-Galaxies semi-analytic model of galaxy formation

Authors: Aswin P. Vijayan, Scott J. Clay, Peter A. Thomas, Robert M. Yates, Stephen M. Wilkins, Bruno M. Henriques

Abstract: We implement a detailed dust model into the L-Galaxies semi-analytical model which includes: injection of dust by type II and type Ia supernovae (SNe) and AGB stars; grain growth in molecular clouds; and destruction due to supernova-induced shocks, star formation, and reheating. Our grain growth model follows the dust content in molecular clouds and the inter-cloud medium separately, and allows gr… ▽ More We implement a detailed dust model into the L-Galaxies semi-analytical model which includes: injection of dust by type II and type Ia supernovae (SNe) and AGB stars; grain growth in molecular clouds; and destruction due to supernova-induced shocks, star formation, and reheating. Our grain growth model follows the dust content in molecular clouds and the inter-cloud medium separately, and allows growth only on pre-existing dust grains. At early times, this can make a significant difference to the dust growth rate. Above $z\sim8$, type II SNe are the primary source of dust, whereas below $z\sim8$, grain growth in molecular clouds dominates, with the total dust content being dominated by the latter below $z\sim6$. However, the detailed history of galaxy formation is important for determining the dust content of any individual galaxy. We introduce a fit to the dust-to-metal (DTM) ratio as a function of metallicity and age, which can be used to deduce the DTM ratio of galaxies at any redshift. At $z\lesssim3$, we find a fairly flat mean relation between metallicity and the DTM, and a positive correlation between metallicity and the dust-to-gas (DTG) ratio, in good agreement with the shape and normalisation of the observed relations. We also match the normalisation of the observed stellar mass -- dust mass relation over the redshift range of $0-4$, and to the dust mass function at $z=0$. Our results are important in interpreting observations on the dust content of galaxies across cosmic time, particularly so at high redshift. △ Less

Submitted 5 July, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 19 pages, 16 figures. Edits: Corrected equation 15 and its parameters. Accepted in MNRAS

Learning the Relationship between Galaxies Spectra and their Star Formation Histories using Convolutional Neural Networks and Cosmological Simulations

Authors: Christopher C. Lovell, Viviana Acquaviva, Peter A. Thomas, Kartheik G. Iyer, Eric Gawiser, Stephen M. Wilkins

Abstract: We present a new method for inferring galaxy star formation histories (SFH) using machine learning methods coupled with two cosmological hydrodynamic simulations. We train Convolutional Neural Networks to learn the relationship between synthetic galaxy spectra and high resolution SFHs from the EAGLE and Illustris models. To evaluate our SFH reconstruction we use Symmetric Mean Absolute Percentage… ▽ More We present a new method for inferring galaxy star formation histories (SFH) using machine learning methods coupled with two cosmological hydrodynamic simulations. We train Convolutional Neural Networks to learn the relationship between synthetic galaxy spectra and high resolution SFHs from the EAGLE and Illustris models. To evaluate our SFH reconstruction we use Symmetric Mean Absolute Percentage Error (SMAPE), which acts as a true percentage error in the low-error regime. On dust-attenuated spectra we achieve high test accuracy (median SMAPE $= 10.5\%$). Including the effects of simulated observational noise increases the error ($12.5\%$), however this is alleviated by including multiple realisations of the noise, which increases the training set size and reduces overfitting ($10.9\%$). We also make estimates for the observational and modelling errors. To further evaluate the generalisation properties we apply models trained on one simulation to spectra from the other, which leads to only a small increase in the error (median SMAPE $\sim 15\%$). We apply each trained model to SDSS DR7 spectra, and find smoother histories than in the VESPA catalogue. This new approach complements the results of existing SED fitting techniques, providing star formation histories directly motivated by the results of the latest cosmological simulations. △ Less

Submitted 9 October, 2019; v1 submitted 25 March, 2019; originally announced March 2019.

Comments: 20 pages, 21 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Morphological evolution and galactic sizes in the L-Galaxies SA model

Authors: Dimitrios Irodotou, Peter A. Thomas, Bruno M. Henriques, Mark T. Sargent, Jessica M. Hislop

Abstract: In this work we update the L-Galaxies semi-analytic model (SAM) to better follow the physical processes responsible for the growth of bulges via disc instabilities (leading to pseudo-bulges) and mergers (leading to classical bulges). We address the former by considering the contribution of both stellar and gaseous discs in the stability of the galaxy, and we update the latter by including dissipat… ▽ More In this work we update the L-Galaxies semi-analytic model (SAM) to better follow the physical processes responsible for the growth of bulges via disc instabilities (leading to pseudo-bulges) and mergers (leading to classical bulges). We address the former by considering the contribution of both stellar and gaseous discs in the stability of the galaxy, and we update the latter by including dissipation of energy in gas-rich mergers. Furthermore, we introduce angular momentum losses during cooling and find that an accurate match to the observed correlation between stellar disc scale length and mass at z ~ 0.0 requires that the gas loses 20% of its initial specific angular momentum to the corresponding dark matter halo during the formation of the cold gas disc. We reproduce the observed trends between the stellar mass and specific angular momentum for both disc- and bulge-dominated galaxies, with the former rotating faster than the latter of the same mass. We conclude that a two-component instability recipe provides a morphologically diverse galaxy sample which matches the observed fractional breakdown of galaxies into different morphological types. This recipe also enables us to obtain an excellent fit to the morphology-mass relation and stellar mass function of different galactic types. Finally, we find that energy dissipation during mergers reduces the merger remnant sizes and allows us to match the observed mass-size relation for bulge-dominated systems. △ Less

Submitted 7 October, 2019; v1 submitted 11 October, 2018; originally announced October 2018.

Cosmic CARNage II: the evolution of the galaxy stellar mass function in observations and galaxy formation models

Authors: Rachel Asquith, Frazer R. Pearce, Omar Almaini, Alexander Knebe, Violeta Gonzalez-Perez, Andrew Benson, Jeremy Blaizot, Jorge Carretero, Francisco J. Castander, Andrea Cattaneo, Sofía A. Cora, Darren J. Croton, Julien E. Devriendt, Fabio Fontanot, Ignacio D. Gargiulo, Will Hartley, Bruno Henriques, Jaehyun Lee, Gary A. Mamon, Julian Onions, Nelson D. Padilla, Chris Power, Chaichalit Srisawat, Adam R. H. Stevens, Peter A. Thomas , et al. (2 additional authors not shown)

Abstract: We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of… ▽ More We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift. △ Less

Submitted 27 July, 2018; v1 submitted 10 July, 2018; originally announced July 2018.

Comments: 16 pages, 8 figures, accepted for publication in MNRAS

Cosmic CARNage I: on the calibration of galaxy formation models

Authors: Alexander Knebe, Frazer R. Pearce, Violeta Gonzalez-Perez, Peter A. Thomas, Andrew Benson, Rachel Asquith, Jeremy Blaizot, Richard Bower, Jorge Carretero, Francisco J. Castander, Andrea Cattaneo, Sofia A. Cora, Darren J. Croton, Weiguang Cui, Daniel Cunnama, Julien E. Devriendt, Pascal J. Elahi, Andreea Font, Fabio Fontanot, Ignacio D. Gargiulo, John Helly, Bruno Henriques, Jaehyun Lee, Gary A. Mamon, Julian Onions , et al. (9 additional authors not shown)

Abstract: We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational… ▽ More We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-σ. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: 20 pages, 14 figures (+5 supplementary figures in the Appendix), accepted for publication in MNRAS

Characterising and Identifying Galaxy Protoclusters

Authors: Christopher C. Lovell, Peter A. Thomas, Stephen M. Wilkins

Abstract: We study the characteristics of galaxy protoclusters using the latest L-galaxies semi-analytic model. Searching for protoclusters on a scale of $\sim 10 \, \mathrm{cMpc}$ gives an excellent compromise between the completeness and purity of their galaxy populations, leads to high distinction from the field in overdensity space, and allows accurate determination of the descendant cluster mass. This… ▽ More We study the characteristics of galaxy protoclusters using the latest L-galaxies semi-analytic model. Searching for protoclusters on a scale of $\sim 10 \, \mathrm{cMpc}$ gives an excellent compromise between the completeness and purity of their galaxy populations, leads to high distinction from the field in overdensity space, and allows accurate determination of the descendant cluster mass. This scale is valid over a range of redshifts and selection criteria. We present a procedure for estimating, given a measured galaxy overdensity, the protocluster probability and its descendant cluster mass for a range of modelling assumptions, particularly taking into account the shape of the measurement aperture. This procedure produces lower protocluster probabilities compared to previous estimates using fixed size apertures. The relationship between AGN and protoclusters is also investigated, and shows significant evolution with redshift; at $z \sim 2$ the fraction of protoclusters traced by AGN is high, but the fraction of all AGN in protoclusters is low, whereas at $z \geqslant 5$ the fraction of protoclusters containing AGN is low, but most AGN are in protoclusters. We also find indirect evidence for the emergence of a passive sequence in protoclusters at $z \sim 2$, and note that a significant fraction of all galaxies reside in protoclusters at $z \geqslant 2$, particularly the most massive. △ Less

Submitted 1 December, 2017; v1 submitted 5 October, 2017; originally announced October 2017.

Comments: 19 pages, 16 figures, submitted to MNRAS

The Cluster-EAGLE project: global properties of simulated clusters with resolved galaxies

Authors: David J. Barnes, Scott T. Kay, Yannick M. Bahe, Claudio Dalla Vecchia, Ian G. McCarthy, Joop Schaye, Richard G. Bower, Adrian Jenkins, Peter A. Thomas, Matthieu Schaller, Robert A. Crain, Tom Theuns, Simon D. M. White

Abstract: We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bahé et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of… ▽ More We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bahé et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of $1.8\times10^{6}\,\mathrm{M}_{\odot}$ and physical softening length of $0.7\,\mathrm{kpc}$. In this paper, we introduce the sample and present the low-redshift global properties of the clusters. We calculate the X-ray properties in a manner consistent with observational techniques, demonstrating the bias and scatter introduced by using estimated masses. We find the total stellar content and black hole masses of the clusters to be in good agreement with the observed relations. However, the clusters are too gas rich, suggesting that the AGN feedback model is not efficient enough at expelling gas from the high-redshift progenitors of the clusters. The X-ray properties, such as the spectroscopic temperature and the soft-band luminosity, and the Sunyaev-Zel'dovich properties are in reasonable agreement with the observed relations. However, the clusters have too high central temperatures and larger-than-observed entropy cores, which is likely driven by the AGN feedback after the cluster core has formed. The total metal content and its distribution throughout the ICM are a good match to the observations. △ Less

Submitted 6 July, 2017; v1 submitted 31 March, 2017; originally announced March 2017.

Comments: 21 pages, 15 figures, MNRAS accepted version

The Hydrangea simulations: galaxy formation in and around massive clusters

Authors: Yannick M. Bahé, David J. Barnes, Claudio Dalla Vecchia, Scott T. Kay, Simon D. M. White, Ian G. McCarthy, Joop Schaye, Richard G. Bower, Robert A. Crain, Tom Theuns, Adrian Jenkins, Sean L. McGee, Matthieu Schaller, Peter A. Thomas, James W. Trayford

Abstract: We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for th… ▽ More We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses Mstar > 10^10 M_Sun this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (Mstar < 10^10 M_Sun) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching. Conversely, the simulations predict an overabundance of massive galaxies in clusters compared to the field that persists to their far outskirts (> 5r_200c). This is caused by a significantly increased stellar mass fraction of (sub-)haloes in the cluster environment, by up to ~0.3 dex even well beyond r_200c. Haloes near clusters are also more concentrated than equally massive field haloes, but these two effects are largely uncorrelated. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Comments: 22 pages, 15 figures (18 pages, 14 figures without appendices); submitted to MNRAS

nIFTy Cosmology: the clustering consistency of galaxy formation models

Authors: Arnau Pujol, Ramin A. Skibba, Enrique Gaztañaga, Andrew Benson, Jeremy Blaizot, Richard Bower, Jorge Carretero, Francisco J. Castander, Andrea Cattaneo, Sofia A. Cora, Darren J. Croton, Weiguang Cui, Daniel Cunnama, Gabriella De Lucia, Julien E. Devriendt, Pascal J. Elahi, Andreea Font, Fabio Fontanot, Juan Garcia-Bellido, Ignacio D. Gargiulo, Violeta Gonzalez-Perez, John Helly, Bruno M. B. Henriques, Michaela Hirschmann, Alexander Knebe , et al. (13 additional authors not shown)

Abstract: We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single ΛCDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Fun… ▽ More We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single ΛCDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs. △ Less

Submitted 13 April, 2017; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: 17 pages, 7 figures

Galaxy formation in the Planck cosmology - IV. Mass and environmental quenching, conformity and clustering

Authors: Bruno M. B. Henriques, Simon D. M. White, Peter A. Thomas, Raul E. Angulo, Qi Guo, Gerard Lemson, Wenting Wang

Abstract: We study the quenching of star formation as a function of redshift, environment and stellar mass in the galaxy formation simulations of Henriques et al. (2015), which implement an updated version of the Munich semi-analytic model (L-GALAXIES) on the two Millennium Simulations after scaling to a Planck cosmology. In this model massive galaxies are quenched by AGN feedback depending on both black ho… ▽ More We study the quenching of star formation as a function of redshift, environment and stellar mass in the galaxy formation simulations of Henriques et al. (2015), which implement an updated version of the Munich semi-analytic model (L-GALAXIES) on the two Millennium Simulations after scaling to a Planck cosmology. In this model massive galaxies are quenched by AGN feedback depending on both black hole and hot gas mass, and hence indirectly on stellar mass. In addition, satellite galaxies of any mass can be quenched by ram-pressure or tidal stripping of gas and through the suppression of gaseous infall. This combination of processes produces quenching efficiencies which depend on stellar mass, host halo mass, environment density, distance to group centre and group central galaxy properties in ways which agree qualitatively with observation. Some discrepancies remain in dense regions and close to group centres, where quenching still seems too efficient. In addition, although the mean stellar age of massive galaxies agrees with observation, the assumed AGN feedback model allows too much ongoing star formation at late times. The fact that both AGN feedback and environmental effects are stronger in higher density environments leads to a correlation between the quenching of central and satellite galaxies which roughly reproduces observed conformity trends inside haloes. △ Less

Submitted 7 November, 2016; originally announced November 2016.

Comments: Submitted to MNRAS; 21 pages, 16 figures

A consistent model for both the HI and stellar mass functions of galaxies

Authors: Hazel Martindale, Peter A. Thomas, Bruno M. Henriques, Jon Loveday

Abstract: Using the L-Galaxies semi-analytic model we simultaneously fit the HI mass function, stellar mass function and galaxy colours. We find good fits to all three observations at z = 0 and to the stellar mass function and galaxy colours at z = 2. Using Markov Chain Monte Carlo (MCMC) techniques we adjust the L-Galaxies parameters to best fit the constraining data. In order to fit the HI mass function w… ▽ More Using the L-Galaxies semi-analytic model we simultaneously fit the HI mass function, stellar mass function and galaxy colours. We find good fits to all three observations at z = 0 and to the stellar mass function and galaxy colours at z = 2. Using Markov Chain Monte Carlo (MCMC) techniques we adjust the L-Galaxies parameters to best fit the constraining data. In order to fit the HI mass function we must greatly reduce the gas surface density threshold for star formation, thus lowering the number of low HI mass galaxies. A simultaneous reduction in the star formation efficiency prevents the over production of stellar content. A simplified model in which the surface density threshold is eliminated altogether also provides a good fit to the data. Unfortunately, these changes weaken the fit to the Kennicutt-Schmidt relation and raise the star-formation rate density at recent times, suggesting that a change to the model is required to prevent accumulation of gas onto dwarf galaxies in the local universe. △ Less

Submitted 27 June, 2016; originally announced June 2016.

Comments: 10 pages, 10 figures, submitted to mnras

Sussing Merger Trees: Stability and Convergence

Authors: Yang Wang, Frazer R. Pearce, Alexander Knebe, Aurel Schneider, Chaichalit Srisawat, Dylan Tweed, Intae Jung, Jiaxin Han, John Helly, Julian Onions, Pascal J. Elahi, Peter A. Thomas, Peter Behroozi, Sukyoung K. Yi, Vicente Rodriguez-Gomez, Yao-Yuan Mao, Yipeng Jing, Weipeng Lin

Abstract: Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not… ▽ More Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots. △ Less

Submitted 5 April, 2016; originally announced April 2016.

Comments: 16 pages, 14 figures, accepted by MNRAS

Iron in galaxy groups and clusters: Confronting galaxy evolution models with a newly homogenised dataset

Authors: Robert M. Yates, Peter A. Thomas, Bruno M. B. Henriques

Abstract: We present an analysis of the iron abundance in the hot gas surrounding galaxy groups and clusters. To do this, we first compile and homogenise a large dataset of 79 low-redshift (|z| = 0.03) systems (159 individual measurements) from the literature. Our analysis accounts for differences in aperture size, solar abundance, and cosmology, and scales all measurements using customised radial profiles… ▽ More We present an analysis of the iron abundance in the hot gas surrounding galaxy groups and clusters. To do this, we first compile and homogenise a large dataset of 79 low-redshift (|z| = 0.03) systems (159 individual measurements) from the literature. Our analysis accounts for differences in aperture size, solar abundance, and cosmology, and scales all measurements using customised radial profiles for the temperature (T), gas density, and iron abundance (Z). We then compare this dataset to groups and clusters in the L-Galaxies galaxy evolution model. Our homogenised dataset reveals a tight T-Z relation for clusters, with a scatter in Z of only 0.10 dex and a slight negative gradient. After examining potential measurement biases, we conclude that at least some of this negative gradient has a physical origin. Our model suggests greater accretion of hydrogen in the hottest systems, via stripping of gas from infalling satellites, as a cause. At lower temperatures, L-Galaxies over-estimates Z in groups, indicating that metal-rich gas removal (via e.g. AGN feedback) is required. L-Galaxies provides a reasonable match to the observed Z in the intracluster medium (ICM) of the hottest clusters from at least z ~ 1.3 to 0.3. This is achieved without needing to modify any of the galactic chemical evolution (GCE) model parameters. However, the Z in intermediate-temperature clusters appears to be under-estimated in our model at z = 0. The merits and problems with modifying the GCE modelling to correct this are discussed. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 19 pages, 14 figures (excluding appendices), submitted to MNRAS

The XMM Cluster Survey: evolution of the velocity dispersion -- temperature relation over half a Hubble time

Authors: Susan Wilson, Matt Hilton, Philip J. Rooney, Caroline Caldwell, Scott T. Kay, Chris A. Collins, Ian G. McCarthy, A. Kathy Romer, Alberto Bermeo-Hernandez, Rebecca Bernstein, Luiz da Costa, Daniel Gifford, Devon Hollowood, Ben Hoyle, Tesla Jeltema, Andrew R. Liddle, Marcio A. G Maia, Robert G. Mann, Julian A. Mayers, Nicola Mehrtens, Christopher J. Miller, Robert C. Nichol, Ricardo Ogando, Martin Sahlén, Benjamin Stahl , et al. (4 additional authors not shown)

Abstract: We measure the evolution of the velocity dispersion--temperature ($σ_{\rm v}$--$T_{\rm X}$) relation up to $z = 1$ using a sample of 38 galaxy clusters drawn from the \textit{XMM} Cluster Survey. This work improves upon previous studies by the use of a homogeneous cluster sample and in terms of the number of high redshift clusters included. We present here new redshift and velocity dispersion meas… ▽ More We measure the evolution of the velocity dispersion--temperature ($σ_{\rm v}$--$T_{\rm X}$) relation up to $z = 1$ using a sample of 38 galaxy clusters drawn from the \textit{XMM} Cluster Survey. This work improves upon previous studies by the use of a homogeneous cluster sample and in terms of the number of high redshift clusters included. We present here new redshift and velocity dispersion measurements for 12 $z > 0.5$ clusters observed with the GMOS instruments on the Gemini telescopes. Using an orthogonal regression method, we find that the slope of the relation is steeper than that expected if clusters were self-similar, and that the evolution of the normalisation is slightly negative, but not significantly different from zero ($σ_{\rm v} \propto T^{0.86 \pm 0.14} E(z)^{-0.37 \pm 0.33}$). We verify our results by applying our methods to cosmological hydrodynamical simulations. The lack of evolution seen in our data is consistent with simulations that include both feedback and radiative cooling. △ Less

Submitted 3 August, 2016; v1 submitted 9 December, 2015; originally announced December 2015.

Comments: Accepted to MNRAS (3 August 2016); Paper: 15 pages, 12 figures; Appendix A: 1 table; Appendix B: 34 Tables; Appendix C: 2 Figures

Sussing Merger Trees: A proposed Merger Tree data format

Authors: Peter A. Thomas, Julian Onions, Dylan Tweed, Andrew J. Benson, Darren Croton, Pascal Elahi, Bruno Henriques, Ilian T. Iliev, Alexander Knebe, Hanni Lux, Yao-Yuan Mao, Mark Neyrinck, Frazer R. Pearce, Vicente Rodriguez-Gomez, Aurel Schneider, Chaichalit Srisawat

Abstract: We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON. We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: Technical report: a proposed merger tree data format. Not submitted to any journal. Source code available at https://bitbucket.org/ProfPAThomas/mergertree

nIFTy Cosmology: Comparison of Galaxy Formation Models

Authors: Alexander Knebe, Frazer R. Pearce, Peter A. Thomas, Andrew Benson, Jeremy Blaizot, Richard Bower, Jorge Carretero, Francisco J. Castander, Andrea Cattaneo, Sofia A. Cora, Darren J. Croton, Weiguang Cui, Daniel Cunnama, Gabriella De Lucia, Julien E. Devriendt, Pascal J. Elahi, Andreea Font, Fabio Fontanot, Juan Garcia-Bellido, Ignacio D. Gargiulo, Violeta Gonzalez-Perez, John Helly, Bruno Henriques, Michaela Hirschmann, Jaehyun Lee , et al. (11 additional authors not shown)

Abstract: We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data… ▽ More We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here. △ Less

Submitted 18 May, 2015; originally announced May 2015.

Comments: 35 pages, 20 figures, accepted for publication in MNRAS

Galaxy formation in the Planck cosmology II. Star formation histories and post-processing magnitude reconstruction

Authors: Sorour Shamshiri, Peter A. Thomas, Bruno M Henriques, Rita Tojeiro, Gerard Lemson, Seb J. Oliver, Stephen Wilkins

Abstract: We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied… ▽ More We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied to data stored in the Millennium data repository. We contrast model SFHs from the Millennium Simulation with observed ones from the VESPA algorithm as applied to the SDSS-7 catalogue. The overall agreement is good, with both simulated and SDSS galaxies showing a steeper SFH with increased stellar mass. The SFHs of blue and red galaxies, however, show poor agreement between data and simulations, which may indicate that the termination of star formation is too abrupt in the models. The mean star-formation rate (SFR) of model galaxies is well-defined and is accurately modelled by a double power law at all redshifts: SFR proportional to $1/(x^{-1.39}+x^{1.33})$, where $x=(t_a-t)/3.0\,$Gyr, $t$ is the age of the stars and $t_a$ is the loopback time to the onset of galaxy formation; above a redshift of unity, this is well approximated by a gamma function: SFR proportional to $x^{1.5}e^{-x}$, where $x=(t_a-t)/2.0\,$Gyr. Individual galaxies, however, show a wide dispersion about this mean. When split by mass, the SFR peaks earlier for high-mass galaxies than for lower-mass ones, and we interpret this downsizing as a mass-dependence in the evolution of the quenched fraction: the SFHs of star-forming galaxies show only a weak mass dependence. △ Less

Submitted 28 September, 2016; v1 submitted 22 January, 2015; originally announced January 2015.

Comments: Accepted version of the paper, to appear in MNRAS. Compared to the original version, contains more detail on the post-processing of magnitudes, including a table of rms magnitude errors. SFHs available on Millennium database http://gavo.mpa-garching.mpg.de/MyMillennium/

Sussing Merger Trees : The Impact of Halo Merger Trees on Galaxy Properties in a Semi-Analytic Model

Authors: Jaehyun Lee, Sukyoung K. Yi, Pascal J. Elahi, Peter A. Thomas, Frazer R. Pearce, Peter Behroozi, Jiaxin Han, John Helly, Intae Jung, Alexander Knebe, Yao-Yuan Mao, Julian Onions, Vicente Rodriguez-Gomez, Aurel Schneider, Chaichalit Srisawat, Dylan Tweed

Abstract: A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for m… ▽ More A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant. △ Less

Submitted 5 October, 2014; originally announced October 2014.

Comments: 15 pages, 18 figures, accepted for publication in MNRAS

Cosmological simulations of galaxy clusters with feedback from active galactic nuclei: profiles and scaling relations

Authors: Simon R. Pike, Scott T. Kay, Richard D. A. Newton, Peter A. Thomas, Adrian Jenkins

Abstract: We present results from a new set of 30 cosmological simulations of galaxy clusters, including the effects of radiative cooling, star formation, supernova feedback, black hole growth and AGN feedback. We first demonstrate that our AGN model is capable of reproducing the observed cluster pressure profile at redshift, z~0, once the AGN heating temperature of the targeted particles is made to scale w… ▽ More We present results from a new set of 30 cosmological simulations of galaxy clusters, including the effects of radiative cooling, star formation, supernova feedback, black hole growth and AGN feedback. We first demonstrate that our AGN model is capable of reproducing the observed cluster pressure profile at redshift, z~0, once the AGN heating temperature of the targeted particles is made to scale with the final virial temperature of the halo. This allows the ejected gas to reach larger radii in higher-mass clusters than would be possible had a fixed heating temperature been used. Such a model also successfully reduces the star formation rate in brightest cluster galaxies and broadly reproduces a number of other observational properties at low redshift, including baryon, gas and star fractions; entropy profiles outside the core; and the X-ray luminosity-mass relation. Our results are consistent with the notion that the excess entropy is generated via selective removal of the densest material through radiative cooling; supernova and AGN feedback largely serve as regulation mechanisms, moving heated gas out of galaxies and away from cluster cores. However, our simulations fail to address a number of serious issues; for example, they are incapable of reproducing the shape and diversity of the observed entropy profiles within the core region. We also show that the stellar and black hole masses are sensitive to numerical resolution, particularly the gravitational softening length; a smaller value leads to more efficient black hole growth at early times and a smaller central galaxy. △ Less

Submitted 2 September, 2014; originally announced September 2014.

Comments: 26 pages, 17 figures. Accepted for publication in MNRAS

The EAGLE project: Simulating the evolution and assembly of galaxies and their environments

Authors: Joop Schaye, Robert A. Crain, Richard G. Bower, Michelle Furlong, Matthieu Schaller, Tom Theuns, Claudio Dalla Vecchia, Carlos S. Frenk, I. G. McCarthy, John C. Helly, Adrian Jenkins, Y. M. Rosas-Guevara, Simon D. M. White, Maarten Baes, C. M. Booth, Peter Camps, Julio F. Navarro, Yan Qu, Alireza Rahmati, Till Sawala, Peter A. Thomas, James Trayford

Abstract: We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massiv… ▽ More We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation. △ Less

Submitted 2 October, 2014; v1 submitted 25 July, 2014; originally announced July 2014.

Comments: Accepted for publication in MNRAS. V2: Minor changes. For images and videos, see http://eagle.strw.leidenuniv.nl/ and http://icc.dur.ac.uk/Eagle/

Journal ref: 2015, MNRAS, 446, 521

Sussing Merger Trees: the influence of the halo finder

Authors: Santiago Avila, Alexander Knebe, Frazer R. Pearce, Aurel Schneider, Chaichalit Srisawat, Peter A. Thomas, Peter Behroozi, Pascal J. Elahi, Jiaxin Han, Yao-Yuan Mao, Julian Onions, Vicente Rodriguez-Gomez

Abstract: Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the construct… ▽ More Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced. △ Less

Submitted 28 April, 2014; v1 submitted 11 February, 2014; originally announced February 2014.

Comments: Accepted in MNRAS. 15 pages, 10 figures, 1 table

Impact of baryons on the cluster mass function and cosmological parameter determination

Authors: Sam J. Cusworth, Scott T. Kay, Richard A. Battye, Peter A. Thomas

Abstract: Recent results by the Planck collaboration have shown that cosmological parameters derived from the cosmic microwave background anisotropies and cluster number counts are in tension, with the latter preferring lower values of the matter density parameter, $Ω_\mathrm{m}$, and power spectrum amplitude, $σ_8$. Motivated by this, we investigate the extent to which the tension may be ameliorated once t… ▽ More Recent results by the Planck collaboration have shown that cosmological parameters derived from the cosmic microwave background anisotropies and cluster number counts are in tension, with the latter preferring lower values of the matter density parameter, $Ω_\mathrm{m}$, and power spectrum amplitude, $σ_8$. Motivated by this, we investigate the extent to which the tension may be ameliorated once the effect of baryonic depletion on the cluster mass function is taken into account. We use the large-volume Millennium Gas simulations in our study, including one where the gas is pre-heated at high redshift and one where the gas is heated by stars and active galactic nuclei (in the latter, the self-gravity of the baryons and radiative cooling are omitted). In both cases, the cluster baryon fractions are in reasonably good agreement with the data at low redshift, showing significant depletion of baryons with respect to the cosmic mean. As a result, it is found that the cluster abundance in these simulations is around 15 per cent lower than the commonly-adopted fit to dark matter simulations by Tinker et al (2008) for the mass range $10^{14}-10^{14.5}h^{-1} \mathrm{M}_\odot$. Ignoring this effect produces a significant artificial shift in cosmological parameters which can be expressed as $Δ[σ_8(Ω_\mathrm{m}/0.27)^{0.38}]\simeq -0.03$ at $z=0.17$ (the median redshift of the $\mathit{Planck}$ cluster sample) for the feedback model. While this shift is not sufficient to fully explain the $\mathit{Planck}$ discrepancy, it is clear that such an effect cannot be ignored in future precision measurements of cosmological parameters with clusters. Finally, we outline a simple, model-independent procedure that attempts to correct for the effect of baryonic depletion and show that it works if the baryon-dark matter back-reaction is negligible. △ Less

Submitted 15 January, 2014; v1 submitted 16 September, 2013; originally announced September 2013.

Comments: 10 pages, 5 figures, Accepted by MNRAS

Sussing Merger Trees: The Merger Trees Comparison Project

Authors: Chaichalit Srisawat, Alexander Knebe, Frazer R. Pearce, Aurel Schneider, Peter A. Thomas, Peter Behroozi, Klaus Dolag, Pascal J. Elahi, Jiaxin Han, John Helly, Yipeng Jing, Intae Jung, Jaehyun Lee, Yao Yuan Mao, Julian Onions, Vicente Rodriguez-Gomez, Dylan Tweed, Sukyoung K. Yi

Abstract: Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were appli… ▽ More Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results. △ Less

Submitted 10 October, 2013; v1 submitted 12 July, 2013; originally announced July 2013.

Comments: 14 pages,11 figures

Modelling Element Abundances in Semi-analytic Models of Galaxy Formation

Authors: Robert M. Yates, Bruno Henriques, Peter A. Thomas, Guinevere Kauffmann, Jonas Johansson, Simon D. M. White

Abstract: We update the treatment of chemical evolution in the Munich semi-analytic model, L-GALAXIES. Our new implementation includes delayed enrichment from stellar winds, supernovae type II (SNe-II) and supernovae type Ia (SNe-Ia), as well as metallicity-dependent yields and a reformulation of the associated supernova feedback. Two different sets of SN-II yields and three different SN-Ia delay-time distr… ▽ More We update the treatment of chemical evolution in the Munich semi-analytic model, L-GALAXIES. Our new implementation includes delayed enrichment from stellar winds, supernovae type II (SNe-II) and supernovae type Ia (SNe-Ia), as well as metallicity-dependent yields and a reformulation of the associated supernova feedback. Two different sets of SN-II yields and three different SN-Ia delay-time distributions (DTDs) are considered, and eleven heavy elements (including O, Mg and Fe) are self-consistently tracked. We compare the results of this new implementation with data on a) local, star-forming galaxies, b) Milky Way disc G dwarfs, and c) local, elliptical galaxies. We find that the z=0 gas-phase mass-metallicity relation is very well reproduced for all forms of DTD considered, as is the [Fe/H] distribution in the Milky Way disc. The [O/Fe] distribution in the Milky Way disc is best reproduced when using a DTD with less than or equal to 50 per cent of SNe-Ia exploding within ~400 Myrs. Positive slopes in the mass-[alpha/Fe] relations of local ellipticals are also obtained when using a DTD with such a minor `prompt' component. Alternatively, metal-rich winds that drive light alpha elements directly out into the circumgalactic medium also produce positive slopes for all forms of DTD and SN-II yields considered. Overall, we find that the best model for matching the wide range of observational data considered here should include a power-law SN-Ia DTD, SN-II yields that take account of prior mass loss through stellar winds, and some direct ejection of light alpha elements out of galaxies. △ Less

Submitted 30 May, 2013; originally announced May 2013.

Comments: 24 pages, 22 figures. Submitted to MNRAS

'Digital' Electron Diffraction - Seeing the Whole Picture

Authors: Richard Beanland, Paul J. Thomas, David I. Woodward, Pamela A. Thomas, Rudolf A. Roemer

Abstract: The advantages of convergent beam electron diffraction for symmetry determination at the scale of a few nm are well known. In practice, the approach is often limited due to the restriction on the angular range of the electron beam imposed by the small Bragg angle for high energy electron diffraction, i.e. a large convergence angle of the incident beam results in overlapping information in the diff… ▽ More The advantages of convergent beam electron diffraction for symmetry determination at the scale of a few nm are well known. In practice, the approach is often limited due to the restriction on the angular range of the electron beam imposed by the small Bragg angle for high energy electron diffraction, i.e. a large convergence angle of the incident beam results in overlapping information in the diffraction pattern. Techniques have been generally available since the 1980s which overcome this restriction for individual diffracted beams, by making a compromise between illuminated area and beam convergence. Here, we describe a simple technique which overcomes all of these problems using computer control, giving electron diffraction data over a large angular range for many diffracted beams from the volume given by a focused electron beam (typically a few nm or less). The increase in the amount of information significantly improves ease of interpretation and widens the applicability of the technique, particularly for thin materials or those with larger lattice parameters. △ Less

Submitted 28 November, 2012; originally announced November 2012.

Comments: 13 pages, 8 figures

Journal ref: Acta Cryst. (2013). A69, 427-434

The inhomogenous reionization of the inter-galactic medium by metal-poor globular clusters

Authors: B. F. Griffen, M. J. Drinkwater, Ilian T. Iliev, P. A. Thomas, Garrelt Mellema

Abstract: We present detailed radiative transfer simulations of the reionization history of the Milky Way by metal-poor globular clusters. We identify potential metal-poor globular cluster candidates within the Aquarius simulation using dark matter halo velocity dispersions. We calculate the local ionization fields via a photon-conserving, three dimensional non-equilibrium chemistry code and allow the model… ▽ More We present detailed radiative transfer simulations of the reionization history of the Milky Way by metal-poor globular clusters. We identify potential metal-poor globular cluster candidates within the Aquarius simulation using dark matter halo velocity dispersions. We calculate the local ionization fields via a photon-conserving, three dimensional non-equilibrium chemistry code and allow the model to propagate through to the present day. The key feature of the model is that globular cluster formation is suppressed if the local gas is ionized. We find that our spatial treatment of the ionization field leads to drastically different numbers and spatial distributions when compared to models where globular cluster formation is simply truncated at a given redshift. We find that it is possible for metal-poor globular clusters to have formed via the dark matter halo formation channel as our secondary model (delayed formation) combined with truncation at z = 10 produces radial distributions statistically consistent with that of the Milky Way metal-poor globular clusters. If globular clusters do indeed form within high-redshift dark matter halos, if only in-part, their contributions to the reionization of the local (i.e. 2^3 h^-3 Mpc^3 centred on the host galaxy) volume and mass by redshift 10 could be as high as 98% and 90%, respectively. In our photon poorest model, this contribution drops to 60% and 50%. The surviving clusters in all models have a narrow average age range (mean = 13.34 Gyr, σ= 0.04 Gyr) consistent with current ages estimates of the Milky Way metal-poor globular clusters. We also test a simple dynamical destruction model and estimate that ~60% of all metal-poor globular clusters formed at high redshift have since been destroyed via tidal interactions with the host galaxy. △ Less

Submitted 26 September, 2012; originally announced September 2012.

Comments: 18 pages, 15 figures, 5 tables, submitted to MNRAS

The Spitzer Extragalactic Representative Volume Survey (SERVS): survey definition and goals

Authors: J. -C. Mauduit, M. Lacy, D. Farrah, J. A. Surace, M. Jarvis, S. Oliver, C. Maraston, M. Vaccari, L. Marchetti, G. Zeimann, E. A. Gonzalez-Solares, J. Pforr, A. O. Petric, B. Henriques, P. A. Thomas, J. Afonso, A. Rettura, G. Wilson, J. T. Falder, J. E. Geach, M. Huynh, R. P. Norris, N. Seymour, G. T. Richards, S. A. Stanford , et al. (59 additional authors not shown)

Abstract: We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a fun… ▽ More We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey. △ Less

Submitted 17 September, 2012; v1 submitted 18 June, 2012; originally announced June 2012.

Comments: 21 pages, 23 figures, published in PASP. This version differs slightly from PASP, mainly due to formatting issues. Figure 2 was also corrected from the previous version

Journal ref: Publications of the Astronomical Society of the Pacific, 2012, Volume 124, issue 917, pp.714-736

The XMM Cluster Survey: Evidence for energy injection at high redshift from evolution of the X-ray luminosity-temperature relation

Authors: Matt Hilton, A. Kathy Romer, Scott T. Kay, Nicola Mehrtens, E. J. Lloyd-Davies, Peter A. Thomas, Chris J. Short, Julian A. Mayers, Philip J. Rooney, John P. Stott, Chris A. Collins, Craig D. Harrison, Ben Hoyle, Andrew R. Liddle, Robert G. Mann, Christopher J. Miller, Martin Sahlen, Pedro T. P. Viana, Michael Davidson, Mark Hosmer, Robert C. Nichol, Kivanc Sabirli, S. A. Stanford, Michael J. West

Abstract: We measure the evolution of the X-ray luminosity-temperature (L_X-T) relation since z~1.5 using a sample of 211 serendipitously detected galaxy clusters with spectroscopic redshifts drawn from the XMM Cluster Survey first data release (XCS-DR1). This is the first study spanning this redshift range using a single, large, homogeneous cluster sample. Using an orthogonal regression technique, we find… ▽ More We measure the evolution of the X-ray luminosity-temperature (L_X-T) relation since z~1.5 using a sample of 211 serendipitously detected galaxy clusters with spectroscopic redshifts drawn from the XMM Cluster Survey first data release (XCS-DR1). This is the first study spanning this redshift range using a single, large, homogeneous cluster sample. Using an orthogonal regression technique, we find no evidence for evolution in the slope or intrinsic scatter of the relation since z~1.5, finding both to be consistent with previous measurements at z~0.1. However, the normalisation is seen to evolve negatively with respect to the self-similar expectation: we find E(z)^{-1} L_X = 10^{44.67 +/- 0.09} (T/5)^{3.04 +/- 0.16} (1+z)^{-1.5 +/- 0.5}, which is within 2 sigma of the zero evolution case. We see milder, but still negative, evolution with respect to self-similar when using a bisector regression technique. We compare our results to numerical simulations, where we fit simulated cluster samples using the same methods used on the XCS data. Our data favour models in which the majority of the excess entropy required to explain the slope of the L_X-T relation is injected at high redshift. Simulations in which AGN feedback is implemented using prescriptions from current semi-analytic galaxy formation models predict positive evolution of the normalisation, and differ from our data at more than 5 sigma. This suggests that more efficient feedback at high redshift may be needed in these models. △ Less

Submitted 4 July, 2012; v1 submitted 24 May, 2012; originally announced May 2012.

Comments: Accepted for publication in MNRAS; 12 pages, 6 figures; added references to match published version