Search

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., 5 pages, 3 figures, accepted to MNRAS Letters, open access version see https://doi.org/10.1093/mnrasl/slad031

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, Comment: 23 pages, 13 figures, Accepted for publication in ApJ

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., 19 pages, 14 figures, accepted for publication in MNRAS

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., 20 pages, 14 figures, accepted for publication in MNRAS

We introduce the Uchuu suite of large high-resolution cosmological N-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion (12 8003) dark matter particles in a box of side-length 2.0 h-1, Gpc$, with particle mass of 3.27 × 108, h-1, M⊙. The highest resolution simulation, Shin-Uchuu, consists of 262 billion (64003) particles in a box of side-length 140 h-1, Mpc, with particle mass of 8.97 × 105, h-1, M⊙. Combining these simulations, we can follow the evolution of dark matter haloes and subhaloes spanning those hosting dwarf galaxies to massive galaxy clusters across an unprecedented volume. In this first paper, we present basic statistics, dark matter power spectra, and the halo and subhalo mass functions, which demonstrate the wide dynamic range and superb statistics of the Uchuu suite. From an analysis of the evolution of the power spectra, we conclude that our simulations remain accurate from the baryon acoustic oscillation scale down to the very small. We also provide parameters of a mass-concentration model, which describes the evolution of halo concentration and reproduces our simulation data to within 5 per cent for haloes with masses spanning nearly eight orders of magnitude at redshift 0 ≤ z ≤ 14. There is an upturn in the mass-concentration relation for the population of all haloes and of relaxed haloes at z ≳ 0.5, whereas no upturn is detected at z < 0.5. We make publicly available various N-body products as part of Uchuu Data Release 1 on the Skies & Universes site. Future releases will include gravitational lensing maps and mock galaxy, X-ray cluster, and active galactic nucleus catalogues. © 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., The Uchuu simulations were carried out on Aterui II supercomputer at Center for Computational Astrophysics, CfCA, of National Astronomical Observatory of Japan, and the K computer at the RIKEN Advanced Institute for Computational Science (Proposal numbers hp180180, hp190161). The Uchuu DR1 effort has made use of the skun6@IAA facility managed by the IAA-CSIC in Spain, this equipment was funded by the Spanish MICINN EU-FEDER infrastructure grant EQC2018-004366-P. The skun@IAA_RedIRIS server was funded by the MICINN grant AYA2014-60641-C2-1-P. The numerical analysis was partially carried out on XC40 at the Yukawa Institute Computer Facility in Kyoto University. TI has been supported by MEXT via the ‘Priority Issue on Post-K computer’ (Elucidation of the Fundamental Laws and Evolution of the Universe), JICFUS, and MEXT via the ‘Program for Promoting Researches on the Supercomputer Fugaku’ (Toward a unified view of the universe: from large scale structures to planets, proposal numbers hp200124). TI thanks the support by MEXT/JSPS KAKENHI Grant Numbers JP17H04828, JP18H04337, JP19KK0344, and JP20H05245. FP, AK, and DM thank the support of the Spanish Ministry of Science and Innovation funding grant PGC2018-101931-B-I00. Parts of this research were conducted under the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. EJ acknowledges financial support from CNRS. FP and EJ want to thank a French–Spanish international collaboration grant from CNRS and CSIC. CVM acknowledges support from FONDECYT through grant 3200918, and he also acknowledges previous support from the Max Planck Society through a Partner Group grant. SAC acknowledges funding from Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET, PIP-0387), Agencia Nacional de Promoción de la Investigación, el Desarrollo Tecnológico y la Innovación (Agencia I+D+i, PICT-2018-03743), and Universidad Nacional de La Plata (G11-150), Argentina., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

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., Comment: 16 pages, 7 figures. Accepted for publication in MNRAS

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., 12 pages, 8 figure, accepted for publication in MNRAS

Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher fraction of active galactic nuclei (AGNs) in interacting and merging galaxies than in 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 AGNs 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.16}^{+0.21}$ and ${1.00}_{-0.31}^{+0.58}$ for X-ray and IR-selected AGNs, 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....

We demonstrate the potential of an emulator-based approach to analyzing galaxy formation models in the domain where constraining data is limited. We have applied the open-source Python package PRISM to the galaxy formation model Meraxes. Meraxes is a semi-analytic model, purposefully built to study the growth of galaxies during the Epoch of Reionization (EoR). Constraining such models is however complicated by the scarcity of observational data in the EoR. PRISM's ability to rapidly construct accurate approximations of complex scientific models using minimal data is therefore key to performing this analysis well. This paper provides an overview of our analysis of Meraxes using measurements of galaxy stellar mass densities; luminosity functions; and color-magnitude relations. We demonstrate the power of using PRISM instead of a full Bayesian analysis when dealing with highly correlated model parameters and a scarce set of observational data. Our results show that the various observational data sets constrain Meraxes differently and do not necessarily agree with each other, signifying the importance of using multiple observational data types when constraining such models. Furthermore, we show that PRISM can detect when model parameters are too correlated or cannot be constrained effectively. We conclude that a mixture of different observational data types, even when they are scarce or inaccurate, is a priority for understanding galaxy formation and that emulation frameworks like PRISM can guide the selection of such data., Comment: 19 pages, 120 figures, accepted for publication in ApJS

We build a theoretical picture of how the light from galaxies evolves across cosmic time. In particular, we predict the evolution of the galaxy spectral energy distribution (SED) by carefully integrating the star formation and metal enrichment histories of semi-analytic model (SAM) galaxies and combining these with stellar population synthesis models which we call mentari. Our SAM combines prescriptions to model the interplay between gas accretion, star formation, feedback process, and chemical enrichment in galaxy evolution. From this, the SED of any simulated galaxy at any point in its history can be constructed and compared with telescope data to reverse engineer the various physical processes that may have led to a particular set of observations. The synthetic SEDs of millions of simulated galaxies from mentari can cover wavelengths from the far UV to infrared, and thus can tell a near complete story of the history of galaxy evolution. \keywords{galaxies: evolution - galaxies: stellar content - galaxies.}, Comment: 5 pages, 3 figures, Proceedings of IAU symposium #341: Challenges in Panchromatic Modelling with Next Generation Facilities

The exact role of quasars during the Epoch of Reionization remains uncertain. With consensus leaning towards quasars producing a negligible amount of ionizing photons, we pose an alternate question: Can quasars indirectly contribute to reionization by allowing ionizing photons from stars to escape more easily? Using the Semi-Analytic Galaxy Evolution model to evolve a galaxy population through cosmic time, we construct an idealized scenario in which the escape fraction of stellar ionizing photons ($f_\mathrm{esc}$) is boosted following quasar wind events, potentially for several dynamical times. We find that under this scenario, the mean value of $f_\mathrm{esc}$ as a function of galaxy stellar mass peaks for intermediate mass galaxies. This mass dependence will have consequences for the 21cm power spectrum, enhancing power at small scales and suppressing it at large scales. This hints that whilst quasars may not directly contribute to the ionizing photon budget, they could influence reionization indirectly by altering the topology of ionized regions., Comment: 5 pages, 3 figures

Virtual Observatories (VOs) are online hubs of scientific knowledge. They encompass a collection of platforms dedicated to the storage and dissemination of astronomical data, from simple data archives to e-research platforms offering advanced tools for data exploration and analysis. Whilst the more mature platforms within VOs primarily serve the observational community, there are also services fulfilling a similar role for theoretical data. Scientific visualization can be an effective tool for analysis and exploration of datasets made accessible through web platforms for theoretical data, which often contain spatial dimensions and properties inherently suitable for visualization via e.g. mock imaging in 2d or volume rendering in 3d. We analyze the current state of 3d visualization for big theoretical astronomical datasets through scientific web portals and virtual observatory services. We discuss some of the challenges for interactive 3d visualization and how it can augment the workflow of users in a virtual observatory context. Finally we showcase a lightweight client-server visualization tool for particle-based datasets allowing quantitative visualization via data filtering, highlighting two example use cases within the Theoretical Astrophysical Observatory., Comment: 10 Pages, 13 Figures, Accepted for Publication in Monthly Notices of the Royal Astronomical Society

We have built PRISM, a "Probabilistic Regression Instrument for Simulating Models". PRISM uses the Bayes linear approach and history matching to construct an approximation ('emulator') of any given model, by combining limited model evaluations with advanced regression techniques, covariances and probability calculations. It is designed to easily facilitate and enhance existing Markov chain Monte Carlo (MCMC) methods by restricting plausible regions and exploring parameter space efficiently. However, PRISM can additionally be used as a standalone alternative to MCMC for model analysis, providing insight into the behavior of complex scientific models. With PRISM, the time spent on evaluating a model is minimized, providing developers with an advanced model analysis for a fraction of the time required by more traditional methods. This paper provides an overview of the different techniques and algorithms that are used within PRISM. We demonstrate the advantage of using the Bayes linear approach over a full Bayesian analysis when analyzing complex models. Our results show how much information can be captured by PRISM and how one can combine it with MCMC methods to significantly speed up calibration processes (>15 times faster). PRISM is an open-source Python package that is available under the BSD 3-Clause License (BSD-3) at https://github.com/1313e/PRISM and hosted at https://prism-tool.readthedocs.io. PRISM has also been reviewed by "The Journal of Open Source Software" (https://doi.org/10.21105/joss.01229)., 28 pages, 13 figures, 1 table. Updated to reflect some changes made to published version

We use three semi-analytical models (SAMs) of galaxy formation and evolution run on the same 1 h(-1) Gpc MultiDark Planck2 cosmological simulation to investigate the properties of [O II] emission line galaxies at redshift z similar to 1. We compare model predictions with different observational data sets, including DEEP2-FIREFLY galaxies with absolute magnitudes. We estimate the [O II] luminosity (L[O II]) of our model galaxies using the public code GET EMLINES, which ideally assumes as input the instantaneous star formation rates (SFRs). This property is only available in one of the SAMs under consideration, while the others provide average SFRs, as most models do. We study the feasibility of inferring galaxies' L[O II] from average SFRs in post-processing. We find that the result is accurate for model galaxies with dust attenuated L[O II] less than or similar to 10(42.2) erg s(-1) (< 5 per cent discrepancy). The galaxy properties that correlate the most with the model L[O II] are the SFR and the observed-frame u and g broad-band magnitudes. Such correlations have r-values above 0.64 and a dispersion that varies with L[O II]. We fit these correlations with simple linear relations and use them as proxies for L[O II], together with an observational conversion that depends on SFR and metallicity. These proxies result in [O II] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used to derive L[O II]. The amplitude of the clustering of model galaxies with L[O II] > 10(40.4) erg s(-1) remains overall unchanged on scales above 1 h(-1) Mpc, independently of the L[O II] computation. © 2020 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, GFandVGPacknowledge support from the University of Portsmouth through the Dennis Sciama Fellowship award. VGP acknowledges support from the European Research Council under the European Union's Horizon 2020 research and innovation programme (grant agreement No 769130). DS is funded by the Spanish Ministry of Economy and Competitiveness (MINECO) under the 2014 Severo Ochoa Predoctoral Training Programme. DS also wants to thank the Mamua Cafe Bar-team for their kind (g)astronomical support. DS and FP acknowledge funding support from the MINECO grant AYA2014-60641-C2-1-P. AO acknowledges support from the Spanish Ministerio de Economia y Competitividad (MINECO) project No. AYA2015-66211-C2-P-2, and funding from the European Union Horizon 2020 research and innovation programme under grant agreement No. 734374. SAC acknowledges funding from Consejo Nacional de Investigaciones Cientificas y Tecnicas (CONICET, PIP0387), Agencia Nacional de Promocion Cientifica y Tecnologica (ANPCyT, PICT-2013-0317), and Universidad Nacional de La Plata (11-G124 and 11-G150), Argentina. CVMacknowledges CONICET, Argentina, for the supporting fellowship. AK is supported by the MINECO and the Fondo Europeo de Desarrollo Regional (FEDER, UE) in Spain through grant AYA2015-63810-P as well as by the MICIU/FEDER through grant number PGC2018-094975-C21. He further acknowledges support from the Spanish Red Consolider MultiDark FPA2017-90566-REDC and thanks Christopher Cross for sailing. ARHS acknowledges receipt of the Jim Buckee Fellowship at ICRAR-UWA. GF, VGP and DS wish to thank La Plata Astronomical Observatory for hosting the MultiDark Galaxies workshop in September 2016, during which this work was started. The authors thank the FIREFLY team and the anonymous referee for providing insightful comments. The analysis of DEEP2 data using the FIREFLY code was done on the Sciama High Performance Compute cluster which is supported by the ICG, SEPNet and the University of Portsmouth (UK). The CosmoSim database used in this paper is a service by the Leibniz-Institute for Astrophysics Potsdam (AIP). The MultiDark database was developed in cooperation with the Spanish MultiDark Consolider Project CSD200900064. The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. (www.gauss-centre.eu) and the Partnership for Advanced Supercomputing in Europe (PRACE, www.praceri.eu) for funding the MultiDark simulation project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (LRZ, www.lrz.de).The authors thank New Mexico State University (USA) and Instituto de Astrof ' isica de Andalucia CSIC (Spain) for hosting the SKIES& UNIVERSES database for cosmological simulation products. This work has benefited from the publicly available software tools and packages: MATPLOTLIB14 2012-2016 (Hunter 2007); PYTHON SOFTWARE FOUNDATION15 19902017, version 2.7., PYTHONBREW16; we use whenever possible in this work a colour-blind friendly colour palette17 for our plots.

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., 15 pages, 8 figures. Accepted for publication in MNRAS

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., 15 pages, 9 figures (3 new figures), accepted on MNRAS

We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), transition galaxies and quiescent galaxies. We study their $UVI$ ($U-V$ versus $V-I$) color gradients to infer their sSFR gradients out to twice effective radii. We show that on average both star-forming and transition galaxies at all masses are not fully quenched at any radii, whereas quiescent galaxies are fully quenched at all radii. We find that at low masses ($M_{\ast} = 10^{9}-10^{10}M_{\odot}$) SFGs both above and below the SFMS ridge generally have flat sSFR profiles, whereas the transition galaxies at the same masses generally have sSFRs that are more suppressed in their outskirts. In contrast, at high masses ($M_{\ast} > 10^{10.5}M_{\odot}$), SFGs above and below the SFMS ridge and transition galaxies generally have varying degrees of more centrally-suppressed sSFRs relative to their outskirts. These findings indicate that at $z\sim~0.5-1.0$ the main galaxy quenching mode depends on its already formed stellar mass, exhibiting a transition from "the outside-in" at $M_{\ast} \leq 10^{10}M_{\odot}$ to "the inside-out" at $M_{\ast} > 10^{10.5}M_{\odot}$. In other words, our findings support that internal processes dominate the quenching of massive galaxies, whereas external processes dominate the quenching of low-mass galaxies., Comment: 19 pages, 12 figures, 2 tables, accepted for publication in ApJ

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., 11 pages, 7 figures. Accepted for publication in ApJ

The WiggleZ Dark Energy Survey measured the redshifts of over 200,000 UV-selected (NUV, Comment: Catalogue available at MNRAS (DOI link below) and also at: https://espace.library.uq.edu.au/view/UQ:3e43575

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 h-1 Gpc)3 and mass resolution of 1.5 × 109 h-1 M⊙) 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 ongoing and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available., Facultad de Ciencias Astronómicas y Geofísicas, Instituto de Astrofísica de La Plata

The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts., Comment: 18 pages, 11 figures, 3tables; accepted for publication in MNRAS

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 1013h−1M⊙ 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., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavily obscured AGN are the AGN most likely to be merger-driven. To better sample this population of AGN, we turn to infrared selection in the CANDELS/COSMOS field. Compared to their lower-luminosity and less obscured X-ray-only counterparts, IR-only AGN (luminous, heavily obscured AGN) are more likely to be classified as either irregular (50$^{+12}_{-12}$% vs. 9$^{+5}_{-2}$%) or asymmetric (69$^{+9}_{-13}$% vs. 17$^{+6}_{-4}$%) and are less likely to have a spheroidal component (31$^{+13}_{-9}$% vs. 77$^{+4}_{-6}$%). Furthermore, IR-only AGN are also significantly more likely than X-ray-only AGN (75$^{+8}_{-13}$% vs. 31$^{+6}_{-6}$%) to be classified either as interacting or merging in a way that significantly disturbs the host galaxy or disturbed though not clearly interacting or merging, which potentially represents the late stages of a major merger. This suggests that while major mergers may not contribute significantly to the fueling of Seyfert luminosity AGN, interactions appear to play a more dominant role in the triggering and fueling of high-luminosity heavily obscured AGN., 13 pages, 7 figures, 1 table, Accepted for publication in the Astrophysical Journal

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., 14 pages, 9 Figures, Accepted MNRAS

With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies., 14 pages, 8 figures, 3 tables. Submitted to MNRAS

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., 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

Successfully reproducing the galaxy luminosity function (LF) and the bimodality in the galaxy distribution requires a mechanism that can truncate star formation in massive haloes. Current models of galaxy formation consider two such truncation mechanisms: strangulation, which acts on satellite galaxies, and active galactic nucleus (AGN) feedback, which predominantly affects central galaxies. The efficiencies of these processes set the blue fraction of galaxies, fblue(L, M), as a function of galaxy luminosity, L, and halo mass, M. In this paper, we use a galaxy group catalogue extracted from the Sloan Digital Sky Survey (SDSS) to determine fblue(L, M). To demonstrate the potential power of these data as a benchmark for galaxy formation models, we compare the results to the semi-analytical model for galaxy formation of Croton et al. Although this model accurately fits the global statistics of the galaxy population, as well as the shape of the conditional LF, there are significant discrepancies when the blue fraction of galaxies as a function of mass and luminosity is compared between the observations and the model. In particular, the model predicts (i) too many faint satellites in massive haloes, (ii) a blue fraction of satellites that is much too low, and (iii) a blue fraction of centrals that is too high and with an inverted luminosity dependence. In the same order, we argue that these discrepancies owe to (i) the neglect of tidal stripping in the semi-analytical model, (ii) the oversimplified treatment of strangulation, and (iii) improper modelling of dust extinction and/or AGN feedback. The data presented here will prove useful to test and calibrate future models of galaxy formation and, in particular, to discriminate between various models for AGN feedback and other star formation truncation mechanisms

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., Comment: 10 pages, 8 figures. Accepted for publication in ApJ; added references

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., 8 pages, 4 figures. Accepted for publication in ApJ

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 Λ cold dark matter 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 two-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 h(−1) Mpc. 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., Monthly Notices of the Royal Astronomical Society, 469 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We use the EAGLE cosmological hydrodynamic simulation suite to study the specific angular momentum of galaxies, $j$, with the aims of (i) investigating the physical causes behind the wide range of $j$ at fixed mass and (ii) examining whether simple, theoretical models can explain the seemingly complex and non-linear nature of the evolution of $j$. We find that $j$ of the stars, $j_{\rm stars}$, and baryons, $j_{\rm bar}$, are strongly correlated with stellar and baryon mass, respectively, with the scatter being highly correlated with morphological proxies such as gas fraction, stellar concentration, (u-r) intrinsic colour, stellar age and the ratio of circular velocity to velocity dispersion. We compare with available observations at $z=0$ and find excellent agreement. We find that $j_{\rm bar}$ follows the theoretical expectation of an isothermal collapsing halo under conservation of specific angular momentum to within $\approx 50$%, while the subsample of rotation-supported galaxies are equally well described by a simple model in which the disk angular momentum is just enough to maintain marginally stable disks. We extracted evolutionary tracks of the stellar spin parameter of EAGLE galaxies and found that the fate of their $j_{\rm stars}$ at $z=0$ depends sensitively on their star formation and merger histories. From these tracks, we identified two distinct physical channels behind low $j_{\rm stars}$ galaxies at $z=0$: (i) galaxy mergers, and (ii) early star formation quenching. The latter can produce galaxies with low $j_{\rm stars}$ and early-type morphologies even in the absence of mergers., Comment: 22 pages (16 without appendices). Accepted for publication in MNRAS. The only difference with previous version is a shorter abstract (to comply with the 250 words limit in MNRAS)

We present an analysis of the pairwise velocity statistics from a suite of cosmological N-body simulations describing the "Running Friedmann-Lema\^itre-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends {\Lambda}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 {\sigma}_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 {\Lambda}CDM which could be readily exploited by future cosmological observations to test and further constrain our understanding of dark energy., Comment: 10 pages, 4 figures. arXiv admin note: text overlap with arXiv:1605.00714

We use EAGLE to quantify the effect galaxy mergers have on the stellar specific angular momentum of galaxies, $j_{\rm stars}$. We split mergers into: dry (gas-poor)/wet (gas-rich), major/minor, and different spin alignments and orbital parameters. Wet (dry) mergers have an average neutral gas-to-stellar mass ratio of $1.1$ ($0.02$), while major (minor) mergers are those with stellar mass ratios $\ge 0.3$ ($0.1-0.3$). We correlate the positions of galaxies in the $j_{\rm stars}$-stellar mass plane at $z=0$ with their merger history, and find that galaxies of low spins suffered dry mergers, while galaxies of normal/high spins suffered predominantly wet mergers, if any. The radial $j_{\rm stars}$ profiles of galaxies that went through dry mergers are deficient by $\approx 0.3$~dex at $r\lesssim 10\,r_{50}$ (with $r_{50}$ being the half-stellar mass radius), compared to galaxies that went through wet mergers. Studying the merger remnants reveals that dry mergers reduce $j_{\rm stars}$ by $\approx 30$\%, while wet mergers increase it by $\approx 10$\%, on average. The latter is connected to the build-up of the bulge by newly formed stars of high rotational speed. Moving from minor to major mergers accentuates these effects. When the spin vectors of the galaxies prior to the dry merger are misaligned, $j_{\rm stars}$ decreases to a greater magnitude, while in wet mergers co-rotation and high orbital angular momentum efficiently spun-up galaxies. We predict what would be the observational signatures in the $j_{\rm stars}$ profiles driven by dry mergers: (i) shallow radial profiles and (ii) profiles that rise beyond $\approx 10\,r_{50}$, both of which are significantly different from spiral galaxies., Accepted for publication in MNRAS. Compared to previous revision, Fig. 6 and A1 are new, and section 3.3 was shortened

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., 15 pages, 12 figures. Accepted for publication in MNRAS

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., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

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., Comment: 15 pages, 9 figures, submitted to MNRAS. arXiv admin note: text overlap with arXiv:1201.5918

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., 9 pages, 8 figures, Accepted to ApJ

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., Comment: 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

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., 12 pages, 9 figures, accepted to MNRAS

Higher order statistics are a useful and complementary tool for measuring the clustering of galaxies, containing information on the non-Gaussian evolution and morphology of large-scale structure in the Universe. In this work we present measurements of the three-point correlation function (3PCF) for 187 000 galaxies in the WiggleZ spectroscopic galaxy survey. We explore the WiggleZ 3PCF scale and shape dependence at three different epochs z = 0.35, 0.55 and 0.68, the highest redshifts where these measurements have been made to date. Using N-body simulations to predict the clustering of dark matter, we constrain the linear and non-linear bias parameters of WiggleZ galaxies with respect to dark matter, and marginalize over them to obtain constraints on σ8(z), the variance of perturbations on a scale of 8 h^−1 Mpc and its evolution with redshift. These measurements of σ_8(z), which have 10–20 per cent accuracies, are consistent with the predictions of the Λ cold dark matter concordance cosmology and test this model in a new way.

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 $\Omega_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$, $\Omega_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 $\Omega_Mh^2$ and $z_{eq}$ of $\pm$10% (assuming a prior $N_{eff}=3$) and on $N_{eff}$ of (+78/-56)% (assuming a prior $\Omega_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, \Omega_Mh^2, N_eff)$ of (3,4,+17/-21)%. This represents a precision approaching our forecasts for the Planck Surveyor., Comment: 12 pages, 7 figures, accepted for publication in MNRAS

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., Accepted in MNRAS. 21 pages, 18 figures. Replacement includes fixed typos and reference updates

We present a suite of cosmological N-body simulations describing the "Running Friedmann-Lema{\"i}tre-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends {\Lambda}CDM with a time-evolving vacuum density, {\Lambda}(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 {\Lambda}CDM, which could be readily exploited by future cosmological observations to test and further constrain the model., Comment: 14 pages, 7 figures

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., Published in MNRAS, 19 pages, 20 figures (both including appendices)

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., Accepted on MNRAS; 24 pages, 11 figures

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, Comment: 16 pages, 10 figures, resubmitted to ApJ after addressing the Referee's comments

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., Comment: 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

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., 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

We have used deep, HST, near-IR imaging to study the morphological properties of the most massive galaxies at high z, modelling the WFC3/IR H-band images of the ~200 galaxies in the CANDELS-UDS field with 1 10^11 M_sun. We have used both single-Sersic and bulge+disk models, have investigated the errors/biases introduced by uncertainties in the background and the PSF, and have obtained formally-acceptable model fits to >90% of the galaxies. Our results indicate that these massive galaxies at 1 2 the compact bulges display effective radii a factor ~4 smaller than local ellipticals of comparable mass. These trends appear to extend to the bulge components of disk-dominated galaxies, and vice versa. We also find that, while such massive galaxies at low z are bulge-dominated, at 1 2 they are mostly disk-dominated. The majority of the disk-dominated galaxies are actively forming stars, but this is also true for many of the bulge-dominated systems. Interestingly, however, while most of the quiescent galaxies are bulge-dominated, we find that a significant fraction (25-40%) of the most quiescent galaxies have disk-dominated morphologies. Thus, while our results show that the massive galaxy population is undergoing dramatic changes at this crucial epoch, they also suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies into present-day giant ellipticals.