Search

Substantial populations of massive quiescent galaxies at $z\ge3$ challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at $z>3$. As a result, studies report that as much as 70\% of quiescent galaxies at $z>3$ may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at $3\lesssim z \lesssim 6$ using JWST NIRCam imaging data. We use empirically-constrained galaxy SED templates to define a region in the $F277W-F444W$ vs. $F150W-F277W$ color plane that captures quiescent galaxies at $z>3$. We apply this color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and identify 44 candidate $z\gtrsim3$ quiescent galaxies. Over half of these sources are newly discovered and, on average, exhibit specific star formation rates of post-starburst galaxies. We derive volume density estimates of $n\sim1-4\times10^{-5}$\,Mpc$^{-3}$ at $3< z <5$, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam's wavelength coverage and sensitivity, this technique provides an efficient tool to search for large samples of these rare galaxies., Comment: 12 pages, 4 figures, 2 tables, submitted to ApJ

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

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

The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting., Comment: 25 pages, 11 figures. Accepted for publication in ApJ

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

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., Comment: 24 pages, 15 figures, submitted to MNRAS, comments are welcome

We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, with each galaxy in the sample classified by three different individuals. We also measure quantitative morphologies using the publicly available codes across all seven NIRCam filters. Using these measurements, we present the fraction of galaxies of each morphological type as a function of redshift. Overall, we find that galaxies at z>3 have a wide diversity of morphologies. Galaxies with disks make up a total of 60\% of galaxies at z=3 and this fraction drops to ~30% at z=6-9, while galaxies with spheroids make up ~30-40% across the whole redshift range and pure spheroids with no evidence for disks or irregular features make up ~20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (~40-50%), while those that are purely irregular increases from ~12% to ~20% at z>4.5. We note that these are apparent fractions as many selection effects impact the visibility of morphological features at high redshift. The distributions of S\'ersic index, size, and axis ratios show significant differences between the morphological groups. Spheroid Only galaxies have a higher S\'ersic index, smaller size, and higher axis ratio than Disk/Irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study and further work with large samples at higher redshift is needed to quantify when these features first formed., Comment: Accepted for publication to ApJL, 24 pages, 14 figures

Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at $z>1$ based on rest-frame NIR F444W images of high resolution (~1.3 kpc at z ~ 1-3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at $z>1$ with spectroscopic redshifts, including the two highest redshift bars at ~2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at $z$ ~ 1.1-2.3 presented in our study have projected semi-major axes of ~2.9-4.3 kpc and projected ellipticities of ~0.41-0.53 in the rest-frame NIR. The barred host galaxies have stellar masses ~ $ 1 \times 10^{10}$ to $2 \times 10^{11}$ $M_{\odot}$, star formation rates of ~ 21-295 $M_{\odot}$ yr$^{-1}$, and several have potential nearby companions. Our finding of bars at $z$ ~1.1-2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8-10 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by z ~ 0., Comment: 16 pages, 5 figures. Accepted for Publication in Astrophysical Journal Letters

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

We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state., Comment: Submitted to ApJ; 12 Pages, 6 Figures

Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the random forests using the merger history information for the simulated galaxies available from IllustrisTNG. The random forests correctly classify $\sim60\%$ of non-merging and merging galaxies across $0.5 < z < 4.0$. Rest-frame asymmetry parameters appear more important for lower redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the random forest classifications match with theoretical Illustris predictions, but are underestimated by a factor of $\sim 0.5$., Comment: 21 pages, 14 figures, submitted to ApJ

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$\nu^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$\nu^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST)., Comment: 18 pages, 16 figures, accepted for publication in MNRAS, Figs. 3, 4, 11, 13, and 14 revised, Figs. 14 and 16 in the original manuscript moved to the Appendix, All Uchuu-$\nu^2$GC data used in this paper are publicly available on our Skies & Universes website https://www.skiesanduniverses.org/Simulations/Uchuu/GalaxyCatalogues/

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

We investigate the structural properties, distribution and abundance of LCDM dark matter subhaloes using the Phi-4096 and Uchuu suite of N-body cosmological simulations. Thanks to the combination of their large volume, high mass resolution and superb statistics, we are able to quantify -- for the first time consistently over more than seven decades in ratio of subhalo-to-host-halo mass -- dependencies of subhalo properties with mass, maximum circular velocity, Vmax, host halo mass and distance to host halo centre. We also dissect the evolution of these dependencies over cosmic time. We provide accurate fits for the subhalo mass and velocity functions, both exhibiting decreasing power-law slopes in the expected range of values and with no significant dependence on redshift. We also find subhalo abundance to depend weakly on host halo mass. We explore the distribution of subhaloes within their hosts and its evolution over cosmic time for subhaloes located as deep as ~0.1 per cent of the host virial radius. Subhalo structural properties are codified via a concentration parameter, cV, that does not depend on any specific, pre-defined density profile and relies only on Vmax. We derive the cV-Vmax relation in the range 7-1500 km/s and find an important dependence on distance of the subhalo to the host halo centre, as already described in Molin\'e et al. (2017). Interestingly, we also find subhaloes of the same mass to be significantly more concentrated into more massive hosts. Finally, we investigate the redshift evolution of cV, and provide accurate fits that take into account all mentioned dependencies. Our results offer an unprecedented detailed characterization of the subhalo population, consistent over a wide range of subhalo and host halo masses, as well as cosmic times. Our work enables precision work in any future research involving dark matter halo substructure., Comment: 17 pages, 15 figures, 5 tables. Comments welcome!

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

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

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 $\Delta V <5000$ km s$^{-1}$) at $0.5

We introduce the Uchuu suite of large high-resolution cosmological $N$-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion ($12800^3$) dark matter particles in a box of side-length 2.0 Gpc/h, with particle mass $3.27 \times 10^{8}$ Msun/h. The highest resolution simulation, Shin-Uchuu, consists of 262 billion ($6400^3$) particles in a box of side-length 140 Mpc/h, with particle mass $8.97 \times 10^{5}$ Msun/h. Combining these simulations we can follow the evolution of dark matter halos and subhalos 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 halos with masses spanning nearly eight orders of magnitude at redshift 00.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 nuclei catalogues., Comment: 22 pages, 14 figures, accepted by MNRAS. We release various $N$-body products as data release 1 on http://skiesanduniverses.org/ such as subsets of simulation particles, matter power spectra, halo/subhalo catalogues, and their merger trees

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

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

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

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

One of the most important parameters in characterizing the Epoch of Reionization, the escape fraction of ionizing photons, $f_\mathrm{esc}$, remains unconstrained both observationally and theoretically. With recent work highlighting the impact of galaxy-scale feedback on the instantaneous value of $f_\mathrm{esc}$, it is important to develop a model in which reionization is self-consistently coupled to galaxy evolution. In this work, we present such a model and explore how physically motivated functional forms of $f_\mathrm{esc}$ affect the evolution of ionized hydrogen within the intergalactic medium. Using the $21$cm power spectrum evolution, we investigate the likelihood of observationally distinguishing between a constant $f_\mathrm{esc}$ and other models that depend upon different forms of galaxy feedback. We find that changing the underlying connection between $f_\mathrm{esc}$ and galaxy feedback drastically alters the large-scale $21$cm power. The upcoming Square Kilometre Array Low Frequency instrument possesses the sensitivity to differentiate between our models at a fixed optical depth, requiring only $200$ hours of integration time focused on redshifts $z = 7.5-8.5$. Generalizing these results to account for a varying optical depth will require multiple $800$ hour observations spanning redshifts $z = 7-10$. This presents an exciting opportunity to observationally constrain one of the most elusive parameters during the Epoch of Reionization., Comment: 14 pages, 8 figures, 2 tables

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)., Comment: 28 pages, 13 figures, 1 table. Updated to reflect some changes made to published version

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., Comment: 20 pages, 8 figures, 7 tables. MNRAS published version

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

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

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

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

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-{\sigma}. 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., Comment: 20 pages, 14 figures (+5 supplementary figures in the Appendix), accepted for publication in MNRAS

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

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

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

We study galactic star-formation activity as a function of environment and stellar mass over 0.59 (9.5)$ at z=1.3 (2.0). This method, when applied to a mock catalog with the photometric-redshift precision ($\sigma_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., Comment: 29 pages, 15 figure, accepted for publication in ApJ

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

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

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

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 {\Lambda}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., Comment: 17 pages, 7 figures

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., Comment: Accepted for publication in MNRAS. Compared to previous revision, Fig. 6 and A1 are new, and section 3.3 was shortened

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

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

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

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

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