Search

We investigate the formation of infant globular cluster (GC) candidates in high-resolution cosmological simulations from the First Billion Years (FiBY) project. By analysing the evolution of the systems in the energy and angular momentum plane, we identify the redshift at which the infant GCs first became gravitationally bound, and we find evidence of radial infall of their gaseous and stellar components. The collapse appears to be driven by internal self-gravity, however, the initial trigger is sourced from the external environment. The phase space behaviour of the infant GCs also allows us to identify some characteristic groupings of objects. Such a classification based on internal properties appears to be reflected in the formation environment: GC candidates that belong to the same class are found in host galaxies of similar morphology, with the majority of the infant GCs located in clumpy, irregular proto-galaxies. Finally, through the inspection of two GC candidates that contain only stars by z = 6, we find that supernova feedback is the main physical mechanism behind their dearth of gas and that the systems subsequently respond with an approximately adiabatic expansion. Such infant GC candidates already resemble the GCs we currently observe in the local Universe., Comment: 17 pages, 12 Figures, Accepted by MNRAS

We computed the luminosity of simulated galaxies of the C-EAGLE project, a suite of 30 high-resolution zoom-in simulations of galaxy clusters based on the EAGLE simulation. The AB magnitudes are derived for different spectral bands, from ultraviolet to infrared, using the simple stellar population modeling based on the E-MILES stellar spectra library. We take into account obscuration due to dust in star forming regions and diffuse interstellar medium. The $g-r$ colour-stellar mass diagram, at z=0.1, presents a defined red sequence, reaching $g-r \simeq 0.8$, 0.05 dex redder than EAGLE at high masses, and a well populated blue cloud, when field galaxies are included. The clusters' inner regions are dominated by red-sequence galaxies at all masses, although a non-negligible amount of blue galaxies are still present. We adopt Bayesian inference to compute the clusters LFs, testing for statistical significance of both single and double Schechter functions. The multicolour LFs at z=0 show a knee luminosity that peaks in the infrared and increases with the cluster's mass. The faint-end is weakly dependent on colour and mass and shows an upturn in the optical, bounded between -1.25 and -1.39, just moderately steeper than the field. The simulations reproduce, within the observational errors, the spectroscopic LFs of the Hercules and Abell 85 clusters, including their faint end upturn. C-EAGLE LFs are in broad agreement with observed LFs taken from SDSS and XXL surveys, up to z=0.67, showing a rather flat faint end when the observational constrains are taken into account., 18 pages, 10 figures, 7 tables. Accepted for publication in MNRAS

We extend the state-of-the-art N-body code PKDGRAV3 with the inclusion of mesh-free gas hydrodynamics for cosmological simulations. Two new hydrodynamic solvers have been implemented, the mesh-less finite volume and mesh-less finite mass methods. The solvers manifestly conserve mass, momentum and energy, and have been validated with a wide range of standard test simulations, including cosmological simulations. We also describe improvements to PKDGRAV3 that have been implemented for performing hydrodynamic simulations. These changes have been made with efficiency and modularity in mind, and provide a solid base for the implementation of the required modules for galaxy formation and evolution physics and future porting to GPUs. The code is released in a public repository, together with the documentation and all the test simulations presented in this work., 18 pages, 14 figures; accepted for publication in MNRAS

Cosmological simulations are crucial tools in studying the Universe, but they typically do not directly match real observed structures. Constrained cosmological simulations, on the other hand, are designed to match the observed distribution of galaxies. Here we present constrained simulations based on spectroscopic surveys at a redshift of z~2.3, corresponding to an epoch of nearly 11 Gyrs ago. This allows us to 'fast-forward' the simulation to our present-day and study the evolution of observed cosmic structures self-consistently. We confirm that several previously-reported protoclusters will evolve into massive galaxy clusters by our present epoch, including the 'Hyperion' structure that we predict will collapse into a giant filamentary supercluster spanning 100 Megaparsecs. We also discover previously unknown protoclusters, with lower final masses than typically detectable by other methods, that nearly double the number of known protoclusters within this volume. Constrained simulations, applied to future high-redshift datasets, represents a unique opportunity for studying early structure formation and matching galaxy properties between high and low redshifts., Submitted: 10 November 2021; Accepted: 28 April 2022 in Nature Astronomy https://www.nature.com/articles/s41550-022-01693-0 32 pages, 9 Figures

We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, ${\gamma}_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of ${\gamma}_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of ${\gamma}_{\rm gNFW}$., Comment: 19 pages, 16 figures, accepted by MNRAS

Background: In Ireland, the paediatric emergency department (PED) is tasked with acute healthcare provision to children and adolescents under 16 years of age. The population > 15 years attending the PED remains undescribed. Aims: The objective of the study is to describe the presentations of patients > 15 years to the PED of a national tertiary academic paediatric hospital. Methods: A retrospective review of electronic records identifying all patients (> 15 years) who presented from January 2014 to December 2015. Patient demographics, presenting complaint, diagnosis, treatment and disposition were recorded. Results: A total of 71,082 patients attended during the study period; of whom, 426 (0.6%) patients were aged 16 to 61 years. Over half were identified as having clearly defined chronic complex conditions. Three hundred and forty-five (89%) patients were known to the hospital paediatric services and under the care of specialist with 131 (34%) patients admitted locally, all of whom had chronic complex conditions (CCC). There was a total of 44 (11%) patients who presented to the PED with de novo issues and had a variety of acute presenting complaints with minor injuries, syncope and chest pain being most common. Conclusion: This is the first Irish study demonstrating a significant population (> 15 years) with medical complexity that requires a suite of services in a paediatric hospital which is accessed in an unscheduled manner through an emergency department. We recommend further research to describe paediatric CCC attending Irish emergency and hospital services.

Rib fractures represent a substantial health burden. Chest injuries contribute to 25% of deaths after trauma and survivors can experience long-standing consequences, such as reduced functional capabilities and loss of employment. Over recent years, there has been an increase in the awareness of the importance of early identification, aggressive pain management and adequate safety netting for patients with chest injuries. Substandard management leads to increased rates of morbidity and mortality. The development of protocols in the emergency department (ED) for management of patients with chest wall injuries has demonstrated reduction of complication rates.Our aim was to develop an evidence-based, multidisciplinary chest injury pathway for the management of patients presenting with rib injury to our ED.Prior to implementation of the pathway in our department, only 39% of patients were documented as having received analgesia and only 7% of discharged patients had documented written verbal advice. There was no standardised method to perform regional anaesthetic blocks. Using quality improvement methods, we standardised imaging modality, risk stratification with a scoring system, analgesia with emphasis on regional anaesthesia blocks and disposition with information leaflets for those discharged.Implementation of the pathway increased rates of documented analgesia received from 39% to 70%. The number of regional anaesthetic blocks performed went from 0% to 60% and the number of patients receiving discharge advice went from 7% to 70%. Compliance of doctors and nurses with the pathway was 63%.Our previous audits showed substandard management of patients with chest injuries in our department. Through this quality improvement project, we were able to improve the quality of care provided to patients attending with rib fractures by increasing rate of analgesia received, regional blocks performed and discharge advice given.

We perform high-resolution cosmological hydrodynamic simulations to study the formation of the first galaxies that reach the masses of $10^{8-9}~h^{-1}~M_\odot$ at $z=9$. The resolution of the simulations is high enough to resolve minihaloes and allow us to successfully pursue the formation of multiple Population (Pop) III stars, their supernova (SN) explosions, resultant metal-enrichment of the inter-galactic medium (IGM) in the course of the build-up of the system. Metals are ejected into the IGM by multiple Pop III SNe, but some of the metal-enriched gas falls back onto the halo after $\gtrsim 100~\rm Myr$. The star formation history of the first galaxy depends sensitively on the initial mass function (IMF) of Pop III stars. The dominant stellar population transits from Pop III to Pop II at $z\sim 12-15$ in the case of power-law Pop III IMF, ${\rm d}n/{\rm d}M \propto M^{-2.35}$ with the mass range $10-500~M_\odot$. At $z\lesssim 12$, stars are stably formed in the first galaxies with a star formation rate of $\sim 10^{-3}$-$10^{-1}~M_\odot/{\rm yr}$. In contrast, for the case with a flat IMF, gas-deprived first galaxies form due to frequent Pop III pair-instability SNe, resulting in the suppression of subsequent Pop II star formation. In addition, we calculate UV continuum, Ly$\alpha$- and H$\alpha$-line fluxes from the first galaxies. We show that the James Webb Space Telescope will be able to detect both UV continuum, Ly$\alpha$ and H$\alpha$ line emission from first galaxies with halo mass $\gtrsim 10^{9}~M_\odot$ at $z \gtrsim 10$., Comment: 14 pages, 15 figures, accepted for publication in MNRAS

We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges' of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the `stellar splashback' feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (mu >> 32 mag/arcsec^2) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid., 12 pages, 9 figures. Accepted for publication in MNRAS

We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization., Comment: 24 pages, 19 figures, accepted for publication in MNRAS

Our location in the Milky Way provides an exceptional opportunity to gain insight on the galactic evolution processes, and complement the information inferred from observations of external galaxies. Since the Milky Way is a barred galaxy, the study of motions of individual stars in the bulge and disc is useful to understand the role of the bar. The Gaia mission enables such study by providing the most precise parallaxes and proper motions to date. In this theoretical work, we explore the effects of the bar on the distribution of higher-order moments --the skewness and kurtosis-- of the proper motions by confronting two simulated galaxies, one with a bar and one nearly axisymmetric, with observations from the latest Gaia data release (GaiaDR2). We introduce the code ASGAIA to account for observational errors of Gaia in the kinematical structures predicted by the numerical models. As a result, we find clear imprints of the bar in the skewness distribution of the longitudinal proper motion $\mu_\ell$ in GaiaDR2, as well as other features predicted for the next Gaia data releases., Comment: 11 pages, 9 figures, accepted in A&A

By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-density contours using the procedure of Montes \& Trujillo (2019), and find that the shape of the stellar mass distribution follows that of the total matter even more closely than observed, although their radial profiles differ substantially. The ratio between stellar and total matter density profiles in circular apertures, shows a slope close to -1, with a small dependence on the cluster's total mass. We propose an indirect method to calculate the halo mass and mass density profile from the radial profile of the intra-cluster stellar mass density., 7 pages, 5 figures, submitted to MNRAS

In the last decade observations have been able to probe the evolution of the galaxy luminosity function, in particular showing a variation of its faint-end with redshift. We employ the data of the Cluster-EAGLE project, a set of cosmological, hydrodynamical zoom-in simulations of 30 galaxy clusters, to study the evolution of the galaxy luminostity functions in clusters with redshift. We compile a catalogue of simulated galaxies’ luminosities in the SDSS bands using the E-MILES spectra database, and taking into account dust attenuation. Stacked luminosity functions present little evolution with redshift of the faint-end slope from z=3.5 to z=0, regardless of the cluster mass.

We investigate higher order symplectic integration strategies within Bayesian cosmic density field reconstruction methods. In particular, we study the fourth-order discretisation of Hamiltonian equations of motion (EoM). This is achieved by recursively applying the basic second-order leap-frog scheme (considering the single evaluation of the EoM) in a combination of even numbers of forward time integration steps with a single intermediate backward step. This largely reduces the number of evaluations and random gradient computations, as required in the usual second-order case for high-dimensional cases. We restrict this study to the lognormal-Poisson model, applied to a full volume halo catalogue in real space on a cubical mesh of 1250 $h^{-1}$ Mpc side and 256$^3$ cells. Hence, we neglect selection effects, redshift space distortions, and displacements. We note that those observational and cosmic evolution effects can be accounted for in subsequent Gibbs-sampling steps within the COSMIC BIRTH algorithm. We find that going from the usual second to fourth-order in the leap-frog scheme shortens the burn-in phase by a factor of at least $\sim30$. This implies that 75-90 independent samples are obtained while the fastest second-order method converges. After convergence, the correlation lengths indicate an improvement factor of about 3.0 fewer gradient computations for meshes of 256$^3$ cells. In the considered cosmological scenario, the traditional leap-frog scheme turns out to outperform higher order integration schemes only at lower dimensional problems, e.g. meshes with 64$^3$ cells. This gain in computational efficiency can help to go towards a full Bayesian analysis of the cosmological large-scale structure for upcoming galaxy surveys., 19 pages, 12 figures, 4 tables, accepted at MNRAS, additional robust mathematical argument supported by numerical tests with longer HMC chains and a solid statistical analysis

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of -15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%., Comment: 45 pages, 21 figures, accepted for publication in the Astrophysical Journal

We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes., 26 pages, 24 figures (16 pages, 14 figures without appendices); accepted by MNRAS (a few extra references added)

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter., 20 pages, 16 Figures. Highlights: Figs 4, 13, 15. To be submitted to MNRAS, comments welcome. Contact: lovell@hi.is

Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole (MBH) and their stellar mass (M∗). The origin of observed positive outliers from this relation with extremely high MBH (> 109M⊙) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations, designed to enable the study of high-mass galaxy formation and evolution in cluster environments. We find 69 MBH(M∗) outliers at z=0, defined as those with MBH107M⊙ and MBH/M∗ > 0.01. This paper focuses on a sample of five extreme outliers, that have been selected based on their MBH and M∗ values, which are comparable to the most recent estimates of observed positive outliers. This sample of five outliers, classified as 'black hole monster galaxies' (BMGs), was traced back in time to study their origin and evolution. In agreement with the results of previous simulations for lower mass MBH(M∗) outliers, we find that these galaxies became outliers due to a combination of their early formation times and tidal stripping. For BMGs with MBH 109 M⊙, major mergers (with a stellar mass ratio of μ > 0.25) at early times (z> 2) precede the rapid growth of their supermassive BHs. Furthermore, the scatter in the relation between MBH and stellar velocity dispersion, σ, correlates positively with the scatter in [Mg/Fe](σ). This indicates that the alpha enhancement of these galaxies, which is closely related to their star formation history, is related to the growth of their central BHs.

We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method., Additional results with respect to v1, new section and new figures added. 25 pages, 1 table 18 figures

Early-type galaxies are considered to be the end-products of massive galaxy formation. Optical spectroscopic studies reveal that massive early-type galaxies formed the bulk of their stars over short timescales (12), followed by passive evolution to the present. However, their optical spectra are insensitive to constrain small episodes of recent star formation, since they are dominated by old stars. Fortunately, this problem can be tackled in the ultraviolet range. While recent studies that make use of ultraviolet absorption lines have suggested the presence of young stars in a few early-type galaxies, the age and mass fractions of young stars and their dependence on galaxy mass, is unknown. Here we report a detailed study of these young stellar populations, from high-quality stacked spectra of 28,663 galaxies from the BOSS survey, analysing optical and ultraviolet absorption lines simultaneously. We find that residual star formation is ubiquitous in massive early-type galaxies, measuring average mass fractions of ~0.5% in young stars in the last 2 Gyr of their evolution. This fraction shows a decreasing trend with galaxy stellar mass, consistent with a down-sizing scenario. We also find that synthetic galaxies from state-of-the-art cosmological numerical simulations significantly overproduce both intermediate and young stellar populations. Therefore, our results pose stringent constraints on numerical simulations of galaxy formation., Comment: Nature Astronomy, December 16, 2019. 21 pages, 8 figures, 2 tables. The published version is available at the following URL: https://www.nature.com/articles/s41550-019-0955-0

For several decades it has been known that stellar bars in disc galaxies can be triggered by interactions, or by internal processes such as dynamical instabilities. In this work, we explore the differences between these two mechanisms using numerical simulations. We perform two groups of simulations based on isolated galaxies, one group in which a bar develops naturally, and another group in which the bar could not develop in isolation. The rest of the simulations recreate 1:1 coplanar fly-by interactions computed with the impulse approximation. The orbits we use for the interactions represent the fly-bys in groups or clusters of different masses accordingly to the velocity of the encounter. In the analysis we focus on bars' amplitude, size, pattern speed and their rotation parameter, ${\cal R}=R_{CR}/R_{bar}$. The latter is used to define fast (${\cal R}1.4$). Compared with equivalent isolated galaxies we find that bars affected or triggered by interactions: (i) remain in the slow regime for longer; (ii) are more boxy in face-on views; (iii) they host kinematically hotter discs. Within this set of simulations we do not see strong differences between retrograde or prograde fly-bys. We also show that slow interactions can trigger bar formation., 12 pages, 7 figures. Accepted for publication in MNRAS

Measurements of the redshift-space galaxy clustering have been a prolific source of cosmological information in recent years. Accurate covariance estimates are an essential step for the validation of galaxy clustering models of the redshift-space two-point statistics. Usually, only a limited set of accurate N-body simulations is available. Thus, assessing the data covariance is not possible or only leads to a noisy estimate. Further, relying on simulated realisations of the survey data means that tests of the cosmology dependence of the covariance are expensive. With these points in mind, this work presents a simple theoretical model for the linear covariance of anisotropic galaxy clustering observations with synthetic catalogues. Considering the Legendre moments (`multipoles') of the two-point statistics and projections into wide bins of the line-of-sight parameter (`clustering wedges'), we describe the modelling of the covariance for these anisotropic clustering measurements for galaxy samples with a trivial geometry in the case of a Gaussian approximation of the clustering likelihood. As main result of this paper, we give the explicit formulae for Fourier and configuration space covariance matrices. To validate our model, we create synthetic HOD galaxy catalogues by populating the haloes of an ensemble of large-volume N-body simulations. Using linear and non-linear input power spectra, we find very good agreement between the model predictions and the measurements on the synthetic catalogues in the quasi-linear regime., 17 pages, 16 figures, 3 tables; modified to match version accepted by MNRAS

We explored a suite of high-resolution cosmological simulations from the First Billion Years Project (FiBY) at $z \geq 6$. All substructures within the simulations have been identified with the SUBFIND algorithm. From our analysis, two distinct groups of objects emerge. We hypothesise that the substructures in the first group, which appear to have a high baryon fraction ($f_{\rm b} \geq 0.95$), are possible infant GC candidates. Objects belonging to the second group have a high stellar fraction ($f_{\rm star} \geq 0.95$) and show a potential resemblance to infant ultra-faint dwarf galaxies. The high baryon fraction objects identified in this study are characterised by a stellar content similar to the one observed in present-day GCs, but they still contain a high gas fraction ($f_{\rm gas} \sim 0.95$) and a relatively low amount of dark matter. They are compact, dense systems. Their sizes are consistent with recent estimates based on the first observations of possible proto-GCs at high redshifts. These types of infant GC candidates appear to be more massive and more abundant in massive host galaxies, indicating that the assembly of galaxies via mergers may play an important role in building several GC-host scaling relations. Specifically, we express the relation between the mass of the most massive infant GC and its host stellar mass as $\log(M_{\rm cl}) = (0.31\pm0.15)\log(M_{\rm *,gal} + (4.17\pm1.06)$. We also report a new relation between the most massive infant GC and the parent specific star formation rate of the form $\log(M_{\rm cl}) = (0.85\pm0.30)\log(sSFR) + \alpha$ that describes the data at both low and high redshift. Finally, we assess the present-day GC mass (GC number) -- halo mass relation offers a satisfactory description of the behaviour of our infant GC candidates at high redshift, suggesting that such a relation may be set at formation., Comment: 17 pages, 11 figures, accepted by A&A

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others., Comment: 23 pages, 11 figures. Replaced to match accepted MNRAS version. Included Kullback-Leibler divergence

Bars are common galactic structures in the local universe that play an important role in the secular evolution of galaxies, including the Milky Way. In particular, the velocity distribution of individual stars in our galaxy is useful to shed light on stellar dynamics, and provides information complementary to that inferred from the integrated light of external galaxies. However, since a wide variety of models reproduce the distribution of velocity and the velocity dispersion observed in the Milky Way, we look for signatures of the bar on higher-order moments of the line-of-sight velocity ($V_{los}$) distribution. We make use of two different numerical simulations --one that has developed a bar and one that remains nearly axisymmetric-- to compare them with observations in the latest APOGEE data release (SDSS DR14). This comparison reveals three interesting structures that support the notion that the Milky Way is a barred galaxy. A high skewness region found at positive longitudes constrains the orientation angle of the bar, and is incompatible with the orientation of the bar at $\ell=0^\circ$ proposed in previous studies. We also analyse the $V_{los}$ distributions in three regions, and introduce the Hellinger distance to quantify the differences among them. Our results show a strong non-Gaussian distribution both in the data and in the barred model, confirming the qualitative conclusions drawn from the velocity maps. In contrast to earlier work, we conclude it is possible to infer the presence of the bar from the kurtosis distribution., Accepted in MNRAS. 14 pages, 5 figures

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance., Comment: 20 pages, 16 figures, replaced to match accepted MNRAS version. Results on parameter errors changed

Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of volumes up to $(100\,\rm Mpc)^3$, while Hydrangea consists of $24$ cosmological simulations of galaxy clusters and their environment. Together they provide a statistically significant sample in the stellar mass range $10^{9.5}\,\rm M_{\odot}-10^{12.3}\,\rm M_{\odot}$, of $16,358$ galaxies. We construct IFS-like cubes and measure stellar spin parameters, $\lambda_{\rm R}$, and ellipticities, allowing us to classify galaxies into slow/fast rotators as in observations. The simulations display a primary dependence of $F_{\rm SR}$ on stellar mass, with a weak dependence on environment. At fixed stellar mass, satellite galaxies are more likely to be SRs than centrals. $F_{\rm SR}$ shows a dependence on halo mass at fixed stellar mass for central galaxies, while no such trend is seen for satellites. We find that $\approx 70$% of SRs at $z=0$ have experienced at least one merger with mass ratio $\ge 0.1$, with dry mergers being at least twice more common than wet mergers. Individual dry mergers tend to decrease $\lambda_{\rm R}$, while wet mergers mostly increase it. However, $30$% of SRs at $z=0$ have not experienced mergers, and those inhabit halos with median spins twice smaller than the halos hosting the rest of the SRs. Thus, although the formation paths of SRs can be varied, dry mergers and/or halos with small spins dominate., Comment: Accepted for publications in MNRAS (20 pages, 17 main body, 1.5 appendix). Changes include analysis of the orbital angular momentum effect on lambdaR and slightly different ellipticity calculation

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $\Lambda$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*, Comment: 21 pages, 17 figures, version accepted for publication in MNRAS

To provide a quantitative cosmological context to ongoing observational work on the formation histories and location of compact massive galaxies, we locate and study a sample of exceptionally compact systems in the Bolshoi simulation, using the dark matter structural parameters from a real, compact massive galaxy (NGC 1277) as a basis for our working criteria. We find that over 80% of objects in this nominal compact category are substructures of more massive groups or clusters, and that the probability of a given massive substructure being this compact increases significantly with the mass of the host structure; rising to ~30% for the most massive clusters in the simulation. Tracking the main progenitors of this subsample back to z=2, we find them all to be distinct structures with scale radii and densities representative of the population as a whole at this epoch. What does characterise their histories, in addition to mostly becoming substructures, is that they have almost all experienced below-average mass accretion since z=2; a third of them barely retaining, or even losing mass during the intervening 10 Gyr., 9 pages, 9 figures

We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range $14.0 \le \log_{10}(M_{\rm 200c}/\mathrm{M_\odot}) \le 15.4$, with their sophisticated sub-grid physics modelling and high numerical resolution (sub-kpc gravitational softening) making them ideal for this purpose. We find that selecting galaxies by their total mass results in a velocity dispersion that is 5-10 per cent higher than the dark matter particles. However, selecting galaxies by their stellar mass results in an almost unbiased ($, 15 pages, 10 figures, accepted for publication in MNRAS

The diffuse soft X-ray emissivity from galactic winds is computed during the Epoch of Reionization (EoR). We consider two analytic models, a pressure-driven wind and a superbubble model, and a 3D cosmological simulation including gas dynamics from the First Billion Years (FiBY) project. The analytic models are normalized to match the diffuse X-ray emissivity of star-forming galaxies in the nearby Universe. The cosmological simulation uses physically motivated star formation and wind prescriptions, and includes radiative transfer corrections. The models and the simulation all are found to produce sufficient heating of the Intergalactic Medium to be detectable by current and planned radio facilities through 21 cm measurements during the EoR. While the analytic models predict a 21 cm emission signal relative to the Cosmic Microwave Background sets in by $z_{\rm trans} \simeq 8 - 10$, the predicted signal in the FiBY simulation remains in absorption until reionization completes. The 21 cm absorption differential brightness temperature reaches a minimum of $\Delta T \simeq -130$ to $-200$ mK, depending on model. Allowing for additional heat from high mass X-ray binaries pushes the transition to emission to $z_{\rm trans} \simeq 10 - 12$, with shallower absorption signatures having a minimum of $\Delta T \simeq -110$ to $-140$ mK. The 21 cm signal may be a means of distinguishing between the wind models, with the superbubble model favouring earlier reheating. While an early transition to emission may indicate X-ray binaries dominate the reheating, a transition to emission as early as $z_{\rm trans} > 12$ would suggest the presence of additional heat sources., Comment: 16 pages. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. The version of record, MNRAS, 2017, volume 471, issue 3, is available online at: https://academic.oup.com/mnras/article/471/3/3632/4002694/Galactic-wind-X-ray-heating-of-the-intergalactic

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

We present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution towards observable galaxies at $z \sim 6$. We focus on three different galaxies which end up in halos with masses $M_{h} = 2.4 \times10^{10}~h^{-1}\; M_{\odot}$ (Halo-10), $1.6 \times10^{11}~h^{-1}\; M_{\odot}$ (Halo-11) and $0.7 \times10^{12}~h^{-1} M_{\odot}$ (Halo-12) at z=6. Our simulations also probe impacts of different sub-grid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that star formation occurs intermittently due to supernova (SN) feedback at z > 10, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no-feedback or lower SF efficiency models can sustain galactic disk for long periods > 10 Myr. The expulsion of gas at the galactic center also affects the inner dark matter density profile. However, SN feedback does not seem to keep the shallow profile of dark matter for a long period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the star formation rates (SFR) and stellar masses of observed Lyman-$��$ emitters (LAEs) at z = 7-8 fairly well given observational uncertainties. In addition, we investigate the effect of UV background radiation on star formation as an external feedback source, and find that earlier reionization extends the quenching time of star formation due to photo-ionization heating, but does not affect the stellar mass at z=6., 14 pages, 13 figures, submitted to ApJ

We explore the cosmological implications of anisotropic clustering measurements in configuration space of the final galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We implement a new detailed modelling of the effects of non-linearities, galaxy bias and redshift-space distortions that can be used to extract unbiased cosmological information from our measurements for scales $s \gtrsim 20\,h^{-1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS with the latest cosmic microwave background (CMB) observations and Type Ia supernovae samples and found no significant evidence for a deviation from the $\Lambda$CDM cosmological model. In particular, these data sets can constrain the dark energy equation of state parameter to $w_{\rm DE}=-0.996\pm0.042$ when assumed time-independent, the curvature of the Universe to $\Omega_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_{\nu} < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=\Omega_{\rm m}(z)^\gamma$ and obtain $\gamma = 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $\gamma=0.55$. We compress the information of our clustering measurements into constraints on the parameter combinations $D_{\rm V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $f\sigma_8(z)$ at the effective redshifts of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and find good agreement with the predictions for these parameters obtained from the best-fitting $\Lambda$CDM model to the CMB data from the Planck satellite. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS., Comment: 20 pages, 16 figures

We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}, Comment: 21 pages, 15 figures, MNRAS accepted version

We extract cosmological information from the anisotropic power spectrum measurements from the recently completed Baryon Oscillation Spectroscopic Survey (BOSS), extending the concept of clustering wedges to Fourier space. Making use of new FFT-based estimators, we measure the power spectrum clustering wedges of the BOSS sample by filtering out the information of Legendre multipoles l > 4. Our modelling of these measurements is based on novel approaches to describe non-linear evolution, bias, and redshift-space distortions, which we test using synthetic catalogues based on large-volume N-body simulations. We are able to include smaller scales than in previous analyses, resulting in tighter cosmological constraints. Using three overlapping redshift bins, we measure the angular diameter distance, the Hubble parameter, and the cosmic growth rate, and explore the cosmological implications of our full shape clustering measurements in combination with CMB and SN Ia data. Assuming a {\Lambda}CDM cosmology, we constrain the matter density to {\Omega}_m = 0.311 -0.010 +0.009 and the Hubble parameter to H_0 = 67.6 -0.6 +0.7 km s^-1 Mpc^-1, at a confidence level (CL) of 68 per cent. We also allow for non-standard dark energy models and modifications of the growth rate, finding good agreement with the {\Lambda}CDM paradigm. For example, we constrain the equation-of-state parameter to w = -1.019 -0.039 +0.048. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. 2016 to produce the final cosmological constraints from BOSS., Comment: 23 pages, 16 figures in the main text, appendix of 6 pages and 5 figures; revision submitted to MNRAS. The data used in this analysis is publicly available at https://www.sdss3.org/science/boss_publications.php

We study the effect of baryons on the inner dark matter profile of the first galaxies using the First Billion Years simulation between z=16-6 before secular evolution sets in. Using a large statistical sample from two simulations of the same volume and cosmological initial conditions, one with and one without baryons, we are able to directly compare haloes with their baryon-free counterparts, allowing a detailed study of the modifications to the dark matter density profile due to the presence of baryons during the first billion years of galaxy formation. For each of the ~ 5000 haloes in our sample we quantify the impact of the baryons using eta, defined as the ratio of dark matter mass enclosed in 100 pc in the baryonic run to its counterpart without baryons. During this epoch of rapid growth of galaxies, we find that many haloes of these first galaxies show an enhancement of dark matter in the halo centre compared to the baryon-free simulation, while many others show a deficit. We find that the mean value of eta is close to unity, but there is a large dispersion, with a standard deviation of 0.677. The enhancement is cyclical in time and tracks the star formation cycle of the galaxy; as gas falls to the centre and forms stars, the dark matter moves in as well. Supernova feedback then removes the gas, and the dark matter again responds to the changing potential. We study three physical models relating the motion of baryons to that of the dark matter: adiabatic contraction, dynamical friction, and rapid outflows. Abridged, see text for full abstract, Comment: accepted in MNRAS

We investigate the environment in which direct-collapse black holes may form by analysing a cosmological, hydrodynamical simulation that is part of the First Billion Years project. This simulation includes the most relevant physical processes leading to direct collapse of haloes, most importantly, molecular hydrogen depletion by dissociation of $H_2$ and $H^-$ from the evolving Lyman-Werner radiation field. We selected a sample of pristine atomic cooling haloes that have never formed stars in their past, have not been polluted with heavy elements and are cooling predominantly via atomic hydrogen lines. Amongst them we identified six haloes that could potentially harbour massive seed black holes formed via direct collapse (with masses in the range of $10^{4-6} M_{sun}$). These potential hosts of direct-collapse black holes form as satellites and are found within 15 physical kpc of proto-galaxies, with stellar masses in the range $10^{5-7} M_{sun}$ and maximal star formation rates of 0.1 Msun/yr over the past 5 Myr, and are exposed to the highest flux of Lyman-Werner radiation emitted from the neighbouring galaxies. It is the proximity to these proto-galaxies that differentiates these haloes from rest of the sample., 10 pages, 8 figures, accepted MNRAS

We use cosmological hydrodynamical simulations to investigate how the inclusion of physical processes relevant to galaxy formation (star formation, metal-line cooling, stellar winds, supernovae and feedback from Active Galactic Nuclei, AGN) change the properties of haloes, over four orders of magnitude in mass. We find that gas expulsion and the associated dark matter (DM) expansion induced by supernova-driven winds are important for haloes with masses M200 < 10^13 Msun, lowering their masses by up to 20% relative to a DM-only model. AGN feedback, which is required to prevent overcooling, has a significant impact on halo masses all the way up to cluster scales (M200 ~ 10^15 Msun). Baryonic physics changes the total mass profiles of haloes out to several times the virial radius, a modification that cannot be captured by a change in the halo concentration. The decrease in the total halo mass causes a decrease in the halo mass function of about 20%. This effect can have important consequences for abundance matching technique as well as for most semi-analytic models of galaxy formation. We provide analytic fitting formulae, derived from simulations that reproduce the observed baryon fractions, to correct halo masses and mass functions from DM-only simulations. The effect of baryonic physics (AGN feedback in particular) on cluster number counts is about as large as changing the cosmology from WMAP7 to Planck, even when a moderately high mass limit of M500 ~ 10^14 Msun is adopted. Thus, for precision cosmology the effects of baryons must be accounted for., 21 pages, 15 figures. Accepted for publication in MNRAS after minor revisions

We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on metal-line cooling and feedback from star formation and active galactic nuclei (AGN). We find that if the sub-grid feedback from star formation is implemented kinetically, the feedback is only efficient if the initial wind velocity exceeds a critical value. This critical velocity increases with galaxy mass and also if metal-line cooling is included. This suggests that radiative losses quench the winds if their initial velocity is too low. If the feedback is efficient, then the star formation rate is inversely proportional to the amount of energy injected per unit stellar mass formed (which is proportional to the initial mass loading for a fixed wind velocity). This can be understood if the star formation is self-regulating, i.e. if the star formation rate (and thus the gas fraction) increase until the outflow rate balances the inflow rate. Feedback from AGN is efficient at high masses, while increasing the initial wind velocity with gas pressure or halo mass allows one to generate galaxy-wide outflows at all masses. Matching the observed galaxy mass function requires efficient feedback. In particular, the predicted faint-end slope is too steep unless we resort to highly mass loaded winds for low-mass objects. Such efficient feedback from low-mass galaxies (M_* << 10^10 Msun) also reduces the discrepancy with the observed specific star formation rates, which are higher than predicted unless the feedback transitions from highly efficient to inefficient just below the observed stellar mass range., 25 pages, 15 figures, accepted by MNRAS

We investigate the evolution of supermassive black holes in the `Evolution and Assembly of GaLaxies and their Environments' (EAGLE) cosmological hydrodynamic simulations. The largest of the EAGLE volumes covers a $(100 \,\rm cMpc)^3$ and includes state-of-the-art physical models for star formation and black hole growth that depend only on local gas properties. We focus on the black hole mass function, Eddington ratio distribution and the implied duty cycle of nuclear activity. The simulation is broadly consistent with observational constraints on these quantities. In order to make a more direct comparison with observational data, we calculate the soft and hard X-ray luminosity functions of the active galactic nuclei (AGN). Between redshifts $0$ and $1$, the simulation is in agreement with data. At higher redshifts, the simulation tends to underpredict the luminosities of the brightest observed AGN. This may be due to the limited volume of the simulation, or a fundamental deficiency of the underlying model. It seems unlikely that additional unresolved variability can account for this difference. The simulation shows a similar `downsizing' of the AGN population as seen in observational surveys., Comment: MNRAS accepted for publication. 19 pages, 3 tables and 14 figures

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 �� to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement., 12 pages, 8 figures, v.2: ApJ, accepted. Model grids are available from http://www.astro.uu.se/~ez/lycan/lycan.html

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods., Comment: 21 pages, 13 figures, 4 tables - submitted to MNRAS