Sergio Contreras

Bio: Sergio Contreras is an academic researcher from Donostia International Physics Center. The author has contributed to research in topics: Galaxy & Physics. The author has an hindex of 15, co-authored 33 publications receiving 763 citations. Previous affiliations of Sergio Contreras include Pontifical Catholic University of Chile.

Quantifying the impact of mergers on the angular momentum of simulated galaxies

Abstract: We use EAGLE to quantify the effect galaxy mergers have on the stellar specific angular momentum of galaxies, jstars. 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 ≥0.3 (0.1–0.3). We correlate the positions of galaxies in the jstars–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 jstars profiles of galaxies that went through dry mergers are deficient by ≈0.3 dex at r ≲ 10 r50 (with r50 being the half-stellar mass radius), compared to galaxies that went through wet mergers. Studying the merger remnants reveals that dry mergers reduce jstars by ≈30 per cent, while wet mergers increase it by ≈10 per cent, 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, jstars decreases by a greater magnitude, while in wet mergers corotation and high orbital angular momentum efficiently spun-up galaxies. We predict what would be the observational signatures in the jstars profiles driven by dry mergers: (i) shallow radial profiles and (ii) profiles that rise beyond ≈10 r50, both of which are significantly different from spiral galaxies.

Cosmic web alignments with the shape, angular momentum and peculiar velocities of dark matter haloes

Abstract: We study the alignment of dark matter haloes with the cosmic web characterized by the tidal and velocity shear fields. We focus on the alignment of their shape, angular momentum and peculiar velocities. We use a cosmological N-body simulation that allows to study dark matter halos spanning almost five orders of magnitude in mass ($10^{9}$-$10^{14}$) $h^{-1}$$M_{\odot}$ and spatial scales of $(0.5$-$1.0)$ $h^{-1}$ Mpc to define the cosmic web. We find that the halo shape presents the strongest alignment along the smallest tidal eigenvector, e.g. along filaments and walls, with a signal that gets stronger as the halo mass increases. In the case of the velocity shear field only massive halos $>10^{12}$ $h^{-1}$$M_{\odot}$ tend to have their shapes aligned along the largest tidal eigenvector; that is, perpendicular to filaments and walls. For the angular momentum we find alignment signals only for halos more massive than $10^{12}$ $h^{-1}$$M_{\odot}$ both in the tidal and velocity shear webs where the preferences are for it to be parallel to the middle eigenvector; perpendicular to filaments and parallel to walls. Finally, the peculiar velocities show a strong alignment along the smallest tidal eigenvector for all halo masses; halos move along filaments and walls. In the velocity shear the same alignment is present but weaker and only for haloes less massive than $10^{12}$ $h^{-1}$$M_{\odot}$. Our results clearly show that the two different algorithms we used to define the cosmic web describe different physical aspects of non-linear collapse and should be used in a complementary way to understand the effect of the cosmic web on galaxy evolution.

The BACCO Simulation Project: Exploiting the full power of large-scale structure for cosmology

Abstract: We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear mass power spectrum over the redshift range $0

The impact of assembly bias on the halo occupation in hydrodynamical simulations.

Abstract: We investigate the variations in galaxy occupancy of the dark matter haloes with the large-scale environment and halo formation time, using two state-of-the-art hydrodynamical cosmological simulations, EAGLE and Illustris. For both simulations, we use three galaxy samples with a fixed number density ranked by stellar mass. For these samples, we find that low-mass haloes in the most dense environments are more likely to host a central galaxy than those in the least dense environments. When splitting the halo population by formation time, these relations are stronger. Hence, at a fixed low halo mass, early-formed haloes are more likely to host a central galaxy than late-formed haloes since they have had more time to assemble. The satellite occupation shows a reverse trend where early-formed haloes host fewer satellites due to having more time to merge with the central galaxy. We also analyse the stellar mass–halo mass relation for central galaxies in terms of the large-scale environment and formation time of the haloes. We find that low-mass haloes in the most dense environment host relatively more massive central galaxies. This trend is also found when splitting the halo population by age, with early-formed haloes hosting more massive galaxies. Our results are in agreement with previous findings from semi-analytical models, providing robust predictions for the occupancy variation signature in the halo occupation distribution of galaxy formation models.

The host dark matter haloes of [O II] emitters at 0.5 < z < 1.5

Abstract: Emission line galaxies (ELGs) are used in several ongoing and upcoming surveys (SDSSIV/eBOSS, DESI) as tracers of the dark matter distribution. Using a new galaxy formation model, we explore the characteristics of [O II] emitters, which dominate optical ELG selections at z 1. Model [O II] emitters at 0.5 < z < 1.5 are selected to mimic the DEEP2, VVDS, eBOSS and DESI surveys. The luminosity functions of model [O II] emitters are in reasonable agreement with observations. The selected [O II] emitters are hosted by haloes with Mhalo ≥ 1010.3h−1M, with ∼90 per cent of them being central star-forming galaxies. The predicted mean halo occupation distributions of [O II] emitters have a shape typical of that inferred for star-forming galaxies, with the contribution from central galaxies, N[O II] cen, being far from the canonical step function. The N[O II] cen can be described as the sum of an asymmetric Gaussian for discs and a step function for spheroids, which plateau below unity. The model [O II] emitters have a clustering bias close to unity, which is below the expectations for eBOSS and DESI ELGs. At z ∼ 1, a comparison with observed g-band-selected galaxy, which is expected to be dominated by [O II] emitters, indicates that our model produces too few [O II] emitters that are satellite galaxies. This suggests the need to revise our modelling of hot gas stripping in satellite galaxies.

PHIBSS: Unified Scaling Relations of Gas Depletion Time and Molecular Gas Fractions

Abstract: This paper provides an update of our previous scaling relations between galaxy-integrated molecular gas masses, stellar masses, and star formation rates (SFRs), in the framework of the star formation main sequence (MS), with the main goal of testing for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ∼1 mm dust photometry, in a large sample of 1444 star-forming galaxies between z=0 and 4. The sample covers the stellar mass range log(M * /M e)=9.0-11.8, and SFRs relative to that on the MS, δMS=SFR/SFR (MS), from 10 −1.3 to 10 2.2. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero-point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t depl , defined as the ratio of molecular gas mass to SFR, scales as (1+z) −0.6 ×(δMS) −0.44 and is only weakly dependent on stellar mass. The ratio of molecular to stellar mass μ gas depends on (* d +´´-) () () z M 1 M S 2.5 0.52 0.36 , which tracks the evolution of the specific SFR. The redshift dependence of μ gas requires a curvature term, as may the mass dependences of t depl and μ gas. We find no or only weak correlations of t depl and μ gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high z.

First results from the TNG50 simulation: the evolution of stellar and gaseous discs across cosmic time

Abstract: We present a new cosmological, magnetohydrodynamical simulation for galaxy formation: TNG50, the third and final instalment of the Illustris TNG project. TNG50 evolves 2 x 2160(3) dark matter particles and gas cells in a volume 50 comoving Mpc across. It hence reaches a numerical resolution typical of zoom-in simulations, with a baryonic element mass of 8.5 x 10(4) M-circle dot and an average cell size of 70-140 pc in the star-forming regions of galaxies. Simultaneously, TNG50 samples similar to 700 (6500) galaxies with stellar masses above 10(10) (10(8)) M-circle dot at z = 1. Here we investigate the structural and kinematical evolution of star-forming galaxies across cosmic time (0 less than or similar to z less than or similar to 6). We quantify their sizes, disc heights, 3D shapes, and degree of rotational versus dispersion-supported motions as traced by rest-frame V-band light (i.e. roughly stellar mass) and by H alpha light (i.e. star-forming and dense gas). The unprecedented resolution of TNG50 enables us to model galaxies with sub-kpc half-light radii and with less than or similar to 300-pc disc heights. Coupled with the large-volume statistics, we characterize a diverse, redshift- and mass-dependent structural and kinematical morphological mix of galaxies all the way to early epochs. Our model predicts that for star-forming galaxies the fraction of disc-like morphologies, based on 3D stellar shapes, increases with both cosmic time and galaxy stellar mass. Gas kinematics reveal that the vast majority of 10(9-11.5) M-circle dot star-forming galaxies are rotationally supported discs for most cosmic epochs (V-rot/sigma > 2-3, z less than or similar to 5), being dynamically hotter at earlier epochs (z greater than or similar to 1.5). Despite large velocity dispersion at high redshift, cold and dense gas in galaxies predominantly arranges in disky or elongated shapes at all times and masses; these gaseous components exhibit rotationally dominated motions far exceeding the collisionless stellar bodies.

A unified multiwavelength model of galaxy formation

Abstract: We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850μm. This was not previously achieved by a single physical model within the ΛCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z = 0 to z = 6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z = 0, as well as far-IR number counts, and far-UV luminosity functions at z ∼ 3 − 6. Discrepancies are however found in galaxy sizes and metallicities at low luminosities, and in the abundance of low mass galaxies at high-z, suggesting the need for a more sophisticated model of supernova feedback.

The SDSS-IV extended baryon oscillation spectroscopic survey: overview and early data

Abstract: In a six-year program started in 2014 July, the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. These observations will be conducted simultaneously with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. In particular, eBOSS will measure with percent-level precision the distance-redshift relation with baryon acoustic oscillations (BAO) in the clustering of matter. eBOSS will use four different tracers of the underlying matter density field to vastly expand the volume covered by BOSS and map the large-scale-structures over the relatively unconstrained redshift range 0.6 0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of d_A(z) to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z = 0.87. A sample of more than 500,000 spectroscopically confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9 2.1; these new data will enhance the precision of dA(z) and H(z) at z > 2.1 by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

Bimodal gas accretion in the HORIZON-MareNostrum galaxy formation simulation

Abstract: The physics of diffuse gas accretion and the properties of the cold and hot modes of accretion on to proto-galaxies between z= 2 and 5.4 is investigated using the large cosmological simulation performed with the ramses code on the MareNostrum supercomputing facility. Galactic winds, chemical enrichment, ultraviolet background heating and radiative cooling are taken into account in this very high resolution simulation. Using accretion-weighted temperature histograms, we have performed two different measurements of the thermal state of the gas accreted towards the central galaxy. The first measurement, performed using accretion-weighted histograms on a spherical surface of radius 0.2Rvir centred on the densest gas structure near the halo centre of mass, is a good indicator of the presence of an accretion shock in the vicinity of the galactic disc. We define the hot shock mass, Mshock, as the typical halo mass separating cold dominated from hot dominated accretion in the vicinity of the galaxy. The second measurement is performed by radially averaging histograms between 0.2Rvir and Rvir, in order to detect radially extended structures such as gas filaments: this is a good proxy for detecting cold streams feeding the central galaxy. We define Mstream as the transition mass separating cold dominated from hot dominated accretion in the outer halo, marking the disappearance of these cold streams. We find a hot shock transition mass of Mshock= 1011.6 M⊙ (dark matter), with no significant evolution with redshift. Conversely, we find that Mstream increases sharply with z. Our measurements are in agreement with the analytical predictions of Birnboim & Dekel and Dekel & Birnboim, if we correct their model by assuming low metallicity(≤10−3 Z⊙)for the filaments, correspondingly to our measurements. Metal enrichment of the intergalactic medium is therefore a key ingredient in determining the transition mass from cold to hot dominated diffuse gas accretion. We find that the diffuse cold gas supply at the inner halo stops at z= 2 for objects with stellar masses of about 1011.1 M⊙, which is close to the quenching mass determined observationally by Bundy et al. However, its evolution with z is not well constrained, making it difficult to rule out or confirm the need for an additional feedback process such as active galactic nuclei.