Showing 1–50 of 97 results for author: Contreras, S

SHAMe-SF: Predicting the clustering of star-forming galaxies with an enhanced abundance matching model

Authors: Sara Ortega-Martinez, Sergio Contreras, Raul Angulo

Abstract: With the advent of several galaxy surveys targeting star-forming galaxies, it is important to have models capable of interpreting their spatial distribution in terms of astrophysical and cosmological parameters. To address this need, we introduce SHAMe-SF, an extension of the subhalo abundance matching (SHAM) technique designed specifically for analyzing the redshift-space clustering of star-formi… ▽ More With the advent of several galaxy surveys targeting star-forming galaxies, it is important to have models capable of interpreting their spatial distribution in terms of astrophysical and cosmological parameters. To address this need, we introduce SHAMe-SF, an extension of the subhalo abundance matching (SHAM) technique designed specifically for analyzing the redshift-space clustering of star-forming galaxies. Our model directly links a galaxy's star formation rate to the properties of its host dark-matter halo, with further modulations based on effective models of feedback and gas stripping. To quantify the accuracy of our model, we show that it simultaneously reproduces key clustering statistics such as the projected correlation function, monopole, and quadrupole of star-forming galaxy samples at various redshifts and number densities. Notably, these tests were conducted over a wide range of scales $[0.6, 30]\hMpc$, using samples from both the TNG300 magneto-hydrodynamic simulation and from a semi-analytical model. SHAMe-SF can also reproduce the clustering of simulated galaxies that fall within the colour selection criteria employed by DESI for emission line galaxies. Our model exhibits several potential applications, including the generation of covariance matrices, exploration of galaxy formation processes, and even placing constraints on the cosmological parameters of the Universe. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 18 pages, 10 figures

CO emission survey of asymptotic giant branch stars with ultraviolet excesses

Authors: Jaime Alonso-Hernández, Carmen Sánchez Contreras, Raghvendra Sahai

Abstract: Context. The transition from the spherically symmetric envelopes around asymptotic giant branch (AGB) stars to the asymmetric morphologies observed in planetary nebulae is still not well understood, and the shaping mechanisms are a subject of debate. Even though binarity is widely accepted as a promising option. Recently, the presence of ultraviolet excesses in AGB stars has been suggested as a po… ▽ More Context. The transition from the spherically symmetric envelopes around asymptotic giant branch (AGB) stars to the asymmetric morphologies observed in planetary nebulae is still not well understood, and the shaping mechanisms are a subject of debate. Even though binarity is widely accepted as a promising option. Recently, the presence of ultraviolet excesses in AGB stars has been suggested as a potential indicator of binarity. Aims. Our main goals are to characterise the properties of the circumstellar envelopes (CSEs) around candidate AGB binary stars, specifically those selected based on their UV excess emission, and to compare these properties with those derived from previous CO-based studies of AGB stars. Methods. We observed the 12CO ($J$=1-0) and 12CO ($J$=2-1) millimetre-wavelength emission in a sample of 29 AGB binary candidates with the IRAM-30 m antenna. We explored different trends between the envelope parameters deduced and compared them with those previously derived from larger samples of AGB stars found in the literature. Results. We derived the average excitation temperature and column density of the CO-emitting layers, which we used to estimate self-consistently the average mass-loss rate and the CO photodissociation radius of our targets. We find a correlation between CO intensity and IRAS 60$μ$m fluxes, revealing a CO-to-IRAS 60$μ$m ratio lower than for AGB stars and closer to that found for pre-planetary nebulae (pPNe). Conclusions. For the first time we have studied the mass-loss properties of UV-excess AGB binary candidates and estimated their main CSE parameters. The different relationships between 12CO and IRAS 60μm, with NUV and FUV are consistent with an intrinsic origin of NUV emission, but potential dominance of an extrinsic process (e.g. presence of a binary companion) in FUV emission. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 30 pages, 13 figures and 6 tables (A&A, in press)

Zooming on the emerging ionized regions of pPNe with ALMA

Authors: C. Sánchez Contreras, D. Tafoya, J. P. Fonfría, J. Alcolea, A. Castro-Carrizo, V. Bujarrabal

Abstract: We report on recent results from our successful and pioneering observational program with ALMA to study emerging ultracom pact HII regions of pre-Planetary Nebulae (pPNe) using mm-wavelength recombination lines (mRRLs) as new optimal tracers. We focus on our study of two poster-child pPNe, namely, M2-9 and CRL618. We reveal the structure and kinematics of the en igmatic inner nebular regions of th… ▽ More We report on recent results from our successful and pioneering observational program with ALMA to study emerging ultracom pact HII regions of pre-Planetary Nebulae (pPNe) using mm-wavelength recombination lines (mRRLs) as new optimal tracers. We focus on our study of two poster-child pPNe, namely, M2-9 and CRL618. We reveal the structure and kinematics of the en igmatic inner nebular regions of these objects with an unprecedented angular resolution down to 20-30mas (~15-30AU). For both targets, the ionized central regions are elongated along the main symmetry axis of the large-scale nebulae, consiste nt with bipolar winds, and show notable axial velocity gradients with expansion velocities of up to ~100km/s. The H30a pr ofiles exhibit time variability, reflecting changes in the physical properties and kinematics on scales of a few years. O ur ongoing analysis employs 3D, non-LTE radiative transfer modeling, providing a detailed description of the innermost la yers of these well known pPNe with exceptional clarity. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 6 pages, 3 figures. Proceedings IAU Symposium no. 384

The Red Rectangle: a thin disk with big grains

Authors: Javier Alcolea, Valentin Bujarrabal, Arancha Castro-Carrizo, Jacques Kluska, Carmen Sanchez Contreras, Hans van Winckel

Abstract: The Red Rectangle is a nebula surrounding the post-AGB star HD 44179. It is the prototype of a particular class of nebulae associated with post-AGB binaries characterised by the presence of stable circumbinary disks in (quasi-)Keplerian rotation. Here we present the results of new high-resolution (20-50 mas) ALMA observations of continuum and line emissions at 0.9 mm. The continuum maps are analys… ▽ More The Red Rectangle is a nebula surrounding the post-AGB star HD 44179. It is the prototype of a particular class of nebulae associated with post-AGB binaries characterised by the presence of stable circumbinary disks in (quasi-)Keplerian rotation. Here we present the results of new high-resolution (20-50 mas) ALMA observations of continuum and line emissions at 0.9 mm. The continuum maps are analysed through a simple model of dust emission, which can reproduce the observational data. We find that most dust emission in the Red Rectangle is concentrated in the central regions of the rotating disk and that the settlement of dust grains onto the equatorial plane is very significant, particularly in comparison with the much larger scale height displayed by the gas distribution. The diameter of the dust-emitting region is about 250 au, with a total width of about 50 au. This region coincides with the warm PDR where certain molecules (like HCN), CI, and CII are presumably formed, as well as probably PAHs. From the spectral index, we confirm the presence in the disk of large grains, with a typical radius of about 0.150 mm, which supports the long-lived hypothesis for this structure. We also confirm the existence of a compact ionised wind at the centre of the nebula, probably emerging from the accretion disk around the companion, for which we derive an extent of about 10 au and a total flux of 8 mJy. We also briefly present the results on molecular lines of 12CO, 13CO, and other less abundant species. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 6 pages, 4 figures, 1 table, Proceedings IAU Symposium No. 384 "Planetary Nebulae: a Universal Toolbox in the Era of Precision Astrophysics"

The very compact dust disk in the Red Rectangle

Authors: V. Bujarrabal, J. Alcolea, A. Castro-Carrizo, J. Kluska, C. Sanchez Contreras, H. Van Winckel

Abstract: We aim to study the dust distribution in the central regions of the Keplerian disk of the Red Rectangle, the prototype of binary post-AGB stars with rotating circumbinary disks, and to compare it with the distribution of relevant molecular gas tracers We present new high-resolution (20 milliarcseconds, mas) ALMA observations of continuum and line emissions at 0.9 mm. The maps have been analyzed… ▽ More We aim to study the dust distribution in the central regions of the Keplerian disk of the Red Rectangle, the prototype of binary post-AGB stars with rotating circumbinary disks, and to compare it with the distribution of relevant molecular gas tracers We present new high-resolution (20 milliarcseconds, mas) ALMA observations of continuum and line emissions at 0.9 mm. The maps have been analyzed by means of a simple model of dust and free-free emissionn that is able to reproduce the continuum data. Resuts: i) We find that most of the dust emission in the Red Rectangle is concentrated in the inner disk regions, with a typical size of 250 AU in diameter and 50 AU in width. ii) The settlement of dust grains onto inner equatorial regions is remarkable when compared with the relatively widespread gas distribution. iii) This region is basically coincident with the warm PDR (photo-dominated region) where CI, CII, and certain molecules such as HCN are presumably formed, as well as probably PAHs (polycyclic aromatic hydrocarbons, whose emission is very strong in this source). iv) We confirm the large size of the grains, with a typical radius ~ 150 mu The opacity of dust at 0.9 mm is deduced to be relatively large, ~0.5. v) We also confirm the existence of a very compact HII region in the center, for which we measure an extent of 10 - 15 mas (~ 10 AU) and a total flux of 7 - 8 mJy at 0.9 mm. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 5 pages, 4 figures, accepted by Astronomy and Astrophysics Letters

On the origin of the evolution of the halo occupation distribution

Authors: Sergio Contreras, Idit Zehavi

Abstract: We use the TNG300 magneto-hydrodynamic simulation and mock catalogues built using subhalo abundance matching (SHAM) to study the origin of the redshift evolution of the halo occupation distribution (HOD). We analyse stellar-mass selected galaxy samples with fixed number densities, spanning the redshift range $0 \le z \le 3$. We measure their halo occupation functions and fit the HOD parameters to… ▽ More We use the TNG300 magneto-hydrodynamic simulation and mock catalogues built using subhalo abundance matching (SHAM) to study the origin of the redshift evolution of the halo occupation distribution (HOD). We analyse stellar-mass selected galaxy samples with fixed number densities, spanning the redshift range $0 \le z \le 3$. We measure their halo occupation functions and fit the HOD parameters to study their evolution over cosmic time. The TNG300 galaxy population strongly depends on the baryonic physics implemented in the simulation. In contrast, the galaxy population predicted by a basic SHAM model without scatter is a direct result of the cosmology of the dark matter simulation. We find that the HOD evolution is similar for both models and is consistent with a previous study of the HOD evolution in semi-analytical models. Specifically, this is the case for the ratio between the characteristic halo masses for hosting central and satellite galaxies. The only HOD parameter whose evolution varies across models is $σ_{\rm logM}$, which contains information about the stellar mass-halo mass relation of the galaxies and does not strongly impact galaxy clustering. We also demonstrate that the dependence on the specific values of the cosmological parameters is small. We conclude that the cosmology of the galaxy sample, i.e. the cosmological hierarchical growth of structure, and not the baryonic physics prescriptions, governs the evolution of the HOD for stellar mass-selected samples. These results have important implications for populating simulated lightcones with galaxies and can facilitate the interpretation of clustering data at different redshifts. △ Less

Submitted 9 August, 2023; v1 submitted 31 May, 2023; originally announced May 2023.

Comments: 12 pages, 8 figures. Accepted in MNRAS

Consistent clustering and lensing of SDSS-III BOSS galaxies with an extended abundance matching formalism

Authors: Sergio Contreras, Jonás Chaves-Montero, Raul E. Angulo

Abstract: Several analyses have shown that LCDM-based models cannot jointly describe the clustering (GC) and galaxy-galaxy lensing (GGL) of galaxies in the SDSS-III BOSS survey, which is commonly known as the 'lensing-is-low problem'. In this work, we show that an extension of Subhalo Abundance Matching, dubbed SHAMe, successfully solves this problem. First, we show that this model accurately reproduces the… ▽ More Several analyses have shown that LCDM-based models cannot jointly describe the clustering (GC) and galaxy-galaxy lensing (GGL) of galaxies in the SDSS-III BOSS survey, which is commonly known as the 'lensing-is-low problem'. In this work, we show that an extension of Subhalo Abundance Matching, dubbed SHAMe, successfully solves this problem. First, we show that this model accurately reproduces the GC and GGL of a mock galaxy sample in the TNG300 hydrodynamic simulation with analogous properties to BOSS galaxies. Then, we switch our attention to observed BOSS galaxies at z=0.31-0.43, and we attempt to reproduce their GC and GGL by evaluating SHAMe on two different simulations: one adopting best-fitting cosmological parameters from Planck and the other from weak gravitational lensing surveys (Low S8), where the amplitude of matter fluctuations is lower for the latter. We find excellent agreement between SHAMe predictions and observations for both cosmologies, indicating that the lensing-is-low problem originates from approximations in previous theoretical descriptions of the data. The main difference between SHAMe results in these cosmologies is the level of galaxy assembly bias, which is approximately 20 and 10% for Planck and Low S8, respectively. These results highlight the dangers of employing oversimplified models to analyse current large-scale structure datasets, and the need for realistic yet flexible descriptions of the galaxy-halo connection. △ Less

Submitted 9 August, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures. Accepted by MNRAS

DES Y3 cosmic shear down to small scales: constraints on cosmology and baryons

Authors: Giovanni Aricò, Raul E. Angulo, Matteo Zennaro, Sergio Contreras, Angela Chen, Carlos Hernández-Monteagudo

Abstract: We present the first analysis of cosmic shear measured in DES Y3 that employs the entire range of angular scales in the data. To achieve this, we build upon recent advances in the theoretical modelling of weak lensing provided by a combination of $N$-body simulations, physical models of baryonic processes, and neural networks. Specifically, we use BACCOemu to model the linear and nonlinear matter… ▽ More We present the first analysis of cosmic shear measured in DES Y3 that employs the entire range of angular scales in the data. To achieve this, we build upon recent advances in the theoretical modelling of weak lensing provided by a combination of $N$-body simulations, physical models of baryonic processes, and neural networks. Specifically, we use BACCOemu to model the linear and nonlinear matter power spectrum including baryonic physics, allowing us to robustly exploit scales smaller than those used by the DES Collaboration. We show that the additional data produce cosmological parameters that are tighter but consistent with those obtained from larger scales, while also constraining the distribution of baryons. In particular, we measure the mass scale at which haloes have lost half of their gas, $\log\,M_{\rm c}=14.38^{+0.60}_{-0.56}\log(h^{-1}{\rm M_{ \odot}})$, and a parameter that quantifies the weighted amplitudes of the present-day matter inhomogeneities, $S_8=0.799^{+0.023}_{-0.015}$. Our constraint on $S_8$ is statistically compatible with that inferred from the Planck satellite's data at the $0.9σ$ level. We find instead a $1.4σ$ shift in comparison to that from the official DES Y3 cosmic shear, because of different choices in the modelling of intrinsic alignment, non-linearities, baryons, and lensing shear ratios. We conclude that small scales in cosmic shear data contain valuable astrophysical and cosmological information and thus should be included in standard analyses. △ Less

Submitted 7 August, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: Accepted for publication on A&A

Journal ref: A&A 678, A109 (2023)

The messy death of a multiple star system and the resulting planetary nebula as observed by JWST

Authors: Orsola De Marco, Muhammad Akashi, Stavros Akras, Javier Alcolea, Isabel Aleman, Philippe Amram, Bruce Balick, Elvire De Beck, Eric G. Blackman, Henri M. J. Boffin, Panos Boumis, Jesse Bublitz, Beatrice Bucciarelli, Valentin Bujarrabal, Jan Cami, Nicholas Chornay, You-Hua Chu, Romano L. M. Corradi, Adam Frank, Guillermo Garcia-Segura, D. A. Garcia-Hernandez, Jorge Garcia-Rojas, Veronica Gomez-Llanos, Denise R. Goncalves, Martin A. Guerrero , et al. (44 additional authors not shown)

Abstract: Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures,… ▽ More Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures, likely shaped by a low-mass companion orbiting the central star at 40-60 AU. The images also reveal a mid-IR excess at the central star interpreted as a dusty disk, indicative of an interaction with another, closer companion. Including the previously known, A-type visual companion, the progenitor of the NGC3132 PN must have been at least a stellar quartet. The JWST images allow us to generate a model of the illumination, ionisation and hydrodynamics of the molecular halo, demonstrating the power of JWST to investigate complex stellar outflows. Further, new measurements of the A-type visual companion allow us to derive the value for the mass of the progenitor of a central star to date with excellent precision: 2.86+/-0.06 Mo. These results serve as path finders for future JWST observations of PNe providing unique insight into fundamental astrophysical processes including colliding winds, and binary star interactions, with implications for supernovae and gravitational wave systems. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 32 pages, 5 figures for the main article. 12 pages 8 figures for the supplementary material

Journal ref: Nature Astronomy, 2022, Vol. 6, p. 1421

The binary and the disk: the beauty is found within NGC3132 with JWST

Authors: Raghvendra Sahai, Valentin Bujarrabal, Guillermo Quintana-Lacaci, Nicole Reindl, Griet Van de Steene, Carmen Sánchez Contreras, Michael E. Ressler

Abstract: The planetary nebula (PN) NGC3132 is a striking example of the dramatic but poorly understood, mass-loss phenomena that (1-8) Msun stars undergo during their death throes as they evolve into white dwarfs (WDs). From an analysis of JWST multiwavelength (0.9-18 micron) imaging of NGC3132, we report the discovery of an asymmetrical dust cloud around the WD central star (CS) of NGC3132, seen most prom… ▽ More The planetary nebula (PN) NGC3132 is a striking example of the dramatic but poorly understood, mass-loss phenomena that (1-8) Msun stars undergo during their death throes as they evolve into white dwarfs (WDs). From an analysis of JWST multiwavelength (0.9-18 micron) imaging of NGC3132, we report the discovery of an asymmetrical dust cloud around the WD central star (CS) of NGC3132, seen most prominently in the 18 micron~image, with a surface-brightness limited radial extent of >~2 arcsec. We show that the A2V star located 1.7 arcsec to CS's North-East (and 0.75 kpc from Earth) is gravitationally-bound to the latter, by the detection of relative orbital angular motion of (0.24+/-0.045) deg between these stars over ~20 yr. Using aperture photometry of the CS extracted from the JWST images, together with published optical photometry and an archival UV spectrum, we have constructed the spectral-energy distribution (SED) of the CS and its extended emission over the UV to mid-IR (0.091-18 micron) range. We find that fitting the SED of the CS and the radial intensity distributions at 7.7, 12.8 and 18 micron with thermal emission from dust requires a cloud that extends to a radius of >~1785 au, with a dust mass of ~1.3 x 10^(-2) M(Earth) and grains that are 70% silicate and 30% amorphous carbon. We propose plausible origins of the dust cloud and an evolutionary scenario in which a system of three stars -- the CS, a close low-mass companion, and a more distant A2V star -- forms a stable hierarchical triple system on the main-sequence but becomes dynamically unstable later, resulting in the spectacular mass-ejections that form the current, multipolar PN. △ Less

Submitted 5 December, 2022; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 23 pages, 5 tables, 8 figures (ApJ, in press) references added, minor revision of Fig 5 caption, minor typos corrected

Consistent and simultaneous modelling of galaxy clustering and galaxy-galaxy lensing with Subhalo Abundance Matching

Authors: Sergio Contreras, Raul E. Angulo, Jonás Chaves-Montero, Simon D. M. White, Giovanni Aricò

Abstract: The spatial distribution of galaxies and their gravitational lensing signal offer complementary tests of galaxy formation physics and cosmology. However, their synergy can only be fully exploited if both probes are modelled accurately and consistently. In this paper, we demonstrate that this can be achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe. Specifically, we use mock c… ▽ More The spatial distribution of galaxies and their gravitational lensing signal offer complementary tests of galaxy formation physics and cosmology. However, their synergy can only be fully exploited if both probes are modelled accurately and consistently. In this paper, we demonstrate that this can be achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe. Specifically, we use mock catalogues built from the TNG300 hydrodynamical simulation to show that SHAMe can simultaneously model the multipoles of the redshift-space galaxy correlation function and galaxy-galaxy lensing, without noticeable bias within the statistical sampling uncertainties of a SDSS volume and on scales r = [0.6-30] Mpc/h. Modelling the baryonic processes in galaxy-galaxy lensing with a baryonification scheme allows SHAMe's range of validity to be extended to r = [0.1-30] Mpc/h. Remarkably, our model achieves this level of precision with just five free parameters beyond those describing the baryonification model. At fixed cosmology, we find that galaxy-galaxy lensing provides a general consistency test but little additional information on galaxy modelling parameters beyond that encoded in the redshift-space multipoles. It does, however, improve constraints if only the projected correlation function is available, as in surveys with only photometric redshifts. We expect SHAMe to have a higher fidelity across a wider range of scales than more traditional methods such as Halo Occupation Distribution modelling. Thus it should provide a significantly more powerful and more robust tool for analysing next-generation large-scale surveys. △ Less

Submitted 10 January, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 14 pages, 7 figures. Accepted by MNRAS

The galaxy formation origin of the lensing is low problem

Authors: Jonas Chaves-Montero, Raul E. Angulo, Sergio Contreras

Abstract: Recent analyses show that $Λ$CDM-based models optimised to reproduce the clustering of massive galaxies overestimate their gravitational lensing by about 30\%, the so-called lensing is low problem. Using a state-of-the-art hydrodynamical simulation, we show that this discrepancy reflects shortcomings in standard galaxy-halo connection models rather than tensions within the $Λ$CDM paradigm itself.… ▽ More Recent analyses show that $Λ$CDM-based models optimised to reproduce the clustering of massive galaxies overestimate their gravitational lensing by about 30\%, the so-called lensing is low problem. Using a state-of-the-art hydrodynamical simulation, we show that this discrepancy reflects shortcomings in standard galaxy-halo connection models rather than tensions within the $Λ$CDM paradigm itself. Specifically, this problem results from ignoring a variety of galaxy formation effects, including assembly bias, segregation of satellite galaxies relative to dark matter, and baryonic effects on the matter distribution. All these effects contribute towards overestimating gravitational lensing and, when combined, explain the amplitude and scale dependence of the lensing is low problem. We conclude that simplistic galaxy-halo connection models are inadequate to interpret clustering and lensing simultaneously, and that it is crucial to employ more sophisticated models for the upcoming generation of large-scale surveys. △ Less

Submitted 19 January, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 16 pages, 12 figures, accepted by MNRAS

The MillenniumTNG Project: Inferring cosmology from galaxy clustering with accelerated N-body scaling and subhalo abundance matching

Authors: Sergio Contreras, Raul E. Angulo, Volker Springel, Simon D. M. White, Boryana Hadzhiyska, Lars Hernquist, Rüdiger Pakmor, Rahul Kannan, César Hernández-Aguayo, Monica Barrera, Fulvio Ferlito, Ana Maria Delgado, Sownak Bose, Carlos Frenk

Abstract: We introduce a novel technique for constraining cosmological parameters and galaxy assembly bias using non-linear redshift-space clustering of galaxies. We scale cosmological N-body simulations and insert galaxies with the SubHalo Abundance Matching extended (SHAMe) empirical model to generate over 175,000 clustering measurements spanning all relevant cosmological and SHAMe parameter values. We th… ▽ More We introduce a novel technique for constraining cosmological parameters and galaxy assembly bias using non-linear redshift-space clustering of galaxies. We scale cosmological N-body simulations and insert galaxies with the SubHalo Abundance Matching extended (SHAMe) empirical model to generate over 175,000 clustering measurements spanning all relevant cosmological and SHAMe parameter values. We then build an emulator capable of reproducing the projected galaxy correlation function at the monopole, quadrupole and hexadecapole level for separations between $0.1\,h^{-1}{\rm Mpc}$ and $25\,h^{-1}{\rm Mpc}$. We test this approach by using the emulator and Monte Carlo Markov Chain (MCMC) inference to jointly estimate cosmology and assembly bias parameters both for the MTNG740 hydrodynamic simulation and for a semi-analytical galaxy formation model (SAM) built on the MTNG740-DM dark matter-only simulation, obtaining unbiased results for all cosmological parameters. For instance, for MTNG740 and a galaxy number density of $n\sim 0.01 h^{3}{\rm Mpc}^{-3}$, we obtain $σ_{8}=0.799^{+0.039}_{-0.044}$ ($σ_{8,{\rm MTNG}} =$ 0.8159), and $Ω_\mathrm{M}h^2= 0.138^{+ 0.025}_{- 0.018}$ ($Ω_{\mathrm{M}} h^2_{\rm MTNG} =$ 0.142). For fixed Hubble parameter ($h$), the constraint becomes $Ω_\mathrm{M}h^2= 0.137^{+ 0.011}_{- 0.012}$. Our method performs similarly well for the SAM and for other tested sample densities. We almost always recover the true amount of galaxy assembly bias within one sigma. The best constraints are obtained when scales smaller than $2\,h^{-1}{\rm Mpc}$ are included, as well as when at least the projected correlation function and the monopole are incorporated. These methods offer a powerful way to constrain cosmological parameters using galaxy surveys. △ Less

Submitted 13 December, 2022; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 18 pages, 12 figures, accepted by MNRAS

The MillenniumTNG Project: An improved two-halo model for the galaxy-halo connection of red and blue galaxies

Authors: Boryana Hadzhiyska, Daniel Eisenstein, Lars Hernquist, Rüdiger Pakmor, Sownak Bose, Ana Maria Delgado, Sergio Contreras, Rahul Kannan, Simon D. M. White, Volker Springel, Carlos Frenk, César Hernández-Aguayo, Fulvio Ferlito, Monica Barrera

Abstract: Approximate methods to populate dark matter halos with galaxies are of great utility to large galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs,… ▽ More Approximate methods to populate dark matter halos with galaxies are of great utility to large galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs, $z=1$ and $z=0$, in the large volume, high-resolution hydrodynamical simulation of the MillenniumTNG project. In a partner study we concentrated on the small-scale, one-halo regime down to $r\sim 0.1{\rm Mpc}/h$, while here we focus on modeling galaxy assembly bias in the two-halo regime, $r\gtrsim 1{\rm Mpc}/h$. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales ($r\sim 20{\rm Mpc}/h$), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10-15\% discrepancy present in the standard halo model prescription is only reconciled when we augment our halo occupation model with a dependence on extrinsic halo properties ("shear" being the best-performing one) rather than intrinsic ones (e.g., concentration, peak mass). We argue that this fact constitutes evidence for two-halo galaxy conformity. Including tertiary assembly bias (i.e. a property beyond mass and "shear") is not an essential requirement for reconciling the galaxy assembly bias signal of LRGs, but the combination of external and internal properties is beneficial for recovering ELG the clustering. We find that centrals in low-mass haloes dominate the assembly bias signal of both populations. Finally, we explore the predictions of our model for higher-order statistics such as nearest-neighbor counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 16 pages, 8 figures; submitted to MNRAS

The MillenniumTNG Project: Refining the one-halo model of red and blue galaxies at different redshifts

Authors: Boryana Hadzhiyska, Lars Hernquist, Daniel Eisenstein, Ana Maria Delgado, Sownak Bose, Rahul Kannan, Rüdiger Pakmor, Volker Springel, Sergio Contreras, Monica Barrera, Fulvio Ferlito, César Hernández-Aguayo, Simon D. M. White, Carlos Frenk

Abstract: Luminous red galaxies (LRGs) and blue star-forming emission-line galaxies (ELGs) are key tracers of large-scale structure used by cosmological surveys. Theoretical predictions for such data are often done via simplistic models for the galaxy-halo connection. In this work, we use the large, high-fidelity hydrodynamical simulation of the MillenniumTNG project (MTNG) to inform a new phenomenological… ▽ More Luminous red galaxies (LRGs) and blue star-forming emission-line galaxies (ELGs) are key tracers of large-scale structure used by cosmological surveys. Theoretical predictions for such data are often done via simplistic models for the galaxy-halo connection. In this work, we use the large, high-fidelity hydrodynamical simulation of the MillenniumTNG project (MTNG) to inform a new phenomenological approach for obtaining an accurate and flexible galaxy-halo model on small scales. Our aim is to study LRGs and ELGs at two distinct epochs, $z = 1$ and $z = 0$, and recover their clustering down to very small scales, $r \sim 0.1 \ {\rm Mpc}/h$, i.e. the one-halo regime, while a companion paper extends this to a two-halo model for larger distances. The occupation statistics of ELGs in MTNG inform us that: (1) the satellite occupations exhibit a slightly super-Poisson distribution, contrary to commonly made assumptions, and (2) that haloes containing at least one ELG satellite are twice as likely to host a central ELG. We propose simple recipes for modeling these effects, each of which calls for the addition of a single free parameter to simpler halo occupation models. To construct a reliable satellite population model, we explore the LRG and ELG satellite radial and velocity distributions and compare them with those of subhalos and particles in the simulation. We find that ELGs are anisotropically distributed within halos, which together with our occupation results provides strong evidence for cooperative galaxy formation (manifesting itself as one-halo galaxy conformity); i.e.~galaxies with similar properties form in close proximity to each other. Our refined galaxy-halo model represents a useful improvement of commonly used analysis tools and thus can be of help to increase the constraining power of large-scale structure surveys. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 14 pages, 10 figures; submitted to MNRAS

The cosmology dependence of the concentration-mass-redshift relation

Authors: Daniel López-Cano, Raúl E. Angulo, Aaron D. Ludlow, M. Zennaro, S. Contreras, Jonás Chaves-Montero, G. Aricò

Abstract: The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological paramet… ▽ More The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological parameters were varied: $σ_{8}$, $Ω_{\mathrm{M}}$, $Ω_{\mathrm{b}}$, $n_{\mathrm{s}}$, $h$, $M_ν$, $w_{0}$ and $w_{\mathrm{a}}$. We characterize the impact of these parameters on concentrations for different halo masses and redshifts. In agreement with previous works, and for all cosmologies studied, we find that there exists a tight correlation between the characteristic densities of dark matter haloes within their scale radii, $r_{-2}$, and the critical density of the Universe at a suitably defined formation time. This finding, when combined with excursion set modelling of halo formation histories, allows us to accurately predict the concentrations of dark matter haloes as a function of mass, redshift, and cosmology. We use our simulations to test the reliability of a number of published models for predicting halo concentration and highlight when they succeed or fail to reproduce the cosmological $c(M,z)$ relation. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 11 pages, 9 figures

The Bacco Simulation Project: Bacco Hybrid Lagrangian Bias Expansion Model in Redshift Space

Authors: Marcos Pellejero-Ibanez, Raul E. Angulo, Matteo Zennaro, Jens Stuecker, Sergio Contreras, Giovanni Arico, Francisco Maion

Abstract: We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We b… ▽ More We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We build the emulator using a forward neural network trained with the simulations of the BACCO project, which covers an 8-dimensional parameter space including massive neutrinos and dynamical dark energy. We show that our emulator provides unbiased cosmological constraints from the monopole, quadrupole, and hexadecapole of a mock galaxy catalogue that mimics the BOSS-CMASS sample down to nonlinear scales ($k\sim0.6$[$h/$Mpc]$^{3}$). This work opens up the possibility of robustly extracting cosmological information from small scales using observations of the large-scale structure of the Universe. △ Less

Submitted 27 July, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 16 pages, 9 figures

Dissecting the central regions of OH 231.8 + 4.2 with ALMA: a salty rotating disk at the base of a young bipolar outflow

Authors: C. Sánchez Contreras, J. Alcolea, R. Rodríguez Cardoso, V. Bujarrabal, A. Castro-Carrizo, G. Quintana-Lacaci, L. Velilla-Prieto, M. Santander-García

Abstract: We present ALMA continuum and molecular line emission maps at $\sim$1 mm of OH 231.8, a well studied bipolar nebula around an AGB star. The excellent angular resolution of our maps ($\sim$20 mas) allows us to scrutinise the central nebular regions of OH 231.8, which hold the clues to unravel how this iconic object assembled its complex nebular architecture. We report, for the first time in this ob… ▽ More We present ALMA continuum and molecular line emission maps at $\sim$1 mm of OH 231.8, a well studied bipolar nebula around an AGB star. The excellent angular resolution of our maps ($\sim$20 mas) allows us to scrutinise the central nebular regions of OH 231.8, which hold the clues to unravel how this iconic object assembled its complex nebular architecture. We report, for the first time in this object and others of its kind, the discovery of a rotating circumbinary disk selectively traced by NaCl, KCl, and H$_2$O emission lines. This represents the first detection of KCl in an oxygen-rich AGB circumstellar envelope. The rotating disk, of radius $\sim$30 au, lies at the base of a young bipolar wind traced by SiO and SiS emission, which also presents signs of rotation at its base. The NaCl equatorial structure is characterised by a mean rotation velocity of $\sim$4 km s$^{-1}$ and extremely low expansion speeds, $\sim$3 km s$^{-1}$. The outflow has a predominantly expansive kinematics characterised by a constant radial velocity gradient of $\sim$65 km s$^{-1}$ arcsec$^{-1}$ at its base. Beyond $r$$\sim$350 au, the gas in the outflow continues radially flowing at a constant terminal speed of $\sim$16 km s$^{-1}$. Our continuum maps reveal a spatially resolved dust disk-like structure perpendicular to the outflow, with the NaCl, KCl and H$_2$O emission arising from the disk's surface layers. Within the disk, we also identify an unresolved point continuum source, which likely represents the central Mira-type star QX Pup enshrouded by a $\sim$3 $R_{\star}$ component of hot ($\sim$1400 K) freshly formed dust. The point source is slightly off-centered from the disk centroid, enabling us for the first time to place constraints to the orbital separation and period of the central binary system, $a$$\sim$20 au and $P_{\rm orb}$$\sim$55 yr, respectively. (abridged). △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 30 pages, 18 figures. Accepted for publication in Astronomy & Astrophysics (on 19/06/2022)

Journal ref: A&A 665, A88 (2022)

Observational identification of a sample of likely recent Common-Envelope Events

Authors: Theo Khouri, Wouter Vlemmings, Daniel Tafoya, Andrés F. Pérez-Sánchez, Carmen Sánchez Contreras, José F. Gómez, Hiroshi Imai, Raghvendra Sahai

Abstract: One of the most poorly understood stellar evolutionary paths is that of binary systems undergoing common-envelope evolution, when the envelope of a giant star engulfs the orbit of a companion. Although this interaction leads to a great variety of astrophysical systems, direct empirical studies are difficult because few objects experiencing common-envelope evolution are known. We present ALMA obser… ▽ More One of the most poorly understood stellar evolutionary paths is that of binary systems undergoing common-envelope evolution, when the envelope of a giant star engulfs the orbit of a companion. Although this interaction leads to a great variety of astrophysical systems, direct empirical studies are difficult because few objects experiencing common-envelope evolution are known. We present ALMA observations towards sources known as water fountains that reveal they had low initial masses ($<4~{\rm M}_\odot$) and ejected a significant fraction of it over less than a few hundred years. The only mechanism able to explain such rapid mass ejection is common-envelope evolution. Our calculations show that the water-fountain sample accounts for a large fraction of the systems in our Galaxy which have just experienced the common-envelope phase. Since water-fountain sources show characteristic fast bipolar outflows, outflows and jets likely play an important role right before, during or immediately after the common-envelope phase. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Comments: Published on Nature Astronomy on Dec 16, 2021: https://rdcu.be/cDlX8

Priors on Lagrangian bias parameters from galaxy formation modelling

Authors: Matteo Zennaro, Raul E. Angulo, Sergio Contreras, Marcos Pellejero-Ibáñez, Francisco Maion

Abstract: We study the relations among the parameters of the hybrid Lagrangian bias expansion model, fitting biased auto and cross power spectra up to $k_{\rm max} = 0.7 \, h \, \mathrm{Mpc}^{-1}$. We consider $\sim 8000$ halo and galaxy samples, with different halo masses, redshifts, galaxy number densities, and varying the parameters of the galaxy formation model. Galaxy samples are obtained through state… ▽ More We study the relations among the parameters of the hybrid Lagrangian bias expansion model, fitting biased auto and cross power spectra up to $k_{\rm max} = 0.7 \, h \, \mathrm{Mpc}^{-1}$. We consider $\sim 8000$ halo and galaxy samples, with different halo masses, redshifts, galaxy number densities, and varying the parameters of the galaxy formation model. Galaxy samples are obtained through state-of-the-art extended subhalo abundance matching techniques and include both stellar-mass and star-formation-rate selected galaxies. All of these synthetic galaxies samples are publicly available at https://bacco.dipc.org/galpk.html. We find that the hybrid Lagrangian bias model provides accurate fits to all of our halo and galaxy samples. The coevolution relations between galaxy bias parameters, although roughly compatible with those obtained for haloes, show systematic shifts and larger scatter. We explore possible sources of this difference in terms of dependence on halo occupation and assembly bias of each sample. The bias parameter relations displayed in this work can be used as a prior for future Bayesian analyses employing the hybrid Lagrangian bias expansion model. △ Less

Submitted 27 June, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 14 pages, 12 figures, galaxy power spectra available at https://bacco.dipc.org/galpk.html

Modelling galaxy clustering in redshift space with a Lagrangian bias formalism and $N$-body simulations

Authors: Marcos Pellejero-Ibanez, Jens Stuecker, Raul E. Angulo, Matteo Zennaro, Sergio Contreras, Giovanni Arico

Abstract: Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describ… ▽ More Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describe redshift-space distortions even down to small scales. This model corresponds to a second-order perturbative Lagrangian bias expansion which is advected to Eulerian space employing a displacement field extracted from $N$-body simulations. Eulerian coordinates are then transformed into redshift space by directly employing simulated velocity fields augmented with nuisance parameters capturing various possible satellite fractions and intra-halo small-scale velocities. We quantify the accuracy of our approach against samples of physically-motivated mock galaxies selected according to either Stellar Mass (SM) or Star Formation Rate (SFR) at multiple abundances and at $z=0$ and $1$. We find our model describes the monopole, quadrupole, and hexadecapole of the galaxy-power spectra down to scales of $k\approx 0.6 [h/$Mpc] within the accuracy of our simulations. This approach could pave the way to significantly increase the amount of cosmological information to be extracted from future galaxy surveys. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: 15 pages, 8 figures

Predicting halo occupation and galaxy assembly bias with machine learning

Authors: Xiaoju Xu, Saurabh Kumar, Idit Zehavi, Sergio Contreras

Abstract: Understanding the impact of halo properties beyond halo mass on the clustering of galaxies (namely galaxy assembly bias) remains a challenge for contemporary models of galaxy clustering. We explore the use of machine learning to predict the halo occupations and recover galaxy clustering and assembly bias in a semi-analytic galaxy formation model. For stellar-mass selected samples, we train a Rando… ▽ More Understanding the impact of halo properties beyond halo mass on the clustering of galaxies (namely galaxy assembly bias) remains a challenge for contemporary models of galaxy clustering. We explore the use of machine learning to predict the halo occupations and recover galaxy clustering and assembly bias in a semi-analytic galaxy formation model. For stellar-mass selected samples, we train a Random Forest algorithm on the number of central and satellite galaxies in each dark matter halo. With the predicted occupations, we create mock galaxy catalogues and measure the clustering and assembly bias. Using a range of halo and environment properties, we find that the machine learning predictions of the occupancy variations with secondary properties, galaxy clustering and assembly bias are all in excellent agreement with those of our target galaxy formation model. Internal halo properties are most important for the central galaxies prediction, while environment plays a critical role for the satellites. Our machine learning models are all provided in a usable format. We demonstrate that machine learning is a powerful tool for modelling the galaxy-halo connection, and can be used to create realistic mock galaxy catalogues which accurately recover the expected occupancy variations, galaxy clustering and galaxy assembly bias, imperative for cosmological analyses of upcoming surveys. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 21 pages, 17 figures, submitted to MNRAS

The cosmological dependence of halo and galaxy assembly bias

Authors: Sergio Contreras, Jonás Chaves-Montero, Matteo Zennaro, Raul E. Angulo

Abstract: One of the main predictions of excursion set theory is that the clustering of dark matter haloes only depends on halo mass. However, it has been long established that the clustering of haloes also depends on other properties, including formation time, concentration, and spin; this effect is commonly known as halo assembly bias. We use a suite of gravity-only simulations to study the dependence of… ▽ More One of the main predictions of excursion set theory is that the clustering of dark matter haloes only depends on halo mass. However, it has been long established that the clustering of haloes also depends on other properties, including formation time, concentration, and spin; this effect is commonly known as halo assembly bias. We use a suite of gravity-only simulations to study the dependence of halo assembly bias on cosmology; these simulations cover cosmological parameters spanning 10$σ$ around state-of-the-art best-fitting values, including standard extensions of the $Λ$CDM paradigm such as neutrino mass and dynamical dark energy. We find that the strength of halo assembly bias presents variations smaller than 0.05 dex across all cosmologies studied for concentration and spin selected haloes, letting us conclude that the dependence of halo assembly bias upon cosmology is negligible. We then study the dependence of galaxy assembly bias (i.e. the manifestation of halo assembly bias in galaxy clustering) on cosmology using subhalo abundance matching. We find that galaxy assembly bias also presents very small dependence upon cosmology ($\sim$ 2$\%$-4$\%$ of the total clustering); on the other hand, we find that the dependence of this signal on the galaxy formation parameters of our galaxy model is much stronger. Taken together, these results let us conclude that the dependence of halo and galaxy assembly bias on cosmology is practically negligible. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: 11 pages, 8 figures. Submitted to MNRAS, comments welcome

Non-universality of the mass function: dependence on the growth rate and power spectrum shape

Authors: Lurdes Ondaro-Mallea, Raul E. Angulo, Matteo Zennaro, Sergio Contreras, Giovanni Aricò

Abstract: The abundance of dark matter haloes is one of the key probes of the growth of structure and expansion history of the Universe. Theoretical predictions for this quantity usually assume that, when expressed in a certain form, it depends only on the mass variance of the linear density field. However, cosmological simulations have revealed that this assumption breaks, leading to 10-20% systematic effe… ▽ More The abundance of dark matter haloes is one of the key probes of the growth of structure and expansion history of the Universe. Theoretical predictions for this quantity usually assume that, when expressed in a certain form, it depends only on the mass variance of the linear density field. However, cosmological simulations have revealed that this assumption breaks, leading to 10-20% systematic effects. In this paper we employ a specially-designed suite of simulations to further investigate this problem. Specifically, we carry out cosmological N-body simulations where we systematically vary growth history at a fixed linear density field, or vary the power spectrum shape at a fixed growth history. We show that the halo mass function generically depends on these quantities, thus showing a clear signal of non-universality. Most of this effect can be traced back to the way in which the same linear fluctuation grows differently into the nonlinear regime depending on details of its assembly history. With these results, we propose a parameterization with explicit dependence on the linear growth rate and power spectrum shape. Using an independent suite of simulations, we show that this fitting function accurately captures the mass function of haloes over cosmologies spanning a vast parameter space, including massive neutrinos and dynamical dark energy. Finally, we employ this tool to improve the accuracy of so-called cosmology-rescaling methods and show they can deliver 2% accurate predictions for the halo mass function over the whole range of currently viable cosmologies. △ Less

Submitted 17 February, 2021; originally announced February 2021.

The BACCO simulation project: biased tracers in real space

Authors: Matteo Zennaro, Raul E. Angulo, Marcos Pellejero-Ibáñez, Jens Stücker, Sergio Contreras, Giovanni Aricò

Abstract: We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advecte… ▽ More We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advected to Eulerian space using the displacement field of numerical simulations. The cosmology-dependence is captured thanks to a cosmology-rescaling algorithm. We show that our emulator is capable of describing the power spectrum of galaxy formation simulations for a sample mimicking that of a typical Emission-Line survey at $z \sim 1$ with an accuracy of $1-2\%$ up to nonlinear scales $k \sim 0.7 h \mathrm{Mpc}^{-1}$. △ Less

Submitted 2 November, 2023; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: 14 pages, 8 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 2, pp.2407-2419, 2023

SHAM through the lens of a hydrodynamical simulation

Authors: Ginevra Favole, Antonio D. Montero-Dorta, M. Celeste Artale, Sergio Contreras, Idit Zehavi, Xiaoju Xu

Abstract: We use the IllustrisTNG100 hydrodynamical simulation to study the dependence of the galaxy two-point correlation function on a broad range of secondary subhalo and galactic properties. We construct galaxy mock catalogues adopting a standard sub-halo abundance matching scheme coupled with a secondary assignment between galaxy colour or specific star formation rate and the following subhalo properti… ▽ More We use the IllustrisTNG100 hydrodynamical simulation to study the dependence of the galaxy two-point correlation function on a broad range of secondary subhalo and galactic properties. We construct galaxy mock catalogues adopting a standard sub-halo abundance matching scheme coupled with a secondary assignment between galaxy colour or specific star formation rate and the following subhalo properties: starvation redshift z$_{\rm starve}$, concentration at infall, overdensity $δ_R^{\rm env}$, tidal anisotropy $α_R$, and tidal overdensity $δ_R$. The last two quantities allow us to fully characterise the tidal field of our subhaloes, acting as mediators between their internal and large-scale properties. The resulting mock catalogues overall return good agreement with the IllustrisTNG100 measurements. The accuracy of each model strongly depends on the correlation between the secondary galaxy and subhalo properties employed. Among all the subhalo proxies tested, we find that z$_{\rm starve}$ and $c_{\rm infall}$ are the ones that best trace the large-scale structure, producing robust clustering predictions for different samples of red/blue and quenched/star-forming galaxies. △ Less

Submitted 3 November, 2021; v1 submitted 26 January, 2021; originally announced January 2021.

Comments: 12 pages, 8 figures

A flexible subhalo abundance matching model for galaxy clustering in redshift space

Authors: Sergio Contreras, Raul Angulo, Matteo Zennaro

Abstract: We develop an extension of subhalo abundance matching (SHAM) capable of accurately reproducing the real and redshift-space clustering of galaxies in a state-of-the-art hydrodynamical simulation. Our method uses a low-resolution gravity-only simulation and it includes orphan and tidal disruption prescriptions for satellite galaxies, and a flexible amount of galaxy assembly bias. Furthermore, it inc… ▽ More We develop an extension of subhalo abundance matching (SHAM) capable of accurately reproducing the real and redshift-space clustering of galaxies in a state-of-the-art hydrodynamical simulation. Our method uses a low-resolution gravity-only simulation and it includes orphan and tidal disruption prescriptions for satellite galaxies, and a flexible amount of galaxy assembly bias. Furthermore, it includes recipes for star formation rate (SFR) based on the dark matter accretion rate. We test the accuracy of our model against catalogues of stellar-mass- and SFR-selected galaxies in the TNG300 hydrodynamic simulation. By fitting a small number of free parameters, our extended SHAM reproduces the projected correlation function and redshift-space multipoles for number densities $10^{-3} - 10^{-2}\, h^{3}{\rm Mpc}^{-3}$, at $z=1$ and $z=0$, and for scales $r \in [0.3 - 20] h^{-1}{\rm Mpc}$. Simultaneously, the SHAM results also retrieve the correct halo occupation distribution, the level of galaxy assembly bias, and higher-order statistics present in the TNG300 galaxy catalogues. As an application, we show that our model simultaneously fits the projected correlation function of the SDSS in 3 disjoint stellar mass bins, with an accuracy similar to that of TNG300 galaxies. This SHAM extension can be used to get accurate clustering prediction even when using low and moderate-resolution simulations. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 14 pages, 12 figures. Submitted to MNRAS

The BACCO Simulation Project: A baryonification emulator with Neural Networks

Authors: Giovanni Aricò, Raul E. Angulo, Sergio Contreras, Lurdes Ondaro-Mallea, Marcos Pellejero-Ibañez, Matteo Zennaro

Abstract: We present a neural-network emulator for baryonic effects in the non-linear matter power spectrum. We calibrate this emulator using more than 50,000 measurements in a 15-dimensional parameters space, varying cosmology and baryonic physics. Baryonic physics is described through a baryonification algorithm, that has been shown to accurately capture the relevant effects on the power spectrum and bisp… ▽ More We present a neural-network emulator for baryonic effects in the non-linear matter power spectrum. We calibrate this emulator using more than 50,000 measurements in a 15-dimensional parameters space, varying cosmology and baryonic physics. Baryonic physics is described through a baryonification algorithm, that has been shown to accurately capture the relevant effects on the power spectrum and bispectrum in state-of-the-art hydrodynamical simulations. Cosmological parameters are sampled using a cosmology-rescaling approach including massive neutrinos and dynamical dark energy. The specific quantity we emulate is the ratio between matter power spectrum with baryons and gravity-only, and we estimate the overall precision of the emulator to be 1-2%, at all scales 0.01 < k < 5 h/Mpc, and redshifts 0 < z < 1.5. We also obtain an accuracy of 1-2%, when testing the emulator against a collection of 74 different cosmological hydrodynamical simulations and their respective gravity-only counterparts. We show also that only one baryonic parameter, namely Mc, which set the gas fraction retained per halo mass, is enough to have accurate and realistic predictions of the baryonic feedback at a given epoch. Our emulator will become publicly available in http://www.dipc.org/bacco. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 13 pages, 11 figures. Comments are welcome

The assembly bias of emission line galaxies

Authors: Esteban Jimenez, Nelson Padilla, Sergio Contreras, Idit Zehavi, Carlton Baugh, Alvaro Orsi

Abstract: The next generation of spectroscopic surveys will target emission-line galaxies (ELGs) to produce constraints on cosmological parameters. We study the large scale structure traced by ELGs using a combination of a semi-analytical model of galaxy formation, a code that computes the nebular emission from HII regions using the properties of the interstellar medium, and a large-volume, high-resolution… ▽ More The next generation of spectroscopic surveys will target emission-line galaxies (ELGs) to produce constraints on cosmological parameters. We study the large scale structure traced by ELGs using a combination of a semi-analytical model of galaxy formation, a code that computes the nebular emission from HII regions using the properties of the interstellar medium, and a large-volume, high-resolution N-body simulation. We consider fixed number density samples where galaxies are selected by either their H$α$, [OIII]$λ5007$ or [OII]$λλ3727-3729$ emission line luminosities. We investigate the assembly bias signatures of these samples, and compare them to those of stellar mass and SFR selected samples. Interestingly, we find that the [OIII]- and [OII]-selected samples display scale-dependent bias on large scales and that their assembly bias signatures are also scale-dependent. Both these effects are more pronounced for lower number density samples. The [OIII] and [OII] emitters that contribute most to the scale dependence tend to have a low gas-phase metallicity and are preferentially found in low-density regions. We also measure the baryon acoustic oscillation (BAO) feature and the $β$ parameter related to the growth rate of overdensities. We find a slight tendency for the BAO peak to shift toward smaller scales for [OII] emitters and that $β$ is scale-dependent at large scales. Our results suggest that ELG samples include environmental effects that should be modelled in order to remove potential systematic errors that could affect the estimation of cosmological parameters. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Comments: 15 pages, 12 figures

Simultaneous modelling of matter power spectrum and bispectrum in the presence of baryons

Authors: Giovanni Aricò, Raul E. Angulo, Carlos Hernández-Monteagudo, Sergio Contreras, Matteo Zennaro

Abstract: We demonstrate that baryonification algorithms, which displace particles in gravity-only simulations according to physically-motivated prescriptions, can simultaneously capture the impact of baryonic physics on the 2 and 3-point statistics of matter. Specifically, we show that our implementation of a baryonification algorithm jointly fits the changes induced by baryons on the power spectrum and eq… ▽ More We demonstrate that baryonification algorithms, which displace particles in gravity-only simulations according to physically-motivated prescriptions, can simultaneously capture the impact of baryonic physics on the 2 and 3-point statistics of matter. Specifically, we show that our implementation of a baryonification algorithm jointly fits the changes induced by baryons on the power spectrum and equilateral bispectrum on scales up to k < 5 h/Mpc and redshifts z<2, as measured in six different cosmological hydrodynamical simulations. The accuracy of our fits are typically 1% for the power spectrum, and for the equilateral and squeezed bispectra, which somewhat degrades to 3% for simulations with extreme feedback prescriptions. Our results support the physical assumptions underlying baryonification approaches, and encourage their use in interpreting weak gravitational lensing and other cosmological observables. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 14 pages, 15 figures. To be submitted to MNRAS. Comments are welcome

Dissecting and Modelling Galaxy Assembly Bias

Authors: Xiaoju Xu, Idit Zehavi, Sergio Contreras

Abstract: Understanding the galaxy-halo connection is fundamental for contemporary models of galaxy clustering. The extent to which the haloes' assembly history and environment impact galaxy clustering (a.k.a. galaxy assembly bias; GAB), remains a complex and challenging problem. Using a semi-analytic galaxy formation model, we study the individual contributions of different secondary halo properties to the… ▽ More Understanding the galaxy-halo connection is fundamental for contemporary models of galaxy clustering. The extent to which the haloes' assembly history and environment impact galaxy clustering (a.k.a. galaxy assembly bias; GAB), remains a complex and challenging problem. Using a semi-analytic galaxy formation model, we study the individual contributions of different secondary halo properties to the GAB signal. These are obtained by comparing the clustering of stellar-mass selected samples to that of shuffled samples where the galaxies are randomly reassigned to haloes of fixed mass and a specified secondary halo property. We explore a large range of internal halo properties and environmental measures. We find that commonly used properties like halo age or concentration amount to only 20-30 per cent of the signal, while the smoothed matter density or the tidal anisotropy can account for the full level of GAB (though care should be given to the specific definition). For the "successful" measures, we examine the occupancy variations and the associated changes in the halo occupation function parameters. These are used to create mock catalogues that reproduce the full level of GAB. Finally, we propose a practical modification of the standard halo occupation distribution model, which can be tuned to any level of assembly bias. Fitting the parameters to our semi-analytic model, we demonstrate that the corresponding mock catalogue recovers the target level of GAB as well as the occupancy variations. Our results enable producing realistic mock catalogues and directly inform theoretical modelling of assembly bias and attempts to detect it in the Universe. △ Less

Submitted 2 April, 2021; v1 submitted 10 July, 2020; originally announced July 2020.

Comments: 22 pages, 16 figures, published by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2021, 502(3), 3242-3263

A flexible modelling of galaxy assembly bias

Authors: Sergio Contreras, Raul Angulo, Matteo Zennaro

Abstract: We use the {\sc Illustris TNG300} magneto-hydrodynamic simulation, the {\sc SAGE} semi-analytical model, and the subhalo abundance matching technique (SHAM) to examine the diversity in predictions for galaxy assembly bias (i.e. the difference in the large scale clustering of galaxies at a fixed halo mass due to correlations with the assembly history and other properties of host haloes). We conside… ▽ More We use the {\sc Illustris TNG300} magneto-hydrodynamic simulation, the {\sc SAGE} semi-analytical model, and the subhalo abundance matching technique (SHAM) to examine the diversity in predictions for galaxy assembly bias (i.e. the difference in the large scale clustering of galaxies at a fixed halo mass due to correlations with the assembly history and other properties of host haloes). We consider samples of galaxies selected according to their stellar mass or star formation rate at various redshifts. We find that all models predict an assembly bias signal of different magnitude, redshift evolution, and dependence with selection criteria and number density. To model these non-trivial dependences, we propose an extension to the standard SHAM technique so it can include arbitrary amounts of assembly bias. We do this by preferentially selecting subhaloes with the same internal property but different {\it individual} large-scale bias. We find that with this model, we can successfully reproduce the galaxy assembly bias signal in either {\sc SAGE} or the {\sc TNG}, for all redshifts and galaxy number densities. We anticipate that this model can be used to constrain the level of assembly bias in observations and aid in the creation of more realistic mock galaxy catalogues. △ Less

Submitted 21 April, 2021; v1 submitted 7 May, 2020; originally announced May 2020.

Comments: Accepted by mnras

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

Authors: Raul E. Angulo, Matteo Zennaro, Sergio Contreras, Giovanni Aricò, Marcos Pellejero-Ibañez, Jens Stücker

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 m… ▽ More 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<z<1.5$ and over scales $10^{-2} < k/(h Mpc^{-1} ) < 5$ for 800 points in an 8-dimensional cosmological parameter space. For an efficient interpolation of the results, we build an emulator and compare its predictions against several widely-used methods. Over the whole range of scales considered, we expect our predictions to be accurate at the 2\% level for parameters in the minimal $Λ$ CDM model and to 3\% when extended to dynamical dark energy and massive neutrinos. We make our emulator publicly available under http://www.dipc.org/bacco △ Less

Submitted 6 July, 2021; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 12 pages; 11 figures. Accepted for publication in MNRAS

Do model emission line galaxies live in filaments at z~1?

Authors: V. Gonzalez-Perez, W. Cui, S. Contreras, C. M. Baugh, J. Comparat, A. J. Griffin, J. Helly, A. Knebe, C. Lacey, P. Norberg

Abstract: Current and future cosmological surveys are targeting star-forming galaxies at $z\sim 1$ with nebular emission lines. We use a state-of-the-art semi-analytical model of galaxy formation and evolution to explore the large scale environment of star-forming emission line galaxies (ELGs). Model ELGs are selected such that they can be compared directly with the DEEP2, VVDS, eBOSS-SGC and DESI surveys.… ▽ More Current and future cosmological surveys are targeting star-forming galaxies at $z\sim 1$ with nebular emission lines. We use a state-of-the-art semi-analytical model of galaxy formation and evolution to explore the large scale environment of star-forming emission line galaxies (ELGs). Model ELGs are selected such that they can be compared directly with the DEEP2, VVDS, eBOSS-SGC and DESI surveys. The large scale environment of the ELGs is classified using velocity-shear-tensor and tidal-tensor algorithms. Half of the model ELGs live in filaments and about a third in sheets. Model ELGs which reside in knots have the largest satellite fractions. We find that the shape of the mean halo occupation distribution of model ELGs varies widely for different large scale environments. To interpret our results, we also study fixed number density samples of ELGs and galaxies selected using simpler criteria, with single cuts in stellar mass, star formation rate and [OII] luminosity. The fixed number density ELG selection produces samples that are close to L[OII] and SFR selected samples for densities above $10^{-4.2}h^{3}{\rm Mpc}^{-3}$. ELGs with an extra cut in stellar mass applied to fix their number density, present differences in sheets and knots with respect to the other samples. ELGs, SFR and L[OII] selected samples with equal number density have similar large scale bias but their clustering below separations of $1h^{-1}$Mpc is different. △ Less

Submitted 14 August, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 20 pages, 17 figures, MNRAS in press, a small fraction of the presented data can be found here: https://github.com/viogp/plots4papers/tree/master/elg_cw_plots

The better half -- Asymmetric star-formation due to ram pressure in the EAGLE simulations

Authors: P. Troncoso-Iribarren, N. Padilla, C. Santander, C. D. P. Lagos, D. García-Lambas, S. Rodriguez, S. Contreras

Abstract: We use the EAGLE simulations to study the effects of the intra-cluster medium (ICM) on the spatially resolved star-formation activity in galaxies. We study three cases of galaxy asymmetry dividing each galaxy in two halves using the plane (i) perpendicular to the $\rm \texttt{velocity}$ direction, differentiating the galaxy part approaching to the cluster center, hereafter dubbed as the "leading h… ▽ More We use the EAGLE simulations to study the effects of the intra-cluster medium (ICM) on the spatially resolved star-formation activity in galaxies. We study three cases of galaxy asymmetry dividing each galaxy in two halves using the plane (i) perpendicular to the $\rm \texttt{velocity}$ direction, differentiating the galaxy part approaching to the cluster center, hereafter dubbed as the "leading half", and the opposite one "trailing half", (ii) perpendicular to the $\rm \texttt{radial}$ position of the satellite to the centre of the cluster, (iii) that maximizes the star-formation rate ($\rm SFR$) difference between the two halves. For (i), we find an enhancement of the $\rm SFR$, star formation efficiency ($\rm SFE$), and interstellar medium pressure in the leading half with respect to the trailing one and normal star-forming galaxies in the EAGLE simulation, and a clear overabundance of gas particles in their trailing. These results suggest that ram pressure (RP) is boosting the star formation by gas compression in the leading half, and transporting the gas to the trailing half. This effect is more pronounced in satellites of intermediate stellar masses $\rm 10^{9.5-10.5} M_{\odot}$, with gas masses above $\rm 10^{9} M_{\odot}$, and located within one virial radius or in the most massive clusters. In (iii) we find an alignment between the velocity and the vector perpendicular to the plane that maximizes the $\rm SFR$ difference between the two halves. It suggests that finding this plane in real galaxies can provide insights into the velocity direction. △ Less

Submitted 3 February, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

3\%-accurate predictions for the clustering of dark matter, haloes and subhaloes, over a wide range of cosmologies and scales

Authors: S. Contreras, R. E. Angulo, M. Zennaro, G. Aricó, M. Pellejero-Ibañez

Abstract: Predicting the spatial distribution of objects as a function of cosmology is an essential ingredient for the exploitation of future galaxy surveys. In this paper we show that a specially-designed suite of gravity-only simulations together with cosmology-rescaling algorithms can provide the clustering of dark matter, haloes, and subhaloes with high precision. Specifically, with only 3 $N$-body simu… ▽ More Predicting the spatial distribution of objects as a function of cosmology is an essential ingredient for the exploitation of future galaxy surveys. In this paper we show that a specially-designed suite of gravity-only simulations together with cosmology-rescaling algorithms can provide the clustering of dark matter, haloes, and subhaloes with high precision. Specifically, with only 3 $N$-body simulations we obtain the power spectrum of dark matter at $z=0$ and $z=1$ to better than 3\% precision for essentially all currently viable values of 8 cosmological parameters, including massive neutrinos and dynamical dark energy, and over the whole range of scales explored, 0.03 < $k/h^{-1}Mpc$ < 5. This precision holds at the same level for mass-selected haloes and for subhaloes selected according to their peak maximum circular velocity. As an initial application of these predictions, we successfully constrain $Ω_{\rm m}$, $σ_8$, and the scatter in subhalo-abundance-matching employing the projected correlation function of mock SDSS galaxies. △ Less

Submitted 22 September, 2020; v1 submitted 9 January, 2020; originally announced January 2020.

Comments: 12 pages, 8 figures. Submitted. Comments are welcome

Cosmological parameter estimation via iterative emulation of likelihoods

Authors: Marcos Pellejero-Ibañez, Raul E. Angulo, Giovanni Aricó, Matteo Zennaro, Sergio Contreras, Jens Stücker

Abstract: The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterativ… ▽ More The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterative Gaussian emulation of the target likelihood function. This requires a minimal number of likelihood evaluations and naturally accommodates for stochasticity in theoretical models. We apply the algorithm to estimate 9 parameters from the monopole and quadrupole of a mock power spectrum in redshift space. We obtain accurate posterior distribution functions with approximately 100 times fewer likelihood evaluations than an affine invariant MCMC, roughly independently from the dimensionality of the problem. We anticipate that our parameter estimation algorithm will accelerate the adoption of more accurate theoretical models in data analysis, enabling more comprehensive exploitation of cosmological observables. △ Less

Submitted 23 January, 2020; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: 12 pages, 5 figures

Modelling the large scale structure of the Universe as a function of cosmology and baryonic physics

Authors: Giovanni Aricò, Raul E. Angulo, Carlos Hernández-Monteagudo, Sergio Contreras, Matteo Zennaro, Marcos Pellejero-Ibañez, Yetli Rosas-Guevara

Abstract: We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons,… ▽ More We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons, such as star formation and AGN feedback. We show how this approach can accurately reproduce the effects of baryons on the matter power spectrum of various state-of-the-art hydro-dynamical simulations (EAGLE, Illustris, Illustris-TNG, Horizon-AGN, and OWLS,Cosmo-OWLS and BAHAMAS), to percent level from very large down to small, highly nonlinear scales, k= 5 h/Mpc, and from z=0 up to z=2. We highlight that, thanks to the heavy optimisation of the algorithms, we can obtain these predictions for arbitrary baryonic models and cosmology (including massive neutrinos and dynamical dark energy models) with an almost negligible CPU cost. Therefore, this approach is efficient enough for cosmological data analyses. With these tools in hand we explore the degeneracies between cosmological and astrophysical parameters in the nonlinear mass power spectrum. Our findings suggest that after marginalising over baryonic physics, cosmological constraints inferred from weak gravitational lensing should be moderately degraded. △ Less

Submitted 17 January, 2020; v1 submitted 19 November, 2019; originally announced November 2019.

Comments: 20 pages, 15 figures. Submitted to MNRAS

A rotating fast bipolar wind and disk system around the B[e]-type star MWC 922

Authors: C. Sánchez Contreras, A. Báez-Rubio, J. Alcolea, A. Castro-Carrizo, V. Bujarrabal, J. Martín-Pintado, D. Tafoya

Abstract: We present interferometric observations with the Atacama Large Millimeter Array (ALMA) of the free-free continuum and recombination line emission at 1 and 3mm of the "Red Square Nebula" surrounding the B[e]-type star MWC922. The unknown distance to the source is usually taken to be d=1.7-3 kpc. The unprecedented angular resolution (up to ~0.02arcsec) and exquisite sensitivity of these data unveil,… ▽ More We present interferometric observations with the Atacama Large Millimeter Array (ALMA) of the free-free continuum and recombination line emission at 1 and 3mm of the "Red Square Nebula" surrounding the B[e]-type star MWC922. The unknown distance to the source is usually taken to be d=1.7-3 kpc. The unprecedented angular resolution (up to ~0.02arcsec) and exquisite sensitivity of these data unveil, for the first time, the structure and kinematics of the emerging, compact ionized region at its center. We imaged the line emission of H30a and H39a, previously detected with single-dish observations, as well as of H51epsilon, H55gamma, and H63delta, detected for the first time in this work. The line emission is seen over a full velocity range of ~180 km/s arising in a region of diameter <0.14arcsec (less than a few hundred au) in the maser line H30a, which is the most intense transition reported here. We resolve the spatio-kinematic structure of a nearly edge-on disk rotating around a central mass of ~10Msun (d=1.7 kpc) or ~18Msun (d=3 kpc), assuming Keplerian rotation. Our data also unveil a fast (~100 km/s) bipolar ejection (a jet?) orthogonal to the disk. In addition, a slow (<15km/s) wind may be lifting off the disk. Both, the slow and the fast winds are found to be rotating in a similar manner to the ionized layers of the disk. This represents the first empirical proof of rotation in a bipolar wind expanding at high velocity (~100 km/s). (abridged) △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: Accecpted in A&A. 21 pages, 18 figures

Journal ref: A&A 629, A136 (2019)

On the prospect of using the maximum circular velocity of halos to encapsulate assembly bias in the galaxy-halo connection

Authors: Idit Zehavi, Stephen E. Kerby, Sergio Contreras, Esteban Jiménez, Nelson Padilla, Carlton M. Baugh

Abstract: We investigate a conceptual modification of the halo occupation distribution approach, using the halos' present-day maximal circular velocity, $\vmax$, as an alternative to halo mass. In particular, using a semi-analytic galaxy formation model applied to the Millennium WMAP7 simulation, we explore the extent that switching to $\vmax$ as the primary halo property incorporates the effects of assembl… ▽ More We investigate a conceptual modification of the halo occupation distribution approach, using the halos' present-day maximal circular velocity, $\vmax$, as an alternative to halo mass. In particular, using a semi-analytic galaxy formation model applied to the Millennium WMAP7 simulation, we explore the extent that switching to $\vmax$ as the primary halo property incorporates the effects of assembly bias into the formalism. We consider fixed number density galaxy samples ranked by stellar mass and examine the variations in the halo occupation functions with either halo concentration or formation time. We find that using $\vmax$ results in a significant reduction in the occupancy variation of the central galaxies, particularly for concentration. The satellites occupancy variation on the other hand increases in all cases. We find effectively no change in the halo clustering dependence on concentration, for fixed bins of $\vmax$ compared to fixed halo mass. Most crucially, we calculate the impact of assembly bias on galaxy clustering by comparing the amplitude of clustering to that of a shuffled galaxy sample, finding that the level of galaxy assembly bias remains largely unchanged. Our results suggest that while using $\vmax$ as a proxy for halo mass diminishes some of occupancy variations exhibited in the galaxy-halo relation, it is not able to encapsulate the effects of assembly bias potentially present in galaxy clustering. The use of other more complex halo properties, such as $\vpeak$, the peak value of $\vmax$ over the assembly history, provides some improvement and warrants further investigation. △ Less

Submitted 11 July, 2019; originally announced July 2019.

Comments: 16 figures, 15 pages, submitted to ApJ

Extensions to the halo occupation distribution model for more accurate clustering predictions

Authors: Esteban Jiménez, Sergio Contreras, Nelson Padilla, Idit Zehavi, Carlton M. Baugh, Violeta Gonzalez-Perez

Abstract: We test different implementations of the halo occupation distribution (HOD) model to reconstruct the spatial distribution of galaxies as predicted by a publicly available semi-analytical model (SAM). We compare the measured two-point correlation functions of the HOD mock catalogues and the SAM samples to quantify the fidelity of the reconstruction. We use fixed number density galaxy samples select… ▽ More We test different implementations of the halo occupation distribution (HOD) model to reconstruct the spatial distribution of galaxies as predicted by a publicly available semi-analytical model (SAM). We compare the measured two-point correlation functions of the HOD mock catalogues and the SAM samples to quantify the fidelity of the reconstruction. We use fixed number density galaxy samples selected according to stellar mass or star formation rate (SFR). We develop three different schemes to populate haloes with galaxies with increasing complexity, considering the scatter of the satellite HOD as an additional parameter in the modelling. We first modify the SAM output, removing assembly bias and using a standard Navarro-Frenk-White density profile for the satellite galaxies as the target to reproduce with our HOD mocks. We find that all models give similar reproductions of the two-halo contribution to the clustering signal, but there are differences in the one-halo term. In particular, the HOD mock reproductions work equally well using either the HOD of central and satellites separately or using a model that also accounts for whether or not the haloes contain a central galaxy. We find that the HOD scatter does not have an important impact on the clustering predictions for stellar mass selected samples. For SFR selections, we obtain the most accurate results assuming a negative binomial distribution for the number of satellites in a halo. The scatter in the satellites HOD is a key consideration for HOD mock catalogues that mimic ELG or SFR selected samples in future galaxy surveys. △ Less

Submitted 10 June, 2019; originally announced June 2019.

How to add massive neutrinos to your $Λ$CDM simulation -- extending cosmology rescaling algorithms

Authors: Matteo Zennaro, Raúl E. Angulo, Giovanni Aricò, Sergio Contreras, Marcos Pellejero-Ibáñez

Abstract: Providing accurate predictions for the spatial distribution of matter and luminous tracers in the presence of massive neutrinos is an important task, given the imminent arrival of highly accurate large-scale structure observations. In this work, we address this challenge by extending cosmology-rescaling algorithms to massive neutrino cosmologies. In this way, a $Λ$CDM simulation can be modified to… ▽ More Providing accurate predictions for the spatial distribution of matter and luminous tracers in the presence of massive neutrinos is an important task, given the imminent arrival of highly accurate large-scale structure observations. In this work, we address this challenge by extending cosmology-rescaling algorithms to massive neutrino cosmologies. In this way, a $Λ$CDM simulation can be modified to provide nonlinear structure formation predictions in the presence a hot component of arbitrary mass, and, if desired, to include non-gravitational modifications to the clustering of matter on large scales. We test the accuracy of the method by comparing its predictions to a suite of simulations carried out explicitly including a neutrino component in its evolution equations. We find that, for neutrino masses in the range $M_ν\in [0.06, 0.3] ~ \mathrm{eV}$ the matter power spectrum is recovered to better than $1\%$ on all scales $k<2~h~\mathrm{Mpc}^{-1}$. Similarly, the halo mass function is predicted at a few percent level over the range $M_{\rm halo} \in [10^{12}, 10^{15}] ~ h^{-1} ~ \mathrm{M}_{\odot}$, and so do also the multipoles of the galaxy 2-point correlation function in redshift space over $r \in [0.1, 200] ~ h^{-1} ~ \mathrm{Mpc}$. We provide parametric forms for the necessary transformations, as a function of $Ω_{\rm m}$ and $Ω_ν$ for various target redshifts. △ Less

Submitted 3 December, 2019; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: 15 pages, 11 figures

Warm CO in evolved stars from the THROES catalogue. II. Herschel/PACS spectroscopy of C-rich envelopes

Authors: J. M. da Silva Santos, J. Ramos-Medina, C. Sánchez Contreras, P. García-Lario

Abstract: This is the second paper of a series making use of Herschel/PACS spectroscopy of evolved stars in the THROES catalogue to study the inner regions of their circumstellar envelopes (CSEs). We analyze the CO emission spectra, including a large number of high-$J$ CO lines (from $J$=14-13 to $J$=45-44), as a proxy for the warm molecular gas in the CSEs of a sample of bright carbon-rich stars spanning d… ▽ More This is the second paper of a series making use of Herschel/PACS spectroscopy of evolved stars in the THROES catalogue to study the inner regions of their circumstellar envelopes (CSEs). We analyze the CO emission spectra, including a large number of high-$J$ CO lines (from $J$=14-13 to $J$=45-44), as a proxy for the warm molecular gas in the CSEs of a sample of bright carbon-rich stars spanning different evolutionary stages from the Asymptotic Giant Branch (AGB) to the young planetary nebulae (PNe) phase. We use the rotational diagram (RD) to derive rotational temperatures ($T_{\rm rot}$) and masses ($M_{\rm H_2}$) of the envelope layers where the CO transitions arise. We also obtain a first order estimate of the mass-loss rates and assess the impact of the opacity correction for a range of characteristic envelope radii. We use multi-epoch spectra for the well studied C-rich envelope IRC+10216 to investigate the impact of CO flux variability on the values of $T_{\rm rot}$ and $M_{\rm H_2}$. PACS sensitivity allowed the study of higher rotational numbers than before, indicating the presence of a significant amount of warmer gas ($\sim$200-900 K) not traceable with lower-$J$ CO observations at sub-mm/mm wavelengths. The masses are in the range $\sim10^{-2}-10^{-5}\,\rm M_{\odot}$, anti-correlated with temperature. For some strong CO emitters we infer a double temperature (warm $T_{\rm rot}\sim$400 K and hot $T_{\rm rot}\sim$820 K) component. From the analysis of IRC+10216, we corroborate that the effect of line variability is perceptible on the $T_{\rm rot}$ of the hot component only, and certainly insignificant on $M_{\rm H_2}$ and, hence, the mass-loss rate. Therefore, the parameters derived from the RD are robust even when strong line flux variability occurs, with the major source of uncertainty in the estimate of the mass-loss rate being the size of the CO-emitting volume. △ Less

Submitted 19 December, 2018; originally announced December 2018.

Comments: accepted in A&A

Journal ref: A&A 622, A123 (2019)

The Effect of Assembly Bias on Redshift Space Distortions

Authors: Nelson Padilla, Sergio Contreras, Idit Zehavi, Carlton Baugh, Peder Norberg

Abstract: We study potential systematic effects of assembly bias on cosmological parameter constraints from redshift space distortion measurements. We use a semi-analytic galaxy formation model applied to the Millennium N-body WMAP-7 simulation to study the effects of halo assembly bias on the redshift space distortions of the galaxy correlation function. We look at the pairwise velocities of galaxies livin… ▽ More We study potential systematic effects of assembly bias on cosmological parameter constraints from redshift space distortion measurements. We use a semi-analytic galaxy formation model applied to the Millennium N-body WMAP-7 simulation to study the effects of halo assembly bias on the redshift space distortions of the galaxy correlation function. We look at the pairwise velocities of galaxies living in haloes with concentrations and ages in the upper and lower quintiles, and find that the velocity differences between these are consistent with those reported for real-space clustering analyses, i.e. samples with higher clustering also exhibit stronger infall pairwise motions. This can also be seen in the monopole and quadrupole of the redshift-space correlation function. We find that regardless of the method of measurement, the changes in the $β$ parameter due to different secondary halo parameters fully tracks the change in the bias Parameter. Hence, assembly bias does not introduce detectable systematics in the inferred logarithmic growth factor. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: 15 pages, 12 figures.Submitted to MNRAS

Through the magnifying glass: ALMA acute viewing of the intricate nebular architecture of OH231.8+4.2

Authors: C. Sánchez Contreras, J. Alcolea, A. Castro-Carrizo, L. Velilla-Prieto, M. Santander-García, G. Quintana-Lacaci, J. Cernicharo

Abstract: We present continuum and molecular line emission ALMA observations of OH 231.8+4.2, a well studied bipolar nebula around an asymptotic giant branch (AGB) star. The high angular resolution (~0.2-0.3 arcsec) and sensitivity of our ALMA maps provide the most detailed and accurate description of the overall nebular structure and kinematics of this object to date. We have identified a number of outflow… ▽ More We present continuum and molecular line emission ALMA observations of OH 231.8+4.2, a well studied bipolar nebula around an asymptotic giant branch (AGB) star. The high angular resolution (~0.2-0.3 arcsec) and sensitivity of our ALMA maps provide the most detailed and accurate description of the overall nebular structure and kinematics of this object to date. We have identified a number of outflow components previously unknown. Species studied in this work include 12CO, 13CO, CS, SO, SO2, OCS, SiO, SiS, H3O+, Na37Cl, and CH3OH. The molecules Na37Cl and CH3OH are first detections in OH 231.8+4.2, with CH3OH being also a first detection in an AGB star. Our ALMA maps bring to light the totally unexpected position of the mass-losing AGB star (QX Pup) relative to the large-scale outflow. QX Pup is enshrouded within a compact (<60 AU) parcel of dust and gas (clump S) in expansion (V~5-7 km/s) that is displaced by 0.6arcsec to the south of the dense equatorial region (or waist) where the bipolar lobes join. Our SiO maps disclose a compact bipolar outflow that emerges from QX Pup's vicinity. This outflow is oriented similarly to the large-scale nebula but the expansion velocities are about ten times lower (~35 km/s). We deduce short kinematical ages for the SiO outflow, ranging from ~50-80 yr, in regions within ~150 AU, to ~400-500 yr at the lobe tips (~3500 AU). Adjacent to the SiO outflow, we identify a small-scale hourglass-shaped structure (mini-hourglass) that is probably made of compressed ambient material formed as the SiO outflow penetrates the dense, central regions of the nebula. The lobes and the equatorial waist of the mini-hourglass are both radially expanding with a constant velocity gradient. The mini-waist is characterized by extremely low velocities, down to ~1 km/s at ~150 AU, which tentatively suggest the presence of a stable structure. (abridged) △ Less

Submitted 3 September, 2018; originally announced September 2018.

The Evolution of Assembly Bias

Authors: Sergio Contreras, Idit Zehavi, Nelson Padilla, Carlton Baugh, Esteban Jiménez, Ivan Lacerna

Abstract: We examine the evolution of assembly bias using a semi-analytical model of galaxy formation implemented in the Millennium-WMAP7 N-body simulation. We consider fixed number density galaxy samples ranked by stellar mass or star formation rate. We investigate how the clustering of haloes and their galaxy content depend on halo formation time and concentration, and how these relationships evolve with… ▽ More We examine the evolution of assembly bias using a semi-analytical model of galaxy formation implemented in the Millennium-WMAP7 N-body simulation. We consider fixed number density galaxy samples ranked by stellar mass or star formation rate. We investigate how the clustering of haloes and their galaxy content depend on halo formation time and concentration, and how these relationships evolve with redshift. At $z=0$ the dependences of halo clustering on halo concentration and formation time are similar. However, at higher redshift, halo assembly bias weakens for haloes selected by age, and reverses and increases for haloes selected by concentration. The variation of the halo occupation with concentration and formation time is also similar at $z=0$ and changes at higher redshifts. In this case, the occupancy variation with halo age stays mostly constant with redshift but decreases for concentration. Finally, we look at the evolution of assembly bias reflected in the galaxy distribution by examining the galaxy correlation functions relative to those of shuffled galaxy samples which remove the occupancy variation. This correlation functions ratio monotonically decreases with larger redshift and for lower number density samples, going below unity in some cases, leading to reduced galaxy clustering. While the halo occupation functions themselves vary, the assembly bias trends are similar whether selecting galaxies by stellar mass or star formation rate. Our results provide further insight into the origin and evolution of assembly bias. Our extensive occupation function measurements and fits are publicly available and can be used to create realistic mock catalogues. △ Less

Submitted 14 February, 2019; v1 submitted 8 August, 2018; originally announced August 2018.

Comments: Accepted by MNRAS, 15 pages, 13 figures

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

Authors: M. Celeste Artale, Idit Zehavi, Sergio Contreras, Peder Norberg

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 environmen… ▽ More 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. △ Less

Submitted 1 August, 2018; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: Accepted for publication in MNRAS, 17 pages, 11 figures

The maser emitting structure and time variability of the SiS lines J=14-13 and 15-14 in IRC+10216

Authors: J. P. Fonfría, M. Fernández-López, J. R. Pardo, M. Agúndez, C. Sánchez Contreras, L. Velilla Prieto, J. Cernicharo, M. Santander-García, G. Quintana-Lacaci, A. Castro-Carrizo, S. Curiel

Abstract: We present new high angular resolution interferometer observations of the v=0 J=14-13 and 15-14 SiS lines towards IRC+10216, carried out with CARMA and ALMA. The maps, with angular resolutions of ~0.25"and 0.55", reveal (1) an extended, roughly uniform, and weak emission with a size of ~0.5", (2) a component elongated approximately along the East-West direction peaking at ~0.13" and 0.17" at both… ▽ More We present new high angular resolution interferometer observations of the v=0 J=14-13 and 15-14 SiS lines towards IRC+10216, carried out with CARMA and ALMA. The maps, with angular resolutions of ~0.25"and 0.55", reveal (1) an extended, roughly uniform, and weak emission with a size of ~0.5", (2) a component elongated approximately along the East-West direction peaking at ~0.13" and 0.17" at both sides of the central star, and (3) two blue- and red-shifted compact components peaking around 0.07" to the NW of the star. We have modeled the emission with a 3D radiation transfer code finding that the observations cannot be explained only by thermal emission. Several maser clumps and one arc-shaped maser feature arranged from 5 to 20R* from the central star, in addition to a thin shell-like maser structure at ~13R* are required to explain the observations. This maser emitting set of structures accounts for 75% of the total emission while the other 25% is produced by thermally excited molecules. About 60% of the maser emission comes from the extended emission and the rest from the set of clumps and the arc. The analysis of a time monitoring of these and other SiS and 29SiS lines carried out with the IRAM 30m telescope from 2015 to present suggests that the intensity of some spectral components of the maser emission strongly depends on the stellar pulsation while other components show a mild variability. This monitoring evidences a significant phase lag of ~0.2 between the maser and NIR light-curves. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 19 pages, 11 figures, 3 tables; accepted for publication in The ApJ

Binarity and Accretion in AGB Stars: HST/STIS Observations of UV Flickering in Y Gem

Authors: R. Sahai, C. Sánchez Contreras, A. Mangan, J. Sanz-Forcada, C. Muthumariappan, M. J. Claussen

Abstract: Binarity is believed to dramatically affect the history and geometry of mass loss in AGB and post-AGB stars, but observational evidence of binarity is sorely lacking. As part of a project to search for hot binary companions to cool AGB stars using the GALEX archive, we discovered a late-M star, Y Gem, to be a source of strong and variable UV and X-ray emission. Here we report UV spectroscopic obse… ▽ More Binarity is believed to dramatically affect the history and geometry of mass loss in AGB and post-AGB stars, but observational evidence of binarity is sorely lacking. As part of a project to search for hot binary companions to cool AGB stars using the GALEX archive, we discovered a late-M star, Y Gem, to be a source of strong and variable UV and X-ray emission. Here we report UV spectroscopic observations of Y Gem obtained with the Hubble Space Telescope that show strong flickering in the UV continuum on time-scales of <~20 s, characteristic of an active accretion disk. Several UV lines with P-Cygni-type profiles from species such as Si IV and C IV are also observed, with emission and absorption features that are red- and blue- shifted by velocities of ~500 km/s from the systemic velocity. Our model for these (and previous) observations is that material from the primary star is gravitationally captured by a companion, producing a hot accretion disk. The latter powers a fast outflow that produces blue-shifted features due to absorption of UV continuum emitted by the disk, whereas the red-shifted emission features arise in heated infalling material from the primary. The outflow velocities support a previous inference by Sahai et al. (2015) that Y Gem's companion is a low-mass main-sequence star. Blackbody fitting of the UV continuum implies an accretion luminosity of about 13 Lsun and thus a mass-accretion rate >5e-7 Msun/yr; we infer that Roche lobe overflow is the most likely binary accretion mode for Y Gem. △ Less

Submitted 8 May, 2018; originally announced May 2018.

High-resolution observations of IRAS 08544-4431. Detection of a disk orbiting a post-AGB star and of a slow disk wind

Authors: V. Bujarrabal, A. Castro-Carrizo, H. Van Winckel, J. Alcolea, C. Sanchez Contreras, M. Santander-Garcia, M. Hillen

Abstract: We are studying a class of binary post-AGB stars that seem to be systematically surrounded by equatorial disks and slow outflows. Although the rotating dynamics had only been well identified in three cases, the study of such structures is thought to be fundamental to the understanding of the formation of nebulae around evolved stars. We present ALMA maps of 12CO and 13CO J=3-2 lines in one of thes… ▽ More We are studying a class of binary post-AGB stars that seem to be systematically surrounded by equatorial disks and slow outflows. Although the rotating dynamics had only been well identified in three cases, the study of such structures is thought to be fundamental to the understanding of the formation of nebulae around evolved stars. We present ALMA maps of 12CO and 13CO J=3-2 lines in one of these sources, IRAS08544-4431. We analyzed the data by means of nebula models, which account for the expectedly composite source and can reproduce the data. From our modeling, we estimated the main nebula parameters, including the structure and dynamics and the density and temperature distributions. We discuss the uncertainties of the derived values and, in particular, their dependence on the distance. Our observations reveal the presence of an equatorial disk in rotation; a low-velocity outflow is also found, probably formed of gas expelled from the disk. The main characteristics of our observations and modeling of IRAS08544-4431 are similar to those of better studied objects, confirming our interpretation. The disk rotation indicates a total central mass of about 1.8 Mo, for a distance of 1100 pc. The disk is found to be relatively extended and has a typical diameter of ~ 4 10^16 cm. The total nebular mass is ~ 2 10^-2 Mo, of which ~ 90% corresponds to the disk. Assuming that the outflow is due to mass loss from the disk, we derive a disk lifetime of ~ 10000 yr. The disk angular momentum is found to be comparable to that of the binary system at present. Assuming that the disk angular momentum was transferred from the binary system, as expected, the high values of the disk angular momentum in this and other similar disks suggest that the size of the stellar orbits has significantly decreased as a consequence of disk formation. △ Less

Submitted 12 February, 2018; originally announced February 2018.

Comments: 10 pages, 7 figures, accepted by A&A

Journal ref: A&A 614, A58 (2018)