Showing 1–50 of 88 results for author: Croton, D J

Using the Star Forming Main Sequence To Explore Quiescent Galaxies Across Cosmic Time

Authors: Tyler Houston, Darren J. Croton, Manodeep Sinha

Abstract: In this letter, we explore the quiescent lives of central galaxies using the SAGE galaxy model and Uchuu dark matter simulation. We ask three questions: (1) How much of a galaxy's life is spent in quiescence? (2) How often do galaxies transition off the main sequence? (3) What is the typical duration of a quiescent phase? We find low and high-mass galaxies spend the highest fraction of their lives… ▽ More In this letter, we explore the quiescent lives of central galaxies using the SAGE galaxy model and Uchuu dark matter simulation. We ask three questions: (1) How much of a galaxy's life is spent in quiescence? (2) How often do galaxies transition off the main sequence? (3) What is the typical duration of a quiescent phase? We find low and high-mass galaxies spend the highest fraction of their lives in quiescence: 45 \pm 19% for log10(Mstar) < 9.0 (3.68 \pm 1.80 Gyr) and 26 \pm 25% for log10(Mstar) > 11.5 (3.46 \pm 3.30 Gyr), falling to 7 \pm 13% for galaxies in-between (0.82 \pm 1.57 Gyr). Low mass galaxies move in and out of quiescence frequently, 2.8 \pm 1.3 times on average, though only for short periods, 1.49 \pm 1.04 Gyr. This can be traced to the influence of supernova feedback on their quite stochastic evolution. Galaxies of higher mass have fewer quiescent periods, ~0.7 \pm 0.9, and their length increases with mass, peaking at 1.97 \pm 2.27 Gyr. However, our high-mass population comprises star-forming and quiescent galaxies with diverging evolutionary paths, so the actual time may be even longer. These high-mass trends are driven by radio mode feedback from supermassive black holes, which, once active, tend to remain active for extended periods. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 5 pages, 3 figures, accepted to MNRAS Letters, open access version see https://doi.org/10.1093/mnrasl/slad031

Dust contribution to the panchromatic galaxy emission

Authors: Dian P. Triani, Darren J. Croton, Manodeep Sinha, Edward N. Taylor, Camilla Pacifici, Eli Dwek

Abstract: We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically… ▽ More We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically motivated attenuation model from the computed dust properties in \dustysage, so each galaxy has a self-consistent set of attenuation parameters based on the complicated dust physics that occurred across the galaxy's assembly history. Then, we explore several dust emission templates to produce infrared spectra. Our results show that a physically-motivated attenuation model is better for obtaining a consistent multi-wavelength description of galaxy formation and evolution, compared to using a constant attenuation. We compare our predictions with a compilation of observations and find that the fiducial model is in reasonable agreement with: (i) the observed $z=0$ luminosity functions from the far-ultraviolet to far-infrared simultaneously, and hence (ii) the local cosmic SED in the same range, (iii) the rest-frame K-band luminosity function across $0 < z < 3$, and (iv) the rest-frame far-ultraviolet luminosity function across $0 < z < 1$. Our model underproduces the far-ultraviolet emission at $z=2$ and $z=3$, which can be improved by altering the AGN feedback and dust processes in \dustysage. However, this combination thus worses the agreement at $z=0$, which suggests that more detailed treatment of such processes is required. △ Less

Submitted 12 December, 2022; originally announced December 2022.

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

Coupling Semi-Analytic and N-body Galaxies (CoSANG) for cosmological stellar halo simulations I- Methods and the structure of dark matter halos

Authors: Shahram Talei, Krista Mccord, Jeremy Bailin, Darren J. Croton, Alexandra Mannings, Michael Sitarz, Annelia Anderson, Brooke Bailey, Alexis Rollins

Abstract: We present CoSANG (Coupling Semi-Analytic and N-body Galaxies), a new hybrid model for cosmological dark matter and stellar halo simulations. In this approach a collisionless model (Gadget3) is used for gravitational interactions while a coupled semi-analytic model (SAGE) calculates baryonic effects at each time-step. This live self-consistent interaction at each time-step is the key difference be… ▽ More We present CoSANG (Coupling Semi-Analytic and N-body Galaxies), a new hybrid model for cosmological dark matter and stellar halo simulations. In this approach a collisionless model (Gadget3) is used for gravitational interactions while a coupled semi-analytic model (SAGE) calculates baryonic effects at each time-step. This live self-consistent interaction at each time-step is the key difference between CoSANG and traditional semi-analytic models that are mainly used for post-processing. By accounting for the gravitational effect of the baryons, CoSANG can overcome some of the deficiencies of pure N-body simulations, while being less computationally expensive than hydrodynamic simulations. Moreover CoSANG can produce stellar halo populations via tagging tracer dark matter particles. We demonstrate the performance and dynamical accuracy of this approach using both controlled test simulations and a set of three cosmological zoom-in simulations of Milky Way (MW) mass halos. We simulate each target halo both without the coupling (dark matter only, hereafter DMO) and with the coupling (CoSANG). We compare the internal structure (subhalo distribution, shape and orientation) of the halos. The following changes are observed in the CoSANG model compared to the DMO model: 1) the total number of subhaloes close to the center of the halo is reduced, 2) the $V_{\mathrm{max}}$ distribution peaks at a lower value and lies below the DMO model, 3) the axis ratio is smaller. The difference between DMO and CoSANG simulations is more significant in early forming halos. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 24 pages, 15 figures, submitted to MNRAS, comments are welcome

Investigating the Effect of Galaxy Interactions on Star Formation at 0.5<z<3.0

Authors: Ekta A. Shah, Jeyhan S. Kartaltepe, Christina T. Magagnoli, Isabella G. Cox, Caleb T. Wetherell, Brittany N. Vanderhoof, Kevin C. Cooke, Antonello Calabro, Nima Chartab, Christopher J. Conselice, Darren J. Croton, Alexander de la Vega, Nimish P. Hathi, Olivier Ilbert, Hanae Inami, Dale D. Kocevski, Anton M. Koekemoer, Brian C. Lemaux, Lori Lubin, Kameswara Bharadwaj Mantha, Stefano Marchesi, Marie Martig, Jorge Moreno, Belen Alcalde Pampliega, David R. Patton , et al. (2 additional authors not shown)

Abstract: Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increa… ▽ More Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increased SFRs in these galaxies. To test this, we measure the SFR enhancement in a total of 2351 (1327) massive ($M_*>10^{10}M_\odot$) major ($1<M_1/M_2<4$) spectroscopic galaxy pairs at 0.5<z<3.0 with $ΔV <5000$ km s$^{-1}$ (1000 km s$^{-1}$) and projected separation <150 kpc selected from the extensive spectroscopic coverage in the COSMOS and CANDELS fields. We find that the highest level of SFR enhancement is a factor of 1.23$^{+0.08}_{-0.09}$ in the closest projected separation bin (<25 kpc) relative to a stellar mass-, redshift-, and environment-matched control sample of isolated galaxies. We find that the level of SFR enhancement is a factor of $\sim1.5$ higher at 0.5<z<1 than at 1<z<3 in the closest projected separation bin. Among a sample of visually identified mergers, we find an enhancement of a factor of 1.86$^{+0.29}_{-0.18}$ for coalesced systems. For this visually identified sample, we see a clear trend of increased SFR enhancement with decreasing projected separation (2.40$^{+0.62}_{-0.37}$ vs.\ 1.58$^{+0.29}_{-0.20}$ for 0.5<z<1.6 and 1.6<z<3.0, respectively). The SFR enhancement seen in our interactions and mergers are all lower than the level seen in local samples at the same separation, suggesting that the level of interaction-induced star formation evolves significantly over this time period. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 23 pages, 13 figures, Accepted for publication in ApJ

Exploring the relation between dust mass and galaxy properties using Dusty SAGE

Authors: Dian P. Triani, Manodeep Sinha, Darren J. Croton, Eli Dwek, Camilla Pacifici

Abstract: We explore the relation between dust and several fundamental properties of simulated galaxies using the Dusty SAGE semi-analytic model. In addition to tracing the standard galaxy properties, Dusty SAGE also tracks cold dust mass in the interstellar medium (ISM), hot dust mass in the halo and dust mass ejected by feedback activity. Based on their ISM dust content, we divide our galaxies into two ca… ▽ More We explore the relation between dust and several fundamental properties of simulated galaxies using the Dusty SAGE semi-analytic model. In addition to tracing the standard galaxy properties, Dusty SAGE also tracks cold dust mass in the interstellar medium (ISM), hot dust mass in the halo and dust mass ejected by feedback activity. Based on their ISM dust content, we divide our galaxies into two categories: ISM dust-poor and ISM dust-rich. We split the ISM dust-poor group into two subgroups: halo dust-rich and dust-poor (the latter contains galaxies that lack dust in both the ISM and halo). Halo dust-rich galaxies have high outflow rates of heated gas and dust and are more massive. We divide ISM dust-rich galaxies based on their specific star formation rate (sSFR) into star-forming and quenched subgroups. At redshift z=0, we find that ISM dust-rich galaxies have a relatively high sSFR, low bulge-to-total (BTT) mass ratio, and high gas metallicity. The high sSFR of ISM dust-rich galaxies allows them to produce dust in the stellar ejecta. Their metal-rich ISM enables dust growth via grain accretion. The opposite is seen in the ISM dust-poor group. Furthermore, ISM dust-rich galaxies are typically late-types, while ISM dust-poor galaxies resemble the early-type population, and we show how their ISM content evolves from being dust-rich to dust-poor. Finally, we investigate dust production from z=3 to z=0 and find that all groups evolve similarly, except for the quenched ISM dust-rich group. △ Less

Submitted 25 February, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 12 pages, 8 figure, accepted for publication in MNRAS

Investigating the Effect of Galaxy Interactions on AGN Enhancement at $0.5<z<3.0$

Authors: Ekta A. Shah, Jeyhan S. Kartaltepe, Christina T. Magagnoli, Isabella G. Cox, Caleb T. Wetherell, Brittany N. Vanderhoof, Antonello Calabro, Nima Chartab, Christopher J. Conselice, Darren J. Croton, Jennifer Donley, Laura de Groot, Alexander de la Vega, Nimish P. Hathi, Olivier Ilbert, Hanae Inami, Dale D. Kocevski, Anton M. Koekemoer, Brian C. Lemaux, Kameswara Bharadwaj Mantha, Stefano Marchesi, Marie Martig, Daniel C. Masters, Elizabeth J. McGrath, Daniel H. McIntosh , et al. (8 additional authors not shown)

Abstract: Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the large… ▽ More Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the largest known sample of major spectroscopic galaxy pairs (2381 with $ΔV <5000$ km s$^{-1}$) at $0.5<z<3.0$ from observations in the COSMOS and CANDELS surveys. We identify X-ray and IR AGN among this kinematic pair sample, a visually identified sample of mergers and interactions, and a mass-, redshift-, and environment-matched control sample for each in order to calculate AGN fractions and the level of AGN enhancement as a function of relative velocity, redshift, and X-ray luminosity. While we see a slight increase in AGN fraction with decreasing projected separation, overall, we find no significant enhancement relative to the control sample at any separation. In the closest projected separation bin ($<25$ kpc, $ΔV <1000$ km s$^{-1}$), we find enhancements of a factor of 0.94$^{+0.21}_{-0.16}$ and 1.00$^{+0.58}_{-0.31}$ for X-ray and IR-selected AGN, respectively. While we conclude that galaxy interactions do not significantly enhance AGN activity on average over $0.5<z<3.0$ at these separations, given the errors and the small sample size at the closest projected separations, our results would be consistent with the presence of low-level AGN enhancement. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ, 28 pages, 15 figures

The origin of dust in galaxies across cosmic time

Authors: Dian P. Triani, Manodeep Sinha, Darren J. Croton, Camilla Pacifici, Eli Dwek

Abstract: We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM… ▽ More We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM by star formation, reheating, inflows and outflows. Our model successfully reproduces the observed dust mass function at redshift z = 0 and the observed scaling relations for dust across a wide range of redshifts. We find that the dust mass content in the present Universe is mainly produced via grain growth in the interstellar medium (ISM). By contrast, in the early Universe, the primary production mechanism for dust is the condensation in stellar ejecta. The shift of the significant production channel for dust characterises the scaling relations of dust-to-gas (DTG) and dust-to-metal (DTM) ratios. In galaxies where the grain growth dominates, we find positive correlations for DTG and DTM ratios with both metallicity and stellar mass. On the other hand, in galaxies where dust is produced primarily via condensation, we find negative or no correlation for DTM and DTG ratios with either metallicity or stellar mass. In agreement with observation showing that the circumgalactic medium (CGM) contains more dust than the ISM, our model also shows the same trend for z < 4. Our semi-analytic model is publicly available at https: //github.com/dptriani/dusty-sage. △ Less

Submitted 12 February, 2020; originally announced February 2020.

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

Journal ref: 2020MNRAS.tmp..415T

The WiggleZ Dark Energy Survey: Final Data Release and the Metallicity of UV-Luminous Galaxies

Authors: Michael J. Drinkwater, Zachary J. Byrne, Chris Blake, Karl Glazebrook, Sarah Brough, Matthew Colless, Warrick Couch, Darren J. Croton, Scott M. Croom, Tamara M. Davis, Karl Forster, David Gilbank, Samuel R. Hinton, Ben Jelliffe, Russell J. Jurek, I-hui Li, D. Christopher Martin, Kevin Pimbblet, Gregory B. Poole, Michael Pracy, Rob Sharp, Jon Smillie, Max Spolaor, Emily Wisnioski, David Woods , et al. (2 additional authors not shown)

Abstract: The WiggleZ Dark Energy Survey measured the redshifts of over 200,000 UV-selected (NUV<22.8 mag) galaxies on the Anglo-Australian Telescope. The survey detected the baryon acoustic oscillation signal in the large scale distribution of galaxies over the redshift range 0.2<z<1.0, confirming the acceleration of the expansion of the Universe and measuring the rate of structure growth within it. Here w… ▽ More The WiggleZ Dark Energy Survey measured the redshifts of over 200,000 UV-selected (NUV<22.8 mag) galaxies on the Anglo-Australian Telescope. The survey detected the baryon acoustic oscillation signal in the large scale distribution of galaxies over the redshift range 0.2<z<1.0, confirming the acceleration of the expansion of the Universe and measuring the rate of structure growth within it. Here we present the final data release of the survey: a catalogue of 225415 galaxies and individual files of the galaxy spectra. We analyse the emission-line properties of these UV-luminous Lyman-break galaxies by stacking the spectra in bins of luminosity, redshift, and stellar mass. The most luminous (-25 mag < MFUV <-22 mag) galaxies have very broad H-beta emission from active nuclei, as well as a broad second component to the [OIII] (495.9 nm, 500.7 nm) doublet lines that is blue shifted by 100 km/s, indicating the presence of gas outflows in these galaxies. The composite spectra allow us to detect and measure the temperature-sensitive [OIII] (436.3 nm) line and obtain metallicities using the direct method. The metallicities of intermediate stellar mass (8.8<log(M*/Msun)<10) WiggleZ galaxies are consistent with normal emission-line galaxies at the same masses. In contrast, the metallicities of high stellar mass (10<log(M*/Msun)<12) WiggleZ galaxies are significantly lower than for normal emission-line galaxies at the same masses. This is not an effect of evolution as the metallicities do not vary with redshift; it is most likely a property specific to the extremely UV-luminous WiggleZ galaxies. △ Less

Submitted 18 October, 2019; originally announced October 2019.

Comments: Catalogue available at MNRAS (DOI link below) and also at: https://espace.library.uq.edu.au/view/UQ:3e43575

Journal ref: Monthly Notices of the Royal Astronomical Society, 2018, Volume 474, p.4151-4168

The 21cm bispectrum during reionization: a tracer of the ionization topology

Authors: Anne Hutter, Catherine A. Watkinson, Jacob Seiler, Pratika Dayal, Manodeep Sinha, Darren J. Croton

Abstract: We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties an… ▽ More We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array. △ Less

Submitted 18 November, 2019; v1 submitted 9 July, 2019; originally announced July 2019.

Comments: 16 pages, 7 figures. Accepted for publication in MNRAS

Evolution of black hole and galaxy growth in a semi-numerical galaxy formation model

Authors: Mackenzie L. Jones, R. C. Hickox, S. J. Mutch, D. J. Croton, A. F. Ptak, M. A. DiPompeo

Abstract: We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit o… ▽ More We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts. △ Less

Submitted 1 July, 2019; originally announced July 2019.

Comments: 9 pages, 8 figures, Accepted to ApJ

Feedback by supermassive black holes in galaxy evolution: impacts of accretion and outflows on the star formation rate

Authors: Mojtaba Raouf, Joseph Silk, Stanislav S. Shabala, Gary A. Mamon, Darren J. Croton, Habib G. Khosroshahi, Ricarda S. Beckmann

Abstract: We describe a physical model of the outflows produced as a result of gas accretion onto a black hole, and the resultant changes to star formation rates and efficiencies in galaxies, using the Radio-SAGE semi-analytic galaxy formation model. We show that the ratio of outflow rate to SFR of galaxies is mainly driven by black hole mass and virial halo mass, and show that the SFR is higher than the ou… ▽ More We describe a physical model of the outflows produced as a result of gas accretion onto a black hole, and the resultant changes to star formation rates and efficiencies in galaxies, using the Radio-SAGE semi-analytic galaxy formation model. We show that the ratio of outflow rate to SFR of galaxies is mainly driven by black hole mass and virial halo mass, and show that the SFR is higher than the outflow rate at low black hole masses. The model consistently reproduces the observed evolution of star formation rate density from z = 6 to z = 0, as well as the trend of the stellar mass - halo mass relations. We show the characteristic growth of massive galaxies influenced by AGN feedback at different redshifts. We find feedback to be prevalent in the most massive galaxy halos, inhibiting the cooling catastrophe. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 15 pages, 8 figures. Accepted for publication in MNRAS

The Three Hundred project: a large catalogue of theoretically modelled galaxy clusters for cosmological and astrophysical applications

Authors: Weiguang Cui, Alexander Knebe, Gustavo Yepes, Frazer Pearce, Chris Power, Romeel Dave, Alexander Arth, Stefano Borgani, Klaus Dolag, Pascal Elahi, Robert Mostoghiu, Giuseppe Murante, Elena Rasia, Doris Stoppacher, Jesus Vega-Ferrero, Yang Wang, Xiaohu Yang, Andrew Benson, Sofía A. Cora, Darren J. Croton, Manodeep Sinha, Adam R. H. Stevens, Cristian A. Vega-Martínez, Jake Arthur, Anna S. Baldi , et al. (12 additional authors not shown)

Abstract: We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic… ▽ More We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 20 pages, 8 figures, 7 tables. MNRAS published version

Journal ref: Cui, W., Knebe, A., Yepes, G., et al.\ 2018, \mnras, 480, 2898

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

Merger history of central galaxies in Semi-Analytic Models of galaxy formation

Authors: Mojtaba Raouf, Habib G. Khosroshahi, Gary A. Mamon, Darren J. Croton, Abdolhosein Hashemizadeh, Ali A. Dariush

Abstract: We investigate the dynamical evolution of galaxies in groups with different formation epochs. Galaxy groups have been selected to be in different dynamical states, namely dynamically old and dynamically young, which reflect their early and late formation times, respectively, based on their halo mass assembly. Brightest galaxies in dynamically young groups have suffered their last major galaxy merg… ▽ More We investigate the dynamical evolution of galaxies in groups with different formation epochs. Galaxy groups have been selected to be in different dynamical states, namely dynamically old and dynamically young, which reflect their early and late formation times, respectively, based on their halo mass assembly. Brightest galaxies in dynamically young groups have suffered their last major galaxy merger typically $\sim 2$ Gyr more recently than their counterparts in dynamically old groups. Furthermore, we study the evolution of velocity dispersion in these two classes and compare them with the analytic models of isolated halos. The velocity dispersion of dwarf galaxies in high mass, dynamically young groups increases slowly in time, while the analogous dispersion in dynamically old high-mass groups is constant. In contrast, the velocity dispersion of giant galaxies in low mass groups decreases rapidly at late times. This increasing velocity bias is caused by dynamical friction, and starts much earlier in the dynamically old groups. The recent {\sc Radio-SAGE} model of galaxy formation suggests that radio luminosities of central galaxies, considered to be tracers of AGN activity, are enhanced in halos that assembled more recently, independent of the time since the last major merger. △ Less

Submitted 20 June, 2018; v1 submitted 6 March, 2018; originally announced March 2018.

Comments: 11 pages, 7 figures. Accepted for publication in ApJ

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

MultiDark-Galaxies: data release and first results

Authors: Alexander Knebe, Doris Stoppacher, Francisco Prada, Christoph Behrens, Andrew Benson, Sofia A. Cora, Darren J. Croton, Nelson D. Padilla, Andrés N. Ruiz, Manodeep Sinha, Adam R. H. Stevens, Cristian A. Vega-Martínez, Peter Behroozi, Violeta Gonzalez-Perez, Stefan Gottlöber, Anatoly A. Klypin, Gustavo Yepes, Harry Enke, Noam I. Libeskind, Kristin Riebe, Matthias Steinmetz

Abstract: We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data… ▽ More We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available. △ Less

Submitted 23 October, 2017; originally announced October 2017.

Comments: 29 pages, 16 figures, 8 tables, accepted for publication in MNRAS. All data incl. the complete galaxy catalogues for all models are publicly available from the CosmoSim database (http://www.cosmosim.org); a selected set of galaxy properties is available via the Skies & Universes website (http://www.skiesanduniverses.org)

Convergence properties of halo merger trees; halo and substructure merger rates across cosmic history

Authors: Gregory B. Poole, Simon J. Mutch, Darren J. Croton, Stuart Wyithe

Abstract: We introduce gbpTrees: an algorithm for constructing merger trees from cosmological simulations, designed to identify and correct for pathological cases introduced by errors or ambiguities in the halo finding process. gbpTrees is built upon a halo matching method utilising pseudo-radial moments constructed from radially-sorted particle ID lists (no other information is required) and a scheme for c… ▽ More We introduce gbpTrees: an algorithm for constructing merger trees from cosmological simulations, designed to identify and correct for pathological cases introduced by errors or ambiguities in the halo finding process. gbpTrees is built upon a halo matching method utilising pseudo-radial moments constructed from radially-sorted particle ID lists (no other information is required) and a scheme for classifying merger tree pathologies from networks of matches made to-and-from haloes across snapshots ranging forward-and-backward in time. Focusing on Subfind catalogs for this work, a sweep of parameters influencing our merger tree construction yields the optimal snapshot cadence and scanning range required for converged results. Pathologies proliferate when snapshots are spaced by $\lesssim{0.128}$ dynamical times; conveniently similar to that needed for convergence of semi-analytical modelling, as established by Benson etal Total merger counts are converged at the level of $\sim{5}$% for friends-of-friends (FoF) haloes of size $n_{\rm p}\gtrsim{75}$ across a factor of 512 in mass resolution, but substructure rates converge more slowly with mass resolution, reaching convergence of $\sim{10}$% for $n_{\rm p}\gtrsim{100}$ and particle mass $m_{\rm p}{\lesssim}10^{9}M_\odot$. We present analytic fits to FoF and substructure merger rates across nearly all observed galactic history ($z{\le}8.5$). While we find good agreement with the results presented by Fakhouri etal for FoF haloes, a slightly flatter dependance on merger ratio and increased major merger rates are found, reducing previously reported discrepancies with extended Press-Schechter estimates. When appropriately defined, substructure merger rates show a similar mass ratio dependance as FoF rates, but with stronger mass and redshift dependencies for their normalisation. △ Less

Submitted 5 September, 2017; originally announced September 2017.

Comments: Accepted for publication in MNRAS. 25 pages, 17 figures. Movies available at http://www.astronomy.swin.edu.au/~gpoole/gbpTrees.html . Code available at https://github.com/gbpoole/gbpCode

Triggering Active Galactic Nuclei in galaxy clusters

Authors: Madeline A. Marshall, Stanislav S. Shabala, Martin G. H. Krause, Kevin A. Pimbblet, Darren J. Croton, Matt S. Owers

Abstract: We model the triggering of Active Galactic Nuclei (AGN) in galaxy clusters using the semi- analytic galaxy formation model SAGE (?). We prescribe triggering methods based on the ram pressure galaxies experience as they move throughout the intracluster medium, which is hypothesized to trigger star formation and AGN activity. The clustercentric radius and velocity distribution of the simulated activ… ▽ More We model the triggering of Active Galactic Nuclei (AGN) in galaxy clusters using the semi- analytic galaxy formation model SAGE (?). We prescribe triggering methods based on the ram pressure galaxies experience as they move throughout the intracluster medium, which is hypothesized to trigger star formation and AGN activity. The clustercentric radius and velocity distribution of the simulated active galaxies produced by these models are compared with that of AGN and galaxies with intense star formation from a sample of low-redshift, relaxed clusters from the Sloan Digital Sky Survey. The ram pressure triggering model that best explains the clustercentric radius and velocity distribution of these observed galaxies has AGN and star formation triggered if $2.5\times10^{-14} < P_{ram} < 2.5\times10^{-13}$ Pa and $P_{ram} > 2P_{internal}$; this is consistent with expectations from hydrodynamical simulations of ram-pressure induced star formation. Our results show that ram pressure is likely to be an important mechanism for triggering star formation and AGN activity in clusters. △ Less

Submitted 20 November, 2017; v1 submitted 18 August, 2017; originally announced August 2017.

Comments: 15 pages, 12 figures. Accepted for publication in MNRAS

The many lives of active galactic nuclei-II: The formation and evolution of radio jets and their impact on galaxy evolution

Authors: Mojtaba Raouf, Stanislav S. Shabala, Darren J. Croton, Habib G. Khosroshahi, Maksym Bernyk

Abstract: We describe new efforts to model radio active galactic nuclei (AGN) in a cosmological context using the SAGE semi-analytic galaxy model. Our new method tracks the physical properties of radio jets in massive galaxies, including the evolution of radio lobes and their impact on the surrounding gas. This model also self consistently follows the gas cooling-heating cycle that significantly shapes star… ▽ More We describe new efforts to model radio active galactic nuclei (AGN) in a cosmological context using the SAGE semi-analytic galaxy model. Our new method tracks the physical properties of radio jets in massive galaxies, including the evolution of radio lobes and their impact on the surrounding gas. This model also self consistently follows the gas cooling-heating cycle that significantly shapes star formation and the life and death of many galaxy types. Adding jet physics to SAGE adds new physical properties to the model output, which in turn allows us to make more detailed predictions for the radio AGN population. After calibrating the model to a set of core observations we analyse predictions for jet power, radio cocoon size, radio luminosity, and stellar mass. We find that the model is able to match the stellar mass--radio luminosity relation at $z\sim0$, and the radio luminosity function out to $z\sim1$. This updated model will make possible the construction of customised AGN-focused mock survey catalogues to be used for large-scale observing programs. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: 15 pages, 12 figures. Accepted for publication in MNRAS

Effect of local environment and stellar mass on galaxy quenching and morphology at $0.5<z<2.0$

Authors: Lalitwadee Kawinwanichakij, Casey Papovich, Ryan F. Quadri, Karl Glazebrook, Glenn G. Kacprzak, Rebecca J. Allen, Eric F. Bell, Darren J. Croton, Avishai Dekel, Henry C. Ferguson, Ben Forrest, Norman A. Grogin, Yicheng Guo, Dale D. Kocevski, Anton M. Koekemoer, Ivo Labbé, Ray A. Lucas, Themiya Nanayakkara, Lee R. Spitler, Caroline M. S. Straatman, Kim-Vy H. Tran, Adam Tomczak, Pieter van Dokkum

Abstract: We study galactic star-formation activity as a function of environment and stellar mass over 0.5<z<2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey. We estimate the galaxy environment using a Bayesian-motivated measure of the distance to the third nearest neighbor for galaxies to the stellar mass completeness of our survey, $\log(M/M_\odot)>9 (9.5)$ at z=1.3 (2.0). This method, when applie… ▽ More We study galactic star-formation activity as a function of environment and stellar mass over 0.5<z<2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey. We estimate the galaxy environment using a Bayesian-motivated measure of the distance to the third nearest neighbor for galaxies to the stellar mass completeness of our survey, $\log(M/M_\odot)>9 (9.5)$ at z=1.3 (2.0). This method, when applied to a mock catalog with the photometric-redshift precision ($σ_z / (1+z) \lesssim 0.02$), recovers galaxies in low- and high-density environments accurately. We quantify the environmental quenching efficiency, and show that at z> 0.5 it depends on galaxy stellar mass, demonstrating that the effects of quenching related to (stellar) mass and environment are not separable. In high-density environments, the mass and environmental quenching efficiencies are comparable for massive galaxies ($\log (M/M_\odot)\gtrsim$ 10.5) at all redshifts. For lower mass galaxies ($\log (M/M)_\odot) \lesssim$ 10), the environmental quenching efficiency is very low at $z\gtrsim$ 1.5, but increases rapidly with decreasing redshift. Environmental quenching can account for nearly all quiescent lower mass galaxies ($\log(M/M_\odot) \sim$ 9-10), which appear primarily at $z\lesssim$ 1.0. The morphologies of lower mass quiescent galaxies are inconsistent with those expected of recently quenched star-forming galaxies. Some environmental process must transform the morphologies on similar timescales as the environmental quenching itself. The evolution of the environmental quenching favors models that combine gas starvation (as galaxies become satellites) with gas exhaustion through star-formation and outflows ("overconsumption"), and additional processes such as galaxy interactions, tidal stripping and disk fading to account for the morphological differences between the quiescent and star-forming galaxy populations. △ Less

Submitted 20 September, 2017; v1 submitted 12 June, 2017; originally announced June 2017.

Comments: 29 pages, 15 figure, accepted for publication in ApJ

The intrinsic characteristics of galaxies on the SFR-stellar mass plane at 1.2<z<4: I. the correlation between stellar age, central density and position relative to the main sequence

Authors: Bomee Lee, Mauro Giavalisco, Katherine Whitaker, Christina C. Williams, Henry C. Ferguson, Viviana Acquaviva, Anton M. Koekemoer, Amber N. Straughn, Yicheng Guo, Jeyhan S. Kartaltepe, Jennifer Lotz, Camilla Pacifici, Darren J. Croton, Rachel S. Somerville, Yu Lu

Abstract: We use the deep CANDELS observations in the GOODS North and South fields to revisit the correlations between stellar mass ($M_*$), star--formation rate (SFR) and morphology, and to introduce a fourth dimension, the mass-weighted stellar age, in galaxies at $1.2<z<4$. We do this by making new measures of $M_*$, $SFR$, and stellar age thanks to an improved SED fitting procedure that allows various s… ▽ More We use the deep CANDELS observations in the GOODS North and South fields to revisit the correlations between stellar mass ($M_*$), star--formation rate (SFR) and morphology, and to introduce a fourth dimension, the mass-weighted stellar age, in galaxies at $1.2<z<4$. We do this by making new measures of $M_*$, $SFR$, and stellar age thanks to an improved SED fitting procedure that allows various star formation history for each galaxy. Like others, we find that the slope of the Main Sequence (MS) of star formation in the $(M_*;SFR)$ plane bends at high mass. We observe clear morphological differences among galaxies across the MS, which also correlate with stellar age. At all redshifts, galaxies that are quenching or quenched, and thus old, have high $Σ_1$ (the projected density within the central 1 kpc), while younger, star-forming galaxies span a much broader range of $Σ_1$, which includes the high values observed for quenched galaxies, but also extends to much lower values. As galaxies age and quench, the stellar age and the dispersion of $Σ_1$ for fixed values of $M_{*}$ shows two different regimes, one, at the low--mass end, where quenching might be driven by causes external to the galaxies; the other, at the high--mass end, where quenching is driven by internal causes, very likely the mass given the low scatter of $Σ_1$ (mass quenching). We suggest that the monotonic increase of central density as galaxies grow is one manifestation of a more general phenomenon of structural transformation that galaxies undergo as they evolve. △ Less

Submitted 28 December, 2017; v1 submitted 7 June, 2017; originally announced June 2017.

Comments: Significant improvements in Section 3.2 and Section 7, but results unchanged. 46 pages, 32 figures, 2 tables, Accepted to ApJ

Do you see what I see? Exploring the consequences of luminosity limits in black hole-galaxy evolution studies

Authors: Mackenzie L. Jones, Ryan C. Hickox, Simon J. Mutch, Darren J. Croton, Andrew F. Ptak, Michael A. DiPompeo

Abstract: In studies of the connection between active galactic nuclei (AGN) and their host galaxies there is widespread disagreement on some key aspects stemming largely from a lack of understanding of the nature of the full underlying AGN population. Recent attempts to probe this connection utilize both observations and simulations to correct for a missed population, but presently are limited by intrinsic… ▽ More In studies of the connection between active galactic nuclei (AGN) and their host galaxies there is widespread disagreement on some key aspects stemming largely from a lack of understanding of the nature of the full underlying AGN population. Recent attempts to probe this connection utilize both observations and simulations to correct for a missed population, but presently are limited by intrinsic biases and complicated models. We take a simple simulation for galaxy evolution and add a new prescription for AGN activity to connect galaxy growth to dark matter halo properties and AGN activity to star formation. We explicitly model selection effects to produce an "observed" AGN population for comparison with observations and empirically motivated models of the local universe. This allows us to bypass the difficulties inherent in many models which attempt to infer the AGN population by inverting selection effects. We investigate the impact of selecting AGN based on thresholds in luminosity or Eddington ratio on the "observed" AGN population. By limiting our model AGN sample in luminosity, we are able to recreate the observed local AGN luminosity function and specific star formation-stellar mass distribution, and show that using an Eddington ratio threshold introduces less bias into the sample by selecting the full range of growing black holes, despite the challenge of selecting low mass black holes. We find that selecting AGN using these various thresholds yield samples with different AGN host galaxy properties. △ Less

Submitted 5 June, 2017; v1 submitted 1 June, 2017; originally announced June 2017.

Comments: 8 pages, 4 figures. Accepted for publication in ApJ

Galaxy And Mass Assembly (GAMA): 'No Smoking' zone for giant elliptical galaxies?

Authors: Habib G. Khosroshahi, Mojtaba Raouf, Halime Miraghaei, Sarah Brough, Darren J. Croton, Simon Driver, Alister Graham, Ivan Baldry, Michael Brown, Matt Prescott, Lingyu Wang

Abstract: We study the radio emission of the most massive galaxies in a sample of dynamically relaxed and un-relaxed galaxy groups from Galaxy and Mass Assembly (GAMA). The dynamical state of the group is defined by the stellar dominance of the brightest group galaxy, e.g. the luminosity gap between the two most luminous members, and the offset between the position of the brightest group galaxy and the lumi… ▽ More We study the radio emission of the most massive galaxies in a sample of dynamically relaxed and un-relaxed galaxy groups from Galaxy and Mass Assembly (GAMA). The dynamical state of the group is defined by the stellar dominance of the brightest group galaxy, e.g. the luminosity gap between the two most luminous members, and the offset between the position of the brightest group galaxy and the luminosity centroid of the group. We find that the radio luminosity of the most massive galaxy in the group strongly depends on its environment, such that the brightest group galaxies in dynamically young (evolving) groups are an order of magnitude more luminous in the radio than those with a similar stellar mass but residing in dynamically old (relaxed) groups. This observation has been successfully reproduced by a newly developed semi-analytic model which allows us to explore the various causes of these findings. We find that the fraction of radio loud brightest group galaxies in the observed dynamically young groups is ~2 times that in the dynamically old groups. We discuss the implications of this observational constraint on the central galaxy properties in the context of galaxy mergers and the super-massive blackhole accretion rate. △ Less

Submitted 15 April, 2018; v1 submitted 28 April, 2017; originally announced April 2017.

Comments: 10 pages, 8 figures. Accepted for publication in ApJ; added references

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

Pairwise velocities in the "Running FLRW" cosmological model

Authors: Antonio Bibiano, Darren J. Croton

Abstract: We present an analysis of the pairwise velocity statistics from a suite of cosmological N-body simulations describing the "Running Friedmann-Lemaître-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends ΛCDM with a time-evolving vacuum energy density. To enforce local conservation of matter a time-evolving gravitational coupl… ▽ More We present an analysis of the pairwise velocity statistics from a suite of cosmological N-body simulations describing the "Running Friedmann-Lemaître-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends ΛCDM with a time-evolving vacuum energy density. To enforce local conservation of matter a time-evolving gravitational coupling is also included. Our results constitute the first study of velocities in the R-FLRW cosmology, and we also compare with other dark energy simulations suites, repeating the same analysis. We find a strong degeneracy between the pairwise velocity and σ_8 at z=0 for almost all scenarios considered, which remains even when we look back to epochs as early as z=2. We also investigate various Coupled Dark Energy models, some of which show minimal degeneracy, and reveal interesting deviations from ΛCDM which could be readily exploited by future cosmological observations to test and further constrain our understanding of dark energy. △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: 10 pages, 4 figures. arXiv admin note: text overlap with arXiv:1605.00714

The evolution of star formation histories of quiescent galaxies

Authors: Camilla Pacifici, Susan A. Kassin, Benjamin J. Weiner, Bradford Holden, Jonathan P. Gardner, Sandra M. Faber, Henry C. Ferguson, David C. Koo, Joel R. Primack, Eric F. Bell, Avishai Dekel, Eric Gawiser, Mauro Giavalisco, Marc Rafelski, Raymond C. Simons, Guillermo Barro, Darren J. Croton, Romeel Dave, Adriano Fontana, Norman A. Grogin, Anton M. Koekemoer, Seong-Kook Lee, Brett Salmon, Rachel Somerville, Peter Behroozi

Abstract: Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2<z<2.1$. A total of… ▽ More Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2<z<2.1$. A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the CANDELS dataset. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected because the Universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim9.5$) grow on average slowly, take a long time to reach their peak of star formation ($\gtrsim 4$ Gyr), and the declining phase is fast ($\lesssim 2$ Gyr). Conversely, high-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim11$) grow on average fast ($\lesssim 2$ Gyr), and, after reaching their peak, decrease the star formation slowly ($\gtrsim 3$ Gyr). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms that operate on different timescales. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: 16 pages, 10 figures, resubmitted to ApJ after addressing the Referee's comments

How to get cool in the heat: comparing analytic models of hot, cold, and cooling gas in haloes and galaxies with EAGLE

Authors: Adam R. H. Stevens, Claudia del P. Lagos, Sergio Contreras, Darren J. Croton, Nelson D. Padilla, Matthieu Schaller, Joop Schaye, Tom Theuns

Abstract: We use the hydrodynamic, cosmological EAGLE simulations to investigate how hot gas in haloes condenses to form and grow galaxies. We select haloes from the simulations that are actively cooling and study the temperature, distribution, and metallicity of their hot, cold, and transitioning `cooling' gas, placing these in context of semi-analytic models. Our selection criteria lead us to focus on Mil… ▽ More We use the hydrodynamic, cosmological EAGLE simulations to investigate how hot gas in haloes condenses to form and grow galaxies. We select haloes from the simulations that are actively cooling and study the temperature, distribution, and metallicity of their hot, cold, and transitioning `cooling' gas, placing these in context of semi-analytic models. Our selection criteria lead us to focus on Milky Way-like haloes. We find the hot-gas density profiles of the haloes form a progressively stronger core over time, the nature of which can be captured by a beta profile that has a simple dependence on redshift. In contrast, the hot gas that will cool over a time-step is broadly consistent with a singular isothermal sphere. We find that cooling gas carries a few times the specific angular momentum of the halo and is offset in spin direction from the rest of the hot gas. The gas loses ~60% of its specific angular momentum during the cooling process, generally remaining greater than that of the halo, and it precesses to become aligned with the cold gas already in the disc. We find tentative evidence that angular-momentum losses are slightly larger when gas cools onto dispersion-supported galaxies. We show that an exponential surface density profile for gas arriving on a disc remains a reasonable approximation, but a cusp containing ~20% of the mass is always present, and disc scale radii are larger than predicted by a vanilla Fall & Efstathiou model. These scale radii are still closely correlated with the halo spin parameter, for which we suggest an updated prescription for galaxy formation models. △ Less

Submitted 24 January, 2017; v1 submitted 15 August, 2016; originally announced August 2016.

Comments: Accepted in MNRAS. 21 pages, 18 figures. Replacement includes fixed typos and reference updates

The Califa and Hipass velocity function for all morphological galaxy types

Authors: S. Bekeraitė, C. J. Walcher, L. Wisotzki, D. J. Croton, J. Falcón-Barroso, M. Lyubenova, D. Obreschkow, S. F. Sánchez, K. Spekkens, P. Torrey, G. van de Ven, M. A. Zwaan, Y. Ascasibar, J. Bland-Hawthorn, R. González-Delgado, B. Husemann, R. A. Marino, M. Vogelsberger, B. Ziegler, the CALIFA collaboration

Abstract: The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy p… ▽ More The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations. △ Less

Submitted 4 August, 2016; originally announced August 2016.

Comments: Accepted to ApJL. Keywords: Galaxies: kinematics and dynamics - galaxies: statistics - galaxies: evolution 7 pages, 3 figures

The growth of disks and bulges during hierarchical galaxy formation. II: metallicity, stellar populations and dynamical evolution

Authors: Chiara Tonini, Simon J. Mutch, J. Stuart B. Wyithe, Darren J. Croton

Abstract: We investigate the properties of the stellar populations of model galaxies as a function of galaxy evolutionary history and angular momentum content. We use the new semi-analytic model presented in Tonini et al. (2016). This new model follows the angular momentum evolution of gas and stars, providing the base for a new star formation recipe, and treatment of the effects of mergers that depends on… ▽ More We investigate the properties of the stellar populations of model galaxies as a function of galaxy evolutionary history and angular momentum content. We use the new semi-analytic model presented in Tonini et al. (2016). This new model follows the angular momentum evolution of gas and stars, providing the base for a new star formation recipe, and treatment of the effects of mergers that depends on the central galaxy dynamical structure. We find that the new recipes have the effect of boosting the efficiency of the baryonic cycle in producing and recycling metals, as well as preventing minor mergers from diluting the metallicity of bulges and ellipticals. The model reproduces the stellar mass - stellar metallicity relation for galaxies above 1e10 solar masses, including Brightest Cluster Galaxies. Model disks, galaxies dominated by instability-driven components, and merger-driven objects each stem from different evolutionary channels. These model galaxies therefore occupy different loci in the galaxy mass-size relation, which we find to be in accord with the Atlas 3D classification of disk galaxies, fast rotators and slow rotators. We find that the stellar populations' properties depend on the galaxy evolutionary type, with more evolved stellar populations being part of systems that have lost or dissipated more angular momentum during their assembly history. △ Less

Submitted 7 December, 2016; v1 submitted 20 June, 2016; originally announced June 2016.

Comments: 15 pages, 9 figures (3 new figures), accepted on MNRAS

Non-linear structure formation in the "Running FLRW" cosmological model

Authors: Antonio Bibiano, Darren J. Croton

Abstract: We present a suite of cosmological N-body simulations describing the "Running Friedmann-Lema{ï}tre-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends ΛCDM with a time-evolving vacuum density, Λ(z), and time-evolving gravitational Newton's coupling, G(z). In this paper we review the model and introduce the necessary analytic… ▽ More We present a suite of cosmological N-body simulations describing the "Running Friedmann-Lema{ï}tre-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends ΛCDM with a time-evolving vacuum density, Λ(z), and time-evolving gravitational Newton's coupling, G(z). In this paper we review the model and introduce the necessary analytical treatment needed to adapt a reference N-body code. Our resulting simulations represent the first realisation of the full growth history of structure in the R-FLRW cosmology into the non-linear regime, and our normalisation choice makes them fully consistent with the latest cosmic microwave background data. The post-processing data products also allow, for the first time, an analysis of the properties of the halo and sub-halo populations. We explore the degeneracies of many statistical observables and discuss the steps needed to break them. Furthermore, we provide a quantitative description of the deviations of R-FLRW from ΛCDM, which could be readily exploited by future cosmological observations to test and further constrain the model. △ Less

Submitted 2 May, 2016; originally announced May 2016.

Comments: 14 pages, 7 figures

Building disc structure and galaxy properties through angular momentum: The DARK SAGE semi-analytic model

Authors: Adam R. H. Stevens, Darren J. Croton, Simon J. Mutch

Abstract: We present the new semi-analytic model of galaxy evolution, DARK SAGE, a heavily modified version of the publicly available SAGE code. The model is designed for detailed evolution of galactic discs. We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, inclu… ▽ More We present the new semi-analytic model of galaxy evolution, DARK SAGE, a heavily modified version of the publicly available SAGE code. The model is designed for detailed evolution of galactic discs. We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, including all channels of star formation and associated feedback, are performed in these annuli. We present the surface density profiles of our model spiral galaxies, both as a function of radius and specific angular momentum, and find the discs naturally build a pseduobulge-like component. Our main results are focussed on predictions relating to the integrated mass--specific angular momentum relation of stellar discs. The model produces a distinct sequence between these properties in remarkable agreement with recent observational literature. We investigate the impact Toomre disc instabilities have on shaping this sequence and find they are crucial for regulating both the mass and spin of discs. Without instabilities, high-mass discs would be systematically deficient in specific angular momentum by a factor of ~2.5, with increased scatter. Instabilities also appear to drive the direction in which the mass--spin sequence of spiral galaxy discs evolves. With them, we find galaxies of fixed mass have higher specific angular momentum at later epochs. △ Less

Submitted 24 June, 2016; v1 submitted 2 May, 2016; originally announced May 2016.

Comments: Published in MNRAS, 19 pages, 20 figures (both including appendices)

The growth of disks and bulges during hierarchical galaxy formation. I: fast evolution vs secular processes

Authors: Chiara Tonini, Simon J. Mutch, Darren J. Croton, J. Stuart B. Wyithe

Abstract: We present a theoretical model for the evolution of mass, angular momentum and size of galaxy disks and bulges, and we implement it into the semi-analytic galaxy formation code SAGE. The model follows both secular and violent evolutionary channels, including smooth accretion, disk instabilities, minor and major mergers. We find that the combination of our recipe with hierarchical clustering produc… ▽ More We present a theoretical model for the evolution of mass, angular momentum and size of galaxy disks and bulges, and we implement it into the semi-analytic galaxy formation code SAGE. The model follows both secular and violent evolutionary channels, including smooth accretion, disk instabilities, minor and major mergers. We find that the combination of our recipe with hierarchical clustering produces two distinct populations of bulges: merger-driven bulges, akin to classical bulges and ellipticals, and instability-driven bulges, akin to secular (or pseudo-)bulges. The model mostly reproduces the mass-size relation of gaseous and stellar disks, the evolution of the mass-size relation of ellipticals, the Faber-Jackson relation, and the magnitude-colour diagram of classical and secular bulges. The model predicts only a small overlap of merger-driven and instability-driven components in the same galaxy, and predicts different bulge types as a function of galaxy mass and disk fraction. Bulge type also affects the star formation rate and colour at a given luminosity. The model predicts a population of merger-driven red ellipticals that dominate both the low-mass and high-mass ends of the galaxy population, and span all dynamical ages; merger-driven bulges in disk galaxies are dynamically old and do not interfere with subsequent evolution of the star-forming component. Instability-driven bulges dominate the population at intermediate galaxy masses, especially thriving in massive disks. The model green valley is exclusively populated by instability-driven bulge hosts. Through the present implementation the mass accretion history is perceivable in the galaxy structure, morphology and colours. △ Less

Submitted 20 April, 2016; v1 submitted 7 April, 2016; originally announced April 2016.

Comments: Accepted on MNRAS; 24 pages, 11 figures

Semi-Analytic Galaxy Evolution (SAGE): Model Calibration and Basic Results

Authors: Darren J. Croton, Adam R. H. Stevens, Chiara Tonini, Thibault Garel, Maksym Bernyk, Antonio Bibiano, Luke Hodkinson, Simon J. Mutch, Gregory B. Poole, Genevieve M. Shattow

Abstract: This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised in a supported format and contain a… ▽ More This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised in a supported format and contain a minimum set of basic halo properties. In this work we present the baryonic prescriptions implemented in SAGE to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling--radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population. △ Less

Submitted 21 January, 2016; v1 submitted 18 January, 2016; originally announced January 2016.

Comments: 15 pages, 9 figures, accepted for publication in ApJS. SAGE is a publicly available codebase for modelling galaxy formation in a cosmological context, available at https://github.com/darrencroton/sage Questions and comments can be sent to Darren Croton: dcroton@astro.swin.edu.au

Galaxy And Mass Assembly (GAMA): The Bright Void Galaxy Population in the Optical and Mid-IR

Authors: S. J. Penny, M. J. I. Brown, K. A. Pimbblet, M. E. Cluver, D. J. Croton, M. S. Owers, R. Lange, M. Alpaslan, I. Baldry, J. Bland-Hawthorn, S. Brough, S. P. Driver, B. W. Holwerda, A. M. Hopkins, T. H. Jarrett, D. Heath Jones, L. S. Kelvin, M. A. Lara-Lopez, J. Liske, A. R. Lopez-Sanchez, J. Loveday, M. Meyer, P. Norberg, A. S. G. Robotham, M. Rodrigues

Abstract: We examine the properties of galaxies in the Galaxies and Mass Assembly (GAMA) survey located in voids with radii $>10~h^{-1}$ Mpc. Utilising the GAMA equatorial survey, 592 void galaxies are identified out to z~0.1 brighter than $M_{r} = -18.4$, our magnitude completeness limit. Using the $W_{\rm{Hα}}$ vs. [NII]/H$α$ (WHAN) line strength diagnostic diagram, we classify their spectra as star formi… ▽ More We examine the properties of galaxies in the Galaxies and Mass Assembly (GAMA) survey located in voids with radii $>10~h^{-1}$ Mpc. Utilising the GAMA equatorial survey, 592 void galaxies are identified out to z~0.1 brighter than $M_{r} = -18.4$, our magnitude completeness limit. Using the $W_{\rm{Hα}}$ vs. [NII]/H$α$ (WHAN) line strength diagnostic diagram, we classify their spectra as star forming, AGN, or dominated by old stellar populations. For objects more massive than $5\times10^{9}$ M$_{\odot}$, we identify a sample of 26 void galaxies with old stellar populations classed as passive and retired galaxies in the WHAN diagnostic diagram, else they lack any emission lines in their spectra. When matched to WISE mid-IR photometry, these passive and retired galaxies exhibit a range of mid-IR colour, with a number of void galaxies exhibiting [4.6]-[12] colours inconsistent with completely quenched stellar populations, with a similar spread in colour seen for a randomly drawn non-void comparison sample. We hypothesise that a number of these galaxies host obscured star formation, else they are star forming outside of their central regions targeted for single fibre spectroscopy. When matched to a randomly drawn sample of non-void galaxies, the void and non-void galaxies exhibit similar properties in terms of optical and mid-IR colour, morphology, and star formation activity, suggesting comparable mass assembly and quenching histories. A trend in mid-IR [4.6]-[12] colour is seen, such that both void and non-void galaxies with quenched/passive colours <1.5 typically have masses higher than $10^{10}$ M$_{\odot}$, where internally driven processes play an increasingly important role in galaxy evolution. △ Less

Submitted 25 August, 2015; originally announced August 2015.

Comments: 22 pages, 12 figures. Accepted for publication in MNRAS

Beyond the Halo: Redefining environment with unbound matter in N-body simulations

Authors: Genevieve M. Shattow, Darren J. Croton

Abstract: Approximately half of the matter in the Universe is "unbound" at z = 0, according to N-body simulations such as the Millennium Run. Here, we use the milli-Millennium simulation to examine the distribution of unbound matter in relation to the dark matter halos which host galaxies. We measure the unbound matter within two types of windows, using a halo dependent radius and a fixed radius at several… ▽ More Approximately half of the matter in the Universe is "unbound" at z = 0, according to N-body simulations such as the Millennium Run. Here, we use the milli-Millennium simulation to examine the distribution of unbound matter in relation to the dark matter halos which host galaxies. We measure the unbound matter within two types of windows, using a halo dependent radius and a fixed radius at several different scales. We also consider the timescales over which a halo can accrete the local unbound matter at z = 2 and z = 0. Finally, we compare the unbound matter to observable properties of galaxies, such as local galaxy count environment and stellar mass. We find that halos at z = 2 can accrete far more of the nearby unbound matter over a Hubble time than halos at z = 0 and that 78% of particles within 5 $R_{vir}$ of a halo at z = 2 will be accreted by z = 0, compared to 36% of particles within 5 $h^{-1}$ Mpc of the halo. We also find that galaxy count environment is closely related to the amount of nearby unbound matter when measured on the same scale. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: 14 pages, 9 Figures, Accepted MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 452 (2): 1779-1791

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

Thinking outside the halo: Tracing the large-scale distribution of diffuse cosmic metals with semi-analytic models

Authors: Genevieve M. Shattow, Darren J. Croton, Antonio Bibiano

Abstract: With the installation of the Cosmic Origins Spectrograph on the Hubble Space Telescope, measurements of the metal content of the low redshift intergalactic medium (IGM) are now available. Using a new grid-based model for diffuse gas coupled to the SAGE semi-analytic model of galaxy formation, we examine the impact of supernova feedback on the pollution of the IGM. We consider different assumptions… ▽ More With the installation of the Cosmic Origins Spectrograph on the Hubble Space Telescope, measurements of the metal content of the low redshift intergalactic medium (IGM) are now available. Using a new grid-based model for diffuse gas coupled to the SAGE semi-analytic model of galaxy formation, we examine the impact of supernova feedback on the pollution of the IGM. We consider different assumptions for the reheating and ejection of gas by supernovae and their dependence on galaxy circular velocity and gas surface density. Where metals are present, we find the most likely metallicity to be $-1.5 < $log$_{10}$(Z/Z$_{\odot}$)$< -1.0$ at $z = 0$, consistent with both observations and more sophisticated hydrodynamic simulations. Our model predicts that the regions of the IGM with the highest metallicities will be near galaxies with M$_{\star} \sim 10^{10.5}h^{-1}$M$_{\odot}$ and in environments of densities $\sim 10 \times$ the mean. We also find that 90% of IGM metals at $z = 0$ are ejected by galaxies with stellar masses less than $10^{10.33}h^{-1}$M$_{\odot}$. △ Less

Submitted 26 March, 2015; originally announced March 2015.

Comments: 12 pages, 9 figures, accepted MNRAS

Modeling the Ages and Metallicities of Early-Type Galaxies in Fundamental Plane Space

Authors: L. A. Porter, R. S. Somerville, J. R. Primack, D. J. Croton, M. D. Covington, G. J. Graves, S. M. Faber

Abstract: Recent observations have probed the formation histories of nearby elliptical galaxies by tracking correlations between the stellar population parameters, age and metallicity, and the structural parameters that enter the Fundamental Plane, radius and velocity dispersion sigma. These studies have found intriguing correlations between these four parameters. In this work, we make use of a semi-analyti… ▽ More Recent observations have probed the formation histories of nearby elliptical galaxies by tracking correlations between the stellar population parameters, age and metallicity, and the structural parameters that enter the Fundamental Plane, radius and velocity dispersion sigma. These studies have found intriguing correlations between these four parameters. In this work, we make use of a semi-analytic model, based on halo merger trees extracted from the Bolshoi cosmological simulation, that predicts the structural properties of spheroid-dominated galaxies based on an analytic model that has been tested and calibrated against an extensive suite of hydrodynamic+N-body binary merger simulations. We predict the radius, sigma, luminosity, age, and metallicity of spheroid-dominated galaxies, enabling us to compare directly to observations. Our model predicts a strong correlation between age and sigma for early-type galaxies, and no significant correlation between age and radius, in agreement with observations. In addition we predict a strong correlation between metallicity and sigma, and a weak correlation between metallicity and radius, in qualitative agreement with observations. We find that the correlations with sigma arise as a result of the strong link between sigma and the galaxy assembly time. Minor mergers produce a large change in radius while leaving sigma nearly the same, which explains the weaker trends with radius. △ Less

Submitted 8 July, 2014; originally announced July 2014.

Comments: 15 pages, 9 figures, submitted to MNRAS. arXiv admin note: text overlap with arXiv:1201.5918

The Gigaparsec WiggleZ Simulations: Characterising scale dependant bias and associated systematics in growth of structure measurements

Authors: Gregory B. Poole, Chris Blake, Felipe A. Marin, Chris Power, Simon J. Mutch, Darren J. Croton, Matthew Colless, Warrick Couch, Michael J. Drinkwater, Karl Glazebrook

Abstract: We use the Gigaparsec WiggleZ (GiggleZ) simulations to characterise galaxy bias and its scale dependence for a range of redshifts and halo masses in a standard $Λ$LCDM cosmology. Assuming bias converges to a scale independent form at large scales, we develop a phenomenological model which fully expresses the mass and redshift dependence of bias and its scale dependence in real or redshift space. W… ▽ More We use the Gigaparsec WiggleZ (GiggleZ) simulations to characterise galaxy bias and its scale dependence for a range of redshifts and halo masses in a standard $Λ$LCDM cosmology. Assuming bias converges to a scale independent form at large scales, we develop a phenomenological model which fully expresses the mass and redshift dependence of bias and its scale dependence in real or redshift space. We then use this to illustrate how scale-dependent bias can systematically skew measurements of the growth-rate of cosmic structure obtained from redshift-space distortion measurements. When data is fit only to scales $k_{\rm max}{\le}0.1$ $[h^{-1} \rm{Mpc}]^{-1}$, we find that these effects are significant only for large biases ($b{\gtrsim}3$) at large redshifts ($z{\gtrsim}1$). However, when smaller scales are incorporated ($k_{\rm max}{\le}0.2$ $[h^{-1} \rm{Mpc}]^{-1}$) to increase measurement precision, the combination of reduced statistical uncertainty and increased scale dependent bias can result in highly significant systematics for most large halos across all redshifts. We identify several new interesting aspects of bias, including a significant large-scale bias boost for small halos at low-redshifts due to substructure effects ($\sim$20\% for Milky Way-like systems) and a nearly redshift-independent halo mass (corresponding to a redshift-space bias of ${\sim}1.5$) for which halo bias has little-or-no scale dependence on scales greater than $3$ $[h^{-1} {\rm Mpc}]$. This suggests an optimal strategy of targeting bias ${\sim}{1.5}$ systems for clustering studies which are dominated by systematic uncertainties in how observed halo (or galaxy) distributions map to their underlying mass distribution, such as cosmological measurements of neutrino masses. Code for generating our fitting formula is publicly available at http://gbpoole.github.io/Poole_2014a_code/ (Abridged) △ Less

Submitted 2 November, 2015; v1 submitted 1 July, 2014; originally announced July 2014.

Comments: 17 Pages, 10 Figures, Accepted for publication in MNRAS

Journal ref: 2015MNRAS.449.1454P

Where do galaxies end? Comparing measurement techniques of hydrodynamic-simulation galaxies' integrated properties

Authors: Adam R. H. Stevens, Marie Martig, Darren J. Croton, Yu Feng

Abstract: Using the suite of high-resolution zoom re-simulations of individual haloes by Martig et al., and the large-scale simulation \emph{MassiveBlack-II}, we examine the differences in measured galaxy properties from techniques with various aperture definitions of where galaxies end. We perform techniques popular in the literature and present a new technique of our own, where the aperture radius is base… ▽ More Using the suite of high-resolution zoom re-simulations of individual haloes by Martig et al., and the large-scale simulation \emph{MassiveBlack-II}, we examine the differences in measured galaxy properties from techniques with various aperture definitions of where galaxies end. We perform techniques popular in the literature and present a new technique of our own, where the aperture radius is based on the baryonic mass profiles of simulated (sub)haloes. For the average Milky-Way-mass system, we find the two most popular techniques in the literature return differences of order 30 per cent for stellar mass, a factor of 3 for gas mass, 40 per cent for star formation rate, and factors of several for gas accretion and ejection rates. Individual cases can show variations greater than this, with the severity dependent on the concentration of a given system. The average difference in integrated properties for a more general galaxy population are not as striking, but are still significant for stellar and gas mass, especially for optical-limit apertures. The large differences that can occur are problematic for comparing results from various publications. We stress the importance of both defining and justifying a technique choice and discourage using popular apertures that use an exact fraction of the virial radius, due to the unignorable variation in galaxy-to-(sub)halo size. Finally, we note that technique choice does not greatly affect simulated galaxies from lying within the scatter of observed scaling relations, but it can alter the derived best-fit slope for the Kennicutt-Schmidt relation. △ Less

Submitted 21 August, 2014; v1 submitted 15 April, 2014; originally announced April 2014.

Comments: Accepted by MNRAS, 18 pages, 16 figures

Journal ref: Monthly Notices of the Royal Astronomical Society 2014 445 (2): 239-255

The Theoretical Astrophysical Observatory: Cloud-Based Mock Galaxy Catalogues

Authors: Maksym Bernyk, Darren J. Croton, Chiara Tonini, Luke Hodkinson, Amr H. Hassan, Thibault Garel, Alan R. Duffy, Simon J. Mutch, Gregory B. Poole, Sarah Hegarty

Abstract: We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance compu… ▽ More We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalogue suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light-cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future. △ Less

Submitted 21 January, 2016; v1 submitted 20 March, 2014; originally announced March 2014.

Comments: 17 pages, 11 figures, 2 tables; accepted for publication in ApJS. The Theoretical Astrophysical Observatory (TAO) is now open to the public at https://tao.asvo.org.au/. New simulations, models and tools will be added as they become available. Contact dcroton@astro.swin.edu.au if you have data you would like to make public through TAO. Feedback and suggestions are very welcome

The WiggleZ Dark Energy Survey: Improved Distance Measurements to z = 1 with Reconstruction of the Baryonic Acoustic Feature

Authors: Eyal A. Kazin, Jun Koda, Chris Blake, Nikhil Padmanabhan, Sarah Brough, Matthew Colless, Carlos Contreras, Warrick Couch, Scott Croom, Darren J. Croton, Tamara M. Davis, Michael J. Drinkwater, Karl Forster, David Gilbank, Mike Gladders, Karl Glazebrook, Ben Jelliffe, Russell J. Jurek, I-hui Li, Barry Madore, D. Christopher Martin, Kevin Pimbblet, Gregory B. Poole, Michael Pracy, Rob Sharp , et al. (4 additional authors not shown)

Abstract: We present significant improvements in cosmic distance measurements from the WiggleZ Dark Energy Survey, achieved by applying the reconstruction of the baryonic acoustic feature technique. We show using both data and simulations that the reconstruction technique can often be effective despite patchiness of the survey, significant edge effects and shot-noise. We investigate three redshift bins in t… ▽ More We present significant improvements in cosmic distance measurements from the WiggleZ Dark Energy Survey, achieved by applying the reconstruction of the baryonic acoustic feature technique. We show using both data and simulations that the reconstruction technique can often be effective despite patchiness of the survey, significant edge effects and shot-noise. We investigate three redshift bins in the redshift range 0.2<$z$<1, and in all three find improvement after reconstruction in the detection of the baryonic acoustic feature and its usage as a standard ruler. We measure model independent distance measures $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ of 1716 $\pm$ 83 Mpc, 2221 $\pm$ 101 Mpc, 2516 $\pm$ 86 Mpc (68% CL) at effective redshifts z = 0.44, 0.6, 0.73, respectively, where $D_{\mathrm V}$ is the volume-average-distance, and $r_{\mathrm s}$ is the sound horizon at the end of the baryon drag epoch. These significantly improved 4.8, 4.5 and 3.4 percent accuracy measurements are equivalent to those expected from surveys with up to 2.5 times the volume of WiggleZ. These measurements are fully consistent with cosmologies allowed by the analyses of the Planck Collaboration and the Sloan Digital Sky Survey.We provide the $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ posterior probability distributions and their covariances. When combining these measurements with temperature fluctuations measurements of Planck, the polarization of WMAP9, and the 6dF Galaxy Survey baryonic acoustic feature, we do not detect deviations from a flat LCDM model. Assuming this model we constrain the current expansion rate to $H_0$ = 67.15 $\pm$ 0.98 kms$^{-1}$Mpc$^{-1}$. Allowing the equation of state of dark energy to vary we obtain $w_\mathrm{DE}$ = -1.080 $\pm$ 0.135. When assuming a curved LCDM model we obtain a curvature value of $Ω_{\mathrm K}$ = -0.0043 $\pm$ 0.0047. △ Less

Submitted 28 July, 2014; v1 submitted 1 January, 2014; originally announced January 2014.

Comments: 20 double column pages, 10 figures, 5 tables

Journal ref: Monthly Notices of the Royal Astronomical Society 441, 3524 (2014)

Measures of Galaxy Environment - III. Difficulties in identifying proto-clusters at z ~ 2

Authors: Genevieve M. Shattow, Darren J. Croton, Ramin A. Skibba, Stuart I. Muldrew, Frazer R. Pearce, Ummi Abbas

Abstract: Galaxy environment is frequently discussed, but inconsistently defined. It is especially difficult to measure at high redshift where only photometric redshifts are available. With a focus on early forming proto-clusters, we use a semi-analytical model of galaxy formation to show how the environment measurement around high redshift galaxies is sensitive to both scale and metric, as well as to clust… ▽ More Galaxy environment is frequently discussed, but inconsistently defined. It is especially difficult to measure at high redshift where only photometric redshifts are available. With a focus on early forming proto-clusters, we use a semi-analytical model of galaxy formation to show how the environment measurement around high redshift galaxies is sensitive to both scale and metric, as well as to cluster viewing angle, evolutionary state, and the availability of either spectroscopic or photometric data. We use two types of environment metrics (nearest neighbour and fixed aperture) at a range of scales on simulated high-z clusters to see how "observed" overdensities compare to "real" overdensities. We also "observationally" identify z = 2 proto-cluster candidates in our model and track the growth histories of their parent halos through time, considering in particular their final state at z = 0. Although the measured environment of early forming clusters is critically dependent on all of the above effects (and in particular the viewing angle), we show that such clusters are very likely (< 90%) to remain overdense at z = 0, although many will no longer be among the most massive. Object to object comparisons using different methodologies and different data, however, require much more caution. △ Less

Submitted 7 June, 2013; originally announced June 2013.

Comments: 12 pages, 9 figures, accepted to MNRAS

The simplest model of galaxy formation I: A formation history model of galaxy stellar mass growth

Authors: Simon J. Mutch, Darren J. Croton, Gregory B. Poole

Abstract: We introduce a simple model to self-consistently connect the growth of galaxies to the formation history of their host dark matter haloes. Our model is defined by two simple functions: the "baryonic growth function" which controls the rate at which new baryonic material is made available for star formation, and the "physics function" which controls the efficiency with which this material is conver… ▽ More We introduce a simple model to self-consistently connect the growth of galaxies to the formation history of their host dark matter haloes. Our model is defined by two simple functions: the "baryonic growth function" which controls the rate at which new baryonic material is made available for star formation, and the "physics function" which controls the efficiency with which this material is converted into stars. Using simple, phenomenologically motivated forms for both functions that depend only on a single halo property, we demonstrate the model's ability to reproduce the z=0 red and blue stellar mass functions. Furthermore, by adding redshift as a second input variable to the physics function we show that the reproduction of the global stellar mass function out to z=3 is improved. We conclude by discussing the general utility of our new model, highlighting its usefulness for creating mock galaxy samples which have a number of key advantages over those generated by other techniques. △ Less

Submitted 3 September, 2013; v1 submitted 9 April, 2013; originally announced April 2013.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS. New accepted version contains an extended discussion, minor corrections to Table 1 and other small changes/additions. Main results and conclusions are unchanged

The Stellar Masses of ~40,000 UV Selected Galaxies from the WiggleZ Survey at 0.3<z<1.0: Analogues of Lyman Break Galaxies?

Authors: Manda Banerji, Karl Glazebrook, Chris Blake, Sarah Brough, Matthew Colless, Carlos Contreras, Warrick Couch, Darren J. Croton, Scott Croom, Tamara M. Davis, Michael J. Drinkwater, Karl Forster, David Gilbank, Mike Gladders, Ben Jelliffe, Russell J. Jurek, I-hui Li, Barry Madore, D. Christopher Martin, Kevin Pimbblet, Gregory B. Poole, Michael Pracy, Rob Sharp, Emily Wisnioski, David Woods , et al. (2 additional authors not shown)

Abstract: We characterise the stellar masses and star formation rates in a sample of almost 40000 spectroscopically confirmed UV luminous galaxies at 0.3<z<1.0 selected from within the WiggleZ Dark Energy Survey. In particular, we match this UV bright population to wide-field infrared surveys such as the near infrared UKIDSS LAS and the mid infrared WISE All-Sky Survey. We find that ~30% of the UV luminous… ▽ More We characterise the stellar masses and star formation rates in a sample of almost 40000 spectroscopically confirmed UV luminous galaxies at 0.3<z<1.0 selected from within the WiggleZ Dark Energy Survey. In particular, we match this UV bright population to wide-field infrared surveys such as the near infrared UKIDSS LAS and the mid infrared WISE All-Sky Survey. We find that ~30% of the UV luminous WiggleZ galaxies are detected at >5sigma in the UKIDSS-LAS at all redshifts. An even more luminous subset of 15% are also detected in the WISE 3.4 and 4.6um bands. We compute stellar masses for this very large sample of extremely blue galaxies and quantify the sensitivity of the stellar mass estimates to various assumptions made during the SED fitting. The median stellar masses are log10(M*/M0)=9.6\pm0.7, 10.2\pm0.5 and 10.4\pm0.4 for the IR-undetected, UKIDSS detected and UKIDSS+WISE detected galaxies respectively. We demonstrate that the inclusion of NIR photometry can lead to tighter constraints on the stellar masses. The mass estimates are found to be most sensitive to the inclusion of secondary bursts of star formation as well as changes in the stellar population synthesis models, both of which can lead to median discrepancies of the order of 0.3dex in the stellar masses. We find that the best-fit M/LK is significantly lower (by ~0.4 dex) than that predicted by simple optical colour based estimators, in particular for the bluer galaxies with younger best-fit ages. The WiggleZ galaxies have star formation rates of 3-10 M0/yr and mostly lie at the upper end of the main sequence of star-forming galaxies at these redshifts. Their rest-frame UV luminosities and stellar masses are comparable to both local compact UV-luminous galaxies as well as Lyman break galaxies at z~2-3.(abridged) △ Less

Submitted 20 February, 2013; originally announced February 2013.

Comments: 25 pages including two appendices, Accepted for publication in MNRAS

The Morphologies of Massive Galaxies at 1 < z < 3 in the CANDELS-UDS Field

Authors: V. A. Bruce, J. S. Dunlop, M. Cirasuolo, R. J. McLure, T. A. Targett, E. F. Bell, D. J. Croton, A. Dekel, S. M. Faber, H. C. Ferguson, N. A. Grogin, D. D. Kocevski, A. M. Koekemoer, D. C. Koo, K. Lai, J. M. Lotz, E. J. McGrath, J. A. Newman, A. van der Wel

Abstract: We have used high-resolution, HST WFC3/IR, near-infrared imaging to conduct a detailed bulge-disk decomposition of the morphologies of ~200 of the most massive (M_star > 10^11 M_solar) galaxies at 1<z<3 in the CANDELS-UDS field. We find that, while such massive galaxies at low redshift are generally bulge-dominated, at redshifts 1<z<2 they are predominantly mixed bulge+disk systems, and by z>2 the… ▽ More We have used high-resolution, HST WFC3/IR, near-infrared imaging to conduct a detailed bulge-disk decomposition of the morphologies of ~200 of the most massive (M_star > 10^11 M_solar) galaxies at 1<z<3 in the CANDELS-UDS field. We find that, while such massive galaxies at low redshift are generally bulge-dominated, at redshifts 1<z<2 they are predominantly mixed bulge+disk systems, and by z>2 they are mostly disk-dominated. Interestingly, we find that while most of the quiescent galaxies are bulge-dominated, a significant fraction (25-40%) of the most quiescent galaxies, have disk-dominated morphologies. Thus, our results suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies. △ Less

Submitted 27 January, 2013; originally announced January 2013.

Comments: Contributed article (4 pages, 2 figures), IAU Symposium 295 on "The intriguing life of massive galaxies", to appear in the proceedings of the IAU XXVIII General Assembly, Beijing, 2012, (Eds: Daniel Thomas, Anna Pasquali, Ignacio Ferreras)

The WiggleZ Dark Energy Survey: Probing the Epoch of Radiation Domination using Large Scale Structure

Authors: Gregory B. Poole, Chris Blake, David Parkinson, Sarah Brough, Matthew Colless, Carlos Contreras, Warrick Couch, Darren J. Croton, Scott Croom, Tamara Davis, Michael J. Drinkwater, Karl Forster, David Gilbank, Mike Gladders, Karl Glazebrook, Ben Jelliffe, Russell J. Jurek, I-hui Li, Barry Madore, D. Christopher Martin, Kevin Pimbblet, Michael Pracy, Rob Sharp, Emily Wisnioski, David Woods , et al. (2 additional authors not shown)

Abstract: We place the most robust constraint to date on the scale of the turnover in the cosmological matter power spectrum using data from the WiggleZ Dark Energy Survey. We find this feature to lie at a scale of $k_0=0.0160^{+0.0041}_{-0.0035}$ [h/Mpc] (68% confidence) for an effective redshift of 0.62 and obtain from this the first-ever turnover-derived distance and cosmology constraints: a measure of t… ▽ More We place the most robust constraint to date on the scale of the turnover in the cosmological matter power spectrum using data from the WiggleZ Dark Energy Survey. We find this feature to lie at a scale of $k_0=0.0160^{+0.0041}_{-0.0035}$ [h/Mpc] (68% confidence) for an effective redshift of 0.62 and obtain from this the first-ever turnover-derived distance and cosmology constraints: a measure of the cosmic distance-redshift relation in units of the horizon scale at the redshift of radiation-matter equality (r_H) of D_V(z=0.62)/r_H=18.3 (+6.3/-3.3) and, assuming a prior on the number of extra relativistic degrees of freedom $N_{eff}=3$, constraints on the matter density parameter $Ω_Mh^2=0.136^{+0.026}_{-0.052}$ and on the redshift of matter-radiation equality $z_{eq}=3274^{+631}_{-1260}$. All results are in excellent agreement with the predictions of standard LCDM models. Our constraints on the logarithmic slope of the power spectrum on scales larger than the turnover is bounded in the lower limit with values only as low as -1 allowed, with the prediction of standard LCDM models easily accommodated by our results. Lastly, we generate forecasts for the achievable precision of future surveys at constraining $k_0$, $Ω_Mh^2$, $z_{eq}$ and $N_{eff}$. We find that BOSS should substantially improve upon the WiggleZ turnover constraint, reaching a precision on $k_0$ of $\pm$9% (68% confidence), translating to precisions on $Ω_Mh^2$ and $z_{eq}$ of $\pm$10% (assuming a prior $N_{eff}=3$) and on $N_{eff}$ of (+78/-56)% (assuming a prior $Ω_Mh^2=0.135$). This is sufficient precision to sharpen the constraints on $N_{eff}$ from WMAP, particularly in its upper limit. For Euclid, we find corresponding attainable precisions on $(k_0, Ω_Mh^2, N_eff)$ of (3,4,+17/-21)%. This represents a precision approaching our forecasts for the Planck Surveyor. △ Less

Submitted 2 November, 2015; v1 submitted 23 November, 2012; originally announced November 2012.

Comments: 12 pages, 7 figures, accepted for publication in MNRAS

Measures of galaxy environment -- II. Rank-ordered mark correlations

Authors: Ramin A. Skibba, Ravi K. Sheth, Darren J. Croton, Stuart I. Muldrew, Ummi Abbas, Frazer R. Pearce, Genevieve M. Shattow

Abstract: We analyze environmental correlations using mark clustering statistics with the mock galaxy catalogue constructed by Muldrew et al. (Paper I). We find that mark correlation functions are able to detect even a small dependence of galaxy properties on the environment, quantified by the overdensity $1+δ$, while such a small dependence would be difficult to detect by traditional methods. We then show… ▽ More We analyze environmental correlations using mark clustering statistics with the mock galaxy catalogue constructed by Muldrew et al. (Paper I). We find that mark correlation functions are able to detect even a small dependence of galaxy properties on the environment, quantified by the overdensity $1+δ$, while such a small dependence would be difficult to detect by traditional methods. We then show that rank ordering the marks and using the rank as a weight is a simple way of comparing the correlation signals for different marks. With this we quantify to what extent fixed-aperture overdensities are sensitive to large-scale halo environments, nearest-neighbor overdensities are sensitive to small-scale environments within haloes, and colour is a better tracer of overdensity than is luminosity. △ Less

Submitted 29 January, 2013; v1 submitted 1 November, 2012; originally announced November 2012.

Comments: 12 pages, 15 figures; published in MNRAS. The mock galaxy catalogues and environment measures used in this paper and Paper I (arXiv:1109.6328) are available here: http://tao.it.swin.edu.au/mock-galaxy-factory/precomputed-mock-catalogues/

Journal ref: 2013, MNRAS, 429, 458

Constraining the last 7 billion years of galaxy evolution in semi-analytic models

Authors: Simon J. Mutch, Gregory B. Poole, Darren J. Croton

Abstract: We investigate the ability of the Croton et al. (2006) semi-analytic model to reproduce the evolution of observed galaxies across the final 7 billion years of cosmic history. Using Monte-Carlo Markov Chain techniques we explore the available parameter space to produce a model which attempts to achieve a statistically accurate fit to the observed stellar mass function at z=0 and z~0.8, as well as t… ▽ More We investigate the ability of the Croton et al. (2006) semi-analytic model to reproduce the evolution of observed galaxies across the final 7 billion years of cosmic history. Using Monte-Carlo Markov Chain techniques we explore the available parameter space to produce a model which attempts to achieve a statistically accurate fit to the observed stellar mass function at z=0 and z~0.8, as well as the local black hole-bulge relation. We find that in order to be successful we are required to push supernova feedback efficiencies to extreme limits which are, in some cases, unjustified by current observations. This leads us to the conclusion that the current model may be incomplete. Using the posterior probability distributions provided by our fitting, as well as the qualitative details of our produced stellar mass functions, we suggest that any future model improvements must act to preferentially bolster star formation efficiency in the most massive halos at high redshift. △ Less

Submitted 29 November, 2012; v1 submitted 11 October, 2012; originally announced October 2012.

Comments: 17 pages, 10 figures, 1 table. Correction to eqn. 10 and minor changes to text. Accepted for publication in MNRAS

Predictions for ASKAP Neutral Hydrogen Surveys

Authors: Alan R. Duffy, Martin J. Meyer, Lister Staveley-Smith, Maksym Bernyk, Darren J. Croton, Barbel S. Koribalski, Derek Gerstmann, Stefan Westerlund

Abstract: The Australian Square Kilometer Array Pathfinder (ASKAP) will revolutionise our knowledge of gas-rich galaxies in the Universe. Here we present predictions for two proposed extragalactic ASKAP neutral hydrogen (HI) emission-line surveys, based on semi-analytic models applied to cosmological N-body simulations. The ASKAP HI All-Sky Survey, known as WALLABY, is a shallow 3 Pi survey (z = 0 - 0.26) w… ▽ More The Australian Square Kilometer Array Pathfinder (ASKAP) will revolutionise our knowledge of gas-rich galaxies in the Universe. Here we present predictions for two proposed extragalactic ASKAP neutral hydrogen (HI) emission-line surveys, based on semi-analytic models applied to cosmological N-body simulations. The ASKAP HI All-Sky Survey, known as WALLABY, is a shallow 3 Pi survey (z = 0 - 0.26) which will probe the mass and dynamics of over 600,000 galaxies. A much deeper small-area HI survey, called DINGO, aims to trace the evolution of HI from z = 0 - 0.43, a cosmological volume of 40 million Mpc^3, detecting potentially 100,000 galaxies. The high-sensitivity 30 antenna ASKAP core (diameter ~2 km) will provide an angular resolution of 30 arcsec (at z=0). Our simulations show that the majority of galaxies detected in WALLABY (87.5%) will be resolved. About 5000 galaxies will be well resolved, i.e. more than five beams (2.5 arcmin) across the major axis, enabling kinematic studies of their gaseous disks. This number would rise to 160,000 galaxies if all 36 ASKAP antennas could be used; the additional six antennas provide baselines up to 6 km, resulting in an angular resolution of 10 arcsec. For DINGO this increased resolution is highly desirable to minimise source confusion; reducing confusion rates from a maximum of 10% of sources at the survey edge to 3%. We estimate that the sources detected by WALLABY and DINGO will span four orders of magnitude in total halo mass (from 10^{11} to 10^{15} Msol) and nearly seven orders of magnitude in stellar mass (from 10^{5} to 10^{12} Msol), allowing us to investigate the process of galaxy formation across the last four billion years. △ Less

Submitted 4 October, 2012; v1 submitted 28 August, 2012; originally announced August 2012.

Comments: 21 pages, accepted for publication in MNRAS, minor updates to published version and fixed links. Movies and images available at http://ict.icrar.org/store/Movies/Duffy12c/