Search

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

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

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

Massive high-redshift quiescent compact galaxies (nicknamed red nuggets) have been traditionally connected to present-day elliptical galaxies, often overlooking the relationships that they may have with other galaxy types. We use large bulge-disk decomposition catalogues based on the Sloan Digital Sky Survey (SDSS) to check the hypothesis that red nuggets have survived as compact cores embedded inside the haloes or disks of present-day massive galaxies. In this study, we designate a "compact core" as the bulge component that satisfies a prescribed compactness criterion. Photometric and dynamic mass-size and mass-density relations are used to show that, in the inner regions of galaxies at z ~ 0.1, there are "abundant" compact cores matching the peculiar properties of the red nuggets, an abundance comparable to that of red nuggets at z ~ 1.5. Furthermore, the morphology distribution of the present-day galaxies hosting compact cores is used to demonstrate that, in addition to the standard channel connecting red nuggets with elliptical galaxies, a comparable fraction of red nuggets might have ended up embedded in disks. This result generalises the inside-out formation scenario; present-day massive galaxies can begin as dense spheroidal cores (red nuggets), around which either a spheroidal halo or a disk are formed later., Comment: 17 pages, 14 figures. Revised to match published MNRAS version; minor editorial changes throughout and updated references, with no change to results. 2016MNRAS.457.1916D

Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionizing photons needed to reionize the Universe, due to their high escape fractions of ionizing photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form, by computing escape fractions in more than 75000 haloes between redshifts 27 and 6 that were extracted from the First Billion Years project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the gas column density in a radius of 10 pc around the stellar populations, causing a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 1e8 solar mass results in low column densities that can be penetrated by radiation from young stars (age < 5 Myr). In haloes with virial mass above 1e8 solar mass supernova feedback is important, but only 30% of the haloes in this mass range have an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. This makes it very hard to predict the escape fraction on the basis of halo properties and results in a highly anisotropic escape fraction. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionization and is something current models of cosmic reionization should strive to take into account., 21 pages, 17 figures, accepted for publication in MNRAS

One of the most promising explanations for the origin of the billion solar mass black holes (BHs) inferred to power quasars at redshifts z > 6 is that supermassive stars (SMSs) with masses > 10,000 solar masses collapse to form the seed BHs from which they grow. Here we review recent theoretical advances which provide support for this scenario. Firstly, given sufficiently high accretion rates of gas into the cores of primordial protogalaxies, it appears that neither the high energy radiation emitted from the stellar surface nor the limited lifetime of SMSs can prevent their growth to masses of up to > 100,000 solar masses. Secondly, recent cosmological simulations suggest that the high fluxes of molecule-dissociating radiation which may be required in order to achieve such high accretion rates may be more common in the early universe than previously thought. We conclude that the majority of supermassive BHs may originate from SMSs at high redshifts., 4 pages, 3 figures; proceedings of "First Stars IV", held in Kyoto, Japan, May 21-25, 2012

We investigate the column density distribution function of neutral hydrogen at redshift z = 3 using a cosmological simulation of galaxy formation from the OverWhelmingly Large Simulations (OWLS) project. The base simulation includes gravity, hydrodynamics, star formation, supernovae feedback, stellar winds, chemodynamics, and element-by-element cooling in the presence of a uniform UV background. Self-shielding and formation of molecular hydrogen are treated in post-processing, without introducing any free parameters, using an accurate reverse ray-tracing algorithm and an empirical relation between gas pressure and molecular mass fraction. The simulation reproduces the observed z = 3 abundance of Ly-A forest, Lyman Limit and Damped Ly-A HI absorption systems probed by quasar sight lines over ten orders of magnitude in column density. Self-shielding flattens the column density distribution for NHI > 10^18 cm-2, while the conversion to fully neutral gas and conversion of HI to H2 steepen it around column densities of NHI = 10^20.3 cm-2 and NHI = 10^21.5 cm-2, respectively., updated to match ApJL accepted version. now includes line list for fiducial model

Feedback from star formation is thought to play a key role in the formation and evolution of galaxies, but its implementation in cosmological simulations is currently hampered by a lack of numerical resolution. We present and test a sub-grid recipe to model feedback from massive stars in cosmological smoothed particle hydrodynamics simulations. The energy is distributed in kinetic form among the gas particles surrounding recently formed stars. The impact of the feedback is studied using a suite of high-resolution simulations of isolated disc galaxies embedded in dark halos with total mass 10^{10} and 10^{12} Msol/h. We focus in particular on the effect of pressure forces on wind particles within the disc, which we turn off temporarily in some of our runs to mimic a recipe that has been widely used in the literature. We find that this popular recipe gives dramatically different results because (ram) pressure forces on expanding superbubbles determine both the structure of the disc and the development of large-scale outflows. Pressure forces exerted by expanding superbubbles puff up the disc, giving the dwarf galaxy an irregular morphology and creating a galactic fountain in the massive galaxy. Hydrodynamic drag within the disc results in a strong increase of the effective mass loading of the wind for the dwarf galaxy, but quenches much of the outflow in the case of the high-mass galaxy., Comment: 15 pages, 11 figures, accepted for publication in MNRAS. Added 4 new figures and made minor textual changes. Simulation videos available at http://www.strw.leidenuniv.nl/DS08/

We use a combination of three large N-body simulations to investigate the dependence of dark matter halo concentrations on halo mass and redshift in the WMAP year 5 cosmology. The median relation between concentration and mass is adequately described by a power-law for halo masses in the range 10^11 - 10^15 Msol/h and redshifts z < 2, regardless of whether the halo density profiles are fit using NFW or Einasto profiles. Compared with recent analyses of the Millennium Simulation, which uses a value of sigma_8 that is higher than allowed by WMAP5, z = 0 halo concentrations are reduced by factors ranging from 23 per cent at 10^11 Msol/h to 16 per cent at 10^14 Msol/h. The predicted concentrations are much lower than inferred from X-ray observations of groups and clusters., 6 pages, 4 figures, Accepted by MNRAS letters. Version 4: Typo fixed